我想按组获取长数据的引导统计信息。

是否有一种简单的方法可以从下表形式的输出中提取 Bootstrap 统计信息(均值和标准误差)？

    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