df_split <- split(testdata, testdata$group)
for (i in 1:7){
model <- metabin(data = df_split[[i]],
event.e = ev.exp, event.c = ev.cont, studlab = df_split[[i]]$author,
n.e = n.exp, n.c = n.cont,
# subgroup = group,
method = "MH", sm = "OR",
common = T, random = F, overall = F , overall.hetstat = F,
test.subgroup = F,)
funnel(model)
title(paste("group",LETTERS[i]))
}
也可以将df_split用做lapply()参数,应该是类似的。
另外tapply()
好像是整合split()
和lapply()
,算是一个函数完成两个操作的