我想按组获取长数据的引导统计信息。
有没有一种简单的方法可以从输出中提取引导统计信息（均值和标准误），如下表所示？
data <- data.frame(list(value = runif(300), group = factor(letters[1:4])))
stat <- function(x, i) c(m1 = mean(x$value[i]))
library(plyr);library(boot)
aa=dlply( data, .(group), function( dat ) boot(dat, stat, R=10) )
group| mean| std.error
a |0.51 |0.035
b |0.56 |0.046
c |0.52 |0.034
d |0.61 |0.017
谢谢，
看一下boot返回的对象str(aa[[1]])，然后您将看到如何提取相关统计信息。
ldply(aa, .fun=function(x) data.frame(mean = mean(x$t),
std.err = sd(x$t)))
group mean std.err
1 a 0.5309598 0.04168062
2 b 0.5194311 0.02583568
3 c 0.5064817 0.02791644
4 d 0.4512118 0.03417612
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句