Skip to main content

Evaluating the relationship between clinical and demographic characteristics of insulin-using people with diabetes and their health outcomes: a cluster analysis application

Abstract

Background

The aim of this study was to determine how clusters or subgroups of insulin-treated people with diabetes, based upon healthcare resource utilization, select social demographic and clinical characteristics, and diabetes management parameters, are related to health outcomes including acute care visits and hospital admissions.

Methods

This was a non-experimental, retrospective cluster analysis. We utilized Aetna administrative claims data to identify insulin-using people with diabetes with service dates from 01 January 2015 to 30 June 2018. The study included adults over the age of 18 years who had a diagnosis of type 1 (T1DM) or type 2 diabetes mellitus (T2DM) on insulin therapy and had Aetna medical and pharmacy coverage for at least 18 months (6 months prior and 12 months after their index date, defined as either their first insulin prescription fill date or their earliest date allowing for 6 months’ prior coverage). We used K-means clustering methods to identify relevant subgroups of people with diabetes based on 13 primary outcome variables.

Results

A total of 100,650 insulin-using people with diabetes were identified in the Aetna administrative claims database and met study criteria, including 11,826 (11.7%) with T1DM and 88,824 (88.3%) with T2DM. Of these 79,053 (78.5%) people were existing insulin users. Seven distinct clusters were identified with different characteristics and potential risks of diabetes complications. Overall, clusters were significantly associated with differences in healthcare utilization (emergency room visits, inpatient admissions, and total inpatient days) after multivariable adjustment.

Conclusions

This analysis of healthcare claims data using clustering methodologies identified meaningful subgroups of patients with diabetes using insulin. The subgroups differed in comorbidity burden, healthcare utilization, and demographic factors which could be used to identify higher risk patients and/or guide the management and treatment of diabetes.

Peer Review reports

Introduction/background

Diabetes is a complex, chronic illness that affects about 30 million people or 9.4% of the population in the United States [1, 2]. The population with diabetes continues to grow and the percentage of adults with diabetes increases with age, where approximately 25% of all adults over 65 years of age have diabetes [3]. Diabetes remains the seventh leading cause of death in the United States with approximately 80,000 death certificates in 2015 listing diabetes as the underlying cause of death [3], although the actual number of diabetes-associated deaths may be higher as diabetes is often underreported as a cause of death [3].

The majority of people with diabetes are classified as having type 2 diabetes mellitus (T2DM), occurring in about 90 to 95% of diagnosed cases [4, 5], with approximately 5% of diagnosed cases classified as type 1 diabetes mellitus (T1DM) [4]. T1DM is characterized by the lack of insulin production, and as such patients require insulin to survive [5]. T2DM is characterized by the body’s resistance to insulin action in addition to a relative insulinopenia [6], and is usually diagnosed in adults. Genetic factors and lifestyle play an important role in disease progression for people with either T1DM or T2DM [4].

People with diabetes visit physician offices and emergency rooms (ER) more frequently than people without diabetes, are more likely to be admitted to the hospital, and more commonly receive home health care [7]. The financial burden of diabetes to individuals and society is estimated to be a total of $327 billion, including $237 billion in direct costs and $90 billion in reduced productivity in 2018 [3, 7]. Healthcare expenditures are about 2.3 times higher for people with diabetes than people without diabetes with an average medical diabetes-related expenditures of approximately $9600 per year in the United States [7].

Although all people with T1DM use insulin, insulin use varies among people with T2DM. Most are initially managed with oral hypoglycemic agents but over time, most will ultimately utilize insulin. The majority of the diabetic population uses injection-based administration with insulin pens or vials and syringes while the remainder uses insulin pumps (continuous subcutaneous insulin infusion) [8, 9]. Since 2005, there has been a significant increase in the use of insulin pens while the use of vials and syringes to deliver insulin has decreased over time [10]. Insulin pens have simplified the administration of insulin, resulting in more accurate and easier delivery of insulin relative to vials and syringes [11]. Pumps deliver rapid-acting insulin throughout the day. Diabetes-related technologies are available to monitor blood glucose and include the standard blood glucose monitor (single reading) and continuous blood glucose monitoring systems (real time and intermittently scanned) [8].

