The ICPHREG Procedure

Hazard Ratios

Consider a dichotomous risk factor variable X that takes the value 1 if the risk factor is present and 0 if the risk factor is absent. The log-hazard function is

\[  \log [\lambda (t|X)] = \log [\lambda _0(t)] + \beta _1 X  \]

where $\lambda _0(t)$ is the baseline hazard function.

The hazard ratio $\psi $ is defined as the ratio of the hazard for those who have the risk factor (X = 1) to the hazard for those who do not have the risk factor (X = 0). The log of the hazard ratio is

\[  \log (\psi ) \equiv \log [\psi (X=1,X=0)] = \log [\lambda (t|X=1)] - \log [\lambda (t|X=0)] = \beta _1  \]

In general, the hazard ratio can be computed by exponentiating the difference of the log-hazard between any two population profiles. This is the approach taken by the HAZARDRATIO statement , so the computations are available regardless of parameterization, interactions, and nestings. However, as shown in the preceding equation for $\log (\psi )$, hazard ratios of main effects can be computed as functions of the parameter estimates. The remainder of this section is concerned with this methodology.

The parameter $\beta _1$ that is associated with X represents the change in the log-hazard from X = 0 to X = 1. So the hazard ratio is obtained by simply exponentiating the value of the parameter that is associated with the risk factor. The hazard ratio indicates how the hazard changes as you change X from 0 to 1. For example, $\psi =2$ means that the hazard when X = 1 is twice the hazard when X = 0.

Suppose the values of the dichotomous risk factor are coded as constants a and b instead of 0 and 1. The hazard when $X=a$ becomes $\lambda (t)\exp (a\beta _1)$, and the hazard when $X= b$ becomes $\lambda (t)\exp (b\beta _1)$. The hazard ratio that corresponds to an increase in X from a to b is

\[  \psi = \exp [(b - a) \beta _1] = [\exp (\beta _1)]^{b-a} \equiv [\exp (\beta _1)]^ c  \]

Note that for any a and b such that $c=b-a=1, \psi =\exp (\beta _1)$. So the hazard ratio can be interpreted as the change in the hazard for any increase of one unit in the corresponding risk factor. However, the change in hazard for some amount other than one unit is often of greater interest. For example, a change of one pound in body weight might be too small to be considered important, whereas a change of 10 pounds might be more meaningful. The hazard ratio for a change in X from a to b is estimated by raising the hazard ratio estimate for a unit change in X to the power of $c=b-a$ as shown previously.

For a polytomous risk factor, the computation of hazard ratios depends on how the risk factor is parameterized. For illustration, suppose that Cell is a risk factor that has four categories: Adeno, Large, Small, and Squamous.

For the effect parameterization scheme (PARAM= EFFECT) with Squamous as the reference group, the design variables for Cell are as follows:

 

Design Variables

Cell

$X_1$

$X_2$

$X_3$

Adeno

1

0

0

Large

0

1

0

Small

0

0

1

Squamous

–1

   –1

–1

The log-hazard for Adeno is

\begin{eqnarray*}  \log [\lambda (t|\mbox{Adeno})] & =&  \log [\lambda _0(t)] + \beta _1 (X_1=1) + \beta _2 (X_2=0) + \beta _3 (X_3=0) \\ & =&  \lambda _0(t) + \beta _1 \end{eqnarray*}

The log-hazard for Squamous is

\begin{eqnarray*}  \log [\lambda (t|\mbox{Squamous})] & =&  \log [\lambda _0(t)] + \beta _1 (X_1=-1) + \beta _2 (X_2=-1)+ \beta _3 (X_3=-1)) \\ & =&  \log [\lambda _0(t)] - \beta _1 - \beta _2 - \beta _3 \end{eqnarray*}

Therefore, the log-hazard ratio of Adeno versus Squamous

\begin{eqnarray*}  \log [\psi (\mbox{Adeno},\mbox{Squamous})] & =&  \log [\lambda (t|\mbox{Adeno})] - \log [\lambda (t|\mbox{Squamous})] \\ & =&  2 \beta _1 + \beta _2 + \beta _3 \end{eqnarray*}

For the reference cell parameterization scheme (PARAM= REF) in which Squamous is the reference cell, the design variables for Cell are as follows:

 

Design Variables

Cell

$X_1$

$X_2$

$X_3$

Adeno

1

0

0

Large

0

1

0

Small

0

0

1

Squamous

0

   0

0

The log-hazard ratio of Adeno versus Squamous is

