Healthcare workers’ willingness to respond following a disaster: a novel statistical approach toward data analysis

Background The willingness of healthcare workers (HCW) to respond is an important factor in the health system’s response capacity during emergencies. Although much research has been devoted to exploring this issue, the statistical methods employed have been predominantly traditional and have not enabled in-depth analysis focused on absenteeism-prone employees during emergencies. The present study employs an innovative statistical approach for modeling HCWs’ willingness to respond (WTR) following an earthquake. Methods A validated questionnaire measuring knowledge, perceptions, and attitudes toward an earthquake scenario was distributed among Israeli HCWs in a hospital setting. Two regression models were employed for data analysis – a traditional linear model, and a quantile regression model that makes it possible to examine associations between explanatory variables across different levels of a dependent variable. A supplementary analysis was performed for selected variables using broken line spline regression. Results Females under the age of forty, and nurses were the most absenteeism-prone sub-groups of employees (showed low WTR) in earthquake events. Professional commitment to care and perception of efficacy were the most powerful predictors associated with WTR across all quantiles. Both marital status (married) and concern for family wellbeing, designated as statistically significant in the linear model, were found to be statistically significant in only one of the WTR quantiles (the former in Q10 and the latter in Q50). Gender and number of children, which were not significantly associated with WTR in the linear model, were found to be statistically significant in the 25th quantile of WTR. Conclusions This study contributes to both methodological and practical aspects. Quantile regression provides a more comprehensive view of associations between variables than is afforded by linear regression alone. Adopting an advanced statistical approach in WTR modeling can facilitate effective implementation of research findings in the field. Electronic supplementary material The online version of this article (10.1186/s12909-019-1561-7) contains supplementary material, which is available to authorized users.


Background
Earthquakes are liable to become mass casualty events [1,2]. Healthcare workers (HCW) in both the pre-hospital and hospital settings are key players in the response to such events. They are required to provide lifesaving medical treatment to an influx of patients, possibly in the face of shortages in resources such as medical supplies and qualified manpower [3][4][5][6]. The State of Israel has contended with numerous mass casualty events caused by terrorism [7], but has little experience in dealing with natural disasters. Israel is located in a seismic risk area, and local experts warn that a powerful earthquake could well occur in the near future [8]. Enhancing the preparedness of its healthcare institutions and staff to deal with such a scenario is therefore seen as being of vital importance.
The clinical staff is directly responsible for fulfilling the core mission of the hospitalproviding quality patient careand is thus one of its key assets [9][10][11]. This is true for routine care and even more so during mass casualty events that call for patient care surge capacity [5,12]. Absenteeism of essential HCWs during such events may result in delays in providing medical aid for casualties, increased morbidity, unnecessary complications, and possibly even loss of life [3,13]. Much research has been devoted to exploring HCWs' willingness to respond (WTR) during and following public health emergencies, and great efforts have been made to reduce absenteeism [14][15][16]. Studies from the field of human resource management suggest that developing individual and collective competencies among core employees and addressing their needs during emergencies can prevent absenteeism and enhance the organization's resilience and capacity to contain such events [17,18].
Three comprehensive literature reviews have sought to identify personal-, professional-, organizational-, and threat-related factors associated with WTR among HCWs following a public health disaster. One review focuses exclusively on pandemic influenza [19] while the other two relate to 'disasters' in general [20,21]. Careful scrutiny of the results reveals conflicting findings regarding most factors. Only a few factors are in broad consensus, e.g. gender. Among the personal-related factors, the findings regarding gender are consistent and indicate that women tend to have a lower WTR than men [19,20]. Concern for family, loved ones, pets, and personal obligations such as childcare and eldercare emerge as powerful barriers against WTR [19][20][21]. Among the professional-related factors, the greater part of the evidence indicates that clinicians have higher WTR than non-clinicians [19], but at the same time that nurses have low WTR compared with other clinicians [19]. Belief in duty to care, knowledge, perception of efficacy, and levels of training are described as facilitators of WTR [19][20][21]. Organizational-related factors are not discussed at length, and it seems that confidence in the employer's ability to meet basic needs of employees (e.g. food, water, etc.) and other concerns is the primary factor in this category [19,21]. The findings regarding threat-related factors demonstrate that WTR rates are lower for events that pose higher personal risk to employees, such as infectious diseases and other biological, chemical, and radiological hazards, as compared with natural disasters and technological or other man-made threats [19,20]. This was also supported in a more recent study [22].
Although an extensive body of literature exists in this field, data analysis approaches have tended to be traditional, usually involving multivariate analysis based on linear or logistic regression. Consequently, the findings merely reflect general associations and cannot lead to identification of specific factors and characteristics associated with distinct levels of WTRinformation that is essential for ameliorating preparedness. WTR is usually measured using an ordinal scale (i.e. Likert scale). The use of linear regression is considered less than optimal when analyzing ordinal data since the latter cannot be assumed to be normally distributed; nevertheless the use of this method is prevalent in WTR modeling [23] and exposes the findings to possible criticism [24]. Another common method is dichotomizing ordinal data in order to apply logistic regression analysis [25,26], but this entrains loss of data, may lead to reduction in statistical power, and is generally not recommended [27]. A more flexible approach toward ordinal data analysis focused on the entire distribution of the dependent variable would provide far more useful and reliable information [28]. A possible candidate is quantile regression (QR) modeling. The use of QR was first introduced in economics [29], but it has since become popular in fields and disciplines such as medicine and health care [30,31], and recently in ecological research and community resilience modeling [32]. QR enables division of the dependent variable into homogeneous segments (much like the structure of original data), and the relationship with independent variables is assessed separately for each of the segments (quantiles). Thus, no loss of data occurs, and a holistic view of the results is made possible. While QR is concerned with modeling the dependent variables, attention should also be paid to the modeling methods used for the independent variables. Effect size measures (e.g. regression coefficient and odds ratio) provide information regarding the direction and strength of association, but fail to identify the exact points across the distribution of continuous explanatory variables at which a statistically significant change occurs in the outcome variable. This is not an issue if the nature of the association is linear, but may constitute a limitation in more complex relationships [33]. Spline regression is a smoothly joined piecewise regression that addresses this issue by performing division of a continuous independent variable into linear segments connected to each other at their edges (called a knot). The model makes it possible to test whether a statistically significant change in the regression line slope occurs above and below a specific knot location. The properties of this regression fit complex distribution patterns, enabling examination of both shape and strength of association [28,34].
Applying these analysis methods in WTR modeling can facilitate characterization of absenteeism-prone sub-groupsi.e. identification of personal characteristics related to HCWs with low levels of WTRas well as in-depth examination of the relationships between professional and organizational factors across different levels of WTR.
The objectives of this study were to identify the personal and professional characteristics of HCWs with low levels of WTR following an earthquake by using a novel statistical approach and help develop strategies to ameliorate their willingness.

