Summary
Introduction Measles outbreaks increased worldwide during the period of 2017-2019. Similarly, Vietnam experienced one of the largest measles outbreaks in recent decades, with various paediatric patients presenting severe complications. In this study, we aim to identify factors associated with death and severity among children with measles admitted to Vietnam National Children’s Hospital (VNCH) in Hanoi, Vietnam between 2017 and 2019.
Method This single-center retrospective cohort study included 2, 072 patients with measles admitted to VNCH from 1/1/2017 to 31/12/2019. Data on epidemiological, clinical characteristics, vaccine status, and outcomes were collected and summarised. We conducted both univariable and multivariable logistic regression analyses to examine the correlations between various characteristics of hospitalized children and mortality.
Findings In total, there were 2, 072 patients, including 1, 297 (62.6%) males and 775 (37.4%) females. The median age was 9 months (interquartile range 7 − 17). 87.3% of cases had not received any measles-containing vaccine (MCV). 30 (1.4%) patients died, with 40% aged less than 9 months. Only 3 among 30 (10%) who died had received at least 1 dose of MCV. Bronchopneumonia was the most common complication, occurring in 1, 413 (68.2%) patients. The following characteristics were significantly associated with mortality in the multivariable analysis: age under 9 months and age from 9 months to 5 years, residing 20 to 200 kilometres from VNCH, and having co-infection with adenovirus or other hospital acquired infections. Age group was also significantly associated with severity in the multivariable analysis.
Interpretation Vietnam continues to face the threat of future measles epidemics, given the burden of hospitalization and the high rate of complications observed in hospitalized patients. This highlights the critical need to maintain high measles vaccine coverage, particularly by targeting the unvaccinated population. To prevent future outbreaks and lower measles incidence, routine immunization needs to be strengthened, and earlier scheduling of MCV1 needs to be further evaluated. The comprehensive analysis of the 2017-19 measles outbreak presented in this study will contribute to informed decision-making regarding appropriate measures to counteract future resurgences of measles in Vietnam.
Funding No specific grant from funding agencies in the public, commercial or not-for-profit sectors supported the submission and publication of this manuscript.
1. Introduction
Measles is considered one of the most contagious diseases for humans with the basic reproduction number of 12 − 18, which means that each person with measles would, on average, infects 12 − 18 other people in a totally susceptible population [1]. This disease is caused by a paramyxovirus virus, manifesting as a febrile rash illness. It is an airborne disease, which can easily spread through susceptible populations, and the only effective measure for preventing this disease is vaccination. Before the invention and widespread use of the measles vaccine in 1963, frequent outbreaks of the disease occurred approximately every 2–3 years, resulting in approximately 2.6 million deaths annually [2]. It was estimated that measles infected over 90% of children before they reached 15, causing over two million deaths and 15,000 to 60,000 cases of blindness worldwide each year [3, 4]. Accelerated immunization activities have significantly reduced the number of measles deaths. From 2000 to 2018, measles vaccination is estimated to have prevented 23.2 million deaths, resulting in a 73% decrease in global measles deaths from 536, 000 in 2000 to 132, 490 in 2016 [3, 5, 6]. Measles outbreaks happen when people who are not immune to the measles virus are infected and spread the disease to unvaccinated or under-vaccinated populations. To control measles and avoid outbreaks and deaths, at both national and sub-national levels, vaccination coverage rates must reach and be maintained at 95% for both the first and second doses. However, global vaccination coverage rates for first-dose measles-containing vaccine (MCV1) have stagnated at 84-85% for over a decade, and while second-dose measles-containing vaccine (MCV2) coverage rates have been gradually increasing, they are only at 71%. Measles vaccination rates remain significantly below the necessary 95% threshold to control measles and prevent outbreaks [5].
From 2017, measles experienced a surge worldwide, resulting in the highest number of reported cases in 23 years in 2019 [8]. World Health Organization (WHO) and the United States Centers for Disease Control and Prevention (CDC) published a report indicating that the number of reported measles cases increased to 869,770 in 2019, the highest number recorded since 1996, with all WHO regions reporting increases. The number of reporting countries with annual measles incidence of less than 5 cases per 1 million population decreased from 125 in 2016 to 85 in 2019 and an estimated 207,500 people worldwide died from measles, a nearly 50% increase since 2016 [5]. Measles cases surged drastically in WHO’s African Region (AFR) by 1, 606%, in the Region of the Americas (AMR) by 19, 739%, and in the European Region (EUR) by 2, 282% since 2016. 9 countries (Central African Republic, Democratic Republic of the Congo, Georgia, Kazakhstan, Madagascar, North Macedonia, Samoa, Tonga, and Ukraine), accounting for 73% of all reported cases worldwide, experienced large outbreaks with reported measles incidence exceeding 500 per 1 million population [6].
Similar to other countries, between 2017 and 2019, Vietnam experienced one of the largest measles outbreaks in recent decades. In 2018, Vietnam had 8,444 cases of typhus-suspected measles with 1,177 confirmed measles cases, twice as many cases as in 2017. Most measles cases were unvaccinated or incompletely vaccinated children, or had no known vaccination history due to lost immunization records. To prevent the outbreak, the Ministry of Health of Vietnam launched an additional measlesrubella vaccination campaign in late 2018 for 4.2 million children aged between 1 and 5 in vulnerable areas in 57 cities and provinces [9]. In 2019, Vietnam reported over 39,000 suspected cases of typhus linked to measles, with more than 7,100 confirmed cases and four deaths [9, 10]. Analysis of the total infections in Hanoi until April 2019 revealed that up to 14% of children under 9 months had measles. This is the group of children who have not yet reached the age of vaccination and the mother has not been vaccinated against measles, so they do not receive any antibodies during pregnancy. Children between 9 months and 5 years old had the highest incidence of the disease, accounting for 42% of total cases. Furthermore, the Center for Preventive Medicine has reported that approximately 15% of children of vaccination age are not received the required two full doses of the measles vaccine [11]. The objective of this study is to analyze epidemiological and clinical information of over 2,000 children with measles admitted to Vietnam National Children’s Hospital (VNCH), the national referral hospital for children in Hanoi, the capital city of Vietnam, between 1 January 2017 and 31 December 2019. We aim to provide an overview of the outbreak in the northern region of Vietnam and assess the potential implications for Vietnam and neighbouring countries.
2. Method
2.1. Study design and participants
A single-center retrospective cohort study including all children admitted with measles admitted to VNCH between 01/01/2017 and 31/12/2019 was undertaken. During this time, individuals with suspected measles were admitted to VNCH, in addition to other hospitals in the region. Admitted children included either self-referrals (walk-in patients) or referrals from other health facilities. Patients with multiple missing data were excluded.
2.2. Ethics statement
The study was approved by the Ethics Committee for Biomedical Research Vietnam National Children’s Hospital (Reference number: 473/BVNTW-HDDD, IRB identifier: IRB-VN01037/IRB00011976/FWA00028418). As the study involved a retrospective analysis of routinely collected data under strict anonymity, patient or caregiver consent was waived.
2.3. Procedures
Following WHO-recommended surveillance standard, the case definition of measles infection was based on laboratory-confirmed diagnosis. Eligible study participants are those who were diagnosed with measles by IgM or PCR test at any time point during their hospital stay, and vaccine-associated illness has been ruled out. All individuals from VNCH fulfilling the clinical case definition acquired information using the measles suspected case inquiry form. The information in the questionnaire was collected by interviewing parents about vaccination status in children’s vaccination records and using paper or electronic medical records to collect other information. Duplicates and unidentified patients are detected and removed by comparing the patient IDs. The date of admission, measles vaccination status according to the vaccination card (if available), diagnosis upon admission, and location of exposure according to the caregiver report, as well as socio-demographic, clinical features, and outcome, were all included. There was no information provided about the rash’s features, dietary status, or whether the case was likely a main or secondary home case.
Furthermore, to facilitate the analysis of factors associated with the risk of severe diseases, we defined severe patients as those requiring any form of oxygen supplementation and respiratory support, namely oxygen supplementation via nasal cannula, continuous positive airway pressure (CPAP), conventional mechanical ventilation, high-frequency oscillatory (HFO) ventilation, or extracorporeal membrane oxygenation (ECMO). To build an anonymised dataset, selected data were extracted and personal information such as name, date of birth, and permanent address were removed.
2.4. Statistical analysis
Univariable and multivariable Firth’s Bias-Reduced logistic regression analysis was used to assess associations between specific characteristics of hospitalized children and clinical outcomes (survived or died). Comparisons were made between three separate locations of exposure, and the Chi-squared or Fisher’s exact test was employed for the comparison of dichotomous and t-test for continuous variables. Odds ratios (ORs) were computed with a 95% confidence interval (CI), and a p-value less than or equal to 0.05 was considered statistically significant.
The selection of variables for the multivariable analysis was based on several factors. Basic variables such as year of hospitalization, vaccination status (patient having at least one dose of MCV), age, and exposure location were included as they are fundamental factors to consider. Additionally, the distance from the hospital was incorporated as it is expected to have an influence on the time taken to reach the hospital and the duration between symptom onset and testing/admission. Patients residing far from VNCH may initially receive treatment at lower-level hospitals, leading to potential delays in obtaining the correct treatment and adequate level of care. In addition to the previously mentioned factors, the duration of stay in the hospital was included as a variable in the analysis. This choice was motivated by the understanding that the duration of hospitalization can serve as an indicator of the severity of measles. Prolonged hospital stays may suggest a more severe manifestation of the disease. Co-infections and underlying conditions of several bodily systems were also added after ensuring that their correlation with other model covariates would not confound the link between the other variables and death. In the model analysing the risk towards severe diseases, the duration between symptom onset and testing was included as a variable based on its statistical significance. Conversely, in the mortality model, the variables “RSV” (respiratory syncytial virus) and “pertussis” (within the co-infection group of variables) were excluded due to their lack of statistical significance. We report the median and interquartile range (IQR) for age groups and the frequency and percentage for categorical variables. Despite the huge number of univariable tests, we did not explicitly apply a multiple-testing correction since the variable selection technique for the multivariable model was only partially based on p-values. The Pandas software library was used to arrange the data, and the R package “logistf” was utilized for the main analysis [12, 13]. For mapping and distance computation from patient’s place of residence to the hospital, the Azure Maps Route and Geopy libraries were used [14, 15].
All conducted statistical analysis can be found within the following GitHub repository.
2.5. Role of the funding source
No specific grant from funding agencies in the public, commercial or not-for-profit sectors supported the publication of this publication. The corresponding author had full access to all the study’s data and was ultimately responsible for the decision to submit it for publication.
3. Results
Within three years from 1/1/17 to 31/12/19, VNCH admitted 2,072 children who met the clinical case definition for measles.
Figure 2 and 3 depict the residence of hospitalized children throughout the study period, the number of measles cases admitted to VNCH each week, and the number of deceased and critical children. The measles outbreak in Vietnam surged in 2018-19, mirroring the trend of VNCH admissions, which increased from low levels in 2017 to over 50 cases per 2 weeks (25 cases per week) during the months of April and May 2019, before declining to pre-epidemic levels by September 2019. Figure 5 illustrates number of measles cases admitted to VNCH and total measles cases in Vietnam, 2017-2019. From 2017 to 2019, VNCH admitted 132, 621, and 1,319 measles cases, accounting for 58%, 27% and 9% of total measles cases in Vietnam.
The demographic and clinical features of the 2,072 children hospitalized with measles are described in Table 1. In 2017, there were 132 admissions and it quintuple in 2018 with 621 measles cases; more than half of the cases are in 2019 with 1,319 (63.7%) cases. 1,297 (62.6%) hospitalized children were male and 775 (37.4%) were female. The median age was 9.6 months, with an interquartile range (IQR) of 7 to 16.9 months; 43.5% were younger than 9 months, 47.5% were aged 9 months to 5 years, and 9% were older than 5 years. Most children (87.3%) were reported not to have received any MCV and only 9.8% received at least one dose. Measles cases were transmitted from the community in 61.8% of cases, from VNCH in 26.4% of cases, and from other hospitals in 11.6% of cases. 69.3% of patients live in the Red River Delta, including Hanoi, and 84.6% live within 200 kilometers of the hospital. In 66% of all instances, the period between the development of the rash and hospital admission was shorter than one week. Furthermore, 611 (29.5%) cases were admitted to the hospital before the onset date due to a different disease. More than half of patients (50.3%) are hospitalized for less than 1 week and 86% for less than three weeks. Most children (76.5%) have no underlying diseases, with only 4.5% having chronic diseases related to gastrointestinal system and 3.5% having chronic diseases of the neurological system. 849 (41%) patients turned severe after hospitialization, with 685 (33.1%) requiring oxygen supplementation via nasal cannula, and 144 patients requiring conventional mechanical ventilation. Complications were reported in most children, the most prevalent of which was bronchopneumonia (68.2%). Middle ear infection and febrile seizure complications were reported in 30 (1.4%) and 10 (0.5%) children, respectively.
Overall, 30 (1.4%) children died and 2,041 (98.6%) were discharged.
Table 2 shows the demographic and clinical characteristics of hospitalized patients with measles by place of exposure, including community, with VNCH or from another health care facilities. These characteristics were significantly associated with exposure location: admission year, age groups, vaccination status, region of residence, clinical outcome, underlying conditions, maximal form of respiratory used, diagnosis upon admission, different types of coinfections, and clinical progression.
Table 3 depicts the relationships between socio-demographic, clinical characteristics and severity (turning severe). In the univariate analysis, the following factors were significantly associated with the likelihood of becoming severe: admission in 2018 and 2019, age 9 months to 5 years and more than 5 years, vaccination at least 1 dose, hospital transmission, living between 20 and 500 kilometres near VNCH, living in the Northeast, Red River Delta (except Hanoi), or Central Vietnam, mortality, being admitted to the hospital before the onset, hospitalized for more than a week, being tested more than 48 hours after the onset, having underlying diseases on the cardiovascular system, gastrointestinal system, kidney and urology system, neurological system. Having other underlying conditions or no underlying diseases is also heavily linked with severity. Moreover, being diagnosed with measles, or bronchopneumonia, having gastroenteritis, middle ear infection, laryngitis, or febrile seizures, having respiratory syncytial virus (RSV), adenovirus, pertussis, or other co-infections are connected to the risk of severity. Children with the following risk factors had significantly higher/lower odds of severity in the multivariable analysis: admission in 2018 (adjusted OR [AOR] 12.92, 95% CI 7.05 − 24.84) and in 2019 (AOR 16.16, 95% CI 8.95 − 30.69) compared with 2017, age under 9 months (AOR 2.60, 95% CI 1.64 − 4.20) and age from 9 months to 5 years (AOR 1.75, 95% CI 1.12 − 2.81) compared with less than 9 months, unvaccinated (AOR 1.56, 95% CI 1.09 − 2.27) compared with vaccinate, hospital transmission (AOR 1.73, 95% CI 1.38 − 2.17) compared to transmission in community, residing 20 to 200 kilometres (AOR 1.46, 95% CI 1.12 − 1.93) and 200 to 500 kilometres from VNCH (AOR 1.75, 95% CI 1.21 − 2.53) compared with less than 20 kilometres, hospitalized 7 to 21 days (AOR 3.50, 95% CI 2.79 − 4.40) and more than 21 days (AOR 14.20, 95% CI 9.15 − 22.04) compared with less than 7 days, having underlying condition in gastrointestinal system (AOR 6.74, 95% CI 1.65 − 28.87). Year of admission, age group, vaccination status, place of exposure, distance from VNCH, region of residence, clinical outcome, duration between onset and admission, duration of stay within VNCH had overall significance (p < 0.001). The association between duration between onset and test and severity had overall significance (p = 0.017 by likelihood ratio test), but not within an individual category. Factors no longer statistically significant in the multivariate analysis: having medical conditions in cardiovascular system, kidney and urology system, or neurological system, absence of underlying medical or other medical conditions, being tested 48 hours after onset and being diagnosed with measles, or bronchopneumonia, having RSV, adenovirus, pertussis or other co-infections. Sex, an unknown number of vaccinations, having inherited metabolic disorders, or the presence of an unknown place of exposure were not significantly associated with mortality in the uni- or multivariable analyses.
Table 4 illustrates the relationships between socio-demographic and clinical factors and mortality. In the univariate analysis, the following factors were significantly associated with fatality: age 9 months to 5 years, community transmission, living 20 to 500 kilometres from VNCH, living in the Northeast, Northwest or Central Vietnam, being admitted to the hospital before the onset or after 2 weeks, hospitalized for more than a week, having no underlying diseases or other underlying con-ditions, requiring a form of respiratory support (oxygen supplementation via nasal cannula, CPAP, conventional mechanical ventilation, HFO ventilation). Furthermore, having adenovirus or other healthcare-related infections, being severe, having sepsis/septic shock sequelae, and being diagnosed with measles, bronchopneumonia, or other diseases on admission are all substantially linked with death. Children with the following risk factors had significantly higher/lower odds of mortality in the multivariable analysis: age under 9 months (adjusted OR [AOR] 0.21, 95% CI 0.07 − 0.70) and age from 9 months to 5 years (AOR 0.19, 95% CI 0.06 − 0.60) compared with less than 9 months, residing 20 to 200 kilometres from VNCH (AOR 4.09, 95% CI 1.12 − 23.19), having adenovirus (AOR 8.29, 95% CI 3.31−20.45) or other hospital acquired infections (AOR 7.20, 95% CI 1.14−40.39). The association between place of exposure and mortality had overall significance (p = 0.001 by likelihood ratio test), but not within an individual category. The duration between onset and admission, as well as duration of stay within the hospital, had overall significance (p = 0.015 and p < 0.001, respectively). Factors such as age less than 9 months, hospital transmission, hospitalization for more than one week, absence of underlying medical or other medical conditions, living 20 to 200 km away from VNCH and being diagnosed with measles, bronchopneumonia, or others on admission in the univariable analysis were no longer statistically significant in the multivariate analysis. Sex, admission year, having bronchopneumonia or the presence of an unknown place of exposure were not significantly associated with mortality in the uni- or multivariable analyses.
4. Dicussions
In this study, we present the epidemiological and clinical characteristics of more than 2,000 children hospitalized with measles in Hanoi, the vast majority during Vietnam’s 2017-19 epidemic. To the best of our knowledge, this study is one of the few that provides data on the number of severe cases. According to the caregiver report or immunization card, 1,808 (87.3%) had not received any MCV. More than half of the individuals in our study reported the community as the likely place of infection exposure, suggesting secondary cases. 69.3% of patients live in the Red River Delta, with 34.5% living in Hanoi, Vietnam’s second most densely populated region, with a population of 2209, 2239, and 2410 per square kilometre in 2017, 2018, and 2019 respectively [16, 17, 18]. Crowded environments and high concentrations of people have known risk factors for disease occurrence and propagation, and Hanoi has been cited as a hotspot of prior epidemics in the country [19, 20]. During the epidemic, children admitted from further distances are most likely referrals from other hospitals due to complications. Despite the high number of hospital admissions and the severity rate, the mortality rate of 1.4% discovered in this study of patients hospitalized with measles is relatively low lower than the global norm. A systematic review and modelling analysis, encompassing various studies, reported an overall pooled mean case fatality rate for measles of 2.2%, with an in-hospital mortality rate ranging from 0.9% to 6.0%, higher than the rate observed in this study [21]. The lower in-hospital mortality rate observed here may be attributed to several factors, including increased measles vaccination coverage, heightened public awareness and recognition, a higher rate of hospitalization, and prompt treatment. Additionally, the fact that the study was conducted at a national-level hospital, where patients with severe symptoms receive appropriate medical care, may have contributed to the lower mortality rate.
The median age of children was under 10 months. Children aged 9 months to 5 years accounted for 12 (60%) of the 30 deaths. To boost MCV1 coverage and reduce cases and deaths in this age range, routine immunization must be strengthened. The strategy is less obvious for the 40% of deaths in infants under the age of 9 months, before their scheduled MCV1. Measles vaccination at 4 to 5 months of age was beneficial in preventing serologically confirmed and definite clinical measles during an outbreak in Guinea-Bissau in 2003-2004 [22]. Vietnam recommends measles vaccination for children from 9 months and in late 2018, 4.2 million children aged 1 to 5 received measles vaccine as part of an additional measles-rubella immunization campaign [9, 23]. In the event of future epidemics, these measures must be swift and effective. The earlier timing of MCV in low-income settings where maternal antibody levels may be low requires further evaluation. Upon admission, 849 (41.0%) children have bronchopneumonia. Most children had complications, with bronchopneumonia, the most common measles complication, observed in two-thirds of cases (68.2%). This prevalence of bronchopneumonia was higher than that reported in Greece (46.2%), Israel (30.5%), China (33%), and Serbia (26.8%), although these studies included adolescent cases, which could have affected the aggregated results [24, 25, 26, 27]. Tu et al. reported that measles-associated bronchopneumonia was associated with younger age, which explains the higher rate of bronchopneumonia in our study [26]. There was no additional information on the pathogens causing pneumonia or gastroenteritis, or the nature of the neurological complications.
There were no statistically significant differences in the risk of severe complications or mortality between the two genders. The odds of death increased by more than three times for children living 20 to 200 kilometres away from the hospital and by more than four times for children living 200 to 500 kilometres away, likely reflecting the increment in travel time to VNCH and referrals from other healthcare institutions due to complications. Children with bronchopneumonia complications had a greater risk of death, as expected and consistent with earlier research [28, 29]. Compared to children who had received one or more doses of MCV, the odds of severity were more than three times greater in unvaccinated children. Vaccine-modified measles is characterized by a milder illness in children with preexisting protection following vaccination [30, 31]. Patients with community transmission have lower risk of severe forms of disease and mortality rates than patients exposed to other sources. The severity rate increased fivefold and the mortality rate increased more than tenfold with adenovirus infection.
It is not a question of if, but when the next measles epidemic will occur in Vietnam and other countries. Reported measles cases have reduced from the last six months of 2020 to 2022, likely due to COVID-19 social distancing measures, reduced domestic and international travel, disturbance in surveillance and reporting, and the phase of the measles epidemic cycle [7]. Furthermore, the Vietnam Ministry of Health listed measles on the list of group B infectious diseases requiring medical isolation as well as expanded vaccination for subjects aged 6-9 months, combined with other reasonable measures to control epidemics from hospitals, limiting the spread of measles in hospitals [32]. However, at the same time, the COVID-19 pandemic has significantly disrupted routine vaccination practices, with estimates of more than 90 million children missing measles vaccine doses by October 2020, and reports of a lost generation of children [33]. MCV1 coverage fell from an average of 85% during 2015-2019 to 84% in 2020 and 81% in 2021, the lowest coverage levels for MCV1 since 2008. The most significant declines in MCV1 coverage during 2019-2021 occurred in the South-East Asia Region (from 94% to 86%) and the Western Pacific Region (from 95% to 91%). MCV2 coverage grew from 63% to 71% between 2015 and 2019, suggesting second-dose introductions in many countries. However, MCV2 coverage remained stable thereafter (72% in 2020 and 71% in 2021), with only four additional countries deploying MCV2 during 2020-2021 [34]. Vaccination coverage in Vietnam declined further to 89% for MCV1 and 85% for MCV2 in 2021 [35]. When COVID-19 limitations are lifted, an increasing number of unvaccinated children who are susceptible to measles creates an environment for measles to resurface [36]. Repeated occurrences of measles outbreaks in Vietnam highlight the significance of re-evaluating measles virus dynamics and prevention and control strategies. Actions are required to reach unimmunised children through catch-up campaigns and to prepare for anticipated outbreaks. Governments must not lose sight of the Measles and Rubella Strategic Framework, 2021-2030 [37].
This study has certain limitations. The patients enrolled in this study are those that were admitted to our hospital, whose illness may be more severe and are not representative of all measles cases in Vietnam. Nonetheless, since patients come from various regions across the country, they still offer insights into the broader measles landscape in Vietnam. Detailed data were partially lacking on the number of MCVs and the place of exposure though the medical files were thoroughly checked, and the percentage of missing data was quite low for almost all variables. Nonetheless, we were able to analyze a large sample of over 2,000 measles cases with extensive clinical information and other socio-demographic characteristics using the measles suspected case inquiry form.
In conclusion, in this study, we report clinical and epidemiological characteristics from a significant measles outbreak that primarily affected unvaccinated children admitted to VNCH with significant number of deaths occurred among infants aged less than nine months. As MCV1 coverage in Vietnam declined to 89% in 2021, routine immunization needs to be strengthened and the earlier timing of MCV1 requires further evaluation to prevent further outbreaks.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Footnotes
↵* Co-first authors
↵** Co-senior authors. Email addresses: h.a.ngo{at}sms.ed.ac.uk (Hoang-Anh Ngo), phucph{at}nch.gov.vn (Phuc H Phan).
The groups within certain categorical variables of the original dataset have either been renamed or merged to enhance understanding upon reading the manuscript. Small changes in wording of the manuscript have also been made. No major transformation to the original dataset for univariate/multivariate analysis or the manuscript has been conducted for this revision.
References
