How do I get the predictive survival probability based on proc lifereg and PHREG

[desireatem]

data Headache;
input Minutes Group Censor @@;
datalines;
11 1 0 12 1 0 19 1 0 19 1 0
19 1 0 19 1 0 21 1 0 20 1 0
21 1 0 21 1 0 20 1 0 21 1 0
20 1 0 21 1 0 25 1 0 27 1 0
30 1 0 21 1 1 24 1 1 14 2 0
16 2 0 16 2 0 21 2 0 21 2 0
23 2 0 23 2 0 23 2 0 23 2 0
25 2 1 23 2 0 24 2 0 24 2 0
26 2 1 32 2 1 30 2 1 30 2 0
32 2 1 20 2 1
;

proc lifereg data=Headache;
class Group;
model Minutes*Censor(1)=Group/Dist=Weibull;
run;

proc phreg data=Headache;

class group;

model Minutes*Censor(1)=Group;
run;

What do I need to do to have the predictive survival probability ?

[Reeza]

Have you tried adding an OUTPUT statement?

Example

https://documentation.sas.com/?docsetId=statug&docsetTarget=statug_lifereg_examples02.htm&docsetVers...

Documentation for OUTPUT statement

https://documentation.sas.com/?docsetId=statug&docsetTarget=statug_lifereg_syntax10.htm&docsetVersio...

---

@desireatem wrote:

data Headache;
input Minutes Group Censor @@;
datalines;
11 1 0 12 1 0 19 1 0 19 1 0
19 1 0 19 1 0 21 1 0 20 1 0
21 1 0 21 1 0 20 1 0 21 1 0
20 1 0 21 1 0 25 1 0 27 1 0
30 1 0 21 1 1 24 1 1 14 2 0
16 2 0 16 2 0 21 2 0 21 2 0
23 2 0 23 2 0 23 2 0 23 2 0
25 2 1 23 2 0 24 2 0 24 2 0
26 2 1 32 2 1 30 2 1 30 2 0
32 2 1 20 2 1
;

proc lifereg data=Headache;
class Group;
model Minutes*Censor(1)=Group/Dist=Weibull;
run;

proc phreg data=Headache;

class group;

model Minutes*Censor(1)=Group;
run;

What do I need to do to have the predictive survival probability ?

---

[desireatem]

How about for proc PHREG predictive probability.

[Reeza · Accepted solution]

Same thing, did you try that?

[desireatem]

I can get the predictive probability using proc lifereg . Now my question is how about for phreg procedure

[Reeza]

Did you try adding an OUTPUT statement as well to your PROC PHREG?