[R] Help: Meta-analysis with metacor

Viechtbauer Wolfgang (STAT) wolfgang.viechtbauer at maastrichtuniversity.nl
Tue Nov 20 14:05:09 CET 2012


There are multiple issues here. First of all, you should simplify your code to:

metacor(cor, n, title="title", complab="comparison", outclab="outcome")

The rest of your call to metacor() is just using the defaults, so no need to repeat all of that in your call and apparently one of these arguments is causing some problems.

Second, your vector of "correlations" (cor) includes values below -1, which is impossible if these are supposed to be Pearson product-moment correlation coefficients. 

Finally, note that the sampling variance of an r-to-z transformed correlation is calculated with 1/(n-3), so you get division by zero when n=3. As far as I can tell, metacor() will handle this by setting the weight for these values to 0. It is certainly unusual to meta-analyze correlations, where some of those correlations are based on only 3 observations.

Best,

Wolfgang

--   
Wolfgang Viechtbauer, Ph.D., Statistician   
Department of Psychiatry and Psychology   
School for Mental Health and Neuroscience   
Faculty of Health, Medicine, and Life Sciences   
Maastricht University, P.O. Box 616 (VIJV1)   
6200 MD Maastricht, The Netherlands   
+31 (43) 388-4170 | http://www.wvbauer.com   


> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
> On Behalf Of Catherine Proulx
> Sent: Monday, November 19, 2012 18:55
> To: r-help at r-project.org
> Subject: [R] Help: Meta-analysis with metacor
> 
> 
> Trying to do a meta-analysis of correlations in R using the meta package;
> have tried several things and keep getting a similar error. Can anyone
> help explain the error?
> > cor<-c(-0.3018, 0.667, -3.8002, -0.607, -0.4885, -3.8002, -0.0701,
> 0.1348, -0.9505, -0.5709, -0.6127, -1.2419, -0.1511, -0.1054)> n<-
> c(3,4,3,3,3,3,16,36,30,9,3,3,30,4)> library(meta)> metacor(cor, n,
> data=NULL, subset=NULL, sm="ZCOR", level=0.95, level.comb=level,
> comb.fixed=TRUE, comb.random=TRUE, hakn=NULL, method.tau="DL",
> tau.preset=NULL, TE.tau=NULL, method.bias="linreg", title="title",
> complab="comparison", outclab="outcome")
> Error in data.frame(subset = NULL, comb.fixed = TRUE, comb.random = TRUE,
> :   arguments imply differing number of rows: 0, 1
> Thanks,
> Catherine
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.




More information about the R-help mailing list