R 计算置信区间的函数

R 计算置信区间的函数,r,function,confidence-interval,R,Function,Confidence Interval,我有两个问题: 首先,我需要一个定性变量Typeclass(Type)=table,具有不同的类别频率,比如 Type Type1 Type2 Type3 150 4900 4 以及我之前在另一个函数中计算的每个类别的统计(类是数据帧): 我想在一个函数中自动计算每个类别的置信区间,该函数返回置信区间上下限的两个向量。置信区间的公式为: df[1,2]+1.96*sqrt(df[1

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

Type
    Type1          Type2            Type3
      150           4900                4  
以及我之前在另一个函数中计算的每个类别的统计(类是数据帧):

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

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)
    }
}

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

Hi Charlotte,尽管有些过时,但此页面的一些代码可能对您有所帮助,尤其是summarySE函数。包“summarySE”不可用(作为R版本3.1.3的二进制包)
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)
    }
}