In 2005–2012, among patients who had any insulin use, only 31.4% had an hemoglobin A1c (HbA1c) < 7%; therefore, as glycemic control is still not attained by most people with diabetes, there is a need for new approaches to identify subgroups of people with diabetes (T1DM or T2DM) who might have risk factors that could impact treatment decisions and targeted disease management [12]. Given the differences across the US population especially in disease severity and utilization of healthcare, we sought insight from a large, real-world database utilizing Aetna’s administrative claims. In this study, we used clustering techniques to identify subgroups of people treated with insulin based upon healthcare resource utilization, select social demographic and clinical characteristics, and diabetes management parameters. We then related these subgroups to health outcomes, including ER visits, hospital admissions, and total inpatient days.

Methods

Study design and data sources

We used a non-experimental, retrospective design in this study utilizing Aetna’s administrative claims data containing membership, eligibility, medical claims, pharmacy claims, laboratory test results, and data derived for Aetna’s care management processes (Health Profile Database [13]) for Aetna fully insured Commercial and Medicare Advantage members, with services dates from 01 January 2015 to 30 June 2018. All of the data used in this study were fully de-identified and Health Insurance Portability and Accountability Act compliant.

Sample selection and patient population

Using the International Classification of Diseases Ninth and Tenth Revision (ICD9-CM v1 and ICD10-CM) diagnostic codes, the study included adults 18 years of age and older who had a diagnosis of T1DM or T2DM who utilized insulin therapy and had Aetna medical and pharmacy coverage for at least 18 continuous months (6 months of coverage prior and 12 months after their index date, defined as either their first insulin prescription fill date or their earliest insulin fill date allowing for 6 months prior coverage; Fig. 1). People with diabetes were excluded from the analyses if any of the following criteria were met during the entire study period: had ≥1 inpatient or outpatient medical claim with a diagnosis in any position of gestational diabetes, steroid-induced diabetes, or metastatic cancer; had indications of hospice use; or were enrolled in Aetna’s Compassionate Care Program. The study was approved by the Sterling Institutional Review Board (Atlanta, Georgia, USA).

Fig. 1
figure 1

Member selection - population funnel

Statistical analysis

Descriptive analyses were performed initially on the overall study population and for T1DM and T2DM separately. After clustering, descriptive analyses were completed for each cluster. Descriptive statistics were generated for all demographic characteristics and utilization measures as applicable to the type of variable. Continuous variables were described using means with standard deviations (SD) or medians with first and third quartiles if data were highly skewed. Categorical and binary variables were described using counts and percentages. Healthcare utilization were reported as means and SD and number of patients with ≥1 visit.

Cluster analysis

A cluster analysis was performed on 13 pre-period variables that were hypothesized to be both related to outcomes and to the use of diabetes technology (blood glucose monitor [BGM], continuous glucose monitoring [CGM], or insulin pumps). Variables were identified by the study team based on clinical opinion and were limited to variables found in the administrative claims database. They included age, number of endocrinologist visits, diabetes complications severity index (DCSI) [14], Charlson comorbidity index (CCI) [15], number of HbA1c tests, number of months on insulin, number of ER visits, total number of inpatient days, number of discordant comorbidities, number of concordant comorbidities, number of medical claims, proportion diabetes-related claims, and estimated income. Concordant comorbidities were defined as those that share the same pathophysiologic risk profile and management plan as with diabetes, whereas discordant conditions are those that are not directly related to diabetes in either their pathophysiology or management [16]. As the underlying pathophysiology differs for T1DM and T2DM, concordant and discordant were defined separately for each type of diabetes (see Supplemental Materials Table S1). Clinically dominant comorbidities, those that are so complex or serious that they tend to eclipse the management of other conditions [16] such as end stage renal failure or dementia, were considered but not found in the study population. Claims were considered diabetes related if they contained a diagnosis for diabetes in any position. All 13 cluster variables (mentioned above) had significant bivariate associations with utilization of diabetes technology and were standardized according to the type of variable. We opted to use variables associated with the utilization of diabetes technology in the cluster analysis to identify unexpected patterns in the characteristics of patients who utilized different types of technology. After clustering, utilization of diabetes technology along with other relevant variables such as diabetes type, insurer, and HbA1c levels (when available, for 31.2% of the population) were reported to help characterize the clusters.

