<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p style="margin-top:0; margin-bottom:0"><br>
</p>
Sorry I was informed by another member that they had to approve my emails each time (even though I have subscribed) so I had taken it out to prevent him from getting notifications. </div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
My understanding is that the "intercept" is equivalent to the univariate model if you were to produce a model using rma.uni for only the first outcome (i.e. total cholesterol). Anything after it would be the "slopes"; The difference from the intercept. </div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
In my example I had given you, the slopes would be HDL and LDL. Therefore, if we take into account HDL and LDL outcomes in a multivariate model. Then, compute the overall intercept for this model using predict(x,newmods=(1,1)) this would give us an "adjusted
intercept". What I mean by adjusted is that we have taken into account the difference between HDL and LDL. Whereas, in a univariate model (rma.uni) this would be without taking into account the difference of HDL and LDL I would refer this to the "unadjusted
model". </div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<span style="font-size: 12pt;"><br>
</span></div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<span style="font-size: 12pt;">I think I may know where the confusion may be coming from. The outcomes I have included in my model are "individual measurements". I had converted these to "cohens D" so that we could include these in a single model (multivariate)
and calculate the covariance between each etc.. In another example, I have combined body mass index (kg/m2), body weight (kg) and waist circumference (cm) by converting them all into "cohens D" so that they can be analysed in a single model (multivariate model).
Then, I can change the order of the intercept to get the overall model for each outcome. Basically, these outcomes are not "groups in a set of data" but individual measurements. </span><br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
Please refer to my presentation I had given to my research centre a few months ago which explains my methodology this show you what I mean. </div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
Hope this makes a little more sense. </div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
I referred to the multivariate model as "trivariate" because it had simply taken into account all "three outcomes". </div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
Kind regards</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
Wasim</div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
<div style="color:rgb(0,0,0)">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> R-sig-meta-analysis <r-sig-meta-analysis-bounces@r-project.org> on behalf of Viechtbauer, Wolfgang (SP) <wolfgang.viechtbauer@maastrichtuniversity.nl><br>
<b>Sent:</b> 18 September 2018 18:14<br>
<b>To:</b> Wasim Iqbal (UG)<br>
<b>Cc:</b> r-sig-meta-analysis@r-project.org<br>
<b>Subject:</b> Re: [R-meta] Inflated confidence intervals</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText">Please always cc the mailing list.<br>
<br>
I don't know what you mean by 'trivariate mean' or 'adjusted mean' (adjusted for what?). If you want the average of the three outcomes, then it would be:<br>
<br>
intrcpt + 1/3 * MeasurementHDL + 1/3 * MeasurementTC<br>
<br>
so in this case:<br>
<br>
-0.2376 + 1/3 * 0.1982 + 1/3 * -0.0148 = -0.1764667<br>
<br>
which you can also get with: predict(resMV, newmods = c(1/3, 1/3)).<br>
<br>
I don't know what you mean by 'Univariate model'. A model that does not distinguish between the three outcomes? Why would this model give an estimate that is equal to the intercept of the model whose output you showed?<br>
<br>
Best,<br>
Wolfgang<br>
<br>
-----Original Message-----<br>
From: Wasim Iqbal (UG) [<a href="mailto:W.Iqbal@newcastle.ac.uk" id="LPlnk81847" class="OWAAutoLink" previewremoved="true">mailto:W.Iqbal@newcastle.ac.uk</a>]
<br>
Sent: Sunday, 16 September, 2018 17:18<br>
To: Viechtbauer, Wolfgang (SP)<br>
Subject: Re: Inflated confidence intervals<br>
<br>
Sorry the difference for both HDL and LDL from the intercept would be....-0.0541. -0.0394 would be just for HDL<br>
<br>
Wasim<br>
________________________________________<br>
From: Wasim Iqbal (UG)<br>
Sent: 16 September 2018 16:13:16<br>
To: Viechtbauer, Wolfgang (SP)<br>
Subject: Re: Inflated confidence intervals <br>
<br>
Sorry I am getting easily confused. Essentially, if you were to work out the trivariate mean (or adjusted mean) for this model would this be the difference of both HDL and LDL from the intercept (i.e. -0.0394)? Therefore, if I was to show this in a table with
the univariate model this would be: <br>
<br>
Univariate model: -0.2376<br>
Trivariate model : -0.0394<br>
<br>
I assumed that we must add the intercept to the slopes of HDL and LDL to get the overall adjusted mean? However, now having second thoughts based on what you said (and my textbooks) by adding the intercept ( i.e. -0.2376+0.1982+-0.0148) I am actually including
the intercept again when it was already included? (silly me). <br>
<br>
Thank you for clarifying. Its the simple things that catch me!<br>
<br>
Kind regards<br>
Wasim<br>
________________________________________<br>
From: Viechtbauer, Wolfgang (SP) <wolfgang.viechtbauer@maastrichtuniversity.nl><br>
Sent: 16 September 2018 14:30:01<br>
To: Wasim Iqbal (UG)<br>
Cc: r-sig-meta-analysis@r-project.org; Gavin Stewart; Chris Seal<br>
Subject: RE: Inflated confidence intervals <br>
<br>
I apparently haven't had enough coffee today, so first a correction on my part:<br>
<br>
predict(resMV, newmods = c(1,0)) gives the estimated (average) outcome for HDL (i.e., -0.2376 + 0.1982 = -0.0394). The coefficient for HDL (i.e., 0.1982) is already the difference between HDL and LDL.<br>
<br>
But you seem to be after something different. Apparently, you want to add the intecept and the two coefficients together, so: -0.2376 + 0.1982 + -0.0148 =~ -0.0541, which indeed you would obtain with predict(resMV, newmods = c(1,1)). But what is the meaning
of this?<br>
<br>
If you want a marginal mean, that is, the average of the three outcomes, then you would want:<br>
<br>
intercept + 1/3 * HDL + 1/3 * TC<br>
<br>
(assuming the intercept corresponds to LDL, as in the output you showed), which you would get with predict(resMV, newmods = c(1/3,1/3)). But maybe I am still misunderstanding.<br>
<br>
Best,<br>
Wolfgang<br>
_______________________________________________<br>
R-sig-meta-analysis mailing list<br>
R-sig-meta-analysis@r-project.org<br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis" id="LPlnk522347" class="OWAAutoLink" previewremoved="true">https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis</a><br>
</div>
</span></font></div>
</div>
</div>
</div>
</body>
</html>