References about powerlmm HOT 1 CLOSED

Hi,

As far as I know there's no article that describe exactly the same calculations. However, several articles use the general formulas of the generalized least squares estimator to get variance of the treatment effect. My calculations are described in the package's technical vignette:.

As you note, there exists several different definitions of both the variance ratio and the standardized effects. Tho, they make no difference for the actual calculations, but they do lead to different raw values. Currently, I use calculations that are more common in my field (clinical psychology). My plan is to add other definitions as well. I the mean time it should be possible to simply enter the corresponding variance components, and avoid using my standardizations. Though, this is not the best solution as you'd have to manually translate the treatment effect to the same units as my definition. The next release will, at least, include the possibility to use unstandardized value for the time:treatment-interaction, which will make custom definitions more easy. And I should probably make this more clear in the documentation.

I suspect the differences to PASS is either caused by the models not being exactly the same, or differences in the degrees of freedom / distribution used when doing the power analysis. Many other implementations simply use normal approximations, and avoid the issue of correct dfs.

Thanks for the feedback!

from powerlmm.