我有一个看起来像:
的数据框架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的自举统计。