Early Enrollment in Diabetes Pay-for-Performance Program Reduced Loss of Life Expectancy in Newly-Diagnosed Patients with Type 2 Diabetes Mellitus

Article information

Diabetes Metab J. 2025;.dmj.2024.0507
Publication date (electronic) : 2025 March 26
doi : https://doi.org/10.4093/dmj.2024.0507
1Department of Public Health, College of Medicine, National Cheng Kung University, Tainan, Taiwan
2Department of Healthcare Administration and Medical Informatics, Kaohsiung Medical University, Kaohsiung, Taiwan
3Department of Statistics and Information Science, Fu Jen Catholic University, New Taipei City, Taiwan
4College of Medicine, Fu Jen Catholic University, New Taipei, Taiwan
5College of Medicine, National Taiwan University, Taipei, Taiwan
6Department of Occupational and Environmental Medicine, National Cheng Kung University Hospital, College of Medicine, National Cheng Kung University, Tainan, Taiwan
Corresponding authors: Li-Jung Elizabeth Ku https://orcid.org/0000-0002-4644-8439 Department of Public Health, College of Medicine, National Cheng Kung University, No. 1, University Road, Tainan 701, Taiwan E-mail: eljku@mail.ncku.edu.tw
Jung-Der Wang https://orcid.org/0000-0002-3176-4500 Department of Public Health, College of Medicine, National Cheng Kung University, No. 1, University Road, Tainan 701, Taiwan E-mail: jdwang121@gmail.com
Received 2024 August 25; Accepted 2024 December 3.

Abstract

Background

Diabetes is associated with reduced lifespan. To explore pay-for-performance (P4P) program and life expectancy (LE), we investigated the impact of interval between diabetes diagnosis and enrollment in P4P program on loss-of-LE among patients with diabetes in Taiwan.

Methods

From diabetes mellitus health database, which collected all newly-diagnosed patients with diabetes by calendar year, we selected patients, aged 40 to 64, with 503,662 in P4P group and 450,071 in non-P4P group, respectively, from 2004 to 2015, and followed them until the end of 2018 using Kaplan–Meier survival analysis. We simulated age-, gender-, and calendar yearmatched referents for each group through Monte Carlo method from Taiwan’s vital statistics. We constructed a restricted cubic spline model on logit-transformed relative survival between each group and its corresponding matched referents, and applied a rolling-over algorithm month-by-month to extrapolate the survival function of each index group to lifetime to estimate the LE, which was subtracted from that of matched referents to obtain the loss-of-LE.

Results

We found stratified analysis by interval showed that the earlier the enrollment, the lower the loss-of-LE, namely, 0.06±0.72 years for interval <1 year, 0.05±0.59 years for interval 1–4 years, 10.01±0.11 years for interval 5–9 years, and 12.77±0.14 years for interval 10–15 years, respectively (P<0.001), compared with 2.60±0.14 years for non-P4P group.

Conclusion

Early enrollment in the P4P program was associated with reduced loss-of-LE, indicating P4P might gain life if implemented early after diabetes diagnosis.

GRAPHICAL ABSTRACT

Highlights

• This study explores if early entry in P4P program extends lifespan in T2DM patients.

• We compared the index cohort with matched referents to estimate LE loss difference.

• Stratified analysis showed earlier P4P enrollment reduced loss of LE.

• Early P4P entry was linked to shorter duration of DM and timely insulin initiation.

INTRODUCTION

Diabetes is a disease of increasing prevalence with 537 million worldwide and 90 million in South-East Asia in 2021 [1]. Years spent with diabetes increased by 156% in men and 70% in women in the USA during 1985–2011 [2]. Type 2 diabetes mellitus was associated with an average 6-year loss of life expectancy (LE) [3], more among whites than among South Asians [4].

Early diagnosis and continuing multidisciplinary care with regular monitoring and timely intervention are the key elements in reducing the disease burden of diabetes [5]. Intensive glucose control starting at the time of diagnosis was associated with a significantly decreased risk of death from any cause [6], while intensive glucose lowering in patients with 10-year median duration of diabetes might increase mortality [7]. It appeared that time interval between diagnosis and target-driven treatment could predict survival.

In August, 1996, a pilot study of a modified “Diabetes Shared Care Network” was launched at Yilan county, Taiwan, and later expanded to three other counties; in 2001, Taiwan implemented a pay-for-performance (P4P) program islandwide (Guaynabo, Puerto Rico) [8], providing multidisciplinary team care to prevent diabetes complications. In a systematic review of the outcome effect of P4P program on diabetes in single-payer national health systems [9], it was associated with a lower risk of mortality in some studies of Taiwan [8,10,11], in one study [12] but not in the other study [13] of England or in one study [14] of Canada. In South Korea, the achievement of the national diabetes quality assessment program indicators was associated with a decreased risk of all-cause mortality [15].

Rarely has any study been done to investigate the impact of P4P program on LE in patients with diabetes. To control potential confounding by different age, gender, and calendar year of care technology for diabetes, we employed age, gender, and calendar year-matched referents simulated from vital statistics of general population to estimate the loss-of-LE as the indicator to determine if early enrollment in P4P program could prolong the survival in patients with diabetes. Namely, we aimed to (1) compare years of life lost associated with the interval between diabetes diagnosis and enrollment in P4P program and (2) quantify these relationships in subgroups of different intervals.

METHODS

Data source and study cohort

Data source

Taiwan introduced a compulsory, single-payer National Health Insurance (NHI) program with universal coverage in 1995. By 2009, approximately 99% of Taiwan’s residents had been enrolled in the NHI program [16]. From the original NHI database, the Ministry of Health and Welfare has provided two main cohort categories: population-based and disease-specific. The latter aimed to facilitate the rapid identification of patients with specific disease. We used the disease-specific diabetes mellitus health database, where patients with diabetes, by ICD-9-CM code 250 or ICD-10-CM code E11, were newly diagnosed based on at least one hospital admission or three or more outpatient visits within a calendar year from 2004 to 2018 [17]. They were interlinked with two other databases, the Taiwan Mortality Registry (2004–2018) and the Taiwan National Vital Statistics (2004–2018), as presented in Fig. 1. The study was conducted after obtaining approval by the Human Research Ethics Committee of National Cheng Kung University (HREC [Exempt] No. 111-493), and the requirement for informed consent was waived.

Fig. 1.

Flow diagram of establishment of subcohorts and estimation for life expectancy (LE) and loss-of-LE. T2DM, type 2 diabetes mellitus; ICD-9 or 10-CM, International Classification of Diseases, 9th or 10th Revision, Clinical Modification; P4P, pay-forperformance. aInterval between diabetes diagnosis and enrollment in P4P program.

Study cohort

