# Improved Ensemble Extreme Learning Machine Regression Algorithm

Meiyi Li, Weibiao Cai, Xingwang Liu

► To cite this version:

Meiyi Li, Weibiao Cai, Xingwang Liu. Improved Ensemble Extreme Learning Machine Regression Algorithm. 10th International Conference on Intelligent Information Processing (IIP), Oct 2018, Nanning, China. pp.12-19, 10.1007/978-3-030-00828-4_2. hal-02197774

HAL Id: hal-02197774

https://inria.hal.science/hal-02197774

Submitted on 30 Jul 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Improved Ensemble Extreme Learning Machine
Regression Algorithm

Meiyi Li¹, Weibiao Cai², Xingwang Liu³

¹ Xiangtan University College of Information Engineering,
Xiangtan University, Xiangtan, China,
1356695709@qq.com

² Xiangtan University College of Information Engineering,
Xiangtan University, Xiangtan, China,
2481596161@qq.com

³ Information Center of Jianglu mechanical & electrical group Co., Ltd.

**ABSTRACT.** Compared with other traditional neural network algorithms, the Extreme Learning Machine (ELM) has the advantages of simple structure, fast learning speed and good generalization performance. However, there are still some shortages that restrict the further development of ELM. For example, the randomly generated input weights, biases and the ill-conditioned appearance of the hidden layer design matrix all affect the generalization performance and robustness of the ELM algorithm model. In order to overcome the adverse affects of both, an improved ensemble extreme learning machine regression algorithm(ECV-ELM) is proposed in this paper. The method first generates multiple sub CV-ELM model through AdaBoost.RT method, and the selects the best set of submodels to integrate. The ECV-ELM algorithm makes use of ensemble learning method to complement each other among sub models, so that the generalization performance and robustness of the algorithm are better than that of the sub model. The results of regression experiments on multiple data sets show that the ECV_ELM algorithm can effectively reduce the influence of the ill conditioned matrix, the random input weight and bias, and has good generalization performance and robustness.

**KEYWORDS:** extreme learning machine, ensemble learning, robustness, CV-ELM algorithm

# 1 Introduction

Artificial neural network (ANN) has been widely used in various fields due to its good learning ability and high speed optimization capability [1-2]. However, the artificial neural network is used to calculate the parameters of the model by the gradient descent method. The gradient descent method will increase the time complexity of the algorithm and cause the algorithm to fall into the local optimal solution easily. Due to these shortcomings, it may take more time to train the network and not necessarily guarantee the best solution. Therefore, searching for high
efficiency and high real-time neural network has become the main research direction of many scholars.

The Extreme Learning Machine (ELM) [3] is a new single-hidden-layer feedforward neural networks (SLFNs) proposed by G.B. Huang et al. in 2006. The method can randomly generate input weights and hidden layer node thresholds, and calculate the output layer weights without the need for iterative operations, so its computation speed is much better than the traditional neural network algorithm. When constructing the model, the traditional neural network algorithm is different from the traditional neural network algorithm. The limit learning function gets the minimum norm of the weight of the output layer while reaching the minimum training error. According to the neural network theory: For feed-forward neural networks, the smaller the training error and norm of the weight of the LM model, the stronger its network generalization ability [4]. Therefore, we can theoretically demonstrate the generalization ability of ELM algorithm. In recent years, ELM has been applied to various fields of various disciplines, such as face recognition [5], classified, regression, image processing [6], ground reconstruction [7], and so on, because of the ability to effectively improve the defects of traditional neural networks.

Although ELM has fast learning speed and good generalization performance, when some column vectors in the hidden layer design array are approximated to linear correlation, that is, the hidden layer design array has multiple collinearity or ill posed. It can result in poor generalization performance and stability by using ordinary least square method to estimate the solution of the ill conditioned matrix. To this end, many scholars have improved the ELM model, for example, M.Y.Li et al. Proposed an improved limit learning machine regression algorithm (CV-ELM) based on conditional index and variance decomposition ratio (CV-ELM) [8], Ceng Lin and others proposed the limit learning machine (PC-ELM) [9] based on the principal component estimation.

