Clinical Knowledge Supported Acute Kidney Injury (AKI) Risk Assessment Model for Elderly Patients

Abstract

From the clinical viewpoint, the statistical approach is still the cornerstone for exploring many diseases. This study was conducted to explore the risk factors related to acute kidney injury (AKI) for elderly patients using the multiple criteria decision-making (MCDM) approach. Ten nephrologists from a teaching hospital in Taipei took part in forming the AKI risk assessment model. The key findings are: (1) Comorbidity and Laboratory Values would influence Comprehensive Geriatric Assessment; (2) Frailty is the highest influential AKI risk factor for elderly patients; and (3) Elderly patients could enhance their daily activities and nutrition to improve frailty and lower AKI risk. Furthermore, we illustrate how to apply MCDM methods to retrieve clinical experience from seasoned doctors, which may serve as a knowledge-based system to support clinical prognoses. In conclusion, this study has shed light on integrating multiple research approaches to assist medical decision-making in clinical practice.

1. Introduction

There are roughly two types of clinical research methods in conventional medical research: observation-based statistical testing and meta-analysis [1]. The former has to collect over one group of data, after controlling for demographic factors and the associated variables (e.g., symptoms) of patients, to examine the premises or consequences of a targeted illness (or disease). The latter stands on the grounds of published homogeneous research to explore hypothesized patterns from the papers that satisfy the inclusion criteria [1,2]. They both belong to the evidence-based statistical approach, and we focus on enhancing the findings of a meta-analysis in this work.

Though the statistical approach has been broadly adopted, certain limitations are rooted in the presumptions of statistics [3]. For instance, most statistical models assume variables to be independent, which is controversial. Besides, researchers have to presume the probabilistic distribution(s) of involved variables. It is difficult to validate those assumptions [4]; those limitations impede researchers to explore complicated relations in their models. Furthermore, there is nearly no room to accommodate the clinical experiences/knowledge from doctors. The research gap between the findings of meta-analysis and clinical prognoses exists.

In clinical research, the statistical approach is still the cornerstone to understanding many illnesses. A well-devised meta-analysis may explore and examine inconclusive disease patterns—with rigorous control and testing—from credible sources. Its contributions are evident. However, owing to the limitations of its modeling methods (i.e., statistics), meta-analysis is insufficient for clarifying the influential relationship among a group of factors without presuming the probabilistic distributions of its variables. In most cases, meta-analysis can merely identify the significance of one or more factors that might be associated with a disease. Therefore, to enrich a meta-analysis’ findings and bridge the gap, this study proposes the multiple criteria decision-making (MCDM) approach [5,6,7] to incorporate clinical knowledge from seasoned doctors.

Unlike traditional statistics, the MCDM approach can support decision-makers (DMs) or experts to make rational judgments by considering multiple factors simultaneously, without presumed probabilistic distributions. This approach hinges upon domain experts’ knowledge/experience. Certain MCDM methods may unveil the convoluted interrelationship of a complex model, such as Decision-Making Trial and Evaluation Laboratory (DEMATEL) [8,9], Analytical Network Process (ANP) [10], and the DEMATEL-based ANP (DANP) [11,12]. The advantages of MCDM methods complement the limitations of meta-analyses; the outcome may even serve as a knowledge-based system to support clinical prognoses.

The present study adopts acute kidney injury (AKI) [13] as an empirical case, which has a high risk of mortality for the elderly. By definition, “AKI is an abrupt decrease in kidney function that includes, but is not limited to, acute renal failure. [14]” Nowadays, AKI is becoming common among elderly patients, who are more vulnerable to kidney failure or dysfunction [15]. In this context, our research purposes are twofold: (1) elaborate the findings from a published meta-analysis based on the clinical experience of doctors and (2) construct a flexible knowledge-based system to support clinical prognoses. The specific aims of this case study are to:

(1)

explore the interrelationships between the risk factors that might induce AKI in elderly patients;

(2)

identify the relative influence of each criterion on the risk assessment of AKI;

(3)

provide precautionary guidance for elderly patients to reduce their potential AKI. risk.

Aside from the plausible advantages, the proposed model has the flexibility to adjust the weighting of doctors’ opinions based on their experience (or expertise in a specific field). This flexibility is yielded from the group decision theory [16]. This advantage contributes to carrying out the policy of a hospital. Evaluating a patient with complicated symptoms or comorbidity usually requires doctors’ participation with different experience levels and varied expertise. The hybrid MCDM approach can fuse experts’ knowledge with a customized weighting policy. This is another contribution that the present study might bring to clinical applications. We will illustrate this idea by an example in Section 4.

