The LIFEREG Procedure

Predicted Values

For a given set of covariates, $\mb{x}$ (including the intercept term), the pth quantile of the log response, $y_ p$, is given by

\[  y_ p = \mb{x}^{\prime } \bbeta + \sigma u_ p  \]

if no offset variable has been specified, or

\[  y_ p = \mb{x}^{\prime } \bbeta + \mr{o} + \sigma u_ p  \]

for a given value o of an offset variable, where $u_ p$ is the pth quantile of the baseline distribution. The estimated quantile is computed by replacing the unknown parameters with their estimates, including any shape parameters on which the baseline distribution might depend. The estimated quantile of the original response is obtained by taking the exponential of the estimated log quantile unless the NOLOG option is specified in the preceding MODEL statement.

The following table shows how $u_ p$ is computed from the baseline distribution $F(u)$:

Distribution

$F(u)$

$u_ p$

Exponential

$1-\exp (-\exp (u))$

$\log ( -\log ( 1 - p ) )$

Generalized Gamma

$\left\{  \begin{array}{lcl} \frac{ \Gamma \left( \delta ^{-2}, \delta ^{-2} \exp (\delta u) \right) }{ \Gamma \left( \delta ^{-2} \right) } & &  \mr{if } \  \delta > 0 \\ 1-\frac{ \Gamma \left( \delta ^{-2}, \delta ^{-2} \exp (\delta u) \right) }{ \Gamma \left( \delta ^{-2} \right) } & &  \mr{if } \  \delta < 0 \\ \end{array} \right.$

$F^{-1}(p)$

Logistic

$1 - ( 1 + \exp ( u ) )^{-1}$

$\log ( p / ( 1 - p ) )$

Log-logistic

$1 - ( 1 + \exp ( u ) )^{-1}$

$\log ( p / ( 1 - p ) )$

Lognormal

$\Phi (u)$

$\Phi ^{-1}(p)$

Normal

$\Phi (u)$

$\Phi ^{-1}(p)$

Weibull

$1-\exp (-\exp (u))$

$\log ( -\log ( 1 - p ))$

For the generalized gamma distribution, $u_ p$ is computed numerically.

The standard errors of the quantile estimates are computed using the estimated covariance matrix of the parameter estimates and a Taylor series expansion of the quantile estimate. The standard error is computed as

\[  \mr{STD} = \sqrt {\mb{z}^{\prime }\mb{Vz}}  \]

where $\mb{V}$ is the estimated covariance matrix of the parameter vector $(\bbeta ^{\prime },\sigma ,\delta )^{\prime }$, and $\mb{z}$ is the vector

\[  \mb{z} = \left[ \begin{array}{c} \mb{x} \\[0.05in] \hat{u}_ p \\[0.05in] \hat{\sigma } \frac{\partial u_ p}{\partial \delta } \\ \end{array} \right]  \]

where $\delta $ is the vector of the shape parameters. Unless the NOLOG option is specified, this standard error estimate is converted into a standard error estimate for $\exp (y_ p)$ as $\exp (\hat{y}_ p)$STD. It might be more desirable to compute confidence limits for the log response and convert them back to the original response variable than to use the standard error estimates for $\exp (y_ p)$ directly. See Example 57.1 for a 90% confidence interval of the response constructed by exponentiating a confidence interval for the log response.

The variable CDF is computed as

\[  \mr{CDF}_ i = F(u_ i)  \]

where the residual is defined by

\[  u_ i=\left( \frac{ y_ i-\mb{x}^{\prime }_ i\mb{b} }{ \hat{\sigma } } \right)  \]

and F is the baseline cumulative distribution function. If the data are interval-censored, then the cumulative distribution function, $\mr{CDF}_ i = F(u_ i)$, is evaluated at the lower interval endpoint.