Este script crea un conjunto de datos 'six' que contiene información sobre síntomas respiratorios (sibilancias) en dos ciudades (Kingston y Portage). Luego utiliza el procedimiento GENMOD para ajustar un modelo de regresión logística (distribución binomial) para predecir la presencia de sibilancias respiratorias en función de la ciudad, la edad y el estado de tabaquismo. La instrucción REPEATED se utiliza para especificar una estructura de correlación intercambiable (type=exch) para manejar las mediciones repetidas en el mismo sujeto (variable 'case').
Análisis de datos
Type : CREATION_INTERNE
Los datos se generan manualmente en el script a través del paso DATA 'six' y la instrucción DATALINES. No se requiere ninguna fuente externa.
1 Bloque de código
DATA STEP Data
Explicación : Crea el conjunto de datos 'six' leyendo datos brutos incluidos en el código. Se utiliza un bucle DO para leer varias observaciones (4 mediciones repetidas por sujeto) de una sola línea de datos fuente. Nota: El código original contiene artefactos ' @code_sas_json...' que parecen ser errores de copiar y pegar que reemplazan el signo estándar de retención de línea ' @code_sas/16.4'.sas.
¡Copiado!
data six;
input case city$ @;
do i=1 to 4;
input age smoke wheeze @;
output;
end;
datalines;
1 portage 9 0 1 10 0 1 11 0 1 12 0 0
...
;
1
DATA six;
2
INPUT case city$ @;
3
DO i=1 to 4;
4
INPUT age smoke wheeze @;
5
OUTPUT;
6
END;
7
DATALINES;
8
1 portage 9 0 110 0 111 0 112 0 0
9
...
10
;
2 Bloque de código
PROC GENMOD
Explicación : Ejecuta un análisis de modelos lineales generalizados. El modelo especifica una distribución binomial para la variable dependiente 'wheeze'. La instrucción REPEATED maneja la correlación entre las observaciones del mismo sujeto ('case') utilizando una estructura de covarianza intercambiable, típica de los análisis GEE.
¡Copiado!
proc genmod data=six;
class case city;
model wheeze(event='1') = city age smoke / dist=bin;
repeated subject=case / type=exch covb corrw;
run;
1
PROC GENMODDATA=six;
2
class case city;
3
model wheeze(event='1') = city age smoke / dist=bin;
4
repeated subject=case / type=exch covb corrw;
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.
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.