【发布时间】：2015-03-27 15:39:12
【问题描述】：
我有两个问题：
首先，我是一个定性变量类型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]])以便函数从表和数据框中提取相应的值。
我必须像这样使用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)
}
}
【问题讨论】：
-
嗨，夏洛特，虽然有些过时，但此页面有一些代码可能对您有帮助 cookbook-r.com/Graphs/Plotting_means_and_error_bars_(ggplot2) 特别是 summarySE 函数。
-
“summarySE”包不可用（作为 R 版本 3.1.3 的二进制包）
标签： r function confidence-interval