- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Here is the code for the model I will be running:
proc glimmix data=ecig method=quad oddsratio;
class ID Year(ref='2018');
model vape(Event='1')=Year /s dist=binary link=logit ddfm=none
chisq;
random intercept/subject=ID;
output out=glmxout predicted(blup ilink)=predprob;
run;
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
You can search on the web for papers on this field like Estimating power in (generalized) linear mixed models: An open introduction and tutorial in R | Beha.... R code is provided in the paper.
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Another paper (Calculating power in analysis of variance. (apa.org)) found on the book Statistical Power Analysis for the Behavioral Sciences | Jacob Cohen | (taylorfrancis.com) also discusses the topic of power calculation in mixed effect models. Of closer revelance to your question is another paper cited in this book: Optimum Sample Size and Number of Levels in a One-Way Random-Effects Analysis of Variance: The Journ.... A table of optimum sample size or number of levels for the random effects model is included in the latter paper. I hope these information helps.