key: cord-1016529-wzhgr09l
authors: Seiti, Hamidreza; Makui, Ahmad; Hafezalkotob, Ashkan; Khalaj, Mehran; Hameed, Ibrahim A.
title: R.Graph: a New Risk-based Causal Reasoning and Its Application to COVID-19 Risk Analysis
date: 2022-01-12
journal: Process Saf Environ Prot
DOI: 10.1016/j.psep.2022.01.010
sha: f5e2c874f565741a195fa0714497ce6c5e052aa9
doc_id: 1016529
cord_uid: wzhgr09l

Various unexpected, low-probability events can have short or long-term effects on organizations and the global economy. Hence there is a need for appropriate risk management practices within organizations to increase their readiness and resiliency, especially if an event may lead to a series of irreversible consequences. One of the main aspects of risk management is to analyze the levels of change and risk in critical variables which the organization's survival depends on. In these cases, an awareness of risks provides a practical plan for organizational managers to reduce/avoid them. Various risk analysis methods aim at analyzing the interactions of multiple risk factors within a specific problem. This paper develops a new method of variability and risk analysis, termed R.Graph, to examine the effects of a chain of possible risk factors on multiple variables. Additionally, different configurations of risk analysis are modeled, including acceptable risk, analysis of maximum and minimum risks, factor importance, and sensitivity analysis. This new method's effectiveness is evaluated via a practical analysis of the economic consequences of new Coronavirus in the electricity industry.

COVID-19 has rapidly led to unprecedented health, economic, and political crises across the globe (Djalante et al., 2020) , becoming a major risk factor for many organizations. Many companies have indicated that the impact of Coronavirus COVID-19 is, or will be, a significant source of uncertainty. According to the OECD Economic Outlook Interim Report (March 2020) , annual global GDP growth projections for 2020 have dropped by half a percentage point to 2.4%, largely due to the coronavirus outbreak. However, a longer-lasting and more intense coronavirus outbreak could even slow global growth to 1.5% (OECD, 2020) . The coronavirus pandemic may eventually fade, as the Ebola, Zika, and SARS viruses have in recent history. However, even if it does, the next as-yet-unknown devastating outbreak is not so much a matter of "if" but "when." How, then, should organizations, societies and governments prepare for similar future possibilities? Even with the knowledge that such events may occur, there is a significant difference between mere awareness and the actual experience; evaluating the correct response to these potential future situations is vital for any company, institution, or country which wishes to remain competitive in the current globalized world ("Embed. Resil. A Guid. to Bus. Implic. COVID-19," 2020) . One solution to increase resiliency is risk-informed development. Risk can be defined as the probability of a certain deviation in achieving a goal, which can be determined by modeling the interacting risk factors (Seiti and Hafezalkotob, 2020a) . Risk factors can be considered as effective factors which may cause variations in predictions (Seiti et al., 2020b) . Risk analysis constitutes a family of approaches to aid top managers in assessing all potential impacts through considering the criticality of various risk factors within systemic procedures. Two significant aspects of specific risk are the probability of a risk occurring (known as stochastic uncertainty) and the variability and changes in predictive consequences (aleatoric uncertainty) (Hillson, 2016) . Therefore, the risk of interactive risk factors on each other can be assumed in three different ways: 1) the effect of a specific risk factor on the probability of an event; 2) the effect of the risk factor on the severity of an event; and 3) the combined effect of the risk factor on the probability and severity of a consequence. For instance, let us assume the probability of a person being imprisoned for a certain period of time. For this, finding a new document in court or a further witness, as a factor, can have an effect on reducing or increase in the length of imprisonment, an increase/decrease in the probability of incarceration, as well as the simultaneous effect on the probability and duration of confinement. One of the most effective methods for modeling and determining the probability and variability in a system is to break it down into smaller components through causality analysis models. These components can then be used to identify a set of factors that affect each other, examining a chain of causes and effects to make better predictions of risk and variability. The motivation for conducting this study is based on the fact that, in some risk analysis problems, the decisionmaker is interested in knowing various change rates. These can relate to the occurrence of events that: have not been previously predicted; have already occurred; or have not yet occurred but may alter relevant predictions and variables. In many cases, either events happen or it is obvious to the decision-maker that the event will happen eventually. In such cases, the models can be assumed to be definite (Hillson, 2016) , thus only variability needs to be considered. Since, in the real world, many variables are continuous, it is necessary to develop a suitable model that can find the percentage of change predicted in the desired variables. The challenge for some of these issues is that the probability of certain events occurring may be very low, meaning that there is little data available for risk analysis. In these cases, therefore, instead of empirical or statistical data, the knowledge of experts in related fields can be employed during the process of evaluating the relationships between the model components (Panula-Ontto, 2016) . There are several causal methods in the literature for risk analysis, whose various characteristics and disadvantages are investigated in Section 2. The motivation of this study is then to overcome J o u r n a l P r e -p r o o f these disadvantages and to find an appropriate solution; a causal mathematical model is proposed, called R.Graph. This model has the ability to estimate variability and risk while considering different scenarios in a causal chain of various factors. It uses data gathered from experts, while its outputs are easily interpretable and explicable by decision-makers. The current paper is arranged in the following order. Section 2 surveys existing literature on causality models and methods of risk analysis, while Section 3 discusses the preliminaries of risk and aggregation operators. Section 4 sets out the R.Graph model. A relevant case study is presented in Section 5, and Section 6 summarizes the findings, and proposes insights for future research.

In this section, the literature survey on causality models and their applications in risk analysis is provided in Section 2.1; the limitations of these methods, and the novel contributions of this paper are then presented in Section 2.2.