We utilized K-means methodology to identify the clusters based on the following 13 variables: age, number of endocrinologist visits, DCSI, CCI, number of HbA1c tests, number of months on insulin, number of ER visits, total number of inpatient days, number of discordant comorbidities, number of concordant comorbidities, number of medical claims, proportion diabetes-related claims, and estimated income. To determine the optimal number of clusters, we used the Jump method [17], a non-parametric method for choosing the number of clusters based on “distortion” which is a measure of within-cluster dispersion. We chose the Jump method as it is a simple and readily available method that has performed well against other competing methods in simulated data analyses [17].

To identify which variables played a key role in the formation of clusters and to obtain simple descriptive rules for clusters, classification and regression trees (CART) were used with cluster assignment as the “outcome” variable and all variables used in clustering as the “predictors.” This analysis was performed using the Ctree function in the R package partykit.

Outcomes of interest, measured in the 12-month follow-up period, included: all-cause ER utilization, all-cause inpatient hospitalization, and total inpatient days. To determine if the clusters were associated with the outcomes measured during the study, multivariable generalized linear models were applied with clusters as the independent variables and other covariates or interactions to account for confounding. The generalized linear models included binomial distribution and logit link function (ER and inpatient hospitalization outcomes) and negative binomial distribution with log link function (total inpatient days outcome). For the multivariable regression analyses, we used a stepwise selection methodology (or backward elimination) with variable removal when p ≥ 0.05. Akaike Information Criterion was also used to determine the best fit model. All statistical analyses were conducted using SAS Enterprise Guide Version 6.1 (SAS Institute Inc., Cary, NC, USA) or R version 3.2.5 (R Foundation for Statistical Computing, Vienna, Austria).

Results

Study population

A total of 100,650 insulin-using people with diabetes met the study criteria, including 11,826 (11.7%) with T1DM and 88,824 (88.3%) with T2DM. Table 1 shows the pre-period demographic characteristics for the total population and differences between people with T1DM and T2DM. The majority of study patients (78.5%) were existing insulin users at the index date, whereas the remaining 21.5% were new insulin initiators. The mean age (SD) was 62.2 (14.2) years with a mean age (SD) of 45.3 (16.3) years for those with T1DM and 64.4 (12.3) years for those with T2DM. There were slightly more males than females (51.6% male versus 48.3% female). The median (interquartile range) household income was $54,143 ($43,184 to $69,606) for the total population with a slightly higher income for those with T1DM of $63,304 ($49,195 to $81,547) compared to T2DM of $53,299 ($42,675 to $68,059).

Table 1 Pre-period1 demographic and clinical characteristics (baseline)

Overall, the mean (SD) CCI score was 1.9 (1.6) and the mean DCSI was 1.1 (1.5), both of which were lower for those with T1DM (Table 1). A total of 60,668 people (60.3%) had one or more HbA1c test performed and the average HbA1c value was 8.8% among those individuals (31.2%) with values in the database. Approximately 6.4% of the population used an insulin pump (37.9% T1DM and 2.2% T2DM) and 67.5% used a pen (38.7% T1DM and 71.4% T2DM). Blood glucose meters were utilized by approximately one-third of the population (33.3%). On the other hand, CGMs were only used by 3.0% overall, where T1DM patients had greater use compared to those with T2DM (19.2 and < 1%, respectively).

The healthcare utilization over the 6-month pre-period is shown in Table 2. The majority of the population had primary care physician (PCP) visits (72.3%), but the proportion with endocrinologist and cardiologist visits was lower at 21.3 and 17.8%, respectively. Inpatient hospitalizations occurred for 14.0% of the population and 18.8% had ER visits. Overall, people with T1DM had lower utilization of all-cause and diabetes-related healthcare services than people with T2DM, with the exception of endocrinologist visits.

Table 2 Pre-period1 healthcare utilization (baseline)

Cluster formation and analysis

We identified seven clusters of people with diabetes, which had distinguishable characteristics and risk factors. Characteristics of the seven identified clusters are shown in Table 3. Seven of the 13 clustering variables studied using the CART analysis were identified as being the most important factors for cluster formation. These included number of endocrinology visits, total inpatient days, concordant comorbidities, number of ER visits, comorbidity burden as measured by CCI and DCSI scores, and percentage of diabetes-related medical claims. Figure 2 illustrates the CART analysis after the clusters were formed.

