[R] metafor estimates using mods and subset do not match
Morris, Kathryn
morrisk10 at xavier.edu
Fri Jul 29 15:44:42 CEST 2016
I am running a meta-analysis using metafor and getting what seem to be conflicting results.
#analysis with species moderator
cropMeta.species<-rma(cropyi, cropvi, data=dataCropMeta, mods=~dataCropMeta$species - 1, method="HE")
The output for the analysis using species as a moderator shows estimates for the seven species in my data set, and makes sense.
#subset analysis with B oleraceae
cropMeta.species.b.oleracea<-rma(cropyi, cropvi, data=dataCropMeta, subset=(species=="B. oleracea"), method="HE")
Then when I start to plot the results and use separate analyses for each species in order to get the coef and variances for plotting the subgroups I get different estimates. The estimates for most subgroups are almost identical in both analyses, but one of them is 5.922 in the analysis with species as a moderator, and 7.0139 in the subset analysis for only that species.
Why are the estimates different? If they should be different, then which analysis is producing the correct estimates?
I'm happy to provide the data and R scripts if that would help.
many thanks,
Kathryn
---------------------------------
Dr. Kathryn Morris
Assistant Professor, Biology
Xavier University
3800 Victory Parkway
Cincinnati, OH
45207
Office: (513) 745-3554
More information about the R-help
mailing list