BookmarkSubscribeRSS Feed
kivanvan
Obsidian | Level 7

In my PH Cox model (using PROC PHREG), one of my continuous covariate (EGFR) violates the PH assumption. Based on the plot of Schoenfeld residual below, I decide to use a heaviside function with the time cut-off at 2 (the 600th rank of time) on this variable.

SASq.png

 

I believe that, with the heaviside function, I still assume the proportional hazard in each time interval. I tried the ZPH option and ASSESS statement, and I got the note: "Model assessment is not available when there are time-dependent covariates. The ASSESS statement is ignored." I further tried to output Schoenfeld residual (see code below), but I got no observation in the output dataset. Does anyone know how to assess the PH assumption in this model?

 

 PROC PHREG DATA= DF;
  CLASS EXPO SEX;
  MODEL TIME*SA(0,)= EXPO AGE SEX  EGFRT1 EGFRT2 / RL=PL;
   EGFRT1=0;
   IF TIME<2 THEN EGFRT1=EGFR;
   EGFRT2=0;
   IF TIME>=2 THEN EGFRT2=EGFR;
  OUTPUT OUT=SCHOEN RESSCH=SCHEXPO SCHAGE SCHSEX SCHEGFRT1 SCHEGFRT2;
RUN;

hackathon24-white-horiz.png

2025 SAS Hackathon: There is still time!

Good news: We've extended SAS Hackathon registration until Sept. 12, so you still have time to be part of our biggest event yet – our five-year anniversary!

Register Now

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 0 replies
  • 1252 views
  • 0 likes
  • 1 in conversation