Ce script utilise la procédure POWER pour déterminer la taille totale de l'échantillon nécessaire pour atteindre une puissance statistique de 0.9 dans un modèle ANOVA à un facteur, avec des moyennes de groupe et un écart-type spécifiés. Il produit également des graphiques illustrant la relation entre la puissance et la taille de l'échantillon.
Analyse des données
Type : CREATION_INTERNE
Les données nécessaires à l'analyse (moyennes des groupes, écart-type, poids, contrastes) sont définies directement dans les instructions de la procédure PROC POWER. Aucune table externe n'est lue.
1 Bloc de code
PROC POWER
Explication :
Calcul de la taille totale de l'échantillon ('ntotal = .') nécessaire pour obtenir une puissance de 0.9 ('power = 0.9') avec un seuil de signification alpha de 0.025. Les paramètres de l'ANOVA (moyennes, écart-type, poids) sont spécifiés.
Calcul de la taille totale de l'échantillon ('ntotal = .') nécessaire pour obtenir une puissance de 0.9 ('power = 0.9') avec un seuil de signification alpha de 0.025. Les paramètres de l'ANOVA (moyennes, écart-type, poids) sont spécifiés.
Copié !
| 1 | PROC POWER; |
| 2 | onewayanova |
| 3 | groupmeans = 35.6 | 33.7 | 30.2 | 29 28 | 25.9 |
| 4 | stddev = 3.75 |
| 5 | groupweights = (2 1 1 1 1) |
| 6 | alpha = 0.025 |
| 7 | ntotal = . |
| 8 | power = 0.9 |
| 9 | contrast = (4 -1 -1 -1 -1) (0 1 1 -1 -1) |
| 10 | (0 1 -1 0 0) (0 0 0 1 -1); |
| 11 | RUN; |
2 Bloc de code
PROC POWER
Explication :
Activation des graphiques ODS. Exécution de PROC POWER en mode 'plotonly' pour générer un graphique de la taille de l'échantillon en fonction de la puissance (axe X défini sur 'power' variant de 0.5 à 0.95).
Activation des graphiques ODS. Exécution de PROC POWER en mode 'plotonly' pour générer un graphique de la taille de l'échantillon en fonction de la puissance (axe X défini sur 'power' variant de 0.5 à 0.95).
Copié !
| 1 | ods graphics on; |
| 2 | |
| 3 | PROC POWER plotonly; |
| 4 | onewayanova |
| 5 | groupmeans = 35.6 | 33.7 | 30.2 | 29 28 | 25.9 |
| 6 | stddev = 3.75 |
| 7 | groupweights = (2 1 1 1 1) |
| 8 | alpha = 0.025 |
| 9 | ntotal = . |
| 10 | power = 0.9 |
| 11 | contrast = (4 -1 -1 -1 -1) (0 1 1 -1 -1) |
| 12 | (0 1 -1 0 0) (0 0 0 1 -1); |
| 13 | plot x=power min=.5 max=.95; |
| 14 | RUN; |
3 Bloc de code
PROC POWER
Explication :
Calcul de la puissance ('power = .') pour une taille d'échantillon fixe ('ntotal = 24') et génération d'un graphique représentant la puissance en fonction de la taille de l'échantillon (axe X défini sur 'n'). Fermeture des graphiques ODS à la fin.
Calcul de la puissance ('power = .') pour une taille d'échantillon fixe ('ntotal = 24') et génération d'un graphique représentant la puissance en fonction de la taille de l'échantillon (axe X défini sur 'n'). Fermeture des graphiques ODS à la fin.
Copié !
| 1 | PROC POWER plotonly; |
| 2 | onewayanova |
| 3 | groupmeans = 35.6 | 33.7 | 30.2 | 29 28 | 25.9 |
| 4 | stddev = 3.75 |
| 5 | groupweights = (2 1 1 1 1) |
| 6 | alpha = 0.025 |
| 7 | ntotal = 24 |
| 8 | power = . |
| 9 | contrast = (4 -1 -1 -1 -1) (0 1 1 -1 -1) |
| 10 | (0 1 -1 0 0) (0 0 0 1 -1); |
| 11 | plot x=n min=24 max=480; |
| 12 | RUN; |
| 13 | |
| 14 | ods graphics off; |
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.
Informations de Copyright : SAS SAMPLE LIBRARY, PRODUCT: STAT, NAME: pwrex01
Le Conseil de l'Expert
Stéphanie
Spécialiste Machine Learning et IA.
« Une puissance de 0.9 est un standard d'excellence, signifiant que vous n'avez que 10 % de risque de manquer un effet réel (Erreur de Type II). Lors de l'interprétation des graphiques, portez une attention particulière à la pente de la courbe : une pente raide indique qu'une petite réduction de l'échantillon pourrait effondrer votre capacité de détection. »