IJPDS International Journal of Population Data Science IJPDS 2399-4908 Swansea University 10.23889/ijpds.v6i1.1414 6:1:1414 S2399490821014142 Population Data Science Estimates of age specific death rates in people with diabetes and associated multimorbidity using Australian administrative pharmaceutical data PurkissShaun Francis1* KeegelTessa12 VallyHassan1 WollersheimDennis1 Department of Public Health, La Trobe University, Bundoora, Victoria, Australia Monash Centre for Occupational and Environmental Health, Monash University, Victoria, Australia Corresponding author: Shaun Francis Purkiss purkiss@latrobe.edu.au

The authors have no conflicts of interest to declare.

 6 1 1414 This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License. This article is available from the IJPDS website at: https://ijpds.org/article/view/1414 Abstract Introduction

Estimating the mortality risk of persons with diabetes can be challenging. Associated conditions such as cardiovascular disease can become the primary cause of mortality and the underlying contribution of diabetes not recorded. Alternative methods to assess mortality risk in people with diabetes would be useful.

 Objective

To evaluate an Australian pharmaceutical database to identify multi-morbidity cohorts associated with diabetes and determine mortality rates in these groups using prescription exchange cessation as a proxy event for death.

 Methods

Australian Pharmaceutical Benefits Scheme data covering the period 2003–14 were used. Persons with diabetes, cardiovascular diseases and dyslipidemia were identified using Anatomic Therapeutic Chemical codes allocated to their recorded dispensed treatments. People with combinations of these conditions were followed and the last recorded prescription exchange used as a proxy event for mortality. Age and gender specific mortality rates and mortality rate ratios for the multi-morbidity cohorts were then calculated from the number of deaths occurring within 10 years.

 Results

346,201 individuals were identified as taking treatments for diabetes, dyslipidemia and cardiovascular conditions in 2004, 86,165 deaths occurred within 10 years of follow up. Overall crude mortality was 26.2/1,000 person years. Age specific mortality rates and rate ratios were calculated for various multi-morbidity groupings. Statin treatments improved the mortality rates associated with diabetes and cardiovascular disease in persons age >54 (Log–Rank <.001).

 Conclusions

Administrative pharmaceutical data can be used to identify persons with diabetes and associated multi-morbidities. Proxy mortality events defined by the cessation of treatment can generate mortality rates, providing an alternative perspective for the assessment of mortality risk.

 Keywords diabetes administrative data multi-morbidity mortality Introduction

Quantifying the mortality risks for people with diabetes can be difficult as diabetes is associated with other chronic diseases that contribute to the risk of death [1–5]. The cause specific mortality for people with diabetes, as recorded on death certificates, may often indicate the primary cause to be from associated conditions such as cardiovascular disease rather than the underlying diabetes [6, 7]. Furthermore, the contribution of diabetes may on occasions not be recorded at all [8, 9]. Accordingly, there are limitations with death certificate data to quantify mortality risk for people with diabetes.

The patterns of mortality for people with diabetes are also changing [1, 10–13]. The prevalence of cardiovascular disease, a major cause of mortality for people with diabetes, is considered to be declining. It is also recognised that the mortality risk of diabetes in elderly people is reducing and approaching similar levels to persons without diabetes [5]. However, cause-specific mortality from other conditions such as cancer may be increasing among people with diabetes [14–16]. Quantifying the relative contributions of age and the associated illnesses and incorporating them into an overall analysis of mortality risk in people with diabetes present challenges [5, 15].

Pharmaceutical administrative data has been recognised as an alternative data source that has utility in the case definition of people with diabetes [17–20]. In this situation, the records of dispensed drugs used in diabetes may act as a proxy clinical diagnosis [19, 20]. This system was developed in Europe using World Health Organization allocated Anatomic Therapeutic Codes to the prescribed treatments established [19, 20]. This method of classification can also be used to identify other chronic diseases associated with diabetes such as cardiovascular disease and dyslipidemia [20]. As a consequence, there is the potential for pharmaceutical data to case define individuals with diabetes and the associated chronic conditions known to have an influence upon mortality [5, 6, 13].