Table 3 Description of people with diabetes identified through cluster analysis
Fig. 2
figure 2

Cluster analysis schematic. Note: Rules within the tree have been transformed back to original un-standardized definitionsAbbreviations: CCI = Charlson Comorbidity Index; DCSI = Diabetes Complications Severity Index; DM = diabetes mellitus; Endo = endocrinologist; ER = emergency room

The seven clusters fell into three main groupings in hierarchical order (Fig. 3): those with endocrinology visits (Clusters 2 and 6) with differing concordant burden profiles; those with high acute care utilization (Clusters 4 and 5) with differences in inpatient days and ER visits; and those with low ER utilization (Clusters 1, 3, and 7) with varying comorbidity burden, diabetes utilization, and complication profiles.

Fig. 3
figure 3

Cluster descriptions

Endocrinology utilizers clusters

Cluster 2 all had at least one endocrinologist visit (mean = 1.8 visits/patient) and the average number of concordant comorbidities was 4.6. Cluster 6 had the youngest mean age (38.3 years) and the second highest utilization of endocrinologists (mean = 0.6 visits/patient) but had few concordant comorbidities (mean = 0.5). Additionally, the DCSI and CCI scores were relatively low.

High acute care utilizers clusters

Those in Cluster 4 had the highest number of inpatient days with nearly all (99.1%) visiting the hospital at least once with the vast majority having an inpatient stay that was diabetes related (95.7%). Cluster 4 also had the highest overall number of medical claims and the highest comorbidity and complications burden (per the CCI and DCSI). Cluster 5 all experienced an ER visit, usually related to diabetes (80.8%).

Low ER utilizers clusters

Those people in Clusters 1, 3, and 7 had few ER and endocrinology visits. Cluster 1 had a relatively high complication and comorbidity burden, second only to DCSI and CCI scores in Cluster 4, but those people in Cluster 1 had no ER utilization. In contrast, Clusters 3 and 7 had relatively low CCI and DCSI scores. Cluster 7 had the lowest overall number of medical claims (median = 3.00) and the lowest CCI scores (mean = 0.30) suggesting that Cluster 7 was the healthiest population, however only 15.2% of those in Cluster 7 had at least one HbA1c test performed. While those in Cluster 7 had nearly no diabetes-related claims (mean = 0.19), more than three-quarters of the claims for Cluster 3 were diabetes related (mean = 0.78).

Other notable differences were observed across clusters for variables not included in the cluster formations. Clusters 6 and 2 had the highest proportion of T1DM patients (69.6 and 13.6%, respectively) and diabetes technology use (pump, CGM, or BGM). Consistent with the young mean age, nearly all patients in Cluster 6 had commercial insurance. The highest proportion of Medicare Advantage members (and oldest mean age) was in Clusters 1 and 4 (> 80%). The highest mean HbA1c value (among the subset of the population with available values) was in Cluster 3 (9.4%). On the other hand, the lowest mean HbA1c of 8.3% was observed in Cluster 4, which also had the highest proportion without prior insulin use in the pre-period and the lowest mean number of months where insulin was filled.

Multivariable modeling analysis

For the multivariable modeling analysis, three key outcomes were chosen including all-cause ER visit (Fig. 4), all-cause inpatient hospitalization (Fig. 5), and total inpatient days (Fig. 6) during the follow-up period. Across all models, clusters were independently associated with the outcomes of interest after controlling for covariates that were potentially related to cluster assignment and/or outcomes of interest. Clusters 1, 2, 4, and 5 had higher odds of an ER visit, whereas Cluster 6 had lower odds compared to Cluster 7. Clusters 2, 4, and 5 also had higher odds of having a hospital admission. If hospitalized, Cluster 4 had significantly increased total inpatient days compared to Cluster 7.

Fig. 4
figure 4

All-cause emergency room visits - model results. Abbreviations: BGM = blood glucose monitor; CCI = Charlson Comorbidity Index; CGM = continuous glucose monitor; DCSI = diabetes complications severity index; ER = emergency room; HbA1c = hemoglobin A1c; IP = inpatient

Fig. 5
figure 5

All-cause inpatient hospitalization - model results. Abbreviations: BGM = blood glucose monitor; CCI = Charlson Comorbidity Index; CGM = continuous glucose monitor; DCSI = diabetes complications severity index; ER = emergency room; HbA1c = hemoglobin A1c; IP = inpatient