This paper is divided into five sections. Section 1 discusses the background information as above. Section 2 provides the research material and the proposed hybrid MCDM model. Section 3 reports the modeling of AKI risk assessments for elderly patients. Section 4 provides discussions with an illustration of how to apply the model and Section 5 concludes this study.

2. Materials and Research Methods

This section explains the materials adopted in modeling. Since the present study aims to enrich meta-analysis findings, a study of AKI in the elderly population from a recently published paper [15] serves as the modeling foundation. Therefore, in Section 3.1, the findings of [15] are briefly highlighted. Besides this, ten experienced doctors took part in this work. Section 2.2 introduces two MCDM methods: the DEMATEL technique [8,9] and the DEMATEL-based ANP (DANP) [11,12]. The first one can indicate the influential relationship of a model with multiple factors (criteria). The second one may explore the relative influence of each factor on supporting clinical prognoses.

2.1. Materials

The meta-analysis article [15] that we referred to mainly discussed whether elderly patients might suffer from a higher risk of AKI because of frailty. It began with 1385 papers that are included in PubMed (566), the Cochrane Library (609), and EMBASE (210). After following the PRISMA (Preferred Reporting Items for Systematic Review and Meta Analysis) research flow [17], it excluded duplicated, unrelated or unsuitable articles. It eventually identified four cohort studies [18,19,20,21] that met the inclusion criteria. Its rigorous research process [15] attempted to overcome the potential issues of meta-analysis (such as low quality of included studies and heterogeneity [22]), which serves as the cornerstone of the present study.

Owing to its research target, the meta-analysis [15] attempted to include studies that controlled for the variable, diagnosed AKI or non-AKI, for those elderly with frailty. The selection process examined the inter-rater reliability and evaluated the degree of variation across studies; both passed the associated statistical thresholds (i.e., the kappa and ${\mathrm{I}}^2$ statistics). Among the four cohort studies, the study duration ranged from two weeks [20] to four years [18]. The total subjects included in each study were between 164 [20] to 533 [19].

Though the four studies adopted different frailty measures (e.g., Clinical Frailty Scores), their assessments on frailty are all grounded on clinically recognized measurements. The meta-analysis from the four papers confirmed the significant association between frailty and AKI among elderly patients.

2.2. Research Methods

In this subsection, we introduce the essentials of the DEMATEL and DEMATEL-based ANP (DANP). DANP is extended from the DEMATEL, which adopts the central idea of ANP to derive the relative influence of each criterion of a hybrid decision model. The required research flows are illustrated in Figure 1.

2.2.1. DEMATEL Technique

The DEMATEL technique was devised to resolve complicated issues, and assumes that all the included factors (criteria) are interrelated. Furthermore, it relies on domain experts’ knowledge (opinions) to derive the influential relationship among a model’s factors. The essential steps to conduct the DEMATEL analysis are as follows.

Step 1: Collect opinions from domain experts to form an initial matrix $A^{\mathrm{I}}$.

As mentioned above, the DEMATEL technique presumes that all the included criteria are interrelated. Therefore, it begins by devising a questionnaire to have all the criteria with questions like: “in your opinion, what is the influence of criterion i on criterion j (the answer ranges from “0 (No Influence),” “1 (Minor Influence),” “2 (Moderate Influence)” to “3 (High Influence),” to “4 (Very High Influence)).” After collecting opinions from domain experts, the averaged result of the influence of criterion i on criterion j will be indicated as $\overline{a_{ij}}$ in $A^{\mathrm{I}}$ ($A^{\mathrm{I}}={\left[\overline{a_{ij}}\right]}_{n\times n}$, for i,j = 1,…,n), located on the i-th row and the j-th column of $A^{\mathrm{I}}$.

Step 2: Normalize $A^{\mathrm{I}}$ to become a matrix $N$.

The normalized matrix $N$ can be obtained by calculating k, where $k=\max \left\{\max {\displaystyle \sum _{j=1}^n\overline{a_{ij}}},\max {\displaystyle \sum _{i=1}^n\overline{a_{ij}}}\right\}$, and the normalized matrix $N=k\times A^{\mathrm{I}}$.

Step 3: Calculate the total-influence matrix $T$.

The total-influence matrix $T$ is defined as $T=N+N^2+\dots +N^w$, which indicates the chain influence by raising the number of w. It can be decomposed as Equation (1):

