Dear Owen, What is your definition of "multivariate analysis"? Do you mean: A meta-regression model with more than one predictor/moderator? In that case, yes, metafor handles that. Usually, this is referred to as "multiple regression" (as opposed to "simple regression" with a single predictor) -- and in the case of a meta-analysis, I guess one could call it "multiple meta-regression".
If you are referring to a model that handles statistical dependencies in the observed outcomes (and hence requires multivariate methods), then you will have to use some other package (e.g., mvmeta). See also: http://stats.stackexchange.com/questions/2358/explain-the-difference-between-multiple-regression-and-multivariate-regression Best, Wolfgang > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of Owen, Branwen > Sent: Friday, March 15, 2013 11:23 > To: r-help@r-project.org > Subject: [R] metafor - multivariate analysis > > > > > 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="D > L",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. ______________________________________________ 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.