[R-meta] Inner|outer model vs multiple random id terms in rma.mv

Divya Ravichandar d|vy@ @end|ng |rom @econdgenome@com
Wed Apr 22 18:51:37 CEST 2020

Hi all,

I am trying to understand why results from running a model of the form
~lvl1|lv2 are not comparable to results of running ~1 | lvl1 ,~ 1 | lvl2

In a simple example case below,results of the 2 models are comparable as

```case <- data.frame(Dataset=
c("a","b","c","d"),Cohort=c("c1","c1","c2","c3"), Tech=
res1 = rma.mv(Effect_size, Standard_error^2, random = list(~ 1 | Dataset,~
1 | Cohort), data=case)
res2=rma.mv(Effect_size, Standard_error^2, random = ~ Dataset | Cohort,
However, when running the 2 model on a more complex example [attached]
markedly different results are obtained with ~ Dataset | Cohort estimating
a pvalue of .02 and list(~ 1 | Dataset,~ 1 | Cohort) estimating a pvalue of
*Divya Ravichandar*
Second Genome

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20200422/ea0cdcad/attachment.html>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: Case.csv
Type: text/csv
Size: 367 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20200422/ea0cdcad/attachment.csv>

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