Analyse de puissance statistique (GLM et T-test)

Le script commence par créer un jeu de données 'sleep' simulant des scores de rappel ('recall') pour différentes durées ('hours'). Il exécute ensuite PROC GLMPOWER pour déterminer la taille d'échantillon nécessaire dans un modèle linéaire général (ANOVA) pour atteindre des puissances de 0.7, 0.8 et 0.9. Enfin, il utilise PROC POWER pour calculer la taille d'échantillon requise pour un test de comparaison de deux moyennes (t-test) avec une différence de moyenne de 3 et un écart-type de 30.

Analyse des données

Type : CREATION_INTERNE

Les données sont générées via une étape DATA avec des instructions DATALINES (valeurs 30 et 33).

1 Bloc de code

DATA STEP Data

Explication :
Création du jeu de données 'sleep' contenant une variable explicative 'hours' (1 et 2) et une variable réponse 'recall' lue depuis les datalines.

Copié !

1	DATA sleep;
2	DO hours = 1 to 2;
3	INPUT recall;
4	OUTPUT;
5	END;
6	DATALINES;
7	30
8	33
9	;
10	RUN;

2 Bloc de code

PROC GLMPOWER

Explication :
Analyse de puissance pour un modèle linéaire général (ici une ANOVA à un facteur 'hours'). Le script demande le calcul de la taille totale de l'échantillon (ntotal = .) pour des puissances cibles de 0.7, 0.8 et 0.9, avec un écart-type supposé de 30. Un graphique de la puissance en fonction de la taille est généré.

Copié !

1	PROC 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;
9	RUN;

3 Bloc de code

PROC POWER

Explication :
Calcul de puissance pour un test de comparaison de deux moyennes (t-test). Paramètres : différence de moyennes attendue de 3, alpha de 0.05, écart-type de 30. Le code demande la taille d'échantillon requise pour les niveaux de puissance spécifiés.

Copié !

1	PROC POWER;
2	twosamplemeans
3	meandiff=3
4	alpha=0.05
5	ntotal = .
6	power = 0.7 0.8 0.9
7	stddev = 30
8	;
9	RUN;

Ce matériel est fourni "tel quel" par We Are Cas. Il n'y a aucune garantie, expresse ou implicite, quant à la qualité marchande ou à l'adéquation à un usage particulier concernant le matériel ou le code contenu dans les présentes. We Are Cas n'est pas responsable des erreurs dans ce matériel tel qu'il existe maintenant ou existera, et We Are Cas ne fournit pas de support technique pour celui-ci.

Retour à la liste