[R-meta] Mismatch between output from sub-group analysis and forest plot

Joao Afonso jot@|on@o @end|ng |rom gm@||@com
Tue Feb 11 18:51:16 CET 2020

Dear all,

I am conducting a meta-analysis on prevalence and incidence data. As there
is plenty of heterogeneity I am doing a sub-group analysis. However I am
finding the output of the same and forest plots produced to have different
values for the pooled estimates. I am using the double arcsine
transformation and then have the data back-transformed to produce the
estimates. The model is the following:

*iies.da.lcmbi<-ies.da %>%  filter(is.na
<http://is.na>(lcmbi)==FALSE)subganal.lcmbi=rma(yi, vi, data=ies.da.lcmbi,
mods=~lcmbi, method="DL")pes.da.lcm=rma(yi, vi, data=ies.da.lcmbi,
mods=~lcmbi=="Records", method="DL")pes.da.records=rma(yi, vi,
data=ies.da.lcmbi, mods=~lcmbi=="Locomotion Scoring Method",
transf=transf.ipft.hm <http://transf.ipft.hm>,
digits=3) print(pes.da.records, digits=3) print(subganal.lcmbi,
digits=3)print(pes.lcmbi, digits=3) *

And to produce the forest plot I use the following code:

*subganal.lcmbi=rma(yi, vi, data=ies.da.lcmbi, mods=~lcmbi, method="DL")*
*pes.summary=metaprop(*nlameanimal*, *ssizeanimal*, author,
*                     sm="PFT",*
*                     method.tau="DL",*
*                     method.ci <http://method.ci>="NAsm",*
*                     byvar=lcmbi,*
*                     tau.common=TRUE,*
*                     tau.preset=sqrt(subganal.lcmbi$tau2))*
*jpeg(file="forest_cow_cow_lcmbi.jpeg", width=11, height=12, units="in",
*       xlim=c(0,1),*
*       pscale=1,*
*       rightcols=FALSE,*
*       leftcols=c("studlab", "effect", "ci"),*
*       leftlabs=c("Study", "Prevalence", "95% C.I."),*
*       text.random="Prevalence of Lameness in British Dairy Cattle",*
*       bylab = "Lameness Data Source",*
*       print.byvar = TRUE,*
*       xlab="Prevalence of Lameness", smlab="",*
*       weight.study="random", squaresize=0.5, col.square="navy",*
*       col.diamond="maroon", col.diamond.lines="maroon",*
*       pooled.totals=FALSE,*
*       comb.fixed=FALSE,*
*       fs.hetstat=10,*
*       print.tau2=TRUE,*
*       print.Q=TRUE,*
*       print.pval.Q=TRUE,*
*       print.I2=TRUE,*
*       digits=2,*
*       sortvar = precision)*

What could I be doing wrong? I will be extremely grateful for any help or

Wishing all a great afternoon,

João Afonso
DVM, MSc Veterinary Epidemiology
PhD Student
Department of Infection and Global Health
University of Liverpool

	[[alternative HTML version deleted]]

More information about the R-sig-meta-analysis mailing list