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

[James Pustejovsky]

Reza,

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:
https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2017-September/000223.html


James

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]]