Fig. 1 illustrates the process of establishing the study cohort. Patients with type 1 diabetes mellitus were excluded from this study. Derived from the sum of newly-diagnosed patients with diabetes, our study population began with 2,040,597 patients with diabetes during 2001–2015. The exclusion criteria were as follows: (1) Diabetes was diagnosed before the end of 2003 (n=315,359); (2) those with missing data in gender (n= 59,164); (3) those with age <40 or ≥65 years (n=651,348); (4) those with death date before diabetes diagnosis (n=3); (5) those with any catastrophic illness before diabetes diagnosis (n=60,990). Among the remaining 953,733 patients, 503,662 were enrolled in the P4P program (internal code in the NHI system: P14xx) between January 1, 2004 and December 31, 2015. As we recruited patients with diabetes into this cohort by calendar year during 2004–2015 and followed them until 2018, patients diagnosed in January 2004 could be followed for 15 years, whereas those diagnosed in December 2015 might be followed for 3 years only.

From mortality registry to loss-of-LE (‘Statistical analysis’ section), we further stratified the P4P group into four subgroups according to the time interval between diabetes diagnosis and enrollment in the P4P program: interval <1, 1–4, 5–9, and 10–15 years. The interval was calculated for an individual from the date of newly-diagnosed diabetes to the date of enrollment in P4P program confirmed by the first appearance of code P1407C in the claim data.

Availability of data and materials

The data that support the findings of this study are available from the Health and Welfare Data Science Center of Taiwan, but restrictions apply to the availability of these data, which were applied to be used exclusively for the current study, and so are not publicly available.

Measures

P4P program

Taiwan’s diabetes P4P program can be summarized according to the enrollment process, enrollment rate, four criteria of performance measure within a calendar year, and data uploading once a quarter. First, only those diagnosed with diabetes (International Classifications of Diseases, 9th Revision, Clinical Modifications [ICD-9-CM] code 250 or ICD-10-CM code E11) by the same doctor in a healthcare facility at least two times within the past 90 days could be enrolled in the P4P program. The primary diagnosis should be diabetes mellitus on the date of enrollment. The enrollment rate has consistently increased from 30% in 2010, 43% in 2015, 50% in 2020 [16] to 60% in 2022. The four criteria of performance include (1) the proportion of the enrollees who had completed four visits, namely once a quarter, within a calendar year; (2) the proportion of enrollees who had glycosylated hemoglobin (HbA1c) <7.0%; (3) the proportion of enrollees who had HbA1c >9.5%; (4) the proportion of enrollees who had low-density lipoprotein (LDL) cholesterol >130 mg/dL at annual evaluation visit. The physicians who ranked in the top 25% would be honored with the ‘Quality Excellence Award’ and rewarded with an annual incentive of Taiwan dollar 1,000 (United States dollar $30) per enrollee. At enrollment visit and annual visit, a physician is required to order exams including body weight, waist circumference, blood pressure, LDL, triglyceride, blood sugar, HbA1c, serum creatinine, urine microalbumin/creatinine ratio, retina and foot examination and upload the results to Health Promotion Administration to qualify for financial incentives. At the quarterly visits in between, the required exams included body weight, waist circumference, blood pressure, LDL, triglyceride, blood sugar, and HbA1c.

Based on the ‘ambulatory care and expenditures by visits’ file, patients with a reimbursement code of ‘P1401C,’ ‘P1407C’ (first enrollment in P4P program in a healthcare facility), ‘P1402C,’ ‘P1403C,’ ‘P1408C,’ ‘P1409C,’ ‘P1410C,’ or ‘P1411C’ in their physician’s order are judged to have been enrolled in diabetes P4P program in Supplementary Methods and Supplementary Tables 1-3.

Comorbidities

Major comorbidities at baseline were defined as those with at least one admission record or at least two outpatient visits for a certain diagnosis within 3 years before the date of diabetes diagnosis (index date). The comorbidities included hypertension (ICD-9 401–402 and 405), hyperlipidemia (ICD-9 272.0–272.4), coronary artery disease (ICD-9 414.8 and 414.9), heart failure (ICD-9 398.91, 402.01, 402.11, 402.91, 404.01, 404.03, 404.11, 404.13, 404.91, 404.93, and 428–428.9), stroke/cerebrovascular disease (ICD-9 430–438), peripheral vascular disease (ICD-9 440–443, 447, and 557), renal disease (ICD-9 016.0, 095.4, 189.0, 189.9, 223.0, 236.91, 250.4, 271.4, 274.1, 283.11, 403.X1, 404.X2, 404.X3, 440.1, 442.1, 447.3, 572.4, 580–588, 591, 642.1, 646.2, 753.12–753.17, 753.19, 753.2, and 794.4), chronic obstructive pulmonary disease (ICD-9 491–494, 496, and 510), liver disease (ICD-9 570, 571, and 572.4), rheumatoid arthritis/collagen deficiency disease (ICD-9 701.0, 710.0–710.9, 714.0–714.9, 720.0–720.9, and 725), gastrointestinal bleeding (ICD-9 456.0–456.2, 530.7, 531–534, 569.84, 569.85, and 578), psychoses (ICD-9 295.00–298.9), depression (ICD-9 300.4, 301.12, 309.1, and 311), dementia (ICD-9 290, 296.2x, 296.3x, 291.1, 291.2, and 294), and cancer (ICD-9 140–239). The corresponding ICD-10 codes were used, beginning in 2016. 50 articles with validation findings of diagnosis codes and related algorithms for health outcomes in Taiwan, including cardiovascular diseases, stroke, renal impairment, malignancy and diabetes etc., have been reported with positive predictive value in the range of 80% to 99% [17].

Other covariates

Covariates analyzed in this study included gender, age, date of diabetes diagnosis, date of enrollment in P4P program, Charlson comorbidity index (CCI) [18], Diabetes Complications Severity Index (DCSI) [19], and catastrophic illnesses [20]. Age was divided into five groups, 40–44, 45–49, 50–54, 55–59, and 60–64 years, respectively. The date of diabetes diagnosis (index date) was defined as the date of the first visit for diabetes; the date of P4P status was defined as the date of first enrollment in the P4P program as evidenced by the presence of a specific code (internal code in the NHI system: P1407C) in the claim data. The general medical status before index date was assessed using a modified version of the CCI, which was the sum for 19 comorbid conditions [18]. DCSI included the following seven categories of complications: cardiovascular, nephropathy, retinopathy, peripheral vascular, neuropathy, cerebrovascular, and metabolic complications [19]. Compared with a simple count of complications, the DCSI performed slightly better and appeared to be a useful tool for predicting mortality and risk of hospitalization [19]. Certificate of catastrophic illness was issued when a patient was diagnosed with one of 30 categories of catastrophic illnesses [20]. The duration of diabetes mellitus was measured from the date of a first visit for diabetes to death or the end of 2018; the duration of enrollment in P4P program was measured from the date of first enrollment to the date of withdrawal from P4P program, death or the end of 2018. Insulin use was defined according the drug prescriptions recorded in the 6 months prior to death or the end of the study.

