# Assessing Dry Weight of Hemodialysis Patients via Sparse Laplacian Regularized RVFL Neural Network with L2,1-Norm

Authors: Xiaoyi Guo; Wei Zhou; Qun Lu; Aiyan Du; Yinghua Cai; Yijie Ding
Journal: BioMed Research International
 (2021)

---

## Abstract

Dry weight is the normal weight of hemodialysis patients after hemodialysis. If the amount of water in diabetes is too much (during hemodialysis), the patient will experience hypotension and shock symptoms. Therefore, the correct assessment of the patient’s dry weight is clinically important. These methods all rely on professional instruments and technicians, which are time-consuming and labor-intensive. To avoid this limitation, we hope to use machine learning methods on patients. This study collected demographic and anthropometric data of 476 hemodialysis patients, including age, gender, blood pressure (BP), body mass index (BMI), years of dialysis (YD), and heart rate (HR). We propose a Sparse Laplacian regularized Random Vector Functional Link (SLapRVFL) neural network model on the basis of predecessors. When we evaluate the prediction performance of the model, we fully compare SLapRVFL with the Body Composition Monitor (BCM) instrument and other models. The Root Mean Square Error (RMSE) of SLapRVFL is 1.3136, which is better than other methods. The SLapRVFL neural network model could be a viable alternative of dry weight assessment.


---

## Body

## 1. Introduction

Fluid overload in patients with chronic renal failure is closely related to poor cardiovascular outcomes [1, 2]. Maintenance of hemodialysis (HD) is the main method for patients with renal failure [3]. However, the accurate assessment of body water volume is still a concern [4]. At present, dry weight has been used as an important indicator to assess the homeostasis of fluids in hemodialysis patients. Medical staff can use the patient’s dry weight to estimate the amount of water needed for dialysis during hemodialysis. The conventional clinical-based dry weight assessment method is time-consuming and labor-intensive [1]. There are already some methods based on bioelectrical impedance analysis (BIA) [5] to determine dry weight, including body composition monitor (BCM) [6] and lung ultrasound (LUS). However, all the above methods require special instruments and professional technicians to complete. Medical staff can use some clinical data to build predictive models [7] to accurately assess dry weight. Currently, machine learning (ML) or deep learning has solved many common clinical problems in medicine, such as brain diseases [8–10], cancer analysis, and diabetes.Some scholars have used artificial neural networks (ANN) to predict the total water volume of hemodialysis patients and have obtained better results than conventional clinical calculation equations [11]. In addition, deep learning methods are also emerging in clinical diagnosis, including pixel-based convolutional neural networks to diagnose skin cancer [12]. In the biological field, microbiology analysis [13], CircRNAs [14], microRNAs, and cancer association prediction [15–17], lncRNA-miRNA association prediction, O-GlcNAcylation site prediction [18], DNA methylation site [19–21], protein remote homology [22], function prediction of proteins [23–29], electron transport proteins [30], breast cancer [31], cell-specific replication [32], osteoporosis diagnoses [33], and drug complex network analysis [34–38].In our previous research, a Multiple Kernel Support Vector Regression (MKSVR) [39] predictor was proposed to assess the dry weight and obtain good predictive performance. Inspired by the previous work and baseline Random Vector Functional Link (RVFL) network [40], we propose a new dry weight assessment model, called Sparse Laplacian regularized RVFL neural network with L2,1-norm (SLapRVFL), which considers the topological relationship between samples and more sparse connections between the input layer and the hidden layer.


## 2. Materials and Methods

### 2.1. Materials

