我想针对长数据获取按组的Bootstrap统计信息。
从下面的输出中提取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
Thanks,
set.seed()
会很有帮助。 - Ricardo Saporta