Continuity of care with physicians and risk of subsequent hospitalization and end-stage renal disease in newly diagnosed type 2 diabetes mellitus patients.

PURPOSE: Effective management for type 2 diabetes mellitus (DM) can slow the progression of kidney outcomes and reduce hospital admissions. Better continuity of care (COC) was found to improve patients' adherence and self-management. This study examined the associations between COC, hospitalization, and end-stage renal disease (ESRD) in DM patients. PATIENTS AND METHODS: In the cohort study, data from 1996 to 2012 were retrieved from the Longitudinal Health Insurance Database, using inverse probability weighted analysis. A total of 26,063 patients with newly diagnosed type 2 DM who had been treated with antihyperglycemic agents were included. COC is to assess the extent to which a DM patient visited the same physician during the study period. This study categorized COC into 3 groups - low, intermediate, and high, - according to the distribution of scores in our sample.
RESULTS: The number of ESRD patients in the high, intermediate, and low COC groups were 92 (22.33%), 130 (31.55%), and 190 (46.12%), respectively, and the mean follow-up periods for the 3 groups were 7.13, 7.12, and 7.27 years, respectively. After using inverse probability weighting, the intermediate and low COC groups were significantly associated with an increased risk of ESRD compared with the high COC group (adjusted hazard ratio (aHR) 1.36 [95% CI, 1.03-1.80] and aHR 1.76 [95% CI, 1.35-2.30], respectively). The intermediate and low COC groups were also significantly associated with the subsequent hospitalization compared with the high COC group (aHR 1.15 [95% CI, 0.99-1.33] and aHR 1.72 [95% CI, 1.50-1.97], respectively).
CONCLUSION: COC is related to ESRD onset and subsequent hospitalization among patients with DM. This study suggested that when DM patients keep visiting the same physician for managing their diseases, the progression of renal disease can be prevented.

Introduction

Diabetes mellitus (DM) is the most common chronic and progressive noncommunicable disease. The World Health Organization estimated the prevalence of DM in patients aged >18 years to have increased from 108 million (4.7%) to 422 million (8.5%) between 1980 and 2014.1 DM is a severe global public health problem and may result in morbidity, mortality, and economic burden on health care systems.2,3

DM complications can lead to end-stage renal disease (ESRD), lower extremity amputation, cardiovascular disease, and retinopathy, and consequently increase the rate of hospitalization.1,4–7 ESRD is a major complication in patients with DM. The United States Renal Data System Annual Data Report has indicated that 44% of patients with newly diagnosed ESRD have DM.8 In the USA, ESRD treatment ($34.3 billion) accounts for 6.3% of the total Medicare budget.9 Thus, ESRD strongly influences individual finances and social economies.

The management of DM can prevent or delay its progression to ESRD, which consequently reduces hospital admissions.2 Studies have proposed that effective DM management can slow the progression of kidney disease, prevent the onset of ESRD,10,11 and reduce the need for hospitalization.12,13 Therefore, it is crucial to understand the management of DM.

DM management is strongly related to continuity of care (COC).14 The COC level refers to the ongoing care partnership between patients and health care providers15 and is important to build physician–patient relationships.16 A long-term physician–patient sustained relationship facilitates enhancing mutual communication, effectively managing diseases, and improving disease outcomes,15,17 particularly in patients with chronic conditions.

Scholars have indicated that the increased continuity of ambulatory care can reduce cardiovascular mortality, cardiovascular events, and health care utilization and costs for patients with newly diagnosed hypertension, DM, and hypercholesterolemia.18 Hussey et al used the Medicare claims database for analyzing data on chronic diseases. They reported that high COC of patients is associated with lower rates of emergency department visits, complications, and episode costs.19 According to our review of relevant literature, no study has yet examined the association between COC and the subsequent development of ESRD. Therefore, this study investigated the association between COC and the risk of ESRD in patients with newly diagnosed type 2 DM.

High COC levels have been reported to be associated with low hospitalization risk13,20,21 and high medication adherence.22,23 Observational studies did not perform random assignment of different groups; therefore, confounding variables may have affected the health outcomes.24 Using traditional statistical methods to reduce bias results in sample size limitations because of the number of variables.24,25 To reduce the effects of confounders, we performed the inverse probability weighted analysis to induce similar covariate distributions between different COC groups. We examined the associations between COC, hospitalization, and ESRD in patients with newly diagnosed type 2 DM.

Methods

Dataset source

