Survival function confidence interval prediction following PROC LIFEREG with Weibull distribution

[TomHsiung]

proc lifereg data = d202_raw order = formatted;
format x1 x2 x3 x4 x5 grpfmt.;
class x2 x3 x4;
model time*censor(0) = x1 x2 x3 x4 x5 / distribution = weibull;
run;

The PROC LIFEREG procedure does tell me the coefficients for variables and the Weibull shape. However, now I would like to draw the survival function curve with a confidence interval. The point estimation is easy to figure out by inputting the point estimations of these variables and the Weibull shape parameter. However, how to figure out the confidence interval of the survival curve? Hitherto, I only used the lower and upper limits of the Weibull shape parameter to do so. Much appreciated.