Direct adjusted survival curve with ICPHREG

[yunfat]

Hi,

I would like to estimate direct adjusted survival curve for interval censored data using ICPHREG.

Unlike PHREG, ICPHREG doesn't have an option DIRADJ.

I calculated the point estimate by simply taking the mean of all the estimated survival probability using covariates in the dataset. However, I have no idea how to calculate SE of the estimated survival probability...

Here is my code;

data cov;
set &DS;
keep age_at_1stvisit SBP HDLCh TG BS UA &VA;
if (SBP=. or HDLCh=. or TG=. or BS=. or UA=.) then delete;
&VA=0; output;
&VA=1; output;
run;

proc icphreg data=&DS ;
class &VA (ref=first);
model (lower, upper) = age_at_1stvisit &VA SBP HDLCh TG BS UA / basehaz=splines (DF=1);
baseline out=a covariates=cov timelist=10 survival=s ;
run;

proc sort data=a;
by &VA;
run;

proc means data=a;
class &VA;
var s ;
output out=b mean=;
run;

proc print data=a;
run;

Thank you in advance!

Mitsuko