This work collects demographic and anthropometric data and bioimpedance spectroscopy (BIS) from historical data (2018-9 to 2019-9) from Wuxi people’s hospital and the northern Jiangsu people’s hospital. This study has been approved by the ethics committees of the hospitals (Nos. KYLLKS201813 and 2018KY-001). The collected patient data meet the following requirements: age greater than 18 years; ESRD for more than three months and maintenance hemodialysis [41]; no heart failure, no metal implants, no pregnancy, no disability, no infection, and no edema and other diseases; and hemodialysis treatment 3 times a week, 4 hours each time. Finally, we obtain a data set of 476 hemodialysis patients. DW is the normal body weight after clinical diabetes. DW is obtained by a clinician under strict clinical supervision using a clinical scoring system (using trial and error method) [42, 43].We choose 7 features, including age, gender (binary feature), systolic blood pressure (SBP), diastolic blood pressure (DBP), body mass index (BMI), heart rate (HR), and years of dialysis (YD) to build our predictive model. Table1 shows the information of the data set. BMI is measured before hemodialysis treatment.Table 1
The information of data set.
FeatureValuer∗Age (years)54.17±14.22-0.2341Gender (males/females)312/164-0.4489BMI22.96±2.950.9558Systolic blood pressure (mmHg)150.64±29.36-0.1739Diastolic blood pressure (mmHg)88.32±19.56-0.1249Heart rate (times/min)73.41±8.920.1862Years of dialysis (years)5.97±3.22-0.1069∗Denotes that each feature correlated with dry weight using Pearson correlation coefficient (r).


### 2.2. Methods

The baseline RVFL was proposed for regression or classification. The schematic diagram of RVFL is shown in Figure1. The basic information of the patient is put into the RVFL neural network model for processing, and the predicted dry weight is the output.Figure 1
Schematic of our proposed method.Suppose, there areN training samples with xi,yi,i=1,2,⋯,N. The output value is yi∈R1×c and the input data is xi∈R1×d. d denotes the dimension of xi. As per Figure 1, RVFL randomly initializes all weights and deviations between the hidden layer and the input layer. These parameters are fixed during the training process and do not need to be tuned. There are connections between the output layer, input layer, and hidden layer. This part of the weight needs to be obtained by training RVFL. The output layer of RVFL is connected to both the input layer and the hidden layer, so as to ensure the nonlinear and linear relationships between the input and the output. The RVFL network with P hidden nodes are formulated as
(1)Hβ=Y,where β denotes the output weight matrix; H is the concatenated matrix, which combines the output of the hidden layer and the input layer; and Y denotes the label matrix. H and β can be represented as
(2)H=H1H2,(3)H1=x11⋯x1d⋮⋱⋮xN1⋯xNdN×d,(4)H2=Ga1x1+b1⋯GaPx1+bP⋮⋱⋮Ga1xN+b1⋯GaPxN+bPN×P,(5)β=β1Tβ2T⋮βd+PTd+P×C.In Equation (4), aj and bj are the weights and bias of the hidden and input layers. C and P are numbers of output and hidden layer nodes. In general, the activation function is a Gaussian function: gx=e−x2. The activation function has a nonlinear approximation effect. To consider the potential linear relationship between the input data and the output value, RVFL adds a direct connection weight between the input layer and the output layer. Therefore, RVFL is a model that contains both linear and nonlinear approximations to improve prediction performance. For optimal β, the RVFL can be formulated as a regularized least-squares:
(6)β∗=argmin12Hβ−Y22+λ2β22,where λ is the parameter of regularization term. The solution of Equation (6) can be found by setting its gradient to 0:
(7)β∗=HTH+λI−1HTY,where I denotes the identity matrix. However, the RVFL network did not consider the topological relationship between samples. For the output node, it must be connected to both the input and the hidden layer.In order to further improve the robustness of RVFL, we propose Sparse Laplacian regularized RVFL neural network with L2,1-norm (SLapRVFL). The objective function is
(8)β∗=argmin12Hβ−Y22+λ12TrHβTLHβ+λ22β2,12,where L∈RN×N denotes the Laplacian matrix. λ1 and λ2 are the coefficients of Laplacian regularization the and L21-norm term, respectively. Laplacian regularization is used to indicate the potential manifold between samples. It can better describe the topological association between samples to improve the generalization ability of the model. Since the third term of β2,12 is not diversified, we convert Equation (8) to
(9)β∗=argmin12Hβ−Y22+λ12TrHβTLHβ+λ22TrβTGβ,where G∈Rd+P×d+P denotes a diagonal matrix whose ith-diagonal element
(10)Gii=12βi2,i=1,2,⋯,d+P.We take the derivative of the formula Equation (10) as
(11a)HTHβ−Y+λ1HTLHβ+λ2Gβ=0,(11b)HTHβ+λ1HTLHβ+λ2Gβ=HTY,(11c)HTH+λ1HTLH+λ2Gβ=HTY,(11d)β=HTH+λ1HTLH+λ2G−1HTY.We use the baseline RVFL solution with Equation (7) as the initial β0. In addition, the Laplacian matrix can be calculate as
(12a)L=D−1/2ΔD−1/2,(12b)Δ=D−S,where D is diagonal matrix, Dii=∑j=1NSij. Similarity matrix S is built by Radial Basis Function (RBF):
(13)Sij=exp−γxi−xj2.The process of SLapRVFL is list in Algorithm1.Algorithm 1. Algorithm of SLapRVFL
Require: Training set xi,yi,i=1,2,⋯,N, test set xjte,j=1,2,⋯,M, the numbers of hidden layer nodes (P), the maximum number of iterations tmax, coefficients of λ1 and λ2;Ensure: The predictive values of yjte,j=1,2,⋯,M(1) Randomly initializing all weights and deviations between the hidden layer and the input layer. Calculating the hidden layer output matrixH (training set)and Laplacian matrix L by Equations (2), (12), and (13);(2) Sett=0, estimate the initial β0 using Equation (7);Repeat(3) Update the diagonal matrixG withGt+1=1/2β1t2⋱1/2βd+Pt2d+P×d+P,(4) Updateβ via Equation (11d);Untilt>tmax;(5) Calculate the hidden layer output matrixHte (test set);(6) Estimateyjte,j=1,2,⋯,M by Yte=Hteβ.