We aim to provide estimates of age specific mortality rates for persons with common comorbidities associated with diabetes, cardiovascular disease and dyslipidemia, using an Australian pharmaceutical data set and cessation of dispensed prescriptions as an implied proxy event for mortality.

 Methods Study design

This study was a longitudinal cohort study using linked administrative pharmaceutical data. Defined cohorts included persons with prevalent diabetes, cardiovascular disease and dyslipidemias alone and in combination. These groups of chronic disease combinations were identified by proxy using the dispensed drugs received. The follow up of individuals was achieved by tracking the records of health service provisions and prescription exchanges within the data set. The date of the last recorded service provision or prescription was used to determine a mortality event and quantify survivals.

 Data source

The pharmacological benefits scheme (PBS) data set published by the Australian Department of Health was employed for this study [21]. The data were selected randomly by government statisticians from the entire current and historical PBS data for the period 2003–14, prior to publication. The information provides an outpatient dispensing perspective of health service users of the Australian universal health-care system which is available to all residents holding a Medicare card [21–23].

This government scheme allows people access to subsidised prescriptions and free health-care in public hospitals. The database provides a comprehensive longitudinal profile of PBS activity for all the people listed and is linked to the Medicare benefits scheme (MBS) data for the same individuals [22, 23].

Supplementary to the local pharmaceutical code of the medication prescribed, associated World Health Organization (WHO) Anatomic Therapeutic Chemical codes were allocated to the class of treatment supplied as an additional field [24].

 Participants

People are de-identified within the PBS data set as part of the confidentialisation methodology employed by the Australian Department of Health. Individuals are assigned a unique identifier (UI) that is maintained across the longitudinal extent of the data. The UI is also employed for database linkage with the MBS data set. Person co-variables listed in the data included age, gender and state of origin.

 Case definition

Persons with diabetes, cardiovascular disorders and dyslipidemias were identified using an ATC system of mapping chronic diseases to the treatments dispensed to individuals (Table 1) [20]. Prevalent chronic disease cases were ascertained by the exchange of at least one prescription of an associated ATC treatment within a calendar year. Persons with prevalent diabetes, cardiovascular disease or dyslipidemia were initially collected. Post data retrieval, persons with multiple co-morbidities were identified by the appearance of their individual UI in the other proxy diagnostic lists at the start of the study. This enabled the construction of multiple chronic disease cohorts from the available combinations of conditions (Table 2).

Chronic condition ATC classification Medication class
Cardiovascular diseases (incl. hypertension) C01, C04A, Anti-hypertensives
C02, C07, Peripheral vasodilators
C08, C09 Beta blocking agents
C02, C07, Calcium channel blockers
C08, C09 Agents acting on the renin-angiotensin system
Diabetes mellitus A10A, A10B, Insulins and analogues
A10X Blood glucose lowering drugs (excl. insulins)
A10A, A10B, Other drugs used in diabetes
Hyperlipidaemia C10 Lipid modifying agents
Multi-morbidity grouping Sex Outcome Mortality rateb (95% CI)
Alive at 10 years Died within 10 years
n Median age (IQR)a n Median age (IQR)a
A, All conditions F 7,560 65 (58–72) 3,704 76 (70–82) 35.9 (34.7, 37.1)
M 7,987 63 (57–70) 4,499 72 (60–84) 40.4 (39.1, 41.6)
B, Diabetes and Cardiovascular F 4,295 62 (53–81) 2,868 79 (73–86) 47.3 (46.0, 48.7)
M 4,060 62 (54–70) 2,750 74 (68–81) 47.6 (46.3, 49.0)
C, Diabetes and Dyslipidemia F 1,517 59 (51–68) 263 70 (62–80) 14.5 (13.8, 15.3)
M 1,663 58 (50–66) 361 68 (60–76) 17.8 (16.9, 18.6)
D, Diabetes alone F 3,968 41 (27-55) 676 72 (60–85) 14.5 (13.7, 15.3)
M 4,158 46 (31–61) 1,093 68 (58–78) 21.4 (20.5, 22.3)
E, Cardiovascular and Dyslipidemia F 34,490 67 (60–74) 11,437 78 (71–85) 25.8 (24.8, 26.8)
M 32,280 63 (56–70) 11,165 74 (68–80) 26.9 (25.9, 27.9)
F, Cardiovascular alone F 64,375 60 (53–67) 24,665 82 (75–89) 29.9 (28.8, 31.0)
M 46,558 58 (57–65) 17,412 75 (67–83) 29.3 (28.2, 30.3)
G, Dyslipidemia alone F 21,863 59 (52–66) 2,259 71 (63–80) 8.9 (8.3, 9.6)
M 25,262 55 (48–62) 3,013 65 (55–75) 10.2 (9.6, 10.9)
All Multi-morbidity groupings Grand Total 260,036 60 (52–68) 86165 75 (62–88) 26.2 (25.1, 27.2)

