The STDRATE Procedure

Attributable Fraction and Population Attributable Fraction

The attributable fraction measures the excess event rate or risk fraction in the exposed population that is attributable to the exposure. That is, it is the proportion of event rate or risk in the exposed population that would be reduced if the exposure were not present. In contrast, the population attributable fraction measures the excess event rate or risk fraction in the total population that is attributable to the exposure.

In the STDRATE procedure, you can compute the attributable fraction by using either indirect standardization or Mantel-Haenszel estimation.

Indirect Standardization

With indirect standardization, you specify a study population that consists of subjects who are exposed to a factor, such as smoking, and a reference population that consists of subjects who are not exposed to the factor. Denote the numbers of events in the study and reference populations by ${\mathcal D}_ s$ and ${\mathcal D}_ r$, respectively.

For the rate estimate, denote the population-times in the study and reference populations by ${\mathcal T}_ s$ and ${\mathcal T}_ r$, respectively. Then the event rates in the two populations can be expressed as the following equations, respectively:

\[  {\hat\lambda }_ s= \,  \frac{\,  {\mathcal D}_ s \, }{{\mathcal T}_ s} \; \; \; \; \;  \mbox{and} \; \; \; \; \;  {\hat\lambda }_ r= \,  \frac{\,  {\mathcal D}_ r \, }{{\mathcal T}_ r}  \]

Similarly, for the risk estimate, denote the numbers of observations in the study and reference populations by ${\mathcal N}_ s$ and ${\mathcal N}_ r$, respectively. Then the event risks in the two populations can be expressed as the following equations, respectively:

\[  {\hat\gamma }_ s= \frac{\,  {\mathcal D}_ s \, }{{\mathcal N}_ s} \; \; \; \; \;  \mbox{and} \; \; \; \; \;  {\hat\gamma }_ r= \frac{\,  {\mathcal D}_ r \, }{{\mathcal N}_ r}  \]

In the next two subsections, $\beta =\lambda $ denotes the rate statistic and $\beta =\gamma $ denotes the risk statistic.

Attributable Fraction with Indirect Standardization

The attributable fraction is the fraction of event rate or risk in the exposed population that is attributable to exposure:

\[  {\mathcal R}_ a = \frac{{\hat\beta }_ s - {\hat\beta }_ r}{{\hat\beta }_ s}  \]

With a standardized mortality ratio ${\mathcal R}_{sm}$, the attributable fraction is estimated by

\[  {\mathcal R}_ a = \frac{{\mathcal R}{sm}-1}{{\mathcal R}_{sm}}  \]

The confidence intervals for the attributable fraction can be computed using the confidence intervals for ${\mathcal R}_{sm}$. That is, with a confidence interval $({\mathcal R}_ l, \;  {\mathcal R}_ u)$ for ${\mathcal R}_{sm}$, the corresponding ${\mathcal R}_ a$ confidence interval is given by

\[  \left( \frac{{\mathcal R}_ l - 1}{{\mathcal R}_ l}, \;  \frac{{\mathcal R}_ u - 1}{{\mathcal R}_ u} \right)  \]

Population Attributable Fraction with Indirect Standardization

The population attributable fraction for a population is the fraction of event rate or risk in a given time period that is attributable to exposure. The population attributable fraction is

\[  {\mathcal R}_{pa} = \frac{ {\hat\beta }_0 - {\hat\beta }_ r }{ {\hat\beta }_0 }  \]

where

\[  {{\hat\beta }_0} = \frac{{\mathcal D}_ s + {\mathcal D}_ r}{{\mathcal T}_ s + {\mathcal T}_ r}  \]

is the combined rate in the total population for the rate statistic and where

\[  {{\hat\beta }_0} = \frac{{\mathcal D}_ s + {\mathcal D}_ r}{{\mathcal N}_ s + {\mathcal N}_ r}  \]

is the combined risk in the total population for the risk statistic.

Denote $\rho = {\mathcal D}_ s / ({\mathcal D}_ s+{\mathcal D}_ r)$, the proportion of exposure among events, then ${\mathcal R}_{pa}$ can also be expressed as

\[  {\mathcal R}_{pa} = \,  \rho \;  \,  \frac{{\mathcal R}_{sm}-1}{{\mathcal R}_{sm}}  \]

where ${\mathcal R}_{sm}$ is the standardized mortality ratio.

An approximate confidence interval for the population attributable rate ${\mathcal R}_{pa}$ can be derived by using the complementary log transformation (Greenland, 2008, p. 296). That is, with

\[  {\mathcal H} = \mbox{log}(1 - {\mathcal R}_{pa})  \]

a variance estimator for the estimated ${\mathcal H}$ is given by

\[  \mbox{Var}({\hat{\mathcal H}}) = \frac{{\mathcal R}_{pa}^{2}}{(1-{\mathcal R}_{pa})^{2}} \left( \frac{\hat{V}}{({\mathcal R}_{sm}-1)^{2}} + \frac{2}{{\mathcal D}_ s \,  ({\mathcal R}_{sm}-1)} + \frac{{\mathcal D}_ r}{{\mathcal D}_ s \,  ({\mathcal D}_ s+{\mathcal D}_ r)} \right)  \]

