Exemple de démarrage pour PROC NLMIXED (nlmgs1)

Ce programme SAS^© utilise la procédure NLMIXED pour modéliser une courbe de croissance logistique. Les données, issues de Draper et Smith (1981) sur la croissance des arbres d'orange, sont créées en ligne via un DATA STEP. Le modèle inclut des paramètres fixes (b1, b2, b3) et un effet aléatoire (u1) au niveau de l'arbre (sujet), pour capturer la variabilité inter-arbres. Le modèle est une distribution normale où la moyenne est définie par une fonction logistique et la variance par s2e.

Analyse des données

Type : CREATION_INTERNE

Les données nommées 'tree' sont créées directement dans le script SAS en utilisant une instruction DATALINES, incorporant les variables 'tree', 'day' et 'y'.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA STEP crée un jeu de données SAS nommé 'tree'. Il utilise l'instruction INPUT pour définir les variables 'tree' (identifiant de l'arbre), 'day' (jour) et 'y' (valeur mesurée). Les données sont directement incluses dans le script après l'instruction DATALINES, ce qui signifie que le jeu de données est créé à partir de données littérales.

Copié !

1	DATA tree;
2	INPUT tree day y;
3	DATALINES;
4	1 118 30
5	1 484 58
6	1 664 87
7	1 1004 115
8	1 1231 120
9	1 1372 142
10	1 1582 145
11	2 118 33
12	2 484 69
13	2 664 111
14	2 1004 156
15	2 1231 172
16	2 1372 203
17	2 1582 203
18	3 118 30
19	3 484 51
20	3 664 75
21	3 1004 108
22	3 1231 115
23	3 1372 139
24	3 1582 140
25	4 118 32
26	4 484 62
27	4 664 112
28	4 1004 167
29	4 1231 179
30	4 1372 209
31	4 1582 214
32	5 118 30
33	5 484 49
34	5 664 81
35	5 1004 125
36	5 1231 142
37	5 1372 174
38	5 1582 177
39	;
40

2 Bloc de code

PROC NLMIXED

Explication :
Ce bloc exécute la procédure NLMIXED pour ajuster un modèle non linéaire à effets mixtes au jeu de données 'tree'.
- L'instruction PARMS initialise les paramètres du modèle : b1, b2, b3 (paramètres fixes de la courbe logistique), s2u (variance des effets aléatoires) et s2e (variance de l'erreur résiduelle).
- Les équations 'num', 'ex', 'den' définissent la forme logistique de la moyenne du modèle.
- L'instruction MODEL spécifie que la variable de réponse 'y' suit une distribution normale avec une moyenne 'num/den' et une variance 's2e'.
- L'instruction RANDOM déclare 'u1' comme un effet aléatoire suivant une distribution normale de moyenne 0 et de variance 's2u', avec 'tree' comme variable sujet, indiquant que les effets aléatoires sont spécifiques à chaque arbre.

Copié !

1	PROC NLMIXED DATA=tree;
2	parms b1=190 b2=700 b3=350 s2u=1000 s2e=60;
3	num = b1+u1;
4	ex = exp(-(day-b2)/b3);
5	den = 1 + ex;
6	model y ~ normal(num/den,s2e);
7	random u1 ~ normal(0,s2u) subject=tree;
8	RUN;
9

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