## 2.1. Materials

This work collects demographic and anthropometric data and bioimpedance spectroscopy (BIS) from historical data (2018-9 to 2019-9) from Wuxi people’s hospital and the northern Jiangsu people’s hospital. This study has been approved by the ethics committees of the hospitals (Nos. KYLLKS201813 and 2018KY-001). The collected patient data meet the following requirements: age greater than 18 years; ESRD for more than three months and maintenance hemodialysis [41]; no heart failure, no metal implants, no pregnancy, no disability, no infection, and no edema and other diseases; and hemodialysis treatment 3 times a week, 4 hours each time. Finally, we obtain a data set of 476 hemodialysis patients. DW is the normal body weight after clinical diabetes. DW is obtained by a clinician under strict clinical supervision using a clinical scoring system (using trial and error method) [42, 43].We choose 7 features, including age, gender (binary feature), systolic blood pressure (SBP), diastolic blood pressure (DBP), body mass index (BMI), heart rate (HR), and years of dialysis (YD) to build our predictive model. Table1 shows the information of the data set. BMI is measured before hemodialysis treatment.Table 1
The information of data set.
FeatureValuer∗Age (years)54.17±14.22-0.2341Gender (males/females)312/164-0.4489BMI22.96±2.950.9558Systolic blood pressure (mmHg)150.64±29.36-0.1739Diastolic blood pressure (mmHg)88.32±19.56-0.1249Heart rate (times/min)73.41±8.920.1862Years of dialysis (years)5.97±3.22-0.1069∗Denotes that each feature correlated with dry weight using Pearson correlation coefficient (r).


## 2.2. Methods