Fig. 6
figure 6

Total inpatient days - model results. Abbreviations: BGM = blood glucose monitor; CCI = Charlson Comorbidity Index; CGM = continuous glucose monitor; DCSI = diabetes complications severity index; ER = emergency room; HbA1c = hemoglobin A1c; IP = inpatient

Discussion and conclusions

The aim of this study was to determine how clusters or subgroups of insulin-treated people with diabetes, based upon healthcare resource utilization, select social demographic/clinical characteristics, and diabetes management parameters, are related to health outcomes including acute care (ER and hospital inpatient) visits and total inpatient days. We did this study to help identify groups of patients that may be amenable to emerging diabetes management technologies. In this study, we identified seven clusters of insulin-treated people with diabetes, which have different patterns of healthcare utilization and diagnosed comorbidities in a large healthcare claims database. The most important factors in defining the clusters were the number of endocrinology visits, total inpatient days, concordant comorbidities, number of ER visits, comorbidity burden as measured by CCI and DCSI scores, and percentage of diabetes-related medical claims. Multivariable modeling showed that these clusters are significantly associated with ER visits, inpatient hospitalizations, and total inpatient days, suggesting that this approach may help identify patients at greater need for targeted disease management efforts at the population level. The clusters also offer providers clinically relevant information regarding treatment decisions for a patient population with diabetes.

Cluster analyses can reveal how variables, in our case administrative claims from people with diabetes, are related in complex datasets. The use of cluster analyses in healthcare decision making is still relatively uncommon but appears to be gaining acceptance [18,19,20,21,22,23]. Our work builds upon a few previously published cluster analyses in diabetes, which focused on readiness of CGM and other diabetes-related devices, self-management patterns in a pediatric population, factors influencing people with diabetes who have poorly controlled conditions, and identifying diabetes phenotypes [23,24,25,26,27]. These previous studies involved smaller numbers of participants from relatively homogeneous populations (e.g., T1DM registry) and/or more controlled conditions (e.g., clinical trial). In contrast, our study used a large healthcare claims database to evaluate whether routinely available data could identify relevant subgroups of insulin users.

Not only did clusters differ with respect to the specific variables used to form them (by design), but also on important other characteristics. Clusters 2 and 6 were formed primarily based on the use of endocrinologists. Not surprisingly, these were the clusters with the highest proportion of people with T1DM as well as utilizers of diabetes technology (pump, BGM, or CGM). Those in Cluster 2, however, had higher comorbidity burden and mean number of HbA1c tests than those in Cluster 6, but there was little difference in mean HbA1c values for these two clusters.

Two clusters were identified with high levels of acute care utilization. Those in Cluster 4 had the highest total inpatient days and everyone in Cluster 5 had an ER visit. These clusters differed, however, in their comorbidity burdens and glycemic control. Interestingly, the lowest mean observed HbA1c value was for Cluster 4, with the highest levels of overall medical utilization (median of 43.0 claims), acute care utilization (99.1% had an inpatient hospitalization), and highest CCI and DCSI scores. These results could suggest that a high burden of comorbidities or diabetes complications and increased interactions with hospitals facilitated more intensive diabetes management. However, because HbA1c values were only available on a subset of the study population (approximately 30%), additional analyses on datasets with more complete HbA1c data are needed to confirm this finding.

Conversely, higher mean HbA1c’s were observed among Clusters 3, 5, and 7 (in order of highest to lowest values). Clusters 3 and 7 differ from Cluster 5 in that they fell into the low utilization grouping (both acute care and overall utilization via number of medical claims) and had among the lowest CCI and DCSI scores. They differed from each other in one key aspect: the proportion of medical claims that were diabetes-related. Approximately three-fourths of the claims for Cluster 3 were related to diabetes, compared to less than 20% in Cluster 7. Because the CCI and DCSI scores are derived from the presence of diagnosis codes in claims data, on one hand it is not surprising for these clusters who have the lowest overall number of medical claims to have the lowest scores due to fewer opportunities to derive those diagnoses. But, on the other hand, the lack of diagnoses of comorbidities in the observed claims or lack of encounters altogether could suggest a healthier underlying population. Either way, the relatively high observed HbA1c values along with the low rates of interactions with healthcare providers suggested suboptimal diabetes self-management.

