Ignacio López De Ullibarri Galparsoro
2012-Jan-14 19:19 UTC
[R] metafor: weights computation in Mantel-Haenszel method
Dear R users, In metafor 1.6-0, the Mantel-Haenszel method is implemented by the rma.mh() function. I have observed that the sum of the weights computed by weights(x) doesn't add to 100% when x is an object of class rma.mh. The consequences of this fact can be clearly seen when a forest diagram is drawn with forest(x), which calls weights(x) (or more precisely, the method weights.rma.mh() defined in the package). Is this, as I suppose, a bug? [[alternative HTML version deleted]]
Viechtbauer Wolfgang (STAT)
2012-Jan-15 15:39 UTC
[R] metafor: weights computation in Mantel-Haenszel method
Yes, I forgot to divide the weights by the sum of the weights (and multiply by
100) for weights.rma.mh(). This bug was pointed out to me a few weeks ago and
this will be changed in the next version of the package.
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) 368-5248 | http://www.wvbauer.com
________________________________________
From: r-help-bounces at r-project.org [r-help-bounces at r-project.org] On
Behalf Of Ignacio L?pez De Ullibarri Galparsoro [ilu at udc.es]
Sent: Saturday, January 14, 2012 8:19 PM
To: r-help at r-project.org
Subject: [R] metafor: weights computation in Mantel-Haenszel method
Dear R users,
In metafor 1.6-0, the Mantel-Haenszel method is implemented by the rma.mh()
function. I have observed that the sum of the weights computed by weights(x)
doesn't add to 100% when x is an object of class rma.mh. The consequences of
this fact can be clearly seen when a forest diagram is drawn with forest(x),
which calls weights(x) (or more precisely, the method weights.rma.mh() defined
in the package).
Is this, as I suppose, a bug?
[[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.