In the observational retrospective cohort study, data from 1996 to 2012 were retrieved from the Longitudinal Health Insurance Database 2005 (LHID 2005), which is provided by the National Health Research Institutes in Taiwan. The LHID 2005 consists of 1 million beneficiaries randomly sampled from the National Health Insurance Research Database (NHIRD) at the end of 2005. The NHIRD includes the enrollees’ demographic data, medical service providers and records, medication use, reimbursement amounts, treatment procedures, and disease diagnoses codes (according to the International Classification of Diseases, Ninth Revision, Clinical Modification; ICD-9-CM). The NHIRD has been widely used in many epidemiological studies in Taiwan.26–29 Gender, age, and average insured payroll-related amount distributions do not differ significantly between the patients in the LHID2005 and NHIRD.30 The NHIRD does not contain the personal information of patients, namely, body weight, height, family history, socioeconomic status, genetic factors, laboratory examination results, smoking, and alcohol consumption. Because LHID dataset is the encrypted secondary data, it does not contain any identified patient information. This study was approved by the TMU-Joint Institutional Review Board (approval number 201204036).

Design and study participants

This study enrolled patients who had received a primary diagnosis of type 2 DM (ICD-9-CM 250.xx, excluding type 1 DM [ICD-9-CM 250.x1 or 250.x3]) and had >3 annual physician visits during the study period under the condition that the visit intervals exceeded 28 days. These patients had been treated with antihyperglycemic agents, such as biguanides, sulfonamides, urea derivatives, α-glucosidase inhibitors, thiazolidinediones, and dipeptidyl peptidase-4 between January 1, 2000, and December 31, 2010. The aforementioned criteria were proposed by American Diabetes Association Clinical Practice Recommendations.31 The patients’ first clinical visit was recognized as the onset of DM (referred as index date hereinafter). The patients who had received dialysis treatment before the index date were excluded.

In this study, we used 1 year as the washout period because patients with newly diagnosed DM could be influenced by psychological factors that might produce misleading results regarding the effects of medical treatment.32 Therefore, we measured the patient’s COC and adherence from the second year of the index date of DM. To ensure that the outcomes were caused by COC and adherence, patients who had ESRD within 2 years of the index date were also excluded. The patients with a history of autoimmune disease or cancer within 1 year of the index date were excluded because these diseases were highly related to renal diseases. The patients with <3 physician visits within 2 years of the index day were excluded because it was inappropriate to use COC indicators for patients with fewer visits.18,22 Patients were followed up from the diagnosis time of type 2 DM to the development of ESRD, death, or the end of the study period (December 31, 2012). The patients were censored based on death or the end of the follow-up period.

COC during follow-up

Considering the variation and higher frequency of physician visits in Taiwan, this study used the COC index developed by Bice and Boxerman for measuring the COC, which is the common measure used in COC-related studies.15,17,19,33,34 The index was only calculated for the COC score on the basis of physician visits to Western-medicine doctors for DM. The COC index reflects the distribution of visits to different physicians and number of visits to each physician. The formula of the COC index is as follows: where N is the total number of physician visits, n is the number of visits to the jth physician, and M is the number of physicians. The COC index value ranges from 0 to 1, with a higher value corresponding to better COC. The COC score of 1 represents the patient visits to the same physician. Consistent with previous studies, we measured DM-related visits’ continuity score from the second to the third year after the first year of the index date and categorized the COC index into 3 equal tertiles (ie, low [0.00–0.43], intermediate [0.43–0.80], and high [0.80–1.00]) according to the distribution of scores across the entire study population because the COC index score lacks inherent clinical relevance.15,21,35,36

Outcome measurements and covariates

The primary outcome was ESRD, defined as patients continuously receiving dialysis treatment for 3 months. The secondary outcome was the first hospitalization caused by DM-related ambulatory care-sensitive condition admissions as defined by the Agency for Healthcare Research and Quality Prevention Quality Indicators, which also state that adequate management and outpatient care can prevent the need for hospitalization.12 Hospitalizations were defined as patients with a hospital stay of >1 day and the main ICD-9-CM diagnosis code for DM with short- or long-term complications (Table S1).