The baseline RVFL was proposed for regression or classification. The schematic diagram of RVFL is shown in Figure1. The basic information of the patient is put into the RVFL neural network model for processing, and the predicted dry weight is the output.Figure 1
Schematic of our proposed method.Suppose, there areN training samples with xi,yi,i=1,2,⋯,N. The output value is yi∈R1×c and the input data is xi∈R1×d. d denotes the dimension of xi. As per Figure 1, RVFL randomly initializes all weights and deviations between the hidden layer and the input layer. These parameters are fixed during the training process and do not need to be tuned. There are connections between the output layer, input layer, and hidden layer. This part of the weight needs to be obtained by training RVFL. The output layer of RVFL is connected to both the input layer and the hidden layer, so as to ensure the nonlinear and linear relationships between the input and the output. The RVFL network with P hidden nodes are formulated as
(1)Hβ=Y,where β denotes the output weight matrix; H is the concatenated matrix, which combines the output of the hidden layer and the input layer; and Y denotes the label matrix. H and β can be represented as
(2)H=H1H2,(3)H1=x11⋯x1d⋮⋱⋮xN1⋯xNdN×d,(4)H2=Ga1x1+b1⋯GaPx1+bP⋮⋱⋮Ga1xN+b1⋯GaPxN+bPN×P,(5)β=β1Tβ2T⋮βd+PTd+P×C.In Equation (4), aj and bj are the weights and bias of the hidden and input layers. C and P are numbers of output and hidden layer nodes. In general, the activation function is a Gaussian function: gx=e−x2. The activation function has a nonlinear approximation effect. To consider the potential linear relationship between the input data and the output value, RVFL adds a direct connection weight between the input layer and the output layer. Therefore, RVFL is a model that contains both linear and nonlinear approximations to improve prediction performance. For optimal β, the RVFL can be formulated as a regularized least-squares:
(6)β∗=argmin12Hβ−Y22+λ2β22,where λ is the parameter of regularization term. The solution of Equation (6) can be found by setting its gradient to 0:
(7)β∗=HTH+λI−1HTY,where I denotes the identity matrix. However, the RVFL network did not consider the topological relationship between samples. For the output node, it must be connected to both the input and the hidden layer.In order to further improve the robustness of RVFL, we propose Sparse Laplacian regularized RVFL neural network with L2,1-norm (SLapRVFL). The objective function is
(8)β∗=argmin12Hβ−Y22+λ12TrHβTLHβ+λ22β2,12,where L∈RN×N denotes the Laplacian matrix. λ1 and λ2 are the coefficients of Laplacian regularization the and L21-norm term, respectively. Laplacian regularization is used to indicate the potential manifold between samples. It can better describe the topological association between samples to improve the generalization ability of the model. Since the third term of β2,12 is not diversified, we convert Equation (8) to
(9)β∗=argmin12Hβ−Y22+λ12TrHβTLHβ+λ22TrβTGβ,where G∈Rd+P×d+P denotes a diagonal matrix whose ith-diagonal element
(10)Gii=12βi2,i=1,2,⋯,d+P.We take the derivative of the formula Equation (10) as
(11a)HTHβ−Y+λ1HTLHβ+λ2Gβ=0,(11b)HTHβ+λ1HTLHβ+λ2Gβ=HTY,(11c)HTH+λ1HTLH+λ2Gβ=HTY,(11d)β=HTH+λ1HTLH+λ2G−1HTY.We use the baseline RVFL solution with Equation (7) as the initial β0. In addition, the Laplacian matrix can be calculate as
(12a)L=D−1/2ΔD−1/2,(12b)Δ=D−S,where D is diagonal matrix, Dii=∑j=1NSij. Similarity matrix S is built by Radial Basis Function (RBF):
(13)Sij=exp−γxi−xj2.The process of SLapRVFL is list in Algorithm1.Algorithm 1. Algorithm of SLapRVFL
Require: Training set xi,yi,i=1,2,⋯,N, test set xjte,j=1,2,⋯,M, the numbers of hidden layer nodes (P), the maximum number of iterations tmax, coefficients of λ1 and λ2;Ensure: The predictive values of yjte,j=1,2,⋯,M(1) Randomly initializing all weights and deviations between the hidden layer and the input layer. Calculating the hidden layer output matrixH (training set)and Laplacian matrix L by Equations (2), (12), and (13);(2) Sett=0, estimate the initial β0 using Equation (7);Repeat(3) Update the diagonal matrixG withGt+1=1/2β1t2⋱1/2βd+Pt2d+P×d+P,(4) Updateβ via Equation (11d);Untilt>tmax;(5) Calculate the hidden layer output matrixHte (test set);(6) Estimateyjte,j=1,2,⋯,M by Yte=Hteβ.


## 3. Results

We test our model on the benchmark data set and obtain the optimal parameters of the predictor through cross-validation. The SLapRVFL network is compared to other machine learning-based models. In addition, the body composition monitor (BCM) device (Fresenius Medical Care, Baden Humboldt, Germany) is also compared with the SLapRVFL network.


### 3.1. Evaluation Measurements

