In reply to Rick Wicklin.

Dear Rick,

It is possible to retrieve the estimated intercepts using the OUTPUT statement together with the LINP option (for example: output out=want LINP).
From the retrieved values one can compute the ‘lambda’ scale parameters for different Weibull distributions as, for example, lambda_1 = exp(Linp_1) (as you also indicate above).

Do you know which parameter estimates to retrieve, using the OUTPUT statement, and from these how to compute the shape parameters ‘a’ for different Weibull distributions?

Thanks!