$$T=N\times (I-N^w)\times {(I-N)}^{-1}$$

(1)

while $w\to \infty$, $N^w\cong {[0]}_{n\times n}$. Therefore, the total-influence matrix $T$ can be approximated and indicated as: $T=N\times {(I-N)}^{-1}={\left[t_{ij}\right]}_{{}_{n\times n}}$. The sum of each row and column of $T$ can form two vectors, termed as $r^C={\left(r_1,\cdots ,r_i,\cdots ,r_n\right)}^{\prime }$ and $d^C=\left(d_1,\cdots ,d_j,\cdots ,d_n\right)$. In the DEMATEL analytic, “$r^C+d^C$” and “$r^C-d^C$” can be adopted to indicate the directional influence relationship among criteria. The details will be explained in Section 4.

Step 4: Form initial unweighted super-matrix $T_C^S$.

Assume that there are m dimensions in $T$. Then, the total-influence matrix $T$ can be partitioned by its dimensions and indicated as an initial unweighted super-matrix $T_C^S$ that comprises $m\times m$ sub-matrices in Equation (2).

$$T_C^S=\begin{array}{cc}& \begin{array}{ccccc}{D}_1 & & D_j & & D_m \end{array}\\ \begin{array}{c}{D}_1 \\ ⋮ \\ D_i \\ ⋮ \\ D_m\end{array}& {\left[\begin{array}{ccccc}{T}_C^{S11}& \cdots & T_C^{S1j}& \cdots & T_C^{S1m}\\ ⋮& & ⋮& & ⋮\\ T_C^{Si1}& \cdots & T_C^{Sij}& \cdots & T_C^{Sim}\\ ⋮& & ⋮& & ⋮\\ T_C^{Sm1}& \cdots & T_C^{Smj}& \cdots & T_C^{Smm}\end{array}\right]}_{n\times n}\end{array},\mathrm{for}1<i,j\le m<n.$$

(2)

For instance, if there were three criteria in $D_1$, the sub-matrix $T_C^{S11}$ would be a $3\times 3$ matrix. $T_C^{S11}$ denotes the corresponding upper-left $3\times 3$ elements of $T_C^S$.

Step 5: Form the dimension matrix $T_D^{}$.

The averaged figures of all elements of each sub-matrix of $T_C^S$ will form a $m\times m$ dimension matrix, which is shown in Equation (3):

$$T_D={\left[\begin{array}{ccccc}\overline{t_C^{S11}}& \cdots & \overline{t_C^{S1j}}& \cdots & \overline{t_C^{S1m}}\\ ⋮& & ⋮& & ⋮\\ \overline{t_C^{Si1}}& \cdots & \overline{t_C^{Sij}}& \cdots & \overline{t_C^{Sim}}\\ ⋮& & ⋮& & ⋮\\ \overline{t_C^{Sm1}}& \cdots & \overline{t_C^{Smj}}& \cdots & \overline{t_C^{Smm}}\end{array}\right]}_{m\times m},\mathrm{for}1<i,j\le m.$$

(3)

Similarly, the sum of each row and column of $T_D^{\mathrm{N}}$ can form two vectors, termed as $r^D={\left(r_1^D,\cdots ,r_i{}^D,\cdots ,r_m^D\right)}^{\prime }$ and $d^D=\left(d_1{}^D,\cdots ,d_j^D,\cdots ,d_m^D\right)$. Then, “$r^D+d^D$” and “$r^D-d^D$” can be adopted to indicate the directional influence relationship among the dimensions. Based on the DEMATEL analytics from Steps 3 and 5, an Internetwork Relationship Map (INRM) can help depict the directional influences among dimensions and criteria.

2.2.2. DEMATEL-Based ANP (DANP)

The DANP method adopts the central idea of the ANP and the transpose of $T_C^S$ (i.e., $W={\left(T_C^S\right)}^{\prime }$) is regarded as an unweighted supermatrix in the ANP. Compared with the ANP’s assumption, the DANP considers the dimensional influences while forming an initial supermatrix.

Step 6: Obtain the initial supermatix of DANP.

First, the dimension matrix $T_D^{}$ from Step 5 should be normalized to become $T_D^{\mathrm{N}}$ (refer to Equation (4)). Then, the initial supermatrix of DANP can be obtained by using $T_D^{\mathrm{N}}$ to adjust the unweighted supermatrix $W$(i.e., $W^{adj}=T_D^{\mathrm{N}}\times W$).

