<div dir="ltr"><div>Hi all,</div><div>I noticed that in Comprehensive Meta Analysis (CMA), there are two options you can select to calculate subgroup average effect sizes. The first approach is assume a common among-study variance component across subgroups and the second one is not assume a common among-study variance component. When you select different approaches, the subgroup average effect sizes, the overall average effect size and their confidence intervals will change accordingly. This has something to do with the tau-squared calculated by different approaches. This<a href="https://www.youtube.com/watch?v=Y7X5ZbfJgDI&t=1714s."> video </a>here explains the reason behind it in detail if you fast forward to 22:40.</div><div><br></div><div>I have attached my dataset here as an example. I want to do a meta-analysis of proportions so the data doesn't have a control group. It includes a moderator called studytype that divides the studies into two subgroups: birth cohort and others.</div><div><br></div><div>Using R, the subgroup average effect sizes I calculated are the same as obtained in CMA. <u>My question is that </u>how to use R to get the same overall average effect size as calculated in CMA using different approaches?</div><div><br></div><div>My code:</div><div>dat=read.csv("data.csv",header=T,sep=",")</div><div><div>es=escalc(xi=cases,ni=total,measure="PLO",data=dat)</div><div><a href="http://overall.es">overall.es</a>=rma(yi,vi,data=es,method="DL",weighted=TRUE)</div><div>pred=predict(<a href="http://overall.es">overall.es</a>,transf=transf.ilogit,digits=5)</div><div>print(pred,digits=5)</div></div><div><br></div><div>The overall average effect size:</div><div> pred   <a href="http://ci.lb">ci.lb</a>   ci.ub   <a href="http://cr.lb">cr.lb</a>   cr.ub</div><div> 0.00042 0.00032 0.00057 0.00013 0.00135 </div><div>The overall average effect size obtained above is derived when you don't do subgrouping. In CMA, if you do subgrouping, the overall effect size will change according to the approach you select. If you assume a common among-study variance, the overall average will become 0.000430 with a 95% CI [0.000307, 0.000602]; if you do not assume that, the overall average will become 0.000453 with a 95% CI [0.000335, 0.000611], which are the values I try to obtain in R<u> (I don't know how to do this, which is why I'm asking this question).</u></div><div><br></div><div>If assuming a common among-study variance component across subgroups:<br></div><div>studytype=rma(yi,vi,data=ies.logit,mods=~studytype,method="DL")</div><div>pred.studytype=predict(studytype,transf=transf.ilogit,digits=5)</div><div>print(pred.studytype[c(1,17)])</div><div>then the subgroup average effect sizes are:</div><div>      pred   <a href="http://ci.lb">ci.lb</a>   ci.ub   <a href="http://cr.lb">cr.lb</a>   cr.ub</div><div>1  0.00034 0.00020 0.00061 0.00008 0.00142 #Birth cohort</div><div>17 0.00049 0.00032 0.00074 0.00012 0.00191 #Others</div><div><br></div><div>If not assuming a common among-study variance component across subgroups:</div><div>birthcohort=rma(yi,vi,data=ies.logit,subset=studytype=="Birth cohort",method="DL")</div><div>others=rma(yi,vi,data=ies.logit,subset=studytype=="Others",method="DL")</div><div>pred.birthcohort=predict(birthcohort,transf=transf.ilogit,digits=5)</div><div>pred.others=predict(others,transf=transf.ilogit,digits=5)</div><div>print(pred.birthcohort,digits=5)</div><div>print(pred.others,digits=5)</div><div>then the subgroup average effect sizes are:</div><div> pred   <a href="http://ci.lb">ci.lb</a>   ci.ub   <a href="http://cr.lb">cr.lb</a>   cr.ub</div><div> 0.00035 0.00016 0.00078 0.00005 0.00274</div><div> pred   <a href="http://ci.lb">ci.lb</a>   ci.ub   <a href="http://cr.lb">cr.lb</a>   cr.ub</div><div> 0.00047 0.00034 0.00065 0.00017 0.00132</div><div><br></div><div>The screen shots I attached here shows the same results I obtained above.</div><div><br></div><div>Please let me know if my question makes sense. It's kind of difficult to explain this very clearly in an email.</div><div><br></div><div>Thank you for your time!</div><div><br></div><div>Cheers,</div><div><br></div><div>Naike</div></div>