 Letter to the Editor Response
 Open
 Published:
Estimation and application of population attributable fraction in ecological studies
Environmental Healthvolume 18, Article number: 52 (2019)
The original article was published in Environmental Health 2019 18:54
Abstract
Estimation of population attributable fraction (PAF) requires unbiased relative risk (RR) by using either Levin’s or Miettinen’s formula, on which decision depends on the available exposure information in reference group, not the types of studies. For ecological studies and studies with aggregated outcomes, once having unbiased RRs, Levin’s and Miettinen’s formulae would provide identical PAF estimates. PAF could also be applied to compare relative burdens of disease between countries across time, which is an additional information in consideration of countrylevel policies.
To the Editor
Population attributable fraction (PAF) is widely used to measure the disease burden attributable to a given risk factor. Concerns on improper estimation of PAF in ecological studies are raised [1] in consideration of potential ecological bias. However, if unbiased relative risks (RR) are available, estimation of PAF from ecological studies is feasible, either by using Levin’s or Miettinen’s formula.
In a recent published ecological study [2], for example, the PAF is referred to the proportion of subjects with lung cancer that would have not occurred if coalfired power plants (the exposure of interest) were absent counterfactually, assuming the same probability of getting lung cancer in the exposed and unexposed groups with the remaining risk factors. Considering a randomly selected country i, the numbers of subjects with and without lung cancer in the exposed and unexposed groups can be presented as Table 1.
where O = observed numbers of subjects with lung cancer; E = expected numbers of subjects with lung cancer had everyone not been exposed in country i.
PAF could also be expressed as either Levin’s or Miettinen’s formula, respectively.
\( \mathrm{PAF}=\frac{P_e\times \left( RR1\right)}{P_e\times \left( RR1\right)+1}\kern0.5em \)(Levin’s formula) [3].
\( ={P}_c\times \frac{RR1}{RR} \) (Miettinen’s formula) [4].
where.
\( {P}_e=\frac{X}{X+Y} \), the proportion of subjects being exposed in the population;
\( \mathrm{RR}=\frac{\frac{a}{X}}{\frac{c}{Y}}=\frac{a\times Y}{X\times c} \), relative risk of lung cancer comparing the subjects among exposed and unexposed groups;
\( {P}_c=\frac{a}{a+c} \), the proportion of exposed cases among lung cancer subjects.
Mathematically, the PAF calculated by using Levin’s formula would be identical to that from Miettinen’s formula, regardless whether P_{e} = P_{c} = 1 or not, as proved below.
Levin’s formula:
Miettinen’s formula:
One major reason that researchers commonly use Levin’s formula to estimate PAF in casecontrol studies and use Miettinen’s formula in cohort studies, respectively, is because of data availability. In casecontrol studies, researchers have sufficient information of the four cells in Table 1. Whereas in cohort studies and studies with aggregated data (e.g., standardized incidence ratio (SIR) and standardized mortality ratio (SMR) studies), most researchers do not have the exposure information among subjects without diseases (i.e. b and d cells in Table 1), and thus P_{e} is not available. Once having sufficient information in observational studies, the use of either Miettinen’s or Levin’s formula would yield an identical PAF estimate.
Furthermore, although both requiring unbiased RRs for estimation of PAF in Levin’s and Miettinen’s formulae, methods to obtain the unbiased RRs are different according to study types, such as using confounding adjustment (mostly from regression models) in casecontrol studies and stratification (e.g., SMR stratified by age and sex) in observational studies [5]. Theoretically, unbiased RRs from either fully adjustment (i.e., no residual confounding) or fully stratification (i.e., no residual confounding within stratum) should be identical and feasible for PAF estimation [6]. World Health Organization and Institute for Health Metrics and Evaluation applied a hybrid method, which age and sexstratified RRs were retrieved by prior metaanalysis or regression models, and summarized in standardized populations to estimate PAFs (i.e., global burden of diseases) [7]. Given fully adjusted and unbiased estimation, RRs derived from ecological studies and other studies with aggregated outcomes would be as valid as those from casecontrol studies, and therefore, are legit for PAF estimation.
Lastly, PAF estimates could be interpreted as relative strength of a relationship between exposure and disease, regardless of the nature of association or causation [8], and subsequently be applied to compare relative burden of diseases across countries/populations. In Lin’s study [2], for example, relative burden of lung cancer between countries across time contribute valuable information in consideration of countrylevel policies.
In conclusion, valid PAF estimation could be achieved from both Levin’s and Miettinen’s formulae with sufficient information in different types of studies with unbiased RR estimation, regardless stratification or adjustment. Comparison of PAFs between countries across time might provide additional information, along with the point estimate per se.
Availability of data and materials
Not applicable.
Abbreviations
 PAF:

Population attributable fraction
 RR:

Relative risks
 SIR:

Standardized incidence ratio
 SMR:

Standardized mortality ratio
References
 1.
Khosravi, A. & Mansournia, M.A., Issues with incorrect computing of population attributable fraction (PAF) in a global perspective on coalfired power plants and burden of lung cancer, Environmental Health, 2019, https://doi.org/10.1186/s1294001904906
 2.
Lin CK, et al. A global perspective on coalfired power plants and burden of lung cancer. Environ Health. 2019;18(1):9.
 3.
Levin ML. The occurrence of lung cancer in man. Acta Unio Int Contra Cancrum. 1953;9(3):531–41.
 4.
Miettinen OS. Proportion of disease caused or prevented by a given exposure, trait or intervention. Am J Epidemiol. 1974;99(5):325–32.
 5.
Rockhill B, Newman B, Weinberg C. Use and misuse of population attributable fractions. Am J Public Health. 1998;88(1):15–9.
 6.
Morgenstern H. Ecologic studies in epidemiology: concepts, principles, and methods. Annu Rev Public Health. 1995;16:61–81.
 7.
Murray CJL, L.A., The Global Burden of Disease: a comprehensive assessment of mortality and disability from diseases, injuries and risk factors in 1990 and projected to 2020. 1996, Cambridge, MA USA: Harvard University Press.
 8.
Flegal KM, Panagiotou OA, Graubard BI. Estimating population attributable fractions to quantify the health burden of obesity. Ann Epidemiol. 2015;25(3):201–7.
Acknowledgements
We sincerely thank Professor HsienHo Lin and Dr. Nicholas Luo for their professional opinion on global burden disease and PAF estimations.
Funding
We declared no funding for the letter.
Author information
Affiliations
Contributions
CKL contributed the idea formation, manuscript preparation. STC contributed manuscript preparation. Both authors read and approved the final manuscript.
Corresponding author
Correspondence to ChengKuan Lin.
Ethics declarations
Ethics approval and consent to participate
Not applicable. No individual data is used in the letter.
Consent for publication
Not applicable. No individual data is used in the letter.
Competing interests
We declared no competing interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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
Received
Accepted
Published
DOI
Keywords
 Population attributable fraction
 Ecological study
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.