I am trying to calculate proportional meta.analysis with clusters and subgroups. The meta-analysis works fine and as expected, however, forest() seems - at least for me - broken:
forest(res, subgroup=T)
Error in x$labels[[i]] : subscript out of bounds
A reproducible example:
# libraries
library(meta)
library(data.table)
# data
dt = data.table(
ID = rep(1:5, each=3),
AE = rep(c("minor","major","deadly"),5),
ni = rep(c(25,19,101,32,50), each=3),
xi = c(2,3,3,4,1,0,29,13,4,7,4,1,7,2,1) )
# meta-analysis & forest plot
res = metaprop(xi, ni, cluster=ID, subgroup=AE, data=dt, tau.common=T)
forest(res)
A similar issue was described earlier for metamean() ( #41 ), but was solved meanwhile.
Hi everyone.
I am trying to calculate proportional meta.analysis with clusters and subgroups. The meta-analysis works fine and as expected, however, forest() seems - at least for me - broken:
A reproducible example:
A similar issue was described earlier for metamean() ( #41 ), but was solved meanwhile.
Any help would be highly appreciated. Best, Felix