Survey instrument
A structured self-administered anonymous questionnaire was designed to assess HCW demographic characteristics, WTR, attitudes, perceptions, and knowledge level regarding an earthquake scenario (see Additional file 1). The questionnaire was a modified version of a validated tool used to measure preparedness components among HCWs in various emergency situations, and was adjusted to an earthquake scenario. Additional information describing the survey design and validation process is detailed in previous studies [23,35,36]. Information pertaining to demographics included gender, age, marital status, number of children residing at home, and professional role. Knowledge regarding different components of hospital response to an earthquake was measured through twelve multiple-choice questions, based on a national standardized operating procedure (SOP) formulated by the Israel Ministry of Health. To understand HCW perceptions, seven questions were asked concerning self-and organizational efficacy in an earthquake scenario (α = .87), using a 7-point Likert-type scale ranging from 1 (I strongly disagree) to 7 (I strongly agree). WTR following an earthquake was assessed using two items that dealt with personal WTR, as well as the perceived willingness of colleagues. Possible barriers (2 items related to concern for family wellbeing and damage to house due to the earthquake) and facilitators (2 items related to professional commitment to care for the injured or ill, and fear of losing place of employment) of WTR were also assessed. The items in this section were also rated on a 7-point Likert-type scale. The questionnaire was pre-tested and piloted with a convenience sample of 27 HCWs who were similar in their characteristics to the study population. Based on the participants' recommendations, minor modifications were made in the final version.

