Análisis de Varianza Desequilibrada con EFFECTPLOT (ICOMEP2)

Nivel de dificultad

Principiante

Publicado el : 14/10/2024

Este script ilustra cómo analizar un diseño experimental de dos factores (fármaco y enfermedad) con interacciones, sobre datos desequilibrados. Utiliza el procedimiento GENMOD para ajustar el modelo y la instrucción EFFECTPLOT para generar varios gráficos: curvas de efectos, interacciones con segmentación (slice), diagramas de caja (boxplots) y gráficos de mosaico.

Análisis de datos

Type : CREATION_INTERNE

Los datos se generan dinámicamente en el paso DATA 'a' a través de la instrucción DATALINES. Simulan un estudio clínico con las variables 'drug', 'disease' y la respuesta 'y'.

1 Bloque de código

DATA STEP Data

Explicación :
Creación de la tabla de datos 'a'. El paso DATA lee los identificadores de medicamento y enfermedad, luego itera para leer hasta 6 medidas 'y' por línea de entrada.

¡Copiado!

1	DATA a;
2	INPUT drug disease @;
3	DO i=1 to 6;
4	INPUT y @;
5	OUTPUT;
6	END;
7	DATALINES;
8	1 1 42 44 36 13 19 22
9	1 2 33 . 26 . 33 21
10	1 3 31 -3 . 25 25 24
11	2 1 28 . 23 34 42 13
12	2 2 . 34 33 31 . 36
13	2 3 3 26 28 32 4 16
14	3 1 . . 1 29 . 19
15	3 2 . 11 9 7 1 -6
16	3 3 21 1 . 9 3 .
17	4 1 24 . 9 22 -2 15
18	4 2 27 12 12 -5 16 15
19	4 3 22 7 25 5 12 .
20	;

2 Bloque de código

PROC GENMOD

Explicación :
Activación de gráficos ODS. Ejecución de PROC GENMOD para modelar la respuesta 'y' (distribución normal por defecto). 'effectplot / obs' muestra el gráfico de efectos con las observaciones. 'effectplot interaction(sliceby=disease)' muestra las interacciones segmentando por la variable enfermedad, con los límites de confianza (clm).

¡Copiado!

1	ods graphics on;
2	PROC GENMOD DATA=a;
3	class drug disease;
4	model y=disease drug disease*drug / d=n;
5	effectplot / obs;
6	effectplot interaction(sliceby=disease) / clm;
7	RUN;

3 Bloque de código

PROC GENMOD

Explicación :
Segundo análisis con visualizaciones complementarias: diagramas de caja (box), interacciones cruzadas con observaciones, gráficos de mosaico (mosaic), e interacciones facetadas por enfermedad (plotby).

¡Copiado!

1	PROC GENMOD DATA=a;
2	class drug disease;
3	model y=drug disease drug*disease / d=n;
4	effectplot box;
5	effectplot interaction(x=drug*disease) / obs;
6	effectplot mosaic;
7	effectplot interaction(plotby=disease);
8	RUN;
9	ods graphics off;

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.

Documentación relacionada

Aucune documentation spécifique pour cette catégorie.