$$T_D^{\mathrm{N}}={\left[\begin{array}{ccccc}\overline{t_C^{S11}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{S1j}}}& \cdots & \overline{t_C^{S1j}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{S1j}}}& \cdots & \overline{t_C^{S1m}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{S1j}}}\\ ⋮& & ⋮& & ⋮\\ \overline{t_C^{Si1}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{Sij}}}& \cdots & \overline{t_C^{Sij}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{Sij}}}& \cdots & \overline{t_C^{Sim}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{Sij}}}\\ ⋮& & ⋮& & ⋮\\ \overline{t_C^{Sm1}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{Smj}}}& \cdots & \overline{t_C^{Smj}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{Smj}}}& \cdots & \overline{t_C^{Smm}}/{\displaystyle \sum _{j=1}^m\overline{t_C^{Smj}}}\end{array}\right]}_{m\times m},\mathrm{for}1<i,j\le m.$$

(4)

Once the adjusted supermatrix $W^{adj}$ is obtained, the influential weights can be reached by raising the number of q in $\underset{q\to \infty }{\lim }{\left(W^{Adj}\right)}^q$ until $W^{adj}$ becomes stable. After normalizing the stable weights, the final DANP influential weights can be obtained.

3. AKI Risk Assessment Model for Elderly Patients

In the present study, the definition of AKI is based on Kidney Disease: Improving Global Outcomes (KDIGO), and the criterion for elderly patients is aged ≥65 years. Due to AKI being a complex disorder related to the interplay of patient-associated factors within the environment where patients live, the dimensions of heterogeneous factors should consider environmental exposures and the inherent risk elements [23]. In each dimension, risk factors could be identified to design and implement preventive interventions or further investigations.

In this context, the preliminary version of the AKI risk assessment model is based on the meta-analysis article [15] and related research (i.e., [18,19,20,21], summarized in

Table 1. Summary of the four studies that are referred to.

Ref	Country	Database Source	Research Subjects	Mean Age of Research Subjects	Duration	Key Findings
[18]	USA	Pubmed	243 AKI 74 non-AKI	57 (AKI) 56 (non-AKI)	2007–2010 (4 years)	AKI was associated with higher frailty (Clinical Frailty Scale).
[19]	Korea	Pubmed	183 (mild frail) 199 (moderate frail) 152 (severe frail)	73.8 (mild frail) 76.5 (moderate frail) 79.0 (severe frail)	2013 (1 year)	The frailest group indicated the highest AKI risk compared with the other groups.
[20]	UK	Pubmed	31 AKI 133 non-AKI	82.6 (AKI) 81.2 (non-AKI)	2 weeks	Severe frailty was associated with AKI significantly (p = 0.01).
[21]	USA	Pubmed	93 (non-frail) 139 (pre-frail) 136 (frail)	68.52 (non-frail) 74.71 (pre-frail) 78.83 (frail)	2013–2014 (2 years)	Frail patients were more likely to induce AKI (p = 0.03).

), and the three dimensions all come from [19] after removing demographic factors. One thing that needs to be mentioned here is that Frailty, (C₃₃), is the primary factor highlighted by [15].

Next, we further consulted on the preliminary version with ten nephrologists from one teaching hospital in Taipei. Based on their clinical experience, they suggested including another two criteria: “Depression (C₁₃)” and “Hepatitis B/C (C₁₅)” for AKI elderly patients. After several rounds of discussions with the doctors, we finalized the 14 criteria in three dimensions, shown in

Table 2. The dimensions and criteria of the AKI (acute kidney injury) risk assessment decision model.

Dimensions	Criteria	Reference
	Diabetes (C11)	[19,20]
	Hypertension (C12)	[19,20]
Comorbidity	Depression (C13)	(from the doctors)
(D1)	Malignancy (C14)	[19,20]
	Hepatitis B/C (C15)	(from the doctors)
	Creatinine (C21)	[19,20]
	Estimated glomerular filtration rate (C22)	[19]
Laboratory Values	Hemoglobin (C23)	[19]
(D2)	Albumin (C24)	[19]
	Na (C25)	[19]
Comprehensive	Activities of daily living (C31)	[19]
Geriatric	Mid-arm circumference (C32)	[19]
Assessments	Frailty (C33)	[15,18,19,20,21]
(D3)	Nutritional assessment (C34)	[19]

.

This study’s purposive sampling was conducted in one medical center in northern Taiwan from April 2020 to June 2020. We recruited well-trained nephrologists who were able to express willingness to participate in the survey.

