El script comienza creando un conjunto de datos 'ossification' a partir de datos integrados (datalines), que representan los resultados de un experimento en camadas (litters) con diferentes tratamientos (PHT, TCPO). Se ajusta un primer modelo GEE con una matriz de correlación de trabajo independiente. Luego, los datos se transforman a un formato binario (un registro por sujeto) en el conjunto de datos 'ossification_b'. Se ajusta un segundo modelo GEE a estos datos transformados, utilizando esta vez una matriz de correlación de trabajo intercambiable. Los resultados de ambos análisis se generan en formato HTML.
Análisis de datos
Type : CREATION_INTERNE
Los datos se generan completamente dentro del script mediante una instrucción 'datalines' en el primer paso DATA.
1 Bloque de código
DATA STEP Data
Explicación : Crea el conjunto de datos 'ossification' leyendo los datos directamente del código a través de 'datalines'. Cada línea de datos representa una camada (litter) y el número de éxitos (t) sobre un número de ensayos (m) para diferentes tratamientos (tx). Se crean variables indicadoras (PHT, TCPO) para modelar los efectos del tratamiento.
Explicación : Ajusta un modelo de ecuaciones de estimación generalizadas (GEE) a los datos. El modelo utiliza una distribución binomial con un enlace logit para modelar la proporción de éxitos (t/m) en función de los tratamientos. Se realiza un análisis de medidas repetidas en las camadas (litters) asumiendo una estructura de correlación de trabajo independiente (type=ind).
¡Copiado!
ods html;
title "*** Ossification Data -- GEE using GENMOD with Independent Working Correlation Matrix***";
proc genmod data=ossification;
class litters tcpo pht / param=ref;
model t/m = tcpo pht tcpo*pht / dist=binomial link=logit;
repeated subject=litters / type=ind;
run;
ods html close;
1
ods html;
2
title "*** Ossification Data -- GEE using GENMOD with Independent Working Correlation Matrix***";
3
PROC GENMODDATA=ossification;
4
class litters tcpo pht / param=ref;
5
model t/m = tcpo pht tcpo*pht / dist=binomial link=logit;
6
repeated subject=litters / type=ind;
7
RUN;
8
ods html close;
3 Bloque de código
DATA STEP Data
Explicación : Transforma el conjunto de datos 'ossification' agregado a un formato binario. Para cada camada, crea 't' observaciones con una variable de respuesta 'y' igual a 1 (éxito) y 'm-t' observaciones con 'y' igual a 0 (fracaso). Este formato es necesario para ciertos tipos de análisis binomiales.
¡Copiado!
data ossification_b;
set ossification;
do i=1 to t;
y = 1;
output;
end;
do i=1 to m-t;
y = 0;
output;
end;
run;
1
DATA ossification_b;
2
SET ossification;
3
DO i=1 to t;
4
y = 1;
5
OUTPUT;
6
END;
7
DO i=1 to m-t;
8
y = 0;
9
OUTPUT;
10
END;
11
RUN;
4 Bloque de código
PROC GENMOD
Explicación : Ajusta un segundo modelo GEE utilizando el conjunto de datos binario 'ossification_b'. El modelo es similar al primero, pero especifica una estructura de correlación de trabajo intercambiable (type=exch), asumiendo una correlación constante entre todas las observaciones dentro de una misma camada.
¡Copiado!
ods html;
title "*** Ossification Data -- GEE using GENMOD with Exchangeable Working Correlation Matrix***";
proc genmod data=ossification_b descending;
class litters tcpo pht / param=ref;
model y = tcpo pht tcpo*pht / dist=binomial link=logit;
repeated subject=litters / type=exch;
run;
ods html close;
1
ods html;
2
title "*** Ossification Data -- GEE using GENMOD with Exchangeable Working Correlation Matrix***";
3
PROC GENMODDATA=ossification_b descending;
4
class litters tcpo pht / param=ref;
5
model y = tcpo pht tcpo*pht / dist=binomial link=logit;
6
repeated subject=litters / type=exch;
7
RUN;
8
ods html close;
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.
SAS y todos los demás nombres de productos o servicios de SAS Institute Inc. son marcas registradas o marcas comerciales de SAS Institute Inc. en los EE. UU. y otros países. ® indica registro en los EE. UU. WeAreCAS es un sitio comunitario independiente y no está afiliado a SAS Institute Inc.
Este sitio utiliza cookies técnicas y analíticas para mejorar su experiencia.
Saber más.
