R-如何将迭代从引导程序保存到列表中



我有一个看起来像:

的数据框架
df <- data.frame(id=c(rep("no.234",20),rep("no.125",15)),y=rnorm(35))

我想在每个ID上使用贝叶斯引导程序,例如1000次,然后使用循环或管道将结果放在列表中。1个ID的过程是:

require(bayesboot)
require(dplyr)
no.234 <- df %>% filter(id=="no.234")
bb <- bayesboot(no.234$y, R=1000)

使用sapply循环浏览所有IDS。

x <- sapply(levels(df$id), 
     FUN=function(x) bayesboot(df$y[df$id==x], statistic=mean, R=1000))

结果将是一个列表,其中每个元素将包含给定ID的自举统计。

最新更新