Thank you very much Wolfgang and Michael for your help and information. Now I know how to get the predicted values.
For information, my "mean" is the effect size for the one-sample-case. Obviously I used misleading labels, but I calculated the same with CMA and getting the same results, so it should be okay. Best wishes, Alma On Friday, December 6, 2013 3:10 PM, Michael Dewey <i...@aghmed.fsnet.co.uk> wrote: At 21:05 05/12/2013, Alma Wilflinger wrote: >Hi, > >I am struggling a bit with creating a forest plot containing the >predicted effect size. As seen in other studies these effect sizes >are shown per study usually as a light grey diamond - which is what >I want to achieve. > >The calls I use are: >iat_result = rma(yi=Mean, vi=Variance_rounded, ni=N, sei=Std_error, >slab=Study_Name, subset=(Country == "AUT"), data=cma_iat, method="HS") Alma You do not need to specify both vi and sei as one is sufficient and you do not need ni as well. I realise that is not the question you asked (which Wolfgang has already answered). >summary.rma(iat_result) > > >#not sure how to use it or if needed >#predict(iat_result) > >forest(iat_result) > > >At the end I am getting the forest plot as is without the predicted values. > >I am not sure if I need the predict function and how to use it? - >the predict function deliveres the same values as already computed >in the rma object. > > >I checked the manual for package metafor but was not able to find >out how to print the predicted values per study. > > >kind regards, Alma > [[alternative HTML version deleted]] Michael Dewey i...@aghmed.fsnet.co.uk http://www.aghmed.fsnet.co.uk/home.html [[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.