 Research
 Open Access
 Published:
Accounting for biases in surveybased estimates of population attributable fractions
Population Health Metrics volume 17, Article number: 19 (2019)
Abstract
Background
This paper discusses best practices for estimating fractions of mortality attributable to health exposures in survey data that are biased by observed confounders and unobserved endogenous selection. Extant research has shown that estimates of population attributable fractions (PAF) from the formula using the proportion of deceased that is exposed (PAF_{pd}) can attend to confounders, whereas the formula using the proportion of the entire sample exposed (PAF_{pe}) is biased by confounders. Research has not explored how PAF_{pd} and PAF_{pe} equations perform when both confounding and selection bias are present.
Methods
We review equations for calculating PAF based on either the proportion of deceased (pd) or the proportion of the entire sample (pe) that receives the exposure. We explore how estimates from each equation are affected by confounding bias and selection bias using hypothetical data and realworld survey data from the National Health Interview Survey–Linked Mortality Files, 1987–2011. We examine the association between cigarette smoking and allcause mortality risk in the US adult population as an example.
Results
We show that both PAF_{pd} and PAF_{pe} calculate the true PAF in the presence of confounding bias if one uses the “weightedsum” approach. We further show that both the PAF_{pd} and PAF_{pe} calculate biased PAFs in the presence of collider bias, but that the bias is more severe in the PAF_{pd} formula.
Conclusion
We recommend that researchers use the PAF_{pe} formula with the weightedsum approach when estimates of the exposureoutcome relationship are biased by endogenous selection.
Background
This paper discusses best practices for estimating the fraction of mortality attributable to health exposures in surveybased data that are biased by both observed confounders and unobserved endogenous selection. Much extant work has reviewed errors in computing population attributable fractions (PAFs) in the presence of confounders [1,2,3,4,5,6], but little work has considered how different formulae for computing PAFs are affected by endogenous selection biases (e.g., collider bias).
Endogenous selection bias can affect estimates of statistical associations in many ways. Conditioning on a collider variable—that is, a variable caused by two other variables that are associated with the exposure and the outcome—can occur through statistical control, stratification of the sample into different groups, or the selection of participants into a study [7,8,9,10,11]. Introducing collider variables through any of these mechanisms can bias estimates of associations between exposure and outcome. In this study, we focus on unobserved endogenous selection—a problem that commonly occurs in health studies through the sampling process of recruiting study participants. Simply put, the likelihood of participation in a health study can be affected by both the exposure and outcome, which can bias estimates of the true association between them.
The most common PAF formulae are based on either the proportion of deceased (pd) in the sample that receives the exposure or the proportion of the entire sample (pe) that receives the exposure [1]. The two main aims of this investigation are to examine the performance of these modelbased methods for calculating PAF in the presence of (1) known and observable confounders of the exposuremortality association and (2) collider bias. We focus on the association between cigarette smoking and allcause mortality risk in the US adult population, which is confounded by other variables and also a likely contributor to unobserved endogenous selection bias in surveybased data of smoking and mortality risk [7].
Methods
We use hypothetical data and realworld survey data to calculate PAF in the presence of confounding and unobserved endogenous selection. In all of our exercises, nonexposed cases are respondents who have never smoked cigarettes and exposed cases are respondents who are current or former smokers. The association of interest is how smoking affects allcause mortality risk. For each exercise, we estimate the fraction of US mortality attributable to cigarette smoking using the PAF_{pd} formula:
where pd is the prevalence of a health exposure among the deceased cases and RR is the mortality risk ratio between the exposed and nonexposed subjects [12]. We also estimate this fraction using the PAF_{pe} formula:
where pe is the prevalence of the exposure among all cases in the sample [6, 13]. For each formula, we adopt a “weightedsum” approach [1, 3, 14, 15], which uses modelbased adjusted estimators of PAF separately for each adjustment level i as well as the distribution of cases by the adjustment levels:
where i indicates the adjustment level (i.e., confounder) and W_{i} indicates the proportion of deaths in adjustment level i. The weightedsum approach is mathematically equivalent to the PAF_{pd} [1, 14]; combining the PAF_{pd} with the weightedsum approach is therefore redundant. Nevertheless, we apply it in all of our exercises to maintain consistency.
Exercise 1: Observed confounding bias in hypothetical data
In our first exercise, we examine PAF estimates from Eqs. (1–4) in the presence of a single confounder, race/ethnicity. For simplicity, we consider race/ethnicity using only two categories, nonHispanic black and nonHispanic white (hereafter black and white). The hypothetical data are composed of 1000 black respondents and 4000 white respondents. Both smoking prevalence (i.e., pe) and mortality risk are higher among black respondents than among white respondents, which confound the smokingmortality association. In these data, black pe is 0.35 compared to white pe of 0.2, and overall mortality risk for black respondents is 0.3 compared to 0.2 for white respondents.
Exercise 2: Unobserved endogenous selection bias in hypothetical data
Our second example uses the same data as before, but presupposes that estimates of the smokingmortality association are biased by differential selection into the sample. We assume that current smokers sampled are relatively more select on health than are nonsmokers. That is, both the nonsmoking and the smoking samples are healthier than the true populations, but the difference between the smoking sample and the smoking population is greater than the difference between the nonsmoking sample and the nonsmoking population. This unobserved process of health selection biases downward the allcause mortality RR estimated in the sample data. When these conditions hold, both PAF estimates will be biased due to the central role of the RRs (see Eqs. 1–4). Moreover, the distribution of deaths by exposure and by adjustment levels, W_{i}, will also be biased. This is because counts of deaths among the exposure group in the sample will be artificially low and, consequently, W_{i} will be incorrect. Thus, PAF_{pd} and PAF_{pe} estimates will remain biased via W_{i} even if our adjusted RRs account for collider bias. Finally, the estimated PAF from the PAF_{pd} formula will be additionally biased, due to the central role of the pd in the calculation of the PAF. That is, the pd in the observed data, like the RRs in the observed sample data, will be downwardly biased because deaths among the smoking sample are underreported.
Exercise 3: PAF estimation with realworld survey data
Finally, we illustrate the points above by analyzing the smokingmortality association in the National Health Interview Survey–Linked Mortality Files (NHISLMF) for years 1987–2009. These data are composed of NHIS waves from 1987 and 1989–2009 that have been linked to official death records at the National Death Index through December 31, 2011 (the 1988 NHIS survey did not contain information about respondents’ smoking behavior). The NHISLMF are designed to form a representative sample of noninstitutionalized US adults [12]. To simplify the example, we limit the analytic sample to contain only US adult black and white men and women aged 40 through 84 at time of interview and whose survival is followed between ages 50 and 84. We extend the example by considering two levels of smoking exposure, “former smoker” and “current smoker,” and by considering three possible confounders of the smokingmortality association: race/ethnicity (i.e., white and black), gender (i.e., men and women), and age group (i.e., 50–59, 60–69, 70–79, and 80–84).
We fit a series of cloglog discretetime survival models to estimate smokingbased differences in US adult mortality risk. First, we fit a baseline model that estimates differences in mortality risks between current, former, and never smokers (reference category). Next, we fit a confounder model that estimates agespecific differences in mortality risks between current, former, and never smokers, adjusting for race/ethnicity and gender as categorical confounders of the smokingmortality association. We also fit models separately for black and white men and women that estimate agespecific RRs for former and current smokers compared to never smokers (i.e., confounderspecific models to be used with the weightedsum approach to calculate PAFs). Finally, we fit a bias model that refits the confounder model by accounting for cohortbased variation in mortality risk and agerelated selection biases in the NHISLMF data.
Participants in health surveys like the NHIS are positively selected on survival, health, and noninstitutional living arrangements [16]. These selection biases tend to grow stronger with increasing age [17]. Thus, older respondents in NHISLMF data are selected on the outcome of interest (i.e., survival) and inclusion in the NHIS sampling frame (i.e., healthy and living in noninstitutionalized housing). Combined, the selective nature of the sample results in collider biases via agerelated selection into the sampling frame and the selective factors associated with age are likely stronger among respondents with health risk factors such as smoking than among healthy respondents [8].
Survival models fitted separately by cohort of entry into the NHIS sample provide evidence consistent with these assumptions about collider biases. For example, the estimated RR between current smokers and never smokers who died at age 70–80 ranges from 1.51 [1.43–1.58 95%CI] among respondents surveyed at age 70–75 to 4.11 [3.02–5.57 95%CI] among respondents surveyed at age 50–55. The bias model is a shared frailty survival model that estimates random effects variation in mortality risk by NHIS respondents’ 5year age cohorts at the time of sampling. Overall, the model fits agespecific mortality risks separately for current, former, and never smokers, adjusting for gender, race/ethnicity, birth year, and random effects for a 5year cohort of entry into the data.
Mortality differences between US adults selfreported to be current, former, and never smokers between ages 50 and 84 are estimated across these three models. We use the adjusted RRs between (1) current smokers and never smokers and (2) former smokers and never smokers, which are estimated from confounderspecific survival models and the weightedsum approach to calculate the PAF for smoking as a cause of death in the US adult black and white populations between ages 50 and 84 for years 1987–2011. For all models, we contrast PAFs calculated from PAF_{pd} with PAFs calculated from PAF_{pe} to examine how each formula is affected by (1) confounders in the estimated smokingmortality association and (2) collider bias.
The NHISLMF data analyzed for the current study are publicuse files made available by the NCHS (https://www.cdc.gov/nchs/datalinkage/mortality.htm). The analytic scripts (Additional file 1) and calculations to generate results (Additional file 2) for Exercise 3 are available in the appendix.
Results
Exercise 1: Observed confounding in hypothetical data
The confounding effect of race/ethnicity on the smokingmortality association is illustrated in Table 1. The allcause mortality RR for smoking when unadjusted for the confounding effects of race/ethnicity is (450/1150)/(650/3850) = 2.32. Alternatively, the RR adjusted for race/ethnicity is 2.23. That is, when we estimate separate RRs for each race/ethnicity sample, we observe
When these race/ethnicspecific RRs for smoking are standardized by the race/ethnic distribution of deaths and the race/ethnic distribution of smoking prevalence, the adjusted RR is 2.23. If one does not account for the confounding effects of race/ethnicity on both mortality risk and the probability of smoking, one would incorrectly estimate the PAF by the following:
 a)
