Article Text
PDF
Abstract
Background While there is a long history of measuring death and disability from injuries, modern research methods must account for the wide spectrum of disability that can occur in an injury, and must provide estimates with sufficient demographic, geographical and temporal detail to be useful for policy makers. The Global Burden of Disease (GBD) 2017 study used methods to provide highly detailed estimates of global injury burden that meet these criteria.
Methods In this study, we report and discuss the methods used in GBD 2017 for injury morbidity and mortality burden estimation. In summary, these methods included estimating cause-specific mortality for every cause of injury, and then estimating incidence for every cause of injury. Non-fatal disability for each cause is then calculated based on the probabilities of suffering from different types of bodily injury experienced.
Results GBD 2017 produced morbidity and mortality estimates for 38 causes of injury. Estimates were produced in terms of incidence, prevalence, years lived with disability, cause-specific mortality, years of life lost and disability-adjusted life-years for a 28-year period for 22 age groups, 195 countries and both sexes.
Conclusions GBD 2017 demonstrated a complex and sophisticated series of analytical steps using the largest known database of morbidity and mortality data on injuries. GBD 2017 results should be used to help inform injury prevention policy making and resource allocation. We also identify important avenues for improving injury burden estimation in the future.
- populations/contexts
- methodology
- descriptive epidemiology
- statistical issues
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.
Statistics from Altmetric.com
INTRODUCTION
The Global Burden of Disease (GBD) study is a comprehensive assessment of population health loss. GBD has expanded in scope since its original release in 1994 (GBD 1990) and was most recently updated in autumn 2018 (GBD 2017).1–7 Each update of the study has provided updated results through the most recent year of data availability as well as increasingly refined detail in terms of locations, age groups and causes. In addition, GBD incorporates new data as well as updated methods for each annual release that represent the expanding complexity of the study. Cumulatively, the increasing volume of data and increasingly sophisticated estimation methods have necessitated near-continual refinements in terms of data processing, statistical modelling, computational storage and processing as well as global collaboration with the over 4000 GBD collaborators in over 140 countries and territories.
Historically, injuries have formed one of the three broad cause groups in the GBD cause hierarchy alongside the other two main groups of health loss (communicable, maternal, neonatal and nutritional diseases; non-communicable diseases). Not surprisingly, there is considerable variation in how morbidity and mortality are estimated across different causes in the GBD hierarchy and study design. The methods for estimating morbidity and mortality from injuries have evolved over time through the most recent release of GBD 2017. Historically, there have been certain challenges in injuries burden estimation, some of which have been addressed and updated over time, and some of which remain as methodological challenges to address as population health measurement develops more sophisticated modelling strategies. For example, methodological challenges that have been identified over the past three decades in population health research have included obtaining data in data-sparse, burden-heavy areas of the world, developing adjustments for ill-defined causes of death, separately estimating cause of injury from the bodily harm that results from an injury event and adjusting for known biases in data, such as underestimation in sexual violence data.3 8 9 Cumulatively, the global injuries research community has developed a wide array of methodological innovations and advancements to overcome many of these challenges, although undoubtedly the science will continue to advance as higher-quality datasets become available, as modelling methods improve and as computational processing power becomes more accessible to population health research groups around the world.
Many studies have been published based on different releases of the GBD study, ranging from studies on intentional injuries in the eastern Mediterranean to detailed assessments of traumatic brain injury and spinal cord injury disability rates on a global scale.10 11 While this array of published GBD injury studies demonstrates a broad spectrum of expert knowledge on specific injuries or specific geographies or both, it is also critical to recognise that population health is a rapidly evolving, collaborative science that has benefited from near-continual improvements even through the current updates being implemented for GBD 2019. As a result, it should benefit the scientific enterprise to focus on publishing the most updated results with perspective on global, demographic and temporal patterns, and on sharing iterative updates on the current state of the science of GBD injuries burden estimation. The goal of this study is to comprehensively review and report methods used for GBD 2017 and associated publications that have gone through extensive collaborator-review and peer-review processes.
METHODS
GBD 2017 study
GBD is predicated on the principle that every case of death and disability in the population should be systematically identified and accounted for in the formulation of global disease and injury burden. On the side of mortality, every death that occurs in the population should have one underlying cause of death which can be assigned to a cause in a mutually exclusive, collectively exhaustive hierarchy of diseases and injuries that can cause death. These data can be used in a method described below to calculate cause-specific mortality rates and years of life lost. For morbidity, every non-fatal case of disease or injury should have an amount of disability assigned for some period of time. These data can be used in a process described below to estimate the incidence, prevalence and years lived with disability. Summing morbidity and mortality from some cause form the burden from that cause, expressed as disability-adjusted life-years (DALY). For causes with known risk factors, some portion of this burden may be explained by exposure to that risk factor. Across causes within some population, it is also a principle of GBD that the sum of all cause-specific deaths should equal all-cause mortality in the population, and that rates of incidence, prevalence, remission and cause-specific mortality can be reconciled with one another such that all death and disability in a population is internally consistent across causes and geographies. As examples, the sum of different types of road injury cases must sum up to overall road injuries, and the sum of deaths from different injuries in a given country must sum up to the estimate of all-injury deaths. The principle of internal consistency extends to populations used in GBD, where every birth, death and net migration must be accounted for in the population estimates which form the denominators of GBD results. While there is immense complexity in the process summarised above, it is important to begin with these core principles which govern the computation processes at the heart of GBD burden estimation. A summarised overview of key GBD 2017 methods is also provided in online supplementary appendix 1.
Supplemental material
GBD study design and hierarchies
GBD study design, including cause-specific methods, is described in a high level of detail in associated publications.2–7 In addition to the injury-focused methods described in this paper, it is important to define hierarchies used in the GBD study design. In particular, GBD 2017 was built around a location hierarchy where different subnational locations (eg, US states, India states, China provinces) which form a composite of a national location (eg, the USA, India, China). National locations are aggregated to form GBD regions, which are then aggregated to form GBD super regions. These designations affect the modelling structure and utilisation of location random effects, processes which are described in more detail later. The country-level and regional-level GBD location hierarchy used in GBD 2017 is provided in online supplementary appendix table 1. In addition to locations, GBD processes are conducted to produce estimates for every one of 22 age groups, male and female sex and across 28 years from 1990 to 2017 (inclusive). Age-standardised, all-age and combined sex results are also computed for each GBD result. Exceptions exist to the rules above, for example, self-harm is not permitted to occur in the 0–6 days (early neonatal) age group in the GBD age hierarchy. There are no sex restrictions placed on any GBD injury causes, although these restrictions exist for other GBD causes, such as cancers like prostate, cervical and uterine being related to one sex.
Supplemental material
GBD injury classification
In the GBD cause hierarchy, injuries are part of the first level of the GBD cause hierarchy, which consists of three broad groups: communicable, maternal, neonatal and nutritional diseases; non-communicable diseases and injuries. Additional levels of the GBD cause hierarchy provide additional detail. The hierarchy of injuries in GBD is provided in table 1. The organisation of the hierarchy has implications both in terms of how results are produced and in terms of analytical and processing steps which are discussed in more detail below. Case definitions including International Classification of Diseases (ICD) codes used to identify injury deaths and cases are provided in table 2.
Global Burden of Disease cause-of-injury hierarchy
Case definitions for cause of injury in GBD 2017
GBD separates the concept of cause of injury from nature of injury. Cause of injury (eg, road injuries, falls, drowning) have historically been used for assigning cause of death as opposed to the ‘nature’ of injury, which more directly specifies the pathology that resulted in death. For example, an individual who falls, fractures his or her hip, undergoes surgery and then develops hospital-acquired pneumonia and dies while hospitalised would still have a fall as the underlying cause of death, regardless of whether sepsis or some other disease process leads to death more proximally in the chain of events. In this individual, the ‘nature’ of injury would have been specified as a hip fracture, since it is the bodily injury that would dictate the disability this person experiences. Since it is evident that a hip fracture is more disabling than a mild skin abrasion, it is important for measuring non-fatal burden to consider both the cause and the nature in the formulation of complete injury burden. A full list of nature of injury is provided in table 3.
GBD nature of injury
Cause-specific mortality and years of life lost
As described above, cause-specific mortality is measured for every cause of injury in the GBD cause hierarchy with the exception of foreign body in the ear and sexual violence, which undergo only non-fatal burden estimation (described in more detail below). GBD adheres to five general principles for measuring cause-specific mortality, which are described in more detail elsewhere but are summarised as follows.12 First, GBD 2017 identifies all available data. For injuries, this includes vital registration (VR), vital registration samples, verbal autopsy (VA), police records and mortuary/hospital data. VR is the preferred data source but is not available in every location in the GBD location hierarchy. Prior VA research has demonstrated that VA is more accurate for certain injury causes than it is for certain diseases.13 Police data undergo additional validity checks to ensure that systematic under-reporting does not occur in comparison to VR data, which is described in more detail in a related publication.6 The second general principle relevant to injury mortality estimation is maximising comparability and quality of the dataset. For the purposes of injury mortality estimation, this process is largely focused on (1) ensuring appropriate accounting for different ICD code versions used for cause of death data classification over time, (2) redistribution of ill-defined causes of death (described in more detail elsewhere) and (3) processing VA studies into usable data that map to the GBD cause hierarchy.8 9 12 The third general principle for injury cause of death models in GBD 2017 is to develop a diverse set of plausible models. This process is conducted via the Cause of Death Ensemble model (CODEm) framework, which is the standard, peer-reviewed cause of death estimation process used extensively in the GBD study. CODEm generates a large set of possible models based on covariates suggested by the modeller based on expert input and literature review (eg, alcohol for road injuries) and then runs every plausible model, which can range into the thousands per cause. These models can be conducted in both rate space and cause fraction space and use an assortment of combinations among the user-selected covariates (table 4). Fourth, the predictive validity of each one of these submodels is tested using test-train holdouts, whereby a specific model is trained on a portion of data and tested on a separate portion to determine out-of-sample predictive validity. Once the submodels are conducted and predictive validity is measured, then an ensemble model is developed out of the submodels. The submodels and the ensemble model are then subject to the fifth principle, which is to choose the best-performing models based on out-of-sample predictive validity. The chosen models may be a single cause model or an ensemble of models. Beyond these processes, which have become automated with expert review in the GBD processing architecture, there is also considerable time required by the analysts, modellers, collaborators and principal investigators who are involved in the GBD study. Such processes also come under expert scrutiny via the GBD Scientific Council and the peer-review process in the annual GBD capstone publications.2–7
Covariates used in GBD cause of death models
Once submodels and ensemble models have been conducted for each cause in the GBD cause hierarchy, a process to correct for cause of death rates to ensure internal consistency is conducted. Specifically, each subcause within some overall cause is rescaled such that, for example, every subtype of road injuries sums to road injuries deaths overall, and then road injuries and other transport injuries sum to equal the overall transport injuries cause. As this cascades to the overall cause hierarchy and the overall all-cause mortality rates, cause-specific mortality across all causes ultimately equals the overall mortality in the population. An example of an injuries cause of death model with vital registration data (Colombia, females) is shown in figure 1. A similar model with relatively less data is shown in figure 2 (Honduras, females). While data are absent in more recent years in Honduras, the model is still able to follow temporal trends, age patterns and broader geographical patterns by harnessing signals from covariate-based fixed effects (eg, alcohol consumption per capita) and location-based random effects (eg, the regional trends in Central Latin America and patterns in neighbouring countries). All cause of death models from GBD 2017 are publicly available for review (https://vizhub.healthdata.org/cod/). Cause-specific deaths are converted to cause-specific mortality rates (CSMRs) using GBD populations. Once CSMRs are established, years of life lost (YLLs) are computed as the product of CSMRs and residual life expectancy at the age of death. The residual life expectancy is based on the lowest observed mortality rate for each age across all populations over 5 million. For example, if a death from road injuries occurs at age 25 and the residual life expectancy is 60 years, then there are 60 YLLs attributed to that death. If the death had occurred at age 50 with a residual life expectancy of 38 years, then 38 YLLs would be attributed. Life tables used for GBD 2017 are provided in related publications. 7
Cause of Death Ensemble model with data points for road injuries in Colombia for females.
Cause of Death Ensemble model with data points for road injuries in Honduras for females
Injury incidence, prevalence and years lived with disability
After cause-specific models for each cause of injury in the GBD cause hierarchy are conducted, the non-fatal estimation process is conducted. An overview of this process is depicted in figure 3. In the first stage, we estimate the incidence of injuries warranting medical care using DisMod-MR 2.1 (abbreviated DisMod). DisMod is a meta-regression tool for epidemiological estimation that uses a compartmental model structure whereby a healthy population may become diseased or injured, at which point the individual either remains a prevalent case, goes into remission or dies. DisMod essentially fits differential equations to reconcile the transitions between these different compartments, so that the final posterior estimate for each epidemiological parameter can be explained in the context of the other parameters. Similar to the principles described in CODEm, DisMod uses all available data, ranging from incidence data to cause-specific mortality rates from the corrected CODEm results, to produce estimates for every age, sex, year and location. For the purposes of injuries, we established our case definition for non-fatal injuries as injuries that require medical care. This is a necessary case definition as we do not want to consider minor stumbles and falls, for example, that led to no actual bodily harm as injuries for GBD, since they would not have any associated disability. These models are conducted only for injury causes as opposed to the nature of injuries references above. Each data input is designated based on type of data—specifically, inpatient data, outpatient data, surveillance data, survey data and literature studies that are population-representative. We model incidence rates for hospital admissions for injuries, so the non-inpatient data sources get adjusted according to their classification so that the model inputs are consistent as injuries that warranted or received inpatient medical care. The coefficients measured by DisMod that were used for adjustment are provided in table 5. Input data for injury cause incidence models included sources identified as part of systematic reviews conducted in past GBD cycles, new sources identified by the GBD collaborator network and new sources of clinical data and other injuries data obtained by the core injuries burden estimation team at the Institute for Health Metrics and Evaluation at the University of Washington. In addition, CSMRs from the corrected CODEm models described above are used in this stage of DisMod modelling. The list of non-fatal injury sources used in GBD 2017 is provided in online supplementary appendix table 2. The completed DisMod models for inpatient incidence for each cause of injury are publicly available at https://vizhub.healthdata.org/epi/.
Supplemental material
Injuries non-fatal estimation flow chart.
Covariates and coefficients used in Global Burden of Disease incidence cause models
Once an incidence cause model is constructed for each cause of injury, an extensive analytical ‘pipeline’ follows which converts injury cause incidence into years lived with disability. First, inpatient incidence is split into inpatient and outpatient incidence using coefficients empirically measured by DisMod. The outpatient coefficients for each injury cause are also included in table 5. Separate pipelines are then conducted for inpatient and outpatient injury incidence—each step below can be considered to have been run for both streams of data, for each cause of injury. After the coefficient is applied, incidence is adjusted by the excess mortality rate measured by DisMod to essentially remove injury cases that died after the injury occurred. Once these deaths are removed from the incidence pool, the resulting steps are applied to these surviving cases of injury. First, each new case of injury is considered to have 47 possible ‘natures’ of injury that can result. These are the types of bodily injury that are considered to be possible outcomes from a given injury cause. The proportion of new cases of injury that would have some nature of injury as the most disabling outcome is determined based on dual-coded clinical data sources where both the cause and nature of injury were included as ICD codes.10 Of note, one limitation of this process is that due to computational demands, it is currently only possible to apportion the most disabling nature of injury for each new case of injury. As such, the probability that each nature of injury is the most disabling nature of injury for some cause of injury is modelled in a Dirichlet regression such that the probabilities sum to 1. In other words, each nature of injury has some probability of being the most disabling injury suffered by the victim of some cause of injury, but if multiple natures of injury occurred, then the less disabling injuries are not captured as part of that injury cause’s disability. This limitation has been recognised as a limitation of GBD injury burden estimation in various peer-reviewed articles and will likely be addressed in future GBD updates as computational efficiency improves.3 10
The probability distributions of each cause-nature are computed separately for each age, sex, year and location. At this point, the analytical stage has the age-specific, sex-specific, year-specific, location-specific incidence of a cause-nature combination, for example, the incidence of road injuries that led to a cervical-level spinal cord injury in males aged 20–24 years in 2017 in Stockholm, Sweden. The next step converts these incidence estimates into short-term and long-term injury incidence estimates, where long-term disability is defined as having a lower functional status 1 year postinjury than at the time of injury. These probabilities were measured using long-term follow-up studies.14–20 For some natures of injury, such as lower extremity amputation, the probability of being a long-term injury is 1. The probabilities of short-term versus long-term injury for each cause-nature combination are used to split the incidence values into short-term and long-term pipelines. The long-term incidence is then converted to prevalence using the ordinary differential equation solver used in DisMod, which also uses as an input excess mortality estimated for certain natures of injury such as traumatic brain injury and spinal cord injury conducted in a previous systematic review and meta-analysis. The short-term incidence is converted to prevalence by multiplying incidence and duration of injury, where duration of injury was either computed directly from follow-up studies or, in the case of unavailable data, estimated by an expert clinical panel involved in previous iterations of the GBD study. Since access to medical treatment is assumed to affect duration of injury and disability, the GBD Healthcare Access and Quality Index is used to estimate the proportion with and without access to medical treatment on a location-specific basis.21 The average duration for short-term injury is therefore calculated as the percentage treated multiplied by treated duration added to the percentage untreated multiplied by the untreated duration. The output from this step is the short-term prevalence of each cause-nature combination. Short-term prevalence is subtracted from long-term prevalence at this stage to avoid double counting the same case of injury. Once short-term and long-term prevalence estimates for each cause-nature are computed, then disability weights as derived by the Salomon et al process are assigned to each injury nature.22 Short-term disability weights by injury nature are shown in table 6, which does not include amputations since we assume they cause only long-term disability. The full list of long-term disability weights by injury nature, location and year are provided in online supplementary appendix table 3, which does not include foreign body in respiratory system, foreign body in gastrointestinal and urogenital system, foreign body in ear and superficial injury of any part of body, since we assume these natures of injury do not cause long-term disability. After disability weights are assigned to each injury case, years lived with disability for each cause of injury are calculated as the prevalence of each health state multiplied by the corresponding disability weight and then summed across natures of injury for each cause to compute years lived with disability (YLDs) for each age, sex, year and location for that injury cause. YLDs then undergo comorbidity adjustment used across the GBD study whereby comorbid cases of disease and injury in the population are simulated and adjusted disability weights are computed. These processes are described in more detail in GBD literature.3 GBD 2017 provided an important methodological update whereby nature of injury results, regardless of cause of injury, could be reviewed in the results from this process; this has enabled more advanced GBD research such as measuring the burden of traumatic brain injury and spinal cord injury, measuring the burden of facial fractures and measuring the burden of hand and finger fractures.10
Supplemental material
Short-term disability weights for each nature of injury
Sexual violence
Sexual violence follows a different analytical pathway than the other causes of injury. This process is shown in figure 4. We used the same study framework as was developed for other injury rates in the GBD 2017 study to estimate the yearly proportion of the population that experienced at least one episode of sexual violence in the past year, using a case definition of any sexual assault including penetrative sexual violence (rape) and non-penetrative sexual violence (other forms of unwanted sexual touching). To inform the sexual violence estimates, we identified data in 93 countries that met the case definition above. This resulted in 263 site-years of data, which mainly were derived from surveys such as Demographic and Health Surveys and Reproductive Health Surveys. Similar to our other injury models, we used DisMod 2.1 to model prevalence. The sexual violence prevalence model used study-level covariates for each type of survey question, for example, we used a study-level covariate to identify surveys that identify penetrative sexual violence only to account for how the overall incidence of sexual violence is greater than this value. This model also used a covariate on alcohol use in litres per capita for each location to help fit the model in data-sparse locations. Once yearly prevalence was measured, sexual violence cases undergo a process by which short-term disability from the physical and psychological harm of sexual violence cases is assigned to each prevalent case; however, long-term sequelae of sexual violence are currently not captured in this process, which has been a known limitation of sexual violence estimation in the GBD framework.
Sexual violence estimation flow chart. HAQI, Healthcare Access and Quality Index.
Disability-adjusted life-years
After estimation of cause-specific mortality and YLLs as well as non-fatal health outcomes estimation including YLDs, DALYs are calculated as the sum of YLLs and YLDs for each cause of injury. YLDs are also calculated for each nature of injury category.
GATHER statement
GBD 2017 adheres to the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER). GATHER is described in more detail in online supplementary appendix 2.
Supplemental material
RESULTS
Results for all GBD 2017 injury estimates are available in associated publications as well as online. Specifically, results by age, sex, year, subnational location and nature of injury can be viewed and downloaded online via the GBD Results Tool (http://ghdx.healthdata.org/gbd-results-tool) and GBD Compare (https://vizhub.healthdata.org/gbd-compare/). These results are available in terms of incidence, prevalence, YLDs, cause-specific mortality, YLLs and DALYs, expressed in counts, rates, and percentages. Analytical code and input datasets are available at http://ghdx.healthdata.org.
CODEm models
Model performance metrics for each injury cause model in GBD 2017 are provided in table 7. Model performance metrics for CODEm models include root mean square error (RMSE) for in-sample tests and out-of-sample tests, percentage of data points that correctly predict the trend in-sample and out-of-sample and percentage of data points that are present within the 95% uncertainty intervals (UIs) of the model fit. RMSE in-sample is generally better than RMSE out-of-sample, which is an expected result that also demonstrates the importance of performing out-of-sample predictive validity tests. While the correct trend is predicted in approximately one in five models, this may also be related to more dynamic temporal trends in injury mortality patterns over time. In general, most data points exist within the 95% UI of the model fit (mean: 98.5% in-sample, 98.0% out-of-sample).
Performance metrics for each cause-of-injury CODEm model
Incidence models
Model performance metrics for each injury cause model in GBD 2017 are provided in table 8. These model performance metrics include in-sample coverage and RMSE of estimated results for cause-specific mortality, excess mortality and incidence. There are no performance metrics for CSMR or excess mortality for foreign body in eyes since we do not estimate mortality from this cause of injury. For incidence, the in-sample coverage average was 55.3% across cause-of-injury models and ranged from a low of 26% in falls to a high of 88% in poisoning by carbon monoxide. Incidence RMSE ranged from a low of 1.04 in pedestrian road injuries to a high of 4.86 in foreign body in eye.
Performance metrics for each cause-of-injury DisMod model
DISCUSSION
Many considerable advancements have been made in the measurement of global injury burden since early versions of the GBD Study. Novel datasets, sophisticated statistical modelling and global collaboration have all facilitated the advancement of injury burden measurement science. Many more advancements in future updates should be possible as larger datasets become available and as computational power allows for more detailed measurement processes. Continued global collaboration will be an integral component. Suggested priority items for the advancement of injury burden estimation are as follows:
First, while much of the global injury burden occurs in low-income and middle-income countries, these countries are frequently the most data-sparse. GBD has rigorously attempted to collect all available data, including police records and verbal autopsy studies and inpatient and outpatient records; however, it is likely that additional data sources in data-sparse countries exist. Parties who are aware of additional data sources that could be used in the GBD estimation framework should consider joining the GBD collaborator network to contribute new sources of data to be used in future estimation updates.
Second, computational and data limitations make it difficult to account for the full disability that might be experienced in the setting of multiple injuries. For example, if an individual sustains a below-neck spinal injury and an upper extremity amputation, the amputation is not directly accounted for in the prevalence or YLD estimate of the injury cause to which this disability is attributed. This problem quickly grows in complexity, as one can imagine an event like a road injury leading to multiple contusions and abrasions, several fractures in different anatomical sites, a mild traumatic brain injury and a spinal cord injury. There are over 3.6 million permutations of injury if one considers only 10 possible natures of injury, making it difficult to quantitatively measure these relationships by cause of injury and by age, sex, year and location. Future research to address this limitation may focus on simulation studies that model the probability of different comorbid injury combinations to better inform disability weight applications.
Third, more data could be used for nature of injury measurement. Traumatic brain injury and spinal cord injury registries, for example, are not currently directly compatible with the GBD injury estimation framework yet provide rich epidemiological information. Future updates to GBD should focus more attention on incorporating data that measure burden of nature of injury in terms of incidence, prevalence or excess mortality. Incorporating these types of data would require a method to be developed such that estimates were internally consistent across cause-nature distributions. While the methods and data required for this update would be complex, they would represent a large increase in the available data that could be used for GBD injuries estimation.
Fourth, measuring the total burden of sexual violence has proven to be a challenging area of estimation in the GBD framework. As noted in the ‘Methods’ section of this paper, one known limitation is how long-term sequelae and conditions may not be adequately accounted for in sexual violence burden estimation. In order to attribute burden from major depressive disorder, anxiety disorders, self-harm and substance use disorders, measuring the relative risk of developing these conditions for victims of sexual violence would allow for population attributable fractions to be calculated and DALYs from these conditions to be attributed to sexual violence. While the premise of this methodological update is relatively simple, currently there are relatively few studies to inform these relative risks, and conducting and adding such studies in the future would be recognised as a major achievement in GBD research as it would allow for more accurate estimation of lifetime disability caused by sexual violence. This effort would moreover represent an important contribution to research surrounding the Sustainable Development Goals related to sexual violence and women’s rights. 23 24
Fifth, non-fatal injuries from conflict and natural disaster are challenging to estimate because of data sparsity in areas that are afflicted by these events. Fatalities are estimated after such events, but there is still considerable injury burden among the population that survives. Since data collection systems and hospitals may also be destroyed in these events, it becomes difficult to collect adequate non-fatal injury data. Global collaboration should also focus on identifying sources of data on non-fatal and fatal injury cases in conflict and natural disaster events.
It will be important to monitor the effects of implementing these priorities as injury measurement science continues to evolve. Global collaborations including the GBD enterprise should monitor performance statistics and utilisation of results by research groups and ministries to track how improvements to injury measurement progress over time. Scientific dialogue and collaboration must be a major focus, and the GBD enterprise is a good forum to support this kind of data sharing. For example, a collaborative effort between researchers in Vietnam and the Institute for Health Metrics on Evaluation on developing a study on Vietnam injury burden following GBD 2017 led to identifying the use of the Vietnam National Injury Survey, which was then added for estimation in GBD 2019. Increasing data collection standardisation efforts should be emphasised as a priority in all countries, particularly countries where data coverage on injuries is sparse. Ongoing dialogue via scientific publications and international conferences should also continue to serve as a forum to discuss data and methodological updates that can continue to refine the science of injuries estimation in GBD.
CONCLUSION
Measuring injuries burden in GBD is a complex scientific endeavour that leverages large amounts of data, a complex analytical framework and a global research network. GBD 2017 included more comprehensive detail of injury burden than any other known efforts to date. GBD 2019 and future updates will continue to add detail and refine methods in the interest of providing injury burden estimates that are robust, accurate and timely. Expanded injury data collection efforts will be a critical component of future injury burden estimation.
What is already known on the subject
Global Burden of Disease (GBD) 2017 provided an extensive peer-reviewed assessment of death and disability.
GBD 2017 methods have been reviewed and updated iteratively as new methods and data become available.
Measuring injury burden in GBD 2017 is complex due to differences in measuring cause of injury versus nature of injury and the temporal difference between them.
What this study adds
This capstone study details key estimation methods that are used for measuring the global burden of injuries as described in related publications in this journal.
More detailed methods descriptions and model performance metrics from GBD 2017 are provided in this study than in related studies.
This study also includes suggested future directions for improving injury burden research.
Acknowledgments
Seyed Aljunid acknowledges the Department of Health Policy and Management, Faculty of Public Health, Kuwait University and International Centre for Casemix and Clinical Coding, Faculty of Medicine, National University of Malaysia for the approval and support to participate in this research project. Alaa Badawi acknowledges support by the Public Health Agency of Canada. Till Bärnighausen acknowledges support by the Alexander von Humboldt Foundation through the Alexander von Humboldt Professor award, funded by the German Federal Ministry of Education and Research. Traolach Brugha received support from NatCen Social Research (http://natcen.ac.uk/) via NHS Digital and Department of Health and Social Care London, for the Adult Psychiatric Morbidity Survey (APMS) programme. Felix Carvalho received support from UID/MULTI/04378/2019 with funding from FCT/MCTES through national funds. Vera M Costa acknowledges support from grant (SFRH/BHD/110001/2015), received by Portuguese national funds through Fundação para a Ciência e Tecnologia (FCT), IP, under the Norma Transitória DL57/2016/CP1334/CT0006. Kebede Deribe is supported by a grant from the Wellcome Trust [grant number 201900] as part of his International Intermediate Fellowship. Tim Driscoll acknowledges that work on occupational risk factors was partially supported by funds from the World Health Organization. Eduarda Fernandes acknowledges support from UID/QUI/50006/2019 with funding from FCT/MCTES through national funds. Yuming Guo acknowledges support from Career Development Fellowships of the Australian National Health and Medical Research Council (numbers APP1107107 and APP1163693). Sheikh Mohammed Shariful Islam is funded by a Fellowship from National Heart Foundation of Australia and Institute for Physical Activity and Nutrition, Deakin University. Mihajlo Jakovljevic acknowledges support by the Ministry of Education Science and Technological Development of the Republic of Serbia through the Grant number OI175014; publication of results was not contingent upon Ministry's censorship or approval. Sudha Jayaraman acknowledges support from: NIH R21: 1R21TW010439-01A1 (PI); Rotary Foundation Global Grant #GG1749568 (PI); NIH P20: 1P20CA210284-01A1 (Co-PI) and DOD grant W81XWH-16-2-0040 (Co-I), during the period of this study. Yun Jin Kim acknowledges support from a grant from the Research Management Centre, Xiamen University Malaysia [grant number: XMUMRF/2018-C2/ITCM/0001]. Kewal Krishan is supported by UGC Centre of Advanced Study (CAS II) awarded to the Department of Anthropology, Panjab University, Chandigarh, India. Mansai Kumar acknowledges support from FIC/ NIH K43 1K43MH114320-01. Amanda Mason-Jones acknowledges support by the University of York. Mariam Molokhia is supported by the National Institute for Health Research Biomedical Research Center at Guy’s and St Thomas’s National Health Service Foundation Trust and King’s College London. Ilais Moreno Velasquez is supported by the Sistema Nacional de Investigación (SNI, Senacyt, Panama). Mukhammad David Naimzada acknowledges support from Government of the Russian Federation (Agreement No – 075-02-2019-967). Duduzile Ndwandwe acknowledges support from Cochrane South Africa, South African Medical Research Council. Stanislav S. Otstavnov acknowledges the support from the Government of the Russian Federation (Agreement No – 075-02-2019-967). Ashish Pathak acknowledges support from Indian Council of Medical Research (ICMR), New Delhi, India (Grant number 2013-1253). Michael R Phillips acknowledges support in part by a grant from the National Natural Science Foundation of China (No.81761128031). Abdallah M. Samy received a fellowship from the Egyptian Fulbright Mission Program. Milena Santric Milicevic acknowledges the support from the Ministry of Education, Science and Technological Development, Republic of Serbia (Contract No. 175087). Seyedmojtaba Seyedmousavi was supported by the Intramural Research Program of the National Institutes of Health, Clinical Center, Department of Laboratory Medicine, Bethesda, MD, USA. Rafael Tabarés-Seisdedos was supported in part by the national grant PI17/00719 from ISCIII-FEDER. Sojib Bin Zaman acknowledges support from an "Australian Government Research Training Program (RTP) Scholarship.”
Footnotes
Funding This study was funded by The Bill and Melinda Gates Foundation. SLJ conducts research for a grant on influenza and RSV which is funded in part by Sanofi Pasteur.
Competing interests Dr. James reports grants from Sanofi Pasteur, outside the submitted work. Dr. Driscoll reports grants from World Health Organisation, during the conduct of the study. Dr Shariful Islam is funded by a Fellowship from National Heart Foundation of Australia and Institute for Physical Activity and Nutrition, Deakin University. Dr. Ivers reports grants from National Health and Medical Research Council of Australia, during the conduct of the study. Dr. Jozwiak reports personal fees from TEVA, personal fees from ALAB, personal fees from BOEHRINGER INGELHEIM, personal fees from SYNEXUS, non-financial support from SERVIER, non-financial support from MICROLIFE, non-financial support from MEDICOVER, outside the submitted work. Walter Mendoza is currently Program Analyst Population and Development at the Peru Country Office of the United Nations Population Fund-UNFPA, which does not necessarily endorses this study. Dr. Rakovac reports grants from World Health Organization, during the conduct of the study. Dr. Sheikh reports grants from Health Data Research UK, outside the submitted work. Dr. Singh reports personal fees from Crealta/Horizon, Medisys, Fidia, UBM LLC, Trio health, Medscape, WebMD, Clinical Care options, Clearview healthcare partners, Putnam associates, Spherix, Practice Point communications, the National Institutes of Health and the American College of Rheumatology, personal fees from Speaker’s bureau of Simply Speaking. Dr. Singh owns stock options in Amarin pharmaceuticals and Viking pharmaceuticals. Dr. Singh serves on the steering committee of OMERACT, an international organization that develops measures for clinical trials and receives arms-length funding from 12 pharmaceutical companies. Dr. Singh serves on the FDA Arthritis Advisory Committee. Dr. Singh is a member of the Veterans Affairs Rheumatology Field Advisory Committee. Dr. Singh is the editor and the Director of the UAB Cochrane Musculoskeletal Group Satellite Center on Network Meta-analysis, outside the submitted work. Dr. Stein reports personal fees from Lundbeck and Sun, outside the submitted work.
Patient consent for publication Not required.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Availability of input data varies by source. Select data are available in a public, open-access repository. Select data are available on reasonable request. Select data may be obtained from a third party and are not publicly available. All results from the study are included in the article or uploaded as supplementary information or are available online.
