[R] Regression and Sub-Groups Analysis in Metafor

Dan Kolubinski kolubind at lsbu.ac.uk
Mon May 30 20:27:52 CEST 2016


I am completing a meta-analysis on the effect of CBT on low self-esteem and
I could use some help regarding the regression feature in metafor.  Based
on the studies that I am using for the analysis, I identified 4 potential
moderators that I want to explore:
- Some of the studies that I am using used RCTs to compare an intervention
with a waitlist and others used the pre-score as the control in a
single-group design.
- Some of the groups took place in one day and others took several weeks.
- There are three discernible interventions being represented
- The initial level of self-esteem varies

Based on the above, I used this command to conduct a meta-analysis using
standarized mean differences:



MetaMod<-rma(m1i=m1, m2i=m2, sd1i=sd1, sd2i=sd2, n1i=n1, n2i=n2,
mods=cbind(dur, rct, int, level),measure = "SMD")



Would this be the best command to use for what I described?  Also, what
could I add to the command so that the forest plot shows a sub-group
analysis using the 'dur' variable as a between-groups distinction?


Also, with respect to the moderators, this is what was delivered:



Test of Moderators (coefficient(s) 2,3,4,5):
QM(df = 4) = 8.7815, p-val = 0.0668

Model Results:

         estimate      se     zval    pval    ci.lb   ci.ub
intrcpt    0.7005  0.6251   1.1207  0.2624  -0.5246  1.9256
dur        0.5364  0.2411   2.2249  0.0261   0.0639  1.0090  *
rct       -0.3714  0.1951  -1.9035  0.0570  -0.7537  0.0110  .
int        0.0730  0.1102   0.6628  0.5075  -0.1430  0.2890
level     -0.2819  0.2139  -1.3180  0.1875  -0.7010  0.1373

---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1



>From this, can I interpret that the variable 'dur' (duration of
intervention) has a significant effect and the variable 'rct' (whether a
study was an RCT or used pre-post scores) was just shy of being
statistically significant?  I mainly ask, because the QM-score has a
p-value of 0.0668, which I thought would mean that none of the moderators
would be significant.  Would I be better off just listing one or two
moderators instead of four?

Much appreciated,
Dan

	[[alternative HTML version deleted]]



More information about the R-help mailing list