There are various methods in the existing literature to analyze these causal chains, which can be divided into two categories: deterministic models including structural models (Bashir et al., 2020) , multi-criteria based models (MCBMs) and probabilistic models, such as cross-impact analysis models (CIAMs) (Kadaifci et al., 2020) , Bayesian models (BMs) (George and Renjith, 2021) , risk-based approaches (RBAs) (Amin et al., 2019) , which are displayed in Fig. 1 . In these approaches, the impacts of various scenarios are measured in different ways. For instance, CIAMs, BAs, and RBAs use probabilistic inference to deal with a variety of interactive risk analyses based on probabilistic input data. They can be used to examine the effects of event probability changes, and then identify requisite actions and interventions to reduce adverse effects (Panula-Ontto & Piirainen, 2018) . In MCBMs and cognitive maps (Liu et al., 2019) , the cause and effect between various interactive concepts are modeled.

Various studies have been conducted to analyze the different aspects of risk analysis in distinct organizations, including Li et al. (2020) , who proposed a causal analysis based on the STAMP model for a safety risk analysis of underground gas pipelines. In this study, CAST analysis was adapted to examine safety flaws, revealing the series of reasons behind decisions made leading up to a catastrophic gas explosion. Yazdi et al. (2020a) augmented a new integrated approachbased on DEMATEL, BWM, and Bayesian network approachesto assess the dependency between risk factors and information sources. BWM was employed to compute relative expert opinion weights, then DEMATEL was mapped into the BN in order to identify critical factors in a dynamic structure. The proposed method was utilized in high-tech safety management. In another study, Yazdi et al. (2020b) introduced an improved solution, termed Pythagorean fuzzy DEMATEL, to evaluate the interrelation of corrective actions within a probabilistic safety analysis of an offshore platform facility. Pythagorean fuzzy numbers were applied to conjoin expert judgment and the DEMATEL method to encompass randomness and uncertainty. developed a fuzzy risk assessment methodology by integrating a fuzzy ANP and interpretive structural modeling, which captured interrelationships and interdependencies between risk factors and risk priorities to avoid data inaccuracy. ISM was used to identify critical risk factors, while fuzzy ANP was utilized to capture the fuzziness of neutral, optimistic, and pessimistic expert opinions, before ranking risk factors. This proposed method was implemented in construction and engineering risk management. A new risk analysis approach was proposed by Huang & Zhang (2020) combining FMEA and pessimistic-optimistic FAD, while accounting for an acceptable risk coefficient. The method was employed in evaluating dangerous goods transportation system risk on a railway. Tran et al. (2016) published a hybrid CIAM and factor analysis model for modeling cost variances in highway project-delivery decisions. Chen et al. (2019) proposed a model for the risk analysis of a multi-reservoir system, employing Monte Carlo simulations, dynamic Bayesian networks, and risk-informed inference. In this paper, Monte Carlo simulations were utilized to provide inputs for DBN, then DBN was built through expert knowledge uncertainty interrelationships. Finally, risk-informed inference provided risk information by a trained DBN. Drakaki et al. (2019) proposed a risk-informed supplier selection based on integrating fuzzy cognitive maps and a risk-based FAD (RFAD). Here, FCM was used to capture between-criteria dependencies, while RFAD was used to rank suppliers. proposed a three-stage approach integrating DEMATEL, ISM, and BN. As a first step, a hierarchical network model was presented combining DEMATEL and ISM methods, thus investigating the coupling relationships between various accident-related factors and BN structure. The hierarchical structure was mapped onto a BN in order to quantify the strength of relationships between accident leading systems, specifying the main resultant systemic cause.

Each available causality analysis method ( Fig.1 ) has certain drawbacks in terms of analyzing the main research problem. For example, existing risk analysis methods, such as fault tree (Kabir, 2017) and event tree (Purba et al., 2020) , only determine the probability of risk occurrence or the probability of occurrence of the desired outcome. They are thus unable to determine the rate of change and the new values of variables due to risk factors. In CIAMs, the effect of a risk factor is modeled only on probability. When the goal is to investigate the impact of a risk factor on the severity and probability of an event, this can be achieved by defining discretely different scenarios (states) and assigning probabilities to each of them. However, increasing states requires a large number of evaluations, causing the complexity of the problem to grow exponentially (Premchaiswadi, 2012) . Moreover, due to a consideration of discrete values, it does not provide a good estimate of the continuous values provided to decision-makers. Bayesian methods can be divided into two categories: discrete and continuous. Discrete methods have the same drawbacks as the CIAM, while continuous methods first require sufficient data to estimate the probability function. This is often unavailable in the problem, and is used less frequently by decision-makers due to its computational complexity. In definite methods such as MCDM, structural models, and existence of natural disasters, such as floods and earthquakes, can be considered an event. In the proposed R.Graph, an event is shown as and in the form of a rectangle.

Each variable or event is called a factor.

If there is a cause and effect relationship between two factors, the parent is the factor that affects the influenced factor.

Arc: This is a directional vector drawn from cause to effect, and shows the causal relationship. It is worth noting that R.Graph edges do not form a loop, because the proposed model is acyclic. The possible states and factors affecting the R.Graph method, as shown in Fig. 2 , are briefly examined in the following.

Since the present paper aims to investigate the extent of changes in the problem variables, and also since the problem is considered to be definite, and the event probability is deemed to be 1, the cause and effect events are represented in this state, according to Fig. 2 -a, that shows which events influence which other events, or are influenced by them.

in Fig. 2 

In this state, according to Fig. 2 -c, the goal is to examine the change in the variable due to the change in the variable .

Proposition 1. Suppose the changes in and are called and , respectively. Now, if is considered as a function of the variable , we will have:

Eq. (8) can be called the "change rate" to understand the changes better. In this study, the change rate of a variable is called the risk of that variable, which is generally shown by . Now, if we write the change rate of (risk of ) according to the change rate of (risk of ), we have:

| indicates the amount of risk due to the existence of the risk in . In this relation, the positive value of | indicates that by increasing the risk of , the risk of also increases (positive correlation), while the negative values indicate that increasing the risk of reduces the risk of (negative correlation).