In Table 2, the three dimensions all come from [19] after removing demographic factors. The estimated glomerular filtration rate (eGFR) is based on the CKD epidemiology collaboration formula [24].

3.1. DEMATEL Analysis

From Steps 1 to 4 in the previous section, we collected the ten doctors’ opinions. In this case, to mitigate the potential risk of outliers that might distort the average result, this study adopted the median value, instead of the arithmetic mean, to form the initial average matrix $A^{\mathrm{I}}$ (

Table 3. Initial average matrix $A^{\mathrm{I}}$.

	C11	C12	C13	C14	C15	C21	C22	C23	C24	C25	C31	C32	C33	C34
C11	0.00	4.00	3.00	3.00	2.00	3.00	3.00	3.00	2.00	2.00	3.00	3.00	3.00	2.00
C12	2.00	0.00	1.00	2.00	1.00	4.00	4.00	2.00	3.00	3.00	2.00	2.00	2.00	2.00
C13	2.00	2.00	0.00	2.00	2.00	2.00	2.00	1.00	2.00	1.00	4.00	3.00	3.50	3.00
C14	2.00	2.00	4.00	0.00	2.50	3.00	3.00	4.00	3.00	3.00	4.00	3.50	4.00	4.00
C15	2.00	1.00	2.00	4.00	0.00	3.00	3.00	3.00	4.00	2.00	4.00	4.00	4.00	4.00
C21	1.00	2.50	1.00	1.00	1.00	0.00	4.00	4.00	3.00	4.00	3.50	3.50	4.00	4.00
C22	1.00	2.50	2.00	3.00	1.50	4.00	0.00	4.00	3.50	4.00	4.00	3.50	4.00	4.00
C23	1.00	1.50	1.00	1.00	1.00	2.00	2.00	0.00	3.00	1.00	3.00	3.00	3.00	2.50
C24	1.00	1.00	1.00	1.00	2.00	1.00	1.00	1.00	0.00	2.00	3.00	3.50	3.00	4.00
C25	1.00	3.00	1.00	1.00	1.00	3.00	3.00	1.00	1.00	0.00	2.00	1.00	2.00	2.00
C31	4.00	4.00	4.00	3.50	4.00	2.00	2.00	2.00	2.00	2.00	0.00	4.00	4.00	4.00
C32	2.00	2.00	2.00	2.00	1.50	2.00	2.00	2.00	2.00	2.00	2.00	0.00	4.00	4.00
C33	2.00	2.00	4.00	3.00	2.00	2.00	2.00	2.00	2.00	2.00	4.00	4.00	0.00	4.00
C34	3.00	3.00	3.00	3.00	3.00	3.00	3.00	3.00	3.00	3.00	4.00	4.00	4.00	0.00

).

Since there are three dimensions (i.e., $m=3$) in this analysis, the dimension matrix $T_D^{}$ is a $3\times 3$ matrix, shown in Appendix A (

Table A1. The dimension matrix $T_D^{}$.

	D1	D2	D3
D1	0.197	0.239	0.307
D2	0.163	0.196	0.266
D3	0.226	0.237	0.309

). The associated directional influence analytics for the dimensions and criteria are summarized in

Table 4. Directional influence analytics from DEMATEL (Decision-Making Trial and Evaluation Laboratory) analysis.

Dimensions	${\mathit{r}}^{\mathit{D}}$	${\mathit{d}}^{\mathit{D}}$	${\mathit{r}}^{\mathit{D}}+{\mathit{d}}^{\mathit{D}}$	${\mathit{r}}^{\mathit{D}}-{\mathit{d}}^{\mathit{D}}$	Criteria	${\mathit{r}}^{\mathit{C}}$	${\mathit{d}}^{\mathit{C}}$	${\mathit{r}}^{\mathit{C}}+{\mathit{d}}^{\mathit{C}}$	${\mathit{r}}^{\mathit{C}}-{\mathit{d}}^{\mathit{C}}$
					C11	3.44	2.40	5.84	1.03
					C12	2.87	2.96	5.83	−0.09
D1	0.74	0.59	1.33	0.16	C13	2.92	2.87	5.80	0.05
					C14	3.97	2.86	6.83	1.11
					C15	3.85	2.44	6.29	1.41
					C21	3.43	3.21	6.63	0.22
					C22	3.85	3.21	7.06	0.64
D2	0.63	0.67	1.30	−0.05	C23	2.43	3.05	5.48	−0.62
					C24	2.42	3.17	5.59	−0.76
					C25	2.18	2.99	5.17	−0.81
					C31	3.98	4.01	8.00	−0.03
