按组提取引导统计信息答案

【问题标题】：extracting bootstrap statistics by group按组提取引导统计信息
【发布时间】：2013-10-16 13:14:17
【问题描述】：

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

有没有一种简单的方法可以以下表的形式从输出中提取引导统计信息（平均值和标准误差）？

    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

谢谢，

【问题讨论】：

您好，在使用随机数生成器生成样本数据时，set.seed() 很有帮助
@RicardoSaporta 感谢您的提示。

标签： r statistics-bootstrap

【解决方案1】：

看一下bootstr(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

【讨论】：

很奇怪。查看getAnywhere(print.boot)，我发现它实际上是在动态计算这些。