Es fácil calcular la potencia para pruebas t pareadas y ANOVA, pero ¿cómo se calcula la potencia para diseños más complicados, como diseños de parcelas divididas o diseños factoriales?
Respuestas
¿Demasiados anuncios?Una vez superados los casos sencillos, como las pruebas t, prefiero utilizar simulaciones. Cuando se hace una simulación, se controlan todas las suposiciones que se hacen y se pueden simular situaciones que pueden no estar en las bonitas rutinas enlatadas.
He aquí una respuesta con un ejemplo de simulación:
Primero convierte tu anova en una regresión. Entonces,
$$ N = \frac{2.84^2}{p(1-p)}\frac{\sigma^2}{MDE^2} $$
N es el tamaño de la muestra, p es la proporción que recibe el tratamiento, $\sigma$ es la desviación estándar de los residuos, y MDE=el efecto mínimo detectable que se busca. 2,84 procede de una confianza alfa del 95% y una potencia del 80% (cuando se tienen muchos grados de libertad. Con menos grados de libertad, tendrá que utilizar la distribución T. Mire la referencia para más detalles).
Cuantos más factores explicativos tengas en tu experimento más $\sigma$ se encogerá. Saber cuánto se encogerá es complicado.
Llega un momento en que las fórmulas analíticas se colapsan y es mejor simular todo el conjunto de datos varias veces mediante algún tipo de método monte Carlo y ajustar el modelo deseado a cada uno de los conjuntos de datos plausibles. La potencia es la proporción de veces que se obtiene el resultado deseado. Al igual que con los métodos analíticos, la calidad depende de los supuestos.
La fuente es de memoria, aproximadamente la siguiente este .