我今天在挣扎着一些基本统计数据。我想按组计算比例的标准误差。
这是我按组比例使用的代码：
tapply(dt$v1, dt$v2, function(x){prop.table(table(x))})
这是支撑标准误差的公式，SE = sqrt(p*(1-p)/n)其中p=proportion和n=total population。
我想在上面的代码中包含此公式，可以吗？
先感谢您！
这不是吗？
tapply(dt$v1, dt$v2, function(x){prop <- prop.table(table(x))
SE <- sqrt(prop*(1-prop)/length(x))})
如果您希望两者都返回，请添加一行：
list(prop=prop, SE=SE)
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句