\begin{eqnarray*}  \lefteqn{\log (\psi (\mbox{Adeno},\mbox{Squamous}))} \\ &  = &  \log [\lambda (t|\mbox{Adeno})] - \log [\lambda (t|\mbox{Squamous})] \\ &  = &  (\log [\lambda _0(t)] + \beta _1 (X_1=1) + \beta _2 (X_2=0) + \beta _3 (X_3=0)) - \\ & &  (\log [\lambda _0(t)] + \beta _1 (X_1=0) + \beta _2 (X_2=0) + \beta _3 (X_3=0)) \\ &  = &  \beta _1 \end{eqnarray*}

For the GLM parameterization scheme (PARAM= GLM), the design variables are as follows:

 

Design Variables

Cell

$X_1$

$X_2$

$X_3$

$X_4$

Adeno

1

0

0

0

Large

0

1

0

0

Small

0

0

1

0

Squamous

0

0

0

1

The log-hazard ratio of Adeno versus Squamous is

\begin{eqnarray*}  \lefteqn{\log (\psi (\mbox{Adeno},\mbox{Squamous}))} \\ &  = &  \log [\lambda (t|\mbox{Adeno})] - \log [\lambda (t|\mbox{Squamous})] \\ &  = &  \log [\lambda _0(t)] + \beta _1 (X_1=1) + \beta _2 (X_2=0) + \beta _3 (X_3=0) + \beta _4 (X_4=0)) -\\ & &  (\log (\lambda _0(t)) + \beta _1 (X_1=0) + \beta _2 (X_2=0) + \beta _3 (X_3=0) + \beta _4(X_4=1)) \\ &  = &  \beta _1 - \beta _4 \end{eqnarray*}

Consider Cell as the only risk factor. The computation of the hazard ratio of Adeno versus Squamous for various parameterization schemes is shown in Table 51.7.

Table 51.7: Hazard Ratio of Adeno to Squamous

 

Parameter Estimates

 

PARAM=

$\hat{\beta _1}$

$\hat{\beta _2}$

$\hat{\beta _3}$

$\hat{\beta _4}$

Hazard Ratio Estimates

EFFECT

0.5772

–0.2115

0.2454

 

$\exp (2 \times 0.5772 -0.2115 + 0.2454) = 3.281$

REF

1.8830

0.3996

0.8565

 

$\exp (1.8830) =3.281 $

GLM

1.8830

0.3996

0.8565

0.0000

$\exp (1.8830) = 3.281 $


The fact that the log-hazard ratio ($\log (\psi )$) is a linear function of the parameters enables the HAZARDRATIO statement to compute the hazard ratio of the main effect even in the presence of interactions and nest effects.

To customize hazard ratios for specific units of change for a continuous risk factor, you can use the UNITS= option in a HAZARDRATIO statement to specify a list of relevant units for each explanatory variable in the model. Estimates of these customized hazard ratios are shown in a separate table. Let $(V_ j,U_ j)$ be a confidence interval for $\log (\psi )$. The corresponding lower and upper confidence limits for the customized hazard ratio $\exp (c\beta _ j)$ are $\exp (cV_ j)$ and $\exp (cU_ j)$, respectively for $c>0$, or $\exp (cU_ j)$ and $\exp (cV_ j)$, respectively for $c<0$.

Let $\mb{e}_ j$ be the jth unit vector—that is, the jth entry of the vector is 1 and all other entries are 0. The hazard ratio for the explanatory variable with regression coefficient $\beta _ j=\mb{e}_ j’\bbeta $ is defined as $\textrm{exp}(\beta _ j)$. In general, a log-hazard ratio can be written as $\mb{h}’\bbeta $ (a linear combination of the regression coefficients), and the hazard ratio $\textrm{exp}({\mb{h}’\bbeta })$ is obtained by replacing $\mb{e}_ j$ with $\mb{h}$.

Point Estimate

The hazard ratio $\textrm{exp}(\mb{e}_ j’\bbeta )$ is estimated by $\textrm{exp}(\mb{e}_ j’\hat{\bbeta })$, where $\hat{\bbeta }$ is the maximum likelihood estimate of the regression coefficients $\bbeta $.

Wald’s Confidence Limits

The $100(1-{\alpha })\% $ confidence limits for the hazard ratio are calculated as

\[  \mr{exp} \left(\mb{e}_ j’ \hat{\bbeta } \pm z_{\alpha /2} \sqrt {\mb{e}_ j’\bSigma _{\hat{\bbeta }}\mb{e}_ j} \right)  \]

where $\bSigma _{\hat{\bbeta }}$ is estimated covariance matrix and $z_{\alpha /2}$ is the $100(1-{\alpha }/2)$ percentile point of the standard normal distribution.