Hello,
I'm using the rmeta package to perform a meta analysis using summary statistics
rather than raw data, and would like to analyze the effects in three different
subgroups of my data. Furthermore, I'd like to plot this on one forest plot,
with corresponding summary weighted averages of the effects displayed beneath
each subgroup.
I am able to generate the subgroup analyses by simply performing 3 separate
meta-analyses with the desired subset of data. However, I can't manage to plot
everything on the same forest plot.
My data and code appear as follows:
data(meta)
subgroup study beta se
subgroupA site1 -0.35 0.12
subgroupA site2 -0.34 0.1
subgroupB site1 -0.28 0.06
subgroupB site2 -0.29 0.07
subgroupC site1 0.34 0.03
subgroupC site2 0.36 0.04
#meta-analysis
A<-meta[c(1,2),]
B<-meta[c(3,4),]
C<-meta[c(5,6),]
Ameta<-meta.summaries(A$beta,A$se,method="fixed",logscale=TRUE,names=A$study)
Bmeta<-meta.summaries(B$beta,A$se,method="fixed",logscale=TRUE,names=B$study)
Cmeta<-meta.summaries(C$beta,A$se,method="fixed",logscale=TRUE,names=C$study)
list(summary(Ameta),summary(Bmeta),summary(Cmeta))
I get my output, but I'm not able to plot the three subgroups (A,B,C) on the
same forest plot using any of the functions I see in rmeta. Any ideas?
Thanks,
Steve Lubitz
Cardiovascular Research Fellow, Brigham and Women's and Massachusetts General
Hospital
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
