Analyse de Variance Déséquilibrée avec EFFECTPLOT (ICOMEP2)

Niveau de difficulté

Débutant

Publié le : 06/06/2023

Ce script illustre comment analyser un plan d'expérience à deux facteurs (médicament et maladie) avec des interactions, sur des données déséquilibrées. Il utilise la procédure GENMOD pour ajuster le modèle et l'instruction EFFECTPLOT pour générer divers graphiques : courbes d'effets, interactions avec découpage (slice), boîtes à moustaches (boxplots) et graphiques en mosaïque.

Analyse des données

Type : CREATION_INTERNE

Les données sont générées dynamiquement dans l'étape DATA 'a' via l'instruction DATALINES. Elles simulent une étude clinique avec les variables 'drug', 'disease' et la réponse 'y'.

1 Bloc de code

DATA STEP Data

Explication :
Création de la table de données 'a'. L'étape DATA lit les identifiants de médicament et de maladie, puis boucle pour lire jusqu'à 6 mesures 'y' par ligne d'entrée.

Copié !

1	DATA a;
2	INPUT drug disease @;
3	DO i=1 to 6;
4	INPUT y @;
5	OUTPUT;
6	END;
7	DATALINES;
8	1 1 42 44 36 13 19 22
9	1 2 33 . 26 . 33 21
10	1 3 31 -3 . 25 25 24
11	2 1 28 . 23 34 42 13
12	2 2 . 34 33 31 . 36
13	2 3 3 26 28 32 4 16
14	3 1 . . 1 29 . 19
15	3 2 . 11 9 7 1 -6
16	3 3 21 1 . 9 3 .
17	4 1 24 . 9 22 -2 15
18	4 2 27 12 12 -5 16 15
19	4 3 22 7 25 5 12 .
20	;

2 Bloc de code

PROC GENMOD

Explication :
Activation des graphiques ODS. Exécution de PROC GENMOD pour modéliser la réponse 'y' (distribution normale par défaut). 'effectplot / obs' affiche le graphique des effets avec les observations. 'effectplot interaction(sliceby=disease)' affiche les interactions en découpant par la variable maladie, avec les limites de confiance (clm).

Copié !

1	ods graphics on;
2	PROC GENMOD DATA=a;
3	class drug disease;
4	model y=disease drug disease*drug / d=n;
5	effectplot / obs;
6	effectplot interaction(sliceby=disease) / clm;
7	RUN;

3 Bloc de code

PROC GENMOD

Explication :
Seconde analyse avec des visualisations complémentaires : diagrammes en boîte (box), interactions croisées avec observations, graphiques en mosaïque (mosaic), et interactions facettées par maladie (plotby).

Copié !

1	PROC GENMOD DATA=a;
2	class drug disease;
3	model y=drug disease drug*disease / d=n;
4	effectplot box;
5	effectplot interaction(x=drug*disease) / obs;
6	effectplot mosaic;
7	effectplot interaction(plotby=disease);
8	RUN;
9	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.

Documentation liée

Scoring Distribué : Déployez vos Modèles Analytiques en RAM

Cet exemple démontre comment télécharger un 'analytic store' (astore) depuis le système de fichie...