The CV-ELM algorithm improves the generalization performance of the extreme learning machine to a certain extent and can ensure good algorithm robustness. However, this algorithm still has defects in some cases, so that the model can not achieve the minimum error. Our analysis considers that there are two main reasons: First, high-dimensional data may obscure the noise components in the data, making the proposed method unable to completely isolate the relevant variables, that is, the algorithm can not reduce the noise completely. Secondly, some initial parameters of CV-ELM, such as input weight and hidden layer bias, are generated randomly, which makes the generalization ability and robustness of the model affected by them. Therefore, this paper introduces the ensemble learning method [10] in the CV-ELM algorithm. This method can train a number of similar learners at the same time, and then extract the best set of learners from these neural networks to integrate [11]. Through regression experiments of multiple data sets, it is proved that this method can achieve good generalization performance and stability.

## 2 Review of ELM and CV-ELM

In this section, we mainly review ELM[3] and CV-ELM[10].
## 2.1 Extreme learning machine

For N arbitrary distinct samples $(x_i, y_i)$, where $x_i = [x_{i1}, x_{i2}, ..., x_{in}]^T \in \mathbb{R}^n$ is an n-dimensional feature of the ith sample, and $y_i = [y_{i1}, y_{i2}, ..., y_{im}]^T \in \mathbb{R}^m$, Then, the output of a feedforward neural network with L hidden nodes and excitation function G(x) can be expressed as

$$f_L(x) = \sum_{i=1}^{L} \beta_i G(a_i \cdot x_i + b_i), a_i \in \mathbb{R}^n, \beta_i \in \mathbb{R}^m, \quad [1]$$

where $a_i = [a_{i1}, a_{i2}, ..., a_{in}]^T$ is the weight connecting the i-th hidden node and the input nodes, and $b_i$ is threshold of the i-th hidden nodes, $\beta_i = [\beta_{i1}, \beta_{i2}, ..., \beta_{im}]^T$ is the weight connecting the i-th hidden node and the output nodes, $a_i \cdot x_i$ denote the inner product of $a_i$ and $x_i$. The excitation function G(x) can choose "Sigmoid", "Sine" or "RBF" and so on.

If this feed forward neural network with L hidden layer nodes and M output layer nodes can approximate this N samples with zero error, then the above N equations can be written compactly as

$$f_L(x) = \sum_{i=1}^{L} \beta_i G(a_i \cdot x_i + b_i) = y_i, i = 1, 2, \dots, L, \quad [2]$$

(2) can be simplified as

$$H\beta = Y \quad [3]$$

where

$$H = \begin{bmatrix} G(a_1, b_1, x_1) & \cdots & G(a_L, b_L, x_1) \\ \vdots & \ddots & \vdots \\ G(a_1, b_1, x_N) & \cdots & G(a_L, b_L, x_N) \end{bmatrix}_{N \times L} \quad [4]$$

$$\beta = \begin{bmatrix} \beta_1^T \\ \vdots \\ \beta_L^T \end{bmatrix}_{L \times M} \qquad Y = \begin{bmatrix} y_1^T \\ \vdots \\ y_N^T \end{bmatrix}_{N \times M} \quad [5]$$

H is called the hidden layer output matrix of the network, and ELM training can be transformed into a problem of solving the least squares solution of output weights. The output weight matrix $\hat{\beta}$ can be obtained from (6)

$$\hat{\beta} = (H^T H)^{-1} H^T Y = H^+ Y \quad [6]$$

Where $H^+$ represents the Moore-penrose generalized inverse of the hidden layer output matrix H.
## 2.2 CV-ELM

The hidden layer matrix H is calculated according to the ELM model, and then the H matrix is separated according to the condition number and variance decomposition ratio to obtain

$$H = (H1, H2) \qquad [7]$$

where H1 is the non-interference data column in the hidden layer matrix H and H2 is the interference data columns.

According to the LS principle, the output weight matrix $\hat{\beta}$ can be obtained from (8)