D3	0.77	0.88	1.66	−0.11	C32	2.87	4.02	6.89	−1.15
					C33	3.39	4.22	7.61	−0.83
					C34	3.98	4.16	8.14	−0.18

and illustrated in Figure 2.

The opinions that comprise $A^{\mathrm{I}}$ (i.e., Table 3) were examined by Equation (5). The confidence level (${\mathrm{CL}}_{A^{\mathrm{I}}}$) for the result is higher than 95% and k is the number of questionnaires (k = 10 in here). In other words, the initial average matrix carries enough confidence to indicate the averaged opinions to forming $A^{\mathrm{I}}$.

$${\mathrm{CL}}_{A^{\mathrm{I}}}=\frac{1}{n(n-1)}{\displaystyle \sum _{i=1}^k{\displaystyle \sum _{j=1}^k\frac{\left|\overline{a_{ij}^k}-\overline{a_{ij}^{k-1}}\right|}{\overline{a_{ij}^k}}\times 100\%}}$$

(5)

From Table 4, D₁ (Comorbidity) would influence D₂ (Laboratory Values) and D₃ (Comprehensive Geriatric Assessments). It is worthwhile to delve into D₃, which indicates that Mid-arm circumference (C₃₂) and Frailty (C₃₃) are influenced by Activities of daily living (C₃₁) and Nutritional assessment (C₃₄).

3.2. DEMATEL-Based ANP (DANP) Weights

Based on the initial average matrix $A^{\mathrm{I}}$ in Table 3 and the six steps in Section 2, we obtained the DANP influential weights for the AKI risk assessment model. The normalized dimensional matrix $T_D^{\mathrm{N}}$ is shown in

Table 5. Normalized dimension matrix $T_D^{\mathrm{N}}$.

	D1	D2	D3
D1	0.27	0.32	0.41
D2	0.26	0.31	0.43
D3	0.29	0.31	0.40

, and the associated influential weight of each criterion is summarized in Table 5. The DEMATEL-adjusted initial supermatrix is in Appendix A (

Table A2. DEMATEL-adjusted initial supermatrix $W^{adj}$.

	C11	C12	C13	C14	C15	C21	C22	C23	C24	C25	C31	C32	C33	C34
C11	0.036	0.053	0.051	0.048	0.049	0.046	0.043	0.047	0.046	0.044	0.055	0.054	0.051	0.054
C12	0.065	0.049	0.059	0.056	0.053	0.064	0.059	0.059	0.054	0.070	0.062	0.063	0.059	0.062
C13	0.059	0.054	0.045	0.065	0.057	0.053	0.055	0.054	0.053	0.051	0.062	0.062	0.069	0.061
C14	0.059	0.061	0.058	0.045	0.066	0.053	0.059	0.054	0.053	0.051	0.059	0.062	0.063	0.061
C15	0.047	0.047	0.052	0.051	0.040	0.046	0.046	0.047	0.054	0.045	0.055	0.051	0.051	0.054
C21	0.068	0.070	0.068	0.065	0.065	0.048	0.068	0.067	0.064	0.076	0.063	0.063	0.063	0.063
C22	0.068	0.070	0.068	0.065	0.065	0.069	0.049	0.067	0.064	0.076	0.063	0.063	0.063	0.063
C23	0.065	0.056	0.059	0.067	0.064	0.067	0.066	0.048	0.061	0.056	0.060	0.060	0.060	0.060
C24	0.062	0.064	0.068	0.064	0.070	0.063	0.065	0.075	0.054	0.058	0.062	0.062	0.062	0.062
C25	0.059	0.062	0.058	0.061	0.058	0.066	0.065	0.056	0.069	0.047	0.059	0.059	0.059	0.059
C31	0.103	0.101	0.106	0.102	0.101	0.103	0.105	0.105	0.101	0.107	0.080	0.095	0.105	0.104
C32	0.103	0.101	0.099	0.099	0.102	0.103	0.102	0.106	0.106	0.096	0.105	0.078	0.105	0.105
C33	0.107	0.106	0.106	0.106	0.105	0.110	0.110	0.110	0.105	0.112	0.109	0.114	0.082	0.109
C34	0.099	0.105	0.102	0.105	0.104	0.109	0.109	0.104	0.113	0.111	0.107	0.113	0.108	0.083

).

In

Table 6. DANP (DEMATEL-based Analytical Network Process) influential weights for the AKI risk assessment model for elderly patients.