Covariates included age, gender, comorbidities, Charlson comorbidity index, number of antihyperglycemic drugs, number of physician visits, and medication adherence. The identified comorbidities (>3 physician visits; defined using ICD-9-CM codes) with potential influence outcomes in the present study included hypertension, dyslipidemia, gout, and chronic kidney disease (Table S1). The Charlson comorbidity index, which is a scoring system for weighting factors on the important concomitant disease, is defined by the ICD-9-CM.37 The number of antihyperglycemic drugs is a proxy indicator for disease severity; patients used more antihyperglycemic drugs corresponding to more severe disease. All covariates variables were defined in the first year of the index date.

Medication adherence was significantly associated with continuity of care.22,23 Therefore, we also considered the effect of the adherence to antihyperglycemic medication and measured medication adherence from the second to the third year after the first year of the index date. Medication adherence was defined as the medication possession ratio (MPR) to estimate the consumption (prescribed dosage) of antihyperglycemic medication. The MPR was dichotomized: patients with an MPR lower than the cutoff point of 80% were identified as nonadherent.10,38

Statistical analysis

The characteristic data of the study participants were first analyzed. The χ2-tests and one-way analysis of variance test were used to examine the associations between patient characteristics and COC tertiles. Because of the imbalance in the distribution of measured baseline covariates among the COC low, intermediate, and high groups, we applied inverse probability weighting analysis to induce similar covariate distribution between different COC groups’ baseline covariate distributions. Inverse probability weighting is based on the propensity score to receive unbiased estimates of average exposure effects.39,40 Prior evidence suggests that the propensity score model should include the confounders or the covariates affecting outcomes.40,41 Therefore, we included the covariates, such as age, gender, hypertension, dyslipidemia, gout, chronic kidney disease, Charlson comorbidity score, number of antihyperglycemic drugs, number of physician visits, and medication adherence in the propensity score model. In addition, we used the absolute standardized difference to assess the balance of baseline covariates among the 3 COC groups in the inverse probability weighting sample. The absolute standardized mean difference value of ≤10% indicates a negligible difference in covariates.42 We examined the associations between COC, hospitalization, and ESRD in patients with newly diagnosed type 2 DM. To investigate the risk of ESRD and hospitalization, we conducted multivariable analyses with Cox proportional hazards regression with the inverse probability weighting by using propensity scores. To further assess the robustness of our results, we also stratified the study population based on age and the adherence level. Predefined subgroups included age ≤65 years and >65 years and medication adherence and nonadherence. The data were managed using SAS version 9.4 (SAS Institute, Inc., Cary, NC, USA). The results were considered statistically significant when p<0.05.

Results

A total of 26,063 patients with newly diagnosed type 2 DM were included in the cohort (Figure 1). Table 1 presents the baseline characteristics and COC of the study patients. The mean age of the patients was 55.82±11.98 years, and male accounted for 54.02% of the population. The prevalence of comorbidities in the patients was 50.26% for hypertension, 24.26% for dyslipidemia, 8.61% for gout, and 1.97% for chronic kidney disease; 5.29% of the patients used >3 prescriptions of antihyperglycemic drugs, and 71.14% exhibited non-adherence with antihyperglycemic drugs. The mean number of physician visits was 16.96±9.69. The mean COC was 0.61±0.30 of the study patients. With increasing age, the COC tended to decrease. The patients with hypertension, gout, or chronic kidney disease showed lower COC score than did those without hypertension, gout, or chronic kidney disease. The patients with a higher Charlson comorbidity score showed low COC scores. The patients with adherence medication showed a higher COC score than did those with nonadherence medication.

Figure 1

Flowchart of patient selection.

Abbreviation: ESRD, end-stage renal disease.

Table 1

Level of continuity of care by patient characteristics before propensity weighting

Characteristic	Overall (n=26,063)		COC score		High COC group (0.80–1.00)	Intermediate COC group (0.43–0.80)	Low COC group (0–0.43)	p-value
	n	%	Mean	SD	n=8,666	n=8,843	n=8,554