$$ \begin{aligned} \hat{\beta} &= (H^T H)^{-1} H^T Y \\ &= ([H1, H2]^T [H1, H2])^{-1} [H1, H2] Y \\ &= \begin{bmatrix} H1^T H1 & H1^T H2 \\ H2^T H1 & H2^T H2 \end{bmatrix}^{-1} [H1, H2] Y \end{aligned} \qquad [8] $$

In order to enhance the generalization performance and stability of the model without destroying the authenticity of the data of the non-interference data, and add a small constant to the diagonal elements of the data matrix of the interfering data. The output weight matrix $\hat{\beta}$ can be obtained from (9)

$$ \hat{\beta} = \begin{bmatrix} H1^T H1 & H1^T H2 \\ H2^T H1 & H2^T H2 + kI \end{bmatrix}^{-1} [H1, H2] Y \qquad [9] $$

where k is a small constant and I is the unit matrix.

CV-ELM algorithm is described as follows:

Known training samples $(x_i, y_i)$, $i = 1, \dots, N$, the number of hidden nodes is L, and the excitation function is $G(x)$.

(1) Random setting of input weights $a_i$ and bias $b_i$, $i = 1, \dots, L$

(2) Computing hidden layer output matrix H

(3) Through the condition number and variance decomposition machine decomposition matrix H, and get $(H1, H2)$

(4) Determining the ridge parameter k

(5) The output matrix $\hat{\beta}$ is calculated by Equation (10)

# 3 Improved Ensemble Extreme learning Machine

The CV-ELM algorithm overcomes the situation that the generalized performance and robustness of the algorithm deteriorate when the hidden layer design array is ill-conditioned. However, the random generation of input weights and incomplete cancellation of noise under high-dimensional data can cause generalization performance and robustness to be poorly processed. Therefore, this paper proposes a CV-ELM regression algorithm based on ensemble learning (ECV-ELM), which makes use of the complementarity between multiple learners, thus making the
ensemble better performance. ECV-ELM overcomes the shortcomings of poor model stability due to the random generation of input weight, bias and incomplete cancellation of CV-ELM noise. It combines the ensemble learning method with the CV-ELM regression algorithm and uses some common methods to select the appropriate sub CV-ELM model, which can further improve the performance of the entire CV-ELM.

It is assumed that the training set and test set are $G = \{(x_i, y_i) | i = 1, 2, \dots, 1\}$, $G' = \{(x_i, y_i) | i = 1, 2, \dots, 1\}$, and $x_i$ is the model input, and the $y_i$ is the output of the model. First, according to the training set G come into the different training sub set $G = \{G_1, \dots, G_T\}$, several different sub CV-ELM models are generated by different training subsets. Then, a part of the excellent sub CV-ELM model is selected according to the training results. Finally, the results are made by means of the average method.

In summary, the proposed ECV-ELM integrated regression algorithm can be summarized as follows:

* Input: Training sample set T

* Output: Integrated CV-ELM regression model

(1) Using training set G to randomly generate T intersecting data sub sets $G = \{G_1, G_2, \dots, G_T\}$, set the activation function of all sub models to $g(x)$, and the number of hidden layer neurons is L;

(2) Initialization t=1;

(3) Determine whether it reaches the number of iterations, that is, $t<=T$; if yes, execute step (4); otherwise, execute step (7);

(4) Using the random function to generate the input weight a and the hidden layer offset b;

(5) The t-th sub CV-ELM model is trained using the randomly generated a, b, and t-th data subsets;

(6) Perform step (3);

(7) Calculate the MSE values of all the submodels. According to the size of the MSE values, select the k best sub CV-ELM models;

$$MSE = \frac{1}{n(Y - H\hat{\beta})^T(Y - H\hat{\beta})} \quad [10]$$

(8) using the simple average method to integrate K sub CV-ELM and get the final model, that is, the ECV-ELM model.

# 4 Experiment and Analysis

