Mesures répétées multi-populations avec PROC CATMOD

Ce script illustre l'utilisation de la procédure CATMOD pour analyser des réponses binaires (0 ou 1) enregistrées lors de quatre essais successifs sur des sujets provenant de trois groupes différents. Il effectue une analyse de mesures répétées en ajustant d'abord un modèle saturé sur les probabilités marginales, puis un modèle avec des effets imbriqués.

Analyse des données

Type : CREATION_INTERNE

Les données sont générées directement dans le script via l'instruction DATALINES. L'instruction INPUT utilise l'opérateur '@@' pour lire plusieurs observations par ligne.

1 Bloc de code

DATA STEP Data

Explication :
Création du jeu de données 'group'. Les variables a, b, c, d représentent les réponses aux 4 essais. 'Group' est l'identifiant du groupe et 'wt' est la pondération (fréquence). L'opérateur '@@' permet de lire plusieurs observations sur une seule ligne physique.

Copié !

1	DATA group;
2	INPUT a b c d Group wt @@;
3	DATALINES;
4	1 1 1 1 2 2 0 0 0 0 2 2 0 0 1 0 1 2 0 0 1 0 2 2
5	0 0 0 1 1 4 0 0 0 1 2 1 0 0 0 1 3 3 1 0 0 1 2 1
6	0 0 1 1 1 1 0 0 1 1 2 2 0 0 1 1 3 5 0 1 0 0 1 4
7	0 1 0 0 2 1 0 1 0 1 2 1 0 1 0 1 3 2 0 1 1 0 3 1
8	1 0 0 0 1 3 1 0 0 0 2 1 0 1 1 1 2 1 0 1 1 1 3 2
9	1 0 1 0 1 1 1 0 1 1 2 1 1 0 1 1 3 2
10	;

2 Bloc de code

PROC CATMOD

Explication :
Exécution de PROC CATMOD. La première partie ajuste un modèle saturé pour des mesures répétées (définies par 'repeated Trial 4') sur les probabilités marginales. La seconde instruction MODEL (mode interactif) teste une structure spécifique où l'effet des essais est imbriqué dans le groupe 3.

Copié !

1	title 'Multiple-Population Repeated Measures';
2	PROC CATMOD DATA=group;
3	weight wt;
4	response marginals;
5	model abcd=Group _response_ Group_response_
6	/ freq;
7	repeated Trial 4;
8	title2 'Saturated Model';
9	RUN;
10
11	model abc*d=Group _response_(Group=3)
12	/ noprofile noparm design;
13	title2 'Trial Nested within Group 3';
14	QUIT;

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, SAS/STAT User's Guide

Retour à la liste