By dividing by , the growth rate or the risk of is obtained, and we have:

is also written as a growth rate, i.e., , Eq. (9) is proven.■ For a real problem, the relationship between the variables is either already known (in which Eq.

(12) can be used to determine the risk of the variables) or unknown. If the relationship between the two variables is unknown, we can suppose the relationship between the two variables is . Using this, it can be shown from Eq. (9) that the relationship between the risks of two variables is as follows:

To correctly estimate the type of relationship between two variables and its coefficients, a sufficient amount of data is needed. However, many real-world problems are caused by lack of accurate knowledge of the function in terms of , since there is often not enough data to estimate the functions of the intended variables. In these cases, as is often applied to predictive intelligence in the case of little data (Deeb, 2015) , it can be better to utilize simpler models such as linear models. Considering linear relationships for many variables, such as the values of profit and cost, can be a rational assumption. In this paper, due to the lack of available data, and since it is assumed that all data are obtained from experts, an approximate linear method for determining the risk of with respect to the risk of is presented in the following.

Definition 4. Let the function be assumed to be linear in terms of . If the unit risk of is specified in terms of risk, then we have:

indicates the amount of risk per unit increase or 100% increase in the risk of , which can be obtained according to experts" opinions. It is worth noting that in this article, all variables are considered continuous and definite, and if the nature of the variable is probabilistic in the real world, its expected value can be entered into the problem.

In this state, according to Fig. 2 -d and assuming the model to be non-probabilistic, it can be said that if occurs as a parent of the variable , grows to a constant and definite amount.

Definition 5. If the risk of the variable due to the occurrence of is , then we have:

(12) ( | )

Definition 6. Suppose that a set of ́ variables and ́ events affects a specific variable, and ́ ́ (Fig. 3) ; now, if the purpose is to investigate the change rate (risk) of the variable in terms of all these factors, assuming all the factors are independent, we have: where represents all parents of , is the risk of the th variable and ( | ) is the change rate (risk) due to changes in or occurrence of . Now, a separate parent can be assumed for the members which can be seen in Fig. 4 . In this state, Eq. (12) can be generally written as follows:

In Eq. (13), the risk of the variable was calculated according to its parent. However, we can determine the risk of the variable relative to the risk of an influential variable.

Definition 7. If we want to define the risk of the variable relative to its influential variable , namely | , we will have:

where represents the variables that are directly or indirectly affected by ; shows the effect of events on that are the parents of , and affects their occurrence. Consequently, from Eq. (14), we have: ( | ) as follows:

In the above relation, shows the effect of events on , which is the parent of , and also shows that affects their occurrences. also indicates variables directly or indirectly affected by . On the other hand, it can be said that if does not affect in any way, we have:

After calculating the risk value for each variable, the new value of each variable can be updated using Definition 9.

Definition 9. Suppose the initial value of the variable is denoted by ; now, if the new value of is displayed as , according to the risk value ( | ), its value is obtained as follows:

The R.Graph data can be displayed through the R.Graph matrix as follows.

Definition 10. Assume that the whole problem includes variables and events; the basic values of the effectiveness of the variables (for 100% risk) and the events can be displayed by the R.Graph matrix which is denoted by . For simplicity, the may be defined as:

where the R.Graph matrix consists of 4 separate sub-matrices: , , and .

These are defined as follows:

where the matrix describes the impact of variables" risks on other variables; the matrix defines the impact of variables on events; the matrix describes the impact of events' risks on variables; and the impact of events" risks on other events are defined by the matrix. In the matrix, shows the risk of the -th variable per 100% increase in risk of the -th variable. Since the graph is acyclic, if adopts a value, we will have .

In the matrix, indicates the risk of the -th variable due to the occurrence of the -th event. Finally, in the and matrices, and indicate how the occurrence ofth event is affected by the occurrence of the -th event and -th variable, respectively. If and take values of one, indicate that the occurrence of the -th event is affected by the -th event and -th variable. If and are zero, it indicates that there is no effect. Here, also because the graph is acyclic, if and , we have and .

Example 1. Assume that the R.Graph matrix corresponds to Fig. 5 ;

The risk values of each variable can be computed as follows:

As mentioned in the introduction, in some cases, the only available data is expert knowledge. However, experts can make errors in judgment, which may lead to ineffective use and transfer of expert knowledge. Human decision-making has been shown to be inconsistent both between and within individuals. There are multiple possible sources for this inconsistency (Funk et al., 2020) : the problem type; the judge's decision criteria; uncertainty in the judge's knowledge; or randomness in judgments. Decision-making by experts may also exhibit random variation in situations with similar alternatives (Funk et al., 2020) ; there is also the possibility of systematic variation, depending on the amount of time available for experts to make choices (Funk et al., 2020) . This is also a factor when decision processes are context-dependent. In small samples with few experts, even minor inconsistency can cause large deviations in terms of commercialization choices (Funk et al., 2020) .

One of the things that should be considered in the R.Graph method is to check the inconsistency of the evaluations. On this basis, after obtaining the results through the R.Graph method, the following condition should be considered for each variable:

where and are the lowest and highest possible risk values that ( | ) can take, respectively.

In the proposed method of R.Graph, if there is a risk in a parent, the risk is passed on to its offspring, which increases the risk values of the downstream variables. Since one of the goals of the R.Graph method is to identify events that affect critical variables and plan to reduce their risk, high risk demands greater efforts in decreasing the risk. Therefore, in some decision-making problems, decision-makers accept some level of risk due to organizational goals and the degree of risk-taking, which is called "acceptable risk" (AR) in the literature (Seiti et al., 2019) . In this paper, the acceptance level of risk and changes in a variable is defined as the percentage change for a one hundred percent risk in the desired variable that can be accepted and compensated by the relevant organization which is not included in the calculations. This can be expressed as:

In the R.Graph method, some of the risk levels of each variable can therefore be considered acceptable. However, since the risk value of one variable depends on the risk values of other variables, the acceptance and compensation of some risks affects the risks of other variables, so the risk values of each variable should be calculated by considering the acceptable risk of the variables. Therefore if the AR value is defined as a percentage between 0 and 100, each risk value can be modified based on Definition 11.

Definition 11. Assume that the risk-taking for a one hundred percent risk of the -th variable is shown by , considering the level of risk-taking, Eq. (13) can be rewritten as follows:

where the values of acceptable risk are considered simultaneously for inputs and outputs, and is the percentage of risks not accepted and entered into the problem. Generally, the matrix of acceptable risk for variables, which is displayed by , can be defined as follows: 

In some decision-making and risk analysis problems, the decision-maker is interested to know the value of the highest and lowest risks that will be experienced in each variable, terming them pessimistic and optimistic values, respectively. Since the risk values of each variable or event on a specific variable are either positive or negative, Proposition 2 can be formulated.

Proposition 2. The maximum and minimum amount of risk of a variable can be determined as follows:

shows the highest value of risk (pessimistic risk) and shows the lowest amount of risk (optimistic risk) of the variable . Also, in these relations, and represent the effects of event risks on variables which are positive and negative, respectively.

Moreover, and indicate positive and negative values, which are defined as follows: 

One of the goals of cause and effect methods, including the proposed R.Graph method, is to identify the events or variables with the most significant impact on other variables, as considered by an organization. Indeed, the organization or managers can identify these critical factors and plan to reduce related negative consequences. The following outlines the prioritization of factors in the proposed R.Graph method.

Proposition 3. Assume that the total number of variables is and the total number of events is ; then, the relative importance of factors is obtained from the following relation:

where shows the weight of the -th variable and shows the weight of -th event among all factors.

Proof. The change rates in the risk of all variables can be written, considering the -th variable, as follows: Fig. 5 , and its binary display, which is indicated by arrows, can be seen in Fig. 6 . To find the effect of the variable on the problem, it is necessary to find the effects of the variable on other variables, and then remove the variable from the problem. This is shown by drawing the yellow lines in Fig. 6 , stating that removal only affects the variable . Now, the weight of the variable can be determined as follows: Fig. 6 . Variable removal and the interaction of factors in Example 1

This section examines how to find the sensitivity of the changes in one variable over other problem variables.

showing that by increasing the constant value per input risks, the problem is more sensitive to risk changes in .

As mentioned in the introduction, the purpose of developing the R.Graph approach is to provide a cause and effect model using expert knowledge. In fact, it aims at analyzing risk and changes in important organizational parameters, based on risk factors with a low probability of occurrence, to provide preventive solutions and relative preparedness. The way of modeling in R.Graph method and its mathematical relations were presented in Section 4.1. Additionally, different cases of the risk analysis were considered, such as considering acceptable risk, optimistic and pessimistic situations, the importance of factors, and sensitivity analysis. The following section aims to provide a group risk analysis framework by experts, based on the proposed R.Graph method, and is summarized as follows:

Step zero. Identifying the set of events and variables affecting each other, and the degree of acceptable risk.

In this step, the organization"s intended variables are considered, examined, and determined according to organizational purposes. In addition to the factors leading to changes in intended variables or other events, the organization also wants to investigate how these effects are specified. In this step, the interrelations between variables and events are determined using expert opinions and by a moderator. Then, a graph of causes and effects is drawn, from which the event-variable matrix (i.e., ) and the event-event matrix (i.e., ) can be created.

Additionally, the experts are asked to determine the effects and risks of events on variables Phase IV: Post-processing

Step 5. Calculating the weights of all variables and events

In this step, Eq. (30) is used to determine the importance of each factor, whether variable or event.

Step 6. Sensitivity analysis

To identify the sensitivity analysis of different factors in their input parameters, considering different values for and , Eqs. (31) and (32) 

The COVID-19 epidemic in Iran is part of the worldwide Coronavirus pandemic, which has had various consequences, including social and economic implications, for the entire country. One of the most critical parts of Iran's economy that has been significantly affected by the Coronavirus is the electricity industry. COVID-19 is projected to reduce production volumes and shut down manufacturing, contracting, and consulting units, leaving manufacturing units temporarily closed and costing a large number of jobs (Asadi, 2020). The volume of executive activities will be reduced, and ongoing projects may also be delayed; with the disruption of trade, imports and exports related to the electricity industry will additionally be severely disrupted (Asadi, 2020). The power grid is unstable, and blackouts occur. Delays in the private sector, due to declining production, will lead to the closure and bankruptcy of manufacturers and contractors in the electricity industry, and ultimately the electricity industry economy will suffer direct and indirect losses due to delays in bank payments, social security, and taxes (Asadi, 2020) . In this study, the aim is thus to investigate the reduction of profitability of Iran's electricity industry due to the Coronavirus outbreak, and its consequences, over a one-year period. 4 categories of events and 13 variables have been identified in the industry which are directly and indirectly affected by the Coronavirus pandemic, as shown in Table 1 . Additionally, the causal relationships of these factors have been identified by relevant experts and drawn into the J o u r n a l P r e -p r o o f R.Graph Diagram (Fig. 8) , which led to the determination of the (Table 2 ) and matrices (Table 3) . The six steps of the proposed methodology (Section 4.2) have been utilized to evaluate the risk of each variable in different scenarios, and risk analysis has been carried out using the R.Graph method. 