The current study demonstrated that even after adjusting for other covariates, cluster assignment was significantly predictive of future outcomes. Specifically, cluster assignment was associated with the likelihood of experiencing an ER or hospital inpatient visit and the total number of inpatient days for those with an admission. These results suggest that the specific combination of variables used in the cluster formations shed additional light onto the risk of untoward outcomes above and beyond traditional risk stratification, for example, based upon parameters including diabetes type, age, and HbA1c.

Furthermore, as these clusters were derived from variables routinely found in healthcare claims data where detailed clinical data are often missing, this approach could aid healthcare payers with population management efforts. We found some clusters utilizing less healthcare resources had higher observed mean HbA1c levels. This finding could suggest population management efforts in diabetes that are targeted at some of the lower healthcare utilizers in efforts to improve glycemic control, which could yield better long-term health outcomes for patients and improved quality metric ratings for providers and payers.

This study has limitations that should be considered. The cluster analysis was based on administrative claims data that rely on diagnostic codes entered by the healthcare provider for billing purposes. As such, they are proxies for clinical outcomes and may be prone to data coding errors or inaccuracies in patient records. Within this particular population, patients may have diagnostic codes for both type 1 and type 2 diabetes during the course of the study period, and while our algorithm to classify patients considered the preponderance of diagnostic codes among all available, in conjunction with treatment patterns common to each type of diabetes, it did not preclude the possibility of misclassification. Additionally, as of October 2015, all claims switched from ICD-9 to ICD-10 (study period January 2015 to June 2018) [28]. This switch possibly could have led to inaccuracies in coding due to unfamiliarity with the new system or mistakes in cross walking codes from ICD-9 to ICD-10 by providers. This potential for error should have had limited impact since condition identification was based on families of codes with multiple codes within a family and not on any single code. The comparisons of HbA1c values were incomplete as only a subset of patients (~ 30%) had values in the database. A number of relevant risk factors, including insulin dosing, diet, and exercise, were not available in the database.

Despite the limitations, this study was based on a well-developed study design and included data from a large number of insulin-using people with either T1DM or T2DM. Evaluating the impact of patient-reported outcomes and more socioeconomic data on cluster formations would be of interest to study in the future.

Conclusion

In conclusion, we demonstrated that clustering analyses of healthcare claims data identified meaningful subgroups of patients that differed in comorbidity burden, healthcare utilization, and demographic factors. These clusters were found to be significantly associated with future outcomes indicating that providers and population health managers may be able to better estimate risk, based upon combinations of specific variables, and modify and/or personalize treatment accordingly.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due to their proprietary nature, having been licensed only for the current study, and their protection under the Health Insurance Portability and Accountability Act (HIPAA). Subject to the data owner’s determination of reasonableness and only after execution of a confidentiality agreement with the data owner, among other requirements, the data is available from the corresponding author upon reasonable request.

Abbreviations

BGM:

blood glucose monitor

CART:

classification and regression trees

CCI:

Charlson comorbidity index

CGM:

continuous glucose monitoring

DCSI:

diabetes complications severity index

ER:

emergency room

HbA1c:

hemoglobin A1c

ICD:

International Classification of Diseases

IP:

inpatient

PCP:

primary care physician

SD:

standard deviation

T1DM:

type 1 diabetes mellitus

T2DM:

type 2 diabetes mellitus

US, USA:

United States/United States of America

