Modelo de Coeficientes Aleatorios con PROC MIXED

El script comienza creando un conjunto de datos llamado 'rc' a través de un paso DATA con datos en línea (datalines). Estos datos simulan un estudio de estabilidad farmacéutica con lotes (Batch) y mediciones (Y) tomadas en diferentes meses (Month). Luego, se ejecutan dos procedimientos MIXED. El primero ajusta un modelo de regresión lineal con coeficientes aleatorios (intersección y pendiente 'Month') para cada lote. El segundo explora una variación del modelo incluyendo una variable de clase adicional 'Monthc'.

Análisis de datos

Type : CREATION_INTERNE

Los datos se generan directamente en el script utilizando un paso DATA y una instrucción DATALINES. El conjunto de datos 'rc' se crea en memoria para la sesión SAS.

1 Bloque de código

DATA STEP Data

Explicación :
Este bloque DATA crea la tabla 'rc'. Lee las variables 'Batch' y 'Month', luego utiliza un bucle 'do' para leer hasta 6 valores de 'Y' para cada registro inicial. Para cada valor 'Y' leído, se genera una nueva observación. La variable 'Monthc' se crea como una copia de 'Month'. Los datos se proporcionan en línea a través de 'datalines'.

¡Copiado!

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

PROC MIXED

Explicación :
Este procedimiento ajusta un modelo mixto sobre los datos 'rc'. 'Batch' se define como una variable de clasificación. El modelo ('model') especifica 'Y' como variable dependiente y 'Month' como efecto fijo. La instrucción 'random' define una intersección ('Int') y una pendiente ('Month') aleatorias para cada nivel de 'Batch', con una estructura de covarianza no estructurada ('type=un'). La opción 's' solicita la visualización de las soluciones para los efectos fijos y aleatorios.

¡Copiado!

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

PROC MIXED

Explicación :
Se realiza un segundo análisis PROC MIXED, añadiendo 'Monthc' a la lista de variables de clasificación. El modelo fijo permanece igual. La instrucción 'random' se modifica para incluir 'Monthc' como efecto aleatorio además de la intersección y 'Month', siempre agrupados por 'Batch'. Esto permite explorar una estructura de modelo aleatoria diferente.

¡Copiado!

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

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