In the initial stage, by determining the R.Graph diagram of Fig. 8 , two sub-matrices of the R.Graph matrix (i.e. and ) were first determined in order to obtain the other two matrices, namely the matrices of the event on variable effect ( ) and the variable on variable effect ( ). Three experts from the field of the electricity industry were recruited to cooperate with the study. The participants were then analyzed demographically. Two experts were held doctorates in related fields, while one held a related postgraduate degree. The decision-makers had high levels of experience in the field, ranging from 5 to 25 years. As per the data collected, the decision-makers involved in this study were the researchers from two research centers (namely from a high-capacity power transmission center, and from a group planning and operating power systems). In order to analyze the risk of each variable, the three experts were asked to specify the impact and risk of potential events on each of the affected variables. They were also asked to estimate the effect (risk) of each variable on other variables, in case of a 100% increase in the effective variable. For example, experts were asked the following question to determine the effect of the event "new safety regulations" on the variable "job difficulty":

How many percentage points do you think new safety regulations will alter job difficulty? In order to determine the effects of variable "personnel medical expenses" on "total cost", the following question was asked:

In both types of questions, participants were requested to indicate an increase in the percentage changes with a positive sign, and a decrease with a negative. The questionnaire used for obtaining the and matrices of each expert can be found in Appendix 1. Next, since the and matrices in the R.Graph method were developed based on the unit value of the changes, after determining the percentage impacts by experts, the evaluation values were converted to unit values in terms of percentage (dividing by 100). Finally, equal weights were considered for each ordered aggregation, and equal weights given to each expert. The and matrices obtained from each expert were aggregated using Eqs. (33) and (34), and can be seen in Tables 4 and 5. In the next stage, an expert in the organization"s policies was asked to determine acceptable risk values for each of the variables in Table 6 . For example, the related expert was asked:

The questionnaire used for obtaining acceptable risk values can be found in Appendix 2. 

At this stage, the risk in different situations was calculated, taking into account the potential risk with and without considering the acceptable risk, and while calculating the maximum and minimum risks. A software code based on Algorithm 1 was developed using MATLAB software; the results can be seen in Table 7 Table 7 . Risk values of each variable in different cases Now, since the risk results in the "without considering" mode are consistent, the analysis proceeds to the next section in order to carry out the weight calculation and sensitivity analysis.

In this section, in order to prioritize all factors for planning risk management and determining preventive measures, the importance of each factor was determined and ranked using Eq. (30), as shown in Table 8 . Also, to investigate problem sensitivity to increases and decreases in the constant value in the risk of various factors, three cases of sensitivity analysis were considered as follows:

Case 1. To determine the risk of each of the variables for each factor, the effect of each of the affecting factors is considered to be 0.1 greater, thus: Table 8 and Fig. 9 . 

An analytical method"s robustness may be defined as an estimation of its capability to remain unaffected by small but deliberate methodological changes across variables (César and Pianetti, 2009 ). The sensitivity analysis in Section 5.3 showed how much a constant change in values of a factor would change the aggregated risks of all variables. This section investigates to what degree small changes in the most influential factors will have an impact on each specific result. Therefore, the most influential factors should first be identified; as can be seen in Table 8 and Figure 9 , the greatest sensitivities are related to factors , , , , and .

Since one of the main outputs of the R.Graph method is the ranking of factors, some dummy values were deliberately entered into each of these factors in order to measure their impacts on the overall ranking of factors. For this purpose, three error levels of 2%, 4%, and 6% were entered into each of the values related to each of the factors. Then the effect of each factor on the results was first examined individually, then the changes entered into a set of influential factors. This means that each of these sets entered 2%, 4%, and 6% errors in their values, respectively, and the result of the overall ranking of all factors in all these cases was obtained, as can be seen in Tables 9-11. In each of these tables, the ranks that differ from the main case are distinguished in blue. Additionally, to check the closeness of these answers with the answers of the main case, Spearman's correlation coefficient (SCC) (Seiti et al., 2020b) was applied, which is a measure of the correlation between the ranking and its values. These SCC results can also be seen in Tables  9-11 . Table 9 . The robustness analysis results at 2% error level 

In this section, the aim is to compare the results obtained from the R.Graph method with other existing methods from the relevant literature. Since the proposed method is definite, it is necessary to examine existing definite approaches in the literature, some of which were previously mentioned in Section 2. In general, though, the results of the proposed R.Graph are not comparable to any of the existing methods, since each one has its own assumptions with different input meanings. However, with some simplifying assumptions, it is possible to compare the EXIT (Panula-Ontto and Piirainen, 2018) and fuzzy cognitive maps (Stylios and Groumpos, 2004 ) methods, which are most similar to the present method, assuming that the inputs of the R.Graph have the same interpretation as the inputs in these methods. Table 12 presents the results of comparing the obtained ranks from the R.Graph method and the ranking results based on the total impact measure of the EXIT method and centrality measure of fuzzy cognitive maps approach. 

As discussed in Section 4, the risk analysis horizon was considered for the case study for a period of one year. However, at the time of writing this article, 6 months have elapsed since the risk analysis was conducted. Thus, Table 13 is provided in order to compare how close the R.Graph risk prediction results obtained were to the changes observed after 6 months, according to available data. Of course, it is worth noting that the values of risk analysis are substantially different from the observed values, due to the existence of the acceptable risk factor. Moreover, using forecasts as inputs in decision-making processes often results in self-predicted outcomesa problem known as self-negating forecasts, or the prophet dilemma (Petropoulos et al., 2020) . However, it can be seen that the predicted and observed values in many cases are relatively close. 

The R.Graph method was developed to provide a framework for risk analysis based on the degree of predicted changes in the variables of interest to an organization, taking into account managerial factors, such as acceptable risk factors. However, it can be a useful tool to investigate the effects of various corrective actions on different variables risks, allowing for appropriate planning to implement preventive measures. In this part, some corrective actions to reduce possible risks in the variables are suggested by the experts, as can be seen in Table 14 . In addition, the interactions of these actions on different variables are depicted in Fig. 10 , in which corrective actions are shown as events.