Study design and data collection
This cross-sectional survey study was conducted from November 2013 through June 2014 in all secondary and tertiary hospitals in Israel (n = 24). The study population included clinical health professionals (i.e. physicians, nurses and paramedical staff) employed in the five hospital departments likely to be the most involved in providing treatment to earthquake casualties: the emergency department; internal medicine division; surgical wards; orthopedics; and operating rooms.
Sample size was determined before study initiation. Based on previous findings indicating that levels of WTR to different emergency situations among HCWs ranged between 58 and 72% [26,37,38], a mean value of 65% participants reporting high levels of WTR (ranking between 5 and 7 on a Likert scale) was assumed for calculating the current sample size. A margin of error of 5% and a 99% confidence interval were taken in account, and accordingly a sample of 604 HCWs was sought. In expectation of a relatively low response rate to the surveys, especially among hospital personnel, it was decided to increase the sample size to 2400 HCWs (anticipating a response rate of 25%).
Distribution of the survey and allocation of participants were carried out in several steps. First, the research team met with hospital emergency managers in the Emergency Division of the Israel Ministry of Health. After the team's objectives and study protocol had been presented, each hospital representative was given 100 questionnaires and asked to distribute twenty to the head nurse in each of the five departments mentioned. The head nurses also received a guideline sheet containing information about the study and detailed instructions regarding how and to whom to circulate the questionnaires (to six physicians, eight nurses, and six paramedical staff employed in their department; selection of specific participants in each category was according to the head nurse's choice). Next, the completed questionnaires were collected by the head nurse within 2 weeks from distribution and returned to the hospital's emergency manager. They were then sent back to the Ministry of Health and collected by the researchers. The importance of taking part in the survey was stressed at all levels by sending designated cover letters to hospital managers, heads of departments and head nurses and by including a personal appeal to the participants in the survey.

Statistical analysis
Data analysis was performed in three main steps. In the first step, a univariate analysis with descriptive statistics was followed by a bivariate analysis conducted using Spearman's correlation, chi-square test, Mann-Whitney U-test, Kruskal-Wallis one-way analysis of variance, and median test to estimate the associations between the dependent variable, namely WTR (defined as the mean score of the two items used to measure this element), and the other survey variables. In the second step. Two multivariate regression models (multivariate linear regression and multivariate quantile regression) were developed in order to identify associations between explanatory variables and WTR: gender, age, number of children residing at home, professional role, marital status, collected 'knowledge' score (number of correct answers), 'efficacy' measure, and four potential barriers and facilitators for WTR.
Linear regression modeling relies on finding a best-fitting regression line that minimizes the sum of squared errors of prediction. This simple and traditional model, which uses the mean values of data to calculate the regression line and is therefore sensitive to outliers, may not be an ideal method for capturing a central tendency of the data. Quantile regression, on the other hand, is a non-parametric method which consists of finding a regression line by minimizing the sum of the absolute residuals using the median rather than the mean value of the data across various quantiles of the dependent variable [39][40][41]. Thus, not only may quantile regression be expected to reflect the data in a more precise manner, it also illuminates associations between variables at different levels of the dependent variables and offers additional insights about the nature of these relationships [32].
Our analysis examined the regression coefficients across quantiles 5, 10, 25, 50 and 75 of the dependent variable (the 90th percentile was excluded due to paucity of observations).
In the third and final step, a supplementary analysis was performed to explore the variance of the WTR measure over the age distribution of participants using a broken line spline regression. To enhance the sensitivity of the results, the analysis was performed separately for females and males. Spline knots were fitted at ages 25,30,35,40,45,50,55, 60 and 64. The model was arrayed in a hierarchical design: the first block included the variables age and profession; the second block added the spline knots and employed a forward stepwise selection strategy (inclusion and exclusion criteria set at p > .05 and p < .10, respectively); and the last block added statistically significant explanatory variables identified in the previous multivariate analysis described above (professional commitment and efficacy). A P value of <.05 was considered statistically significant.
Data were analyzed using IBM Statistical Package for the Social Sciences (version 23.0; IBM corp, Armonk, NY) and STATA software (version 12.1; StataCorp, College Station, TX).

