Dear all, I use R 3.0 for Windows. I ask how escalc(measure="SMCC") [metafor package] mathematically calculate yi and vi when only change score and SD change score are provided. I used the example (repoted below) posted by Qiang Yue at: http://r.789695.n4.nabble.com/using-metafor-for-meta-analysis-of-before-after-studies-escalc-SMCC-td4667233.html but it is unclear for me the formula used to derive y1 and v1. I read the documentation of metafor package and it is all very well described, but I ask if possible for the formula used by escalc(measure="SMCC") to mathematically calculate yi and vi. Unfortunatel, I have no free access to the paper quoted in metafor package. Beginning example: fMRS author year n mean_r sd_r mean_s sd_s r 1 Tom 2006 9 0 0 0.12 0.03 0 0 2 Jack 2012 6 0 0 0.23 0.05 0 0 3 Zhu 2013 8 0 0 0.18 0.05 0 0> dat_SMCC=escalc(measure="SMCC",data=fMRS,ni=n,m1i=mean_s,m2i=mean_r,sd1i=sd_s,sd2i=sd_r,ri=r,append=TRUE)> dat_SMCCauthor year n mean_r sd_r mean_s sd_s r yi vi 1 Tom 2006 9 0 0 0.12 0.03 0 3.6108 0.8354 2 Jack 2012 6 0 0 0.23 0.05 0 3.8674 1.4131 3 Zhu 2013 8 0 0 0.18 0.05 0 3.1975 0.7640 Sincerely -- Mario Petretta Department of Translational Medical Sciences Naples University Federico II Italy
Dear Mario, You can always just inspect the code: escalc.default 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 petretta at unina.it > Sent: Tuesday, November 19, 2013 11:29 > To: r-help at r-project.org > Subject: [R] metafor escalc(measure="SMCC") > > Dear all, > > I use R 3.0 for Windows. > > I ask how escalc(measure="SMCC") [metafor package] mathematically > calculate yi and vi when only change score and SD change score are > provided. > > I used the example (repoted below) posted by Qiang Yue at: > > http://r.789695.n4.nabble.com/using-metafor-for-meta-analysis-of-before- > after-studies-escalc-SMCC-td4667233.html > > but it is unclear for me the formula used to derive y1 and v1. I read > the documentation of metafor package and it is all very well > described, but I ask if possible for the formula used by > escalc(measure="SMCC") to mathematically calculate yi and vi. > Unfortunatel, I have no free access to the paper quoted in metafor > package. > > Beginning example: > > fMRS > author year n mean_r sd_r mean_s sd_s r > 1 Tom 2006 9 0 0 0.12 0.03 0 0 > 2 Jack 2012 6 0 0 0.23 0.05 0 0 > 3 Zhu 2013 8 0 0 0.18 0.05 0 0 > > > > dat_SMCC=escalc(measure="SMCC",data=fMRS,ni=n,m1i=mean_s,m2i=mean_r,sd1i=s > d_s,sd2i=sd_r,ri=r > ,append=TRUE) > > dat_SMCC > author year n mean_r sd_r mean_s sd_s r yi vi > 1 Tom 2006 9 0 0 0.12 0.03 0 3.6108 0.8354 > 2 Jack 2012 6 0 0 0.23 0.05 0 3.8674 1.4131 > 3 Zhu 2013 8 0 0 0.18 0.05 0 3.1975 0.7640 > > > > Sincerely > > -- > Mario Petretta > Department of Translational Medical Sciences > Naples University Federico II > Italy > > ______________________________________________ > 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.
Many thanks to Wolfgang Viechtbauer for the prompt and clear answer. However, I'm unable to understand what .cmicalc mathematically does in the code line: cmi <- .cmicalc(mi) I look in the metafor documentation and in R (?.cmical and ??.cmicalc) but I have no result. Please, can I have further explanation on this point? Sorry for the trouble Sincerely Mario Petretta Department of Translational Medical Sciences Naples University Federico II Italy>Message: 50 Date: Tue, 19 Nov 2013 11:48:33 +0100 From: "Viechtbauer Wolfgang (STAT)" <wolfgang.viechtbauer at maastrichtuniversity.nl> To: "petretta at unina.it" <petretta at unina.it>, "r-help at r-project.org" <r-help at r-project.org> Subject: Re: [R] metafor escalc(measure="SMCC") Message-ID: <077E31A57DA26E46AB0D493C9966AC730D9925568A at UM-MAIL4112.unimaas.nl> Content-Type: text/plain; charset="us-ascii" Dear Mario, You can always just inspect the code: escalc.default 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 petretta at unina.it > Sent: Tuesday, November 19, 2013 11:29 > To: r-help at r-project.org > Subject: [R] metafor escalc(measure="SMCC") > > Dear all, > > I use R 3.0 for Windows. > > I ask how escalc(measure="SMCC") [metafor package] mathematically > calculate yi and vi when only change score and SD change score are > provided. > > I used the example (repoted below) posted by Qiang Yue at: > > http://r.789695.n4.nabble.com/using-metafor-for-meta-analysis-of-before- > after-studies-escalc-SMCC-td4667233.html > > but it is unclear for me the formula used to derive y1 and v1. I read > the documentation of metafor package and it is all very well > described, but I ask if possible for the formula used by > escalc(measure="SMCC") to mathematically calculate yi and vi. > Unfortunatel, I have no free access to the paper quoted in metafor > package. > > Beginning example: > > fMRS > author year n mean_r sd_r mean_s sd_s r > 1 Tom 2006 9 0 0 0.12 0.03 0 0 > 2 Jack 2012 6 0 0 0.23 0.05 0 0 > 3 Zhu 2013 8 0 0 0.18 0.05 0 0 > > > > dat_SMCC=escalc(measure="SMCC",data=fMRS,ni=n,m1i=mean_s,m2i=mean_r,sd1i=s > d_s,sd2i=sd_r,ri=r > ,append=TRUE) > > dat_SMCC > author year n mean_r sd_r mean_s sd_s r yi vi > 1 Tom 2006 9 0 0 0.12 0.03 0 3.6108 0.8354 > 2 Jack 2012 6 0 0 0.23 0.05 0 3.8674 1.4131 > 3 Zhu 2013 8 0 0 0.18 0.05 0 3.1975 0.7640 > > > > Sincerely > > -- > Mario Petretta > Department of Translational Medical Sciences > Naples University Federico II > Italy >
Many thanks to Wolfgang Wiechtbauer for the explanation.
Sincerely
Mario Petretta
Message: 31
Date: Thu, 21 Nov 2013 19:08:14 +0100
From: "Viechtbauer Wolfgang (STAT)"
<wolfgang.viechtbauer at maastrichtuniversity.nl>
To: "petretta at unina.it" <petretta at unina.it>,
"r-help at r-project.org"
<r-help at r-project.org>
Subject: Re: [R] metafor escalc(measure="SMCC")
Message-ID:
<077E31A57DA26E46AB0D493C9966AC730D99255ED7 at
UM-MAIL4112.unimaas.nl>
Content-Type: text/plain; charset="us-ascii"
.cmicalc is a non-exported function. You can see the code with:
getAnywhere(.cmicalc)
Best,
Wolfgang
[Nascondi Testo quotato]
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On Behalf Of petretta at unina.it
Sent: Thursday, November 21, 2013 13:39
To: r-help at r-project.org
Subject: [R] metafor escalc(measure="SMCC")
Many thanks to Wolfgang Viechtbauer for the prompt and clear answer.
However, I'm unable to understand what .cmicalc mathematically does in the
code line:
cmi <- .cmicalc(mi)
I look in the metafor documentation and in R (?.cmical and ??.cmicalc) but
I have no result. Please, can I have further explanation on this point?
Sorry for the trouble
Sincerely
Mario Petretta
Department of Translational Medical Sciences
Naples University Federico II
Italy
Message: 50 Date: Tue, 19 Nov 2013 11:48:33 +0100
From: "Viechtbauer Wolfgang (STAT)"
<wolfgang.viechtbauer at maastrichtuniversity.nl>
To: "petretta at unina.it" <petretta at unina.it>,
"r-help at r-project.org"
<r-help at r-project.org>
Subject: Re: [R] metafor escalc(measure="SMCC")
Message-ID:
<077E31A57DA26E46AB0D493C9966AC730D9925568A at UM-
MAIL4112.unimaas.nl>
Content-Type: text/plain; charset="us-ascii"
Dear Mario,
You can always just inspect the code:
escalc.default
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 petretta at unina.it
Sent: Tuesday, November 19, 2013 11:29
To: r-help at r-project.org
Subject: [R] metafor escalc(measure="SMCC")
Dear all,
I use R 3.0 for Windows.
I ask how escalc(measure="SMCC") [metafor package] mathematically
calculate yi and vi when only change score and SD change score are
provided.
I used the example (repoted below) posted by Qiang Yue at:
http://r.789695.n4.nabble.com/using-metafor-for-meta-analysis-of-before-
after-studies-escalc-SMCC-td4667233.html
but it is unclear for me the formula used to derive y1 and v1. I read
the documentation of metafor package and it is all very well
described, but I ask if possible for the formula used by
escalc(measure="SMCC") to mathematically calculate yi and vi.
Unfortunatel, I have no free access to the paper quoted in metafor
package.
Beginning example:
fMRS
author year n mean_r sd_r mean_s sd_s r
1 Tom 2006 9 0 0 0.12 0.03 0 0
2 Jack 2012 6 0 0 0.23 0.05 0 0
3 Zhu 2013 8 0 0 0.18 0.05 0 0
dat_SMCC=escalc(measure="SMCC",data=fMRS,ni=n,m1i=mean_s,m2i=mean_r,sd1i=s
d_s,sd2i=sd_r,ri=r
,append=TRUE)
dat_SMCC
author year n mean_r sd_r mean_s sd_s r yi vi 1 Tom 2006 9 0
0 0.12 0.03 0 3.6108 0.8354
2 Jack 2012 6 0 0 0.23 0.05 0 3.8674 1.4131
3 Zhu 2013 8 0 0 0.18 0.05 0 3.1975 0.7640
Sincerely
--
Mario Petretta
Department of Translational Medical Sciences
Naples University Federico II
Italy
______________________________________________
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.
--
Mario Petretta
Department of Translational Medical Sciences
Naples University Federico II
Italy