Patterns of computed tomography utilisation in injury management: latent classes approach using linked administrative data in Western Australia

Purpose Whilst computed tomography (CT) imaging has been a vital component of injury management, its increasing use has raised concern regarding ionising radiation exposure. This study aims to identify latent classes (underlying patterns) of CT use over a 3-year period following the incidence of injury and factors predicting the observed patterns. Method A retrospective observational cohort study was conducted in 21,544 individuals aged 18 + years presenting to emergency departments (ED) of four tertiary public hospitals with new injury in Western Australia. Mixture modelling approach was used to identify latent classes of CT use over a 3-year period post injury. Results Amongst injured people with at least one CT scan, three latent classes of CT use were identified including a: temporarily high CT use (46.4%); consistently high CT use (2.6%); and low CT use class (51.1%). Being 65 + years or older, having 3 + comorbidities, history with 3 + hospitalisations and history of CT use before injury were associated with consistently high use of CT. Injury to the head, neck, thorax or abdomen, being admitted to hospital after the injury and arriving to ED by ambulance were predictors for the temporarily high use class. Living in areas of higher socio-economic disadvantage was a unique factor associated with the low CT use class. Conclusions Instead of assuming a single pattern of CT use for all patients with injury, the advanced latent class modelling approach has provided more nuanced understanding of the underlying patterns of CT use that may be useful for developing targeted interventions. Supplementary Information The online version contains supplementary material available at 10.1007/s00068-023-02303-y.


