Cómo probar una hipótesis de riesgos proporcionales Con SAS

Modelos de riesgos ayudan a los estadísticos y otros investigadores a evaluar el riesgo a través de la comprensión de las variables asociadas con el riesgo. Un tipo de modelo de riesgo, el modelo de riesgos proporcionales de Cox, es una forma simplificada y específica de este modelo que se basa en una suposición fuerte: las variables asociadas con el riesgo de peligros están relacionados de una manera multiplicativa (que se multiplican en el modelo y la adición ). Antes de que pueda utilizar con éxito un modelo de riesgos proporcionales de Cox, debe comprobar esta hipótesis. Usando SAS, un potente paquete de software estadístico, se puede probar fácilmente esta suposición en sus datos.

Instrucciones

1 Ingrese los datos en SAS. Haga clic en "Archivo" y "Importar datos". Siga los pasos que aparecen y se pueden abrir su archivo de datos.

2 Tipo de "proc = PHREG de datos; yourdata", donde "yourdata" es la matriz de datos que se va a probar. Esto le indica a SAS los datos que va a trabajar con.

3 El doble de las variables en los datos, multiplicando cada una con el registro de tiempo. Para cada variable, utilice el comando "newvar = oldvar

log (tiempo);". Por ejemplo, si tiene una variable para el sueldo llamado "salario", introducir el comando "newsalary = log salario (tiempo);". Hacer esto para todas las variables del conjunto de datos.

4 Tipo "* censor de tiempo del modelo (0) = yourmodel", donde "yourmodel" es el conjunto de variables que desee incluir en el modelo. Este conjunto de variables debe incluir todas las variables originales en sus datos, así como las nuevas variables que se multiplicaron por el registro de tiempo.

5 Preparar el criterio de proporcionalidad para las nuevas variables. Tipo "proportionality_test: newvariables de prueba;", donde "newvariables" es el conjunto de todas las nuevas variables creadas por la multiplicación de las variables de edad por el registro de tiempo. Por ejemplo, si tiene dos variables en su estudio, el salario y la edad, tendrá dos nuevas variables, y newsalary newage. Por lo tanto, se debe introducir "proportionality_test: Prueba newsalary newage;".

6 Ejecutar la prueba. Tipo de "correr;".

7 Observe los valores de p en la salida. Los valores de p se encuentran en la columna entitiled "Pr> chisq". Si hay valores de p inferior a 0,05, la prueba falla, es decir, la asunción de riesgos proporcionales de falla.