[datatable-help] random factors in rma.mv (metafor)
liso
lholper at gmx.ch
Sat Feb 10 15:07:29 CET 2018
Hi,
I have a question regarding the use of rma.mv (metafor) in a multivariate
meta-analysis.
In particular, I have outcomes that vary not only between trials but also
between subjects. In other words, there are outcomes that have been measured
twice in one subject using different methods (mRNA, protein). This is true
for only a few studies, while others only report one method.
Hereby, I would like to deal with it?
So far, I have one random variable 'trial'.
Please find attached the code and a test dataset.
https://www.dropbox.com/s/gcn7e7e89g5nfr9/test.xlsx?dl=0
Thank you very much in advance.
Lisa
> rm(list = ls())
> library(metafor)
> library(readxl)
> dat1 <-
> read_excel("Documents/00_Data/12_MetaAnalysis/Data/test.xlsx",sheet =
> "BD_sorted")
>
>
> yi <- dat1$SMD
> V <- dat1$Variance
> W<-dat1$NoPatients
> res <- rma.mv(yi, V, W,mods = ~ Gen - 1, random =~ Gen | trial,
> struct="UN", data=dat1)
> res
Multivariate Meta-Analysis Model (k = 53; method: REML)
Variance Components:
outer factor: trial (nlvls = 11)
inner factor: Gen (nlvls = 4)
estim sqrt k.lvl fixed level
tau^2.1 0.2331 0.4828 5 no NDUFA1
tau^2.2 0.0009 0.0303 10 no NDUFS1
tau^2.3 0.0147 0.1212 17 no NDUFV1
tau^2.4 0.0595 0.2440 21 no NDUFV2
rho.NDUFA rho.NDUFS rho.NDUFV1 rho.NDUFV2 NDUFA NDUFS NDUFV1
NDUFV2
NDUFA1 1 0.5679 0.9493 0.9989 - no no
no
NDUFS1 0.5679 1 0.2804 0.6060 1 - no
no
NDUFV1 0.9493 0.2804 1 0.9335 1 4 -
no
NDUFV2 0.9989 0.6060 0.9335 1 2 4 5
-
Test for Residual Heterogeneity:
QE(df = 49) = 119.9290, p-val < .0001
Test of Moderators (coefficient(s) 1:4):
QM(df = 4) = 11.8456, p-val = 0.0185
Model Results:
estimate se zval pval ci.lb ci.ub
GenNDUFA1 0.4482 0.3377 1.3272 0.1845 -0.2137 1.1100
GenNDUFS1 -0.1999 0.0641 -3.1211 0.0018 -0.3255 -0.0744 **
GenNDUFV1 0.0145 0.0779 0.1857 0.8527 -0.1382 0.1671
GenNDUFV2 -0.0109 0.1272 -0.0854 0.9319 -0.2601 0.2383
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
--
Sent from: http://r.789695.n4.nabble.com/datatable-help-f2315188.html
More information about the datatable-help
mailing list