Introduction
Globally, injury is one of the leading causes of mortality in adults [1].Computed tomography (CT) imaging is a vital component of diagnosis and management for most injuries [2], especially in the emergency department (ED) since CT offers high sensitivity and specificity in the evaluation of trauma patients.A Western Australian study showed that the rate of CT examinations in EDs doubled from 58 to 105 per 1000 ED presentations between 2003 and 2015, with injury being the second most likely reason for CT use [3].
The increasing volume of these examinations has raised concerns for patients [4], and is receiving increased attention from health care providers, regulators and the media [4,5] due to its association with healthcare expenditure and unnecessary exposure to ionising radiation [6].Although CT scans comprise a small proportion of all diagnostic radiological procedures, a report from the USA in 2009 found that CT and nuclear medicine accounted for 36% of the total radiation exposure and 75% of the medical radiation exposure of the population [7].The UK review in 2011 shows that typical radiation effective dose for common CT scans in adults increased between 20% (head) to 400% (high-resolution chest) compared to the dose observed in the review in 2003 [8,9].Risks are accentuated in trauma patients who are potentially more radiosensitive due to their relatively younger age compared with patients presenting for medical conditions [10], notwithstanding evidence showing that the average age of trauma patients is increasing [11].On Extended author information available on the last page of the article average, trauma patients received a mean effective dose of 22.7 mSv which is almost ten times the average annual background radiation dose (2.4 mSv) [10].Therefore, monitoring long-term CT use is important for understanding the appropriateness of current imaging practice.Assessment of change in use of CT scanning incorporating characterisation of patients with different patterns of use could help policy-makers construct a comprehensive strategy to evaluate population risk.The finding may also influence clinicians employing CT scanning to consider the individual risks of excessive ionizing radiation to their patients.
Current studies generally focus on examining the use of CT during a single life event, such as an ED presentation [3,[12][13][14][15] or hospital admission [16,17], over short periods.Patients who had substantial increase in CT use following the life event may not be captured in the short periods of observation.Availability of large linked administrative data at the individual level has provided opportunity to use advanced analytic methods such as trajectory modelling techniques to identify unobserved patterns (latent classes) within the population over time [18,19].Patterns are defined directly from the data rather than using arbitrarily pre-specified thresholds [19].In this study, we used mixture modelling of linked administrative data in Western Australia (WA) to empirically identify discrete underlying patterns of CT use over a 3-year period following an initial injury requiring assessment in an ED and factors predicting these patterns.

Methods
This was a retrospective observational cohort study reported following the Reporting of studies conducted using Observational Routinely-Collected health data (RECORD) statement [20].

Data source
The study used de-identified individual level health administrative data linked by the WA Data Linkage System and the Australian Institute for Health and Welfare pertaining to adults aged 18 years and over who had any hospitalisation (except for pregnancy), ED presentation or CT scan in WA between 1 January 2003 and 31 December 2016 [21].For the cohort, the following data were available: (i) WA Public Hospital ED data collection (EDDC) (2003-2016) which provided details on all ED presentations at all public hospital EDs in WA; (ii) WA hospital morbidity data collection (HMDC) data which contained information related to inpatient care in all hospitals (public and private) in WA between 2003 and 2016; (iii) Information on all CT scans performed in WA was sourced from: [1] Picture Archival communication System (PACS) records from 2003 to 2015 which contained records of CT undertaken in all tertiary public hospitals and the majority of secondary public hospitals in WA; and [2] Medicare Benefits Schedule (MBS) claim items from 2005 to 2015 for CT scans for subsidised under Australia's universal health insurance scheme Medicare and undertaken for private patients in private and public hospitals and community-based radiology practices.(iv) The WA Death Registry of Births, Deaths and Marriages records from 2004 to 2015 containing information on all deaths in including the age, date and cause of death.Further details of each datasets are published elsewhere [22].

Study population
This study included patients aged 18 years or older who was: (1) presenting to any of the four tertiary (teaching or major) hospital EDs in WA; and (2) with a first-time (incident) injury in 2012, defined as no history of injury recorded in both ED and HMDS data over a 3-year lookback period prior to their injury record in 2012.
The restriction to tertiary hospital EDs was made to ensure a consistent cohort could be captured and that initial CT scanning in the ED could be comprehensively identified.This was necessary because (i) private hospital EDs are not included in the ED data collection and (ii) the PACS data does not capture all secondary (district or regional) hospitals.Injured patients were identified using the International Statistical Classification of Diseases, tenth revision, Australian modification (ICD-10-AM) code S00-T14 recorded in the diagnosis field of the ED data [23].

Use of CT over three-year following the incident injury
All CT scans from the date of ED presentation for the incident injury to 3 years post injury or dead whichever comes first were included.CT scans recorded on the same day and the same anatomical area were counted as one to avoid over counting.The number of CT scans was summed by quarter over the three-year follow-up period.Amongst individuals who had at least one CT scan over the follow-up period, Box-Cox transformation was applied to produce more normally distributed data for use in modelling latent classes of CT use [24].Individuals with no CT scans over the whole study period were grouped as "no CT use" group and excluded from trajectory modelling for computational efficiency (i.e. to improve convergence and time).

Potential predictors
The following variables were captured: sex; age at incident injury ED presentation (classified into 18-24 years, 25-64 years, 65-79 years and 80 + years) and Indigenous status (used only for adjustment of confounding in models, not reported due to ethics approval conditions).Socioeconomic status, in quintiles, was derived from postcode of residence at time of ED presentation using the socioeconomic index for areas (SEIFA) index of relative socioeconomic disadvantage [25].Patient's residential postcode at incident ED presentation was classified into major cities, inner regional areas, outer regional areas, remote and very remote areas according to the Accessibility and Remoteness Index of Australia (ARIA) [26].External cause of injury was sourced from ED data and was categorised into four main groups after consultation with local ED physicians.The groups were transport/pedestrian, fall, force (blunt force, cut/ pierced or stabbed, shot by weapon and contact with machinery) and others (bite or sting, contact burn, contact with fire or flame, exposure or poisoning by chemicals, electrocution, other cause and unknown) and published elsewhere [22].Anatomical area of injury was classified as head, neck, thorax, abdomen, extremity and multiple injuries.Patients were classified as either admitted to hospital or discharged following the incident ED presentation since those who died in ED were excluded from the cohort.Severity of injury was evaluated using the ICD-10-AM-based Injury Severity Score (ICISS) which has been evaluated internationally and in Australia [23,[27][28][29] and classified into mild (survival > 99%), moderate (survival > 94%) and severe/very severe (survival < 94%).Details of how severity scores were derived have been published previously [15].Date of ED presentation was classified into weekday and weekend/public holiday.Mode of arrival was classified into three groups: private transport, ambulance (including air ambulance) and other.ED presentation time was classified into three-time blocks: day (8:00-15:59), evening (16:00-23:59) and night (0:00-7:59).Length of stay was the number of days between admitted date and discharged date of the hospitalisation associated with the incident ED presentation with those not admitted given the value of 0. History of comorbidity was based on the number of comorbidities captured in HMDC records within 5 years prior to the incident injury, classified as none, 1-2 and 3 + comorbidities using the Multipurpose Australian Comorbidity Scoring system [30] using ICD-AM-10 across all diagnostic fields.History of CT scanning was classified as having at least one CT scan recorded in the year prior to the index injury.Number of ED presentations and hospitalisations related to injury during the follow-up period were also captured to adjust for the potential impact of the new event.

Statistical analysis
Descriptive statistics were used to describe socio-demographic and clinical characteristics of the study cohort at baseline-the time presenting to tertiary ED for the incident injury.We followed the three-step approach [31] including: (1) determining number of classes; (2) assigning the most likely class membership to the study cohort; and (3) a subsequent separate analysis using a multinomial logistic regression model to identify factors predicting different CT use classes.
For identifying latent classes, we utilised finite mixture models (FMMs), originally from Pearson and Henrici [32].Generically, the approach involves endogenously, and probabilistically, splitting the sample into a finite number of discrete latent classes.Within each class, observations are relatively homogenous, but potentially heterogeneous across them.The probabilistic splitting of the sample is usually achieved by employing multinomial logit (MNL) techniques.Within each class, the same statistical model applies, but are characterised by differing parameters of that density.Importantly, in this way, the same covariates can have differing effects on the outcome variable within each class.The optimal number of classes is usually determined by information criteria and/or entropy [33].Posterior probabilities are used to predict class membership of individuals (based on the maximum probability rule).
The defined classes are "labelled" to reflect the typical composition of the class members with respect to the outcome variable ( (i) total number of CT scans within a class to reflect total burden of CT use in each class, (ii) average number of CT scans in an individual within each class to capture burden at individual level, and (iii) time between the first and last CT to indicate level of intensity in exposing to CT scan as suggested by previous publication [34].These indicators were used to support interpretation of the classes.
In addition to describing the classes by simple descriptive statistics, multinomial logistic regression model, a more nuanced approach, was used to the explain the posterior probabilities-which drive class membership-by a range of predictor variables.No CT use group was acted as the reference (i.e.control group) in the model to estimate the relative risk of being classified into other CT use class across a range of the predictors.The analysis was conducted using the lcmm package under R version 4.0 [35] and STATA MP Version 16 [36].

Cohort characteristics
A total of 21,544 individuals aged 18 years and older were identified with an incident injury presenting to a tertiary level ED in 2012.The majority of the study cohort were male (58.2%), aged less than 45 years old (60.7%), and lived in a major city (93.5%) (Table 1).People with a socio-economic status of least disadvantaged accounted for 43.1%.

Latent classes of CT use over three-year follow-up
Based on entropy and information criteria, the model with three classes was preferred, with patterns of CT use presented in Fig. 1 (details model output is presented in Appendix 1).The three classes were characterised descriptively according to the pattern of CT use as follows: class 1 -termed " consistently high use" (2.6% of patients having at least one CT scan over the 3-year follow-up ~ 0.8% of the whole cohort) had a consistently high (but slightly diminishing) use of CT across the study period with no peaks or troughs; class 2-"low use" (51.1% amongst those with at least one CT scan over the 3-year follow-up ~ 16.9% of the cohort) had no CT use at the beginning and a slight increase over the follow-up time; and class 3-"temporarily high use" (46.3% of people with at least one CT scan over the 3-year follow-up ~ 14.6% of the whole cohort) in which the CT use was very high at the start, declining to zero at the end of first year followed by an increase in CT use after 2.5 years.Whilst class 1, "consistently high use", accounted for nearly 1% of the study population, it accounted for 10% of the total CT use in the cohort (Table 2).Notably the average number of CTs per individual was the highest (7.9 CT scans [95%CI 7.29; 8.45]) of all classes.In contrast, class 2, the "low use" and class 3, the "temporarily high use", had an average of 1.7 CT scans [95% CI 1.64; 1.73] and 2.2 CT scans [95% CI 2.16; 2.28], respectively (Table 2).As noted in Table 2, 67% of people presenting to a WA tertiary ED with an incident injury had no CT scanning undertaken within three years of the date of presentation.

Factors predicting classes of CT use
Baseline characteristics of the patients in each latent class are presented in Table 3.There were significant differences in age and sex distribution between no CT use and the CT use classes.The majority of the no CT use group were male (60.3%) and aged under 45 years (70.1%).In contrast, class 1-the consistently high use of CT, was characterised by equal distribution between sexes and a high proportion of people aged 65 + years (63.0%).Both class 2 and class 3 had a similar distribution across the age groups, whilst class 3 had a significantly higher proportion of males.In terms of clinical characteristics, individuals in class 1 and class 3 had a significantly higher proportion of head injuries, hospital admissions for the incident injury, injuries due to falls and moderate or severe injuries than those in either class 2 or the After adjustment for all observed baseline demographic and clinical characteristics, factors significantly associated with all three CT use classes were mild to older age (30 + years), having 3 + comorbidities and history of having CT scan.Number of ED presentations and hospitalisations related to injury during follow-up period were also associated with all 3 CT scanning classes.However, the highest magnitude of association with all these factors was observed in class 1, "consistently high use".(Table 4).
Whilst there was no factor associated with both class 1 and 2 (except for severe/very severe injury) or both class 2 and 3, several factors were associated with both class 1 and 3 including head injury, multiple injury, moderate level of injury and length of stay of the index hospital admission.However, the magnitude of the association was different.Multiple injury and moderate level of injury had the highest association with class 1, whilst head injury had the highest association with class 3 and length of stay was equally associated with both classes 1 and 3.
There was no factor uniquely associated with class 1.In contrast, neck injury, thorax injury, abdomen injury, traffic injuries and injury arrived by ambulance were uniquely associated with class 3. Living in an area of moderate to highest socio-economic disadvantage was only associated with class 2.

Discussion
Our study found a third of patients with injury have at least one CT scan at any time over the 3-year follow-up.Using an advanced classification approach, instead of assuming a single pattern of CT use applied to all injury patients, this study has provided more nuanced understanding of the underlying patterns of CT use that may be useful for developing targeted interventions.Whilst we normally expect that CT scanning would be concentrated during an acute episode of injury (i.e. in ED and hospital), only half of the patients with an incident injury were found to follow that pattern.The remainder (of those who had a CT) had  a substantial number of CT scans performed sometime after the index injury.This study found that nearly 10% of the total CT scans were performed on only 1% of the injured subjects.This group had an average of nearly 8 CT scans done within two and a half years.A recent metaanalysis found that cancer risks increased rapidly for radiation exposures above 55 mSv (which corresponds to the amount patients may get from 3 or more CT scans) [29].
For each additional CT scan, the risk of cancer increased by 0.16 (95% CI 0.13-0.19)[30].Literature data also indicate that multiple CT scans within a short period of time may represent potentially avoidable CT scan procedures [31].Therefore, identification of the consistently high-use group and factors associated with persisting high use in our study can contribute to raising awareness in clinical practice about the potential risk of excessive radiation exposure in this sub-population.Caution is needed when requesting additional CT scans for the high-risk group to minimise the risk of unnecessary radiation exposure.
Although CT scanning has become the screening test of choice for most injuries, its increasing use has been particularly concerning as injury patients tend to be relatively younger, hence, more radiosensitive [39].On average, the cumulative effective dose is 2.6 milli-Sieverts (mSv) per injured patient per ED visit [40].Studies found that the increasing use of CT scans in patients with injury does not correspond with improved patient outcomes [14,15,41].However, findings from previous studies may have overlooked and underestimated the use of CT because of the restriction to a short emergent period (i.e.ED presentation) which failed to account for different patterns/trajectories of CT use over the course of injury management.Our study found only half of injured patients who had CT scanning had a peak of CT use occurring at the time of injury.This will enable the development of targeted interventions aimed at improvements in the efficiency of use of CT scanning for injured patients determining which patients should be the focus of radiation dose reduction strategies.
An interesting finding in our study is that moderate to highest disadvantage SEIFA were only significantly associated with low CT use class.Theoretically socio-economic factor should not affect access to CT scan whilst patients are in public hospitals which are fully funded by the State government.However, for CT scan performed out of hospitals, although they are subsidised under Medicare Benefits Scheme funded by the Federal government, there is an outof-pocket payment for the cost above the Medicare reimbursement level.Since CT scanning is a relative expensive diagnostic imaging modality, the out-of-pocket cost can be a significant barrier for patients with lower socio-economic status.Patients with lower socio-economic status often live in more remote areas (i.e. less accessible areas) where lack of transport may also be a considerable barrier to access CT scanning.This is in line with our previous study which only examined the use of CT scanning in hospitals and suggested that accessibility (i.e.living in rural/remote areas) is a factor driving a reduction of CT scanning use.In this study, both CT scanning use in hospitals and out-of-hospital were captured, we provided further evidence that patients with lower socio-economic status rather than less accessibility (i.e.living in rural/remote areas) is a predictor of low CT use Western Australia after adjusting for other demographic and clinical characteristics.This finding is in line with many previous studies which indicated the association between low socio-economic status and low use of health care services [42,43], such as specialist and the use of CT scanning [44].
The use of the advanced latent class mixture modelling applied in our study has advantages over the conventional cross-sectional approach to better capture high-need-highcost-high-dose patients to improve efficient use and safety of healthcare resources and health outcomes.The latent class mixture modelling approach can further refine the prediction of individuals with a similar pattern, and in our study has successfully identified the top 1% of high CT use with an average of 8 CT scans per person over 3 years post injury.Using the cross-sectional approach obfuscates distinct The high use of CT scanning in the consistent high use class may result in a small but considerable excess cancer risk.However, the risk of missing life-threatening injuries may outweigh the small long-term risk for cancer from imaging tests.A previous study suggested the median effective doses ranged from 2.1 mSv for a head CT to 31 mSv for a multiple abdomen and pelvis CT [45].Another study in Western Australia examining 34 common CT scanning protocols found that the mean effective doses can be as low as 0.4 mSv for CT of the sinuses and as high as 31.2mSv for CT of the whole body angiography [46].About a third of CT protocols had mean effective doses greater than 10 mSv [46].Our results are consistent with a recent systematic review investigating high-use-high-cost patients [47].Our study provides further evidence to warrant intervention to avoid unnecessary CT scans, especially targeting at-risk subpopulations.
Several limitations should be considered when interpreting our findings.Whilst this study used linked whole-of-population health administrative data capturing the use of CT scanning both inside and outside hospital settings, CT scanning used in several secondary hospitals was not fully recorded.Due to incomplete records of CT use in secondary hospitals, this study was limited to injuries presenting to tertiary EDs in Western Australia.This limits generalising our findings to the whole population.Whilst the use of CT following the injury event was captured in all tertiary and most of secondary hospitals as well as private providers, the use of CT may be underestimated as the PACS data did not include all secondary hospitals.In addition, there is no information about the indication for the CT scan.Hence, it is possible that subsequent CT scans may not all be related to injury.This also prevented us from including the indication for the CT examination as a factor in predicting the CT use classes, although other proxy measures of clinical characteristics (i.e.comorbidity and history of hospital admission) were included, in addition to injury group and severity of injury.We cannot access the impact of the new injury on the use of CT scan.However, since the patients were followed up in a short time [3-year period], we assumed the chance of multiple new injuries occurred is minimal.We have included number of ED presentation and hospitalisations related to injury during the follow-up period to account for the potential effect.However, the rich source of longitudinal data has enabled us to fully explore the patterns of CT scanning use post injury and identify subpopulation with high use and risk due to the exposure to potentially large radiation doses.In addition, our study has broadly examined various factors associated with being classified in different classes of CT use that is informative for future interventions/ policies targeting reduction of the unnecessary medical radiation exposure.

Conclusion
Our broad exploration of the patterns of CT use 3-year postincident injury will provide valuable information to assist with interpretation of findings in the current literature as * Low cell count is not allowed to present SEIFA Socio-economic Index for areas, ARIA Accessibility and Remoteness Index of Australia, PH Public holiday, IQR the interquartile range  well as to design future studies investigating the use of CT scans.It was found that the patients with the consistently high use of CT scan accounted for only 2.6% in the CT use groups but consumed substantial proportion of CT scans (10%).This study has also illustrated the use of a latent class mixture model in identifying distinctive and interpretable patterns of the CT use in patients with an incidence of injury.This method facilitates the identification of associated factors (within each class) which can be useful to design future policies/intervention targeting the high CT use sub-population to minimise the risk of exposure to medical radiation.

Fig. 1
Fig. 1 Predicted class trajectory amongst injury patients with at least one CT scan over 3 years

Table 1
Cohort characteristics at the baseline 2012 no CT use group.Class 1 had the highest proportion of people with 3 + comorbidities (70.0%) compared with classes 2 and 3 (38.6%and 35.1%, respectively).

Table 1
(continued) SEIFA Socio-economic Index for areas, ARIA Accessibility and Remoteness Index of Australia, PH Public holiday, IQR the interquartile range $ Force includes blunt force, cut/pierced or stabbed, shot by weapon and contact with machinery $$ Others/unknown includes bite or sting, contact burn, contact with fire or flame, exposure or poisoning by chemicals, electrocution, other cause and unknown

Table 2
Characteristics of CT use classes identified from the finite mixture models amongst people with at least one CT scan over the 3-year follow-up period

Table 3
Characteristics of patients in No CT use and in each CT use class