Introduction
The world is facing a devastating pandemic of a novel coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Singhal, Reference Singhal2020). The virus emerged and the outbreak was first identified in the disease epicentre in Wuhan, Hubei province of China, in late November 2019. The World Health Organization (WHO) 2020 declared the outbreak to be a Public Health Emergency of International Concern (PHEIC) on 30th January 2020 and recognized it as a pandemic on 11th March 2020 (World Health Organization [WHO], 2020a).
As of August 6th 2020, a total of 18,354,342 confirmed cases of COVID-19 have been reported in 213 countries and territories, resulting in a death toll of 696 147. For instance, more than 11,302,021 people have since recovered (World Health Organization [WHO], 2020a). Egypt has been experiencing widespread transmission of the virus in the community for several months, with an increasing number of COVID-19 outbreaks in healthcare facilities (World Health Organization [WHO], 2020a) (Figure 1).
With the rapidly increasing numbers of infections, COVID-19 has become a considerable strain for all health care settings. So far, we lack a rapid point of care test for SARS-CoV-2, and the diagnosis of COVID-19 is limited by the use of reverse transcription polymerase chain reaction (RT-PCR) which is not widely available or more time-consuming, thus not practicable at primary health care (PHC) settings. Moreover, expanding testing supplies and capabilities is lagging behind what is needed to curb such unprecedented pandemic scale. Clinical prediction scores support the assessment of patients in the PHC setting to determine the need for further diagnostic and therapeutic steps (Sheridan, Reference Sheridan2020).
Early in the epidemic, the Egyptian Ministry of Health and Population (MoHP) introduced a self-assessment score as remote care solution to support emergency services in the fight against COVID-19 (Figure S1). It combines specific clinical and epidemiological features to determine an individual’s risk of having COVID-19. The MoHP has released the score for the public in the local language as a free online tool supported by a disease control hotline. Users can access the tool through media outlets and the MoHP official website. Based on their input, patients can have telephone consultations with a healthcare professional to determine the need for testing or be guided for home self-care (The Egyptian Ministry of Health and Population (MoHP), 2020a).
We evaluated the MoHP COVID-19 prediction score in an analysis of 818 consecutive patients with suspected COVID-19 from mid-February until mid-July 2020.
Methods
Study setting, design and population
The study was conducted at a PHC clinic in Alexandria between February and July 2020. Patients presenting with acute respiratory symptoms were enrolled consecutively and evaluated in order to determine the initial features which may help to distinguish probable COVID-19 cases from other respiratory problems. All patients were evaluated with chest imaging and blood testing. Detection of SARS-CoV-2 in nasopharyngeal swabs by RT-PCR test was done at a central reference laboratory.
Due to variations of the disease occurrence over time, study sample size estimations were not performed, and alternatively, sample size was determined by practical convenience.
The case definition of a COVID-19 suspected or confirmed case was set by the MoHP and was based on the encounter of elements of the epidemiological history as well as clinical symptoms (Table 1) (The Egyptian Ministry of Health and Population (MoHP), 2020b).
Scoring tools
To assess the validity of the MoHP self-assessment score, we applied it for all the enrolled suspected COVID-19 patients and tested its performance versus the finally obtained diagnosis. The score includes nine questions with each criterion scored with a different number of points. Q1: Fever (≥ 38.0°C) (2 points), Q2: Severe cough (2 points), Q3: Severe sore throat (1 point), Q4: Vomiting or diarrhea (0 point), Q5: Chronic disease (diabetes mellitus, hypertension, ischemic heart disease, renal disease, liver disease, autoimmune disease, etc….) (1 point), Q6: Travel in or outside Egypt (Sharm Elshiek, or Europe, China or any country) (5 points), Q7: Direct contact to a respiratory infection case (care, face to face) (4 points), Q8: Visiting health care place that had received confirmed case (3 points), Q9: Being a health care worker (HCW) (2 points).
A numeric score is calculated by counting the number of criteria met at the initial presentation (score range 0–20 points). A cut-off of 6 points and more was highly suggestive for COVID-of infection and was set for calling the disease control hotline, whereas a score of 4 necessitated home isolation and a score of 5 was required before consulting a health care physician.
Statistical analysis
The collected data were reviewed for accuracy and integrity and entered into computer software. For all statistical analyses, we used MS Excel 2016 and a statistical software package (IBM® SPSS® Statistics Base, IBM Corp: Armonk, NY; version 21.0). We afterward summarized all variables. Continuous variables are presented as the mean ± standard deviation (SD). Categorical variables are expressed as numbers with proportions, n (%). Chi-Square test, Fisher’s exact tests, Student’s t-tests, multivariate logistic regression and ROC curve analyses were performed for group comparisons.
The association between patient’s variables (clinical and sociodemographic) and having a SARS-CoV-2 positive PCR test was examined using a univariate analysis. We run this step three times to infer the association in early epidemic phase, late epidemic phase and throughout the epidemic. The variables associated with SARS-CoV-2 positive PCR test at the P < 0.05 significance level on univariate analyses and deemed potentially useful for the clinico-epidemiological prediction of the disease were selected for multivariable analysis. The final list of eligible variables differed in relation to the epidemic phase. The selected variables were entered as covariates to develop a multivariable logistic regression equation by backward (Wald) step-wise elimination, with SARS-CoV-2 positive PCR test being the outcome variable.
Testing the validity and accuracy of the prediction scores
After testing the MoHP scoring system (thereafter score 1), we created two new scores named ‘score 2’ and ‘score 3’. In score 2, vomiting or diarrhea symptoms in Q4 worth 1 point, and score 3 additionally included patient age [20 years (0 point), 20–35 (1 point), 36–55 (2 points) and 56–99 (3 points)] and smoking status [non-smoker/ex-smoker (0 point) and current smoker (1 point)]. The aggregate of these weighted variables was expressed as a total score (prediction index) for each patient individually. A receiver operating characteristic (ROC) curve was plotted with the total score as the test variable and SARS-CoV-2 PCR test as the state variable. The area under the ROC curve (AUC) was used to assess its overall predictive performance, sensitivity, specificity, positive and negative predictive values. We identified the optimal cut-off (Youden index) of the total score from the ROC curve as the point corresponding to the best trade-off between sensitivity and specificity.
We assessed the goodness-of-fit of all the prediction models using the Hosmer and Lemeshow test.
Results
Characteristics of the study cohort
Of a total of 818 suspected COVID-19 patients presenting with acute respiratory symptoms were enrolled in the study, 390 (47.7%) were confirmed to have SARS-CoV-2 by RT-PCR. The mean age of SARS-CoV-2 positive subjects was 43.0 ± 17.0 years, with women being more affected than men (53.8% vs. 46.8% respectively) (Table 2). The total scores met at presentation differed significantly (t = −8.5, P < 0.001) between SARS-CoV-2 positive (mean ± SD = 9.3 ± 3.9, standard error = 0.2) and SARS-CoV-2 negative (mean ± SD = 7.1 ± 3.2, standard error = 0.15) patients (Table 2).
DM = diabetes mellitus; HT = hypertension; IHD = ischemic heart disease; PCR = polymerase chain reaction; SD = standard deviation.
t, t value for student t-test.
p significant at < 0.05.
Performance of the prediction scores
The validation sample was split up in an early and late phase of the epidemic as driven by the results of data analysis. The performance of the three tested scores (1, 2 or 3) was comparable regarding the prediction of SARS-CoV-2 infection when tested early in the epidemic (February–April), late in the epidemic (May–July) as well as throughout the epidemic period (February–July) (Figure 2). A cut-off of 6.0 points met in the initial score (score 1) was associated with correct RT-PCR-based SARS-CoV-2 detection in nasopharyngeal swabs and presence of symptomatic COVID-19 in 296 out of 390 patients (sensitivity = 75.9%, specificity = 42.3%, positive predictive value = 54.5%, negative predictive value = 65.8%), whereas lower scores were associated with other respiratory conditions (bronchitis, pneumonia, common cold). When we adopted another cut-off point (8.5) with a better trade-off between sensitivity and specificity, more accurate score performance was obtained (sensitivity = 69.5%, specificity = 63.1%, positive predictive value = 63.2%, negative predictive value = 69.4%). Collectively, the tested prediction scores had a sensitivity range of 60.0%–89.3%, a specificity range of 48.6%–65.7% and corresponding positive and negative predictive value ranges of 52.5–65.2% and 57.6%–88.0%, respectively. Likewise, the interrater reliability of the initial and the derived scores versus the diagnostic PCR test used for confirming the diagnosis of COVID-19 showed minimal degrees of agreement, with kappa values ranging between 0.204 and 0.388 (P < 0.001). The AUC ranged between 0.641 (95% CI = 0.589–0.692) and 0.769 (95% CI = 0.721–0.817), (P < 0.001). Overall, the performance of the triage scores was more accurate in predicting SARS-CoV-2 infection later in the epidemic (Figure 2).
All the derived models showed appropriate goodness-of-fit (Hosmer and Lemeshow P > 0.05).
Predictors of COVID-19 infection
On univariate analysis, the variables associated with the presence of COVID-19 infection in patients presenting with acute respiratory symptoms at the P < 0.05 significance level differed according to the epidemic phase (Tables 2 and 3). Step-wise elimination generated epidemic phase wise predictor variables in the final multivariable equations. Early in the epidemic, history of travel in or outside Egypt [OR (95% CI) = 12.1 (5.0–29.4)] and direct contact with a respiratory infection case [OR (95% CI) = 3.5 (2.1–5.6)] were strong predictors of COVID-19 infection. On the other hand, having severe cough [OR (95% CI) = 0.29 (0.18–0.48)] or visiting health care place that had received COVID-19 confirmed cases [OR (95% CI) = 0.73 (0.19–0.74)] reduced the likelihood of having COVID-19. Later in the epidemic, having fever [OR (95% CI) = 28.1 (3.7–214.0)], severe cough [OR (95% CI) = 3.7 (1.4–9.8)], severe sore throat [OR (95% CI) = 2.9 (1.5–5.8)], direct contact with a respiratory infection case [OR (95% CI) = 4.3 (1.7–10.9)], and being a HCW [OR (95% CI) = 5.8 (2.8–11.9)] increased the odds of having COVID-19 infection. Throughout the epidemic period, older age [OR (95% CI) = 1.01 (1.0–1.02)], history of travel [OR (95% CI) = 10.6 (5.5–20.4)], direct contact with a respiratory infection case [OR (95% CI) = 3.1 (2.1–4.4)], and being a HCW [OR (95% CI) = 2.2 (1.4–3.4)] increased the risk of having SARS-CoV-2 infection, whereas having severe cough was identified as a negative disease predictor [OR (95% CI) = 0.50 (0.34–0.73)] (Table 4).
DM = diabetes mellitus; HT = hypertension; IHD = ischemic heart disease; PCR = polymerase chain reaction; SD = standard deviation.
t, t value for student t-test.
p significant at < 0.05.
C.I. = Confidence interval; df = degree of freedom.
p significant at < 0.05.
Discussion
Based on the evaluation of the initial data of 818 patients presenting with acute respiratory symptoms, we assessed the diagnostic performance of the MoHP score as a simple self-assessment tool for COVID-19 prediction. The score combines consistently available clinico-epidemiological elements and thus it is ideal for triage in PHC settings.
During the COVID-19 pandemic and under severe time constraints, the MoHP score was released for the public as a virtual assessment tool for selecting symptomatic COVID-19 cases from rather unspecific general or respiratory symptoms and prioritize those with a higher epidemiological probability of the disease for further medical attention. Based on personal symptoms, exposure and risk factors, individuals can find the right care by either assurance, home self-care, consulting a physician or immediately calling the disease control centre to get a SARS-CoV-2 swab. However, since the tool was created very early in the epidemic with little accumulation of relevant clinical, there was no clear or definitive data available regarding the most specific frequent symptoms encountered in SARS-CoV-2 infection, the tool lacked some important symptoms frequently encountered by COVID-19 patients including altered taste or smell, headache, confusion, difficult breathing, myalgia and fatigue (Harapan et al., Reference Harapan, Itoh, Yufika, Winardi, Keam, Te, Megawati, Hayati, Wagner and Mudatsir2020, Lüers et al., Reference Lüers, Klußmann and Guntinas-Lichius2020, Cummings et al., Reference Cummings, Baldwin, Abrams, Jacobson, Meyer, Balough, Aaron, Claassen, Rabbani, Hastie, Hochman, Salazar-Schicchi, Yip, Brodie and O’donnell2020). It is worth noting that difficult breathing may have not been considered in the scoring system because symptoms usually begin as mild in all patients with occasionally dyspnea manifesting in severe cases (Cascella et al., Reference Cascella, Rajnik, Cuomo, Dulebohn and di Napoli2020). Although the score included gastrointestinal symptoms, a common presentation in COVID-19 disease (Hajifathalian et al., Reference Hajifathalian, Mahadev, Schwartz, Shah, Sampath, Schnoll-Sussman, Brown, Carr-Locke, Cohen and Sharaiha2020, Goyal et al., Reference Goyal, Choi, Pinheiro, Schenck, Chen, Jabri, Satlin, Campion, Nahid, Ringel, Hoffman, Alshak, Li, Wehmeyer, Rajan, Reshetnyak, Hupert, Horn, Martinez, Gulick and Safford2020), it ignored its scoring. However, when a history of vomiting or diarrhea worth 1 point in our proposed score 2 and score 3, this did not result in an appreciable change in the performance of the initial score.
Although relying on respiratory symptoms was not appropriate to identify individuals with COVID-19 early in the disease, epidemiological risk factors such as travel history and direct exposure to a symptomatic case were strong predictors of SARS-CoV-2 infection. This suggests that the clinical differentiation of COVID-19 from other respiratory illnesses at the early stages of the disease is limited. One reason could be that the symptoms included in the scoring system were relatively unusual although they appear to be much more common in other countries, a difference which may reflect geographic variation or differential reporting. Another possibility is that travelers returning from the disease epicentres in Europe and the hardest-hit regions (for our study cohort; China, Italy, Spain, Germany, USA and KSA) at that time were more attentive of their symptoms and sought medical advice immediately before they develop overt symptoms.
Late in the epidemic, the risk attributable to travel disappeared from the model which reflects international travel bans during this period. Conversely, the disease started to affect more HCWs. The risk increased by 5.8-fold compared to the early epidemic phase (30.2% vs. 10.4% respectively, P < 0.001). This could also reflect that they had been tested earlier and more frequently.
From the beginning of the epidemic, older age has been identified as an important risk factor for disease susceptibility, severity and outcome as explained by the age-related abundance of the angiotensin converting enzyme 2 (ACE2) receptors that may facilitate viral invasion (Guan et al., Reference Guan, Ni, Hu, Liang, Ou, He, Liu, Shan, Lei, Hui, Du, Li, Zeng, Yuen, Chen, Tang, Wang, Chen, Xiang, Li, Wang, Liang, Peng, Wei, Liu, Hu, Peng, Wang, Liu, Chen, Li, Zheng, Qiu, Luo, Ye, Zhu and Zhong2020, Huang et al., Reference Huang, Wang, Li, Ren, Zhao, Hu, Zhang, Fan, Xu, Gu, Cheng, Yu, Xia, Wei, Wu, Xie, Yin, Li, Liu, Xiao, Gao, Guo, Xie, Wang, Jiang, Gao, Jin, Wang and Cao2020, Livingston and Bucher, Reference Livingston and Bucher2020, McMichael et al., Reference Mcmichael, Currie, Clark, Pogosjans, Kay, Schwartz, Lewis, Baer, Kawakami, Lukoff, Ferro, Brostrom-Smith, Rea, Sayre, Riedo, Russell, Hiatt, Montgomery, Rao, Chow, Tobolowsky, Hughes, Bardossy, Oakley, Jacobs, Stone, Reddy, Jernigan, Honein, Clark and Duchin2020). In the present report, advancing age was associated with one-fold higher odds of being COVID-19 positive throughout the epidemic period. Nevertheless, including patient age in score 3 did not improve its accuracy in predicting SARS-CoV-2 infection.
Likewise, smoking neither increased the risk of COVID-19 nor improved the accuracy of the initial prediction score. Indeed, smoking as a risk factor for SARS-CoV-2 infection is under discussion although pooled analysis linked it to increased disease severity and death in hospitalized COVID-19 patients (World Health Organization [WHO], 2020b, Grundy et al., Reference Grundy, Suddek, Filippidis, Majeed and Coronini-Cronberg2020, Engin et al., Reference Engin, Engin and Engin2020). Furthermore, although COVID-19 infections are relatively more likely to be seen in men (Park, Reference Park2020), the risk was comparable in both sexes throughout the epidemic phases. So, it may be difficult to identify who will more likely develop the illness by using some established risk factors such as age, sex, smoking and comorbidities. Nevertheless, as the COVID-19 pandemic evolves, the accumulation of relevant clinical and other health information will likely lead to the development of a more robust risk prediction model and increase the success of assessment strategies to support medical decision making.
The structured questions in the MoHP score did not include a focus on change in symptoms or assessing clinical deterioration in mild disease, which might have a role in the clinical assessment.
Given its modest accuracy and to avoid misclassifying patients, the tested prediction tool could be considered as one component within a broader clinical evaluation of patients presenting with acute respiratory symptoms since for instance, it gives attention to a trajectory of symptoms and patient characteristics such as age, exposure and comorbid conditions.
So far, there is no reliable risk assessment tool for predicting COVID-19 in daily practice. Wynants et al. (Reference Wynants, van Calster, Collins, Riley, Heinze, Schuit, Bonten, Dahly, Damen, Debray, de Jong, de Vos, Dhiman, Haller, Harhay, Henckaerts, Heus, Kreuzberger, Lohmann, Luijken, Ma, Martin, Andaur Navarro, Reitsma, Sergeant, Shi, Skoetz, Smits, Snell, Sperrin, Spijker, Steyerberg, Takada, Tzoulaki, van Kuijk, van Royen, Verbakel, Wallisch, Wilkinson, Wolff, Hooft, Moons and van Smeden2020) conducted a systematic review of several COVID-19 prediction models ranging from rule-based scoring systems to advanced machine learning models that have been rapidly proposed and published in response to the pandemic. Most of these models were found of limited quality, at high risk of bias and inadequate in making precise predictions about disease risk and outcomes (Wynants et al., Reference Wynants, van Calster, Collins, Riley, Heinze, Schuit, Bonten, Dahly, Damen, Debray, de Jong, de Vos, Dhiman, Haller, Harhay, Henckaerts, Heus, Kreuzberger, Lohmann, Luijken, Ma, Martin, Andaur Navarro, Reitsma, Sergeant, Shi, Skoetz, Smits, Snell, Sperrin, Spijker, Steyerberg, Takada, Tzoulaki, van Kuijk, van Royen, Verbakel, Wallisch, Wilkinson, Wolff, Hooft, Moons and van Smeden2020). It is worth noting that the currently available data suggest suboptimal test performance for the gold standard test “RT-PCR” since it detected only the SARS-CoV-2 virus in 63% of nasal swabs and 32% of pharyngeal swabs in patients with the known disease (Wang et al., Reference Wang, Xu, Gao, Lu, Han, Wu and Tan2020).
The proposed tool is well-suited to carry out a preliminary assessment of suspected patients and help them to get timely treatment and quarantine suggestion. Since COVID-19 is changing the face of health care delivery, the adaptation of this tool could be useful in virtual primary care or community settings. The predictors identified in included models should be considered as candidate predictors for new models. However, to accommodate an evolving COVID-19 pandemic, the proposed model will need to be recalibrated and refit over time. Widespread dissemination, tracking of its utilization, and direct integration with the digital health record can further improve its utility.
This type of self-assessment tool provides immediate benefit to the patients and health care providers as we face anticipated increased demand and limited resources. They undoubtedly have the potential to (i) facilitate patient triage and preventing unnecessary in-person visits during the COVID-19 pandemic to prevent patient exposure to pathogens in their route to clinic visits and in waiting rooms, (ii) avoid unnecessary testing of a large proportion of patients who are indeed COVID negative, (iii) reduce personal protective use by clinic staff and liberate front-line staff to care for sicker patients. However, it remains unclear how many patients might be harmed if the tool falsely predicted that self-care was sufficient, and thus limit its use in the context of a pandemic (Jehi et al., Reference Jehi, Ji, Milinovich, Erzurum, Rubin, Gordon, Young and Kattan2020, Wynants et al., Reference Wynants, van Calster, Collins, Riley, Heinze, Schuit, Bonten, Dahly, Damen, Debray, de Jong, de Vos, Dhiman, Haller, Harhay, Henckaerts, Heus, Kreuzberger, Lohmann, Luijken, Ma, Martin, Andaur Navarro, Reitsma, Sergeant, Shi, Skoetz, Smits, Snell, Sperrin, Spijker, Steyerberg, Takada, Tzoulaki, van Kuijk, van Royen, Verbakel, Wallisch, Wilkinson, Wolff, Hooft, Moons and van Smeden2020).
Conclusion
An updated implementation of this home-grown tool holds promise and provides new insights into the rapid diagnosis, the timely isolation and treatment of COVID-19, which are currently the known most effective response. Given its potential implications in terms of global health and in improving the COVID-19 response at the PHC level, the proposed tool could be shared with other health organizations and authorities in other countries. The tool can be deployed in a web-based or smartphone contexts as an interactive clinical assessment digital tool or a Self-Checker paired with resources and recommended actions. In this regard, it can be customized to fit the needs of both the individual and the community. Based on user’s responses, the tool can provide valuable information on seeking medical care specific to his area. Moreover, it can track geographic trends by creating heat map zones in cases when users share their zip codes. A further credible use of this tool is to provide passengers entering from risky countries to check their symptoms daily for 14 days. This will ultimately help in combating the spread of the disease globally.
Limitations of the study
The major limitation of this work is the single-centre evaluation of only a limited number of patients representing a convenience sample. Early in the epidemic, the MoHP score was introduced to predict whether a user of the tool should seek a particular level of care. However, we evaluated how this model predicts the probability of a current COVID-19 diagnosis. Indeed, the introduced MoHP tool predicts three triage levels (the probability of COVID-19 infection, the requirement for home isolation and consultation). Yet, these outcomes were not evaluated in the main analysis, which revolved around the distinction between COVID-19 positive and negative outcomes only. It is noteworthy that the introduced self-assessment tool is supposed to be used by individuals in the general population, while we evaluated its performance in individuals that have already presented themselves at PHC. Thus, the selected individuals do not seem to match the research question of interest. The fact that the tool was not evaluated in the target population might result in a bias (self-selection bias) that cannot be corrected for the analysis phase of this study. We must also acknowledge that the clinical motivation for the different cut-off was not described due to limited access to the baseline data used for creating the score. Moreover, the application of the proposed scores on the identification of COVID-19 cases with atypical symptoms was not confirmed which limits their use in clinical practice.
While the preliminary results were promising, it is obvious that the MoHP score needs to be thoroughly updated and validated in further large-scale studies conducted across different populations in larger cohorts of patients to exploit our findings and gain more reliable data regarding its diagnostic yield. For instance, the score may be used for the assessment of patients with suspected COVID-19, but with some caution undertaken in this context. We did not keep track of the disease outcomes in the study cohort and therefore we could not assess the appropriateness of our tool’s recommendations. Like the majority of the currently available COVID-19 diagnostic models, we used viral nucleic acid test results as the gold standard, which may have unacceptable false negative rates (Wynants et al., Reference Wynants, van Calster, Collins, Riley, Heinze, Schuit, Bonten, Dahly, Damen, Debray, de Jong, de Vos, Dhiman, Haller, Harhay, Henckaerts, Heus, Kreuzberger, Lohmann, Luijken, Ma, Martin, Andaur Navarro, Reitsma, Sergeant, Shi, Skoetz, Smits, Snell, Sperrin, Spijker, Steyerberg, Takada, Tzoulaki, van Kuijk, van Royen, Verbakel, Wallisch, Wilkinson, Wolff, Hooft, Moons and van Smeden2020). The strengths are its simplicity, immediate availability as well as wide applicability due to simple components.
Supplementary material
For supplementary material accompanying this paper visit https://doi.org/10.1017/S1463423621000736
Ethical consideration
Ethical review board and the ethics committee
The study was approved by the institutional review board and the ethics committee of the High Institute of Public Health affiliated with Alexandria University, Egypt
Ethical considerations
The study was conducted in accordance with the international ethical guidelines and of the Helsinki’s Declaration (2013). Data sheets were coded by number to ensure anonymity and confidentiality of the participants’ data.
This article does not contain any studies with animals performed by any of the authors.
Consent to participate
Informed written consent was obtained from each participant after explaining the aim and concerns of the study.
Conflicts of interest
All authors declare that they have no conflict of interest.
Funding
Authors own work.
Consent for publication
All authors approved the manuscript for publication.
Availability of supporting data
All data are fully available without restriction by the corresponding author at [email protected] the and through the public data repository “Harvard Dataverse” at https://dataverse.harvard.edu/dataverse
Acknowledgement
We would like to thank all the participants for accepting to participate in the study.
Author contributions
EWAW: Conceptualization, developed the theoretical framework and study design, conducted the data base search, data extraction, development and validation of the data collection tools, took the lead for overall direction and planning, supervised the study implementation, data curation, analysis and interpretation of data, major contribution to writing, revised and approved final version of the manuscript
MM: Study implementation, recruitment of the study participants, data collection, interpretation of data, contributed to writing of the manuscript, revised and approved final version of the manuscript.