Participants
Of the 2400 questionnaires distributed to HCWs in the hospitals, a total of 852 were completed and returned, yielding a response rate of 35% (response rates ranged from 11 to 93% across hospitals). The demographic and professional characteristics of the participants are detailed in Table 1. The majority of participants were nurses, female, and married (78% of all nurses were female). The mean age was 43, and most participants had 1-3 children residing with them at the time of the survey.
Knowledge was measured through twelve multiple choice questions. An average knowledge score was calculated based on the participants' number of correct answers. Participants' responses are detailed in Table 1. The median knowledge score of participants was relatively low: MD = 6 (IQR = 5-8). When examined per profession, nurses had the highest knowledge score, namely MD = 7 (IQR = 6-8), versus MD = 6 (IQR = 4-8) for physicians and MD = 5 (IQR = 3-6.5) for paramedical staff (p < .001). Four knowledge questions had a correct response rate lower than 50%. These items related to: appropriate actions to protect oneself during an earthquake (39%); indicating the person authorized to evacuate a department following an earthquake (8%); familiarity with the treatment protocol of a common earthquake-related injury (32%); and communication procedures with external institutions or organizations (e.g. EMS services) following an earthquake (26%).
Willingness to respond to an earthquake scenario was measured through two items (personal WTR and perceived WTR of colleagues). Seventy six percent of participants reported they had "high WTR" (rating between 5 and 7 on a Likert scale). Only 72% of nurses reported high WTR, significantly lower compared with 82% among physicians and 80% among paramedical staff (p = .02).
'Concern for family wellbeing' and 'Professional commitment to care for the injured or ill' were reported by the largest number of participants as factors that could strongly affect the decision to report to work (91 and 92% of participants, respectively, rated these items between 5 and 7 on a Likert scale), followed by "Concern for house possibly damaged in the quake" (67%) and "Fear of losing place of employment" (34%). A significantly higher proportion of nurses (94%) noted "Concern for family wellbeing" as a possible barrier to WTR compared with physicians (84%) and paramedical staff (90%) (p = .001). Similar results were demonstrated for "Concern for house possibly damaged in the quake" -70% of nurses indicated this as a possible barrier, compared with only 57% among physicians and 65% among paramedical staff (p = .01). No statistically significant differences were demonstrated for the remaining two items (p > .05).

Multivariate analysis
Two multivariate regression models were employed to identify statistically significant associations between study variables and WTR. The linear regression results revealed six statistically significant predictors. Age (per year), beta (standardized regression coefficient) = .01, 95% CI: (.01, .02); efficacy score, beta = .24, 95% CI: (.16, .31); professional commitment, beta = .61, 95% CI: (.53, .68); being a nurse (compared with paramedical staff), beta = −.46, 95% CI: (−.70, −.22); being married (compared with being single), beta = −.25, 95% CI: (−.49, −.03); and 'concern for family's wellbeing' , beta = −.11, 95% CI: (−.18, −.03) had significant predictive value for WTR (see Table 2). The quantile regression results indicated statistically significant associations regarding eight predictors (see Table 2): gender, age, number of children, being a nurse, being married, efficacy score, concern for family's wellbeing, and professional commitment. Gender was statistically significant only in Q25. Age was statistically significant across Q5 to Q50, the effect size Out of bracketsmean ± SD; In bracketsmedian, minimum-maximum being strongest in Q5 and decreasing gradually thereafter. Number of children was statistically significant only in Q25. Being a nurse was statistically significant in Q10 to Q75, the strongest effect size occurring in Q10. Being married was significant only in Q10. Concern for family's wellbeing was statistically significant but only in Q50. Efficacy score and professional commitment were statistically significant across all quantiles of WTR. Figure 1 presents the statistically significant covariates (found to be statistically significant in more than one quantile) and the change in their effect size across various quantiles of WTR.

Supplementary analysisage and WTR (spline regression)
To gain in-depth understanding of the relationship between the statistically significant demographic variable 'age' and WTR, an additional analysis was conducted using a linear regression model with spline. This analysis was conducted separately for females and for males. The results are described in Table 3 and Fig. 2.
The results of the analysis conducted for the female participants revealed that the fitted spline curve was statistically significant (p < .001). The slope took a significant turn at one knot: age > 40 (beta = .05, p = .01), meaning that WTR increases significantly among females older than forty. The results for the males group did not reveal any statistically significant spline knots, as they were all excluded from the regression equation during a stepwise selection process. This indicates a rather linear relationship between age and WTR among males.

