key: cord-0852400-cyyez7rr authors: Shi, Xianliang; Li, Jiangning; Huang, Anqiang; Song, Shaohua; Yang, Zaili title: Assessing the Outbreak Risk of Epidemics Using Fuzzy Evidential Reasoning date: 2021-04-17 journal: Risk Anal DOI: 10.1111/risa.13730 sha: 0ff7da16f20516bce6f7a91f1af7e6332ec64942 doc_id: 852400 cord_uid: cyyez7rr Epidemic diseases (EDs) present a significant but challenging risk endangering public health, evidenced by the outbreak of COVID‐19. Compared to other risks affecting public health such as flooding, EDs attract little attention in terms of risk assessment in the current literature. It does not well respond to the high practical demand for advanced techniques capable of tackling ED risks. To bridge this gap, an adapted fuzzy evidence reasoning method is proposed to realize the quantitative analysis of ED outbreak risk assessment (EDRA) with high uncertainty in risk data. The novelty of this article lies in (1) taking the lead to establish the outbreak risk evaluation system of epidemics covering the whole epidemic developing process, (2) combining quantitative and qualitative analysis in the fields of epidemic risk evaluation, (3) collecting substantial first‐hand data by reviewing transaction data and interviewing the frontier experts and policymakers from Chinese Centers for Disease Control and Chinese National Medical Products Administration. This work provides useful insights for the regulatory bodies to (1) understand the risk levels of different EDs in a quantitative manner and (2) the sensitivity of different EDs to the identified risk factors for their effective control. For instance, in the case study, we use real data to disclose that influenza has the highest breakout risk level in Beijing. The proposed method also provides a potential tool for evaluating the outbreak risk of COVID‐19. Epidemic diseases (EDs) can spread faster and wider than even after their outbreaks in today's community, where more and more megacities of over 10 million population appear and people communicate strengthen the capacity of ED controlling. However, the outbreak of COVID-19 causes over 1.4 million confirmed cases and over 82,000 death toll (by April 08, 2020), and the numbers are still increasing fast in 2020. The nation-wide outbreak of epidemic in December 2017 also questions the effectiveness of the current running system, with a shortage of 310,000 packets of the special drug "Duffy" merely in its capital, Beijing. Consequently, it is direly urgent to develop a new proactive risk assessment system to assess ED's outbreak risk and take corresponding prevention actions. One of the innovations of this study is to use primary data obtained through our collaboration with Chinese National Medical Products Administration (CNMPA) to develop a quantitative risk assessment method for evaluating the outbreaks of EDs. Infectious disease outbreaks can spread rapidly, causing enormous losses to individual health, national economies, and social wellbeing. Through the early detection of an infectious disease outbreak, a small outbreak can potentially be contained at the local level, thereby reducing adverse impacts. Early detection has been and remains the current narrative of infectious disease surveillance. However, it is challenging to assess the outbreak risk of EDs, given that it is affected by a variety of factors, both qualitative and quantitative, such as nature of the pathogen, people's communication style, medical treatment level, weather and natural environment. Besides, limitations to the accessibility to data on certain occasions, some factors cannot be scaled precisely but described in a fuzzy way. For example, people's hygiene awareness is often described by "very good", "good", "fair", "bad" and "very bad". Moreover, domain experts that include all the senior managers of the Ministry of Health, the State Drug Administration, Chinese Center for Disease Control and Prevention (CCDC) and the designated (by the healthcare authorities) experts in the field of public medicine healthcare argue that they have low/medium confidence to make a reliable evaluation without the aid of advanced risk tools that can address the inherent high uncertainty in data. Obviously, the assessment of ED's outbreak risk involves the integration of qualitative and quantitative factors and requires to deal with fuzzy and incomplete information, where a fuzzy evidential reasoning (ER) approach fits in well. The aim of this article is (1) to propose a fuzzy ER assessment approach to the ED's outbreak risk, (2) to apply the proposed approach to assess the outbreak risk of four kinds of major EDs in Beijing. The remainder of this article is organized as follows. Section 2 reviews the relevant literature on emergent risk assessment methods. Section 3 introduces the evidence reasoning approach and proposes a new index system for the assessment and elaborates the ER based EDRA. Section 4 conducts an empirical study related to ED risk of Beijing using the primary data from the experts of the highest healthcare authorities in China. Section 5 concludes the article with implications and limitations. There are few studies on risk assessment of ED's outbreak in megacities in the current literature. This section describes the relevant literature of risk assessment in the areas of "humanitarian logistics," "disaster operations management," and "outbreak of infectious diseases." Some studies address the network designing and inventory problems of humanitarian logistics based on cost analysis such as monetary cost, time cost, social cost, etc. For example, by considering relief time, Ahmadi, Seifi, and Tootooni (2015) propose a multidepot location-routing model to determine the locations of local depots and routing for last mile distribution after an earthquake. By minimizing facility placement, logistics, and deprivation costs, Loree & Aros-Vera (2018) propose a model to determine the location of points of distribution and inventory allocation in Post-Disaster Humanitarian Logistics. However, some others argued that humanitarian should highlight "equity" (e.g., Gutjahr & Fischer, 2018) . Once there is a disaster, relief commodities are to be supplied to the victims of a natural disaster, aid the organizations that not only have to take the total degree of demand satisfaction into account, but also require relief goods to be distributed as equally as possible among the affected population. Ideally, no region or population group should be disadvantaged (Farris, 2010; Yitzhaki & Schechtman 2013) . Substantial studies on humanitarian logistics have been published, which mainly focused on transportation and allocation of disaster relief resources, taking into account the goal of equity (e.g., Huang, Jiang, Yuan, & Zhao, 2015; Lin, Batta, Rogerson, Blatt, & Flanigan, 2011; Ransikarbum & Mason, 2016; Tomasini & Van Wassenhove, 2009) . For example, a fuzzy multicriteria, multiperiod linear programming model for relief distribution was presented in Tzeng, Cheng, & Huang (2007) . They measured the equity of distribution of relief materials to demand nodes by cumulating the least satisfaction score among demand points for each item and for each time period. Moreover, Sun, DePuy, & Evans (2014) dealt with a patient-hospital allocation model for the case of an influenza pandemic. Equity was taken into account by considering the maximal travel distance to the assigned hospital as a second objective function in a biobjective optimization model. The risk assessment of the outbreak of infectious diseases in this article also involved the principle of "equity," focusing on the infectious diseases affecting public health. Research about disaster operations management mainly focuses on "natural disaster," "war or terrorism," but rarely on the management of ED. Disaster operations management can be divided into resource/facilities management, response management, and risk management (e.g., Al-Dahash, Kulatunga, & Amaratunga, 2014; Tsadikovich, Kamble, & Elalouf , 2019) . At first, several studies aimed to develop and increase the ability and facilities to cope with disasters and to mitigate some of the effects in order to minimize the consequences of the disaster (e.g., Al-Dahash, Kulatunga, & Amaratunga, 2014; Harding, 2007) . Furthermore, the case of the Katrina hurricane was considered by Baker & Refsgaard (2007) to identify successful strategies that enabled institutions to respond effectively at an appropriate scale. Finally, the article by Carreño, Cardona, & Barbat (2007) proposed a risk management index. A group of indicators were brought together to measure risk management performance and effectiveness. Organizational and developmental capacity and institutional actions were taken to reduce vulnerability and losses in a given area. Such factors were reflected by these indicators to prepare for a crisis and to recover efficiently from disasters. This index was designed to assess risk management performance. Four public policies, include the identification of risk, risk reduction, disaster management, and governance and financial protection were constructed by the proposed risk management index. Although there are few studies on the risk management of ED, research on disaster operations management provides theoretical sup-port for the construction of the ED risk assessment model in this article. Regarding the outbreak of infectious diseases, there has been more studies in the medical field, focusing on the factors affecting the transmission of pathogens, propagation mode, and the control of infectious diseases (Lal et al., 2018) . Some found that global environmental changes, especially climate change and human exploitation of productive ecosystems (e.g., Tilman, Balzer, Hill, & Befort, 2011) had important implications for infectious disease risk (Cutler, Fooks, & Van, 2010; Murray & Daszak, 2013) . As for the propagation mode, the parasites are primarily spread through contaminated drinking or recreational water; however, infection in humans may arise through contaminated food, contact with animals, especially livestock or infected individuals (Cacciò, Thompson, McLauchlin, & Smith, 2005; Savioli, Smith, & Thompson, 2006) . With regard to diseases control, a few studies identified that the environmental and sociodemographic exposures can help develop disease control priorities in high-risk areas (Lal, Cornish, Fearnley, Glass, & Kirk, 2015; Yoder & Beach, 2010) . In addition to these theoretical studies, practitioners have also made many outstanding contributions to the investigation of the outbreak of infectious diseases. Innovative governance structures have been established to promote early detection. Disease surveillance networks have been formed, such as the World Health Organization (WHO), Global Outbreak Alert and Response Network (GOARN), combining human and technical resources around the world to rapidly identify, confirm, and respond to outbreaks. Cross-border regional disease surveillance networks have been established across the globe, connecting epidemiologists, scientists, ministry officials, health workers, border officers, and community members to engage in activities, such as training, capacity-building, and multidisciplinary research. Agreements have been instituted, setting legal mandates around surveillance activities, such as the IHR (2005), which call for all the WHO member states to build, improve, and strengthen their capacity to prevent, detect, and respond to infectious diseases outbreaks that can have global spread (Steele, Orefuwa, & Dickmann, 2016) . To the best of our knowledge, the risk assessment of ED's outbreak is still an open task, owing to the high complexity and requirements of integrating qualitative and quantitative factors and dealing with uncertainty in data. Fuzzy ER presents an effective solution to dealing with the above requirements. Therefore, this article aims to use a fuzzy ER approach to develop a new risk assessment model for evaluating the ED's outbreak in megacities. As far as the risk assessment of ED's outbreak is concerned, some factors cannot be precisely measured but often can be described by subjective judgements using linguistic terms. Such linguistic terms are then modeled by with fuzzy membership functions to address their discretization. Linear (e.g., triangular and trapezoidal) fuzzy membership functions are widely used to describe the linguistics variables because of their simplification (Wang, 1997) . A triangular fuzzy number is a fuzzy set with three parameters (a, b, c) , where a is the membership function's left intercept (lower boundary) with a grade of 0 and c is the membership function's right intercept (upper boundary) with a grade equal to 0, while b is the value of having the highest membership with a grade equal to 1 (Pam, Li, Wall, Yang, & Wang, 2013) . Fuzzy logic provides experts with a wide scope instead of a precise point to define the risk grade at which a key risk indicator (KRI) may contribute to a specific epidemic disease. Fig. 1 visually depicts the proposed methodology composed of three major steps. The first step is to establish a fuzzy link-based KRI hierarchy, which can offer a generic platform for carrying out the risk assessment of a particular place under the outbreak of any specific ED. The developed generic framework can be appropriately modified with reference to the disease's specifications, involving less or more KRIs. The second step is to analyze the risk level of a targeted area facing the impact of the specific disease. The evaluation of KRIs is a dynamic process depending on different EDs. All the estimations of the identified KRIs can then be aggregated using an ER approach to obtain an assessment result. It will then be measured by a risk threshold indicating whether the criticality of the outbreak under an ED is acceptable. The risk index value is obtained through aggregating the marginal risk estimation of each KRI. Otherwise, all the acceptable risk assessment results obtained from different impact factors and from different conditions will be aggregated to assess the overall risk level of an ED. The overall risk level can then be used as a benchmark to improve the level of prevention and monitoring of infectious diseases in CCDC, as well as a theoretical reference for the government to regulate and control drug reserve (Comfort, 2002) . Details about ER can be found in the book (Lee & Yang, 2018) providing a comprehensive review. This section extracts the KRIs for the risk assessment from the four stages of the ED's outbreak process, including prevention stage, propagation stage, resistance stage, and post-control stage. We obtained the final indicator system in Table I by interviewing 10 experts. Out of these experts, six are from public healthcare authorities, three from the Chinese Food and Drug Administration, three from the Chinese Center for Disease Control and Prevention (CCDC). The other four are the domain frontline experts from Chinese PLA General Hospital and People's Hospital of Peking University, two top-level medical institutions in China. We first drafted a set of indicators referring to the existing studies and some national regulations like National Regulations on Emergency Response to Public Health Emergencies. Then, we presented the draft to the experts and asked them to make adjustment to the draft including adding new ones, alter the existing or drop the unsuited. The feedbacks from the experts are aggregated and sent back to the experts again for second round improvement. After two rounds, we finalized the indicator system by dropping those rejected by more than six experts. Appendix 2 presents the mapping from indicators in Table I to corresponding references. It is notable that the KRIs in Table I , the value of R denotes the capacity of inhibiting diseases, thus a lower value of R implies higher outbreak risk. The real meaning of risk in practice is varied and specific. Therefore, it is necessary to measure the risk of using KRI level descriptions of different units in common space. Triangular fuzzy membership functions developed on a [0−1] utility domain have been widely used to define different risk grades/descriptors (Yang, Bonsall, & Wang, 2009 ), which include four to five linguistic terms ( Table I) . The typical membership functions for four to five linguistic variables are defined and characterized (Table II) . There is some flexibility in the definition of membership functions to accommodate different risk situations. However, it is worth noting that any change in the defined member functions requires very cautious reasons from domain experts (Yang, Ng, & Wang, 2014) . New definitions of memberships need to be verified through empirical tests before being applied in practice. Once the fuzzy memberships are defined, the grades of the lowest level KRIs will need to be transformed and presented by the ones of the top level KRI so as to aggregate the risk estimations of all the KRIs on the same plate. Fuzzy similarity calculation can be used to model the transformation process, as follows (Yang et al., 2014) : (1) where x covers the domain [0,1] of the fuzzy memberships μ A i (x) and μ A j (x) of the lower level KRI grades A i and the upper level KRI grades A j , i, j ∈ (4 or 5); a i j presents the similarity degree to which A i belongs to A j . To keep the completeness of risk estimation a i j is normalized into b ij , as follows: Therefore, all KRI scores at the lowest level can be converted and presented through the top KRI scores (Good, Average, Fair, Poor). An example of transforming KRI "R-P1-I2" to "R" can be found in Fig. 2 . The belief degrees assigned against each link in the figure is calculated by Equations (1) and (2). In the face of different infectious diseases, the risk of a particular region will vary. Therefore, disease is identified in a region. When identifying specific diseases in a region, the risk level can be assessed by considering two types of estimates. One of them is to estimate the risk of all KRIs at the lowest level (level I), and the other is to assign weights to all relevant KRIs. For a specific disease, not all KRIs within the hierarchy have the same importance/weight. Some of them may be less (or even less) influential than others. The weights of all KRIs are dynamic and affected by identified diseases, so a set of KRI weights is needed for each region against a specific disease. Here Analytic Hierarchy Process (AHP) (Saaty, 1998) can be repeatedly used to assign the weights w, of KRIs (Yang, Ng, & Wang, 2013) . After analyzing the weight allocation, the KRI risk assessment for each region-disease pair becomes simple and clear. To present assessors' uncertainty when estimating KRI I i , (4 or 5), a set of probabilities p i is introduced correspondingly. For instance, p i of the KRI "R-P1-I1" can be estimated as {p 1 Full, p 2 Partial, p 3 Poor, p 4 Barely, p 5 No}. Such estimation can be transformed through β i j in Equation (2) and the pathways in Fig. 2 and expressed based on the four linguistics grades of R using Equation (3): where p i ( n i=1 p i ≤ 1) indicates the probabilities assigned to the grades in the lower-level indicator; β i j means the normalized similarities between the grades of lower and upper levels (see Equation (2)); β i j represents the transformed probabilities assigned to the grades in the upper-level indicator. Assume that w i and w j indicate the relative weights associated with the lower and upper-level indicators, respectively. Then, where θ i represents the importance of the ith indicator at the lowest level in the process of synthesizing all the relevant lower-level indicators to their common (the upper level jth) upper-level indicator. Equations. (3) and (4) can be repeatedly used to transform the risk estimation and importance of each indicator at the lowest level to the top-level R. We note that the sum of all the θ i associated with the "R" is equal to one. As the risk estimations and weights of all the KRIs at the lowest level are transformed to their counterparts θ k and β k j at the "R" with respect to one region-disease pair (k = (1, 2, …, L), where L means the total number of all the KRIs at the lowest level. To capture the nonlinear relationship between all the KRIs at the lowest level, the ER approach (Yang & Xu, 2002a ) is used to combine all β k j (transformed from each KRI at the lowest level) and generate a risk estimation of the region-disease pair. First, it is required to transform β k j into basic probability masses using Equations (5)-(8) (Yang & Xu, 2002a) : where each m k j is a degree to which each KDI supports the final synthesized estimation D (inhibition level in Table I D . Then, the overall combined belief degree in D j is generated as follows: {D j} : where β j indicates the normalized belief degree assigned to D j in the final synthesized estimation D and β D represents the normalized remaining belief degree unassigned to any D j . To facilitate the calculations involved in Equations (5)-(15), IDS Software has been developed through a user-friendly interface by Yang and Xu (2002b) . The risk estimation result for the region-disease pair can then be expressed as {β 1 Highest, β 2 High, β 3 Low, β 4 Lowest, β D Un-known}. To measure the acceptance of such risk result, D j needs to be given utility values for a crisp risk result R C and β D require to be assigned back to β 1 and β 4 for the possible best risk R B and the possible worst risk R W . Consequently, where each U j (j = 1, 2, 3, 4) represents the utility values of D j , which can be calculated using a centroid defuzzification method as the set of {0.1, 0.4, 0.6, 0.9} from Table II . Such defuzzied utility values can be used as the criteria to define the three risk levels (Yang, 2001) , as follows: When R C < 0.4, the risk level of the regiondisease pair is unacceptable at a normal situation. Additional risk control measures should be developed and adopted to the KRIs with the highest risk contributions. Although this strategy seems to be debatable on the perspective of cost-benefit analysis, it is actually reasonable given that cost is not the first concern in a pandemic. When R C ≥ 4, the risk level is acceptable. When the risk levels of all the region-disease pairs related to a terminal become acceptable, the risk level of the terminal can be obtained through synthesizing all the risk levels of the involved region-disease pairs for risk improvement purposes. β j and β D obtained from Equations (14) and (15) are based on a specific region-disease pair by an expert. In practice, the risk level of a certain area with respect to all the possible diseases and disease response capacity in specific regions are likely to be assessed by multiple experts. Therefore, the risk estimations from multiple experts and all the regiondisease pairs are synthesized using the ER approach. If there are different risk results from multiple experts with respect to a specific pair, then β k j in Equation (5) means the estimation from the kth expert. The weight θ k indicates the role/importance that the expert plays compared to the other experts. If all the region-disease pairs associated with a particular area under different disease are combined to calculate the risk level of the facility, then β k j in Equation (5) represents the estimation from the kth region-disease pair. The weight θ k is determined by taking into account both likelihood and consequence through a simple risk matrix approach of each disease. Similarly, if the risk estimations of all the facilities are synthesized to calculate the risk level of a certain area, then β k j in Equation (5) represents the estimation from the kth type of disease under investigation. The corresponding R C can be used a benchmark to monitor the risk control of an area under the outbreak of diseases in a longitude study, as well as to priories all the region-disease pairs to identify the risk vulnerability for optimal risk improvement. In order to further analyze and verify the key factors affecting the outbreak of EDs, such as "which factor should attract the most attention by healthcare authorities," we will conduct a sensitivity analysis. Decision makers are often interested in examining how much changes in attribute weights and alternative performance on an attribute can affect the overall ranking of alternatives. On top of score ranges, there are three other types of sensitivity analyses in ER and its associated computing software IDS: Change Weight, Change Input Data, and Trade Off Analysis. In this article, in order to obtain the key means to suppress the outbreak of infectious diseases, the method of change weight is selected to conduct sensitivity analysis. At present, the scope of drug and medical device reserve in China is mainly general accidents and disasters. However, in recent years, with the change of natural environment and ecosystem, as well as the increase of international exchanges of personnel and materials, the incidence and scale of new infectious diseases and infectious diseases in overseas epidemic areas have been continuously increasing. The incidence of unconventional epidemics and new infectious diseases is increasing, and epidemics and diseases beyond the control of the use of existing drugs and medical devices occur. At the same time, due to the shortage of supply, there is no domestic production capacity or insufficient production capacity, in the emergency response process, the event handlers often require the use of imported products on the premise of having domestic products. The inaccurate forecasting of epidemic trend and scale leads to the inability of reserve varieties and quantities to reflect emergency needs. Through sensitivity analysis, the key factors affecting the suppression of infectious disease outbreaks in the region can be obtained, which will provide useful insights for CCDC and other relevant departments in disease prevention, response and control decision-making. As the capital of China, Beijing is the political, economic, and cultural center of China, which has a dense population and the number of permanent residents ranks 1st in the country. Taking Beijing's public health emergency drug reserve as the object of empirical research has important theoretical value and practical significance for emergency measures of public health emergencies in megacities. The aforementioned expert panel recommend four diseases worth careful investigation in light of their background knowledge, including influenza, human infection with highly pathogenic avian influenza (Avian Influenza), hand, foot and mouth disease (HFMD), and tuberculosis. Subsequently, we take Beijing as an example to explore the risk of the outbreak of the four infectious diseases. Firstly, based on the risk assessment system in Section 3, a questionnaire survey was conducted with the panel of experts. On the one hand, the scientific nature of the assessment system will be tested by the survey data, on the other hand, the data will be substituted for evidence reasoning method to make the risk assessment. Finally, we will conduct a sensitivity analysis on the key indicators affecting the risk of infectious disease outbreak. The indicators used in this study are evaluated by a subset of perceptual items in the outbreak of a disaster. First, we adopted indicators that had been validated by previous studies. Second, when indicators had not been well documented in the literature, we develop new indicators based on our observations and face-to-face interviews with the senior managers in the Ministry of Health, CNMPA, and CCDC. The survey instruments included questions on the basic information of respondents and evaluation of four main infectious diseases in four stages. In the prevention stage, we evaluate the risk in terms of effective vaccines (Lopes, Okino, Fernando, Pavani, & Montassier, 2018; Lucero et al., 2019; Wiederman, 2017) , personal hygiene habits (Pehlivan et al., 2011; Zivich, Gancz, & Aiello, 2018) , and pathogens (Xue et al., 2019; Yanagihara, Matsumoto, Aoki, & Sato, 2019) ; in the propagation stage, in terms of transmission speed (Ding & Wang, 2019) , human beings' resistance and mortality rate (Emam, Livernois, Paibomesai, Atalla, & Mallard, 2019; Rao & Ayres, 2017) ; in the resistance stage, in terms of diagnosis and control of diseases, storage of drugs/vaccines, effectiveness of the treatment, and possibility of death (Arji et al., 2019; Odedra, Lalloo, Kennedy, Llewellyn, & McCarthy, 2019; Vaughn, Salas-Wright, Cohen, & Holzer, 2019) ; and in the post-control stage, in terms of the level of drug stocks available, the ability to regulate the supply of drugs and targeted science popularization work (Spencer et al., 2019; Tong et al., 2016) . On the four stages affecting the outbreak of infectious diseases, all indicators were computed on a five-point Likert scale and are listed in Appendix 1, with "1" indicating "strongly disagree" and "5" indicating "strongly agree." Through reliable research collaboration, we managed to survey CCDC, Chinese Drug Administration and the Ministry of Health using anonymous questionnaire in the period from 1st Oct to 31st Oct, 2018. Four weeks were given to each respondent, and every two weeks was considered as one period. We then sent reminding emails to the respondents when the first period passed. By the end of the fourth week, 76 respondents had responded to the survey. On the basis of a comprehensive check of the returned answers, 55 questionnaires were valid, translating to a response rate of 72.3%. The samples cover all the senior managers of the Ministry of health, the State Drug Administration, the CCDC, and the designated (by the healthcare authorities) experts in the field of public healthcare. The good response rate and comprehensive coverages of the stakeholders ensure the representativeness of the sample data. Profiles of the respondents and the organizations are summarized in Table III . According to Armstrong & Overton (1977) , nonresponse bias should be tested based on the differences of measured items between early and late responses. In this study, 19 answers were collected at the end of the first period and 36 questionnaires were returned during the second period. By comparing the late responses obtained after a reminder email with early responses on the measured items and demographic variables, the results of t-tests indicate no significant differences between the responses in the two periods, suggesting that nonresponse bias was not a concern in this study. At the same time, we test the reliability of the data by using the Cronbach's alpha, and the result of this test is 0.786, showing a good reliability of the dataset. In addition, to ensure the suitability of collected data for the factor analysis, a Kaiser-Meyer-Olkin (KMO) measure of sampling adequacy and Bartlett's test were conducted (Eltayeb, Zailani, & Ramayah, 2011) . The overall result of the KMO and Bartlett's test for the dataset is 0.873 and the minimum value of a single variable reached 0.753, which is for sphericity, showing a satisfactory result when p < 0.000. The results indicate that the correlation between variables were strong and the data set structure was robust. Principal components analysis (PCA) is employed by using SPSS on 14 measured variables with VARIMAX rotation to extract the crucial factors. An eigenvalue greater than 1 was used to determine the number of factors (Churchill, 1979) . For simplicity, only loadings of 0.50 or higher were extracted (Hair, Anderson, Tatham, & Black, 1998) . Therefore, these 14 variables were loaded on four factors (shown as Table IV ). The cumulative variance explained by the four factors was 0.857 and the alpha value of each factor group was higher than the suggested threshold of 0.70 (Skipper & Hanna, 2009 ), which indicates a good construct validity. Moreover, it shows that the factor model produced by factor analysis was largely consistent with the theoretical one. Next, the fuzzy link base hierarchy in Table I is directly employed in this section (Step 1). The risk assessment of the four diseases starts with the analysis of the KRI weights with respect to each region-disease pair and the risk estimation of each KRI at the lowest level ( Step 2). AHP is first used to assign the four parameter (P1-P4) weights in the influenza analysis in Table V . In a similar way, the relative weights of all other KRIs at the lowest level are assigned (Table V) . Next, the probabilistic risk estimation of each KRI using the defined linguistics grades is evaluated in the context of influenza by the panel of experts mentioned in Section 3-Step 1 (Table VI) . We adopt a three-round Delphi algorithm to narrow the range of experts' idea, and then average them into a unified one. Using Equations (1)-(3), such probabilistic risk estimations can be transformed and expressed by the four linguistic grades of the R at the top level (Table VI) . Based on Equations (4)-(15), all the transformed risk estima-tions of the KRIs can be synthesized using the ER approach. IDS software package is used to facilitate the synthesis process with the final combined risk estimation result of the influenza obtained as R(influenza) = {0.3533 Highest; 0.3432 High; 0.2653 Low; 0.0382 Lowest}. The above result can be interpreted as that, given the risk input of each KRI in Table VI , the risk level of the influenza pair is estimated as the highest with 35.33% degree of belief, high with 34.32% degree of belief, low with 26.53% belief degree, the lowest with 3.82% degree of belief. This result is visually presented in Fig. 3 . The above results obviously indicate that, to a rather large extent, the KRIs have been assessed as the highest or the second highest grade. For example, R-P1-I1 has been assessed as "Full" with a belief of 30%; R-P2-I1 has been evaluated to a significantly larger extent as "Good" with 30%. Since the risk level of a region-disease pair is determined by the risk performance of each basic indicator, the top risk level should be evaluated as "Highest" to a large extent. This is harmonious with the results obtained above as the risk level of influenza has been assessed as "Highest" to the extent of 35.33%. The acceptance of the risk level could be measured by Rc using Equation (16) and the result is as follows: R C = 0.20402 < 0.4 implies that the risk level of influenza is worth special vigilance. Similarly, the risk levels and their crisp utility values of the other three diseases can be obtained (Table VII) . The results indicate that the risk levels of the other three infectious diseases are needed heightened vigilance, while the smaller the inhibitory effect on influenza, the greater the risk of outbreak (can see in Fig. 4 ). Fig. 5 shows the scoring of four infectious diseases in four stages, which indicates that in the prevention stage and resistance stage, the differences of four infectious diseases are small, while in the propagation and post-control stages, the scoring differences are large. In order to further analyze the key factors to inhibit the risk of infectious diseases outbreak, the sensitivity analysis of the indicators is made in the ensuing section. We invited the aforementioned 10 experts to provide the ranges of the weights of KRIs, considering each of them has at least 8 years of serving the governmental departments or medical institutes on the infectious disease. We formulate and adopt the overlapping range in the sensitivity analysis: the range of the weight of R-P2 (indicators in the propagation stage) is [0.07, 0.61], while R-P4 (indicators in the post-control stage) [0.2, 0.8]. The results of sensitivity analysis are shown in Figs. 6 and 7, from which the following implications can be drawn: First, as shown in Fig. 6 , the weights of indicators in the transmission stage impose a significant impact on the results. Along with the weight increasing from 0.07 to 0.61, tuberculosis and avian influenza transpose their rankings in terms of inhibition level, which implies the ranking of outbreak risk change reversely; besides, influenza and HFMD drop a lot though their rankings remain unchanged. This indicates that the weight in the transmission stage should be determined with much cautiousness. To achieve this goal, integration of expert opinions, like in this study, could be a good option. Second, as depicted in Fig. 7 , the weight of R-P4 (indicators in the post-control stage) has no significant impact on the risk estimation of influenza, HFMD and tuberculosis, but relatively bigger influence on avian influenza. With the weight increasing from 0.2 to 0.8, of the first three diseases, the rankings keep unchanged and the inhibition levels gain only slight increase. In terms of avian influenza, the ranking keep unchanged, though the inhibition level gains bigger increase. This implies that the empirical result is robust on the weight of R-P4. The article reviews the current status of risk control process to identify the need of developing a quantitative risk assessment method. Following the development of a novel risk assessment method, a case study and sensitivity analysis are undertaken to validate the results. The unique contributions of this study comprise (1) constructing a new hierarchical system of indicators for outbreak risk evalu- ation, based on previous studies and the first-hand practical experience of the Chinese authoritative infectious disease control institutions like Ministry of Health, CNMPA and CCDC; (2) proposing a general methodology of outbreak risk evaluation based on fuzzy evidential reasoning, which could be easily extended to other areas or diseases. The proposed system of indicators thoroughly covers key elements in successive stages of fighting a disease including prevention stage, transmission stage and resistance stage, involving human, disease and preventive measures. By adopting fuzzy evidential reasoning, the proposed EDRA methodology is enabled to integrate the valuable expertise and quantitative analysis for more successful disease outbreak risk evaluation. The development of a quantitative EDRA methodology capable of realizing the integrity of risk identification, analysis, control and security improvement will enjoy substantial benefits, including a consistent regime which addresses all the aspects of infectious disease in a certain area through an integrated and quantitative manner; a rational basis that prevention control and regulatory requirements are in proportion to the severity of the risks; a proactive approach, enabling threats that have not yet given rise to tragedies to be properly considered. To further appreciate the benefits, the applications of the developed methodology in more real cases will be carried out. The method will be tailored to deal with specific requirements of CCDC or other healthcare authorities in a wide context (i.e., Other public health events in a certain area) in order to reinforce its validity and flexibility. To conclude, this study offers substantial value to scholars, planners and policymakers in the design, planning and evaluation of effective risk control measures and strategies, thus reducing the outbreak risk of infectious disease in megcities. Moreover, the developed model is designed to deal with EDs in a generic manner. It therefore shows potential significance for tackling the risk analysis of COVID-19. Within this context, the model can be used at different scales, referring to a community, a city, and a regional/country. The uniqueness of this model lies in the fact that it covers all the stages of the pandemic development process, which enables it to evaluate the risk of the countries suffering from COIVD-19, as well as, the outbreak risk of COIVD-19 in the other countries where the large-scale outbreak of COVID-19 has not occurred yet. Furthermore, while this model can aid risk analysts and policymakers to make rational decision based on the final overall risk evaluation score, the score against each specific KPI can also provide useful insights into the response system against COIVD-19 and facilitates the exploration of the deficiency of the system. It offers an actionable method to enhance the effects of fighting the pandemic. Furthermore, although the world has put much effort to fight COIVD-19, there is a scarcity of interdisciplinary studies on how COIVD-19 reacts with other social systems. For example, how the logistics system can adapt to large-scale epidemic outbreak. Hence, by focusing on COVID-19, the proposed model can be tailored to develop effective logistics strategies and planning in tackling the outbreak risk caused by sudden, large-scale epidemic outbreak and to enhance the resilience of emergency logistics chains (e.g., medical and food). It can strive to achieve two interactive analysis results to provide insightful solutions to COVID-19: outbreak risk assessment of COVID-19 and logistics adaptive responses. It can aid to investigate and define the key factors influencing the outbreak risk (e.g., both likelihood and consequence) of COVID-19 in a country/society and to develop a novel dynamic multiobjective optimization model for logistics planning to support societies with respect to their dynamic COVID-19 risk levels. Although this study contributes to both academics and practice, there are still some limitations. For example, by interviews and questionnaires, the indicators in this study are selected by surveying infectious disease institutions located in Beijing, subsequently more indicators should be considered when it comes to other areas. This study is supported by the National Natural Science Foundation of China (Grant number 71831001), the international cooperation project with Liverpool John Moores University "Emergency logistics response to outbreak risk of COVID-19: Phase 1 risk analysis" and Beijing Logistics Informatics Research Base. The purpose of this questionnaire is to assess the risk factors for EDs and to select four representative infectious diseases, including plague, influenza, SARS, avian influenza, and hand, foot and mouth disease. Each question represents a type of indicator and has five options, and the more left option represents the higher the score value. Prevention stage • There is the effective corresponding vaccine (the higher the score, the greater the impact on inhibiting the epidemic) • Good personal hygiene awareness (washing hands, washing clothes, etc.) has a significant preventive effect on the spread of infectious dis-eases (the higher the score, the greater the impact on suppressing the epidemic) • The pathogenicity of the pathogen itself (the higher the score, the greater the impact on the epidemic) • Pathway of pathogen transmission and speed of transmission (the higher the score, the greater the impact on the epidemic) • The level of immunity/resistance of the population to the disease (the higher the score, the smaller the impact on the epidemic) • The possibility of death and serious sequelae after infection in healthy people (the higher the score, the greater the harm of the disease) • Have relevant stocks of drugs/vaccine for the disease (the higher the score, the stronger the ability to control the disease) • The possibility of death and serious sequelae after infection in healthy people (the higher the score, the weaker the ability to control the disease) • Have good diagnosis and control ability for the disease (the higher the score, the stronger the ability to control the disease) • The effectiveness of the treatment (the higher the score, the stronger the ability to control the disease) • There is a comprehensive prevention and control plan for the disease (the higher the score, the stronger the ability to control the disease) • The level of drug stocks available for this type of infectious disease (the higher the score, the stronger the ability to control the disease) • The ability of the government to regulate the supply of drugs for this type of infectious disease (mainly including market scheduling, capacity control, availability of reserves) (the higher the score, the stronger the ability to control the disease) • Targeted science popularization work for this type of infectious disease (the higher the score, the greater the impact of the measure on the suppression of the epidemic) A humanitarian logistics model for disaster relief operation considering network failure and standard relief time: A case study on San Francisco district Effect of hand hygiene on infectious disease risk in the community setting: a meta-analysis Evaluation of the system of disaster management resulting from war operations and terrorism in Iraq Fuzzy logic approach for infectious disease diagnosis: A methodical evaluation, literature and classification Estimating nonresponse bias in mail surveys Institutional development and scale matching in disaster response management Unravelling cryptosporidium and giardia epidemiology A disaster risk management performance index A paradigm for developing better measures of marketing constructs Managing intergovernmental responses to terrorism and other extreme events Public health threat of new, reemerging, and neglected zoonoses in the industrialized world Control strategy design for the anti-HBV mathematical model Green supply chain initiatives among certified companies in Malaysia and environmental sustainability: Investigating the outcomes. Resources Conservation and Recycling Genetic and epigenetic regulation of immune response and resistance to infectious diseases in domestic ruminants The Gini index and measures of inequality Nuanced risk assessment for emerging infectious diseases Equity and deprivation costs in humanitarian logistics Multivariate data analysis with reading Man-made disaster and development: The case of Iraq Modeling multiple humanitarian objectives in emergency response to largescale disasters Cryptosporidiosis: A disease of tropical and remote areas in Australia A Bayesian spatio-temporal framework to identify outbreaks and examine environmental and social risk factors for infectious diseases monitored by routine surveillance Multi-criteria decision making in maritime studies and logistics: Applications and cases A logistics model for emergency supply of critical items in the aftermath of a disaster Inactivated infectious bronchitis virus vaccine encapsulated in chitosan nanoparticles induces mucosal immune responses and effective protection against challenge Points of distribution location and inventory management model for Post-Disaster Humanitarian Logistics Plant-based vaccine candidate against infectious bursal disease: An alternative to inactivated vaccines for breeder hens Law of the People's Republic of China on the prevention and control of infectious diseases. National Health of Commission of The People's Republic of China Safety and effectiveness of apheresis in the treatment of infectious diseases: A systematic review A subjective approach for ballast water risk estimation The comparison of socio-economic conditions and personal hygiene habits of neuro-Behçet's disease and multiple sclerosis patients Goal programming-based post-disaster decision making for integrated relief distribution and early-stage network restoration Resistance and tolerance defenses in cancer: Lessons from infectious diseases The analytic hierarchy process Giardia and Cryptosporidium join the 'neglected diseases initiative' Minimizing supply chain disruption risk through enhanced flexibility Systemic peptide mediated delivery of an siRNA targeting α-syn in the CNS ameliorates the neurodegenerative process in a transgenic model of Lewy body disease Drivers of earlier infectious disease outbreak detection: A systematic literature review Multi-objective optimization models for patient allocation during a pandemic influenza outbreak Regulations on emergency response to public health emergencies Food demand and the sustainable intensification of agriculture Humanitarian logistics Beauvericin counteracted multi-drug resistant candida albicans by blocking ABC transporters Controlled information spread for population preparedness in disaster operations management Multi-objective optimal planning for designing relief delivery systems Recent trends in cooperativeness among participants in the national survey of drug use and health Assessing public health risk due to extremely high temperature events: climate and social parameters A subjective methodology for safety analysis of safety requirements specifications How effective are influenza vaccines International Health Regulations Aquaporin-4 deficiency reduces TGF-β1 in mouse midbrains and exacerbates pathology in experimental Parkinson's disease Nationwide surveillance of bacterial respiratory pathogens conducted by the surveillance committee of Japanese Society of Chemotherapy, the Japanese Association for Infectious Diseases, and the Japanese Society for clinical microbiology in 2014: General view of the pathogens' antibacterial susceptibility Rule and utility based evidential reasoning approach for multi-attribute decision analysis under uncertainties Nonlinear information aggregation via evidential reasoning in multi-attribute decision analysis under uncertainty On the evidential reasoning algorithm for multiple attribute decision analysis under uncertainty Use of fuzzy evidential reasoning in maritime security assessment A new risk quantification approach in port facility security assessment. Transportation Research Part A: Policy and Practice Prioritizing security vulnerabilities in ports More than a dozen alternative ways of spelling Gini Cryptosporidium surveillance and risk factors in the United States Effect of hand hygiene on infectious diseases in the office workplace: A systematic review (2006) R-P4-I4 Targeted science popularization work Aiello, Coulborn, & Perez (2008)