Aggregating the probability of smoking to be (1150/5000) = 0.23,
 b)
Aggregating the probability of smoking among decedents to be (450/1100) = 0.41, and
 c)
Aggregating the RR associated with smoking to be (450/1150)/(650/3850) = 2.32.
As a result, estimates of the PAF for smoking, irrespective of the formula used, would be biased by not attending to the confounding effects of race/ethnicity:
The actual PAF shown in the counterfactual example above is (1100 − 856)/1100 = 0.222
Thus, by failing to account for (1) the higher prevalence of smoking among black respondents and (2) the higher mortality risks among black respondents, we would incorrectly inflate the RR associated with smoking and misattribute numerous deaths to smoking as a cause of mortality in the population. As such, it is necessary to identify the RR by accounting for confounders in model estimates, and then use this confounderadjusted RR to calculate PAFs [18]. It has been argued that only the PAF_{pd} formula can accurately estimate the PAF when using confounderadjusted RRs [3,4,5,6]. Yet, as others have noted, one can use the PAF_{pe} equation with the confounderadjusted RR to derive the true PAF [1, 15]. To do so, one needs to first estimate separate PAFs for each confounder group (i.e., each adjustment level i), and then standardize these confounderspecific PAFs by the distribution of deaths across groups (i.e., W_{i}).
To illustrate, when we estimate separate PAFs for black and white respondents, we see for black:
and for white:
To estimate the total PAF, we further attend to the distribution of deaths across groups. That is, we simply weight the confounderspecific PAFs by the proportion of total deaths occurring in the confounder groups (i.e., W_{i}) [18]. The proportion of the total deaths that occurred among black respondents = (300/1100) = 0.273 and the proportion of total deaths that occurred among white respondents = (800/1100) = 0.727. When we weight the confounderspecific PAFs by the proportion of deaths in the two groups, W_{i}, we retrieve the true overall PAF:
This shows that the weightedsum approach can calculate the true PAF regardless if one uses the PAF_{pd} or PAF_{pe} formula. So long as (1) unobservable confounders or unobservable selection do not induce bias, and (2) one attends to observable confounders of the smokingmortality association, one can use adjusted RRs with either PAF_{pd} or PAF_{pe} and the weightedsum approach to calculate the PAF for smokingrelated mortality in the sample [18].
Exercise 2: Unobserved endogenous selection in hypothetical data
In the next exercise, we extend the previous example to consider sample data that are biased by unobserved selection, causing underestimation of mortality risk in the smoking population. To simplify matters, let us assume that the prevalence of smoking is the same in both the sample and population so that the only change pertains to q_{x} for smokers in the sample. The new information about population parameters is presented in Table 2 below.
The mortality probabilities for the nonsmoking populations equal those in the sample data (0.231 among blacks and 0.156 among whites). Smoking prevalence is also the same (pe_{NHB} = 0.35 and pe_{NHW} = 0.20). However, we now see discrepancies in the mortality risks for the smoking populations (0.500 in the white population vs. 0.375 in the white sample, and 0.500 in the black population vs. 0.429 in the black sample). These, in turn, affect the RRs for smoking (e.g., 2.17 vs. 1.86 for black and 3.21 vs. 2.40 for white), the pds (0.538 vs. 0.500 for black and 0.444 vs. 0.375 for white), and the W_{i} (e.g., 0.265 of population deaths are among blacks vs. 0.273 of sample deaths).
The confounderspecific PAFs using both the PAF_{pd} and PAF_{pe} formulae are as follows (estimates might be slightly different due to rounding):
nonHispanic black:
nonHispanic white:
Standardizing these confounderspecific PAFs by the distribution of deaths, W_{i}, we use the weightedsum approach to calculate the true PAF:
We see that the PAFs in the sample data underestimate the true PAF in the population (0.222 vs. 0.301), and this bias is the same in the PAF_{pd} and PAF_{pe} formulae. The discrepancy arises from one’s inattention to (unobservable) endogenous selection bias in the sample data, resulting in biased sample estimates of the mortality RRs associated with smoking as well as biased W_{i} in the sample.
Imagine that we had accounted for unobservable selection bias in our survival models and correctly identified the RRs for smoking for both the black and white samples. Even though the adjusted RRs would be correct in our survival models, the counts of deaths in the sample data would remain biased. Consequently, the pd values in the sample stay at 0.50 and 0.375, and the proportion of deaths occurring among blacks and whites stay at 0.273 and 0.727, respectively. As a result, if we were to calculate the PAF using the adjusted RRs with PAF_{pd}, we would find
The confounderspecific PAFs are biased (i.e., 0.270 estimated vs. 0.290 actual for blacks and 0.258 estimated vs. 0.306 actual for whites) even when using the adjusted RRs. Furthermore, when we use the weightedsum approach and standardize these PAFs by W_{i}, we add another source of bias because the distribution of deaths in each confounder group is biased as well: total PAF = (0.270*0.273) + (0.258*0.727) = 0.263. Yet, were we to follow conventional wisdom [3,4,5,6] and use the adjusted RR with the PAF_{pd} for the entire sample, we would estimate the same biased PAF:
Thus, even if we accurately accounted for selection bias in our survival models and estimated an unbiased RR (e.g., by fitting frailty models that account for selection bias in the smoking RR [19]), the PAF calculated from the PAF_{pd} formula will still be biased. In this case, a biased 0.263 is estimated for the sample when the true PAF in the population is 0.301 (a bias on the proportionate scale of 12.6%: (0.263 − 0.301)/0.301).
If we calculate the PAF using the confounder and selectionadjusted RRs with the PAF_{pe} formula, we find
We see that the confounderspecific PAFs are unbiased. Only when we standardize these PAFs by the distribution of deaths, W_{i}, do we introduce slight bias in the total PAF = (0.290*0.273) + (0.306*0.727) = 0.302 (a bias on the proportionate scale of − 0.3%: (0.301 − 0.302)/0.302). Thus, when we account for selection bias in our survival models and estimate unbiased adjusted RRs, the PAF calculated from PAF_{pe} will be biased, but only via W_{i}. By using the PAF_{pe} equation, we avoid bias in estimates from the pd and dramatically reduce the overall bias in the PAF estimate (0.3% vs. 12.6%).
To recap, when sample data are biased by unobserved selection, both the PAF_{pd} formula and the PAF_{pe} formula will calculate a biased PAF—even if researchers adjust for selection bias in the data. However, the PAF_{pd} formula is far more affected by the bias than is the PAF_{pe} formula because bias is introduced in both the pd and W_{i}. Conversely, estimates of the confounderspecific PAF from the PAF_{pe} equation are not biased, but some bias is introduced in the weightedsum approach via W_{i}. Theoretically, one could completely eliminate bias by identifying the true RR (i.e., attend to both observable confounders and unobservable selection biases) and standardizing the PAFs by the true distribution of deaths for each adjustment level (i.e., use population data to estimate W_{i}).
Exercise 3: PAF estimation with realworld survey data
For the final exercise, we calculate PAF for smoking as a cause of US adult mortality in the NHISLMF data, which are biased by confounding (i.e., age, race/ethnicity, and gender) and likely biased by endogenous selection (i.e., likelihood of sample inclusion depends on health). Table 3 shows agespecific mortality risks between years 1987 and 2011 for NHIS respondents who are current, former, and never smokers. The pd for former smokers (0.352) combined with the pd for current smokers (0.338) indicates that nearly 70% of the deceased NHIS sample had been exposed to smoking.
From the sample data in Table 3, we calculate the unadjusted RRs:
Because we are calculating a PAF for twolevels of an exposure, former smokers and current smokers, the PAF formulae change slightly [18, 20]:
We see that if we did not consider age, race/ethnicity, or gender as confounders of the smokingmortality association in these NHISLMF data, we would estimate about 23% of US black and white adult deaths between ages 50 and 85 for years 1987–2011 were attributable to cigarette smoking.
Average RRs for current smoking estimated from cloglog discrete time hazard models are presented in Table 4, and overall PAFs estimated from the PAF_{pe} and PAF_{pd} formula are included as well.
The baseline model estimates mortality risks for former and current smokers relative to never smokers that match the RRs observed in Table 3 (i.e., 1.49 and 1.52, respectively). Using these RRs, we estimate the same 0.232 PAF for smoking as a cause of US adult mortality, regardless if we estimate the PAF from the PAF_{pe} formula or the PAF_{pd} formula. The confounder model estimates agespecific RRs for former and current smokers relative to never smokers while controlling for confounding by gender and race/ethnicity. The age patterns in the RRs for current smokers suggest that the mortality consequences of smoking significantly decline with age. For example, current smokers are estimated to have about 2.6 to 2.7 times the mortality risk as never smokers in agegroups 50–59 and 60–69, but only about 1.2 times the mortality risk in agegroup 80–84. When using these confounderadjusted and agespecific RRs for smoking, we estimate a 0.247 PAF for smoking as a cause of US adult mortality.
Finally, the estimated agespecific RRs from the bias model are significantly larger than the agespecific RRs from the confounder model, especially at older ages. Although the smokingmortality relationship attenuates with age, it is substantially less than the attenuation observed in the confounder model. Using these confounder and selectionadjusted RRs, we calculate a PAF of 0.289 from the PAF_{pd} formula and a PAF of 0.326 from the PAF_{pe} formula. This is the only case in which we observe different PAF values depending on the formula used. This is because the PAF_{pd} formula remains biased by pd and likely underestimates the amount of mortality attributable to cigarette smoking in the US adult population. In this case, the PAF estimated from the PAF_{pd} formula is likely additionally biased by − 11.3% over the PAF_{pe} (0.289 − 0.326)/0.326) because it does not fully account for collider bias in estimates of the smokingmortality association in the NHISLMF data.
Discussion
Betweengroup differences in mortality (e.g., smokers and nonsmokers) estimated from survey data are often biased by unobserved endogenous selection [8, 10]. These biases can distort research findings and lead to incorrect conclusions and misguided policy recommendations. Researchers should therefore be wary of collider biases and, when possible, adjust estimates to account for them. Relatedly, researchers should be wary of how these biases affect PAF calculations. In this paper, we demonstrated that the PAF_{pd} formula is far more sensitive to collider bias than the PAF_{pe} formula. Results from both our hypothetical examples and realworld illustration using the NHISLMF show the PAF_{pd} formula calculated severely biased estimates of the PAF for smoking as a cause of mortality. As such, if estimates of the exposureoutcome association are likely biased by endogenous selection, researchers should consider calculating PAFs using the PAF_{pe} formula with the weightedsum approach. The main challenge to using the weightedsum approach is the data required to scale estimates by W_{i}, which increase with the number of confounders in the model. In addition, the weightedsum approach may not be appropriate in small samples because estimates of W_{i} are unreliable [1].
The findings are important for researchers aiming to estimate the mortality burden of exposures that may induce collider bias in sample data. For example, estimates from the NHISLMF data indicate that widening educational disparities in US adult mortality have greatly increased deaths attributable to low educational attainment [21]. Yet, estimates of the educationmortality association in the NHISLMF data may be biased by mortality and health selection across age [22]. Deaths attributable to low education in the USA may, in fact, be underestimated by not accounting for collider bias in PAF calculations. Also, researchers have reported discrepant PAFs for obesity as a cause of US mortality. For example, Flegal et al. [5] review PAF values indicating 2–15% of adult deaths are attributable to high BMI. The discrepancies likely reflect the extent to which researchers attend to confounder and collider biases in model estimates and how these biases affect PAF calculations. While Flegal et al. ([5] p. 203) consider the PAF_{pe} to be “the invalid formula” and PAF_{pd} to be the “formula appropriate for use with adjusted relative risks when confounding exists,” their review did not consider how the PAF formulae were affected by collider bias. Results here indicate that the PAF_{pd} is, in fact, the formula that calculates more biased estimates when relative risks are adjusted for confounding and selection biases.
Conclusion
Many studies have addressed best practices for calculating and interpreting PAFs for causes of mortality [1, 3, 5, 6, 20, 23,24,25]. In this paper, we extend these discussions to consider how unobserved endogenous selection bias (e.g., collider bias) distorts calculations of PAFs in the PAF_{pd} and PAF_{pe} formulae. Prior research has highlighted the importance of confounding bias in PAF calculations, but it has not considered how collider bias may affect PAF calculations. We used both hypothetical and realworld data on the smokingmortality relationship to explore these considerations. Results from our examples demonstrate that both the PAF_{pd} and PAF_{pe} formulae can equally attend to observable confounders and accurately calculate PAFs via the weightedsum approach [1, 3, 18]. Yet, the PAF_{pe} formula via the weightedsum approach is preferred to the PAF_{pd} formula if RR estimates for the exposure are biased from endogenous selection. In contrast to conventional wisdom that recommends using the PAF_{pd} formula with adjusted RRs [3, 5, 6], we conclude by recommending the use of the PAF_{pe} formula with the weightedsum approach when using RRs adjusted for both confounding bias and selection bias.
Availability of data and materials
The datasets supporting the conclusions of this article are available as public use files of the NHISLMF data (https://www.cdc.gov/nchs/datalinkage/mortalitypublic.htm). The analytic scripts are available as “Additional file 1” and the PAF calculations are made available as excel files as “Additional file 2.”
Abbreviations
 95%CI:

95% confidence interval
 NHISLMF:

National Health Interview Survey–Linked Mortality Files
 PAF:

Population attributable fraction
 PAF_{pd} :

Population attributable fraction estimated from the formula using the proportion of deceased that is exposed
 PAF_{pe} :

Population attributable fraction estimated from the formula using the proportion of sample that is exposed
 pd:

The proportion of deceased exposed
 pe:

The proportion of sample exposed
 q _{x} (NS):

Probability of death among nonsmokers in population
 q _{x} (S):

Probability of death among smokers in population
 q _{ x } :

Probability of death
 RR:

Risk ratio
 W _{ i } :

The proportion of total deaths in each adjustment level
References
 1.
Benichou J. A review of adjusted estimators of attributable risk. Stat Methods Med Res. 2001;10(3):195–216.
 2.
Walter SD. Attributable risk in practice. Am J Epidemiol. 1998;148(5):411–3.
 3.
Darrow LA, Steenland NK. Confounding and bias in the attributable fraction. Epidemiology. 2011;22(1):53–8.
 4.
Flegal KM, Graubard BI, Williamson DF. Methods of calculating deaths attributable to obesity. Am J Epidemiol. 2004;160(4):331–8.
 5.
Flegal KM, Panagiotou OA, Graubard BI. Estimating population attributable fractions to quantify the health burden of obesity. Ann Epidemiol. 2015;25(3):201–7.
 6.
Rockhill B, Newman B, Weinberg C. Use and misuse of population attributable fractions. AJPH. 1998;88(1):15–9.
 7.
Schooling CM, Yeung SLA. “Selection bias by death” and other ways collider bias may cause the obesity paradox. Epidemiology. 2017;28(2):e16–7.
 8.
Greenland S. Quantifying biases in causal models: classical confounding vs colliderstratification bias. Epidemiology. 2003;14(3):300–6.
 9.
Elwert F, Winship C. Endogenous selection bias: the problem of conditioning on a collider variable. Annu Rev Sociol. 2014;40:31–53.
 10.
Flanders WD, Eldridge RC, McClellan W. A nearly unavoidable mechanism for collider bias with indexevent studies. Epidemiology. 2014;25(5):762–4.
 11.
Snoep JD, Morabia A, HernándezDíaz S, Hernán MA, Vandenbroucke JP. Commentary: A structural approach to Berkson’s fallacy and a guide to a history of opinions about it. Int J Epidemiol. 2014;43(2):515–21.
 12.
National Center for Health Statistics (NCHS). Office of Analysis and Epidemiology, Publicuse Linked Mortality File. Hyattsville; 2015. Available at the following address: http://www.cdc.gov/nchs/data_access/data_linkage/mortality.htm
 13.
Levin ML. The occurrence of lung cancer in man. Acta Unio Int Contra Cancrum. 1953;9:531–41.
 14.
Miettinen OS. Proportion of disease caused or prevented by a given exposure, trait or intervention. Am J Epidemiol. 1974;99:325–32.
 15.
Gefeller O. Comparison of adjusted attributable risk estimators. Stat Med. 1992;11(16):2083–91.
 16.
Keyes KM, Rutherford C, Popham F, Martins SS, Gray L. How healthy are survey respondents compared with the general population? Using surveylinked death records to compare mortality outcomes. Epidemiology. 2018;29(2):299–307.
 17.
Mendes de Leon CF. Aging and the elapse of time: a comment on the analysis of change. J Gerontol Ser B Psychol Sci Soc Sci. 2007;62(3):S198–202.
 18.
Bruzzi P, Green SB, Byar DP, Brinton LA, Schairer C. Estimating the population attributable risk for multiple risk factors using casecontrol data. Am J Epidemiol. 1985;122(5):904–14.
 19.
Vaupel JW, Manton KG, Stallard E. The impact of heterogeneity in individual frailty on the dynamics of mortality. Demography. 1979;16(3):439–54.
 20.
Hanley JA. A heuristic approach to the formulas for population attributable fraction. J Epidemiol Community Health. 2001;55(7):508–14.
 21.
Krueger PM, Tran MK, Hummer RA, Chang VW. Mortality attributable to low levels of education in the United States. PLoS One. 2015;10(7):e0131809.
 22.
Lynch SM. Cohort and lifecourse patterns in the relationship between education and health: a hierarchical approach. Demography. 2003;40(2):309–31.
 23.
Greenland S, Robins JM. Conceptual problems in the definition and interpretation of attributable fractions. Am J Epidemiol. 1988;128(6):1185–97.
 24.
Poole C. A history of the population attributable fraction and related measures. Ann Epidemiol. 2015;25(3):147–54.
 25.
Greenland S. Concepts and pitfalls in measuring and interpreting attributable fractions, prevented fractions, and causation probabilities. Ann Epidemiol. 2015;25(3):155–61.
Acknowledgements
We thank Bruce Link and Dan Powers for helpful comments and contributions to earlier works related to this paper, and to the referees for helpful comments and suggestions.
Consent to publication
Not applicable
Funding
We thank the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD)funded University of Colorado Population Center (Award Number P2C HD066613) for the development, administrative, and computing support. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NICHD or the National Institutes.
Author information
Affiliations
Contributions
Dr. ER and Dr. RM conceived of the paper together. Dr. RM wrote the bulk of the original text and carried out the analyses. Dr. ER provided invaluable comments and suggestions that guided the analyses, and also edited and rewrote much of the text. Both authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1.
This is a Stata dofile containing commands to fit cloglog discrete time survival models.
Additional file 2.
This is a Microsoft Excel file containing race/ethnic and genderspecific PAF estimates used to estimate PAF_{pe} and PAF_{pd} in Table 4.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Masters, R., Reither, E. Accounting for biases in surveybased estimates of population attributable fractions. Popul Health Metrics 17, 19 (2019). https://doi.org/10.1186/s1296301901966
Received:
Accepted:
Published:
Keywords
 Attributable fractions
 Selection bias
 Confounding bias
 Mortality