COCI (mean, SD)					0.96 (0.06)	0.60 (0.11)	0.27 (0.10)	<0.001
Age, years (mean, SD)	55.82	11.98			55.12 (11.72)	55.96 (11.96)	56.38 (12.21)	<0.001
≤54	12,426	47.68	0.622	0.30	4,307 (49.70)	4,173 (47.19)	3,946 (46.13)	<0.001
55–64	7,072	27.13	0.617	0.29	2,385 (27.52)	2,417 (27.33)	2,270 (26.54)
65–74	4,844	18.59	0.597	0.30	1,509 (17.41)	1,668 (18.86)	1,667 (19.49)
≥75	1,721	6.60	0.561	0.30	465 (5.37)	585 (6.62)	671 (7.84)
Gender								0.013
Female	11,985	45.98	0.609	0.29	3,875 (44.71)	4,140 (46.82)	3,970 (46.41)
Male	14,078	54.02	0.614	0.30	4,791 (55.29)	4,703 (53.18)	4,584 (53.59)
Hypertension	13,098	50.26	0.601	0.30	4,169 (48.11)	4,426 (50.05)	4,503 (52.64)	<0.001
Dyslipidemia	6,323	24.26	0.612	0.29	2,085 (24.06)	2,175 (24.60)	2,063 (24.12)	0.661
Gout	2,243	8.61	0.593	0.30	703 (8.11)	741 (8.38)	799 (9.34)	0.010
Chronic kidney disease	514	1.97	0.557	0.30	135 (1.56)	174 (1.97)	205 (2.40)	<0.001
Charlson comorbidity score (mean, SD)	0.42	0.70			0.37 (0.64)	0.42 (0.70)	0.46 (0.75)	<0.001
Charlson comorbidity score								<0.001
0	17,652	67.73	0.620	0.30	6,069 (70.03)	5,969 (67.50)	5,614 (65.63)
1–2	7,976	30.6	0.598	0.30	2,508 (28.94)	2,723 (30.79)	2,745 (32.09)
≥3	435	1.67	0.512	0.29	89 (1.03)	151 (1.71)	195 (2.28)
Number of antihyperglycemic drugs (mean, SD)	1.48	0.64			1.47 (0.62)	1.48 (0.64)	1.49 (0.67)	0.084
Number of antihyperglycemic drugs								0.014
1	15,198	58.31	0.613	0.30	5,049 (58.26)	5,183 (58.61)	4,966 (58.05)
2	9,487	36.4	0.615	0.30	3,204 (36.97)	3,201 (36.20)	3,082 (36.03)
≥3	1,378	5.29	0.581	0.30	413 (4.77)	459 (5.19)	506 (5.92)
Number of physician visits, (mean, SD)	16.96	9.69			15.76 (8.55)	17.17 (9.56)	17.97 (10.71)	<0.001
Medication adherence								<0.001
Adherence (MPR ≥80)	7,523	28.86	0.664	0.29	3,088 (35.63)	2,505 (28.33)	1,930 (22.56)
Nonadherence (MPR <80)	18,540	71.14	0.591	0.30	5,578 (64.37)	6,338 (71.67)	6,624 (77.44)

Abbreviations: COC, continuity of care; COCI, COC index; MPR, medication possession ratio.

Among 26,063 patients in the overall cohort, 8,666 (33.25%), 8,843 (33.93%), and 8,554 (32.82%) were in the high, intermediate, and low COC groups, respectively. Before applying inverse probability weighting, the patients in the low COC group were older (56.38 [12.21]) than those in the high COC group (55.12 [11.72]) and intermediate COC group (55.96 [11.96]). The lower the COC group is, the higher the proportions of comorbidities (ie, hypertension, gout, and chronic kidney disease) and Charlson comorbidity scores are. Compared with the high and intermediate COC groups, the low COC group had a higher number of physician visits, a higher number of antihyperglycemic drugs, and a higher proportion of medication nonadherence.

After applying inverse probability weighting, the 3 COC groups were well balanced in all covariates (Figure 2). The mean age, gender proportions, hypertension, dyslipidemia, gout, and chronic kidney disease among the 3 COC groups were comparable. Moreover, the mean Charlson comorbidity score, proportions of the number of antihyperglycemic drugs, and medication adherence among the 3 COC groups became similar. On examining the standardized differences among the 3 COC groups (with the high COC group as the reference), the values for all covariates were found to be <10% (Table S2).

Figure 2

Level of continuity of care by patient characteristics before and after propensity weighting.

Abbreviations: COC, continuity of care; IPTW, inverse probability of treatment weighting; MPR, medication possession ratio.

The number of ESRD patients in the high, intermediate, and low COC groups were 92 (22.33%), 130 (31.55%), and 190 (46.12%), and the mean follow-up periods for the 3 groups were 7.13, 7.12, and 7.27 years, respectively. The number of hospitalization patients for the high, intermediate, and low COC groups were 351 (24.39%), 422 (29.33%), and 666 (46.28%), and the mean follow-up periods were 7.01, 6.95, and 6.99 years, respectively (Table S3).

