计算置信区间的函数

时间:2015-03-27 15:39:13

标签: r function confidence-interval

我有两个问题: 首先,我是一个定性变量类型class(Type)=table,具有不同的类别频率,比如

Type
    Type1          Type2            Type3
      150           4900                4  

以及我之前在另一个函数中计算的每个类别的统计量(类是数据框):

df
     Var1           Freq
1     Type1            0.1800
2     Type2            0.1130
3     Type3            0.2500

我想在一个函数中自动计算每个类别的置信区间,该函数返回置信区间上限和下限的两个向量。置信区间的公式为: df[1, 2] + 1.96*sqrt(df[1, 2]*(1 - df[1, 2])/t[[1]])以便函数从表和数据框中提取相应的值。 我必须像这样使用循环吗?

ci <- function(stat, Type){
  p <- tapply(stat, Type, function(x) sum(x)/length(x))
  df <- as.data.frame(as.table(p))
  t <- table(Type)
    for (i in nrow(df)) {
      for (j in t)
    }
}

1 个答案:

答案 0 :(得分:0)

summarySE是包Rmisc中的函数。通过此功能,您可以计算不同类别的CI。当前版本1.5可用于R 3.4.0