Skip to content
Prev 334091 / 398513 Next

How do I print predicted effect sizes in forest plot?

Viechtbauer Wolfgang (STAT) 
One more thing ... You used the command:

iat_result = rma(yi=Mean, vi=Variance_rounded, ni=N, sei=Std_error, slab=Study_Name, subset=(Country == "AUT"), data=cma_iat, method="HS")

This probably does not do what you want it to do. First of all, if you specify vi, there is no need to specify sei (or vice-versa). One is sufficient. But more crucially, I assume 'Mean' is what it says it is - a mean of a certain variable X. And I assume that 'Variance_rounded' is the variance of said variable X. But vi is used to specify the *sampling variance* of yi (or sei is used to specify the standard error), which, for a mean, is the variance divided by N (and the standard error is the SD divided by the square root of N):

http://en.wikipedia.org/wiki/Standard_error_of_the_mean#Standard_error_of_the_mean

So, my hunch is that you are not supplying the right information to the rma() function.

Best,
Wolfgang

Thread (5 messages)

Alma Wilflinger How do I print predicted effect sizes in forest plot? Dec 5 Viechtbauer Wolfgang (STAT) How do I print predicted effect sizes in forest plot? Dec 6 Viechtbauer Wolfgang (STAT) How do I print predicted effect sizes in forest plot? Dec 6 Michael Dewey How do I print predicted effect sizes in forest plot? Dec 6 Alma Wilflinger How do I print predicted effect sizes in forest plot? Dec 9