The Kaplan–Meier analysis showed a significant difference in the ESRD survival and hospitalization survival among the COC groups (Figure 3). After using inverse probability weighting and adjusting for competing risk of death, the intermediate and low COC groups were significantly associated with an increased risk of ESRD compared with the high COC group (adjusted hazard ratio [aHR] for intermediate COC group was 1.36 [95% CI, 1.03–1.80] and low COC group was 1.76 [95% CI 1.35–2.30]). The intermediate and low COC groups were also significantly associated with the subsequent hospitalization before and after adjustment, compared with the high COC group (adjusted HRs for intermediate and low COC groups were 1.15 [95% CI, 0.99–1.33] and 1.72 [1.50–1.97]). The results of inverse probability weighting were very similar to those obtained without inverse probability weighting (Table 2).

Figure 3

Kaplan–Meier curves showing ESRD survival and hospitalization survival among COC groups.

Abbreviations: COC, continuity of care; ESRD, end-stage renal disease.

Table 2

Risk of ESRD and hospitalization by level of continuity of care, with adjustment for competing risk of death

Continuity of care	ESRD, HR (95% CI)				Hospitalization, HR (95% CI)
	Before propensity weighting		After propensity weighting		Before propensity weighting		After propensity weighting
	Crude	Adjusteda	Crude	Adjusteda	Crude	Adjusteda	Crude	Adjusteda
High COCI	1 (reference)	1 (reference)	1 (reference)	1 (reference)	1 (reference)	1 (reference)	1 (reference)	1 (reference)
Intermediate COCI	1.40 (1.05–1.85)	1.32 (1.00–1.76)	1.38 (1.04–1.82)	1.36 (1.03–1.80)	1.17 (1.01–1.36)	1.12 (0.97–1.31)	1.13 (0.98–1.32)	1.15 (0.99–1.33)
Low COCI	1.98 (1.52–2.58)	1.79 (1.36–2.34)	1.77 (1.36–2.31)	1.76 (1.35–2.30)	1.88 (1.64–2.15)	1.74 (1.52–2.00)	1.68 (1.47–1.93)	1.72 (1.50–1.97)

Note:

Adjusted for age, gender, hypertension, dyslipidemia, gout, chronic kidney disease, Charlson comorbidity score, number of antihyperglycemic drugs, number of physician visits, and medication adherence.

Abbreviations: COCI, continuity of care index; ESRD, end-stage renal disease; HR, hazard ratio.

Table 3 showed the interactive effects of hospitalization and COC level on the risk of ESRD. We treated the nonsubsequent hospitalization and high COC groups as the reference group. Compared to the reference group, the patients with subsequent hospitalization and low COC had the highest risk of ESRD, after using inverse probability weighting with an adjusted HR of 13.06 (95% CI, 9.29–18.35).

Table 3

Interactive effects of hospitalization and COC level on the risk of ESRD

Hospitalization	Continuity of care	Before propensity weighting, HR (95% CI)		After propensity weighting, HR (95% CI)
		Crude	Adjusteda	Crude	Adjusteda
Nonhospitalization	High COCI	1 (reference)	1 (reference)	1 (reference)	1 (reference)
	Intermediate COCI	1.29 (0.91–1.84)	1.25 (0.88–1.77)	1.23 (0.87–1.74)	1.22 (0.87–1.72)
	Low COCI	1.79 (1.29–2.49)	1.66 (1.19–2.32)	1.57 (1.13–2.18)	1.56 (1.12–2.16)
Hospitalization	High COCI	10.28 (6.70–15.78)	9.89 (6.41–15.25)	8.91 (5.82–13.65)	9.37 (6.11–14.36)
	Intermediate COCI	12.29 (8.30–18.19)	11.14 (7.47–16.60)	13.07 (8.97–19.06)	12.74 (8.74–18.59)
	Low COCI	14.68 (10.42–20.68)	12.77 (8.89–18.16)	13.07 (9.30–18.37)	13.06 (9.29–18.35)
p-value for trend		<0.001	<0.001	<0.001	<0.001

Note:

Adjusted for age, gender, hypertension, dyslipidemia, gout, chronic kidney disease, Charlson comorbidity score, number of antihyperglycemic drugs, number of physician visits, and medication adherence.

Abbreviations: COCI, continuity of care index; ESRD, end-stage renal disease; HR, hazard ratio.