The 10-fold cross-validation (10-CV) is employed to evaluate the robustness of methods. Root Mean Square Error (RMSE),R square, correlation coefficient (R), Bland–Altman analysis, and Empirical Cumulative Distribution Plot (ECDP) [44] are all used in our study. To evaluate the agreement of two different methods, the Bland–Altman analysis usually can obtain whether the two methods can be substituted for each other (equivalence). Evaluating the agreement of the two methods can answer the question, “Can these two methods replace each other?”


### 3.2. Selection of Optimal Parameters

To get the optimal parameters of the predictive method, we obtain them through a grid search method. The parameters that need to be determined include the numbers of hidden layer nodesP, maximum iterations, and coefficients of λ1 and λ2. For the numbers of hidden layer nodes P, we fix the iterations, λ1 and λ2. Setting the maximum number as 50, λ1=1 and λ2=1. The value of P is from 10 to 140 with step of 10. The results are shown in Figure 2. From 10 to 100, the more neurons in the hidden layer, the lower the RMSE. Since then, RMSE has gradually increased. So, we get the lower RMSE under P=100.Figure 2
The RMSE under different numbers of hidden layer nodes (SLapRVFL network).Next,P=100, λ1=1, and λ2=1. We gradually increase the number of iterations from 1 to 100 (shown in Figure 3). After the number of iterations reaches 10, the RMSE value drops to a minimum and slightly oscillates within a certain value. In our study, maximum number of iterations is 10.Figure 3
The RMSE of iterations on the training set.Then, we use the better number of hidden layer nodes and iterations to search for the bestλ1 and λ2. The search range of parameters is from 2−5 to 20 (with step of 20.5). Figure 4 shows the results of different parameters. When λ1 and λ2 are 2−3 and 2−2.5, RMSE is the lowest.Figure 4
The RMSE under differentλ1 and λ2.


### 3.3. Comparison to Other Predictive Models and BCM

To evaluate our model, SLapRVFL is compared with our previous work of Multiple Kernel Support Vector Regression (MKSVR) [39], Multikernel Ridge Regression (MKRR), Linear Regression (LR), Artificial Neural Network based on Back Propagation algorithm (ANN with BP), and BCM measuring instrument. Clinical dry weight is our reference standard (also the regression target value of the prediction model). The comparisons are listed in Table 2, which shows that SLapRVFL achieves best performance of RMSE (1.3136). Although the ECDP median value (peak) of MKSVR (0.0082) is more close to zero, Figure 5 shows that SLapRVFL has the least bias and much less tails than MKSVR (smaller width). The RMSE of BCM is 1.9694, which is larger than SLapRVFL.Table 2
Comparison on existing methods via 10-fold cross-validation.
MethodRR squaredRMSEEmpirical cumulative distribution plotHighest valueLowest valueMedian valueBCM∗0.94730.91371.96943.2235-6.2776-0.9863LR∗0.94030.93081.43354.2524-4.40140.1418ANN (BP)∗0.93980.92951.47947.3661-4.74470.1324MKRR∗0.93990.92891.50154.9227-4.26040.1104MKSVR∗0.94120.93211.38174.3962-4.12730.0082RVFL0.93890.93001.38286.7004-4.35570.0704SLapRVFL (our method)0.96320.95011.31363.1940-3.50660.1014∗The results are from previous work on MKSVR [39].Figure 5
Folded empirical cumulative distribution plot between different methods.


### 3.4. Bland–Altman Analysis