Discussion
Mitigating absenteeism is one of the most important challenges facing healthcare organizations in their struggle to improve response capabilities in the face of emergency and disaster situations. The present study sheds new light on the association of factors influencing WTR following an earthquake and discusses advanced statistical methods that might be applied for their identification and estimation.
The general high WTR rate (rating between 5 and 7 on a Likert scale) we found (76%) agrees with the rate reported in an earlier study, namely 78% WTR following an earthquake scenario; that study examined the response of hospital employees to various emergencies and disasters in the U.S. [42]. The WTR rate reported by nurses (70%) in the present study is higher than that cited in a previous Israeli study, in which only 57% of nurses indicated that they would report for duty following an earthquake [37]. This difference may be attributable to the inclusion of nursing students in the earlier study and its relatively small sample size. The WTR rate reported for nurses in the present study is also higher than the actual rate of Japanese nurses reporting for work following the Great East Japan earthquake in 2011 (48%) [15]; however, the   radiological hazard posed by that event may have affected the nurses' decision to report. A comprehensive examination of our results reveals statistically significant differences between HCWs with respect to WTR and other elements of earthquake preparedness and response depending on profession. Nurses exhibited the highest levels of knowledge and efficacy, but also the lowest levels of WTR (compared with physicians). Being a nurse (in comparison with paramedical staff ) was found to be a significant predictor of WTR, especially in the case of participants with lower levels of WTR. Knowledge was not a significant predictor of WTR. These are interesting findings that call into question theories that link knowledge and efficacy as determinants of behavioral changerepresented in this study by participants' WTR following an earthquake [43,44]. To some degree our findings also stand in contradiction with several empirical studies that report a positive association between knowledge levels, self-efficacy, and participation in education and training programs (knowledgeand efficacy-enhancing activities), on one hand, and WTR following an earthquake [37] or other menacing scenarios [19,20,26], on the other. Our results agree in part with two studies conducted in the US and one in the UK, which found that physicians were more likely to be willing to report for duty during a catastrophic event than nurses [45][46][47]. However, only one of the latter three studies directly examined earthquakes, making it difficult to compare their results with our own. One possible explanation for this disparity might be the higher incidence of nurses pointing to family wellbeing and safety of the house as possible barriers to WTR (as compared with physicians and paramedical staff). Concern for family and loved ones emerges as one of the most powerful barriers to WTR reported in the literature. This is confirmed by two literature reviews that examined factors associated with healthcare personnel WTR following disaster situations [20] and during an influenza public health emergency [19]. (Concern for damage to the house was not generally discussed in these earlier studies, as that barrier is specific to earthquake scenarios.) The fact that nurses reported such barriers more frequently may be cultural; it may reflect the fact that nursing is a profession with female dominance and that women are still perceived (even by themselves) as the ones primarily responsible for dealing with household and family matters.
A strong positive association was found between the participants' efficacy and professional commitment to care for the injured and ill and their WTR; the two variables were also identified as significant predictors of WTR in both regression models and across all WTR levels (quantiles). These factors seem to overshadow other potential barriers, such as concern for the safety of the home and fear of losing one's place of employment, neither of which was found to be statistically significant in any of the models. While professional obligation is perceived almost as a moral imperative by HCWs and has been identified as a powerful factor influencing WTR even in scenarios that pose a personal risk for the worker (e.g. infectious diseases) [19,20], it also lies at the heart of a complex conflict with personal and family obligations (e.g. child or elder care) [21]. This conflict is reflected in the present finding that 'number of The results for the males group did not reveal any statistically significant spline knots, thus, they were all excluded from the regression equation during the stepwise selection process p values set in boldface indicate statistical significance (p<0.05) children' and 'concern for family wellbeing' are statistically significant predictors of WTR (although the quantile model did not identify them as significant at the lowest levels of WTR). The supplementary analysis also revealed that young females (< 40 of age), and nurses (both females and males) were the most prone to absenteeism groups in regard to WTR following an earthquake. Here too tensions between professional and personal considerations could be at play, since in Israel's family oriented society most women in this age group have childcare obligations. Indeed, as noted above, rates of concern for family and the home were high among this group. The present findings should be taken into account when designing programs and activities aimed at improving staff preparedness for earthquake situations, both at the stage of pinpointing 'at-risk' groups among those employees whose participation is regarded as vital and in the selection of contents for seminars, drills, etc. For example, our results revealed that staff had low-to-medium levels of knowledge about their hospital's standard operating procedure following an earthquake. Areas of acute shortage of knowledge were also identified; it is recommended that these areas be targeted in future activities. Special emphasis should be placed on enhancing the employees' sense of efficacy and professional commitment in the context of duty to care following an earthquake, but at the same time they should be encouraged to voice their concerns regarding potential barriers and conflicts [23]. The healthcare system should do everything in its power to minimize these conflicts and strengthen the moral grounds for duty to care by providing employees with solutions during emergencies, such as childcare centers in the hospital compound and other family support measures [48,49]. Planning for an appropriate workforce to address these issues should be carried out well in advance. This is also noted in the World Health Organization recommendations for developing a public health response to emergencies [50].
To the best of our knowledge, this study is the first to suggest the use of advanced regression models (quantile and spline regression) to analyze data related to knowledge, attitudes, and beliefs among healthcare personnel in the context of emergency and disaster situations. Employing regression equations that perform analysis of variables by dividing them into more homogeneous sections enables assessment of complex (nonlinear) relationships and can increase the sensitivity of the results. Such a capability can be particularly valuable for organizational heads (e.g. hospital directors) as it provides a prognostic tool enabling the identification of vulnerable sub-groups of a certain population and determining the factors that may influence these groups [32,[51][52][53][54]. This can help target preparedness efforts, and improve their efficiency.
In the present analysis, being a nurse (in comparison with paramedical staff ) and age were found to be significantly associated across almost all levels (quantiles) of WTR, in contradistinction with other personal characteristics (e.g. gender and number of children), which were significantly associated only in the 25th quantile of WTR (a relatively low level of WTR) and did not even reach statistical significance in the traditional linear model. Being married (in comparison with being single) and concern for family wellbeing were singled out as statistically significant predictors in the linear model, but in the quantile model the two were significantly associated only at Q10 (low level of WTR) and Q50 (medium level of WTR), respectively. Efficacy and professional commitment to care were found to be statistically significant factors associated with all levels of WTR. An additional analysis using broken line spline regression was performed to determine the point across the participants' age distribution at which a statistically significant change occurred in WTR levels (females under/above forty). The use of quantile and broken line spline regressions provides a more comprehensive view of the relationships between variables and is especially important when the source of interest lies at the ends of the distribution of the dependent variable (characterizing, for instance, HCWs with low WTR) [51,52].
The present study demonstrates that the adoption of robust non-parametric regression analysis methods for exploring and characterizing a diverse population group such as healthcare employees provides insights that would not have been revealed had a traditional analysis been used. The implications of these findings go far beyond the present context of emergency and disaster preparedness. The use of a statistical approach that enables in-depth investigation of a particular variable of interest at all levels of the distribution makes it easier to translate research findings into practical and effective measures in the field. In this study, our ultimate purpose was to make vital information available to the heads of relevant authoritieshospital directors, department heads and emergency officersin order to support decision and policy making aimed at improving preparedness and response capabilities in a future earthquake scenario. However, this approach can also be utilized to explore other issues related to organizational culture and the behavior of healthcare personnel, whether during routine operation or in emergencies.
The present study has certain limitations. One is that it relies on self-reported data regarding intention to respond following an earthquake. Another is that the analysis does not encompass administrative and ancillary staff, who also have an important role to play in the hospital's response to earthquake events and who may present their own unique characteristics. The relatively low response rate in the current study may also pose a certain limitation, as it has implications for potential nonresponse bias, which may affect generalizability of study results. However, this was anticipated and taken under consideration while planning the data collection process and calculating sample size. Further research is needed along two fronts, first to explore and characterize the entire hospital workforce, and second to examine the proposed statistical approach in the context of other emergency and disaster scenarios, each posing its own distinct threats and risks to hospital personnel.

Conclusions
The present study shows that the use of quantile regression and broken line spline regression analysis for modeling healthcare workers' knowledge, attitudes, beliefs and intended behavior in the context of a future disaster provides a more comprehensive and sensitive view of the data than can be achieved by traditional regression analysis. These advanced analysis strategies can promote our understanding of the issues investigated and facilitate effective implementation of research findings in the field. We recommend adoption of the statistical approach presented in this study for analysis of issues related to the organizational behavior of healthcare personnel. In the context of an anticipated earthquake scenario, this can help healthcare institutions improve their preparedness efforts, ultimately enhancing response capabilities.

Additional file
Additional file 1: An English version of the survey instrument (DOC 67 kb) Abbreviations HCW: Health care workers; QR: Quantile regression; WTR: Willingness to respond