Skip to content
Prev 80 / 5632 Next

[R-meta] Robust function when rows have NA values

Could you provide fully reproducible examples to illustrate the problems you are running into? robust() should work just fine even if there are missings and I^2 is also computed correctly. For example:

library(metafor)
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
dat$yi[1] <- NA
res <- rma(yi, vi, data=dat)
res
robust(res, dat$trial)

Best,
Wolfgang