if (FALSE) {
require(metafor)
dat <- escalc(measure = "OR", ai = tpos, bi = tneg, ci = cpos, di = cneg,
data = dat.bcg)
# fit meta-regression
robma_dat <- data.frame(
logOR = dat$yi,
se = sqrt(dat$vi),
ablat = dat$ablat,
alloc = dat$alloc
)
fit <- NoBMA.reg(~ ablat + alloc, data = robma_dat,
seed = 1, algorithm = "ss", parallel = TRUE)
residuals(fit)
}
Run the code above in your browser using DataLab