Publicado el :
Estadística CREATION_INTERNE

Análisis de potencia estadística (GLM y T-test)

Este código también está disponible en: Deutsch English Français
En espera de validación
El script comienza creando un conjunto de datos 'sleep' que simula puntuaciones de recuerdo ('recall') para diferentes duraciones ('hours'). Luego ejecuta PROC GLMPOWER para determinar el tamaño de muestra necesario en un modelo lineal general (ANOVA) para alcanzar potencias de 0.7, 0.8 y 0.9. Finalmente, utiliza PROC POWER para calcular el tamaño de muestra requerido para una prueba de comparación de dos medias (t-test) con una diferencia de media de 3 y una desviación estándar de 30.
Análisis de datos

Type : CREATION_INTERNE


Los datos se generan mediante un paso DATA con instrucciones DATALINES (valores 30 y 33).

1 Bloque de código
DATA STEP Data
Explicación :
Creación del conjunto de datos 'sleep' que contiene una variable explicativa 'hours' (1 y 2) y una variable de respuesta 'recall' leída desde las datalines.
¡Copiado!
1DATA sleep;
2 DO hours = 1 to 2;
3 INPUT recall;
4 OUTPUT;
5 END;
6 DATALINES;
730
833
9;
10RUN;
2 Bloque de código
PROC GLMPOWER
Explicación :
Análisis de potencia para un modelo lineal general (aquí un ANOVA de un factor 'hours'). El script solicita el cálculo del tamaño total de la muestra (ntotal = .) para potencias objetivo de 0.7, 0.8 y 0.9, con una desviación estándar supuesta de 30. Se genera un gráfico de la potencia en función del tamaño.
¡Copiado!
1PROC GLMPOWER DATA=sleep;
2 class hours;
3 model recall = hours;
4 power
5 stddev = 30
6 ntotal = .
7 power = 0.7 0.8 0.9;
8 plot x=power min=0.7 max=0.95;
9RUN;
3 Bloque de código
PROC POWER
Explicación :
Cálculo de potencia para una prueba de comparación de dos medias (t-test). Parámetros: diferencia de medias esperada de 3, alfa de 0.05, desviación estándar de 30. El código solicita el tamaño de muestra requerido para los niveles de potencia especificados.
¡Copiado!
1PROC POWER;
2 twosamplemeans
3 meandiff=3
4 alpha=0.05
5 ntotal = .
6 power = 0.7 0.8 0.9
7 stddev = 30
8 ;
9RUN;
Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.