Modèle à Coefficients Aléatoires avec PROC MIXED

Le script commence par créer un jeu de données nommé 'rc' via une étape DATA avec des données en ligne (datalines). Ces données simulent une étude de stabilité pharmaceutique avec des lots (Batch) et des mesures (Y) prises à différents mois (Month). Ensuite, deux procédures MIXED sont exécutées. La première ajuste un modèle de régression linéaire avec des coefficients aléatoires (intercept et pente 'Month') pour chaque lot. La seconde explore une variation du modèle en incluant une variable de classe supplémentaire 'Monthc'.

Analyse des données

Type : CREATION_INTERNE

Les données sont générées directement dans le script à l'aide d'une étape DATA et d'une instruction DATALINES. Le jeu de données 'rc' est créé en mémoire pour la session SAS.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA crée la table 'rc'. Il lit les variables 'Batch' et 'Month', puis utilise une boucle 'do' pour lire jusqu'à 6 valeurs de 'Y' pour chaque enregistrement initial. Pour chaque valeur 'Y' lue, une nouvelle observation est générée. La variable 'Monthc' est créée comme une copie de 'Month'. Les données sont fournies en ligne via 'datalines'.

Copié !

1	DATA rc;
2	INPUT Batch Month @;
3	Monthc = Month;
4	DO i = 1 to 6;
5	INPUT Y @;
6	OUTPUT;
7	END;
8	DATALINES;
9	1 0 101.2 103.3 103.3 102.1 104.4 102.4
10	1 1 98.8 99.4 99.7 99.5 . .
11	1 3 98.4 99.0 97.3 99.8 . .
12	1 6 101.5 100.2 101.7 102.7 . .
13	1 9 96.3 97.2 97.2 96.3 . .
14	1 12 97.3 97.9 96.8 97.7 97.7 96.7
15	2 0 102.6 102.7 102.4 102.1 102.9 102.6
16	2 1 99.1 99.0 99.9 100.6 . .
17	2 3 105.7 103.3 103.4 104.0 . .
18	2 6 101.3 101.5 100.9 101.4 . .
19	2 9 94.1 96.5 97.2 95.6 . .
20	2 12 93.1 92.8 95.4 92.2 92.2 93.0
21	3 0 105.1 103.9 106.1 104.1 103.7 104.6
22	3 1 102.2 102.0 100.8 99.8 . .
23	3 3 101.2 101.8 100.8 102.6 . .
24	3 6 101.1 102.0 100.1 100.2 . .
25	3 9 100.9 99.5 102.2 100.8 . .
26	3 12 97.8 98.3 96.9 98.4 96.9 96.5
27	;
28

2 Bloc de code

PROC MIXED

Explication :
Cette procédure ajuste un modèle mixte sur les données 'rc'. 'Batch' est défini comme une variable de classification. Le modèle ('model') spécifie 'Y' comme variable dépendante et 'Month' comme effet fixe. L'instruction 'random' définit un intercept ('Int') et une pente ('Month') aléatoires pour chaque niveau de 'Batch', avec une structure de covariance non structurée ('type=un'). L'option 's' demande l'affichage des solutions pour les effets fixes et aléatoires.

Copié !

1	PROC MIXED DATA=rc;
2	class Batch;
3	model Y = Month / s;
4	random Int Month / type=un sub=Batch s;
5	RUN;

3 Bloc de code

PROC MIXED

Explication :
Une seconde analyse PROC MIXED est effectuée, ajoutant 'Monthc' à la liste des variables de classification. Le modèle fixe reste le même. L'instruction 'random' est modifiée pour inclure 'Monthc' comme effet aléatoire en plus de l'intercept et de 'Month', toujours groupé par 'Batch'. Cela permet d'explorer une structure de modèle aléatoire différente.

Copié !

1	PROC MIXED DATA=rc;
2	class Batch Monthc;
3	model Y = Month / s;
4	random Int Month Monthc / sub=Batch s;
5	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.

Informations de Copyright : S A S S A M P L E L I B R A R Y

Retour à la liste