Using multi-year national survey cohorts for period estimates: an application of weighted discrete Poisson regression for assessing annual national mortality in US adults with and without diabetes, 2000–2006
© The Author(s). 2017
Received: 20 April 2016
Accepted: 9 December 2016
Published: 15 December 2016
Monitoring national mortality among persons with a disease is important to guide and evaluate progress in disease control and prevention. However, a method to estimate nationally representative annual mortality among persons with and without diabetes in the United States does not currently exist. The aim of this study is to demonstrate use of weighted discrete Poisson regression on national survey mortality follow-up data to estimate annual mortality rates among adults with diabetes.
To estimate mortality among US adults with diabetes, we applied a weighted discrete time-to-event Poisson regression approach with post-stratification adjustment to national survey data. Adult participants aged 18 or older with and without diabetes in the National Health Interview Survey 1997–2004 were followed up through 2006 for mortality status. We estimated mortality among all US adults, and by self-reported diabetes status at baseline. The time-varying covariates used were age and calendar year. Mortality among all US adults was validated using direct estimates from the National Vital Statistics System (NVSS).
Using our approach, annual all-cause mortality among all US adults ranged from 8.8 deaths per 1,000 person-years (95% confidence interval [CI]: 8.0, 9.6) in year 2000 to 7.9 (95% CI: 7.6, 8.3) in year 2006. By comparison, the NVSS estimates ranged from 8.6 to 7.9 (correlation = 0.94). All-cause mortality among persons with diabetes decreased from 35.7 (95% CI: 28.4, 42.9) in 2000 to 31.8 (95% CI: 28.5, 35.1) in 2006. After adjusting for age, sex, and race/ethnicity, persons with diabetes had 2.1 (95% CI: 2.01, 2.26) times the risk of death of those without diabetes.
Period-specific national mortality can be estimated for people with and without a chronic condition using national surveys with mortality follow-up and a discrete time-to-event Poisson regression approach with post-stratification adjustment.
KeywordsComplex survey Surveillance Mortality Survival analysis Discrete survival time Poisson regression
National surveillance of incidence, prevalence, and mortality is key to guiding and evaluating progress in chronic disease control and prevention. The prevalence of a chronic disease like diabetes can be affected by increasing incidence among persons without the disease as well as decreasing mortality among persons with the disease. Good estimates of mortality among persons with a chronic disease improve understanding of secular changes in prevalence, incidence, and mortality and their relationships. Since diabetes is often not recorded on death certificates as a direct, underlying, or contributing cause of death, the impact of diabetes on deaths in the United States population could be underestimated [1, 2]. The linkage of nationally representative surveys that include baseline disease status with mortality follow-up provides the opportunity to examine all-cause and cause-specific mortality among persons with diabetes or other chronic conditions.
From the policymaking and resource allocation perspectives, a cross-sectional estimate of mortality by calendar period (e.g., year) is highly desirable. Analyses of mortality follow-up data typically use survival approaches to examine the association between risk factors and death. In these analyses, the data are analyzed as a cohort covering the entire follow-up period, and the hazard of death is estimated for the cohort. However, this approach does not permit estimation of the hazard of death across time periods, nor does it provide valid annual or other calendar period estimates.
By following the conceptual framework of age-period-cohort analysis (APC) as represented by the Lexis diagram, multi-year cohort data can be decomposed into discrete time-to-event data and aggregated by calendar period [3, 4]. Calendar period all-cause mortality rates can be calculated by simply using the total number of deaths divided by the total person-years in each calendar period. Poisson regression, a generalized linear model, is appropriate for modeling unadjusted and adjusted mortality rates of multiple periods . Discrete Poisson regression yields identical estimates to the piecewise exponential model, which is another alternative to the Cox proportional hazards model . Nevertheless, most discrete time-to-event studies use aggregated group data and categorized independent variables ; we are not aware of previous publications using discrete Poisson regression applied to multi-year mortality follow-up data from national sample surveys.
In this study, to increase the awareness of estimating cross-sectional period mortality using multi-year national survey mortality follow-up data, we describe the construction of discrete survival time data in detail and demonstrate our approach from data preparation to data analysis. With diabetes as an example, we use population-weighted Poisson regression to model discrete survival time and estimate annual all-cause mortality by diabetes status. The US National Health Interview Survey (NHIS) 1997–2004 with mortality follow-up up to 2006 was used to illustrate this approach; US mortality estimates from the National Vital Statistics System (NVSS) were compared for validation purposes.
Continuous time-to-event data
Demonstration of continuous and discrete time-to-event survival formats
Part I: Continuous time-to-event format
Date of birth
Date of interview
Data of death or censored
Discretized time period
Age at interview (years)
(age at t0i)
Calendar year at interview
(year at t0i)
(1 = died)
(0 = lived)
(1 = Male)
(0 = Female)
Part II: Discrete time-to-event format
Person (i) & time period (j)
Date of birth
Date of interview
Data of death or censored
Discretized time period
Year of follow-up
Age at the beginning of year of follow-up
Calendar year during follow-up
(year at ij)
(1 = died)
(1 = Male)
(0 = Female)
Discrete time-to-event survival data
To calculate a period-specific mortality rate, we divided the continuous survival times into discrete calendar years. Since interviews did not all take place on the first day of the survey year, to make sure the survival time was allocated correctly we added an individual-specific partial time period (t_exti), calculated as the difference between the interview date and the first day of the year, then we divided the extended survival time [(t1−t0) + t_exti] into years. That is, each person’s total continuous survival time was discretized into multiple records, one for each calendar year. An individual’s survival time for a given calendar year was between 0 and 1 year, and the survival time in the first year was [1–t_exti]. In the analysis, age and calendar year were treated as time-varying (i.e., time-dependent) covariates. The age during each discrete period was assigned as the age on the first day of that calendar year.
With each participant contributing multiple discrete person-years during the follow-up, the sum of a person’s discretized annual person-years is equal to the total continuous survival time of that person (Fig. 1 – 1b). For example in Table 1 – Part II, person A was interviewed on 02/06/2000 at age 70.3 years and contributed 3.5 person-years of follow-up. That is, person A contributed 0.9 person-years in 2000 with age of 70.2 years at the beginning of year 2000, contributed 1 person-year in each of years 2001 and 2002 with age of 71.2, and 72.2 years at the beginning of those years, respectively, and in 2003, person A contributed 0.6 person-years before dying on 08/08/2003 at age 73.2 years at the beginning of the year. In this way, the continuous time-to-event survival records of participants have been decomposed into discrete survival time with time-varying age.
US national health interview survey and mortality follow-up
We used the NHIS mortality follow-up data to demonstrate our approach. The NHIS, conducted by the Centers for Disease Control and Prevention’s National Center for Health Statistics (NCHS), is an annual ongoing nationally representative cross-sectional household interview survey of US non-institutionalized civilians of all ages . The sampling plan covers the 50 states and the District of Columbia, and follows a multistage area probability design that permits the representative sampling of households and non-institutional group quarters. The annual response rate of NHIS is approximately 80% of the eligible households in the sample . All information about sex, race/ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, and others), and diabetes status was self-reported. Participants were classified as having diabetes if they answered “yes” to the question “Other than during pregnancy, have you EVER been told by a doctor or other health professional that you have diabetes or sugar diabetes?”
For our analysis, we selected the NHIS 1997 to 2004 surveys as the baseline, with mortality follow-up up through 2006; the NHIS 2005 to 2006 surveys were used to obtain the demographic distribution for the post-stratification reweighting of those two years. We included 307,280 adults aged 18 to 84 years from the NHIS 1997 to 2004 (range: 30,141 to 35,437 per year) and followed them through 2006. We excluded 15,882 (range: 1,723 to 2,642 per year) respondents because of insufficient identifying data to create a death status record, which yielded a final mortality follow-up sample of 242,397 (range: 29,076 to 29,193 per year) adults. The mortality follow-up sampling weights provided by NCHS accounted for excluded respondents.
Diabetic death was defined as a death with an associated International Classification of Diseases, 10th Revision (ICD-10) code of E10-E14. All-cause with diabetes death was defined as a person with diabetes who died of any cause. The total weighted person-time was used as the denominator for mortality calculation. We also estimated mortality by self-reported diagnosed diabetes at baseline. To validate our findings empirically, we compared the all-cause and diabetic mortality rates from NHIS with mortality rates from the NVSS, the fundamental source of US cause-of-death information. Mortality rates from the NVSS were directly calculated as number of death (all-cause or diabetic death coded as E10 to E14) divided by total population using structured query language from CDC WONDER by following the step-by-step instruction on the WONDER website (http://wonder.cdc.gov/mortSQL.html).
To reduce potential selection bias due to respondents being healthier than non-respondents, we excluded each individual’s first two years of follow-up. The final analytical discrete time-to-event data set included adults aged 20 years or older during the years 2000 to 2006.
Here, the natural logarithm of the expected value of the event, log(d), with an offset of natural logarithm of follow-up time, log(pt), is a linear combination of independent covariates, X i , with regression parameters β i ,.
Poisson regression provided the estimate of mortality for each calendar year/period. We used the robust error variances estimation approach to minimize over-dispersion  and the polynomial function of calendar time to smooth year-to-year variation in mortality rates [6, 13]. To smooth the variation in mortality due to low mortality rates in some age subgroups, the age at the beginning of a calendar year was defined as a continuous variable with polynomial terms (quadratic polynomial). The mortality rates in our study were estimated by the predictive margins of the regression coefficients from the Poisson model.
Adjusted sampling weights for the discrete time-to-event data
The age of sampled participants in each survey cohort increased with the year of follow-up and those multi-year survey cohorts also overlapped time periods. Without accounting for the demographic discrepancy between the participants from different cohorts and the US population at each specific year, the demographic distribution of a discrete period after the baseline year would not represent the demographic distribution of the US population at that specific year or time-period, and the total crude mortality of the US population would be biased toward the older population. In order to correct for these issues, we adjusted the sample weights using a post-stratification procedure in which sampled units were divided into subgroups based on age, sex, and race/ethnicity; we used the nationally representative weighted size of each subgroup of NHIS 2000 to 2006 at interview to estimate the US population size. The analysis weights for the discrete time-to-event data were reweighted proportionally. The adjusted analysis weights thus sum to the US population size within each subgroup. The sum of the analysis weights equaled the total non-institutionalized US population for each calendar year.
We used Stata 13.1 (StataCorp LP, College Station, Texas) to account for the complex multistage sampling design and to produce weighted estimates and 95% confidence intervals (CI).
For all comparisons we used a two-sided statistical test with significance defined as p value (p) <0.05 or a 95% CI that did not include the null value. The ggplot2 package of R was used to produce graphics .
All-cause mortality rates of US adults aged 20 to 84 years (per 1,000 person-years and 95% CI) by calendar year using discretized survival time, NVSS and NHIS follow-up
National Vital Statistics System
All-cause deaths, n
All-cause mortality (M0)a
NHIS mortality follow-up with the post-stratification sampling weights
All-cause deaths, n
Sample adults, n
All-cause mortality (M1)b
8.81 (8.00, 9.62)
8.45 (7.75, 9.14)
8.35 (7.77, 8.93)
8.39 (7.89, 8.89)
7.98 (7.54, 8.42)
8.11 (7.65, 8.58)
7.94 (7.56, 8.33)
NHIS mortality follow-up with original sampling weights
All-cause mortality (M2)b
8.99 (8.17, 9.82)
8.80 (8.08, 9.52)
8.67 (8.07, 9.28)
8.83 (8.31, 9.36)
8.41 (7.95, 8.87)
8.58 (8.10, 9.06)
8.36 (7.95, 8.77)
(M2 – M1)
0.18 (−0.97, 1.33)
0.35 (−0.65, 1.35)
0.32 (−0.53, 1.17)
0.44 (−0.29, 1.17)
0.43 (−0.19, 1.05)
0.47 (−0.18, 1.12)
0.42 (−0.15, 0.99)
To show the importance of post-stratification reweighting, we compared the NHIS follow-up estimates with the results from NVSS and the NHIS estimates that used the original weights (Table 2). Mortality estimates using the original sampling weights without post-stratification adjustment were higher than the mortality estimates using adjusted sampling weights, because of the aging of cohorts during the follow-up. Mortality in each year from the NVSS was within the 95% CIs of mortality rates from the NHIS using the adjusted sampling weights. The average annual decrease in crude mortality (per 1,000 person-years) was 0.12 for both the NHIS and the NVSS. The correlation of NHIS and the NVSS mortality was 0.94. Age-sex-race/ethnicity-adjusted mortality decreased 2.6% per year (p < 0.001).
Diabetica and all-cause with diabetesb mortality rates (per 1,000 person-years and 95% CI) of US adults aged 20 to 84 years by calendar year using discretized survival time data, NVSS and NHIS follow-up
National Vital Statistics System
Diabetic death, n
Diabetic mortality (M0)c
NHIS mortality follow-up: total adults with and without diabetes
Diabetic deaths, n
Diabetic mortality (M1)d
0.26 (0.14, 0.37)
0.29 (0.18, 0.41)
0.42 (0.28, 0.57)
0.34 (0.23, 0.44)
0.26 (0.18, 0.34)
0.30 (0.21, 0.38)
0.25 (0.18, 0.33)
All-cause with diabetes death, n
All-cause with diabetes mortality (M2)e
1.78 (1.43, 2.13)
1.94 (1.59, 2.30)
1.80 (1.51, 2.09)
1.80 (1.55, 2.05)
1.63 (1.42, 1.83)
1.70 (1.51, 1.88)
1.66 (1.48, 1.83)
NHIS mortality follow-up: adults with diagnosed diabetes
Diabetic deaths, n
Diabetic mortality (M3)f
3.50 (1.48, 5.51)
4.04 (2.08, 6.01)
7.03 (4.34, 9.72)
5.76 (3.85, 7.67)
4.04 (2.57, 5.51)
4.41 (2.94, 5.88)
3.06 (2.09, 4.03)
All-cause deaths, n
All-cause mortality (M4)g
35.7 (28.4, 42.9)
39.5 (32.4, 46.5)
36.1 (30.4, 41.8)
35.6 (30.7, 40.5)
31.7 (27.8, 35.6)
33.0 (29.4, 36.5)
31.8 (28.5, 35.1)
9.80 (9.05, 10.56)
10.13 (9.39, 10.86)
19.39 (18.76, 20.02)
16.29 (15.76, 16.82)
12.62 (12.20, 13.03)
13.37 (13.00, 13.75)
9.75 (9.40, 10.10)
NHIS mortality follow-up: adults without diagnosed diabetes
Diabetic death, n
Diabetic mortality (M5)h
0.09 (0.03, 0.15)
0.10 (0.03, 0.17)
0.08 (0.03, 0.13)
0.05 (0.02, 0.08)
0.05 (0.02, 0.09)
0.07 (0.04, 0.11)
0.10 (0.05, 0.16)
Period mortality among persons with chronic conditions such as diabetes is an important surveillance indicator of disease prevention and control. However, since chronic disease status is not reported in many vital statistics registries, it is often not possible to use vital statistics data to estimate mortality of persons with and without the condition. This presents a particular limitation for diabetes-related death statistics because diabetic death is often not recorded on US death certificates as a direct underlying or contributing cause of death, and diabetic deaths in the US population could be underestimated by solely using death certificate information [1, 2]. Assembly of national cohorts by linking national survey data with vital statistics provides a potential remedy to the data gap, but requires specific methods to permit estimation of period effects. In this study, we described the use of weighted discrete Poisson regression to estimate national mortality rates by diabetes status using a complex sample survey, and we validated this approach using mortality registry estimates. Our study showed that all-cause mortality of US adults estimated by the NHIS mortality follow-up decreased 2.6% per year from 2000 to 2006, which was similar to mortality rates estimated using the NVSS mortality registry data. Meanwhile all-cause mortality of US adults with diabetes decreased 3.7% per year during the same period [1, 2].
The method that is most often used to analyze mortality cohort data is the Cox proportional hazards regression model, which is useful for analyzing the data from the association or cause-effect relationship perspective. However, it is cumbersome to use this method to calculate hazard rates for a large number of combinations of predictors. Alternatively, parametric survival models can be more convenient for predicting, but cannot deal easily with time-varying covariates .
Age-period-cohort (APC) analysis provides a third option. If a vital statistics registry includes complete information on disease status, the APC method can be used to estimate the annual/period mortality among persons with and without diabetes [4, 6, 11, 16–18]. In the US, diabetes status is not recorded in the national vital statistics registry system. So we cannot apply this method directly. However, US nationally representative survey mortality follow-up data provide information on both diabetes status and death status. The APC model and life table framework can be applied to these data.
The APC analysis has been applied in demography, social science, and disease surveillance research using cross-sectional registration or survey data for a long time . The data are usually cross-sectional and grouped for data analyses. One of the major purposes of these studies was to separate the age, period, or cohort effects using cross-sectional data . In our study, we applied the concept and analytic framework of this widely used APC model. Compared to traditional APC models, our study had several differences. First, our study used longitudinal national complex survey mortality follow-up data. Second, the purpose was to estimate period mortality, which is a sum of the age and cohort effects. Finally, to account for the aging of the cohort during follow-up, we post-stratified the aggregated multiple segments from different survey cohorts using the US population structure at each period.
Both Poisson and logistic regression can be used for discretized time-to-event data analysis. Efron combined the logistic regression with discrete time-to-event survival time by 1-month intervals and obtained direct estimates of the hazard rates . A polynomial or spline model can be used to smooth out the random variation/noise. This partial logistic regression gives good estimates when the discrete time interval is small. Nevertheless, a Poisson regression that accounts for person-time of follow-up gives more accurate hazard rate estimates for longer discrete time intervals than a logistic regression. Poisson regression has been used frequently to compare mortality rates among different categories of cohorts in epidemiological studies and is a convenient alternative to Cox proportional hazards regression especially when the proportional hazards assumptions are not met . Early studies on the analysis of cohort survival data showed that Poisson regression is a straightforward and intuitive approach for directly estimating the hazard rates while incorporating time scale as a covariate in the model [16, 22]. We were interested in annual (or longer) time periods rather than monthly or daily periods and thus discrete Poisson regression was chosen for our analysis.
To obtain valid national estimates from a complex sample survey, it is critical to use proper statistical methods to account for the sample design and sampling weights. Our study shows that in later years, the distribution of age in the follow-up cohort shifted to the right; thus without post-stratification reweighting, the overall mortality rates combining all ages would have been overestimated. Using the US population as the standard population for post-stratification re-weighting yielded all-cause and diabetic mortality estimates that were similar to the national registry estimates. Our study demonstrated that discrete Poisson regression with post-stratification is a feasible approach for estimating annual mortality for the US population with and without diabetes.
The major limitation of our approach is the amount of time needed to discretize and analyze a large sample with long follow-up time. Poisson regression using complex sample data is computationally time-consuming with large discretized person-time datasets because data cannot be collapsed over covariates to account for the design-based analysis of complex sample data. Estimation based on a small number of events can create problems with model convergence. Without careful programming and reweighting, the results can be biased. In addition, the NHIS mortality data represented deaths among the civilian non-institutionalized population with person-year as the denominator, whereas mortality data from the NVSS represented deaths among the entire US population with the whole population at risk as the denominator. Thus, mortality rates from the two systems might have subtle differences.
To demonstrate our approach, we used self-reported diabetes. While any self-reported condition is subject to recall error, the self-report of diabetes is considered a valid measure of diagnosed diabetes . Although it is recognized as being non-sensitive, it has been shown to be highly specific . Another source of bias may arise from the lack of information about diabetes status between the baseline interview and death or censoring. Even though the rate of remission from diabetes to non-diabetes is likely small , the lack of information on incident cases would likely lead to an overestimation of diabetes duration. Furthermore, if incident cases have a higher mortality rate than non-cases and a lower mortality rate than prevalent cases, then lacking this information on incidence could lead to an overestimation of mortality rates for the populations both with and without diabetes. Future analyses with information with multiple follow-up visits could quantify the impact of this bias. We demonstrated that weighted discrete Poisson regression is an efficient applicable approach to estimate period mortality from the national mortality follow-up data. To our knowledge, there has been no similar report, though all the steps of this approach are well established. Several reasons could explain the scant usage of the discrete Poisson regression approach, including lack of data availability, lack of its inclusion as part of biostatistics educational curricula, the computing time required to analyze discrete time-to-event data, and the complex sampling design of national surveys, which further complicates using this approach. However, the increasing availability of more powerful statistical software and computing capabilities permits a revisitation of this method for the analysis of national survey mortality follow-up data.
We conclude that combining national follow-up cohorts from multiple survey years and analyzing them using population weighted discrete Poisson regression can yield annual national mortality rates by disease status.
Centers for disease control & prevention
Wide-ranging online data for epidemiologic research
National center for health statistics
National health interview survey
National vital statistics system
A special thanks to the women and men who participated in the study. We would also like to thank all of the staff involved in the US National Health Interview Survey for the study design, data collection and data dissemination. We thank Dr. Giuseppina Imperatore and Dr. Elizabeth Luman for their contributions to this study.
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
This paper was not presented anywhere outside of the CDC.
Availability of data and materials
All the data are available online: 1) CDC WONDER (http://wonder.cdc.gov/mortSQL.html); 2) NHIS Public-Use file: http://www.cdc.gov/nchs/nhis/nhis_questionnaires.htm; and 3) NHIS Public-Use Linked Mortality Files: http://www.cdc.gov/nchs/data-linkage/mortality-public.htm.
YJC conceived the study, acquired and analyze the data, and draft the manuscript. YJC and TJT designed the study. YJC, TJT, EWG, and DBR revised and interpreted the manuscript critically for important intellectual content. YJC, TJT, EWG, and DBR read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
The NHIS mortality follow-up and CDC Wonder data used in this study were approved public data sets. We did not merging any of the data sets in such a way that individuals might be identified, and did not enhance the public data sets with any identifiable or potentially identifiable data.
Ethics approval and consent to participate
This is a secondary analysis of public-use national survey data. The content and conduct of the US NHIS were subject to an institutional review board. All individuals participating in the original survey were required to sign an informed consent document.
Open AccessThis 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.
- US Centers for Disease Control and Prevention. 2014 US National Diabetes Statistics Report Atlanta (GA): Centers for Disease Control and Prevention; [updated May 15, 2015; cited 2015 February 19, 2016.]. Available from: http://www.cdc.gov/diabetes/data/statistics/2014StatisticsReport.html.
- Miniño AM, Murphy SL, Xu J, Kochanek KD. Deaths: final data for 2008. Natl Vital Stat Rep. 2011;59(10):1–126.Google Scholar
- Frost WH. The age selection of mortality from tuberculosis in successive decades. 1939. Am J Epidemiol. 1995;141(1):4–9. discussion 3.View ArticlePubMedGoogle Scholar
- Keiding N. Statistical Inference in the Lexis Diagram. Philos T Roy Soc A. 1990;332(1627):487–509.View ArticleGoogle Scholar
- Laird N, Olivier D. Covariance analysis of censored survival data using log-linear analysis techniques. J Am Stat Assoc. 1981;76(374):231–40. PMID: WOS:A1981LT69100002.View ArticleGoogle Scholar
- Hansen MB, Jensen ML, Carstensen B. Causes of death among diabetic patients in Denmark. Diabetologia. 2012;55(2):294–302.View ArticlePubMedGoogle Scholar
- Gail MH, Graubard B, Williamson DF, Flegal KM. Comments on ‘Choice of time scale and its effect on significance of predictors in longitudinal studies’ by Michael J. Pencina, Martin G. Larson and Ralph B. D'Agostino, Statistics in Medicine 2007; 26:1343–1359. Stat Med. 2009;28(8):1315–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Korn EL, Graubard BI, Midthune D. Time-to-event analysis of longitudinal follow-up of a survey: choice of the time-scale. Am J Epidemiol. 1997;145(1):72–80.View ArticlePubMedGoogle Scholar
- Pencina MJ, Larson MG, D’Agostino RB. Choice of time scale and its effect on significance of predictors in longitudinal studies. Stat Med. 2007;26(6):1343–59.View ArticlePubMedGoogle Scholar
- US Centers for Disease Control and Prevention. About the National Health Interview Survey Hyattsville (MD): US Centers for Disease Control and Prevention; 2016 [updated October 8, 2016; cited 2016 February 19, 2016.]. Available from: http://www.cdc.gov/nchs/nhis/about_nhis.htm.
- Frome EL. The analysis of rates using Poisson regression models. Biometrics. 1983;39(3):665–74.View ArticlePubMedGoogle Scholar
- Zou G. A modified poisson regression approach to prospective studies with binary data. Am J Epidemiol. 2004;159(7):702–6.View ArticlePubMedGoogle Scholar
- Carter D, Signorino C. Back to the future: Modeling time dependence in binary data. Polit Anal. 2010;18(3):271–92.View ArticleGoogle Scholar
- Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2009.View ArticleGoogle Scholar
- Reid N. A Conversation with Sir David Cox. Stat Sci. 1994;9(3):17.View ArticleGoogle Scholar
- Carstensen B, Kristensen JK, Ottosen P, Borch-Johnsen K. Steering Group of the National Diabetes R. The Danish National Diabetes Register: trends in incidence, prevalence and mortality. Diabetologia. 2008;51(12):2187–96.View ArticlePubMedGoogle Scholar
- O’Brien RM. Age period cohort characteristic models. Soc Sci Res. 2000;29(1):123–39. PMID: WOS:000085452600006.View ArticleGoogle Scholar
- Bell FC, Miller ML. Life tables for the United States social security area 1900–1200 2005 [cited 2016 11/07]. Available from: https://www.ssa.gov/oact/NOTES/as120/LifeTables_Body.html.
- Vandeschrick C. The Lexis diagram, a misnomer. Demogr Res. 2001;4(3):97–124.View ArticleGoogle Scholar
- Keyes KM, Utz RL, Robinson W, Li G. What is a cohort effect? Comparison of three statistical methods for modeling cohort effects in obesity prevalence in the United States, 1971–2006. Soc Sci Med. 2010;70(7):1100–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Efron B. Logistic Regression, Survival Analysis, and the Kaplan-Meier Curve. J Am Stat Assoc. 1988;83(402):414–25.View ArticleGoogle Scholar
- Breslow NE, Lubin JH, Marek P, Langholz B. Multiplicative models and cohort analysis. J Am Stat Assoc. 1983;78(381):1–12.View ArticleGoogle Scholar
- Jackson JM, DeFor TA, Crain AL, Kerby TJ, Strayer LS, Lewis CE, et al. Validity of diabetes self-reports in the Women’s Health Initiative. Menopause. 2014;21(8):861–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Centers for Disease Control Prevention. National diabetes statistics report: estimates of diabetes and its burden in the United States, 2014: HHS CDC; 2014 [cited 2016 11/08]. Available from: http://www.cdc.gov/diabetes/pubs/statsreport14/national-diabetes-report-web.pdf.
- Gregg EW, Chen H, Wagenknecht LE, Clark JM, Delahanty LM, Bantle J, et al. Association of an intensive lifestyle intervention with remission of type 2 diabetes. JAMA. 2012;308(23):2489–96.View ArticlePubMedPubMed CentralGoogle Scholar