- Research article
- Open Access
- Open Peer Review
Time series modelling to forecast prehospital EMS demand for diabetic emergencies
BMC Health Services Research volume 17, Article number: 332 (2017)
Acute diabetic emergencies are often managed by prehospital Emergency Medical Services (EMS). The projected growth in prevalence of diabetes is likely to result in rising demand for prehospital EMS that are already under pressure. The aims of this study were to model the temporal trends and provide forecasts of prehospital attendances for diabetic emergencies.
A time series analysis on monthly cases of hypoglycemia and hyperglycemia was conducted using data from the Ambulance Victoria (AV) electronic database between 2009 and 2015. Using the seasonal autoregressive integrated moving average (SARIMA) modelling process, different models were evaluated. The most parsimonious model with the highest accuracy was selected.
Forty-one thousand four hundred fifty-four prehospital diabetic emergencies were attended over a seven-year period with an increase in the annual median monthly caseload between 2009 (484.5) and 2015 (549.5). Hypoglycemia (70%) and people with type 1 diabetes (48%) accounted for most attendances. The SARIMA (0,1,0,12) model provided the best fit, with a MAPE of 4.2% and predicts a monthly caseload of approximately 740 by the end of 2017.
Prehospital EMS demand for diabetic emergencies is increasing. SARIMA time series models are a valuable tool to allow forecasting of future caseload with high accuracy and predict increasing cases of prehospital diabetic emergencies into the future. The model generated by this study may be used by service providers to allow appropriate planning and resource allocation of EMS for diabetic emergencies.
The prevalence of diabetes is increasing  and is becoming one of the most significant health issues of the developed world . In Australia 1.24 million people are registered as having diabetes  and an additional 500,000 are estimated to have undiagnosed diabetes . Between 2013 and 2016, this number increased by 4.3% per year with the increase in the prevalence of type 2 diabetes greater (4.4% increase per year) than that of type 1 diabetes (0.6% increase per year) . The impact of diabetes on the health system is widespread and includes both acute and chronic complications. Patients experiencing acute glycemic complications, most frequently severe hypoglycemia and hyperglycemic crises (diabetic ketoacidosis and hyperglycemic hyperosmolar state), will often seek emergency medical assistance in the community from the prehospital Emergency Medical Services (EMS) . When combined with the rising demand for prehospital EMS  , the increasing diabetes prevalence will require careful resource planning.
Health forecasting, the prediction of health or disease episodes and signaling of future events, is increasingly recognized as a valuable tool to facilitate health service provision and resource allocation. Health forecasting is often based on a time series, a sequence of data points collected at successive, equally spaced time intervals, which may be characterised by trend, seasonality, cyclicality and randomness . A time series provides a statistical means of temporal trend analysis and prediction of future events based on observed values. There are various forecasting methods in the health literature, including historical averaging , smoothing techniques [10, 11], linear regression  and autoregressive integrated moving average (ARIMA) modelling [13, 14]. ARIMA modelling has demonstrated successful prediction of a range of specific health disease events [15, 16] as well as in aggregate caseload, such as hospital emergency department [9–11, 13] and prehospital EMS  and has been recognised for its simplicity and ease of administration .
Temporal trends in diabetic emergencies have been investigated at the hospital level   but trends of demand for prehospital EMS for acute diabetic emergencies have not been reported. The aims of this study were 1) to quantify the temporal trends in utilization of prehospital EMS for acute diabetic emergencies, 2) to model the temporal variation of prehospital diabetic emergencies and 3) to use the model to make short term predictions of future EMS demand for diabetic emergencies.
A time series analysis on monthly case rates of diabetic emergencies attended by Ambulance Victoria (AV) between January 2009 and December 2015 was conducted. AV is a two-tiered, prehospital EMS system and the sole provider of prehospital EMS for the state of Victoria, Australia. Individuals of all ages receiving prehospital care from AV during the study period with a documented final primary assessment of “hyperglycemia” or “hypoglycemia” were included. The final primary assessment, as assigned by the attending paramedic, is defined as the main problem at the time the patient is discharged from EMS care. No blood glucose level threshold parameters were imposed, however AV uses a BGL < 4 mmol/L (<72 mg/dl) to treat for hypoglycemia and does not specify a glycemic threshold for hyperglycemia. Every case attended by AV is recorded by the attending paramedic, using the VACIS®, an electronic patient care record and integrated data warehouse . Diabetes type, based on patient/bystander self-report, was classified as type 1 diabetes, type 2 diabetes or unspecified diabetes type/status. In this study de-identified data was used with no ability to distinguish repeat callers, thus repeat attendances were treated as individual cases. Data spanning 4 calendar months (September 2014 – December 2014) were unavailable due to lapse in electronic data collection linked to industrial action. The Monash Health Human Research Ethics Committee approved this study.
Descriptive analysis was conducted by tabulating total annual case count and median monthly case count for each year. Categorical variables (gender, diabetes type, emergency type) are reported as absolute number and percentage, and differences between subgroups analysed using χ2 test. Age was summarised as median with interquartile range (IQR), and differences between subgroups analysed using Kruskal-Wallis test.
The autoregressive integrated moving average (ARIMA) modelling process was conducted using monthly case counts during the study period, initially with total case count and then separately for hypoglycemia/hyperglycemia and for male/female. ARIMA models make use of previous observations to make predictions of future values using lag parameter values, under the assumption that the pattern will persist. Lags of the differenced series, termed Auto Regressive (AR), indicate the strength of relationship between incidence rates in consecutive months and lags of the forecast errors, termed Moving Average (MA), check for dependence of monthly incidence on current and past model residuals. A differencing term, D, is applied to make the data stationary when the time series displays a long term trend and a seasonal term, S, is incorporated when the time series displays a seasonal pattern, and Seasonal ARIMA (SARIMA) modelling is used.
The data was divided into two sets; training data and validation data. This was because models are expected to perform well on the data from which they were derived, so post sample validation on a prospective ‘new’ dataset was performed. Model generation was based on the dataset; January 2009 - August 2014 (training dataset) and model validation was based on the dataset; January 2015 - December 2015 (validation dataset). The Box-Jenkins approach was used to fit the models . Annual seasonality was apparent in the plotted monthly caseload (Fig. 2), thus SARIMA modelling with a 12-month seasonality term was used. The Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) were plotted to examine stationarity and lags and to assist in identification of the order the MA and AR terms of each model.
Different formulations of the AR and MA terms were modelled. The following measures of prediction accuracy were calculated; Mean Absolute Error (MAE), Mean Square Error (MSE) and Mean Absolute Percentage Error (MAPE), defined as:
where Oi, Pi and n are the observed and predicted counts for month “i” and the number of observations, respectively and where a lower prediction error indicates a better model fit. Each model was used to predict values for the training dataset and the validation dataset. The most parsimonious model with lower MAE, MSE and MAPE values for each condition (overall caseload, hypoglycaemia, hyperglycaemia, female and male) was selected. The final selected model for each condition was then used to generate predicted monthly caseload values until the December 2017. The predicted values for the training data were based on a one-step forecast while the predicted values for the validation dataset and the 2017 forecast were based on dynamic forecasting, commencing on August 2014, the final month of training dataset. Compared to one-step forecasting, which bases forecasts on observed values of the preceding month, dynamic forecasts are based on predicted values of the preceding month, thus are susceptible to errors accumulating over time.
The chosen SARIMA model for overall caseload was then used to generate 1-month forecasts, 3-month forecasts and 12-month forecasts, and the MAE, MSE and MAPE were calculated for each, using dynamic forecast predictions. In addition, the chosen SARIMA model was compared to other modelling techniques for overall caseload. The models generated for comparison were; a) SARIMA, b) SARIMA + time trend, c) ARIMA + seasonality, d) exponential smoothing and e) linear time trend + seasonality. Modelling was performed using one-step forecasts and model accuracy was compared using the three accuracy measures; MAE, MSE and MAPE. All analyses were performed using Stata software version 14.0 (StataCorp, Texas, USA) and the level of significance was set at 5%.
During the 7-year study period 41,454 prehospital diabetic emergencies were attended by AV, equating to an annual mean (±SD) of 5922 ± 508 cases. Table 1 displays the total annual case load for each year, median age and the proportional distribution of gender, type of diabetes and type of emergency. The median [IQR] monthly attendance over the study period was 506 [481, 534], with an increasing trend in median monthly attendance per year between 2009 (484.5 [447.5, 540.5]) and 2015 (549.5 [513.5, 581.5]) (Fig. 1). Gender distribution (approximately 55% male) remained relatively stable and median [IQR] age decreased very slightly from 60 [40, 76] years to 59 [39, 75] years over the study period (p = 0.11). Overall, hypoglycemia accounted for 70% of attendances however, over the study period, the proportion of attendances for hypoglycemia decreased (75% to 63%) and the proportion of hyperglycemia increased (25 to 37%) (p < 0.001). Attendances to people with type 1 diabetes and type 2 diabetes accounted for 48% and 39% of overall attendances, respectively. The relative proportion of attendances to persons with type 1 diabetes and type 2 diabetes decreased (53% to 44%) and increased (35 to 42%), respectively (p < 0.001).
The monthly caseload of prehospital diabetic emergencies (Fig. 2) showed possible seasonality, with peaks apparent around December-January and troughs around April-May, as well as an overall increasing trend, suggesting a differencing term be incorporated in the model. The partial autocorrelation plot showed spikes at approximately month 12 and month 24 and the autocorrelation plot demonstrated an initial significant spike with a rapid decline and a peak at 12 months, both indicating the inclusion of a seasonality term (Fig. 3). The modelling process (Table 2) found the SARIMA (0,1,0) (0,1,0,12) term provided the best fit for the overall case counts, generating a MAE of 23.0, a MSE of 758.4 and a MAPE of 4.2% for the validation dataset. The predicted values of the final model (Fig. 4a, heavy purple line) demonstrate the increasing trend seen in the earlier data and predict a monthly caseload of 743 by the end of 2017, representing a 17% increase in monthly caseload from December 2015.
The modelling process found the SARIMA (0,1,1) (0,1,1,12) term provided the best fit for cases of hypoglycemia, generating a MAE of 25.5, a MSE of 941.2 and a MAPE of 7.6% for the validation dataset. The predicted values of the final model (Fig. 4a, heavy blue line) suggest a stable trend of EMS-attended hypoglycemia cases in the short-term.
The modelling process found the SARIMA (0,1,1) (0,1,1,12) provided the best fit for cases of hyperglycemia, generating a MAE of 12.9, a MSE of 254.3 and a MAPE of 6.3% for the validation dataset. The predicted values of the final model (Fig. 4a, heavy orange line) suggest a rise in the number of EMS-attended cases of hyperglycemia, with a monthly caseload of 290 predicted by the end of 2017, representing a 45% increase from December 2015.
The modelling process found the SARIMA (2,1,0) (2,1,0,12) provided the best fit for female cases generating a MAE of 14.8%, a MSE of 360.3 and a MAPE of 6.1%. The predicted values of the final model (Fig. 4b, green line) suggest a rise in caseload of 24% between December 2015 and December 2017.
The modelling process found the SARIMA (1,1,1) (1,1,1,12) provided the best fit for male cases of diabetic emergencies generating a MAE of 18.9%, a MSE of 618.4 and a MAPE of 6.0%. The predicted values of the final model (Fig. 4b, golden line) suggest a very small rise in caseload of less than 1% between December 2015 and December 2017.
Comparison of forecast horizon
The SARIMA (0,1,0,12) model was used to generate 1-month, 3-month and 12-month forecasts (Table 3). The shorter forecast horizon (1 month) demonstrated greatest accuracy (MAPE of 7.1%), and the 12-month projection demonstrated slightly greater accuracy than the 3-month projection, with a MAPE of 8.6% and 9.0%, respectively.
Comparison of time series models
When compared to other types of time series models, SARIMA (0,1,0,12) (one-step) generated a MAPE of 7.3%, out-performing b) SARIMA+ time trend (MAPE = 9.1%), d) exponential smoothing (MAPE = 8.8%) and e) linear time trend + seasonality (MAPE = 7.6%), and was comparable to c) ARIMA + seasonality (MAPE = 7.2%) (Table 4).
Utilization of prehospital EMS for acute diabetic emergencies is increasing with increased prevalence of diabetes. After separate analysis of the type of emergency, the majority of this increase appears to be due to cases of hyperglycemia. The SARIMA modelling process was able to model the monthly incidence of prehospital diabetic emergencies with good accuracy and predicts an increased caseload in the short term.
The temporal pattern of EMS demand for diabetic emergencies exhibited a long term rising trend with a seasonal component, where peaks were apparent around December-January (summer months) and troughs around April-May (autumn months). This is consistent with documented seasonal fluctuation in HbA1c levels showing lower levels in the summer months [23, 24]. Lower HbA1c levels are associated with an increased risk of severe hypoglycemia in people with type 1  and type 2 diabetes . Seasonal fluctuation has also been reported in hospital attendance and admission for hyperglycemia in people with type 2 diabetes . The recognition of these season-related patterns may assist short-term organizational resource allocation, such as the increased service provision in the peak months.
While knowledge of seasonal fluctuations are useful for short-term service provision and rostering, the underlying increase of EMS use for diabetic emergencies throughout the study period demonstrates a need to address service provision longer-term. Cases of hypoglycemia and people with type 1 diabetes accounted for the majority of attendances, however, their relative proportions declined, with greater increases in attendances to cases of hyperglycemia and people with type 2 diabetes. This shifting demographic in demand for prehospital EMS has significant implications for the system in light of the types of patients who require transport to hospital. Local transport rates for people presenting with hyperglycemia are approximately 90%, much higher than those presenting with hypoglycemia (40%) . Given transport to hospital consumes EMS resources for longer periods of time, and the types of patients who generally require hospital transport is increasing, consideration to long-term resource planning and sustainable alternatives to hospital transport are warranted. The forecast increases in EMS-attended diabetic emergencies highlights the importance of increasing routine detection of diabetes at the primary-care level and before emergency intervention is required as well as improved community-based management of diabetes, to avert some of the future caseload.
The International Diabetes Federation  reports that for the Western Pacific region (which includes Australia), the number of people with diabetes is estimated to increase from 153 million in 2015 to 215 million in 2040, representing a 1.6% annual increase in community prevalence of diabetes over the time period. Our model predicts a monthly caseload of 743 by the end of 2017, representing a 17% increase from December 2015, and an annual increase of approximately 8%, which far outweighs the predicted increase in community diabetes prevalence. Furthermore, between 2013 and 2015, the Australian National Diabetes Services Scheme register reported annual increases of prevalence of people with type 1 and type 2 diabetes of 0.6% and 4.4%, respectively. Over this same time period, there was an increase in EMS attendances for diabetic emergencies to people with type 1 and type 2 diabetes, suggesting the increases in EMS caseload may be partially explained by the increasing prevalence of diabetes, with the future predicted health burden likely to extend to prehospital EMS.
As predicted, using the 1-month forecast horizon resulted in greater accuracy than the 3-month and 12-month forecast horizon, however, the 12-month horizon performed better than the 3-month horizon. This could be explained by the inclusion of the 12-month seasonality term to the SARIMA models. Across model comparison showed that the selected SARIMA model resulted in a lower MAPE than SARIMA+ time trend, d) exponential smoothing and linear time trend + seasonality and a MAPE comparable to ARIMA + seasonality. While both SARIMA and ARIMA + seasonality are conceptually similar, SARIMA was selected due to simplicity in formulation and execution. It should be highlighted, however, that this is not a formal comparison of various time series models but rather a selection of the most appropriate model for the dataset.
The SARIMA (0,1,0) (0,1,0,12) provided the best fit for the overall caseload, generating a MAPE of 4.2%. Given the lack of published health forecast models in the prehospital field, and the variability of modelling approaches and methods of evaluation, comparison of our model’s performance is difficult. The MAE, MSE and MAPE were used in this study, however various other indices, such as route mean squared error (RMSE)  or goodness of fit criterion  have been used to evaluate performance of SARIMA models of emergency department caseloads. One study  developed an ARIMA model for forecast of daily rates of attendance at an emergency department and found a seasonal ARIMA (0,1,1)(1,0,1) yielded a MAPE of 4.8%, comparable to our results.
The use of time series modelling in health care to predict future events is increasing. Although there are many methods of time series prediction modelling, the ARIMA/SARIMA method has some distinct advantages. A study  comparing simple regression methods with ARIMA methods to predict acute hospital presentations argued that in the acute hospital environment, where trends can change suddenly and unpredictably, simple regression methods were not as effective forecasting tools as ARIMA, which can accommodate for autocorrelated data. In disease surveillance, ARIMA modelling was found to be favorable when compared to Bayesian processes to predict weekly incidence of Dengue fever . Prediction models were based on weekly data from 2001–2006, and validated models on data between 2007 and 2008. Authors found that while the Bayesian K-H model provided marginally better prediction performance, this was outweighed by the relative ease of execution of ARIMA modelling.
The VACIS dataset captures all prehospital cases in the state of Victoria, Australia. The study is further strengthened by the prolonged timeframe of the data capture period (84 months), considered to be in the range of required observations for optimization (n = 50–100) . Limitations of this study include the four months of missing data (Sept-Dec 2014), which impacted ability to generate 12-month forecasts for the validation period (2015) as well as the inherent limitations of forecast modelling. ARIMA models reflect and extend on past patterns. Thus, forecast accuracy is inversely related to the length of time of the forecast, whereby accuracy decreases with increasing length of forecast horizon as vulnerability to environmental or resourcing changes not accounted for in the model increases. Forecasting is also reliant on the reliability of the health data and robustness of the forecasting technique and accuracy is improved by updating the models as more data becomes available. A further limitation is that only cases listing hypoglycemia or hyperglycemia in the primary assessment were included, possibly underestimating the magnitude of the issue (as cases listing hypoglycemia or hyperglycemia as a secondary assessment were not included). The generalizability of the findings to other regions of the world may also be impacted by the local health care model and the socioeconomic and cultural background of those regions.
In conclusion, ARIMA time series models forecast future prehospital caseload of diabetic emergencies with high accuracy. The results of this study demonstrate the requirement of appropriate health services planning and resource provision into the future with increases in both diabetes prevalence and pressures on prehospital emergency medical services anticipated. Further research regarding the effects of specific variables such as community prevalence of diabetes and climate factors, possible precipitants of diabetic emergencies and potential preventative measures to ease demand are required.
Emergency medical services
Autoregressive integrated moving average
Seasonal autoregressive integrated moving average
Auto correlation function
Partial autocorrelation function
Mean absolute error
Mean square error
Mean absolute percentage error
Magliano DJ, Peeters A, Vos T, Sicree R, Shaw J, Sindall C, Haby M, Begg SJ, Zimmet PZ. Projecting the burden of diabetes in Australia - what is the size of the matter? Aust N Z J Public Health. 2009;33(6):540–3.
Tabish SA. Is Diabetes Becoming the Biggest Epidemic of the Twenty-first Century? Int J Health Sci. 2007;1(2):V-VIII.
National Diabetes Service Scheme: Australian Diabetes Map. 2016. Available from http://www.diabetesmap.com.au/#/. Last updated 31/3/16. Accessed 15 Apr 17.
Diabetes Australia [Homepage on the internet] Canberra: Diabetes Australia: Last update 2014, cited Jan 2016. Available from https://www.diabetesaustralia.com.au/diabetes-in-australia.
Villani M, Nanayakkara N, Ranasinha S, Tan C, Smith K, Morgans A, Soldatos G, Teede H, Zoungas S. Utilisation of emergency medical services for severe hypoglycaemia: An unrecognised health care burden. J Diabetes Complicat. 2016;30(6):1081–6.
Lowthian JA, Cameron PA, Stoelwinder JU, Curtis A, Currell A, Cooke MW, McNeil JJ. Increasing utilisation of emergency ambulances. Aust Health Rev. 2011;35(1):63–9.
Toloo S, FitzGerald G, Aitken P, Ting J, Tippett V, Chu K. Emergency Health Services: Demand and Service Delivery Models. Monograph 1: Literature Review and Activity Trends. In: Queensland University of Technology. 2011.
Soyiri IN, Reidpath DD. An overview of health forecasting. Environ Health Prev Med. 2013;18(1):1–9.
Schweigler LM, Desmond JS, McCarthy ML, Bukowski KJ, Ionides EL, Younger JG. Forecasting Models of Emergency Department Crowding. Acad Emerg Med. 2009;16(4):301–8.
Boyle J, Jessup M, Crilly J, Green D, Lind J, Wallis M, Miller P, Fitzgerald G. Predicting emergency department admissions. Emerg Med J. 2012;29(5):358–65.
Champion R, Kinsman LD, Lee GA, Masman KA, May EA, Mills TM, Taylor MD, Thomas PR, Williams RJ. Forecasting emergency department presentations. Aust Health Rev. 2007;31(1):83–90.
Takase M, Carlin J. Modelling seasonal variations in presentations at a paediatric emergency department. Hiroshima J Med Sci. 2012;61(3):51–8.
Sun Y, Heng BH, Seow YT, Seow E. Forecasting daily attendances at an emergency department to aid resource planning. BMC Emerg Med. 2009;9(1):1–9.
Earnest A, Chen MI, Ng D, Sin LY. Using autoregressive integrated moving average (ARIMA) models to predict and monitor the number of beds occupied during a SARS outbreak in a tertiary hospital in Singapore. BMC Health Serv Res. 2005;5:36.
Wah W, Das S, Earnest A, Lim LKY, Chee CBE, Cook AR, Wang YT, Win KMK, Ong MEH, Hsu LY. Time series analysis of demographic and temporal trends of tuberculosis in Singapore. BMC Public Health. 2014;14(1):1–10.
Medina DC, Findley SE, Guindo B, Doumbia S. Forecasting Non-Stationary Diarrhea, Acute Respiratory Infection, and Malaria Time-Series in Niono, Mali. PLoS ONE. 2007;2(11), e1181.
Channouf N, L’Ecuyer P, Ingolfsson A, Avramidis AN. The application of forecasting techniques to modeling emergency medical system calls in Calgary, Alberta. Health Care Manag Sci. 2007;10(1):25–45.
Earnest A, Tan SB, Wilder-Smith A, Machin D. Comparing Statistical Models to Predict Dengue Fever Notifications. Comput Math Methods Med. 2012;2012:6.
Lipska KJ, Ross JS, Wang Y, Inzucchi SE, Minges K, Karter AJ, Huang ES, Desai MM, Gill TM, Krumholz HM. National trends in US hospital admissions for hyperglycemia and hypoglycemia among Medicare beneficiaries, 1999 to 2011. J Am Med Assoc Int Med. 2014;174(7):1116–24.
Mungrue K, Honore A. Patterns of Hyperglycemia in Patients With Type 2 Diabetes Mellitus With or Without Hypertension, Requiring Emergency Care and Hospitalization. J Endocrinol Metab. 2015;5(5):291–8.
Cox S, Martin R, Somaia P, Smith K. The development of a data-matching algorithm to define the 'case patient'. Aust Health Rev. 2013;37(1):54–9.
Box GE and Jenkins GM. Time Series Analysis: Forecasting and Control, 2nd ed. San Francisco: Holden-Day; 1976.
Kim YJ, Park S, Yi W, Yu K-S, Kim TH, Oh TJ, Choi J, Cho YM. Seasonal Variation in Hemoglobin A1c in Korean Patients with Type 2 Diabetes Mellitus. J Korean Med Sci. 2014;29(4):550–5.
Liang WW. Seasonal Changes in Preprandial Glucose, A1C, and Blood Pressure in Diabetic Patients. Diabetes Care. 2007;30(10):2501–2.
Tsujimoto T, Yamamoto-Honda R, Kajio H, Kishimoto M, Noto H, Hachiya R, Kimura A, Kakei M, Noda M. Seasonal variations of severe hypoglycemia in patients with type 1 diabetes mellitus, type 2 diabetes mellitus, and non-diabetes mellitus: clinical analysis of 578 hypoglycemia cases. Medicine. 2014;93(23), e148.
Lipska KJ, Warton EM, Huang ES, Moffet HH, Inzucchi SE, Krumholz HM, Karter AJ. HbA1c and Risk of Severe Hypoglycemia in Type 2 Diabetes. Diabetes Care. 2013;36(11):3535–42.
IDF. International Diabetes Federation. IDF Diabetes Atlas, 7th edn. Brussels, Belgium: International Diabetes Federation, 2015. 2015.
Farmer RD, Emami J. Models for forecasting hospital bed requirements in the acute sector. J Epidemiol Community Health. 1990;44(4):307–12.
MV is supported by APA scholarship from Monash University, BdC is supported by a National Heart Foundation Future Leader Fellowship. SZ is supported by a NHMRC Senior Research Fellowship.
This work received no specific source of grant or funding.
Availability of data and materials
Dataset to be made available upon request.
MV: study design, literature review and synthesis, data collection, data analysis, critical discussion, drafting and revision of the manuscript, AE: study concept and design, data analysis, critical discussion and revision of the manuscript, NN: critical revision of the manuscript, KS: data collection, critical revision of the manuscript, BdC: study conception and design, interpretation of the data, critical revision of the manuscript, SZ: study conception and design, critical revision of the manuscript. The authors MV and AE had full access to the data and take responsibility for the integrity of the data and accuracy of the analysis. All authors read and approved the final manuscript.
S Zoungas reports past participation in advisory boards and/or receiving honoraria from Amgen Australia Pty Ltd, AstraZeneca Pty Ltd/Bristol-Myers Squibb Australia Pty Ltd, Janssen-Cilag Pty Ltd, Merck Sharp & Dohme (Australia) Pty Ltd, Novartis Australia Pty Ltd, Sanofi Pty Ltd, Servier Laboratories Pty Ltd and Takeda Australia Pty Ltd for work unrelated to this paper as well as Monash University undertaking contract work for AstraZeneca Pty Ltd/Bristol-Myers Squibb Australia Pty Ltd. The remaining authors declare that they have no competing interests.
Ethics approval and consent to participate
The Monash Health Human Research Ethics Committee approved this study. Reference number: 12197Q.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.