Statistical analysis

Estimation of life expectancy, loss of life expectancy, and life-year gained

By interlinking the Taiwan National Mortality Registry and the reimbursement database of NHI for 503,662 subjects in P4P group and 450,071 subjects in non-P4P group, we applied the Kaplan–Meier method to estimate the survival function in the follow-up period of 15 years, namely 2004–2018, and then extrapolated to lifetime by a semiparametric method [21]. Namely, we applied Monte Carlo methods to obtain the lifetime survival function from age-, gender- and calendar yearmatched referents simulated from the National Vital Statistics database. With premature mortality caused by diabetes, the relative survival between each group and its corresponding reference population ranged between 0 and 1. Next, we performed logit transformations of the relative survivals obtained from the above methods and constructed a restricted cubic splines model to assure linear beyond the last knot. We used the fitted model to extrapolate the first month beyond 15-year follow-up, which was generally accurate and was treated as a new ‘observed value.’ Month-by-month, we repeatedly updated the data of the logit-transformed function by dropping the first month and adding the new ‘observed value’ as the last value, and re-fitted the restricted cubic splines model to obtain updated model for predicting a new value in the immediate next month. Through the above month-by-month rolling-over algorithm, we extrapolated the survival of the index group down until below 1% to stop (Fig. 2). Thereby, we estimated the lifetime survival, and the area under the survival curve was the LE of the patients in each index subcohort. We subtracted the LE of each index cohort from the corresponding age-, gender- and calendar year-matched reference population to obtain the loss-of-LE. The standard errors of LE and loss-of-LE were generated through a bootstrap method by implementing the extrapolation process with data simulated by repeatedly sampling with replacement from the real data set 100 times. The rolling extrapolation algorithm of iSQoL2 R package (http://sites.stat.sinica.edu.tw/isqol) used in this study is based on the estimated survival function derived from the observed survival times of all patients with diabetes in non-P4P group, P4P group and four subgroups of P4P group, respectively, during the followup period [21]. The survival times of patients with diabetes who were still alive at the end of the follow-up period, namely the end of 2018, were treated as censored, regardless of various lengths of their follow-up. Specifically, the Kaplan–Meier method was employed to estimate the survival function from the survival data. A key advantage of the Kaplan–Meier curve is its ability to incorporate these censored data effectively.

Fig. 2.

This figure shows estimation of lifetime survival function in non-pay-for-performance (P4P) group (A), P4P group (B), and its four subgroups (C, D, E, F), compared with its corresponding reference population, respectively. We subtracted the life expectancy (LE) of each index cohort from that of its corresponding age-, gender- and calendar year-matched reference population to obtain the loss-of-LE. LE and loss-of-LE of each index cohort were shown numerically in its corresponding panel.

In addition, we used this longitudinal cohort data to determine the incidence of adverse outcomes, including myocardial infarction, stroke, and end-stage renal disease (ESRD) (‘Comorbidities’ part of ‘Measures’ section), and catastrophic illnesses, DCSI and insulin use (‘Other covariates’ part of ‘Measures’ section) until death or the end of 2018, which was presented in Tables 1 and 2. We performed descriptive analyses where the continuous variables were presented as means and standard deviations, and dichotomized variables as percentages. The differences between the two compared groups were tested by t-test for continuous variables and chi-square test for categorical variables. Test for a linear trend was also performed for stratifications of age strata and for strata of intervals preceding enrollment in P4P program after diabetes diagnosis. Descriptive statistical analyses and survival analyses were performed using the SAS software version 9.4 (SAS Institute Inc., Cary, NC, USA). The above rolling extrapolations to estimate LEs and loss-of-LEs were carried out by the iSQoL2 R package developed by Hwang and Wang [22].

Life expectancy and loss of life expectancy in patients with diabetes with catastrophic illness, comorbidities, and insulin use

Stratification by P4P status, interval between diabetes diagnosis and enrollment in P4P program and diabetes complications severity index

RESULTS

Table 3 compares the characteristics of the study cohort before diabetes diagnosis, including 503,662 subjects in P4P group and 450,071 subjects in non-P4P group; P4P group was stratified into four subgroups based on interval: 164,711 for interval less than 1 year, 163,313 for interval 1–4 years, 134,620 for interval 5–9 years, and 41,018 for interval 10–15 years, respectively. Mean age was 53.0 to 54.0 years, and 43.3%–45.7% were women across non-P4P group and P4P subgroups. Significant P value occurred among four subgroups of P4P group across comorbidities, except depression.

Distribution of study cohort by gender, age, Charlson comorbidity index, and comorbidities

In Table 1, when P4P group was stratified by interval between diabetes diagnosis and enrollment in P4P program, the result showed that loss-of-LE was 0.06±0.72 years for interval less 1 year, 0.05±0.59 years for interval 1–4 years, 10.01±0.11 years for interval 5–9 years, and 12.77±0.14 years for interval 10–15 years, respectively. When the difference in loss was measured between with and without comorbidity, namely catastrophic illness, stroke, ESRD, or myocardial infarction. There was significant difference (P<0.001) among the four subgroups of P4P group where the difference in loss was large in intervals less than 1 and 1–4 years, and became smaller and smaller from interval 5–9 years to interval 10–15 years, suggesting that someone did not have one comorbidity, while the other comorbidity or more could have developed, leading to large loss-of-LE. In terms of insulin use, subgroup of interval 5–9 years suffered greater loss-of-LE, namely, 2.34±1.96 vs. 0.02±0.48 years without insulin and 13.21±0.20 vs. 7.61±0.30 years with insulin, compared with subgroups of interval <5 years.

Table 2 summarizes the loss-of-LE in each subcohort after stratification by scores of DCSI. We found that P4P group enrolled within 5 years after diabetes diagnosis did show higher proportions of DCSI 1–2 (18.4% vs. 14.6%) with nearly no loss-of-LE, but lower proportion of DCSI ≥3 (67.4% vs. 76.5%) with less loss-of-LE (2.33 years vs. 10.43 years), compared with those enrolled later. Furthermore, those enrolled after 5 years of diabetes diagnosis showed longer diabetes duration accompanied with shorter P4P care. Moreover, a large loss-of-LE was associated with DCSI ranging from 0 to 2 in those patients.

Fig. 2 shows estimation of lifetime survival function in non-P4P group, P4P group and its four subgroups, compared with its corresponding reference population. Loss-of-LE was 2.60 years for non-P4P group, 3.55 years for all P4P group, 0.06 year for interval less than 1 year, 0.05 year for interval 1–4 years, 10.01 years for interval 5–9 years, and 12.77 years for interval 10–15 years, respectively.

Fig. 3 demonstrates the validation of survival function by comparing the result of extrapolation from the end of the 8th year to the 15th year with that of the actual survival rate based on Kaplan–Meier’s estimate throughout 15 years of follow-up. The relative bias assuming the Kaplan–Meier’s estimate as the gold standard was 2.5%.

Fig. 3.

This figure demonstrates the validation of survival function by comparing the result of extrapolation from the end of the 8th year to the end of the 15th year with that of the actual survival rate based on Kaplan–Meier’s (K-M) estimate throughout 15 years of follow-up. The relative bias assuming the K-M estimate as the gold standard was 2.5%.

In Supplementary Table 4, female showed a higher percentage (44.7% vs. 43.3%) in P4P group than in non-P4P group. The mean age of diabetes diagnosis was 53.6 years for P4P group and 53.7 years for non-P4P group, respectively. For CCI score, P4P group vs. non-P4P group is 53.2% vs. 51.9% for score 0, 31.3% vs. 31.5% for score 1, and 10.8% vs. 11.3% for score 2. Significant difference occurred across comorbidities except coronary artery disease, rheumatology disorders and psychosis.

In Supplementary Table 5, during 15-year follow-up, mean duration of diabetes was 8.16 years for non-P4P group and 9.16 years for P4P group; mean duration of P4P status was 5.34 years. P4P group (n=503,662) was associated with significantly more loss-of-LE (3.55±0.25 years vs. 2.60±0.14 years, P<0.001), compared with non-P4P group (n=450,071). In terms of difference in loss-of-LE, P4P group was associated with less difference in loss-of-LE in patients with catastrophic illness (P<0.001), but more difference in loss-of-LE in patients with myocardial infarction (P<0.001).

DISCUSSION

We found loss-of-LE was 0.06 year for interval less than 1 year, 0.05 year for interval 1–4 years, 10.01 years for interval 5–9 years, and 12.77 years for interval 10–15 years, respectively (P<0.001), compared with 2.60 years for the non-P4P group. Although our study seemed to provide evidence for the long-term or lifetime survival benefit for P4P program in diabetes care, our results of the survival benefit from P4P program could not be directly considered as causal. We have following arguments to corroborate the association between early enrollment in P4P program and reduced loss-of-LE. First, we used islandwide diabetes mellitus health database. This made our big study cohort highly representative of Taiwan’s whole diabetes population. The linkage with Taiwan Mortality Registry assured the comprehensiveness of their survival states. Second, high risk of all-cause mortality was associated with early age at diagnosis of diabetes [3] and intensive glucose lowering in patients with 10-year median duration of diabetes [7], while 10 years after the cessation of randomized intervention, a legacy effect was associated with intensive glucose control in patients with newly-diagnosed type 2 diabetes mellitus [6]. As diabetic complications increased with diabetes duration [23], elevated HbA1c [24], high blood pressure [25] and dyslipidemia [26], our findings (Table 1) on the probable impact of early enrollment in P4P program by stratification analysis corroborated these previous reports. Taiwan’s P4P program was target-driven multidisciplinary team care to empower patients for lifestyle modification [27], which corroborated the benefits of multifactorial risk reduction and improved survival as demonstrated in the Steno-2 study [28] with a median gain of 7.9 years in life [29]. Third, P4P subgroup of interval less than 1 year seemed to have a lower prevalence of comorbidities at baseline, suggesting an unintended consequence of P4P program, namely the exclusion of high-risk patients. Hsieh et al. [30] reported that patients who were male, older than 75 years, and those with DCSI scores or CCI scores larger than 2 were more likely to be excluded from enrollment, which could be gradually mitigated by increasing health literacy and growing enrollment rates. Therefore, early enrollment in Taiwan’s P4P program could achieve strict adherence, regular monitoring and timely intervention, and result in reduced loss-of-LE in the P4P group. Fourth, since we controlled for the potential confounding factors of gender, age, interval, major comorbidities (through catastrophic illness, stroke, ESRD, myocardial infarction, and DCSI), and insulin use in stratification analysis, the above factors could not explain the estimated difference in loss-of-LE between non-P4P and P4P subgroups (Table 1). Furthermore, in terms of comorbidities at baseline, non-P4P group and subgroup of interval 1–4 years had similar distributions of prevalence rates; among those involving target organ damages, namely coronary artery disease, heart failure, stroke, peripheral vascular disease and renal disease, the differences in prevalence across non-P4P group and four subgroups of P4P group were all less than or equal to 1.2% (Table 3), which corroborated the hypothesis of beneficial impact of early enrollment on LE. Fifth, in Table 1, the subgroup of interval 5–9 years had higher rates of catastrophic illnesses (17% vs.14%), stroke (10% vs. 6%), ESRD (3% vs. 2%), and myocardial infarction (4% vs. 3%), compared with subgroup of interval 1–4 years, and, to a greater degree, catastrophic illnesses (17% vs.11%), stroke (10% vs. 5%), ESRD (3% vs. 1%), and myocardial infarction (4% vs. 2%), compared with non-P4P group. Furthermore, the diabetes duration of subgroup of interval 5–9 years was 2 years longer than those of subgroup of interval 1–4 years and non-P4P group: 10.79±2.44 years vs. 8.05±3.08 and 8.16±3.50 years. Similarly, the difference in loss-of-LE among patients receiving insulin versus no insulin in subgroup of interval 5–9 years was higher than those of subgroup of interval 1–4 years and non-P4P group (10.87±1.98 years vs. 7.59±0.54 and 10.77±0.22 years). All these facts implied that early enrollment in the P4P program was associated with shorter duration of diabetes, timely insulin use [31,32], and fewer complications, corroborating Yoo et al.’s [33] observations. Sixth, as the validation of our extrapolation was demonstrated by an algorithm of repeatedly fitting restricted cubic spline models and rolling-over month-by-month from the end of the 8th year to the 15th year (Fig. 3) [34,35], we anticipated that the same method using 15 years of real data extrapolated to lifetime would be reasonably accurate. In fact, our approach of external additive hazards using the general population mortality as a splint for lifetime extrapolation was considered as relatively accurate in a recent review of statistical methods by van Oostrum et al. [36]. Thus, we tentatively concluded that the early enrollment (namely <5 years) in Taiwan’s P4P program of diabetes care was associated with a significant reduced loss-of-LE, which deserved further corroboration.

P4P program benefited the patients with diabetes with disability or depression. Among the disabled patients with diabetes, the risk of death was 0.41 times lower in those enrolled in the P4P program than in those who were not, suggesting that enrollment in the P4P program was beneficial in mitigating the risk of death [37]. Moreover, patients with diabetes with depressive symptoms enrolled in the diabetes P4P program had lower expenses for depression-related care than those not enrolled in the program, indicating a positive spillover effect of P4P program [38].

Lastly and importantly, stratification analysis of P4P subgroups revealed that loss-of-LE for subgroups of interval 5–9 and 10–15 years were 10.0 and 12.8 years, respectively, in contrast to almost no loss-of-LE in those who were enrolled in P4P within less than 5 years after diagnosis (Table 1). It corroborated that early diagnosis of diabetes and prompt enrollment in P4P program were essential to robustly implement target-driven multidisciplinary approach to ensure regular monitoring and timely intervention to reduce disease burden and prolong life in patients with diabetes.

There are some limitations in this study that must be acknowledged. First, this is not a randomized control trial; therefore, some unmeasured confounders cannot be ruled out despite stratification analysis. Second, limiting P4P enrollees to those aged 40 to 64 years, this study cannot be generalized to all patients with diabetes. Third, as there were no detailed data of HbA1c, blood pressure, LDL level, body mass index, cigarette smoking, physical activity, and educational level in claim data of NHI database, we were unable to make any additional pathophysiological inference. Future studies are warranted to dig into detailed mechanisms.

In conclusion, early enrollment in the P4P program was associated with reduced loss-of-LE, indicating P4P program might help patients with diabetes gain life if implemented early and faithfully after diabetes diagnosis.

SUPPLEMENTARY MATERIALS

Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2024.0507.

Supplementary Table 1.

Components of initial enrollment visit (package P1401C: 1,845 NTD)

dmj-2024-0507-Supplementary-Table-1.pdf
Supplementary Table 2.

Continuing care visit (package P1402C: 875 NTD): potential components of continuing care visit

dmj-2024-0507-Supplementary-Table-2.pdf
Supplementary Table 3.

Annual evaluation visit (package P1403C: 2,245 NTD): potential components of annual evaluation visit

dmj-2024-0507-Supplementary-Table-3.pdf
Supplementary Table 4.

Distribution of study cohort by gender, age, Charlson comorbidity index, and comorbidities

dmj-2024-0507-Supplementary-Table-4.pdf
Supplementary Table 5.

Life expectancy and loss of life expectancy in patients with diabetes with catastrophic illness, comorbidities, and insulin use

dmj-2024-0507-Supplementary-Table-5.pdf

Notes

CONFLICTS OF INTEREST

The diabetes P4P program in Taiwan is a national policy administered by the National Health Insurance to assure quality of care for patients with diabetes. Although one of the authors (Yu-Ching Chen) has patients who have joined the P4P program since 2001, there are no relevant conflicts of interest to disclose for all authors.

AUTHOR CONTRIBUTIONS

Conception or design: Y.C.C., J.D.W., L.J.E.K.

Acquisition, analysis, or interpretation of data: W.M.W., B.J.L., J.D.W., L.J.E.K.

Drafting the work or revising: all authors.

Final approval of the manuscript: all authors.

FUNDING

This work was partially funded by National Science and Technology Council (NSTC 112-2627-M-006-002).

ACKNOWLEDGMENTS

We are grateful to Health Data Science Center, National Cheng Kung University Hospital for providing administrative and technical support.

References

1. International Diabetes Federation. IDF Diabetes Atlas. 10th edth ed. Brussels: IDF; 2021.
2. Gregg EW, Zhuo X, Cheng YJ, Albright AL, Narayan KM, Thompson TJ. Trends in lifetime risk and years of life lost due to diabetes in the USA, 1985-2011: a modelling study. Lancet Diabetes Endocrinol 2014;2:867–74.
3. Emerging Risk Factors Collaboration. Life expectancy associated with different ages at diagnosis of type 2 diabetes in high-income countries: 23 million person-years of observation. Lancet Diabetes Endocrinol 2023;11:731–42.
4. Wright AK, Kontopantelis E, Emsley R, Buchan I, Sattar N, Rutter MK, et al. Life expectancy and cause-specific mortality in type 2 diabetes: a population-based cohort study quantifying relationships in ethnic subgroups. Diabetes Care 2017;40:338–45.
5. Chan JC, Lim LL, Wareham NJ, Shaw JE, Orchard TJ, Zhang P, et al. The Lancet commission on diabetes: using data to transform diabetes care and patient lives. Lancet 2021;396:2019–82.
6. Holman RR, Paul SK, Bethel MA, Matthews DR, Neil HA. 10-Year follow-up of intensive glucose control in type 2 diabetes. N Engl J Med 2008;359:1577–89.
7. Action to Control Cardiovascular Risk in Diabetes Study Group, Gerstein HC, Miller ME, Byington RP, Goff DC Jr, Bigger JT, et al. Effects of intensive glucose lowering in type 2 diabetes. N Engl J Med 2008;358:2545–59.
8. Chen YC, Lee CT, Lin BJ, Chang YY, Shi HY. Impact of pay-for-performance on mortality in diabetes patients in Taiwan: a population-based study. Medicine (Baltimore) 2016;95e4197.
9. Gupta N, Ayles HM. Effects of pay-for-performance for primary care physicians on diabetes outcomes in single-payer health systems: a systematic review. Eur J Health Econ 2019;20:1303–15.
10. Hsieh HM, Chiu HC, Lin YT, Shin SJ. A diabetes pay-for-performance program and the competing causes of death among cancer survivors with type 2 diabetes in Taiwan. Int J Qual Health Care 2017;29:512–20.
11. Kung FP, Tsai CF, Lu CL, Huang LC, Lu CH. Diabetes pay-for-performance program can reduce all-cause mortality in patients with newly diagnosed type 2 diabetes mellitus. Medicine (Baltimore) 2020;99e19139.
12. McKay AJ, Gunn LH, Vamos EP, Valabhji J, Molina G, Molokhia M, et al. Associations between attainment of incentivised primary care diabetes indicators and mortality in an English cohort. Diabetes Res Clin Pract 2021;174:108746.
13. Kontopantelis E, Springate DA, Ashworth M, Webb RT, Buchan IE, Doran T. Investigating the relationship between quality of primary care and premature mortality in England: a spatial whole-population study. BMJ 2015;350:h904.
14. Thavam T, Devlin RA, Thind A, Zaric GS, Sarma S. The impact of the diabetes management incentive on diabetes-related services: evidence from Ontario, Canada. Eur J Health Econ 2020;21:1279–93.
15. Huh JH, Kwon S, Kim GO, Kim BY, Ha KH, Kim DJ. Long-term effectiveness of the national diabetes quality assessment program in South Korea. Diabetes Care 2023;46:1700–6.
16. Chen TT, Oldenburg B, Hsueh YS. Chronic care model in the diabetes pay-for-performance program in Taiwan: benefits, challenges and future directions. World J Diabetes 2021;12:578–89.
17. Huang YT, Wei T, Huang YL, Wu YP, Chan KA. Validation of diagnosis codes in healthcare databases in Taiwan, a literature review. Pharmacoepidemiol Drug Saf 2023;32:795–811.
18. Elixhauser A, Steiner C, Harris DR, Coffey RM. Comorbidity measures for use with administrative data. Med Care 1998;36:8–27.
19. Young BA, Lin E, Von Korff M, Simon G, Ciechanowski P, Ludman EJ, et al. Diabetes complications severity index and risk of mortality, hospitalization, and healthcare utilization. Am J Manag Care 2008;14:15–23.
20. National Health Insurance Administration, Ministry of Health and Welfare. Catastrophic illnesses [Chinese]. Available from: https://www.nhi.gov.tw/ch/cp-6086-caf5f-2957-1.html (cited 2025 Feb 20).
21. Hwang JS, Hu TH, Lee LJ, Wang JD. Estimating lifetime medical costs from censored claims data. Health Econ 2017;26:e332–44.
22. Hwang JS, Wang JD. Integrating health profile with survival for quality of life assessment. Qual Life Res 2004;13:1–14.
23. Zoungas S, Woodward M, Li Q, Cooper ME, Hamet P, Harrap S, et al. Impact of age, age at diagnosis and duration of diabetes on the risk of macrovascular and microvascular complications and death in type 2 diabetes. Diabetologia 2014;57:2465–74.
24. Stratton IM, Adler AI, Neil HA, Matthews DR, Manley SE, Cull CA, et al. Association of glycaemia with macrovascular and microvascular complications of type 2 diabetes (UKPDS 35): prospective observational study. BMJ 2000;321:405–12.
25. UK Prospective Diabetes Study Group. Tight blood pressure control and risk of macrovascular and microvascular complications in type 2 diabetes: UKPDS 38. BMJ 1998;317:703–13.
26. Colhoun HM, Betteridge DJ, Durrington PN, Hitman GA, Neil HA, Livingstone SJ, et al. Primary prevention of cardiovascular disease with atorvastatin in type 2 diabetes in the Collaborative Atorvastatin Diabetes Study (CARDS): multicentre randomised placebo-controlled trial. Lancet 2004;364:685–96.
27. Lee YJ, Shin SJ, Wang RH, Lin KD, Lee YL, Wang YH. Pathways of empowerment perceptions, health literacy, self-efficacy, and self-care behaviors to glycemic control in patients with type 2 diabetes mellitus. Patient Educ Couns 2016;99:287–94.
28. Gaede P, Lund-Andersen H, Parving HH, Pedersen O. Effect of a multifactorial intervention on mortality in type 2 diabetes. N Engl J Med 2008;358:580–91.
29. Gaede P, Oellgaard J, Carstensen B, Rossing P, Lund-Andersen H, Parving HH, et al. Years of life gained by multifactorial intervention in patients with type 2 diabetes mellitus and microalbuminuria: 21 years follow-up on the Steno-2 randomised trial. Diabetologia 2016;59:2298–307.
30. Hsieh HM, Tsai SL, Mau LW, Chiu HC. Effects of changes in diabetes pay-for-performance incentive designs on patient risk selection. Health Serv Res 2016;51:667–86.
31. Tsai ST, Pathan F, Ji L, Yeung VT, Chadha M, Suastika K, et al. First insulinization with basal insulin in patients with type 2 diabetes in a real-world setting in Asia. J Diabetes 2011;3:208–16.
32. Chan WB, Chen JF, Goh SY, Vu TT, Isip-Tan IT, Mudjanarko SW, et al. Challenges and unmet needs in basal insulin therapy: lessons from the Asian experience. Diabetes Metab Syndr Obes 2017;10:521–32.
33. Yoo H, Choo E, Lee S. Study of hospitalization and mortality in Korean diabetic patients using the diabetes complications severity index. BMC Endocr Disord 2020;20:122.
34. Liao CT, Hsieh TH, Shih CY, Liu PY, Wang JD. Cost-effectiveness of percutaneous coronary intervention versus medical therapy in patients with acute myocardial infarction: real-world and lifetime-horizon data from Taiwan. Sci Rep 2021;11:5608.
35. Lien WC, Wang WM, Wang F, Wang JD. Savings of loss-of-life expectancy and lifetime medical costs from prevention of spinal cord injuries: analysis of nationwide data followed for 17 years. Inj Prev 2021;27:567–73.
36. van Oostrum I, Ouwens M, Remiro-Azocar A, Baio G, Postma MJ, Buskens E, et al. Comparison of parametric survival extrapolation approaches incorporating general population mortality for adequate health technology assessment of new oncology drugs. Value Health 2021;24:1294–301.
37. Kuo WY, Tsai WC, Kung PT. Participation and outcomes among disabled and non-disabled people in the diabetes pay-for-performance program. Healthcare (Basel) 2023;11:2742.
38. Sung MC, Chung KP, Cheng SH. Impact of a diabetes pay-for-performance program on nonincentivized mental disorders: a panel study based on claims database analysis. Cost Eff Resour Alloc 2023;21:41.

Article information Continued

Fig. 1.

Flow diagram of establishment of subcohorts and estimation for life expectancy (LE) and loss-of-LE. T2DM, type 2 diabetes mellitus; ICD-9 or 10-CM, International Classification of Diseases, 9th or 10th Revision, Clinical Modification; P4P, pay-forperformance. aInterval between diabetes diagnosis and enrollment in P4P program.

Fig. 2.

This figure shows estimation of lifetime survival function in non-pay-for-performance (P4P) group (A), P4P group (B), and its four subgroups (C, D, E, F), compared with its corresponding reference population, respectively. We subtracted the life expectancy (LE) of each index cohort from that of its corresponding age-, gender- and calendar year-matched reference population to obtain the loss-of-LE. LE and loss-of-LE of each index cohort were shown numerically in its corresponding panel.

Fig. 3.

This figure demonstrates the validation of survival function by comparing the result of extrapolation from the end of the 8th year to the end of the 15th year with that of the actual survival rate based on Kaplan–Meier’s (K-M) estimate throughout 15 years of follow-up. The relative bias assuming the K-M estimate as the gold standard was 2.5%.

Table 1.

Life expectancy and loss of life expectancy in patients with diabetes with catastrophic illness, comorbidities, and insulin use

Study cohort Non-P4P group (n=450,071) P4P group (n=503,662)a
P value
Interval <1 year (n=164,711) Interval 1–4 years (n=163,313) Interval 5–9 years (n=134,620) Interval 10–15 years (n=41,018)
Female sex 194,837 (43.3) 72,597 (44.1) 73,786 (45.2) 60,227 (44.7) 18,746 (45.7)
Age, mean±SD, yr 53.7±6.6 53.5±6.6 54.0±6.5 53.6±6.4 53.0±6.5
Duration of diabetes, mean±SD, yr 8.16±3.50 7.86±3.55 8.05±3.08 10.79±2.44 13.44±1.43 <0.001
Duration of P4P status, mean±SD, yr 0±0 7.69±3.56 5.25±3.18 3.62±2.24 1.93±1.44 <0.001
LE, mean±SE, yr 26.44±0.14 29.79±0.51 28.87±0.59 19.20±0.10 16.84±0.12 <0.001
Loss-of-LE, mean±SE, yr 2.60±0.14 0.06±0.72 0.05±0.59 10.01±0.11 12.77±0.14 <0.001
Catastrophic illness 51,686 (11.0) 18,251 (11.0) 22,349 (14.0) 23,460 (17.0) 7,043 (17.0)
 Loss-of-LE, mean±SE
  Without disease (%) 0.11±0.04 (89.0) 0.11±0.77 (89.0) 0.05±0.69 (86.0) 6.82±0.24 (83.0) 12.43±0.19 (83.0)
  With disease (%) 16.14±0.12 (11.0) 11.76±0.37 (11.0) 12.05±0.40 (14.0) 13.97±0.07 (17.0) 13.73±0.16 (17.0)
 Difference in loss, mean±SEb
  With and without 16.02±0.12 11.65±0.78 12.00±0.75 7.15±0.25 1.30±0.27 <0.001
Stroke 20,890 (5.0) 7,998 (5.0) 10,264 (6.0) 12,801 (10.0) 4,202 (10.0)
 Loss-of-LE, mean±SE
  Without disease (%) 2.41±0.12 (95.0) 0.01±0.59 (95.0) 0.23±0.73 (94.0) 9.53±0.14 (90.0) 12.70±0.13 (90.0)
  With disease (%) 9.49±0.26 (5.0) 9.20±0.67 (5.0) 8.63±0.54 (6.0) 11.91±0.12 (10.0) 12.04±0.32 (10.0)
 Difference in loss, mean±SEb
  With and without 7.08±0.30 9.19±0.92 8.39±0.78 2.38±0.19 –0.66±0.36 <0.001
ESRD 5,469 (1.0) 2,656 (2.0) 3,150 (2.0) 4,197 (3.0) 1,394 (3.0)
 Loss-of-LE, mean±SE
  Without disease (%) 3.02±0.16 (99.0) 0.26±0.83 (98.0) 0.03±0.98 (98.0) 9.69±0.13 (97.0) 12.48±0.15 (97.0)
  With disease (%) 14.04±0.48 (1.0) 13.10±0.64 (2.0) 13.71±0.49 (2.0) 14.82±0.16 (3.0) 15.12±0.21 (3.0)
 Difference in loss, mean±SEb
  With and without 11.02±0.49 12.84±0.88 13.68±1.12 5.13±0.21 2.65±0.26 <0.001
Myocardial infarction 7,778 (2.0) 3,275 (2.0) 4,358 (3.0) 5,183 (4.0) 1,721 (4.0)
 Loss-of-LE, mean±SE
  Without disease (%) 3.10±0.15 (98.0) 0.01±0.49 (98.0) 0.04±0.71 (97.0) 9.93±0.14 (96.0) 12.63±0.15 (96.0)
  With disease (%) 10.00±0.46 (2.0) 12.06±0.61 (2.0) 10.02±0.82 (3.0) 12.74±0.17 (4.0) 13.13±0.33 (4.0)
 Difference in loss, mean±SEb
  With and without 6.90±0.49 12.05±0.79 9.98±0.97 2.81±0.22 0.50±0.35 <0.001
Insulin use 85,149 (19.0) 59,744 (36.0) 63,941 (39.0) 64,384 (48.0) 20,495 (50.0)
 Loss-of-LE, mean±SE
  Without insulin (%) 0.01±0.90 (81.0) 0.01±1.89 (64.0) 0.02±0.48 (61.0) 2.34±1.96 (52.0) 11.62±0.61 (50.0)
  With insulin (%) 10.78±0.20 (19.0) 8.56±0.51 (36.0) 7.61±0.30 (39.0) 13.21±0.20 (48.0) 14.47±0.64 (50.0)
 Difference in loss, mean±SEb
  With and without 10.77±0.22 8.55±0.53 7.59±0.54 10.87±1.98 2.85±0.85 <0.001

Values are presented as number (%) unless otherwise indicated. P value was measured for comparison of difference in loss-of-LE among four subgroups of P4P group.

P4P, pay-for-performance; SD, standard deviation; LE, life expectancy; SE, standard error; ESRD, end-stage renal disease.

a

Interval between diabetes diagnosis and enrollment in P4P program,

b

Difference in loss=(loss-of-LE±SE with disease or insulin)–(loss-of-LE±SE without disease or insulin).

Table 2.

Stratification by P4P status, interval between diabetes diagnosis and enrollment in P4P program and diabetes complications severity index

Study cohort No. (%) (n=953,733) Age, mean±SD, yr Duration
LE, mean±SE, yr Loss-of-LE, mean±SE, yr P value
Diabetes, mean±SD, yr P4P status, mean±SD, yr
Non-P4P group 450,071 <0.001
 DCSI 0 135,061 (30.0) 52.8±6.6 7.04±3.23 0±0 28.23±0.15 1.58±0.15
 DCSI 1–2 94,950 (21.1) 53.2±6.5 8.15±3.28 0±0 27.57±0.78 1.98±0.78
 DCSI ≥3 220,060 (48.9) 54.4±6.5 8.84±3.57 0±0 21.97±0.27 6.47±0.27
P4P group 503,662 <0.001
 DCSI 0 65,550 (13.0) 52.7±6.6 7.19±3.13 4.41±3.01 30.04±0.50 0.08±0.50
 DCSI 1–2 88,359 (17.6) 52.8±6.5 8.24±3.29 4.99±3.32 30.02±0.60 0.11±0.62
 DCSI ≥3 349,753 (69.4) 54.0±6.5 9.76±3.41 5.60±3.64 20.79±0.16 8.03±0.16
Interval <1 yeara 164,711 <0.001
 DCSI 0 27,652 (16.8) 52.6±6.7 6.20±2.85 6.04±2.88 29.85±0.04 0.10±0.01
 DCSI 1–2 33,468 (20.3) 52.8±6.6 7.17±3.16 7.01±3.17 29.44±1.34 0.44±1.33
 DCSI ≥3 103,591 (62.9) 53.9±6.5 8.52±3.66 8.35±3.66 27.63±0.49 1.25±0.48
Interval 1–4 yearsa 163,313 <0.001
 DCSI 0 23,218 (14.2) 53.2±6.6 6.52±2.52 3.77±2.61 29.29±0.31 0.11±0.51
 DCSI 1–2 30,043 (18.4) 53.3±6.5 7.37±2.83 4.60±2.94 29.55±1.21 0.18±1.21
 DCSI ≥3 110,052 (67.4) 54.3±6.5 8.56±3.12 5.73±3.23 26.30±0.47 2.33±0.47
Interval 5–9 yearsa 134,620 <0.001
 DCSI 0 12,028 (8.9) 52.2±6.4 9.47±2.49 2.58±2.14 21.82±1.12 8.43±1.12
 DCSI 1–2 19,687 (14.6) 52.5±6.4 10.09±2.44 3.06±2.21 22.84±1.14 7.34±1.09
 DCSI ≥3 102,905 (76.5) 53.9±6.4 11.07±2.35 3.85±2.20 18.41±0.17 10.43±0.18
Interval 10–15 yearsa 41,018 <0.001
 DCSI 0 2,652 (6.5) 51.3±6.3 12.95±1.53 1.41±1.39 16.52±0.56 14.47±0.56
 DCSI 1–2 5,161 (12.6) 51.5±6.3 13.11±1.50 1.58±1.42 16.97±0.33 14.03±0.34
 DCSI ≥3 33,205 (81.0) 53.4±6.5 13.53±1.39 2.03±1.43 13.01±0.10 15.29±0.11

P value was measured for comparison of loss-of-LE among three subgroups based on DCSI in non-P4P group, P4P group and its four subgroups, respectively.

P4P, pay-for-performance; SD, standard deviation; LE, life expectancy; SE, standard error; DCSI, diabetes complications severity index.

a

Interval between diabetes diagnosis and enrollment in P4P program.

Table 3.

Distribution of study cohort by gender, age, Charlson comorbidity index, and comorbidities

Study cohort Non-P4P group (n=450,071) P4P group (n=503,662)a
P value
Interval <1 year (n=164,711) Interval 1–4 years (n=163,313) Interval 5–9 years (n=134,620) Interval 10–15 years (n= 41,018)
Sex <0.001
 Female 194,837 (43.3) 72,597 (44.1) 73,786 (45.2) 60,227 (44.7) 18,746 (45.7)
 Male 255,234 (56.7) 92,114 (55.9) 89,527 (54.8) 74,393 (55.3) 22,272 (54.3)
Age, yr 53.7±6.6 53.5±6.6 54.0±6.5 53.6±6.4 53.0±6.5 <0.001
 40–44 55,232 (12.3) 21,096 (12.8) 18,352 (11.2) 16,078 (11.9) 5,483 (13.4)
 45–49 81,869 (18.2) 30,813 (18.7) 28,635 (17.5) 25,379 (18.9) 8,432 (20.6)
 50–54 108,766 (24.2) 40,294 (24.5) 39,321 (24.1) 33,431 (24.8) 10,721 (26.1)
 55–59 112,098 (24.9) 39,971 (24.3) 41,882 (25.6) 34,034 (25.3) 9,203 (22.4)
 60–64 92,106 (20.5) 32,537 (19.8) 35,123 (21.5) 25,698 (19.1) 7,179 (17.5)
Charlson comorbidity index <0.001
 0 233,428 (51.9) 92,869 (56.4) 86,847 (53.2) 66,964 (49.7) 21,363 (52.1)
 1 141,947 (31.5) 47,710 (29.0) 51,004 (31.2) 45,516 (33.8) 13,565 (33.1)
 2 50,935 (11.3) 16,399 (10.0) 17,783 (10.9) 15,781 (11.7) 4,399 (10.7)
 3 15,114 (3.4) 4,924 (3.0) 5,147 (3.2) 4,297 (3.2) 1,169 (2.9)
 4 5,295 (1.2) 1,771 (1.1) 1,655 (1.0) 1,375 (1.0) 399 (1.0)
 ≥5 3,352 (0.7) 1,038 (0.6) 877 (0.5) 687 (0.5) 123 (0.3)
Comorbidities
 Hypertension 169,548 (37.7) 52,313 (31.8) 61,927 (37.9) 55,906 (41.5) 16,099 (39.3) <0.001
 Hyperlipidemia 113,487 (25.2) 32,450 (19.7) 39,376 (24.1) 32,094 (23.8) 8,826 (21.5) <0.001
 Coronary artery disease 14,455 (3.2) 4,254 (2.6) 5,250 (3.2) 5,142 (3.8) 1,550 (3.8) <0.001
 Heart failure 9,678 (2.2) 3,177 (1.9) 3,393 (2.1) 2,945 (2.2) 802 (2.0) 0.001
 Stroke 20,401 (4.5) 6,733 (4.1) 6,551 (4.0) 6,140 (4.6) 1,645 (4.0) <0.001
 Peripheral vascular disease 4,784 (1.1) 1,502 (0.9) 1,639 (1.0) 1,423 (1.1) 410 (1.0) <0.001
 Renal disease 14,032 (3.1) 4,873 (3.0) 4,489 (2.8) 3,272 (2.4) 873 (2.1) <0.001
 COPD 18,021 (4.0) 6,607 (4.0) 7,069 (4.3) 6,208 (4.6) 1,843 (4.5) <0.001
 Liver disease 45,352 (10.1) 13,334 (8.1) 14,959 (9.2) 13,540 (10.1) 4,276 (10.4) <0.001
 Rheumatoid arthritis 5,100 (1.1) 1,698 (1.0) 1,910 (1.2) 1,554 (1.2) 489 (1.2) <0.001
 Gastrointestinal bleeding 38,363 (8.5) 12,344 (7.5) 13,538 (8.3) 11,374 (8.5) 3,334 (8.1) <0.001
 Psychosis 4,939 (1.1) 1,853 (1.1) 1,973 (1.2) 1,399 (1.0) 366 (0.9) <0.001
 Depression 8,081 (1.8) 2,703 (1.6) 3,036 (1.9) 2,275 (1.7) 592 (1.4) 0.134
 Dementia 3,448 (0.8) 1,211 (0.7) 1,328 (0.8) 906 (0.7) 212 (0.5) <0.001

Values are presented as number (%) or mean±standard deviation. P value was measured for comparison of each variable among four subgroups of P4P group.

P4P, pay-for-performance; COPD, chronic obstructive pulmonary disease.

a

Interval between diabetes diagnosis and enrollment in P4P program.