Me pregunto si hay algún método para calcular el tamaño de la muestra en modelos mixtos. Estoy usando lmer
R para ajustar los modelos (tengo pendientes e intercepciones aleatorias).
r
mixed-model
lme4-nlme
power-analysis
Nikita Kuznetsov
fuente
fuente
Respuestas:
El
longpower
paquete implementa los cálculos del tamaño de la muestra en Liu y Liang (1997) y Diggle et al (2002). La documentación tiene un código de ejemplo. Aquí hay uno, usando lalmmpower()
función:Compruebe también
liu.liang.linear.power()
que " realiza el cálculo del tamaño de muestra para un modelo mixto lineal"Liu, G. y Liang, KY (1997). Cálculos de tamaño de muestra para estudios con observaciones correlacionadas. Biometrics, 53 (3), 937-47.
Diggle PJ, Heagerty PJ, Liang K, Zeger SL. Análisis de datos longitudinales. Segunda edicion. Oxford Ciencia Estadística Serires. 2002
Editar: Otra forma es "corregir" el efecto de la agrupación. En un modelo lineal ordinario, cada observación es independiente, pero en presencia de agrupamiento, las observaciones no son independientes, lo que puede considerarse como tener menos observaciones independientes: el tamaño efectivo de la muestra es menor. Esta pérdida de efectividad se conoce como el efecto de diseño :
fuente
Para cualquier cosa más allá de las simples pruebas de 2 muestras, prefiero usar la simulación para el tamaño de la muestra o los estudios de potencia. Con las rutinas preempaquetadas, a veces puede ver grandes diferencias entre los resultados de los programas en función de los supuestos que están haciendo (y es posible que no pueda descubrir cuáles son esos supuestos, y mucho menos si son razonables para su estudio). Con la simulación, controlas todos los supuestos.
Aquí hay un enlace a un ejemplo:
https://stat.ethz.ch/pipermail/r-sig-mixed-models/2009q1/001790.html
fuente