References

  1. American Diabetes Association. Statistics about Diabetes. https://www.diabetes.org/resources/statistics/statistics-about-diabetes. Accessed 19 August 2019.

  2. American Diabetes Association. 1. Introduction: Standard of Medical Care in Diabetes – 2019. Diab Care 2019;42(Suppl. 1):S1–S2. https://0-doi-org.brum.beds.ac.uk/10.2337/dc19-SINT01. https://care.diabetesjournals.org/content/42/Supplement_1. Accessed 19 August 2019.

  3. Centers for Disease Control and Prevention. National Diabetes Statistics Report, 2017. Atlanta: Centers for Disease Control and Prevention, U.S. Dept of Health and Human Services; 2017.

    Google Scholar 

  4. American Diabetes Association. 2. Classification and Diagnosis of Diabetes: Standard of Medical Care in Diabetes – 2019. Diab Care. 2019;42(Suppl. 1):S13–28. https://0-doi-org.brum.beds.ac.uk/10.2337/dc19-S002.

    Article  Google Scholar 

  5. American Diabetes Association. 9. Pharmacologic approaches to glycemic treatment: Standard of Medical Care in Diabetes – 2019. Diab Care. 2019;42(Suppl. 1):S90–S102. https://0-doi-org.brum.beds.ac.uk/10.2337/dc19-S009.

    Article  Google Scholar 

  6. Schwartz SS, Epstein S, Corkey BE, Grant SFA, Gavin JR III, Aguilar RB. The time is right for a new classification system for diabetes: rationale and implications of the β-cell–centric classification schema. Diabetes Care. 2016;39(2):179–86. https://0-doi-org.brum.beds.ac.uk/10.2337/dc15-1585.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. American Diabetes Association. Economic costs of diabetes in the U.S. in 2017. Diabetes Care. 2018;41:917–28. https://0-doi-org.brum.beds.ac.uk/10.2337/dci18-0007.

    Article  PubMed Central  Google Scholar 

  8. American Diabetes Association. 7. Diabetes Technology:Standard of Medical Care in Diabetes – 2019. Diab Care. 2019;42(Suppl. 1):S71–80. https://0-doi-org.brum.beds.ac.uk/10.2337/dc19-S007.

    Article  Google Scholar 

  9. Fry A. Insulin delivery device technology 2012: where are we after 90 years? J Diabetes Sci Technol. 2012;6(4):947–53. https://0-doi-org.brum.beds.ac.uk/10.1177/193229681200600428.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Perez-Nieves M, Jiang D, Eby E. Incidence, prevalence, and trend analysis of the use of insulin delivery systems in the United States (2005 to 2011). Curr Med Res Opin. 2015;31(5):891–9. https://0-doi-org.brum.beds.ac.uk/10.1185/03007995.2015.1020366.

    Article  CAS  PubMed  Google Scholar 

  11. Spollett GR. Improved disposable insulin pen devices provide an alternative to vials and syringes for insulin administration. Diabetes Spectrum. 2012;25(2):117–22. https://0-doi-org.brum.beds.ac.uk/10.2337/diaspect.25.2.117.

    Article  Google Scholar 

  12. Selvin E, Parrinello CM, Daya N, Bergenstal RM. Trends in insulin use and diabetes control in the U.S.: 1988-1994 and 1999-2012. Diabetes Care. 2016;39(3):e33–5. https://0-doi-org.brum.beds.ac.uk/10.2337/dc15-2229.

    Article  PubMed  Google Scholar 

  13. Hanchak NA, Murray JF, Hirsch A, McDermott PD, Schlackman N. USQA health profile database as a tool for health plan quality improvement. Manag Care Q. 1996;4:58–69. Available at: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/10157263. .

  14. Young BA, Lin E, Von Korff M, et al. Diabetes complications severity index and risk of mortality, hospitalization, and healthcare utilization. Am J Manag Care. 2008;14(1):15–23.

    PubMed  PubMed Central  Google Scholar 

  15. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83. https://0-doi-org.brum.beds.ac.uk/10.1016/0021-9681(87)90171-8.

    Article  CAS  PubMed  Google Scholar 

  16. Piette JD, Kerr EA The impact of comorbid chronic conditions on diabetes care. Diabetes Care 2006;29(3):725–731. 10.2337/diacare.29.03.06.dc05–2078.

  17. Sugar CS, James GM. Finding the number of clusters in a dataset. J Am Stat Assoc. 2003;98(463):750–63. https://0-doi-org.brum.beds.ac.uk/10.1198/016214503000000666.

    Article  Google Scholar 

  18. Cornell JE, Pugh JA, Willliams JW, et al. Multimorbidity clusters: clustering binary data from multimorbidity clusters: clustering binary data from a large administrative medical database. Appl Multivar Res. 2017;12(3):163–82. https://0-doi-org.brum.beds.ac.uk/10.22329/amr.v12i3.658.

    Article  Google Scholar 

  19. Dilts D, Khamalah J, Plotkin A. Using cluster analysis for medical resource decision making. Med Decis Making. 1995;15(4):333–47. https://0-doi-org.brum.beds.ac.uk/10.1177/0272989X9501500404.

    Article  CAS  PubMed  Google Scholar 

  20. Newcomer SR, Steiner JF, Bayliss EA. Identifying subgroups of complex patients with cluster analysis. Am J Manag Care. 2011;17(8):e324–32.

    PubMed  Google Scholar 

  21. Liao M, Li Y, Kianifard F, Obi E, Arcona S. Cluster analysis and its application to healthcare claims data: a study of end-stage renal disease patients who initiated hemodialysis. BMC Nephrol. 2016;17(1):25. https://0-doi-org.brum.beds.ac.uk/10.1186/s12882-016-0238-2.

    Article  PubMed  PubMed Central  Google Scholar 

  22. van den Berge MJC, Free RH, Arnold R, de Kleine E, Hofman R, van Dijk JMC, et al. Cluster analysis to identify possible subgroups in tinnitus patients. Front Neurol. 2017;8:115. https://0-doi-org.brum.beds.ac.uk/10.3389/fneur.2017.00115.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ahlquist E, Storm P, Käräjämäki A, et al. Novel subgroups of adult-onset diabetes and their association with outcomes: a data-driven cluster analysis of six variables. Lancet Diabetes Endocrinol. 2018;6(5):361–9. https://0-doi-org.brum.beds.ac.uk/10.1016/S2213-8587(18)30051-2.

    Article  Google Scholar 

  24. Shahreza FA, Hazar N. Using routine data to categorize poor control diabetic patients: an application of cluster analysis technique. Iran J Public Health. 2017;46(1):120–7.

    Google Scholar 

  25. Tanenbaum ML, Adams RN, Lanning MS, Hanes SJ, Agustin BI, Naranjo D, et al. Using cluster analysis to understand clinician readiness to promote continuous glucose monitoring adoption. J Diabetes Sci Technol. 2018;12(6):1108–15. https://0-doi-org.brum.beds.ac.uk/10.1177/1932296818786486.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Tanenbaum ML, Adams RN, Lanning MS, et al. From wary wearers to d-embracers: personas of readiness to use diabetes devices. J Diabetes Sci Technol. 2018;12(6):1101–7. https://0-doi-org.brum.beds.ac.uk/10.1177/1932296818793756.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Rohan JM, Rohan JM, Delamater A, et al. Identification of self-management patterns in pediatric type 1 diabetes using cluster analysis. Pediatr Diabetes. 2011;12(7):611–8. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1399-5448.2010.00752.x.

    Article  PubMed  PubMed Central  Google Scholar 

  28. US Department of Health and Human Services. HIPAA administrative simplification: modifications to medical data code set standards to adopt ICD-10-CM and ICD10-PCS. [FR Doc. E9–743.]. Fed Reg. 2014;79(149):45128 https://www.govinfo.gov/content/pkg/FR-2014-08-04/pdf/2014-18347.pdf.

    Google Scholar 

