我有两个问题:
首先,我是一个定性变量类型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)
}
}