The PHREG Procedure

The Frailty Model

You can use the frailty model to model correlations between failures of the same cluster. The hazard rate for the jth individual in the ith cluster is

\[  \lambda _{ij}(t) = \lambda _0(t) \mr{e}^{\bbeta '\bZ _{ij}(t) + \gamma _ i}  \]

where $\lambda _0(t)$ is an arbitrary baseline hazard rate, $\bZ _{ij}$ is the vector of (fixed-effect) covariates, $\bbeta $ is the vector of regression coefficients, and $\gamma _ i$ is the random effect for cluster i.

Frailties are the exponential transformations of the random components, and the frailty model can be written as

\[  \lambda _{ij}(t) = \lambda _0(t) \mr{e}^{\gamma _ i}\mr{e}^{\bbeta '\bZ _{ij}(t)}  \]

The random components $\gamma _1,\ldots ,\gamma _ s$ (alternatively, the frailties $\mr{e}^{\gamma _1},\ldots ,\mr{e}^{\gamma _ s}$) are assumed to be independent and identically distributed. Modeling is based on the random effects rather than on the frailties.

Two frailty distributions are available in PROC PHREG: gamma and lognormal. Use the DIST= option in the RANDOM statement to choose the distribution. Let $\theta $ be an unknown parameter. The frailty distributions are listed in Table 73.11.

Table 73.11: Frailty Distributions

Frailty Option

Distribution

Density f($\gamma _ i)$

Mean and Variance

DIST=GAMMA

$\mr{e}^{\gamma _ i} \sim G\left(\frac{1}{\theta },\frac{1}{\theta }\right)$

$ \frac{\exp (\frac{\gamma _ i}{\theta }) \exp \left(-\frac{\exp (\gamma _ i)}{\theta } \right)}{\theta ^{\frac{1}{\theta }}\Gamma (\frac{1}{\theta })} $

E($\mr{e}^{\gamma _ i}$)=1   V($\mr{e}^{\gamma _ i}$)=$\theta $

DIST=LOGNORMAL

$\gamma _ i \sim N(0, \theta )$

$ \frac{1}{\sqrt (2\pi \theta )}\exp \left(-\frac{\gamma _ i^2}{2\theta }\right) $

E($\gamma _ i$)=0    V($\gamma _ i$)=$\theta $


The unknown parameter $\theta $ is a dispersion parameter. Each frailty distribution has a central tendency of 1 (the gamma frailty has a mean of 1, and the lognormal frailty has a median of 1). Thus, you can infer that individuals in cluster i with frailty $\mr{e}^{\gamma _ i}>1$ (or $\mr{e}^{\gamma _ i}<1$) tend to fail at a faster (or slower) rate than they fail under an independence model.

PROC PHREG estimates the regression coefficients $\bbeta $, the random effects $\gamma _1,\ldots ,\gamma _ s$, and the dispersion parameter $\theta $. The RANDOM statement in PROC PHREG enables you to fit a shared frailty model by a penalized partial likelihood approach . If you also specify the BAYES statement, PROC PHREG performs a Bayesian analysis of the shared frailty model .

If the RANDOM statement is specified, any ASSESS, BASELINE, and OUTPUT statements are ignored. Also ignored are the COVS options in the PROC PHREG statement and the following options in the MODEL statement: BEST=, DETAILS, HIERARCHY=, INCLUDE=, NOFIT, PLCONV=, SELECTION=, SEQUENTIAL, SLENTRY=, SLSTAY=, TYPE1, and TYPE3(ALL, LR, SCORE). Profile likelihood confidence intervals for the hazard ratios are not available for the frailty model analysis.