Download references

Acknowledgments

The authors wish to acknowledge Debra Hidayetoglu for providing medical writing services on behalf of Covance, Inc.

Funding

This work was supported by Eli Lilly and Company.

Author information

Authors and Affiliations

Authors

Contributions

ELE, EM, BDB, and KS conceptualized the study; AE and IL conducted the analyses; ELE and AE wrote the manuscript; and EM, BDB, KS, and IL critically commented on the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Elizabeth L. Eby.

Ethics declarations

Ethics approval and consent to participate

This study complies with relevant guidelines and regulations on observational studies. The study was approved by the Sterling Institutional Review Board (Atlanta, Georgia, USA). A waiver of informed consent was granted by the Sterling Institutional Review Board. All of the data used in this study were fully de-identified and Health Insurance Portability and Accountability Act compliant.

Consent for publication

Not applicable.

Competing interests

ELE, EM, IL, and BDB are employees and shareholders of, and serve as authors for Eli Lilly and Company and AE and KS are employees of, and serve as authors for CVS Health Clinical Trial Services LLC.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1 Table S1.

Classification of top 20 baseline comorbidities.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Eby, E.L., Edwards, A., Meadows, E. et al. Evaluating the relationship between clinical and demographic characteristics of insulin-using people with diabetes and their health outcomes: a cluster analysis application. BMC Health Serv Res 21, 669 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12913-021-06603-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12913-021-06603-0

Keywords