Dimensions	Dimensional Weights	Criteria	DANP Weights (%)
		C11	4.92
		C12	6.00
Comorbidity	27.53%	C13	5.83
(D1)		C14	5.80
		C15	4.98
		C21	6.45
Laboratory		C22	6.46
Values	31.31%	C23	6.06
(D2)		C24	6.34
		C25	5.99
Comprehensive		C31	10.05
Geriatric	41.17%	C32	10.04
Assessments		C33	10.59 (1st)
(D3)		C34	10.48 (2nd)

, D₃ (Comprehensive Geriatric Assessments) shows the highest influential weight (i.e., 41.17%), and Frailty (C₃₃) is the most influential criterion. This finding echoes the conclusion of [15]. Though the result is similar to the previous research [15,18,20], the DEMATEL analysis further identifies the source factors (dimensions/criteria) of AKI for the elderly patients. In Figure 1, Comorbidity (D₁) and Laboratory Values (D₂) both would influence D₃ (Comprehensive Geriatric Assessments). Besides, within D₃ (Comprehensive Geriatric Assessments), the DEMATEL analysis suggests that C₃₁ (Activities of daily living) and C₃₄(Nutritional assessment) may influence Frailty (C₃₃).

In other words, to mitigate the AKI risk level of an elderly patient, adequate enhancement of daily activities and nutrition might offer another type of treatment; those two criteria would improve an elderly patient’s overall frailty. This finding is complementary to the result of [15]. Under certain controllable circumstances, elderly patients may improve their frailty before surgery that might cause AKI complications. Thus, the DANP model may help improve patients’ awareness of daily activities and nutrition in practice.

4. Discussions

Healthcare decisions are complex and often involve trade-offs among multiple conflicting objectives. Using structured and explicit approaches to medical decisions involving multiple criteria could improve decision-making using a set of techniques. MCDM methods are widely used in other sectors and there has been an increase in healthcare applications [25].

Clinical examples of MCDM applications include prioritizing patients for non-urgent surgery [26], disease diagnosis and classification [27], antibiotic-resistant diseases for R&D [28], supporting patients and physicians in selecting treatments [29], and weighing up the benefits and risks of new medicines to support licensing decisions [30]. From the clinical viewpoint, MCDM models are devised to serve as decision aids. After referring to an MCDM model’s outcome, DMs could deliberate on which might be the most beneficial option before making their final decision.

In Table 6, we transformed ten doctors’ opinions into a DANP-based AKI risk assessment model. Once an elderly patient was under evaluation, he/she could be graded from 1 (Low Risk) to 10 (High Risk) on each criterion by doctors. The overall AKI risk level of an elderly patient can be aggregated by adopting the DANP weights (in Table 5), also called the simple additive weighting (SAW) method. This model may be deployed to accumulate clinical prognoses to associate the overall AKI risk level with doctors’ evaluations. For instance, if there were three elderly patients (i.e., A, B, and C) who need to be assessed, their associated grades on each criterion are assumed in

Table 7. DANP-based AKI risk assessment model with three assumed cases.

DANP			C11	C12	C13	C14	C15	C21	C22	C23	C24	C25	C31	C32	C33	C34	Risk Level
Weights (%)			4.92	6.00	5.83	5.80	4.98	6.45	6.46	6.06	6.34	5.99	10.05	10.04	10.59	10.48
	A		9.00	8.00	8.00	2.00	1.00	6.00	5.00	4.00	6.00	7.00	8.00	8.00	9.00	6.00	6.50
	B		2.00	7.00	8.00	2.00	1.00	4.00	3.00	4.00	5.00	7.00	5.00	8.00	7.00	5.00	5.16
	C		7.00	6.00	2.00	2.00	1.00	5.00	2.00	3.00	2.00	4.00	7.00	5.00	5.00	5.00	4.25

.

Table 7 shows that the three patients’ AKI risk levels are ranked as: $A\succ B\succ C$ (i.e., A has the highest AKI risk), by using the DANP-based decision model. Though the model can be applied to show the ranking of patients’ AKI risk level, it is insufficient to suggest a threshold that could be regarded as dangerous. However, suppose a hospital may adopt this approach to accumulate records and doctors’ prognoses; it may collect cases for further analytics. In that case, this approach might conclude with a meaningful threshold or interval as a warning indicator. For instance, if only case A (in Table 7) was diagnosed as dangerous by doctors, the model may keep one record categorizing the risk level 6.50 as dangerous. This approach provides a systematic way to collect AKI-related measures for elderly patients.

