Simulating Cluster RCTs
If you can find an existing package that does the kind of model you want, that would certainly be recommended. You could also look at the power analysis section of the GLMM FAQ https://bbolker.github.io/mixedmodels-misc/glmmFAQ.html#power-analysis For scalar/intercept-only random effects, the theta parameters are the among-group standard deviations (scaled by the residual standard deviation, for linear mixed models). For random-slopes models it's (alas) more complicated -- the theta parameters are the elements of the (scaled) Cholesky factor of the variance-covariance matrix -- but you could use the functions in ?vcconv to convert from variance-covariance matrices to Cholesky factors ..
On Fri, Jun 1, 2018 at 4:54 PM, Ulf K?ther <ukoether at uke.de> wrote:
Hi Kevin, I think that a more appropriate way to achieve your goal is to have a look at a library that explicitly deals with the simulation of data. Especially, you should look into the "simstudy" package, which is developed by Keith Goldfield, which is, as I recall it, clearly able to do what you want: The package: https://cran.r-project.org/web/packages/simstudy/index.html Goldfield's blog about generating data using simstudy: https://www.rdatagen.net/ His posts about clustered data: https://www.rdatagen.net/page/clustered/ Good luck, Ulf Am 01.06.2018 um 14:53 schrieb Kevin E. Thorpe:
Hi All. Apologies if this is the wrong list but after searching I have not found what I am looking for. I would like (mainly for teaching) to be able to simulate cluster RCTs with continuous and binary outcomes. It appears to me that simulate.merMod in lme4 may be one way to do this. Unfortunately, I am a bit lost in the help file. For one thing, in the newparams argument, I have no idea what theta is. Also, in the simulations I would like to be able to specify an ICC to control the clustering effect as well as the usual things (e.g. mean response in the control group, treatment effect, number of clusters and cluster sizes). I would be most appreciative for any guidance or examples of this type of simulation, either with simulate.merMod or other approaches. Thank you in advance for your time. Kevin
--
_____________________________________________________________________ Universit?tsklinikum Hamburg-Eppendorf; K?rperschaft des ?ffentlichen Rechts; Gerichtsstand: Hamburg | www.uke.de Vorstandsmitglieder: Prof. Dr. Burkhard G?ke (Vorsitzender), Prof. Dr. Dr. Uwe Koch-Gromus, Joachim Pr?l?, Martina Saurin (komm.) _____________________________________________________________________ SAVE PAPER - THINK BEFORE PRINTING _______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models