To further assess the robustness of our results, we stratified the study population by age and medication adherence. Among patients aged <65 years with nonadherence, intermediate and low COC groups were significantly associated with the increased risk of ESRD, compared with the high COC group. No significant association was present between COC and ESRD occurrence in patients aged >65 years with adherence. In addition, the low COC group was significantly associated with subsequent hospitalization in all subgroups compared with the high COC group (Figure 4).

Figure 4

Multivariable stratified analyses for the association between continuity of care and ESRD and hospitalization after propensity weighting.

Note: aAdjusted for age, gender, hypertension, dyslipidemia, gout, chronic kidney disease, Charlson comorbidity score, number of antihyperlycemic drugs, number of visits, and adherrence level.

Abbreviations: COCI, continuity of care index; ESRD, end-stage renal disease; HR, hazard ratio; MPR, medication possession ratio.

Discussion

This large cohort study using inverse probability weighted analysis revealed that the low COC level was associated with the increased risk of ESRD and subsequent hospitalization in patients with newly diagnosed type 2 DM. With regard to the interactive effects of the hospitalization and COC level, the subsequent hospitalization and low COC level was associated with the highest risk of the ESRD. In addition, the results revealed the inverse association between COC values and comorbidities. Through subgroup analyses, low COC and intermediate COC levels were associated with the increased risk of ESRD in patients aged <65 years or nonadherence of medication. The low COC level was associated with the subsequent hospitalization in all subgroups.

Our results suggested that the low COC level leads to renal progression for patients with DM. A possible explanation is that COC indirectly affected the therapeutic relationship, mutual trust, understanding, communication quality, and satisfaction between patients and health care providers.15,16,43 A qualitative study interviewing 25 patients with DM indicated that relational continuity between patients and physicians is determined by whether patients trust the technical abilities of usual providers.44 Therefore, patients with DM without relational continuity might have poor quality of glycemic control,14 which consequently negatively affects renal function. As in previous studies, glycemic control was associated with renal progression45 and ESRD11 in patients with DM in this study.

Compared with the high COC level, the intermediate and low COC levels were associated with the increased subsequent hospitalization for patients with DM; these results are similar to those in previous studies.13,20,21 Hussey et al examined patients with congestive heart failure, chronic obstructive pulmonary disease, or type 2 DM, and reported that the patients with lower odds of hospitalization and complications have higher levels of care continuity.19 Recently, a study on the diagnosis of dementia Medicaid beneficiaries indicated that a low COC level is associated with high rates of hospitalization, emergency department visits, and testing and higher health care costs.16

We noticed that the low COC level is associated with ESRD and subsequent hospitalization in patients aged <65 years with DM or medication nonadherence. The possible explanation is that patients with early-onset DM have an aggressive disease phenotype than do those with late-onset DM and have a higher risk of the development and progression of complications.46,47 Another study on type 2 DM indicated that early-onset patients had a 20% higher risk of microalbuminuria compared with usual onset patients.48 In addition, Chen and Cheng reported that adherence to antidiabetic medications was positively associated with the degree of COC in patients with newly diagnosed DM.22 This positive association between COC and medication adherence may be because patients and physicians share information and the common objective of treating the disease.49 Therefore, a decreased COC level and adherence might strongly affect disease management, which consequently led to poor health care outcomes. A previous study has indicated that nonadherence to antihyperglycemic medications significantly increase the risk of ESRD.10

The advantages of the current study include the utilization of the LHID2005, which has a large sample size from a comprehensive nationwide database, which enabled our study patients to be highly representative of the entire population of Taiwan. To eliminate the possibility of selection bias, we included only patients newly diagnosed with DM between 2000 and 2010. In addition, we applied inverse probability weighting analysis to induce similar covariate distribution between different COC groups, reducing the effects of confounding in this observational study.

The study has several limitations. The claims data did not include biochemical data on glucose levels and renal function, such as fasting glucose, glycated hemoglobin, serum creatinine, and estimated glomerular filtration rate. Therefore, we could not assess the potential effect by renal function and control of blood glucose on the relationship between COC and the subsequent development of ESRD and hospitalization. However, we used the number of antihyperglycemic drugs as a proxy indicator for disease severity.10 In addition, we lacked details on several potential confounders, namely, level of physicians, body weight, height, family history, socioeconomic status, genetic factors, laboratory examination results, lifestyle, smoking, and alcohol consumption. The aforementioned factors may contribute to outcomes. To reduce these biases, we used inverse probability weighting and multivariable analysis, which were widely used in previous studies that have used claims data to adjust for potential confounders.16,28,29 Exclusion of patients who have less than 3 physician visits from the second to the third year after the first year of the index date might limit the representativeness of the study subjects. Finally, the finding could be different for treatment of type 2 DM worldwide because Taiwan’s health resource may be different from other countries. The generalizability of our findings to DM patients around the world might be limited. Thus, caution needs to be taken because of heath resource differences.