where ${\hat{V}}$ is a variance estimate for $\mbox{log}({\mathcal R}_{sm})$.

Mantel-Haenszel Estimation

With Mantel-Haenszel estimation, you specify one study population that consists of subjects who are exposed to a factor and another study population that consists of subjects who are not exposed to the factor. Denote the numbers of events in the exposed and nonexposed study populations by ${\mathcal D}_1$ and ${\mathcal D}_2$, respectively.

For the rate estimate, denote the population-times in the two populations by ${\mathcal T}_1$ and ${\mathcal T}_2$, respectively. Then the event rates in the two populations can be expressed as the following equations, respectively:

\[  {\hat\lambda }_1= \,  \frac{\,  {\mathcal D}_1 \, }{{\mathcal T}_1} \; \; \; \; \;  \mbox{and} \; \; \; \; \;  {\hat\lambda }_2= \,  \frac{\,  {\mathcal D}_2 \, }{{\mathcal T}_2}  \]

Similarly, for the risk estimate, denote the numbers of observations in the two populations by ${\mathcal N}_1$ and ${\mathcal N}_2$, respectively. Then the event risks in the two populations can be expressed as the following equations, respectively:

\[  {\hat\gamma }_1= \frac{\,  {\mathcal D}_1 \, }{{\mathcal N}_1} \; \; \; \; \;  \mbox{and} \; \; \; \; \;  {\hat\gamma }_2= \frac{\,  {\mathcal D}_2 \, }{{\mathcal N}_2}  \]

In the next two subsections, $\beta =\lambda $ denotes the rate statistic and $\beta =\gamma $ denotes the risk statistic.

Attributable Fraction with Mantel-Haenszel Estimation

The attributable fraction is the fraction of event rate or risk in the exposed population that is attributable to exposure:

\[  {\mathcal R}_ a = \frac{{\hat\beta }_1 - {\hat\beta }_2}{{\hat\beta }_1}  \]

Denote the rate or risk ratio by ${\mathcal R}= {\hat\beta }_1 / {\hat\beta }_2$. The attributable fraction is given by

\[  {\mathcal R}_ a = \frac{{\mathcal R}-1}{\mathcal R}  \]

The confidence intervals for the attributable fraction can be computed using the confidence intervals for the rate or risk ratio ${\mathcal R}$. That is, with a confidence interval $({\mathcal R}_ l, \;  {\mathcal R}_ u)$ for ${\mathcal R}$, the corresponding ${\mathcal R}_ a$ confidence interval is given by

\[  \left( \frac{{\mathcal R}_ l - 1}{{\mathcal R}_ l}, \;  \frac{{\mathcal R}_ u - 1}{{\mathcal R}_ u} \right)  \]

For Mantel-Haenszel estimation, you can use the Mantel-Haenszel rate or risk ratio to estimate ${\mathcal R}$.

Population Attributable Fraction with Mantel-Haenszel Estimation

The population attributable fraction for a population is the fraction of event rate or risk in a given time period that is attributable to exposure. The population attributable fraction is

\[  {\mathcal R}_{pa} = \frac{ {\hat\beta }_0 - {\hat\beta }_2 }{ {\hat\beta }_0 }  \]

where

\[  {{\hat\beta }_0} = \frac{{\mathcal D}_1 + {\mathcal D}_2}{{\mathcal T}_1 + {\mathcal T}_2}  \]

is the combined rate in the total population for the rate statistic and where

\[  {{\hat\beta }_0} = \frac{{\mathcal D}_1 + {\mathcal D}_2}{{\mathcal N}_1 + {\mathcal N}_2}  \]

is the combined risk in the total population for the risk statistic.

Denote the proportion of exposure among events as $\rho = {\mathcal D}_1 / ({\mathcal D}_1+{\mathcal D}_2)$. Then ${\mathcal R}_{pa}$ can also be expressed as

\[  {\mathcal R}_{pa} = \,  \rho \;  \,  \frac{{\mathcal R}-1}{\mathcal R}  \]

where ${\mathcal R}= {\hat\beta }_1 / {\hat\beta }_2$ is the rate or risk ratio.

An approximate confidence interval for the population attributable rate ${\mathcal R}_{pa}$ can be derived by using the complementary log transformation (Greenland, 2008, p. 296). That is, with

\[  {\mathcal H} = \mbox{log}(1 - {\mathcal R}_{pa})  \]

a variance estimator for the estimated ${\mathcal H}$ is given by

\[  \mbox{Var}({\hat{\mathcal H}}) = \frac{{\mathcal R}_{pa}^{2}}{(1-{\mathcal R}_{pa})^{2}} \left( \frac{\hat{V}}{({\mathcal R}-1)^{2}} + \frac{2}{{\mathcal D}_1 \,  ({\mathcal R}-1)} + \frac{{\mathcal D}_2}{{\mathcal D}_1 \,  ({\mathcal D}_1+{\mathcal D}_2)} \right)  \]

where ${\hat{V}}$ is a variance estimate for $\mbox{log}({\mathcal R})$.

For Mantel-Haenszel estimation, you can use the Mantel-Haenszel rate or risk ratio to estimate ${\mathcal R}$.