Abstract
At present, multi-attribute group decision making (MAGDM) researches are concerned with obtaining the alternatives’ ranking order based on the gathered group fuzzy preference matrix, as a result, it suffers from the group experts’ strong subjectivity fatigue without approaches to obtain and use the group intelligence. This paper proposes a new definition of group intelligence based on sub-attribute’s weight coefficients and expected intervals of the MAGDM. Group intelligence generating methods are provided as well. Additionally, a kind of group intelligence-based decision approach is introduced that can aid in reducing subjective fatigue in humans. In particular, a test function and a kind of medical MAGDM problem addressing Traditional Chinese medical (TCM) dysphagia rehabilitation treatment are employed to verify the correctness and feasibility of the proposed methods.
Keywords
Introduction
Currently, MAGDM [33, 34] approaches are more concerned with the decision-making results than the method employed to obtain and use the group intelligence based on the sub-attribute’s weight coefficients and expected intervals of MADM problems. Practically, traditional MAGDM methods require every expert to give his preferred judgment preference for each decision-making alternative in order to obtain the alternatives’ ranking order; however, it is important to note that these preferences are influenced by the experts’ strong subjectivity and fatigue.Therefore, it is necessary to gain further insight into group intelligence, which will lay a solid foundation for group intelligence generating techniques.
MAGDM problems [8, 24] are employed in many complex decision making areas [7, 16]. Additionally, the expert’s judgment preferences are always complex, and as a result, there have been many studies conducted on group experts’ judgment preferences gathering mechanisms [22, 25]. For example, the new determination model [2, 30] was proposed to obtain the weights of different attributes and the group consistency index. Additionally, methods on calculating the fuzzy decision matrix in the fuzzy environments were introduced [11, 38]. To describe preferences in the MAGDM, linguistic preference [17], interval preference matrices [3, 9], and intuitionistic fuzzy preference [4] have been discussed in regards to their impact on decision-making results. In addition, to obtain the group experts’ preference, many operators have been proposed [12, 18]. However, most of these studies focused on providing methods for solving MAGDM problems without addressing the characteristics of their sub-attributes (each sub-attribute’s weight coefficient and expected interval). Furthermore, traditional MAGDM studies require each expert to give his preferred judgment preference (such as, a utility vector, ordering relationship value vectors, etc.) for each decision-making alternative [19]. Meanwhile, the gathered group fuzzy preferences matrix cannot describe the group intelligence of the specific MAGDM problem visually, it suffers from the group intelligence’s definition, and its generating approaches aid computers in making an automatic decision.
One such MAGDM problem is the rehabilitation treatment of dysphagia. Dysphagia refers to a dysfunction that occurs when food is transported from the mouth to the stomach. This dysfunction of the lower jaw, lips, tongue, throat, and esophagus can cause serious consequences, such as aspiration pneumonia and asphyxia. The anatomical structures of dysphagia are shown in Fig. 1.
Based on the dysphagia diagnostic results (the mouth period, tongue period, pharynx period, and larynx period) [20, 21], dysphagia rehabilitation treatment [13, 29] usually incorporates TCM techniques [1, 32] (acupuncture, massage, muscle physiotherapy, etc.) to promote muscle function (soft palate, epiglottis capacity, tongue, etc.) and to stimulate the swallowing nervous system (Fengfu acupoint, Lian Quan acupoint, Yifeng acupoint, etc.). Practically, dysphagia rehabilitation programs involve multiple treatment objectives (multiple muscle and nervous systems) as well as a variety of therapeutic attributes (acupuncture, massage, muscle physiotherapy, etc.); however, it is difficult to quantify the relationship between each treatment objective and therapeutic attribute. Therefore, the dysphagia rehabilitation treatment is a kind of MAGDM problem. Additionally, to ensure the safety and effectiveness of rehabilitation, the current treatment program drafting needs doctor groups to discuss all treatment alternatives for each patient. Furthermore, it is very difficult to build and learn the rehabilitation experts’ treatment decision-making intelligence, where the scientific basis of the treatment decision-making should be improved. Therefore, it is necessary to build a quantitative model of the group intelligence and propose its generating methods. Such a case can greatly benefit from an improved interactive generating method for group intelligence.
In regards to some decision processes, especially for TCM clinical problems without a structured description of the objective function, group experts are often required to provide experience and to assist with the decision making. Therefore, it is particularly important to provide group intelligence generating methods and decision strategies based on group intelligence. With this motivation, based on either the characteristics of the sub-attributes or the weight coefficient and expected interval of each sub-attribute in the MAGDM problems, a new definition of group intelligence is proposed to demonstrate the generating process of group intelligence in this paper. In addition, a group intelligence generating algorithm and a group intelligence-based decision making approach are introduced, which can help gradually learn [23, 24] essentials in a group’s decision intelligence, reduce human subjective fatigue, as well as make the decisions more scientific and objective.
In this study, firstly, we proposes the definition and generating methods of group preferences in Section 2, which are followed by the group intelligence-based decision making approach. Secondly, the applications consisting of test functions and traditional Chinese medicine (TCM) dysphagia rehabilitation treatment issues are provided in Section 3. Furthermore, the performance of the proposed methods is discussed in Section 4. In the end, the conclusions and the future perspectives are shown in Section 5.
Group intelligence generating methods
Tradition MADM methods often base the decision solution on utility function and probability approaches. However, the mathematical descriptions of an individual person’s preferred forms judgment preference are always complexes, which include utility vectors and so on [19]. The fuzzy MADM approach is always used to solve MAGDM problems.
The fuzzy MAGDM methods tend to focus on the selection of alternatives, ordering, and evaluation issues based on a group fuzzy preferences matrix, which is incapable of obtaining weight coefficients and expected intervals of the sub-attributes. In contrast to group fuzzy preferences matrix, a definition of group intelligence is proposed to describe an alternative to the sub-attribute weight coefficient and expected interval in MAGDM problems. Therefore, the group intelligence generating algorithm is designed under the Genetic Algorithms (GA) framework.
Background of fuzzy MAGDM methods [12, 19]
In traditional fuzzy MAGDM methods, it is important to convert the person’s different complex judgment preferences into the fuzzy complementary judgment matrix , where , and k is the number of decision-makers. The fuzzy complementary judgment matrices are gathered into the group fuzzy preferences matrix . Additionally, based on the gathered group fuzzy preferences matrix, the key step for obtaining the ranking order of the alternatives is shown as follows:
In order to create a new document, do the following:
Set k decision-makers, ξ sub-attributes corresponding to each alternative and n alternatives in MAGDM problems. In addition,
(1) Conversion methods translate experts’ complex judgment preference of a preferred form [19] to a fuzzy complementary judgment matrix
(a) The utility vector
(b) The ordering relationship value vector (where ) can be converted to a fuzzy complementary judgment matrix:
(c) The reciprocal judgment matrix can be converted into a fuzzy complementary judgment matrix to achieve unification:
Where, φ′ is the maximum possible scaling value of the reciprocal judgment matrix [19]. In particular, the scaling value is 9 in this paper; thus, c = 81.
(2) This section details gathering judgment preferences to the group fuzzy preferences matrix.
(a) The ordered weight averaging (OWA) operator [24] is described in the following.
Firstly, Q (γ) is the fuzzy quantification operator defined by Equation 4.
Where, 0 < γ ≤ 1, (α, β) is designed as the decision principle of “majority,” “half at least,” or “as much as possible,” which corresponds to (0.3, 0.8), (0, 0.5), and (0.5, 1), respectively.
Secondly, the OWA operator weight that denotes the lth decision-maker’s weight is expressed as:
Additionally, the OWA operator weight is constrained by β = 1.
(b) The gathered group fuzzy preferences matrix is described in the following.
We designed the group fuzzy preferences matrix
Based on the OWA operator weight vector, denotes the superior degree that the ith alternative is better than the jth alternative in the group judgment, expressed as:
(3) This section details the ranking order of the alternative.
(a) The weight vector of the OWA operator for the qth alternative is calculated as follows.
Based on Equation (5), the qth alternative’s OWA operator weight vector is obtained.
(b) The superior degree indicating that the ithalternative is more excellent than the other alternatives, is expressed as:
where is the ranked qth element in set { | j = 1, 2, ⋯, n}.
Equation 8 yields the ranking result of the ithalternative.
(c) The group consistency deviation index is defined as:
Where, c k corresponds to the authority of each expert. For simplification, c k = 1/k in this paper.
In order to thoroughly explore the sub-attribute characteristics (sub-attribute weight coefficient and expected interval) in fuzzy MAGDM, a new definition of group intelligence is proposed. Further, in order to reduce expert subjective fatigue in the decision making process, we designed group intelligence generating algorithms, which embeds learning metrics into the decision process. Furthermore, we provide the complex analysis of the algorithms.
(1) Group intelligence generating
where correspond to each sub-attribute’s expected interval vector (the interval’s upper and lower limits); λ ξ = = φ (λ1, λ2 … λ ξ ) corresponds to each sub-attribute’s priority vector (weight coefficient); and ξ is the number of sub-attributes.
where Equation (12) represents the Euclidean distance between all of the actual and upper limit values of the jth decision alternative’s sub-attributes, υ
ji
corresponds to the value of the jth decision alternative’s ith sub-attribute, υ
i
represents the decision alternative value corresponding to the article attributes.
where Equation (13) represents the Euclidean distance between all of the upper and lower limit values of the jth decision alternative’s sub-attributes, is the ith sub-attribute’s expected lower limit value, is the ith sub-attribute’s expected upper limit value, λ i corresponds to the ith sub-attribute’s weight coefficient, and i is the number of sub-attributes.
According to Equation (8), the group intelligence generating problems can be converted into a kind of optimization statement, and they can bedescribed by:
where Ij is the jth decision alternative’s satisfaction degree, d j is the superior degree index that is indicated in Equation (8), and j is the number of decision alternatives. The expected upper limit value and the expected lower limit value are unknown parameters, and i is the number of sub-attributes. Simultaneously, .
Key components of the group intelligence generating methods are described as follows: (a) Obtain every expert’s judgment preference (such as
The flowchart of group intelligence generating algorithm in regards to MAGDM problems is shown in Fig. 2.
Based on Fig. 2, parameters of group intelligence can be identified by GA. Therefore, the complexity analysis for the proposed group intelligence generating algorithms is demonstrated in detail.
1) Time complexity analysis:
The elements of vector λ ξ and are parameters to be optimized. Firstly, we specify the parameters of the GA as follows. t is described as the number of evolutionary generations (we set t = 10). χ (t) is designed as the population of the tth generation. ℓ denotes the population size (we set ℓ=40). q corresponds to the chromosome’s length, and the individual in the population was coded by the binary-encoded method. In addition, σ c is designed as the crossover rate (we set σ c = 0.9). σ m is designed as the mutation rate (we set σ m = 0.9). Ψ is designed as the fitness function (Equation (14)). Thus, the basic statements involved in the algorithms can be summarized as follows. In addition, π indicates the time performance.
Based on the GA, set the initializing parameters (σ
c
, σ
m
, t, Ψ); {t = 0; give the initiating population χ (t); // initiating population randomly (λ
ξ
and (ξ is the number of sub-attributes) in the Equation (10)); obtain the multi-objective fitness index; for (j = 1; j < t+1; j++) { execute the selection operations; execute the crossover operations; execute the mutation operations; achieve chromosome groups; // calculating the fitness index (Equation (12)) of the next generation; } }
According to statements ➂–➅, the time complexity of the algorithms is π [t × (ℓ × q)].
2) Space complexity analysis:
Space complexity refers to the cost of the storage space in computers when the algorithm is executing. In this sense, the space complexity of a group intelligence generating algorithm is π (ℓ × q).
Group intelligence-based MADM methods
Based on group intelligence generating algorithms, group intelligence-based decision making methods automatically calculate each overall value of alternatives through linear programming. The key step of the decision-making process is as follows:
(1) The decision-making matrix is expressed as:
Where, g ij is the ith sub-attribute value of the jth alternative.
(2) Based on harmonious parameters (Equation 10), as described in Equation 14, the standard matrix is written as:
Where, it is a normalized processing procedure and can be written as: .
(3) Thus, the comprehensive alternative value of z
j
is expressed as:
Where, λ i (i = 1, 2, ⋯ ξ) is described as the ith sub-attribute’s weight in λ ξ (Equation 10).
Test functions
The quality of our technique was measured by comparing the results from traditional fuzzy MAGDM approaches [19] to results from the proposed methods. Unlike the traditional fuzzy MAGDM approaches, the proposed method allows experts to only give their judgment preference for parts of the alternatives and have the computer calculate the ranking results of all the alternatives, in order to reduce the human’s subjectivity fatigue.
Consider a MAGDM problem that consists of four decision-makers and five decision-making alternatives, with each alternative containing three sub-attributes {δ1, δ2, ⋯ δ3 } tabulated out in Table 1.
The decision-makers wish obtain the best alternative whose each sub-attribute’s value is as large as possible. In addition, the expected sub-attribute values are (5,8,2). The group consistency deviation index is constrained by D (ω) ≤ 2.0 (Equation (9)).
In this experiment, in order to generate the group fuzzy preferences matrix, each of the 4 decision-makers are assigned a fuzzy complementary judgment matrix, represented as P1, P2, P3, and P4:
(a) Computing the OWA operator.
The OWA operator is obtained under the “majority” principle (0.3, 0.8). According to Equation (5), the OWA operator weight is obtained by:
(b) Gathering the group fuzzy preferences matrix.
According to Equations (6–10), the following are obtained: the superior degree index and the group consistency deviation index
As compared to the traditional MAGDM approaches, the proposed group intelligence-based decision methods show the convenience and flexibility of the decision process, the group intelligence generating methods towards the same decision problem above are described as follows: The interval of the unknown group intelligence’s parameters () are given, which is expressed as:
(b) Without requiring the decision-makers to evaluate all five alternatives, the fuzzy complementary preference matrices for the first three alternatives are provided by the four decision makers, which are expressed as:
(c) The superior degree is obtained based on Equation (8):
(d) According to Equation (10), the group consistency deviation index is obtained:
D = 0.2048 < 2.00
Additionally, the group intelligence generating problem can be into a kind of optimization statement (Equation 12).
GA are used to identify the parameters of group intelligence, such as each sub-attribute’s expected upper and lower limits of interval vector and sub-attribute’s weight vector λ ξ of the group intelligence in Equation (10), where ζ = 3.
The fitness function is:
where d j (j = 1, 2, 3) is shown in Equation (19).
There are nine unknown parameters (, , , , , , λ1, λ2, λ3]) that need to be identified by the GA. We set the specific parameters (10 generations and 40 individual species coded by the binary encoded method), and the initial individual species are set based on Equation (18) randomly. We then obtain the weight vector of the sub-attributes (0.22, 0.69, 0.09)
Figures 3–5 represent the parameter evolution curves corresponding to parameters λ ξ , , of the group intelligence in ten generations.
Finally, the group intelligence is obtained, expressed as:
Then, based on the proposed MAGDM methods shown in Equations 15–17 in Section 2.3, the satisfaction degree index is obtained for all five alternatives: (0.4256, 0.9832, 0.7372, 0.9516, 0.3032). Therefore, the ranking order of the alternative decisions is determined to be (2, 4, 3, 1, 5), thereby verifying the results obtained by traditional methods.
(1) Group Intelligence generating
In regards to a patient with a diagnosis based on the ‘mouth period’, a traditional Chinese swallowing dysfunction rehabilitation medicine problem is considered, consisting of eight expert rehabilitation doctors and ten decision-making alternatives, where each alternative contains five sub-attributes {δ1, δ2, ⋯ δ5 } that correspond to the treatment time (minutes) of physical therapy, muscle massage, medicament-induction, functional training, and acupuncture, respectively. Similarly, the group preference consistency deviation degree is constrained by D (ω) ≤ 2.0.
In addition, we give the optimization interval of the unknown group intelligence’s parameters (), which is expressed as:
Several decision-making alternatives are provided in Table 2.
(a) In regards to the ten alternatives, the Utility vectors for the eight decision makers are expressed as:
(b) Similarly, based on Equation (8), the superior degree index is obtained:
(c) According to Equation (10), the group consistency deviation index is obtained:
D = 1.02 < 2.00
Similarly, the group intelligence generating problem can be into a kind of optimization statement (Equation 12).
GA are used to identify the parameters of group intelligence, such as each sub-attribute’s expected upper and lower limits of interval vector and sub-attribute’s weight vector λ ξ of the group intelligence in Equation (10), where ζ = 5.
The fitness function is;
where d j (j = 1, 2, ⋯ 10) is shown in Equation (22).
There are fifteen unknown parameters (, , , , , , , , , , λ1, λ2, λ3, λ4, λ5]) that need to be identified by the GA. We set the specific parameters (10 generations and 40 individual species coded by the binary encoded method), and the initial individual species are set based on Equation (20) randomly. Figures 6–8 represent the parameter evolution curves corresponding to parameters λ. Vmax and vmin of the group intelligence in ten generations, respectively.
Finally, the group preferences are obtained, expressed as:
Therefore, based on the proposed MADM methods (presented by Equations 15–17 in Section 3.2), the ranking order of all ten alternatives is obtained: (3,10,1,4,9,5,6,2,7,8).
(2) Test of Medical Statistics
Thirty samples of the same symptom can be used as convincing statistical data in the medical statistics field [8]. Thirty samples with a diagnosis in the “mouth period” in the Beijing Zhongguancun Hospital physical rehabilitation department were used to verify the correctness of the proposed MADM method.
In this experiment, three expert doctors were employed to give their expert consultation solution based on every patient, while the computer will give its treatment solution based on the proposed MADM method.
Table 3 demonstrates the treatment strategies comparison between the proposed group intelligence-based decision making method and the experts’ consultation (which is expressed in the Appendix). Figure 9 illustrates a comparison of the time in minutes of the five treatment attributes (physical treatment, massage, etc.) between the proposed method and the expert consultation solution.
The consistency index [28] of the diagnosis and treatment is expressed as:
where φ′ is the result of the proposed MADM method, φ′ is the result of the expert consultation, |φ′ - φ|/φ is the relative error of the treatment solution, and β ∈ [0 1]. If β = 1, then the results of the proposed MADM method and the expert consultation are completely consistent.
Based on Table 3, it can be concluded that the average consistency index is 0.95 in the comparison between the proposed MADM method and the consultation of doctors.
In order to reflect the differences between the traditional fuzzy MAGDM and the group intelligence-based decision making methods proposed in this paper, a novel definition of group intelligence was presented. A Genetic Algorithm (GA) was employed to realize the interactive generating algorithms for the group intelligence.
Based on the experiment described in Section 3.1, the traditional fuzzy MAGDM methods sacrifice flexibility and analysis speed by requiring a judgment preference for each alternative. In contrast to traditional approaches, our proposed decision method utilizes only the first three alternatives to generate the group preferences (Figs. 3–5), which can be used to determine the ordering of the remaining two alternatives. With the same decision results, the proposed decision methods (Section 2.4) based on group intelligence generating are capable of gradually learning the group’s judgment preference in decision-making and reducing their subjective fatigue.
In addition, a traditional Chinese medicine rehabilitation treatment problem was designed in Section 3.2. To overcome the traditional Chinese rehabilitation medicine example, wherein expert doctors’ treatment preferences lack quantitative descriptions, the proposed methods were used to obtain a treatment weight and time (Figs. 6–9) of different TCM techniques (i.e., physical therapy, massage, and acupuncture) for a patient with dysphagia.
In the end, the proposed decision methods based on the group intelligence did not decrease the decision accuracy.
Conclusions
In order to provide beneficial methods for determining group intelligence in MAGDM, innovative definitions of group intelligence based on the weight coefficient and the expected interval of each sub-attribute were proposed. The group intelligence generating algorithms were explicitly proposed, which embedded learning metrics into the decision process, thus helping reduce subjective fatigue in the humans. In addition, a complexity analysis was presented that helped ensure the accuracy of the MAGDM results. In the end, a test function and a kind of medical MAGDM problem addressing Traditional Chinese medical (TCM) dysphagia rehabilitation treatment were used to verify the correctness and feasibility of the proposed methods.
Nonetheless, it should be acknowledged that this research is rudimentary, and as such, there are some challenges for the algorithm proposed in this paper. For example, (1) a large number of the parameters should be optimized when the number of attributes is excessive; thus, a heavy searching burden is suffered. Therefore, in future work, machine learning methods, such as Principal Component Analysis (PCA) and Deep learning (DL), and so on, can be used to extract the sub-attribute’s feature to optimize the generating process of group intelligence. Additionally, (2) it is very boring that each expert must give his judgment preference, such as a utility vector and so on, in the decision process. Deep research can be conducted on the ‘human’s affective—human’s preference’ mapping approaches based on affect computing methods, and the harmonious of the interaction can be improved.
Footnotes
Appendix
Treatment strategies comparison between the proposed method and the experts’ consultation
| Patient with | The group intelligence based decision making methods | Consultation of expert doctors | ||||||||
| diagnosis | Physical | Muscle | Medicament- | Functional | Acupuncture | Physical | Muscle | Medicament- | Functional | Acupuncture |
| ‘mouth- | therapy | massage | induction | training | (min) | therapy | massage | induction | training | (min) |
| period’ | (min) | (min) | (min) | (min) | (min) | (min) | (min) | (min) | ||
| 1 | 12 | 11 | 16 | 19 | 19 | 12 | 12 | 16 | 18 | 19 |
| 2 | 12 | 12 | 17 | 18 | 18 | 12 | 12 | 17 | 19 | 18 |
| 3 | 10 | 13 | 15 | 20 | 20 | 11 | 13 | 16 | 19 | 20 |
| 4 | 11 | 12 | 18 | 15 | 19 | 11 | 13 | 19 | 16 | 17 |
| 5 | 12 | 14 | 17 | 20 | 18 | 12 | 12 | 17 | 20 | 16 |
| 6 | 9.5 | 12 | 17 | 18 | 18 | 10 | 11 | 17 | 18 | 16 |
| 7 | 8.5 | 15 | 17 | 16.5 | 18 | 9 | 13 | 16 | 16 | 18 |
| 8 | 9.5 | 14 | 15 | 19.5 | 16 | 10 | 12 | 14 | 20 | 16 |
| 9 | 9 | 13 | 16 | 16 | 16 | 9 | 13 | 15 | 16 | 15 |
| 10 | 11.5 | 15 | 16 | 18 | 18 | 11 | 14 | 15 | 18 | 16 |
| 11 | 10.5 | 12 | 19 | 18 | 20 | 10 | 12 | 12 | 18 | 20 |
| 12 | 13 | 10 | 21 | 18 | 20 | 13 | 10 | 20 | 18 | 20 |
| 13 | 12.5 | 11 | 20 | 15.5 | 22 | 13 | 10 | 20 | 15 | 22 |
| 14 | 12.5 | 11 | 19 | 17 | 18 | 13 | 10 | 20 | 16 | 16 |
| 15 | 11 | 14 | 17 | 17 | 18 | 11 | 14 | 15 | 16 | 19 |
| 16 | 12 | 13 | 16 | 20 | 16 | 12 | 12 | 15 | 20 | 16 |
| 17 | 11 | 16 | 19 | 21.5 | 18 | 10 | 16 | 18 | 20 | 18 |
| 18 | 10.5 | 15 | 19 | 21 | 18 | 11 | 14 | 20 | 20 | 18 |
| 19 | 9.5 | 16 | 20 | 19.5 | 20 | 10 | 16 | 20 | 20 | 18 |
| 20 | 12 | 14 | 17 | 18 | 16 | 12 | 14 | 18 | 18 | 15 |
| 21 | 9 | 15 | 15 | 17 | 15 | 9 | 16 | 16 | 18 | 15 |
| 22 | 11 | 14 | 17 | 18 | 16 | 11 | 14 | 18 | 19 | 15 |
| 23 | 9.5 | 12 | 17 | 18 | 18 | 10 | 12 | 18 | 18 | 18 |
| 24 | 10 | 11 | 15 | 20.5 | 16 | 10 | 10 | 15 | 20 | 15 |
| 25 | 12.5 | 16 | 15 | 18 | 16 | 12 | 15 | 15 | 18 | 15 |
| 26 | 11.5 | 11 | 21 | 16.5 | 20 | 11 | 10 | 20 | 18 | 20 |
| 27 | 11 | 15 | 20 | 17.5 | 22 | 11 | 14 | 20 | 18 | 22 |
| 28 | 10 | 12 | 21 | 18 | 20 | 10 | 12 | 20 | 18 | 20 |
| 29 | 12.5 | 12 | 18 | 15.5 | 18 | 13 | 12 | 18 | 16 | 18 |
| 30 | 11 | 15 | 19 | 16 | 18 | 11 | 16 | 20 | 16 | 18 |
Acknowledgments
This study is supported by the Open Research Project under Grant from SKLMCCS (20150103), the Beijing outstanding talent training Project (2015000020124G041), the Fundamental Research Funds for the Central Universities (YS1404) and the National Natural Science Foundation of China (61603023).
The authors Jie Chen and Chong Su contributed equally to this work and should be considered co-first authors.
