


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

