Understanding multivariable analysis

Hazards Regression

Survival

Log-rank test

Atopic dermatitis and osteoporosis

Survival Analysis

Log-rank test

P <0.0001, indicates a significant difference in the patient's survival in the 4 clinical stages.

Yes, this agrees with the plot and statistical test. The survival rates are different for patients in various clinical stages.

Stage 2 – 3.739. Death risk is 3.739 times higher in Stage 2 compared to Stage 1.

Stage 3 – 3.659. The risk of death is 3.659 times higher in Stage 3 compared to Stage 1.

Stage 4 – 13.851. The risk of death is 13.851 times higher in Stage 4 compared to Stage 1.

This agrees with the previous findings.

Stage 2 – 3.848. The risk of death is 3.848 times higher in Stage 2 compared to Stage 1.

Stage 3 – 4.115. The risk of death is 4.115 times higher in Stage 3 compared to Stage 1.

Stage 4 – 15.893. The risk of death is 15.893 times higher in Stage 4 compared to Stage 1.

The risk of death is 1.070 times higher among patients who don’t have insurance covers.

The risk of death is 1.459 times higher for Black patients.

The risk of death is 1.436 times higher for Hispanic patients.

Increase in the Annual hospital patient volume by increasing the risk of death by 0.981.

The confounders increase the risk ratios between clinical stage and time to mortality.

Code for Survival

*Upload the dataset and run proc import statement*;

*Kaplan-Meier Method for time to mortality; 

ods graphics on;

proclifetest data=work.mortality plots=s(atrisk=0 to 42 by 6);

timedx_lastcontact_death_months*mortality(0); label dx_lastcontact_death_months='Months from Diagnosis';

run;

ods graphics off;

*Kaplan-Meier Method for time to mortality by clinical stage (1, 2, 3, and 4); 

ods graphics on;

proclifetest data=work.mortality plots=s(nocensoratrisk=0 to 36 by 6) maxtime=36;

timedx_lastcontact_death_months*mortality(0); label dx_lastcontact_death_months='Months from Diagnosis';

strataClin_Stage; label Clin_Stage='Clinical Stage';

run;

ods graphics off;

*Univariate Cox Proportional Hazards Model for Clinical Stage*; 

procphreg data=work.mortality; class clin_stage(ref='1') /param=reference;

modeldx_lastcontact_death_months*mortality(0) = clin_stage /rl;

run;

*Multivariable Cox Proportional Hazards Model for Clinical Stage*; 

procphreg data=work.mortality; class clin_stage(ref='1') /param=reference;

modeldx_lastcontact_death_months*mortality(0) = clin_stagenoinsur black hispanic volume /rl;

run;

Information on Study Data, including Data Dictionary:

NAME: Cancer Mortality Data (Mortality.csv)

SIZE: 200 observations, 8 variables

DATA DICTIONARY: