Publié le :
Statistique CREATION_INTERNE

Exemple 4 pour PROC GLM - Analyse de Covariance

Ce code est également disponible en : Deutsch English Español
En attente de validation
Le script crée d'abord un jeu de données 'DrugTest' contenant l'effet de trois médicaments (A, D, F). Ensuite, il utilise PROC GLM pour modéliser l'effet du médicament et de la mesure avant traitement sur la mesure après traitement. Les moyennes des moindres carrés (lsmeans) sont calculées pour comparer les effets des médicaments et sont sauvegardées dans la table 'adjmeans', qui est ensuite affichée. Finalement, le script utilise ODS Graphics et une seconde passe de PROC GLM pour visualiser le modèle d'analyse de covariance ajusté, y compris un graphique des moyennes avec des limites de confiance.
Analyse des données

Type : CREATION_INTERNE


Les données sont créées directement dans le script à l'aide d'une instruction DATA et de 'datalines'.

1 Bloc de code
DATA STEP Data
Explication :
Ce bloc crée le jeu de données 'DrugTest' en lisant les données intégrées via 'datalines'. Il définit trois variables : 'Drug' (caractère), 'PreTreatment' (numérique) et 'PostTreatment' (numérique).
Copié !
1DATA DrugTest;
2 INPUT Drug $ PreTreatment PostTreatment;
3 DATALINES;
4A 11 6 A 8 0 A 5 2 A 14 8 A 19 11
5A 6 4 A 10 13 A 6 1 A 11 8 A 3 0
6D 6 0 D 6 2 D 7 3 D 8 1 D 18 18
7D 8 4 D 19 14 D 8 9 D 5 1 D 15 9
8F 16 13 F 13 10 F 11 18 F 9 5 F 21 23
9F 16 12 F 12 5 F 12 16 F 7 1 F 12 20
10;
11RUN;
2 Bloc de code
PROC GLM Data
Explication :
Cette procédure de modèle linéaire généralisé (GLM) effectue une analyse de covariance. 'PostTreatment' est la variable dépendante, modélisée par l'effet de 'Drug' (variable de classe) et 'PreTreatment' (covariable). L'option 'solution' demande l'affichage des estimations des paramètres. L'instruction 'lsmeans' calcule les moyennes ajustées pour chaque niveau de 'Drug', avec les erreurs standard, les p-values pour les différences, la matrice de covariance, et sauvegarde les résultats dans un nouveau dataset nommé 'adjmeans'.
Copié !
1PROC GLM DATA=DrugTest;
2 class Drug;
3 model PostTreatment = Drug PreTreatment / solution;
4 lsmeans Drug / stderr pdiff cov out=adjmeans;
5RUN;
3 Bloc de code
PROC PRINT
Explication :
Ce bloc de code affiche le contenu du jeu de données 'adjmeans', qui a été créé par l'instruction 'lsmeans' dans l'étape PROC GLM précédente.
Copié !
1PROC PRINT DATA=adjmeans;
2RUN;
4 Bloc de code
PROC GLM
Explication :
Ce bloc active la sortie graphique ODS pour visualiser les résultats. La procédure GLM est de nouveau appelée, cette fois avec l'option 'plot=meanplot(cl)' pour générer un graphique des moyennes ajustées avec des limites de confiance pour la variable 'Drug'. L'instruction 'lsmeans' avec l'option 'pdiff' calcule les différences par paires entre les moyennes des moindres carrés. ODS Graphics est ensuite désactivé.
Copié !
1ods graphics on;
2 
3PROC GLM DATA=DrugTest plot=meanplot(cl);
4 class Drug;
5 model PostTreatment = Drug PreTreatment;
6 lsmeans Drug / pdiff;
7RUN;
8 
9ods 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 : S A S S A M P L E L I B R A R Y