aData are medians (IQR) in 2004.

bExpressed per 1,000 person years.

 Proxy mortality definition

The most recent date of a record from either the PBS or the MBS dataset was used to define a mortality event. Mortality was assessed by implication provided the event occurred at least 1 year before the end of the longitudinal extent of either dataset on the last day of 2014. The time to this proxy mortality event was calculated in days from the start of 2004. Censored proxy mortality events included all persons that had a last service or prescription exchange during 2014. These data therefore had the potential to include occasional individuals who may have died during that year.

 Data analysis

Persons within allocated multiple morbidity cohorts were identified for the calendar year 2004. Sorting of individuals into single and multiple comorbidity combination cohorts was performed in Excel 2013 (Microsoft, Redmond, WA, USA) and the time to event analysis utilised SPSS statistical software (IBM Corp. Released 2017. IBM SPSS Statistics for Windows, Version 25.0. Armonk, NY: IBM Corp). Log-rank tests were used for comparisons between multiple comorbidities in survival analyses.

 Age specific mortality rates and rate ratio

Age specific mortality rates defined by 10-year age groupings and gender for each of the multiple comorbidity groups were calculated using survival tables generated with SPSS and expressed per 1,000 person years. Mortality rate ratios used published age specific mortality figures from the Australian Bureau of Statistics (2010) for contemporaneous years as the denominator values below the PBS mortality estimates.

 Ethics

Ethics approval for the study was obtained from the La Trobe University ethics committee (S17-198). The construction of this manuscript followed the guidelines recommended in the Reporting of studies conducted using observational routinely-collected health data (RECORD) statement [25].

 Results Participant identification

This study identified 346,201 individuals as having received treatments for diabetes (n = 51,422), cardiovascular disease including hypertension (n = 280,105) and dyslipidemia (n =169,323) during the 2004 calendar year. The mean age of this cohort of 183,940 females (53%) and 162,261 males was 63.3 years (SD 15.1). A total of 86,165 deaths occurred during the follow up period between 2004 and 2014 (33.1%). No causes for the deaths are presented in these results. This analysis displays only the association of chronic conditions, inferred by medications received by persons, ante-mortem, with death implied by the cessation of medical services or dispensed prescription exchanges (Table 2).

 Mortality of multimorbidity cohorts

Within the study population, persons were categorised into seven combinations of multi-morbidity (A-G) and are presented in Table 2. The median age of persons that survived and died within these multi-morbidity cohorts as well as the crude mortality rate for combined age groups are presented in Table 2. Crude mortality overall was 26.2 (CI, 25.1, 27.2) per 1,000 person years.

