我今天在努力学习一些基本的统计数据。我想按组计算比例的标准误差。
这是我按比例逐个使用的代码:
tapply(dt$v1, dt$v2, function(x){prop.table(table(x))})
这是针对SE = sqrt(p*(1-p)/n)
和p=proportion
的标准错误n=total population
的公式。
我想在上面的代码中包含这个公式,是否可能?
提前谢谢!
答案 0 :(得分:2)
这不会是:
tapply(dt$v1, dt$v2, function(x){prop <- prop.table(table(x))
SE <- sqrt(prop*(1-prop)/length(x))})
如果您想要两者都返回,请添加一行:
list(prop=prop, SE=SE)