The aim now, then, is to rank these actions based on their suitability in reducing the risks of the entire problem. It should be noted that the effects of corrective actions are examined even when there is no acceptable risk in the problem, because in this case, the effect of corrective actions on the absolute risk is modeled. It is worth noting that any corrective action which has a positive effect on improving the risk value of one variable may have negative effects on the risk values of other variables. If the percentage of improvement and worsening in the risk value of variable due to taking action are denoted with and , the modified amount of risk assuming corrective action can be obtained as follows: (36) where denotes the amount of risk without considering any acceptable risk value, and is the modified risk value. Also, we have:

Finally, the total effect of -th corrective action on the whole problem, which is shown by , can be defined as: The total impacts of all suggested actions are presented in Table 14 .  Table 8 shows the importance of each factor on the total risks of problem variables, and also their ranking. It can be seen that the Coronavirus pandemic is the most important of all the factors that could have already been predicted, because all the values of the risks were due to its occurrence. It is also observed that the highest importance among the variables is related to the power branch sales revenue, and the lowest concern belongs to the variable of degree of work difficulty. Decision-makers and managers in the organization can thus design and implement preventive plans according to the importance of each individual factor, in order to reduce its risk and its consequences.

 For many problems, modeling with traditional techniques is difficult or even impossible. This is especially true when certain statistical data related to the model and its parameters are not available, and/or it is difficult to extract the relationships between the components of the modeled system from the quantitative prediction models such as time series models. In such cases, opinions provided by knowledgeable individuals are the best available data. The R.Graph method can be seen as a tool for modeling and simulating systems, using data obtained from experts in the field, thus providing an appropriate, interpretable framework for decision-makers and system modelers. In this method, the intended data can be collected and aggregated through interviews, workshops, or surveys, or by using the Delphi method and questionnaire, or a combination of these methods. Additionally, experts can be asked to vote on entries anonymously using an online questionnaire, then to discuss the results directly.

 In the proposed R.Graph method, information about indirect interactions between system components is extracted based on direct interactions. Indeed, direct interactions are input data, and the model"s utility lies in analyzing indirect interactions. This is especially useful, since, in a complex network with a variety of factors, incidental interactions and their causal chains can be complicated and lengthy. Additionally, the proposed R.Graph method can identify the importance of a system component that may not appear to be related to another element, such as when some seemingly essential components may be canceled or neutralized by network interactions.

 Generally speaking, interpretability of an algorithm refers to how causes and effects may be observed within a system, allowing for the ability to predict outcomes (Richard Gall, 2018) . Explicability, on the other hand, refers to how the internal mechanics of a model can be expressed in simpler human terminology (Richard Gall, 2018) . Algorithms based on the concept of causality, such as the proposed R.Graph method, have the ability to be easily interpreted, and the results can be explained to decision-makers. Therefore, its results are more reliable for managers and decision-makers, as they can participate in all stages of analysis.

 In Section 5.6, the results of the R.Graph method were compared with two deterministic methods in the literature that are similar to this method, and the results were relatively close. However, it should be noted that the R.Graph method is logically different from other existing approaches. The most important differences between R.Graph and other causal methods can be seen in: the difference of system inputs (considering change instead of scenario); the difference of interactions between factors (how to consider the relationships between factors); assuming a definite nature instead of a probabilistic one; the display of various factors; and considering a static nature. Indeed, the proposed R.Graph risk concept was developed assuming possible variations in model variables, so that it can be used in cases where the decision-maker is interested in investigating the effect of different events on a series of specific variables, or is interested in examining the rate of effect due to unforeseen events. Other advantages of the R.Graph method include: considering acceptable risk levels for desired organizational parameters; determining the highest and lowest risk values of each variable; determining the weight and ranking of factors; and presenting an appropriate method for model sensitivity analysis. A selection of causal models have been compared with the proposed R.Graph from different perspectives, and the findings are depicted in Table 15 . 1) The first simplification assumption is that in the proposed model, events are assumed to be definite and the probability of occurrence was not considered for each event. This is because definite consideration of the occurrence of events and their impact is based on the assumption that it is known that the event in question will occur sooner or later; the purpose is to determine its impact. However, in the real world and especially in risk management over longer periods of time, this is not the case; the probability of each event should be considered for each variable.

2) The second simplification assumption is that the relationships between the variables are considered to be unknown and statistical data are unavailable. Since it was assumed that the data was obtained through experts, the risk of each variable was considered in terms of the linear sum of changes in the influential variables, as well as in terms of influential events which affect a variable independently. However, the process to obtain the risk of each variable in a general case was described in Section 4. In general, how much data is needed to reliably obtain the risk of one variable based on other variables depends on the degree of relationship between the variables and the number of variables: this can involve the use of statistical methods such as regression analysis (Fox, 1997) .

3) The third simplifying assumption was that the initial and subjective weight was not considered for any of the risk factors (i.e., for both events and variables the same initial weight was considered). The initial weight can be easily determined by extending Eq. (30) on this basis.

 An important point to be noted is that in this study, when collecting data from experts, attention was paid to the quality of experts and not the quantity, although the number and access to experts in this field were limited. Therefore, due to the lack of sufficient statistical data, it was not possible to analyze the dimensions related to gender and age; only the quality and subjective importance of each expert was aggregated by the HWA operator and entered into the problem. In group decision-making, there is no specific rule for determining the minimum number of experts needed to achieve a reliable result, but if the decision-maker believes in the high reliability of data from statistical methods, more experts guarantee a better and more reliable result.

 In this section, the advantages and disadvantages of the proposed R.Graph method in risk analysis and management are briefly mentioned, and its differences and similarities with existing methods were discussed. It should be mentioned that the development of the proposed method and its innovations, with all its positives and negatives, is a step toward a longer-term goal and has been developed to provide new ideas for researchers in the relevant field and decisionmakers in their future decision-making issues.