Bland–Altman plot is a useful tool to evaluate the agreement between predictive methods and clinical DW. In Table3 and Figure 6, SLapRVFL, MKSVR, LR, ANN (BP), MKRR, and BCM are analyzed via Bland-Altman difference plot. SLapRVFL achieves the smallest range of 95% confidence interval (-0.1133 to 0.2866) and standard deviation (2.2202). In addition, the number (ratio) of outside agreement interval for predictive models is all less than 24 (5%) predictive samples. These results of models are clinically acceptable. SLapRVFL achieves least number (20) of the outside agreement interval in Table 3. As shown in Figure 6, two red horizontal dotted lines (upper and lower) denote the upper and lower limits of the 95% agreement limit, respectively. The middle blue solid line is the average value of the difference (between measurement methods and clinical DW). While one measurement method and clinical method can be considered as a better agreement, they can be substituted for each other (equivalence). If 95% of the points of the data set are in the agreement range, the measurement method (predictive model) is clinically acceptable. The results of the evaluation show that SLapRVFL can help clinicians assess DW with low cost.Table 3
Bland–Altman plot analysis for different models.
ModelDifferences with DW (%)Limits of agreement (%)MeanSD95% confidence intervalLower limitUpper limitNumber (ratio) of outside agreement intervalBCM∗-1.82322.7466-2.0706 to -1.5759-7.20663.560130/476 (6.30%)LR∗0.00022.4269-0.2184 to 0.2187-4.75664.756921/476 (4.41%)ANN (BP)∗0.11522.5139-0.1112 to 0.3416-4.81195.042422/476 (4.62%)MKRR∗-0.08012.5007-0.3053 to 0.1451-4.98144.821223/476 (4.83%)MKSVR∗-0.26382.3372-0.4743 to -0.05329-4.84464.317122/476 (4.62%)SLapRVFL (our method)0.08672.2202-0.1133 to 0.2866-4.26504.438320/476 (4.20%)∗The results are from previous work on MKSVR [39].Figure 6
Bland–Altman plot analysis.
(a)
ANN(b)
LR(c)
MKRR(d)
MKSVR(e)
BCM(f)
SLapRVFL


## 3.1. Evaluation Measurements

The 10-fold cross-validation (10-CV) is employed to evaluate the robustness of methods. Root Mean Square Error (RMSE),R square, correlation coefficient (R), Bland–Altman analysis, and Empirical Cumulative Distribution Plot (ECDP) [44] are all used in our study. To evaluate the agreement of two different methods, the Bland–Altman analysis usually can obtain whether the two methods can be substituted for each other (equivalence). Evaluating the agreement of the two methods can answer the question, “Can these two methods replace each other?”


## 3.2. Selection of Optimal Parameters

To get the optimal parameters of the predictive method, we obtain them through a grid search method. The parameters that need to be determined include the numbers of hidden layer nodesP, maximum iterations, and coefficients of λ1 and λ2. For the numbers of hidden layer nodes P, we fix the iterations, λ1 and λ2. Setting the maximum number as 50, λ1=1 and λ2=1. The value of P is from 10 to 140 with step of 10. The results are shown in Figure 2. From 10 to 100, the more neurons in the hidden layer, the lower the RMSE. Since then, RMSE has gradually increased. So, we get the lower RMSE under P=100.Figure 2
The RMSE under different numbers of hidden layer nodes (SLapRVFL network).Next,P=100, λ1=1, and λ2=1. We gradually increase the number of iterations from 1 to 100 (shown in Figure 3). After the number of iterations reaches 10, the RMSE value drops to a minimum and slightly oscillates within a certain value. In our study, maximum number of iterations is 10.Figure 3
The RMSE of iterations on the training set.Then, we use the better number of hidden layer nodes and iterations to search for the bestλ1 and λ2. The search range of parameters is from 2−5 to 20 (with step of 20.5). Figure 4 shows the results of different parameters. When λ1 and λ2 are 2−3 and 2−2.5, RMSE is the lowest.Figure 4
The RMSE under differentλ1 and λ2.


## 3.3. Comparison to Other Predictive Models and BCM

To evaluate our model, SLapRVFL is compared with our previous work of Multiple Kernel Support Vector Regression (MKSVR) [39], Multikernel Ridge Regression (MKRR), Linear Regression (LR), Artificial Neural Network based on Back Propagation algorithm (ANN with BP), and BCM measuring instrument. Clinical dry weight is our reference standard (also the regression target value of the prediction model). The comparisons are listed in Table 2, which shows that SLapRVFL achieves best performance of RMSE (1.3136). Although the ECDP median value (peak) of MKSVR (0.0082) is more close to zero, Figure 5 shows that SLapRVFL has the least bias and much less tails than MKSVR (smaller width). The RMSE of BCM is 1.9694, which is larger than SLapRVFL.Table 2
Comparison on existing methods via 10-fold cross-validation.
MethodRR squaredRMSEEmpirical cumulative distribution plotHighest valueLowest valueMedian valueBCM∗0.94730.91371.96943.2235-6.2776-0.9863LR∗0.94030.93081.43354.2524-4.40140.1418ANN (BP)∗0.93980.92951.47947.3661-4.74470.1324MKRR∗0.93990.92891.50154.9227-4.26040.1104MKSVR∗0.94120.93211.38174.3962-4.12730.0082RVFL0.93890.93001.38286.7004-4.35570.0704SLapRVFL (our method)0.96320.95011.31363.1940-3.50660.1014∗The results are from previous work on MKSVR [39].Figure 5
Folded empirical cumulative distribution plot between different methods.