Mortality rates tended to be higher in persons with cardiovascular disease (Groups A, B, E and F). The lowest crude mortality rates occurred amongst female persons treated for diabetes alone, dyslipidemia alone and diabetes with dyslipidemia (C, D and G). However, the female diabetes alone group was approximately twenty years younger than the other groups, presenting a potential age-related survival advantage (table 2). Of the individuals that died within ten years of observation from 2004, similar average ages at the time of death were noticed for all multi-morbidity subgroupings (median age 75, IQR 62–88 years).

Age specific mortality for females and males are shown in Table 3. These demonstrate a general trend for lower mortality rates amongst females when compared with males in similar age groupings. Mortality rates were age related and increased in both gender subgroups and multi-morbidity subgroupings (Table 3).

Age Interval (years) A All conditions B Diabetes and cardiovascular conditions C Diabetes and dyslipidemia D Diabetes alone
Female
0–14 0 0 0 3.1 (‒)
15–24 26.8 (−)b 0.0 (−) 39.1 (−) 2.4 (1.9, 3.0)
25–34 14.8 (11.9, 17.7) 9.3 (8.3, 10.4) 5.5 (−) 3.5 (3.3, 3.8)
35–44 6.1 (5.2, 6.9) 7.9 (7.3, 8.5) 5.3 (4.3, 6.1) 4.1 (3.9, 4.3)
45–54 7.1 (6.8, 7.4) 9.1 (8.9, 9.3) 5.7 (5.2, 6.2) 6.5 (6.1, 6.8)
55–64 13.0 (12.8, 13.2) 14.3 (14.1, 14.6) 8.0 (7.7, 8.2) 14.3 (13.9, 14.7)
65–74 22.2 (21.9, 22.4) 30.9 (30.8, 31.0) 14.7 (14.3, 15.1) 27.5 (27.3, 27.8)
75–84 48.8 (48.6, 49.1) 67.9 (67.9, 68.0) 36.1 (35.1, 37.1) 63.4 (62.9, 63.8)
85–100 96.9 (96.5, 97.3) 137.4 (136.4, 138.4) 68.0 (56.8, 79.1) 155.2 (148.0, 162.4)
Males
0–14 0.7 (−)
15–24 0 (−) 0 (−) 34.8 (20.4, 49.2) 3.1 (2.9, 3.2)
25–34 13.9 (11.2, 16.7) 11.5 (10.3, 12.6) 20.1 (20.1, 20.1) 7.2 (6.9, 7.5)
35–44 9.5 (9.2, 9.8) 5.8 (5.1, 6.5) 4.1 (3.2, 5.0) 7.5 (7.3, 7.7)
45–54 12.4 (12.2, 12.5) 12.9 (12.6, 13.1) 6.1 (5.8, 6.5) 7.4 (7.1, 7.6)
55–64 16.8 (16.6, 17.0) 20.2 (20.0, 20.5) 13.1 (12.9, 13.3) 22.1 (22.0, 22.2)
65–74 33.4 (33.2, 33.6) 39.6 (39.6, 39.6) 15.0 (14.5, 15.5) 33.0 (32.8, 33.2)
75–84 63.0 (62.6, 63.3) 76.1 (75.8, 76.4) 37.3 (34.6, 40.1) 71.0 (70.4, 71.7)
85–100 124.9 (124.3, 125.5 154.0 (152.2, 155.8) 70.2 (69.2, 81.4) 173.9 (166.9, 180.9)

aExpressed per 1,000 Person years.

b(−) Insufficient data to calculate CI.

The mortality of persons by multi-morbidity sub-groupings and age sub-groupings showed better survival for persons treated with statins across age groupings more than 35–44 years (Table 3).

 Age specific mortality rate ratio of multi-morbidity groups

The mortality rate ratio gradually lowered with increasing age for both genders but was generally higher in males as compared to females (Table 4). The mortality rate ratio in both genders approached unity with increasing age towards the 85–100 age groupings (Figure 1).

Age groups Multi-morbidity grouping AIHW Type 2 diabetes (estimateda) [1]
All conditions Diabetes and cardiovascular Diabetes and dyslipidaemia Diabetes alone
Female
35–44 7.8 (6.6, 9.0) 10.1 (8.5, 11.7) 6.7 (5.6, 7.8) 5.3 (4.5, 6.1) 6.8b
45–54 4.0 (3.4, 4.6) 5.1 (4.3, 5.9) 3.2 (2.7, 3.7) 3.6 (3.0, 4.2) 2.5
55–64 3.3 (2.8, 3.8) 3.7 (3.1, 4.3) 2.0 (1.7, 2.3) 3.6 (3.0, 4.2) 2.1
65–74 2.2 (1.8, 2.6) 3.1 (2.6, 3.6) 1.5 (1.2, 1.7) 2.8 (2.4, 3.2) 1.9
75–84 1.5 (1.3, 1.7) 2.1 (1.8, 2.4) 1.1 (0.9, 1.3) 2.0 (1.7, 2.3) 1.3
85–100 0.6 (0.5, 0.7) 0.8 (0.7, 0.9) 0.4 (0.3, 0.5) 0.9 (0.7, 1.0) 1.0
Males
35–44 6.7 (5.6, 7.8) 4.1 (3.4, 4.8) 2.9 (2.4, 3.4) 5.3 (4.5, 6.1) 5.7b
45–54 4.2 (3.5, 4.9) 4.4 (3.7, 5.1) 2.1 (1.8, 2.4) 2.5 (2.1, 2.9) 2.4
55–64 2.5 (2.1, 2.9) 3.1 (2.6, 3.6) 2.0 (1.7, 2.3) 3.3 (2.8, 3.8) 2.0
65–74 2.1 (1.8, 2.4) 2.5 (2.1, 2.9) 0.9 (0.8, 1.0) 2.1 (1.8, 2.4) 1.8
75–84 1.3 (1.1, 1.5) 1.6 (1.3, 1.9) 0.8 (0.7, 0.9) 1.5 (1.3, 1.7) 1.5
85–100 0.6 (0.5, 0.7) 0.8 (0.7, 0.9) 0.4 (0.3, 0.5) 0.9 (0.8, 1.0) 1.1

aEstimates extracted from figure data points.

bAIHW age group <45 years.

 Multi-morbidity combinations include A, diabetes, cardiovascular disease and dyslipidemia; B, diabetes and cardiovascular disease; C, diabetes and dyslipidemia; D, diabetes only; E, cardiovascular disease and dyslipidemia; F, cardiovascular disease only; G Dyslipidemia only.

A comparison of mortality rate ratios from this study using PBS data for diabetes (types 1 and 2) and similar AIHW data for type 2 diabetes demonstrates reasonable correlation and accuracy for females. However, the PBS derived mortality ratio for males aged 45-75 are higher than those reported in the AIHW study (Figure 2) [1].

 Discussion

People with diabetes are considered to have a shorter life expectancy as compared to those without diabetes [1]. This increased risk of mortality is largely attributable to cardiovascular related causes, a common co-morbidity associated with diabetes [1, 3, 5]. This propensity of individuals with diabetes to develop multiple co-morbidities complicates the understanding of the impact of diabetes upon survival. In a report on deaths among Australian people with diabetes for the period 2009–14, it was recognised that 56% of people did not have diabetes recorded on the death certificate as a cause of death [1]. As a consequence, the influence of diabetes on mortality may be biased by these death registration anomalies.

This study demonstrates the utility of a pharmaceutical based system to identify individuals with diabetes and associated conditions and determine the extent to which the combination of these conditions can influence survival. The quasi-experimental design of the study although using data that is retrospective, follows people forwards in time, and defines death by proxy using the last date of supply for a prescription or a medical service. As a consequence, the identification of people with diabetes is made ante-mortem which contrasts with mortality data arising from death certificates.

The proxy mortality event marker used in this study has been previously applied in clinical experiments and was demonstrated to be of particular utility in persons age >45 [26]. Furthermore, the proxy marker has been useful in the assessment of long-term outcomes following surgical procedures [27, 28]. In this study, reasonable correlation and accuracy of mortality rate ratios was shown in females with similar data from the AIHW [1]. However, males in the age groups 45–74 demonstrated higher PBS derived mortality rate ratios in general. The reason for this disparity is unclear but death registration anomalies reported in AIHW data may be contributory.

The survival estimates using PBS derived data have demonstrated consistent trends. The survival of females compared to males is generally better across all age groupings and multi-morbidity combinations. The addition of statin treatments for dyslipidemia improves the mortality rate and mortality rate ratio of persons with diabetes alone and persons with diabetes and associated cardiovascular disease. This effect is seen across both genders and within the more elderly age groupings.

There are limitations in PBS and the AIHW derived mortality estimates. PBS derived mortality estimates are based upon the time of cessation of medical services or prescription exchanges. As a consequence, it is an implied event and mortality is assumed by proxy. Erroneous mortality events may occur with PBS estimates, when for instance, a person leaves Australia, or discontinues medical treatments for whatever reason. AIHW reports in contrast, determine mortality by probabilistic linkage with mortality databases. This process will likely introduce small inherent errors which when combined with the death registration anomalies will likely reduce the accuracy of AIHW mortality estimates.

PBS data is generated from a network of community pharmacies from a wide geographic area with good rural coverage [21, 22]. In contrast, AIHW reports acknowledge that the population coverage may be lower in remote areas of Australia and as a consequence their data may not fully represent this portion of society [1]. PBS and AIHW derived mortality statistics also exclude persons who are able to manage their diabetes without medication with presumed lifestyle changes [1]. As a consequence, this cohort is not included in the calculation of mortality from both data sources.

A consistent trend identified in both AIHW data and PBS derived mortality rate ratios is the improvement with age for both genders [1]. The mortality of the most elderly cohorts with diabetes have similar mortality rate ratios as compared to the general population. This is likely the consequence of the increasing values of the denominator, the rate of mortality for the general population with age having a greater influence on the mortality rate ratio overall. In contrast, the higher mortality rate ratios noticed in the younger age groupings are likely influenced by the lower mortality rates of these cohorts in the general population.

PBS mortality data provides an alternative perspective of diabetes risk. These data could be used alone or in combination with other data using techniques of triangulation to provide improved mortality estimates overall [29].

 Conclusions

The mortality analyses presented in this study identified people with multiple chronic diseases using the medications they received. The categories and mortality events, also defined by implication from medication profiles focused upon the conditions of diabetes, cardiovascular disease with hypertension and dyslipidemia alone and in combination. The potential to assess other chronic disease combinations can be considered, and may be of use in actuarial analyses, generation of life expectancy tables and the prediction of health-care utilisation costs. The technique also has the potential to explore many combinations of other chronic conditions providing they have a therapeutic signature. As a consequence, bespoke mortality risk profiles could be generated for individual persons based upon their profile of dispensed medication.

We consider that the medication received by an individual and recorded within pharmaceutical administrative data can identify chronic disease multimorbidity cohorts and may be a useful tool to estimate mortality risk of persons with combinations of diabetes, cardiovascular disease and dyslipidemia.

 Acknowledgements

We are grateful for the authorisation of the study by the La Trobe research ethics committee. (Approval number S17-177).

 References Australian Institute of Health and Welfare 2017. Deaths among people with diabetes in Australia, 2009–2014. Cat. no. CVD 79. Canberra: AIHW. Available from: https://www.aihw.gov.au/getmedia/c164ec7c-fc66-4991-8e53-62eba274ead4/aihw-cvd-79.pdf.aspx?inline=true Sijbrands EJ, Tornij E, Homsma SJ. Mortality risk prediction by an insurance company and long-term follow-up of 62,000 men. PloS one. 2009 May 6;4(5):e5457. 10.1371/journal.pone.0005457 Mathers CD, Lopez AD, Murray CJ. The burden of disease and mortality by condition: data, methods and results for 2001. Global burden of disease and risk factors. 2006 Jan 1;45:88. Leal J, Gray AM, Clarke PM. Development of life-expectancy tables for people with type 2 diabetes. European heart journal. 2009 Apr 1;30(7):834–9. 10.1093/eurheartj/ehn567 Harding, J. L., Shaw, J. E., Peeters, A., Guiver, T., Davidson, S., & Magliano, D.J. (2014). Mortality trends among people with type 1 and type 2 diabetes in Australia: 1997. 10.2337/dc14-0096 McEwen LN, Kim C, Haan M, Ghosh D, Lantz PM, Mangione CM, Safford MM, Marrero D, Thompson TJ, Herman WH. Diabetes reporting as a cause of death: results from the Translating Research Into Action for Diabetes (TRIAD) study. Diabetes care. 2006 Feb 1;29(2):247–53. 10.2337/diacare.29.02.06.dc05- Goldacre MJ. Cause-specific mortality: understanding uncertain tips of the disease iceberg. Journal of Epidemiology & Community Health. 1993 Dec 1;47(6):491–6. 10.1136/jech.47.6. Cheng WS, Wingard DL, Kritz-Silverstein D, Barrett-Connor E. Sensitivity and specificity of death certificates for diabetes: as good as it gets?. Diabetes care. 2008 Feb 1;31(2):279–84. 10.2337/dc07-1327 Australian Institute of Health and Welfare 2015. Leading cause of premature mortality in Australia fact sheet: diabetes. Cat. no. PHE 200. Canberra: AIHW. ISBN 978-1-74249-800-3 (PDF) https://www.aihw.gov.au/getmedia/c5e584b2-06bd-4d80-9300-6581e640836f/phe200-diabetes.pdf.aspx Zimmet PZ, Alberti KG. The changing face of macrovascular disease in non-insulin-dependent diabetes mellitus: an epidemic in progress. The Lancet. 1997 Jul 1;350:S1–4. 10.1016/S0140-6736(97)90020-9 Wild SH, Roglic G, Green A, Sicree R, King H. Global prevalence of diabetes: estimates for the year 2000 and projections for 2030: response to Rathman and Giani. Diabetes care. 2004 Oct 1;27(10):2569. 10.2337/diacare.27.5.1047 Gregg EW, Sattar N, Ali MK. The changing face of diabetes complications. The lancet Diabetes & endocrinology. 2016 Jun 1;4(6):537–47. 10.1016/s2213-8587(16)30010-9 Austalian Institute of Health and Welfare. Australian Burden of Disease Study: Impact and Causes of Illness and Death in Australia 2015. Available online: https://www.aihw.gov.au/reports/burden-of-disease/burdendisease-study-illness-death-2015-summary/contents/summary (accessed on 13 October2019). Diabetes: the silent pandemic and its impact on Australia (2012). Edited by Associate Professor Jonathan Shaw, Associate Director, Baker IDI Heart and Diabetes Institute and Stephanie Tanamas, Epidemiologist, Baker IDI Heart and Diabetes Institute, with input from Diabetes Australia and Juvenile Diabetes Research Foundation (JDRF). Available from https://static.diabetesaustralia.com.au/s/fileassets/diabetes-australia/e7282521-472b-4313-b18e-be84c3d5d907.pdf Harding JL, Shaw JE, Peeters A, Cartensen B, Magliano DJ. Cancer risk among people with type 1 and type 2 diabetes: disentangling true associations, detection bias, and reverse causation. Diabetes care. 2015 Feb 1;38(2):264–70. 10.2337/dc15-er04a Giovannucci E, Harlan DM, Archer MC, Bergenstal RM, Gapstur SM, Habel LA, Pollak M, Regensteiner JG, Yee D. Diabetes and cancer: a consensus report. Diabetes care. 2010 Jul 1;33(7):1674–85. 10.3322/caac.20078 Clarke P, Leal J, Kelman C, Smith M, Colagiuri S. Estimating the cost of complications of diabetes in Australia using administrative health-care data. Value in health. 2008 Mar 1;11(2):199–206. 10.1111/j.1524-4733.2007.00228.x Asghari S, Courteau J, Carpentier AC, Vanasse A. Optimal strategy to identify incidence of diagnostic of diabetes using administrative data. BMC medical research methodology. 2009 Dec 1;9(1):62. 10.1186/1471-2288-9-62 Chini F, Pezzotti P, Orzella L, Borgia P, Guasticchi G. Can we use the pharmacy data to estimate the prevalence of chronic conditions? a comparison of multiple data sources. BMC Public Health. 2011 Dec;11(1):1–8. 10.1186/1471-2458-11-688 Huber CA, Szucs TD, Rapold R, Reich O. Identifying patients with chronic conditions using pharmacy data in Switzerland: an updated mapping approach to the classification of medications. BMC public health. 2013 Dec;13(1):1030. 10.1186/1471-2458-13-1030 Mellish L, Karanges EA, Litchfield MJ, Schaffer AL, Blanch B, Daniels BJ, Segrave A, Pearson SA. The Australian Pharmaceutical Benefits Scheme data collection: a practical guide for researchers. BMC research notes. 2015 Dec;8(1):634. 10.1186/s13104-015-1616-8 Moles RJ, Stehlik P. Pharmacy practice in Australia. The Canadian journal of hospital pharmacy. 2015 Sep;68(5):418. 10.4212/cjhp.v68i5.1492 Benrimoj SI, Frommer MS. Community pharmacy in Australia. Australian Health Review. 2004;28(2):238–46. 10.1071/AH040238 World Health Organization. (2006). WHO Collaborating Centre for Drug Statistics Methodology: ATC classification index with DDDs and Guidelines for ATC classification and DDD assignment. Oslo, Norway: Norwegian Institute of Public Health. Benchimol EI, Smeeth L, Guttmann A, Harron K, Moher D, Petersen I, Sørensen HT, von Elm E, Langan SM, RECORD Working Committee. The REporting of studies Conducted using Observational Routinely-collected health Data (RECORD) statement. PLoS medicine. 2015 Oct 6;12(10):e1001885. 10.1371/journal.pmed.1001885 Purkiss S, Keegel T, Vally H, Wollersheim D. Estimates of age-specific death rates and mortality risk using administrative pharmaceutical data. International Journal of Population Data Science. 2019 Nov 7;4(3). 10.23889/ijpds.v4i3.1165 Purkiss S, Keegel T, Vally H, Wollersheim D. Long-term survival following successful abdominal aortic aneurysm repair evaluated using Australian administrative data. ANZ Journal of Surgery. 2020 Mar;90(3):339–44. 10.1111/ans.15598 Purkiss S, Keegel T, Vally H, Wollersheim D. Cervical re-explorations and proxy survival following parathyroidectomy for primary hyperparathyroidism using Australian administrative data. ANZ Journal of Surgery. 2020 May;90(5):872–6. 10.1111/ans.15238 Morse JM. Approaches to qualitative-quantitative methodological triangulation. Nursing research. 1991 Mar 1;40(2):120–3. 10.1097/00006199-199103000-00014 Abbreviations AIHW Australian Institute of Health and Welfare ATC Anatomic Therapeutic Chemical Codes CVS Cardiovascular System IQR Interquartile Range MBS Medicare Benefits Scheme PBS Pharmaceutical Benefits Scheme SPSS Statistical Package for the Social Sciences UI Unique Identifier WHO World Health Organization