Risks and disruptions, especially unforeseen ones, have led many companies today to be concerned about the consequences on the complex, global business environment. The Coronavirus pandemic represents such a specific disruptive case. Epidemic/pandemic outbreaks lead to long-term disruptions, which are unpredictable in terms of time, severity, and scale. Even with certain known, low-probability, high-impact events, many traditional risk analysis methods have proven insufficient in assessing their risks and outcomes. The present paper presented a new risk analysis method to measure the risk of interactive factors in the causal chain which can be fully explainable for decision-makers. Its application was studied in a risk analysis of Coronavirus on the Iranian electricity industry. The advantages of the proposed model included: determining the amount of change in variables; accounting for preferred organizational parameters, such as acceptable risk; determining the importance and ranking of factors; analyzing sensitivity; and interpretability and explicability for decisionmakers.

Simplified assumptions were made in the development of the R.Graph method in order to better comprehend this new model and its applicability. Consequently, new assumptions can be introduced as future research topics developing upon the current study. For instance, the R.Graph model has been developed assuming that the problem in hand is definite, so its extension into a probabilistic model to capture randomness can provide useful future research. Moreover, linearity is one of the main assumptions in the R.Graph method, but the linearity assumption is not always established in certain risk analysis problems. It is also suggested to employ supervised (Osarogiagbon et al., 2020) or unsupervised learning (Michau and Fink, 2021) methods in cases where sufficient data is available for statistical interpretation. Another issue not discussed in the R.Graph method is how to combine the information obtained from experts with methods to consider the reliability of their evaluations ; this can be considered as another future research area. Finally, other areas of research might pursue the development of the R.Graph model in a dynamic mode (Wang et al., 2021) , while considering other uncertainties (Jianxing et al., 2021) in input data.

Risk assessment of ex-post transaction cost in construction projects using structural equation modeling

A bibliometric review of process safety and risk analysis

Explainable Artificial Intelligence (XAI): Concepts, taxonomies, opportunities and challenges toward responsible AI. Information Fusion

Modeling and analyzing factors affecting project delays using an integrated social network-fuzzy MICMAC approach

Risk analysis for real-time flood control operation of a multi-reservoir system using a dynamic Bayesian network

Robustness evaluation of the chromatographic method for the quantitation of lumefantrine using Youden's test

What to do with "small" data? 15AD

Building resilience against biological hazards and pandemics: COVID-19 and its implications for the Sendai Framework

Supplier selection problem in fuzzy environment considering risk factors

Applied regression analysis, linear models, and related methods

Individual inconsistency and aggregate rationality: Overcoming inconsistencies in expert judgment at the technical frontier

Evolution of Safety and Security Risk Assessment methodologies to use of Bayesian Networks in Process Industries

Can a risk have 100% probability?

Future scenarios: The BASICS computational method

Railway dangerous goods transportation system risk assessment: An approach combining FMEA with pessimistic-optimistic fuzzy information axiom considering acceptable risk coefficient

Process system failure evaluation method based on a Noisy-OR gate intuitionistic fuzzy Bayesian network in an uncertain environment

A new 2-additive Choquet integral based approach to qualitative cross-impact analysis considering interaction effects

An overview of fault tree analysis and its application in model based dependability analysis

Analysis on accident-causing factors of urban buried gas pipeline network by combining DEMATEL, ISM and BN methods

A CAST-based causal analysis of the catastrophic underground pipeline gas explosion in Taiwan. Engineering failure analysis

Fault tree analysis for integrated and probabilistic risk analysis of drinking water systems

Some hybrid weighted averaging operators and their application to decision making

A novel approach to fuzzy cognitive map based on hesitant fuzzy sets for modeling risk impact on electric power system

Using fuzzy analytic network process and ISM methods for risk assessment of public-private partnership: A China perspective

Unsupervised transfer learning for anomaly detection: Application to complementary operating condition transfer

The world economy at risk n.d. www.oecd.org › berlin› Interim-Economic-Assessment

Review and analysis of supervised machine learning algorithms for hazardous events in drilling operations

