Dear Metafor users, I'm conducting a metaanalysis of prevalence of a particular
behaviour based on someone elses' code. I've been labouring under the
impression that this:
summary(rma.1<-rma(yi,vi,mods=cbind(approxmeanage,interviewmethodcode),data=mal,method="DL",knha=F,weighted=F,intercept=T))
is doing the multivariate analysis that i want, but have read that multivariate
analysis can't be done in metafor.
this is the output:
Mixed-Effects Model (k = 22; tau^2 estimator: DL)
logLik Deviance AIC BIC
18.7726 -37.5452 -27.5452 -22.0899
tau^2 (estimate of residual amount of heterogeneity): 0.0106
tau (sqrt of the estimate of residual heterogeneity): 0.1031
Test for Residual Heterogeneity:
QE(df = 18) = 1273.9411, p-val < .0001
Test of Moderators (coefficient(s) 2,3,4):
QM(df = 3) = 11.0096, p-val = 0.0117
Model Results:
estimate se zval pval ci.lb ci.ub
intrcpt 0.4014 0.1705 2.3537 0.0186 0.0671 0.7356 *
continent -0.0206 0.0184 -1.1200 0.2627 -0.0568 0.0155
approxmeanage 0.0076 0.0091 0.8354 0.4035 -0.0102 0.0254
interviewmethodcode -0.0892 0.0273 -3.2702 0.0011 -0.1426 -0.0357 **
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
my questions are
1. what is this line of code?
2. if it isn't multivariate analysis, will i have to use the mvmeta instead.
thanks very much for any help
Branwen
.<http://r.789695.n4.nabble.com/metafor-multivariate-analysis-td4661233.html>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
