[R-meta] Replicating robumeta package output by metafor package in R

James Pustejovsky jepu@to @end|ng |rom gm@||@com
Tue Nov 5 04:45:01 CET 2019


Both packages use the same method to calculate standard errors and
covariances for meta-analytic average effect size estimates and
meta-regression coefficients. However, the packages implement different
estimation methods for between-study heterogeneity and they use different
weighting schemes when estimating meta-regression coefficients. Some
details and an example can be found here:


On Mon, Nov 4, 2019 at 3:15 PM Reza Norouzian <rnorouzian using gmail.com> wrote:

> I was wondering if there might be a way to replicate the output from
> library robumenta by library metafor in R?
> Both package state that they use the same estimation method for robust
> variance estimation.
> *Here is reproducible data and the code I tried with no success:*
> d <- read.csv("https://raw.githubusercontent.com/izeh/m/master/e.csv",
> h = T) # DATA
> library(robumeta)
> robu(y~factor(time), data = d, studynum = study.name, var = SD^2, small =
> F)
> library(metafor)
> robust(rma(y~factor(time), sei = SD, data = d), cluster = d$id, adjust = F)
>         [[alternative HTML version deleted]]
> _______________________________________________
> R-sig-meta-analysis mailing list
> R-sig-meta-analysis using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis

	[[alternative HTML version deleted]]

More information about the R-sig-meta-analysis mailing list