This section analyzes the CV-ELM regression method (ECV-ELM) based on ensemble learning proposed in the previous section. In order to better carry out experimental analysis, the time and prediction results of standard ELM algorithm,
CV-ELM algorithm and ECV-ELM algorithm are compared. The experiment uses 5 regression analysis data sets from UCI database and LIACC, which are Balloon data set, California House data set, Cloud dataset, Strike data set, Bodyfat data set. There are huge differences between the input attributes and the number of samples in these five data sets, which can better analyze the performance of the algorithm, as shown in Table 1. The number of hidden layer nodes required for each dataset in the ECV-ELM neutron model is shown in table 2. All experiments in this chapter are run on Windows7 64 bit operating system and Matlab 2016 environment in 3.30GHz i5-4590 CPU, 4G RAM.

<table><caption>**Table 1.** Regression Analysis Dataset</caption><thead><tr><th rowspan="2">Data set</th><th rowspan="2">Attributes</th><th colspan="2">samples</th></tr><tr><th>Training</th><th>Testing</th></tr></thead><tbody><tr><td>Cloud</td><td>9</td><td>599</td><td>1499</td></tr><tr><td>Strike</td><td>6</td><td>416</td><td>209</td></tr><tr><td>Balloon</td><td>1</td><td>1334</td><td>667</td></tr><tr><td>Bodyfat</td><td>14</td><td>168</td><td>84</td></tr><tr><td>House</td><td>8</td><td>8000</td><td>12640</td></tr></tbody></table>

<table><caption>**Table 2.** Number of required hidden layer nodes in each dataset</caption><thead><tr><th>Data set</th><th>Cloud</th><th>Strike</th><th>Balloon</th><th>Bodyfat</th><th>House</th></tr></thead><tbody><tr><th>Hidden nodes</th><td>80</td><td>110</td><td>100</td><td>30</td><td>100</td></tr></tbody></table>

Table 3 shows the comparison of test time and training time of ELM, ECV-ELM and CV-ELM on multiple data sets. Table 4 shows the comparison of RMSE of ELM, ECV-ELM and CV-ELM on multiple data sets. Table 5 shows the comparison of DEV of ELM, ECV-ELM and CV-ELM on multiple data sets. Among them, the activation functions of the standard ELM, the CV-ELM model and the ECV-ELM model all use the Sigmoid function. In the ECV-ELM model, the number of training sub models is T = 20, the number of integrated sub models is k = 10, and the subtraining set selected by the submodel is three quarters of the randomly selected training set.
<table><caption>**Table 3.** Comparison of training and testing time of ELM,ECV-ELM,CV-ELM</caption><thead><tr><th colspan="2">Data set</th><th>Cloud</th><th>Strike</th><th>Balloon</th><th>Bodyfat</th><th>House</th></tr></thead><tbody><tr><td rowspan="2">ELM</td><td>Training</td><td>0.0115</td><td>0.0156</td><td>0.0293</td><td>0.0037</td><td>0.0836</td></tr><tr><td>Testing</td><td>0.0193</td><td>0.0107</td><td>0.0168</td><td>0.0034</td><td>0.0755</td></tr><tr><td rowspan="2">ECV-ELM</td><td>Training</td><td>0.4368</td><td>0.4368</td><td>0.6240</td><td>0.3705</td><td>2.002</td></tr><tr><td>Testing</td><td>0.1832</td><td>0.1248</td><td>0.1014</td><td>0.0732</td><td>0.9844</td></tr><tr><td rowspan="2">CV-ELM</td><td>Training</td><td>0.0375</td><td>0.0502</td><td>0.0506</td><td>0.0250</td><td>0.1732</td></tr><tr><td>Testing</td><td>0.0502</td><td>0.0310</td><td>0.0380</td><td>0.0299</td><td>0.1114</td></tr></tbody></table>

<table><caption>**Table 4.** Comparison of testing RMSE of ELM,ECV-ELM,CV-ELM</caption><thead><tr><th rowspan="2">Data set</th><th colspan="3">RMSE</th></tr><tr><th>ELM</th><th>ECV-ELM</th><th>CV-ELM</th></tr></thead><tbody><tr><td>Cloud</td><td>0.0793</td><td><strong>0.0666</strong></td><td>0.0685</td></tr><tr><td>Strike</td><td>0.4950</td><td><strong>0.4511</strong></td><td>0.4640</td></tr><tr><td>Balloon</td><td>0.0939</td><td><strong>0.0897</strong></td><td>0.0898</td></tr><tr><td>Bodyfat</td><td>0.0385</td><td><strong>0.0265</strong></td><td>0.0316</td></tr><tr><td>House</td><td>0.2626</td><td><strong>0.2530</strong></td><td>0.2542</td></tr></tbody></table>