Based on the group-decision theory [16], the current model can be adjusted using doctors’ experience to re-calculate the DANP weights. Take the ten doctors involved in this study as an example; they can be grouped into three categories based on their clinical experience: (1) higher than 15 years, (2) between 10 to 15 years, and (3) less than ten years.

A hospital may follow its policy to give the associated weighting to each category. Assume that the three types were assigned 50%, 30%, and 20%, respectively; the re-calculated (experience-weighted) DANP weight for each criterion will be changed accordingly.

Table 8. The number of the ten doctors in three categories.

Groups	Clinical Experience	Number of Doctors	Weighting
G1	>15 years	1	50%
G2	10 to 15 years	5	30%
G3	<10 years	4	20%

shows the ten doctors’ numbers in each category and

Table 9. Comparison of the original DANP weights and the re-calculated (experience-weighted) ones.

Original	C11	C12	C13	C14	C15	C21	C22	C23	C24	C25	C31	C32	C33	C34
DANP (%)	4.92	6.00	5.83	5.80	4.98	6.45	6.46	6.06	6.34	5.99	10.05	10.04	10.59	10.48
(Rank)	(14th)	(9th)	(11th)	(12th)	(13th)	(6th)	(5th)	(8th)	(7th)	(10th)	(3rd)	(4th)	(1st)	(2nd)
Re-calculated DANP (%)	4.96	5.89	6.00	5.23	4.58	6.41	6.41	6.03	6.42	6.20	10.28	10.16	10.87	10.57
(Rank)	(13th)	(11th)	(10th)	(12th)	(14th)	(6th)	(7th)	(9th)	(5th)	(8th)	(3rd)	(4th)	(1st)	(2nd)

compares the re-calculated DANP weights with the original ones (refer to Table 6).

The re-calculated (experience-weighted) DANP weights reveal the same ranking order for the top four criteria (i.e., $C_{33}\succ C_{34}\succ C_{31}\succ C_{32}$); however, the remaining criteria ranking orders are not entirely the same. The experience-weighted method can be applied to follow a hospital’s policy, which is more flexible in practice.

In this illustration, we use assumed cases to illustrate how to adopt the DANP model to evaluate patients. Once an elderly patient was regarded as risky by using the model, the DEMATEL analysis findings may further guide to improving their overall AKI risk level. For instance, in Table 7, case A’s score on Frailty (C₃₃) was 9.00 (very high). Based on the findings of DEMATEL analytics (refer to Figure 2 and Table 4), Activities of daily living (C₃₁) and Nutritional assessment (C₃₄) both would influence C₃₃. Thus, case A may consult medical professionals proactively in those two criteria to mitigate their AKI risk level. As suggested by [15], the increased AKI rate among elderly patients worsens with higher frailty. This MCDM approach proposes an overall AKI risk assessment model and provides precautionary guidance to reduce their potential AKI risk.

The MCDM approach incorporates doctors’ clinical knowledge to form a concise AKI evaluation model for elderly patients. Its short-term drawback might be lacking evidence to validate its predictive ability. However, once a hospital collects elderly patients’ scores on those criteria and organizes progress reports consistently, the model’s predictability may be further examined by statistics.

The study population is selected voluntarily; that is, the potential participants were chosen via referrals by the nephrologists, which would introduce selection bias potentially. Hawthorne effect is inevitable since the participants made a conscious decision to be in the selected hospital. Voluntary bias which could be viewed as coming from a particular sample could contain only those participants who are actually willing to participate in the survey and who participate and find the topic particularly interesting and are more likely to volunteer for that study, the same as those who are expected to be evaluated on a positive level [31].

Finally, we evaluated only one tertiary hospital, which might have characteristics that differ from those of the general population. The study’s external validity and generalization should be further considered. Future studies using a random sampling of hospitals over a broader range of regions would make the research more discursive.

5. Conclusions

To conclude, the present study proposes a complementary approach to enrich meta-analysis research findings [15] regarding AKI risk assessment for elderly patients. Its contributions are threefold.

(1)

Identified the influential relationships among the AKI risk assessment model for elderly patients.

(2)

Obtained the influential weights of the AKI evaluation criteria from ten experienced doctors.

(3)

Proposed a flexible method (experience-based weighting) to follow a hospital’s policy to form a decision support system (or termed as a knowledge-based system).

(4)

This study has shed light on integrating multiple research approaches (e.g., statistics and MCDM) to assist medical decision-making in practice.