In conclusion, this study offers evidence demonstrating that COC is related to ESRD onset and subsequent hospitalization among patients with DM. We observed that COC has significant adverse effects, especially in patients aged <65 years or nonadherence of medication. Therefore, to reduce the risk of ESRD, we recommended that clinicians should continuously follow-up on the conditions of patients with DM and build trust and understanding between each other. In addition, when patients with type 2 DM should return to the same doctor and adhere to drug therapy, they can effectively reduce the subsequent hospitalization and prevent the onset of ESRD.

Table S1

ICD-9-CM codes for ACSC admissions and comorbidities

Code description	Codes
ACSC admissions
Type 2 DM with short-term complications (ketoacidosis, hyperosmolarity, or coma)	25010, 25012, 25020, 25022, 25030, 25032
Type 2 DM with long-term complications (renal, eye, neurological, circulatory, or complications not otherwise specified)	25040, 25042, 25050, 25052, 25060, 25062, 25070, 25072, 25080, 25082, 25090, 25092
Hypertension	401, 402, 403, 404, 405
Dyslipidemia	272, 273
Gout	274
Chronic kidney disease	016.0, 095.4, 189, 223, 236.9, 250.4, 271.4, 274.1, 283.11, 403–404, 440.1, 442.1, 447.3, 572.4, 581–584, 586–588, 591, 642.1, 646.2, 753, 794.4

Abbreviations: ACSC, ambulatory care sensitive condition; DM, diabetes mellitus; ICD-9-CM, International Classification of Diseases, Ninth Revision, Clinical Modification.

Table S2

Level of continuity of care by patient characteristics after propensity weighting

Characteristic	High COC group (0.80–1.00)	Intermediate COC group (0.43–0.80)	Low COC group (0–0.43)	Absolute standardized mean difference
				Intermediate vs high COCI group	Low vs high COCI group
Age, years (mean, SD)	55.84 (11.69)	55.83 (11.99)	55.76 (12.30)	−0.0008	−0.0067
Gender, %
Female	46.16	46.01	45.87	−0.3009	−0.5819
Male	53.84	53.99	54.13	0.3009	0.5819
Hypertension, %	50.19	50.29	50.02	0.2000	−0.3400
Dyslipidemia, %	24.16	24.25	24.28	0.2101	0.2801
Gout, %	8.49	8.61	8.59	0.4291	0.3578
Chronic kidney disease, %	1.89	1.97	1.99	0.5814	0.7249
Charlson comorbidity score, (mean, SD)	0.41 (0.67)	0.42 (0.70)	0.41 (0.71)	0.0088	0.0019
Number of antihyperglycemic drugs, (mean, SD)	1.48 (0.62)	1.48 (0.64)	1.48 (0.66)	0.0089	0.0018
Number of physician visits, (mean, SD)	16.88 (10.42)	16.99 (9.71)	16.54 (10.06)	0.0107	−0.0330
Medication adherence, %
Adherence (MPR ≥80)	29.69	28.90	30.17	−1.7360	−1.0481
Nonadherence (MPR <80)	70.31	71.10	69.83	1.7360	−1.0481

Abbreviations: COC, continuity of care; COCI, COC index; MPR, medication possession ratio.

Table S3

Number of ESRD and hospitalization patients in the high, intermediate, and low COC groups

COC	Follow-up (year)		ESRD		NonESRD		Follow-up (year)		Hospitalization		Non-hospitalization
	Mean	SD	n	%	n	%	Mean	SD	n	%	n	%
High COC group	7.13	3.03	92	22.33	8,574	33.43	7.01	3.02	351	24.39	8,315	33.77
Intermediate COC group	7.12	3.05	130	31.55	8,713	33.97	6.95	3.06	422	29.33	8,421	34.20
Low COC group	7.27	3.02	190	46.12	8,364	32.61	6.99	3.09	666	46.28	7,888	32.03

Abbreviations: COC, continuity of care; ESRD, end-stage renal disease.