<table><caption>**Table 5.** Comparison of testing DEV of ELM,ECV-ELM,CV-ELM</caption><thead><tr><th rowspan="2">Data set</th><th colspan="3">Dev</th></tr><tr><th>ELM</th><th>ECV-ELM</th><th>CV-ELM</th></tr></thead><tbody><tr><td>Cloud</td><td>0.0067</td><td><strong>8.3245E-04</strong></td><td>0.0054</td></tr><tr><td>Strike</td><td>0.1403</td><td><strong>4.3562E-03</strong></td><td>0.0334</td></tr><tr><td>Balloon</td><td>0.0168</td><td><strong>9.6536E-06</strong></td><td>1.0863E-04</td></tr><tr><td>Bodyfat</td><td>0.0230</td><td><strong>2.8888E-04</strong></td><td>0.0163</td></tr><tr><td>House</td><td>0.0099</td><td><strong>7.9302E-05</strong></td><td>0.0037</td></tr></tbody></table>
# 5 Conclusions

In this paper, a CV-ELM regression algorithm based on ensemble learning is proposed. This algorithm uses ensemble learning method to overcome the disadvantage of poor model stability caused by CV-ELM random input weight, bias and incomplete noise elimination. Through different regression data sets, the performance of ECV-ELM algorithm is analyzed. It is concluded that although the time cost of ECV-ELM algorithm is longer than that of CV-ELM, the generalization performance and robustness of the algorithm have been greatly improved.

Compared to the CV-ELM algorithm, the ECV-ELM algorithm proposed in this section is trained by training multiple CV-ELM sub models, and then using ensemble learning methods to make multiple CV-ELM learners complementation, thus improving the generalization ability and robustness of the algorithm.

# 6 References

1. Han F, Ling Q H, Huang D S. An improved approximation approach incorporating particle swarm optimization and a priori information into neural networks[J]. Neural Computing & Applications, 2010, 19(2):255-261.

2. Han F, Ling Q H. A New Learning Algorithm for Function Approximation By Incorporating A Priori Information Into Feedforward Neural Networks[J]. Neural Computing & Applications, 2008, 17(5-6):433-439.

3. Huang G B, Zhu Q Y, Siew C K. Extreme learning machine: Theory and applications[J]. Neurocomputing, 2006, 70(1-3):489-501.

4. Zhang Min, Zeng Xinmiao, Ma Changchun. An online learning algorithm based on cluster-based extreme learning machine[J]. Computer Engineering and Applications, 2014, 50(11): 188-191.

5. He B, Xu D, Rui N, et al. Fast Face Recognition Via Sparse Coding and Extreme Learning Machine[J]. Cognitive Computation, 2014, 6(2):264-277.

6. An L, Yang S, Bhanu B. Efficient smile detection by Extreme Learning Machine[J]. Neurocomputing, 2015, 149(PA):354-363.

7. Zhou Z H, Zhao J W, Cao F L. Surface reconstruction based on extreme learning machine[J]. Neural Computing & Applications, 2013, 23(2):283-292.

8. Mei-yi LI, Wei-biao CAI, Qing-shuai SUN. Extreme Learning Machine For Regression Based on Condition Number and Variance Decomposition Ratio[C]// 2018 International Conference on Mathematics and Artificial Intelligence.

9. ZENG Lin, ZHANG Xiaogang, BU Zhanyu et al. Extreme learning machine based on principal components estimation[J]. CEA, 2016, 52(4): 110-114.

10. Hansen L K, Liisberg C, Salamon P. Ensemble methods for handwritten digit recognition[C]// Neural Networks for Signal Processing. IEEE, 2002:333-342.

11. Tang W, Zhou Z H. Bagging-Based Selective Clusterer Ensemble[J]. Journal of Software, 2005, 16(4).