## 3.4. Bland–Altman Analysis

Bland–Altman plot is a useful tool to evaluate the agreement between predictive methods and clinical DW. In Table3 and Figure 6, SLapRVFL, MKSVR, LR, ANN (BP), MKRR, and BCM are analyzed via Bland-Altman difference plot. SLapRVFL achieves the smallest range of 95% confidence interval (-0.1133 to 0.2866) and standard deviation (2.2202). In addition, the number (ratio) of outside agreement interval for predictive models is all less than 24 (5%) predictive samples. These results of models are clinically acceptable. SLapRVFL achieves least number (20) of the outside agreement interval in Table 3. As shown in Figure 6, two red horizontal dotted lines (upper and lower) denote the upper and lower limits of the 95% agreement limit, respectively. The middle blue solid line is the average value of the difference (between measurement methods and clinical DW). While one measurement method and clinical method can be considered as a better agreement, they can be substituted for each other (equivalence). If 95% of the points of the data set are in the agreement range, the measurement method (predictive model) is clinically acceptable. The results of the evaluation show that SLapRVFL can help clinicians assess DW with low cost.Table 3
Bland–Altman plot analysis for different models.
ModelDifferences with DW (%)Limits of agreement (%)MeanSD95% confidence intervalLower limitUpper limitNumber (ratio) of outside agreement intervalBCM∗-1.82322.7466-2.0706 to -1.5759-7.20663.560130/476 (6.30%)LR∗0.00022.4269-0.2184 to 0.2187-4.75664.756921/476 (4.41%)ANN (BP)∗0.11522.5139-0.1112 to 0.3416-4.81195.042422/476 (4.62%)MKRR∗-0.08012.5007-0.3053 to 0.1451-4.98144.821223/476 (4.83%)MKSVR∗-0.26382.3372-0.4743 to -0.05329-4.84464.317122/476 (4.62%)SLapRVFL (our method)0.08672.2202-0.1133 to 0.2866-4.26504.438320/476 (4.20%)∗The results are from previous work on MKSVR [39].Figure 6
Bland–Altman plot analysis.
(a)
ANN(b)
LR(c)
MKRR(d)
MKSVR(e)
BCM(f)
SLapRVFL


## 4. Discussion

Due to the limitations of clinical and BCM measurement (more time and cost), this study uses a machine learning method to assess the dry weight of hemodialysis patients. Based on the basic RVFL, we propose a sparse Laplace regularized RVFL network (SLapRVFL) model. SLapRVFL is compared not only with other machine learning methods (such as LR, MKRR, ANN with BP, and MKSVR) but also with BCM equipment (commonly used in hospitals). The RMSE and Bland–Altman analysis of the model are better than the BCM instrument. It is proven that the predictive model driven by data can provide reference for clinical dry weight assessment.BCM requires the patient’s information on weight (before hemodialysis) and height. It is a portable, inexpensive, and noninvasive technology that has been used to measure DW [45, 46]. For the Bland–Altman analysis, SLapRVFL achieves the least number (20) of outside agreement interval. However, BCM has 30/476 (6.30%) points (ratio) of the outside agreement interval. Obviously, our method has better agreement with the clinical method.


## 5. Conclusions

To further improve the robustness of RVFL, we introduce sparse Laplacian regular term with L2,1-norm. In the training process, the graph topology information and the sparse weight matrix (output) are employed to improve the robustness of the RVFL. In fact, our work provides a new idea for assessing patients’ dry weight. Not only that, in the fields of biology [47–57], pharmacy [58], and medicine [12, 59, 60], machine learning methods have helped solve many analysis tasks. In future research, we will consider collecting more samples, introducing more patient personal information, and building a predictor based on a deep learning model to more accurately assess the dry weight of hemodialysis patients.



---

*Source: 6627650-2021-02-05.xml*