topic Re: Prroc phreg output statement in Statistical Procedures

Prroc phreg output statement

Yasu — Fri, 20 Jan 2023 13:21:24 GMT

Hi,

I'm applying two different cox models to the same dataset.

The linear predictor (XBETA) is different, but the survivor function estimate is exactly the same between two models.

Does this happen by chance?

Below are codes and outputs:

PROC PHREG DATA=&DATA.;
CLASS Z/PARAM=REF REF=FIRST;
MODEL TIME*EVENT(0)=Z P1_3 P2_3 / RL;
OUTPUT OUT=A SURVIVAL=SURV XBETA=XBETA;
RUN;

PROC PHREG DATA=&DATA.;
CLASS Z/PARAM=REF REF=FIRST;
MODEL TIME*EVENT(0)=Z SP1_3 SP2_3 / RL;
OUTPUT OUT=B SURVIVAL=SURV XBETA=XBETA;
RUN;

Re: Prroc phreg output statement

PaigeMiller — Fri, 20 Jan 2023 14:20:38 GMT

So the difference in the two models is the two different model statements:

MODEL TIME*EVENT(0)=Z P1_3 P2_3 / RL;
MODEL TIME*EVENT(0)=Z SP1_3 SP2_3 / RL;

Since you know, and we don't, what is different between P1_3 and SP1_3 and what is different between P2_3 and SP2_3, perhaps you could tell us what the differences are. In particular, if there is a mathematical function that gives an exact relationship between P2_3 and SP2_3 and an exact relationship between P1_3 and SP1_3, then the fact that the survivor function estimates are identical would be expected.

Re: Prroc phreg output statement

Yasu — Fri, 20 Jan 2023 14:45:02 GMT

Hi Paige,

Thank you for your reply.

They are imaginary variables that are normally distributed. Please imagine that P1_3 and P2_3 are office systolic and diastolic blood pressure, and SP1_3 and SP2_3 are home systolic and diastolic blood pressure.

Thanks,

Yasu

Re: Prroc phreg output statement

PaigeMiller — Fri, 20 Jan 2023 14:57:22 GMT

Were they created to have some mathematical relationship between the office variable and the home variable?

Re: Prroc phreg output statement

Yasu — Sat, 21 Jan 2023 01:25:39 GMT

They were simulated so that they have weak correlations among them.