El script comienza creando un conjunto de datos SAS© a partir de datos brutos incluidos directamente en el código. Estos datos, inicialmente en formato ancho, se reestructuran a formato largo para facilitar el análisis de medidas repetidas en el tiempo. Posteriormente, se utilizan dos procedimientos estadísticos, PROC SURVEYLOGISTIC y PROC GENMOD, para ajustar modelos logit acumulativos. Este enfoque permite modelar la variable de respuesta ordinal (el nivel de dolor) teniendo en cuenta la correlación entre las observaciones de un mismo sujeto en diferentes momentos.
Análisis de datos
Type : CREATION_INTERNE
Los datos se generan completamente dentro del script. Un paso DATA lee las observaciones a través de una instrucción 'datalines' y las transforma de un formato ancho (una fila por sujeto) a un formato largo (una fila por sujeto por punto de medición).
1 Bloque de código
DATA STEP Data
Explicación :
Este bloque DATA STEP crea la tabla 'Shoulder_tip_pain'. Lee los datos integrados a través de 'datalines', asigna un identificador único 'subject_id' a cada paciente y luego utiliza un bucle para transformar las 6 mediciones de dolor (t1-t6) de un formato ancho a un formato largo, con una fila por punto temporal ('time').
Este bloque DATA STEP crea la tabla 'Shoulder_tip_pain'. Lee los datos integrados a través de 'datalines', asigna un identificador único 'subject_id' a cada paciente y luego utiliza un bucle para transformar las 6 mediciones de dolor (t1-t6) de un formato ancho a un formato largo, con una fila por punto temporal ('time').
¡Copiado!
| 1 | DATA Shoulder_tip_pain; |
| 2 | INPUT trt $ gender $ age t1-t6; |
| 3 | subject_id = _n_; |
| 4 | array tt t1-t6; |
| 5 | DO over tt; |
| 6 | y = tt; |
| 7 | time = _i_; |
| 8 | OUTPUT; |
| 9 | END; |
| 10 | DATALINES; |
| 11 | y f 64 1 1 1 1 1 1 |
| 12 | y m 41 3 2 1 1 1 1 |
| 13 | y f 77 3 2 2 2 1 1 |
| 14 | y f 54 1 1 1 1 1 1 |
| 15 | y f 66 1 1 1 1 1 1 |
| 16 | y m 56 1 2 1 1 1 1 |
| 17 | y m 81 1 3 2 1 1 1 |
| 18 | y f 24 2 2 1 1 1 1 |
| 19 | y f 56 1 1 1 1 1 1 |
| 20 | y f 29 3 1 1 1 1 1 |
| 21 | y m 65 1 1 1 1 1 1 |
| 22 | y f 68 2 1 1 1 1 2 |
| 23 | y m 77 1 2 2 2 2 2 |
| 24 | y m 35 3 1 1 1 3 3 |
| 25 | y m 66 2 1 1 1 1 1 |
| 26 | y f 70 1 1 1 1 1 1 |
| 27 | y m 79 1 1 1 1 1 1 |
| 28 | y f 65 2 1 1 1 1 1 |
| 29 | y f 61 4 4 2 4 2 2 |
| 30 | y f 67 4 4 4 2 1 1 |
| 31 | y f 32 1 1 1 2 1 1 |
| 32 | y f 33 1 1 1 2 1 2 |
| 33 | n f 20 5 2 3 5 5 4 |
| 34 | n f 50 1 5 3 4 5 3 |
| 35 | n f 40 4 4 4 4 1 1 |
| 36 | n m 54 4 4 4 4 4 3 |
| 37 | n m 34 2 3 4 3 3 2 |
| 38 | n f 34 3 4 3 3 3 2 |
| 39 | n m 56 3 3 4 4 4 3 |
| 40 | n f 82 1 1 1 1 1 1 |
| 41 | n m 56 1 1 1 1 1 1 |
| 42 | n m 52 1 5 5 5 4 3 |
| 43 | n f 65 1 3 2 2 1 1 |
| 44 | n f 53 2 2 3 4 2 2 |
| 45 | n f 40 2 2 1 3 3 2 |
| 46 | n f 58 1 1 1 1 1 1 |
| 47 | n m 63 1 1 1 1 1 1 |
| 48 | n f 41 5 5 5 4 3 3 |
| 49 | n m 72 3 3 3 3 1 1 |
| 50 | n f 60 5 4 4 4 2 2 |
| 51 | n m 61 1 3 3 3 2 1 |
| 52 | ; |
2 Bloque de código
PROC FORMAT
Explicación :
Este procedimiento define dos formatos personalizados. El formato '$abc' se aplica a la variable 'trt' para etiquetar los grupos de tratamiento ('Active', 'Placebo'), y el formato '$xyz' se aplica a 'gender' para etiquetar el sexo ('Female', 'Male').
Este procedimiento define dos formatos personalizados. El formato '$abc' se aplica a la variable 'trt' para etiquetar los grupos de tratamiento ('Active', 'Placebo'), y el formato '$xyz' se aplica a 'gender' para etiquetar el sexo ('Female', 'Male').
¡Copiado!
| 1 | PROC FORMAT; |
| 2 | value $abc 'y' = 'Active' |
| 3 | 'n' = 'Placebo'; |
| 4 | value $xyz 'f' = 'Female' |
| 5 | 'm' = 'Male'; |
| 6 | RUN; |
3 Bloque de código
PROC SURVEYLOGISTIC
Explicación :
Este bloque utiliza PROC SURVEYLOGISTIC para ajustar un modelo GEE con un enlace logit acumulativo. La instrucción 'cluster subject_id' especifica que las observaciones están correlacionadas dentro de cada sujeto. El modelo evalúa el efecto del tratamiento, el sexo, la edad y el tiempo sobre la respuesta ordinal 'y'.
Este bloque utiliza PROC SURVEYLOGISTIC para ajustar un modelo GEE con un enlace logit acumulativo. La instrucción 'cluster subject_id' especifica que las observaciones están correlacionadas dentro de cada sujeto. El modelo evalúa el efecto del tratamiento, el sexo, la edad y el tiempo sobre la respuesta ordinal 'y'.
¡Copiado!
| 1 | ods html; |
| 2 | PROC SURVEYLOGISTIC DATA=Shoulder_tip_pain; |
| 3 | class trt gender; |
| 4 | model y = trt gender age time / link=clogit ; |
| 5 | cluster subject_id; |
| 6 | FORMAT trt $abc. gender $xyz. ; |
| 7 | title1 '*** Results from fitting a GEE Cumulative Logit Model ***'; |
| 8 | title2 '*** to the Shoulder Tip Pain Data in Lumley (1996) ***'; |
| 9 | RUN; |
| 10 | ods html close; |
4 Bloque de código
PROC GENMOD
Explicación :
Este bloque utiliza PROC GENMOD para realizar un análisis GEE alternativo o de confirmación. Especifica una distribución multinomial con un enlace logit acumulativo. La instrucción 'repeated' declara la naturaleza de las medidas repetidas por sujeto, con una estructura de correlación de trabajo independiente ('type=ind').
Este bloque utiliza PROC GENMOD para realizar un análisis GEE alternativo o de confirmación. Especifica una distribución multinomial con un enlace logit acumulativo. La instrucción 'repeated' declara la naturaleza de las medidas repetidas por sujeto, con una estructura de correlación de trabajo independiente ('type=ind').
¡Copiado!
| 1 | PROC GENMOD DATA=Shoulder_tip_pain; |
| 2 | class subject_id trt gender; |
| 3 | model y = trt gender age time / dist=mult link=clogit type3; |
| 4 | repeated subject=subject_id / type=ind; |
| 5 | FORMAT trt $abc. gender $xyz. ; |
| 6 | title3 '*** Results using PROC GENMOD ***'; |
| 7 | 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.