AXIOM Method for Cross-Impact Modeling and Analysis (Master's thesis)

EXIT: An alternative approach for structural cross-impact modeling and analysis

Forecasting: theory and practice

Fuzzy probability based event tree analysis for calculating core damage frequency in nuclear power plant probabilistic safety assessment

Interpretability: Two concepts that could help restore trust in AI n

R-numbers, a new risk modeling associated with fuzzy numbers and its application to decision making

A New Risk-Based Fuzzy Cognitive Model and Its Application to Decision-Making

A novel linguistic approach for multi-granular information fusion and decision-making using risk-based linguistic D numbers

Developing the modified R-numbers for risk-based fuzzy information fusion and its application to failure modes, effects, and system resilience analysis (FMESRA)

Modeling complex systems using fuzzy cognitive maps

How should relative changes be measured

A hybrid cross-impact approach to predicting cost variance of project delivery decisions for highways

Using VIKOR with structural equation modeling for constructing benchmarks in the Internet industry

An approach to evaluate the methods of determining experts" objective weights based on evolutionary game theory

FMDBN: A first-order Markov dynamic Bayesian network classifier with continuous attributes. Knowledge-Based Systems

Reliabilities analysis of evacuation on offshore platforms: A dynamic Bayesian Network model

Improved DEMATEL methodology for effective safety management decision-making

A novel extension of DEMATEL approach for probabilistic safety analysis in process systems

Identifying critical risks in Sponge City PPP projects using DEMATEL method: A case study of China

System modeling and control with resource-oriented Petri nets

in Section I, and fill the response sheet in Section II. All personal information will remain fully confidential, and will not be shared with anyone

The present investigation is supported under the Iran University of Science and Technology, Tehran, Iran, Postdoctoral Fellowship granted by the Iran National Science Foundation (INSF), and Iran's National Elites Foundation, Iran (Grant No. 99003878).

 According to Table 8 , the effect of increasing or decreasing the values of each factor on the risk of its children can be examined, which is termed the amount of sensitivity analysis of that factor. It can be observed that the highest sensitivities to changing values are seen in power branch sales revenue in variables. Moreover, it is also observed that the problem is very sensitive to changes in the new safety regulations, in addition to the Coronavirus pandemic. It can also be seen that, for each factor, the following is apparent:Generally, it can be said that the problem sensitivity to increases or decreases in a certain value are opposites in the R.Graph method, thus: To investigate the degree to which the model ranking results are robust subject to minor changes, three different scenarios of change in the values of the influencing factors were carried out. The results are presented in Tables 9 to 11 ; it can be seen that the Spearman's correlation coefficient values changed from 0.898 to 0.916 in worst and best cases. This indicates slight changes in the ranking results, and an acceptable degree of robustness within the model. In Section 5.7, a framework for evaluating the effectiveness of corrective measures was presented by the R.Graph method through examining which measures had the highest priority in reducing the overall risk of the problem. It was observed that periodic corona tests were the first priority. The proposed framework can be a useful tool for ranking various factors, allowing for appropriate planning to implement preventive measures. It permits decision-makers to effectively determine policies and interventions affecting system outputs, and to examine their impact on reducing adverse effects or increasing.

The COVID-19 epidemic in Iran is part of the worldwide Coronavirus pandemic, which has had various consequences, including social and economic implications, for the entire country. One of the most critical parts of Iran's economy that has been significantly affected by the Coronavirus is the electricity industry. COVID-19 is projected to reduce production volumes and shut down manufacturing, contracting, and consulting units, leaving manufacturing units temporarily closed and costing a large number of jobs. The volume of executive activities will be reduced, and ongoing projects may also be delayed; with the disruption of trade, imports and exports related to the electricity industry will additionally be severely disrupted. The power grid is unstable, and blackouts occur. Delays in the private sector, due to declining production, will lead to the closure and bankruptcy of manufacturers and contractors in the electricity industry, and ultimately the electricity industry economy will suffer direct and indirect losses due to delays in bank payments, social security, and taxes. The purpose of this questionnaire is to evaluate the impact of COVID-19 on the intended variables in the electricity industry. We would be grateful if you would respond to the questions Authors 

The purpose of this questionnaire is to evaluate the impact of COVID-19 on specific variables in the electricity industry, which is determined by the percentage of variability (percentage decrease or percentage increase), and the risk of that variable.

Note that if a factor or change in one variable causes a positive change in another variable, please indicate it with a positive sign, and use a negative sign if it leads to a decrease. For example:

Your comment: + 50%, which means that you think it would increase by 50%.Therefore, you are asked to determine the effect of COVID-19 on each of these factors, according to the instructions mentioned.

How many percentage points do you think new safety regulations will alter job difficulty? How many percentage points do you think new safety regulations will alter environmental health costs? How many percentage points do you think new safety regulations will alter levels of stress? How many percentage points do you think staff infections will alter number of key personnel? How many percentage points do you think staff infections will alter personnel medical expenses? How many percentage points do you think first-degree relatives being infected will alter number of key personnel? How many percentage points do you think first-degree relatives being infected will alter personnel medical expenses? How many percentage points do you think the coronavirus pandemic will alter the percentage of receivables? How many percentage points do you think the coronavirus pandemic will alter energy sales revenue? How many percentage points do you think the coronavirus pandemic will alter power branch sales revenue? How many percentage points do you think a one hundred percent increase in personnel medical expenses will alter total cost? How many percentage points do you think a one hundred percent increase in environmental health costs will alter total cost? How many percentage points do you think a one hundred percent increase in job difficulty will alter total personnel efficiency? How many percentage points do you think a one hundred percent increase in levels of stress will alter job total personnel efficiency? How many percentage points do you think a one hundred percent increase in number of key personnel will alter total personnel efficiency? How many percentage points do you think a one hundred percent increase in total personnel efficiency will alter project delay rate? How many percentage points do you think a one hundred percent increase in project delay rate will alter total income? How many percentage points do you think a one hundred percent increase in percentage of receivables will alter total income? How many percentage points do you think a one hundred percent increase in energy sales revenue will alter total income? How many percentage points do you think a one hundred percent increase in power branch sales revenue will alter total income? Appendix 2. Questionnaire B Dear respondents Greetings.The COVID-19 epidemic in Iran is part of the worldwide Coronavirus pandemic, which has had various consequences, including social and economic implications, for the entire country. One of the most critical parts of Iran's economy that has been significantly affected by the Coronavirus is the electricity industry. COVID-19 is projected to reduce production volumes and shut down manufacturing, contracting, and consulting units, leaving manufacturing units temporarily closed and costing a large number of jobs. The volume of executive activities will be reduced, and ongoing projects may also be delayed; with the disruption of trade, imports and exports related to the electricity industry will additionally be severely disrupted. The power grid is unstable, and blackouts occur. Delays in the private sector, due to declining production, will lead to the closure and bankruptcy of manufacturers and contractors in the electricity industry, and ultimately the electricity industry economy will suffer direct and indirect losses due to delays in bank payments, social security, and taxes. The purpose of this questionnaire is to evaluate the impact of COVID-19 on the intended variables in the electricity industry. We would be grateful if you would respond to the questions in Section I, and fill the response sheet in Section II. All personal information will remain fully confidential, and will not be shared with anyone.

Section I General expert information Please select only one option 1. What is your highest qualification level? 

The purpose of this section is to determine the risk acceptance of each influence and influential factor in the electricity industry.

The acceptable levels of risk and changes in a specific variable means the percentage change due to the occurrence or change in the desired variable which can be accepted and/or compensated by the relevant organization and thus not included in their calculations. This denotes the degree of risk acceptable to the organization, which can be expressed as a percentage between 0 and 100. We need your judgment to define the acceptable risk for each variable, so please fill out the following table, based on your organizational policy, using values between 0 and 100. 

☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: