Introduction
Intellectual disability (ID) is a life-long condition that impacts individuals in key areas of health, wellbeing, education, employment, citizenship, community participation and economic sustainability (Thompson et al., Reference Thompson, Bradley, Buntinx, Schalock, Shogren, Snell, Wehmeyer, Borthwick-Duffy, Coulter and Craig2009; McKenzie et al., Reference Mckenzie, Milton, Smith and Ouellette-Kuntz2016; Scior et al., Reference Scior, Hamid, Hastings, Werner, Belton, Laniyan, Patel, Groce and Kett2016; Reynolds et al., Reference Reynolds, Caldwell, Boonchaisri, Ragon and Palmer2022). The American Association on Intellectual and Developmental Disabilities, Diagnostic Statistical Manual – 5th edition (DSM-5; 2013) and International Classification of Diseases – 10th edition (ICD-10) define ID as limitations in intellectual functioning and adaptive behaviour (World Health Organization, 1978; Schalock et al., Reference Schalock, Borthwick-Duffy, Bradley, Buntinx, Coulter, Craig, Gomez, Lachapelle, Luckasson, Reeve, Shogren, Snell, Spreat, Tasse, Thompson, Verdugo-Alonso, Wehmeyer and Yeager2010; American Psychiatric Association, 2013). Intellectual functioning refers to an individual's general mental capacity, such as learning, reasoning, problem-solving, etc. Adaptive behaviour is defined in terms of functioning in three areas of conceptual (e.g. language, numeracy and self-direction), social (e.g. social problem solving, interpersonal skills and social responsibility) and practical skills (e.g. personal care, vocational skills, community participation and healthcare) that are learned and performed by people in their everyday lives.
The DSM-5 and ICD-10 identify ID among individuals based on an intellectual quotient of less than 70 on cognitive test batteries and a below-average (or lower) adaptive functioning quotient reported by adaptive behaviour measures. Although these international organisations provide similar diagnostic frameworks for the identification of ID, changes in diagnostic practices (e.g. the tightening and broadening of the diagnostic criteria for ID) have been observed over the past few decades. These changes could have had an impact on prevalence estimates for ID over the years (Mckenzie et al., Reference Mckenzie, Milton, Smith and Ouellette-Kuntz2016). Despite changes in diagnostic practices, the identification of ID continues to inform service planning and allocation of resources through structural supports, interventions and accommodations. International goals for such efforts are outlined by the United Nations' sustainable development goals (SDGs) 2030 (United Nations-Department of Economic and Social Affairs Disability, 2018). SDGs 1, 2, 3, 4, 8, 10, 11 and 17 are directly applicable to individuals with ID (United Nations-Department of Economic and Social Affairs, 2019). These SDGs promote inclusive and equitable education and healthcare. The SDGs also promote employment, poverty alleviation, hunger eradication and provide social opportunities (United Nations-Department of Economic and Social Affairs Disability, 2018).
Several intervenable prenatal, perinatal and postnatal risk factors that interfere with normal brain development can result in ID among newborns (Huang et al., Reference Huang, Zhu, Qu and Mu2016; Nemerimana et al., Reference Nemerimana, Chege and Odhiambo2018). Prenatal factors during pregnancy such as advanced parental age, exposure to alcohol and tobacco, chronic maternal illnesses (e.g. diabetes, hypertension, asthma, etc.), certain maternal infections (e.g. rubella, chickenpox, syphilis, etc.) and maternal malnutrition increases the risk for ID (Nemerimana et al., Reference Nemerimana, Chege and Odhiambo2018). Perinatal factors such as birth complications, low birth weight, preterm birth and exposure to perinatal infections are other risk factors for the development of ID. Postnatal risk factors include lead and mercury poisoning, severe child malnutrition, central nervous system malignancies, etc. Additionally, inequalities in prenatal, perinatal and postnatal child developmental outcomes pertaining to ID for low- and middle-income countries (as per their gross national income) have been reported (Walker et al., Reference Walker, Wachs, Grantham-Mcgregor, Black, Nelson, Huffman, Baker-Henningham, Chang, Hamadani and Lozoff2011; de Graaf et al., Reference De Graaf, Steegers and Bonsel2013). Several of the risk factors that contribute to ID could be prevented by providing quality healthcare services during pregnancy, delivery and soon after birth to mothers and children (Sutherland et al., Reference Sutherland, Yelland and Brown2012). In addition, healthcare services are known to vary according to the individual-level per-capita income and country-level socio-economic position as measured using highest attained education (Stolk et al., Reference Stolk, Kulik, Le Rutte, Jacobson, Richardus, de Vlas and Houweling2016; Nusselder et al., Reference Nusselder, Rubio Valverde, Bopp, Brønnum-Hansen, Deboosere, Kalediene, Kovács, Leinsalu, Martikainen and Menvielle2021). As these healthcare needs may not be adequately met in several countries, primary studies on the distribution of ID and their meta-analyses indicate an increasing trend of prevalence of ID in low- and middle-income countries compared to high-income nations (Carulla et al., Reference Carulla, Reed, Vaez-Azizi, Cooper, Leal, Bertelli, Adnams, Cooray, Deb and Dirani2011; Maulik et al., Reference Maulik, Mascarenhas, Mathers, Dua and Saxena2011).
Studies published in the past decade on the worldwide prevalence of ID either have reported a descriptive range of rates estimated at 1–3% (Maulik et al., Reference Maulik, Mascarenhas, Mathers, Dua and Saxena2011; McKenzie et al., Reference Mckenzie, Milton, Smith and Ouellette-Kuntz2016) or have focused on a specific country (Arora et al., Reference Arora, Nair, Gulati, Deshmukh, Mohapatra, Mishra, Patel, Pandey, Das, Divan, Murthy, Sharma, Sapra, Aneja, Juneja, Reddy, Suman, Mukherjee, Dasgupta, Tudu, Das, Bhutani, Durkin, Pinto-Martin, Silberberg, Sagar, Ahmed, Babu, Bavdekar, Chandra, Chaudhuri, Dada, Dass, Gourie-Devi, Remadevi, Gupta, Handa, Kalra, Karande, Konanki, Kulkarni, Kumar, Maria, Masoodi, Mehta, Mohanty, Nair, Natarajan, Niswade, Prasad, Rai, Russell, Saxena, Sharma, Singh, Singh, Sumaraj, Suresh, Thakar, Parthasarathy, Vyas, Panigrahi, Saroch, Shukla, Rao, Silveira, Singh and Vajaratkar2018), or certain age groups (Olusanya et al., Reference Olusanya, Davis, Wertlieb, Boo, Nair, Halpern, Kuper, Breinbauer, De Vries, Gladstone and Halfon2018, Reference Olusanya, Wright, Nair, Boo, Halpern, Kuper, Abubakar, Almasri, Arabloo, Arora, Backhaus, Berman, Breinbauer, Carr, De Vries, Del Castillo-Hegyi, Eftekhari, Gladstone, Hoekstra, Kancherla, Mulaudzi, Kakooza-Mwesige, Ogbo, Olsen, Olusanya, Pandey, Samms-Vaughan, Servili, Shaheen, Smythe, Wertlieb, Williams, Newton, Davis and Kassebaum2020). ID has been associated with high mortality rates among the older populations, referred to as a mortality disadvantage (Landes, Reference Landes2017; Reppermund and Walker, Reference Reppermund and Walker2021). Similarly, several risk factors associated with ID, such as lack of optimal nutrition and inadequate access to child and maternal health services, are significant challenges faced by countries with lower socio-economic development (Graham, Reference Graham2005; Simkiss et al., Reference Simkiss, Blackburn, Mukoro, Read and Spencer2011; Banks et al., Reference Banks, Kuper and Polack2017). Therefore, there is a need to delineate detailed prevalence rates among individuals with ID across various age groups and countries of varying socio-economic development to inform national and global policies. These policies could target services to quantify better the disease burden and the resources required to address the needs and rights of individuals with ID as mandated by the SDGs. To this end, this paper aims to report worldwide trends in prevalence estimates and the annual percentage change (APC) for ID based on the Global Burden of Diseases (GBD) 2019 study (Vos et al., Reference Vos, Lim, Abbafati, Abbas, Abbasi, Abbasifard, Abbasi-Kangevari, Abbastabar, Abd-Allah and Abdelalim2020), across sex, age, socio-economic development (i.e. low, low-middle, middle, high-middle and high) and for all regions listed in the GBD 2019 study across 30 years (1990–2019).
Methods
Data source
Individuals with ID and related conditions (i.e. ICD-10 codes F70–79 or ICD-9 codes 317–319) were included in this study (World Health Organization, 1993). While the GBD 2019 study (Vos et al., Reference Vos, Lim, Abbafati, Abbas, Abbasi, Abbasifard, Abbasi-Kangevari, Abbastabar, Abd-Allah and Abdelalim2020) provides details on the methodology, some relevant details about the data sources and data processing for ID are presented here. Those with ID were added to the dataset based on ICD-10 and legacy ICD-9 codes (World Health Organization, 1978, 1993). The prevalence rates for ID were calculated from IQ survey data, cognitive function data, international educational attainment data, relevant registries and other survey data (Vos et al., Reference Vos, Lim, Abbafati, Abbas, Abbasi, Abbasifard, Abbasi-Kangevari, Abbastabar, Abd-Allah and Abdelalim2020). They combine data from standardised tests for IQ from various sources using analytical methods described earlier. Data were then modelled using Disease Modelling – Meta-regression (Dismod-MR 2.1) and distribution fitted using mean and prevalence. Uncertainty intervals (UI) were calculated based on 1000 draws of simulations. This expressed the relevant uncertainty from random probability distributions of input sources, variations in sex, age, diagnostic uncertainties, data manipulations and choice of models used in these estimations.
Socio-demographic information
This study sourced data for annual prevalence rates and estimated the number of individuals with ID for specific age groups, sexes and socio-demographic index (SDI) from 1990 to 2019 based on the GBD 2019 study (Vos et al., Reference Vos, Lim, Abbafati, Abbas, Abbasi, Abbasifard, Abbasi-Kangevari, Abbastabar, Abd-Allah and Abdelalim2020). Separate information was extracted for both females and males for their global prevalence rates and numbers. Global prevalence rates and numbers were separately extracted for the five age groups that included 0–9, 10–24, 25–49, 50–69 and 70+ years. Regions listed in the GBD 2019 study were grouped according to their SDI levels in 2019. SDI is a composite made of several variables that predict a region's social and economic conditions (Global Burden of Diseases Collaborative Network, 2020). These socio-economic conditions play a key role in predicting the health of the relevant regions. Thus, the values for low, low-middle, middle, high-middle and high SDI regions were extracted. All regions were also separately extracted for a more detailed perspective of their prevalence rates. The results across each region are attached in the online Supplementary Appendix.
Data analyses
This study complies with the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) statement (Stevens et al., Reference Stevens, Alkema, Black, Boerma, Collins, Ezzati, Grove, Hogan, Hogan and Horton2016). Global prevalence rates of ID are presented in the tables as represented in a per 100 000 of the population (as a standard practice of reporting prevalence in epidemiological studies) and in text in the Results section as a percentage or per 100 of the population (by dividing the prevalence rate per 100 000 by a 1000) for ease of reading. Similarly, the number of individuals is presented in the tables as n × 1000 and in millions (by dividing the number by 1 000 000) in the text within the Results section for each category. Trends in prevalence rates over the three decades (i.e. one time point for each year) of this study were assessed by calculating the APC for the overall global population and each age range within the SDI regions separately. The prevalence rate, age and SDI were extracted for the APC calculations. Following previously used methodology, models were fitted for each row where the prevalence rates were regressed against the year (Kim et al., Reference Kim, Fay, Feuer and Midthune2000). The (eβ − 1) × 100 indicates the APC, and the (eβ ±1.96s.e. − 1) × 100 indicates the 95% confidence intervals (CIs) of the APC, where β and s.e. denote the coefficient and standard error. APC was considered statistically significant, where the 95% CI did not include a zero. A positive APC indicates an increasing trend, whereas a negative APC indicates a decreasing trend. Consistent with recommended practices for inequality estimation (Mackenbach et al., Reference Mackenbach, Martikainen, Menvielle and De Gelder2016), absolute and relative inequalities were calculated by comparing the average prevalence rates of the higher SDI regions (high and high-middle SDI) with lower SDI regions (low and low-middle SDI). Absolute inequalities were calculated as a difference of the average prevalence rates for each year of the higher SDI regions from the lower SDI regions. Similarly, the relative inequalities were calculated as a ratio of the average prevalence rates for each year of the higher SDI v. lower SDI regions. Further trends in absolute and relative inequalities in prevalence rates were mapped across the years. Global heat maps for prevalence rates and APCs were also created. All analyses were carried out using R-studio (version 1.2.1335).
Results
ID prevalence by sex, age and SDI
Globally, 92.8 million individuals or approximately 1.74% of the population were estimated to have ID in 1990 (Table 1). There was a steady increasing trend in the global number of individuals through 2019, with an estimated 107.6 million individuals with ID. But the proportion of individuals with ID continued to decrease to 1.39% of the global population in 2019, with an APC of −0.80 (−0.88 to −0.72). There were slightly higher estimates for the number and proportion of males v. females with ID. The 2019 data indicate 54.9 million males (1.42%) v. 52.7 million females with ID (1.37%).
UI, uncertainty intervals.
a Prevalence rates are per 100 000.
A considerable reduction was observed in the number (and prevalence rate) of individuals with ID across the increasing age categories. In 1990 only 0.98 million individuals (0.49% of the population) had ID among those aged 70+ years compared to 29 million (2.38% of the population) in the 0–9 years age group. Similarly, in 2019 there were only 1.92 million individuals with ID in the 70+ years age group (0.41% of the population) when compared to 25.8 million (1.98% of the population) in the 0–9 years age group.
In 1990, there was a relatively higher number and prevalence rate (38.2 million and 3.38%) of individuals with ID in the low-middle SDI regions compared to the low SDI (13.8 million and 2.62%, respectively), middle SDI (27.43 million and 1.60%), high-middle SDI (9.57 million and 0.83%) and high SDI (3.78 million and 0.46%) regions. This was true for 2000, 2010 and 2019, where the highest number and proportion of individuals with ID were in the low-middle SDI regions. The global distribution of the prevalence rate of ID (Fig. 1) shows the regional skews with regions in Asia and Africa skewing towards higher prevalence rates. Detailed country-specific values for prevalence rate in 2019 and APC for 30 years are given in the online Supplementary Appendix.
ID prevalence in each SDI stratified by age
Table 2 shows the number of individuals with ID in the five SDI regions separated by the five age groups in 1990, 2000, 2010 and 2019. In the low SDI regions, the 0–9 years age group had the highest number in 2010 (8.67 million), while in all other age groups, the highest number was observed in 2019. In the low-middle SDI regions, the highest number aged 0–9 years were found in 1990 and 2000 with 13.04 million individuals in both those time points. The highest number in the 10–24 years age group was found in 2010 with 15.16 million individuals. For the three older age groups, the highest number of individuals were all found in 2019.
The prevalence rates for all SDI regions by the five age groups (Fig. 2) show an overall decreasing trend over three decades. For all the age groups, the low-middle SDI and low SDI regions had the highest prevalence rates while the other regions were ordered with the next highest being the middle SDI regions, followed by the high-middle SDI region and the high SDI regions. Low SDI regions had an increasing trend until 2004 for the 0–69 year age groups, while the 70+ year age group had an overall increasing trend.
Among the low SDI regions, the APC values were negative for individuals with ID aged 0–9 years, 10–24 years and 25–49 years (Table 2). The APC value was close to 0 for individuals aged 50–69 years, with a CI including 0. In addition, for individuals aged 70 years and above, the APC value of 0.97 indicates a significantly increasing trend. For the other four SDI regions (low-middle, middle, high-middle and high SDI) over the 30 years, the APC values were all negative, and their CIs did not include the value of 0 for all age groups, indicating a decreasing trend of prevalence rates for ID in each age group.
The absolute inequalities estimate between low and high SDI regions showed marked differences according to age groups, with differences higher among younger age groups than older age groups across the three decades (Fig. 3). The 50–69 and 70+ years age groups had increasing absolute inequalities, while the other age groups had a decreasing trend. On the relative scale, while inequalities between low and high SDI regions increased for all age-groups, there was marked increase from three-fold to over five-fold among the 70+ years age group over the three decades.
Discussion
This study quantified the global trends in the prevalence of individuals with ID. In 2019, an estimated 107.62 million individuals with ID formed about 2% of the global population. While there was a global trend of increasing number of individuals with ID, the prevalence rate continued to decrease across the three decades in this study. Consequently, absolute inequalities reduced or remained largely stagnant over the three decades. However, relative inequalities increased suggesting the need for more targeted services for people with ID in low SDI regions particularly in older age groups. The most notable findings were that regions with less socio-economic development had a far greater prevalence rate of ID, and the prevalence rate of ID reduced with age.
There was a slightly higher number and prevalence rate for males v. females, and it agrees with earlier studies that assessed the prevalence rates of individuals with ID (Boyle et al., Reference Boyle, Boulet, Schieve, Cohen, Blumberg, Yeargin-Allsopp, Visser and Kogan2011; McKenzie et al., Reference Mckenzie, Milton, Smith and Ouellette-Kuntz2016). The global data for the number of individuals with ID separated into age groups show that the younger ages had passed peak numbers in 2000 or 2010. Among the 25+ years age group, the numbers continued to rise in the three decades through 2019. This suggests a demographic shift in the number of individuals with ID, where the population overall gets older due to lower number of children globally being diagnosed with ID over time. This aligns with global population trends, where the number of individuals with ID follows the overall declining fertility rates (Pantazis and Clark, Reference Pantazis and Clark2018; Da Silva et al., Reference Da Silva, De Oliveira, Dionizio, De Santana, Bahadori, Dias, Ribeiro, Gomes, Ferreira, Ferreira, De Moraes, Silva, Barnabe, De Araujo, Santana and Monteiro2021). The more alarming trend in the age groups is the precipitous decrease in prevalence rates with increasing age. This suggests a significantly higher mortality over the lifespan among those with ID v. others (Glover et al., Reference Glover, Williams, Heslop, Oyinlola and Grey2017; O'Leary et al., Reference O'Leary, Cooper and Hughes-Mccormack2018; Cooper et al., Reference Cooper, Allan, Greenlaw, Mcskimming, Jasilek, Henderson, Mccowan, Kinnear and Melville2020; Das-Munshi et al., Reference Das-Munshi, Chang, Bakolis, Broadbent, Dregan, Hotopf, Morgan and Stewart2021; Doyle et al., Reference Doyle, O'Sullivan, Craig and Mcconkey2021). There was almost a five-fold reduction in the 70+ year old category v. the 0–9 year old category. An earlier systematic review found higher mortality rates among individuals with ID (Landes, Reference Landes2017; O'Leary et al., Reference O'Leary, Cooper and Hughes-Mccormack2018; Reppermund and Walker, Reference Reppermund and Walker2021). These included studies were restricted to high SDI countries due to the absence of studies in other regions.
The analyses by SDI categories found a significant difference in the prevalence rates of ID which, was 5–6 times higher in the low- and low-middle income regions v. the high SDI regions. This increase was seen across time, where the more developed regions had progressively lower prevalence rates. Comparisons across time points indicate that the low SDI regions had the highest overall number of individuals with ID in 2019. Other SDI regions had their highest numbers in progressively earlier time periods, with the high SDI group peaking at or earlier than 1990. This large difference in rates could be attributable to the increased exposure to risk factors that include a lack of maternal and child health care, malnutrition, more births with proportionally higher genetic illnesses, etc. (Graham, Reference Graham2005; Maulik et al., Reference Maulik, Mascarenhas, Mathers, Dua and Saxena2011; Simkiss et al., Reference Simkiss, Blackburn, Mukoro, Read and Spencer2011; Banks et al., Reference Banks, Kuper and Polack2017). Thus, there are structural inequalities based on the socio-economic development of regions that predispose a greater number of individuals to ID in the lower SDI regions. The trends of overall increases in relative inequalities show the widening gap between the high and low SDI regions, and the relatively unchanged or increasing absolute inequalities emphasise the need for greater actions needed for care of those with ID in the lower SDI regions. These findings necessitate greater levels of resource-intensive care to assist with activities of daily living for individuals with ID in regions with a higher prevalence of ID (Global Burden of Diseases Collaborative, 2017).
The study found a global trend of decreasing prevalence rates of ID for all age groups except for the 50+ year groups in the low SDI regions. It is important to note that while all regions had similar APC reductions for the younger age groups, the trend across the age ranges suggests a pattern of greater reductions for the older age groups. This trend could be due to systemic changes in longevity that happened earlier in the higher SDI regions, followed later by the lower SDI regions (The Global Burden of Disease Child and Adolescent Health Collaboration, Reference Reiner, Olsen, Ikeda, Echko, Ballestreros, Manguerra, Martopullo, Millear, Shields, Smith, Strub, Abebe, Abebe, Adhena, Adhikari, Akibu, Al-Raddadi, Alvis-Guzman, Antonio, Aremu, Asgedom, Asseffa, Avila-Burgos, Barac, Bärnighausen, Bassat, Bensenor, Bhutta, Bijani, Bililign, Cahuana-Hurtado, Malta, Chang, Charlson, Dharmaratne, Doku, Edessa, El-Khatib, Erskine, Ferrari, Fullman, Gupta, Hassen, Hay, Ilesanmi, Jacobsen, Kahsay, Kasaeian, Kassa, Kebede, Khader, Khan, Khan, Khang, Khubchandani, Kinfu, Kochhar, Kokubo, Koyanagi, Defo, Lal, Kumsa, Larson, Leung, Mamun, Mehata, Melku, Mendoza, Mezgebe, Miller, Moges, Mohammed, Mokdad, Monasta, Neupane, Nguyen, Ningrum, Nirayo, Nong, Ogbo, Olagunju, Olusanya, Olusanya, Patton, Pereira, Pourmalek, Qorbani, Rafay, Rai, Ram, Ranabhat, Renzaho, Rezai, Ronfani, Roth, Safiri, Sartorius, Scott, Shackelford, Sliwa, Sreeramareddy, Sufiyan, Terkawi, Topor-Madry, Tran, Ukwaja, Uthman, Vollset, Weldegwergs, Werdecker, Whiteford, Wijeratne, Yonemoto, Yotebieng, Zuhlke, Kyu, Naghavi, Vos, Murray and Kassebaum2019).
Limitations
The wide UI ranges across all regions suggest a scarcity of worldwide data on the prevalence of ID (Eurostat, 2015; Friedman et al., Reference Friedman, Parrish and Fox2018). Although this study includes data from three decades, the direct comparison of the 70+ year group v. the youngest should be taken with caution as they are not the same cohort at different times. But the trend of decreasing prevalence of ID in all age groups suggests that this may be a worst-case scenario, as the current younger cohort possibly has a lower prevalence rate than the corresponding 70+ year group at their corresponding younger age.
Implications
The worldwide higher mortality among those with ID calls for global action to help achieve SDGs specifically on providing equitable health and wellbeing for all by 2030. This finding suggests the need for considerable initiatives for preparing global healthcare systems to cater to the needs of those with ID. The inequalities estimated within this study highlight services for people with ID need to be scaled according to regional disadvantage, otherwise inequalities in social and health outcomes will escalate over time. This is especially true since many of these deaths are due to preventable causes which include respiratory and cardiovascular conditions (O'Leary et al., Reference O'Leary, Cooper and Hughes-Mccormack2018).
The progressively greater prevalence of ID among lower SDI regions suggests a higher propensity for these regions to have risk factors associated with ID. As risk factors range from prenatal, perinatal and postnatal factors, it is important to concentrate efforts to provide maternal support to help avert these risk factors for ID. This is especially important in the lower SDI regions, especially South Asia, where the prevalence of ID is very high. Thus, there is a need for quantifying the causal variables to calculate the attributable risk of ID to produce policies and interventions that reduce the cases of preventable ID. This could contribute towards a greater equity between various regions in the world.
Conclusions
In conclusion, our investigation provides detailed estimates of the prevalence of ID across different regions over three decades. The study found a decreasing trend of prevalence rates of ID, despite an increased number of individuals diagnosed with ID primarily observed in the low and low-middle SDI regions. Remarkable differences in ID prevalence rates and trends found between higher and lower SDI areas show significant inequities. Similarly, significant decreases in the prevalence rates among the older age groups suggest the lack of comprehensive and tailored healthcare provision both within and outside of clinical settings.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S2045796022000701.
Data
All data used in this paper came from openly available dataset hosted by Global Burden of Diseases collaborative. The code for analyses and detailed extraction from the data set are given in detail in the online Supplementary Appendices.
Financial support
This study was supported by Singapore's Ministry of Education (MOE-DEV 05/19 CM). The funding source played no role in the acquisition, analyses, interpretation or publication of this manuscript. All decisions for this work were based on the decisions of the authors alone.
Conflict of interest
The authors affirm that there are no conflicts of interest.
Ethical standards
Approval with an ethics board was not required for this anonymised open dataset.