Ejemplo de inicio para PROC NLMIXED (nlmgs1)

Este programa SAS^© utiliza el procedimiento NLMIXED para modelar una curva de crecimiento logístico. Los datos, de Draper y Smith (1981) sobre el crecimiento de árboles de naranjo, se crean en línea mediante un DATA STEP. El modelo incluye parámetros fijos (b1, b2, b3) y un efecto aleatorio (u1) a nivel del árbol (sujeto), para capturar la variabilidad entre árboles. El modelo es una distribución normal donde la media se define por una función logística y la varianza por s2e.

Análisis de datos

Type : CREATION_INTERNE

Los datos llamados 'tree' se crean directamente en el script SAS utilizando una instrucción DATALINES, incorporando las variables 'tree', 'day' y 'y'.

1 Bloque de código

DATA STEP Data

Explicación :
Este bloque DATA STEP crea un conjunto de datos SAS llamado 'tree'. Utiliza la instrucción INPUT para definir las variables 'tree' (identificador del árbol), 'day' (día) e 'y' (valor medido). Los datos se incluyen directamente en el script después de la instrucción DATALINES, lo que significa que el conjunto de datos se crea a partir de datos literales.

¡Copiado!

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 Bloque de código

PROC NLMIXED

¡Copiado!

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

Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.

Información de copyright : S A S S A M P L E L I B R A R Y

Retour à la liste