Хотя это не предназначено, вы можете переопределить логические внутри rma объекты, которые указывают, что вы установили модель «только для перехвата». Для иллюстрации:
library(metafor)
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
res <- rma(yi, vi, data=dat)
funnel(res)
Итак, вы получите стандартный воронкообразный график. Теперь заставьте int.only быть FALSE. Затем вы получаете воронкообразный график остатков:
res$int.only <- FALSE
funnel(res)
Другой способ сделать это - вручную извлечь остатки и соответствующие отклонения и передать их непосредственно в функцию funnel():
ei <- resid(res)
vei <- diag(vcov(res, type="resid"))
funnel(ei, vei, xlab="Residual Value")