Downsizing and enhancing broad learning systems by feature augmentation and residuals boosting

Recently, a broad learning system (BLS) has been theoretically and experimentally confirmed to be an efficient incremental learning system. To get rid of deep architecture, BLS shares the same architecture and learning mechanism of the well-known functional link neural networks (FLNN), but works in broad learning way on both the randomly mapped features of original features of data and their randomly generated enhancement nodes. As such, BLS often requires a huge heap of hidden nodes to achieve the prescribed or satisfactory performance, which may inevitably cause both overwhelming storage requirement and overfitting phenomenon. In this study, a stacked architecture of broad learning systems called D&BLS is proposed to achieve enhanced performance and simultaneously downsize the system architecture. By boosting the residuals between previous and current layers and simultaneously augmenting the original input space with the outputs of the previous layer as the inputs of current layer, D&BLS stacks several lightweight BLS sub-systems to guarantee stronger feature representation capability and better classification/regression performance. Three fast incremental learning algorithms of D&BLS are also developed, without the need for the whole re-training. Experimental results on some popular datasets demonstrate the effectiveness of D&BLS in the sense of both enhanced performance and reduced system architecture.


Introduction
As everyone may know well, over the past decade, deep learning systems have earned great success in many application fields [1,2] and have been drawing more and more attentions in today's academic and industrial communities.Typical deep learning systems include deep restricted Boltzmann machines (RBM) [3], deep belief networks (DBN) [4,5] and deep convolutional neural networks (CNN) [6].Deep learning systems exhibit their overwhelming performance by taking their complicated deep structures and deep learning methods.However, they often suffer from very time-consuming training because a huge number of hyperparameters in the corresponding complicated deep structures are involved.Besides, updating their deep structure becomes an extraordinarily tedious task, due to the need for the whole re-training.To get rid of complicated deep architecture and effectively avoid very time-consuming training (and even re-training), as an alternative to deep learning systems, broad learning system (BLS in brevity) [7] has been recently invented to organize the well-known functional link neural networks (FLNN) [8,9] in broad learning way.Up to date, several variants of BLS have been developed.A typical example is fuzzy BLS [10] which is based on Takagi-Sugeno-Kang fuzzy systems [11][12][13].Theoretical and experimental evidences have revealed BLS's excellent classification/regression performance, with fast wide learning capability.
Since FLNN has both its simple structure of single-layer feedforward neural networks and its fast learning which can be realized quickly using its analytical solution, BLS takes FLNN as its basic structure and then expands it in broad learning way.In contrast to FLNN, BLS has two distinc-tive discrepancies in the structure sense: (1) BLS transforms the set of the original input features of data into randomly mapped features and then take sparse autoencoders (SAE) to optimize them into sparse mapped features [14] so as to capture the intrinsic correlations of the input features.(2) BLS generates the enhancement nodes using nonlinear activation functions with randomly generated weights on all the sparse mapped features.All the sparse mapped features and the enhancement nodes are linked to the output layer, where the weights between the hidden nodes and the output layer are determined analytically by ridge regression [15], without the need for iterative updates of the weights.More importantly, without the need for re-training the whole network, BLS can train its structure in broad learning way for the incremental increase of randomly mapped features, enhancement nodes and input data, thereby avoiding very time-consuming training especially for big data.BLS has been theoretically proved to be a universal approximator [16].A large amount of experimental evidences has indicated that BLS indeed has at least comparably and even much more satisfactory classification/regression performance and strong generalization capability, in contrast to the corresponding deep learning system.However, because the mapped features and enhancement nodes are randomly generated, BLS often needs a huge heap of enhancement nodes to achieve the prescribed performance, which may inevitably cause both overwhelming storage requirement and overfitting phenomenon, therefore, how to downsize a BLS system and simultaneously keep the strong capability of the whole system is becoming an urgent demand.
In this study, we attempt to tackle with this issue by stacking several lightweight BLS sub-systems with both feature augmentation and residuals boosting into its stacked structure called D&BLS.The basic idea can be clearly stated as follows.
1. To avoid overfitting phenomenon of the general BLS and simultaneously downsize its structure, D&BLS is proposed to stack several lightweight BLS sub-systems.As a deep classifier, the whole stacked structure and learning of D&BLS has its novelty in the sense of two joint stacking ways: feature augmentation and residuals boosting, which indeed guarantees enhanced performance of D&BLS.2. Using the bootstrap strategy on feature nodes of each lightweight BLS sub-system, enhancement nodes are generated and hence each lightweight BLS sub-system is built to maintain powerful capability in both fast training and good generalization.3. When D&BLS has its fixed depth, three incremental learning algorithms are designed for three incremental cases to endow D&BLS with the ability of fast learn-ing in broad expansion of each layer without the whole re-training process.4. Experimental results on six classification datasets demonstrate that D&BLS can reduce the number of enhancement nodes of BLS to about 76% or less with the same or better performance.On ten regression datasets, D&BLS can achieve smaller error with comparable number of enhancement nodes.On the MNIST dataset, three incremental algorithms of D&BLS can reach the performance comparable to that of D&BLS with one-shot construction and the training time of three incremental algorithms is much less than that of D&BLS with the one-shot construction, In addition, D&BLS only needs much less running time than BLS when they both share the same sizes.In general, D&BLS with some (i.e., from 2 to 4) lightweight sub-systems can obtain promising learning ability.
The rest of this paper is organized as follows.Section "On BLS" reviews BLS in brevity.Section "On D&BLS" proposes the structure of D&BLS, gives the learning algorithm of D&BLS, and discusses its computational complexity.Besides, three incremental algorithms of D&BLS are also designed in this section.Experimental results about image/non-image dataset classification, regression and incremental learning are presented in section "Experimental results".And the last section concludes this paper.

On BLS
In this section, we give a brief review about BLS's structure and three incremental learning algorithms.

Framework of BLS
In [7], BLS is conceptually composed of four basic parts: (1) original input; (2) feature nodes; (3) enhancement nodes; and (4) output layer.BLS begins with the original input as the input to certain feature mapping algorithms [7], which then generates the feature nodes by taking the sparse autoencoder to slightly fine-tune the random features.As such, all the mapped feature nodes actually form an efficient representation way of the original data.Accordingly, BLS generates a series of enhancement nodes by consecutively acting certain activation function with random parameters on all the mapped feature nodes, then all mapped features generated an enhancement node.Finally, BLS links both feature nodes and enhancement nodes to the output layer with each output weight.As one of the most famous merits, all the output weights can be trained analytically by the well-known ridge regression method [15].Without loss of generality, here we only consider the multi-input and single-output case.Suppose the training dataset X contains the total T original input samples with D dimensions and its output vector Y corresponds the total T target outputs, i.e., X ∈ R T ×D , Y ∈ R T ×1 .We assume BLS has the total n groups of feature nodes in which the ith group of feature nodes takes the mapping function φ i , i 1, 2, . . ., n and contains the total d i feature nodes and the total m groups of enhancement nodes in which the jth group of enhancement nodes takes the activation function ξ j , j 1, 2, . . ., m and contains the total q j enhancement nodes.As such, the ith group of feature nodes may be expressed as: where the weight matrix W ei ∈ R D×d i and the bias matrix β ei ∈ R T ×d i are randomly generated according to a certain distribution, such as normal distribution.According to [7], sparse autoencoder is taken to slightly fine-tune the random feature matrix M i into sparse and compact feature matrix.
Using alternating direction method of multiplier (ADMM), the weight matrix tuned by sparse autoencoder can be solved by the following iterative steps: where Mi is the desired sparse feature matrix, Ŵ is the sparse autoencoder solution, Ŵ 0 , o 0 and u 0 are zero matrices, λ is the given regularization parameter.In addition, r > 0 and S is the soft thresholding operator which is defined as: After a certain number of iterations, the sparse feature matrix can be obtained by: We collect all the groups of feature nodes to form the group M n [M 1 , M 2 , . . ., M n ], which will be used to generate the enhancement nodes.Similarly, the jth group of enhancement nodes may be expressed as: where the weight matrix W h j ∈ R ( n i 1 d i )×q j and the bias matrix β h j ∈ R T ×q j are also randomly generated according to a certain distribution.We collect all the groups of enhancement nodes E j to form the group In this way, the output vector of such a BLS becomes: ( 6) Using ridge regression [15], the pseudoinverse of U m n can be solved as follows: where λ is the given regularization parameter, h n i 1 d i + m j 1 q j is the total number of hidden nodes of BLS.As a result, such a BLS is built quickly, without the need for iterative updates at a very slow learning speed, likewise in BP neural network learning [18].

Incremental learning of BLS
Except for fast training of BLS, another famous merit exists in its simple yet fast incremental learning methods, which indeed reveals its very applicability in online scenarios and/or practical situations where the more promising performance is expected.According to [7], without re-training the whole network structure, incremental learning of BLS can be realized from three cases.

Increment of enhancement nodes
In this case, BLS is expanded by adding the (m + 1)th group of enhancement nodes.That is to say, we have After denoting U ex E m+1 , we can decompose the pseudoinverse of U m+1 n as, where Therefore, the new output weight vector Since W m n , U m n and U m n † have been obtained in advance, W m+1 n can be quickly obtained according to Eqs. ( 12) and (13).

Increment of feature nodes
In this case, BLS is expanded by adding the (n + 1)th group of feature nodes.As such, we have, where are the additional outputs of the enhancement nodes corresponding to the (n + 1)th group of feature nodes.After denoting where A and B can be obtained according to Eq. ( 12).Therefore, the new output weight vector W m n+1 is, Obviously, since W m n , U m n and U m n † have been obtained in advance, W m n+1 can be quickly obtained according to Eqs. (12) and (17).

Increment of input data
In this case, BLS is fed with the additional training data { X ex , Y ex } .The additional feature node M ex generated from X ex is The output of enhancement nodes corresponding to M ex is Please note, both W ei , β ei (i 1,…,n) and W h j , β h j (j 1,…,m) are initialized by taking their corresponding values in Eqs.(1) and (5).Let and U ex M ex |E ex , the pseudoinverse of U m n ex can be decomposed as, where Therefore, the new output weight vector Obviously, since W m n , U m n and U m n † have been obtained in advance, W m n ex can be quickly obtained according to Eqs. (22) and (23).
While BLS shares the fast training originating from the flexible yet random construction of both feature nodes and enhancement nodes, it indeed requires a large amount of feature nodes and especially enhancement nodes to achieve satisfactory performance, which may cause an overwhelming storage requirement and overfitting phenomenon.For example, for MNIST dataset, according to [7], the built BLS needs 11,000 enhancement nodes.Therefore, in the next section, we will deepen BLS by stacking its structures into the proposed classifier so as to downsize the structure of BLS and simultaneously share the promising advantages of BLS.

On D&BLS
In this section, by means of feature augmentation, we first state the deep structure of D&BLS and its learning algorithm and then derive three incremental learning algorithms of D&BLS for three incremental cases.Their computational complexities are also discussed.

Structure and learning of D&BLS
Given D&BLS with L layers, the original training dataset X ∈ R T ×D and its target output vector Y ∈ R T ×1 .The parameters of the ith lightweight sub-system include the total n i groups of feature nodes in which each group contains the total d i feature nodes, the number of the selected feature nodes p i that are chosen to generate enhancement nodes, and the total m i groups of enhancement nodes in which each group contains the total q i enhancement nodes.For convenience, we use Figure 2 illustrates its deep structure.According to Fig. 2, D&BLS essentially stacks several lightweight BLS subsystems into its deep structure simultaneously by two ways: (1) augmenting the original input space with the outputs of previous layer as the inputs of current layer; (2) boosting the residual outputs between the original outputs and the outputs of all the previous layers as the target outputs of current layer.According to stacked generalization principle [17], the use of feature augmentation guarantees enhanced performance of D&BLS, whereas the use of residual outputs maintains good local approximation of current BLS sub-system for subtle output differences remaining in current layer.In general, a lightweight BLS sub-system takes not so many feature nodes and enhancement nodes, in contrast to the BLS system for the same training dataset.To stack lightweight BLS sub-systems well, the deep structure of D&BLS embodies our three reliable considerations as follows.
(1) According to the adopted feature augmentation, the ith lightweight sub-system has its augmented training dataset X i , i.e., X i X, i 1 where Y i−1 corresponds to the augmented feature and denotes the output vector of the (i − 1)th sub-system.In other words, the augmented training dataset X i contains the output information of previous sub-system.Please note, each system always has only one augmented feature whose values are the outputs of previous layer in D&BLS.
(2) Inspired by the idea of residuals boosting [19][20][21], each target output of sub-system in current layer is taken to be the residual value between the corresponding target output and the corresponding outputs of the sub-systems of all the previous layers.As such, the ith sub-system has its target output vector Y i : Obviously, according to stacked generalization principle, the above stacking way of both feature augmentation and Fig. 2 Structure of D&BLS residuals boosting can surely enhance the performance of D&BLS with the increase of the number of sub-systems.In other words, with the prescribed performance, a downsized structure can be expected in the sense of the total number of hidden nodes, in contrast to the general BLS.
(3) When combining all the lightweight sub-systems together, diversity between them should be emphasized to achieve an effective ensemble.To do so, D&BLS only randomly boostraps [22] some feature nodes to generate an enhancement node.Suppose the ith sub-system selects p i feature nodes to generate an enhancement node.After denoting . ., M c p i as the random selected feature nodes, where c 1 , c 2 , . . ., c p i ∈ [1, n i ] denote random integers, then an enhancement node E j can be generated as, Once feature nodes and enhancement nodes of each lightweight BLS system are fixed according to the above strategy, this sub-system can be quickly constructed by the same learning algorithm of BLS as in the last section.
D&BLS begins with the construction of the first lightweight.BLS sub-system on the original training dataset X and its target output vector Y.And then D&BLS obtains the output vector Y 1 of the first lightweight BLS subsystem.After generating X 2 X|Y 1 and Y 2 Y − Y 1 , respectively, as the inputs and target outputs of the second lightweight BLS sub-system, the second sub-system is built in the same way as in the above.This procedure is repeated until the maximum number of sub-systems (i.e., layers) or the prescribed performance of D&BLS is achieved.As a result, the final output vector of D&BLS on the original training dataset can be expressed as, In summary, the learning algorithm is given in Algorithm 1.To visualize the effectiveness of Algorithm 1, a simple twoclass classification dataset contained 200 training samples (100 positive samples) is taken for didactic experimentation.The dataset is generated from the model make_moons of Scikit-learn [29] package with 0 noise, as shown in Fig. 3a.Four D&BLSs with different number of sub-systems are taken and the deeper one is deepened from the shallower one.The deepest one is a four-layer D&BLS: (6 × 6, 9, 1 × 3)-(5 × 6, 9, 1 × 2)-(5 × 6, 9, 1 × 2)-(5 × 6, 9, 1 × 2), its whole structure becomes (126, 9).Here the average results of the ten runs are reported.Figure 3b-e shows the decision boundary of all D&BLSs.
From Fig. 3, the decision boundary of D&BLS becomes more complex to match the distribution of training samples as D&BLS becomes deeper.It indicates that Algorithm 1 is effective and the learning ability of D&BLS can be enhanced by adding more sub-systems with appropriate parameters.
The computational complexity of Algorithm 1 can be analyzed in terms of D&BLS with L sub-systems.Without loss of generality, according to the same strategy of BLS in [7], Fig. 3 The moon dataset and corresponding decision boundary of four D&BLSs with different number of sub-systems here we suppose ξ in step 12 is a sigmoid activation function [23] and the number of iterations of sparse autoencoder taken in step 6 is K.
Thus the computational complexity of step 2 and 3 can be calculated to be O(2T ) ≈ O(T ).The computational burden from step 4 to step 7 about the generation of feature nodes takes, In general, d i < T , d i < D and 2 < K can be assured, so the computational complexity mentioned above will be reduced

to O(K T Dn i d i ). Then step 8 requires O(T n i d i ).
The computational burden of the generation of enhancement nodes from step 9 to step 13 takes O(m i (T p i + p i q i + T q i p i + 3T q i )) ≈ O(T m i q i p i ) Then steps 14 and 15 require O(2T m i q i ) ≈ O(T m i q i ).The calculation of the pseudoinverse of U m i n i i in step 16 requires Hence, with the total L sub-systems, the computational complexity of Algorithm 1 becomes As we can see, the most time-consuming parts for each sub-system in D&BLS are the generation of feature nodes, the generation of enhancement nodes and the calculation of pseudoinverse.Similarly, we can calculate the computational complexity of BLS for comparison.Suppose BLS contains d × n feature nodes and q × m enhancement nodes, its computational complexity is For easy comparison, here we suppose each sub-system of D&BLS has the total n/L ×d feature nodes and m/L ×q enhancement nodes.As such, the computational complexity of Algorithm 1 becomes In general, 1 L L i 1 p i < nd and (nd + mq) ≤ T can be assured.Thus the computational complexity is obviously less than that of BLS.

Incremental learning for D&BLS
To achieve enhanced performance, it seems that we may expand D&BLS by increasing a lightweight BLS subsystem successively.However, too many lightweight BLS sub-systems often cause overfitting phenomenon.Our experimental evidences demonstrate that the total number of layers (i.e., sub-systems) should generally take a small integer (i.e., from 2 to 4).Therefore, a feasible strategy is to develop its incremental algorithms with the fixed number of layers.What is more, we should still consider how to expand D&BLS for the case of incremental input data.Here we develop three incremental algorithms of D&BLS for three incremental cases, without the need of re-training the whole classifier.

Increment of enhancement nodes
In this case, each sub-system in D&BLS is expanded by adding m (i.e., a small integer) groups of enhancement nodes in which each group contains q enhancement nodes.
The additional enhancement nodes are generated from a portion of the original feature nodes.After denoting n i |U ex can be calculated by Eqs. ( 11) and (12).Then the output weight vector can be determined by and hence the output of this sub-system after expansion becomes . The incremental algorithm for this case is summarized in Algorithm 2.
Below let us analyze the computational complexity of Algorithm 2 in terms of L sub-systems.Because the training of each sub-system in this incremental case only deals with both the generation of the incremental enhancement nodes and the calculation of the corresponding pseudoinverse, we mainly observe the corresponding steps in Algorithm 2.
That is to say, the computational complexity from step 3 to step 7 takes O(T m q p i ), while steps 8 and 9 require O (2T m q) ≈ O(T m q).What is more, because the number of the additional enhancement nodes is generally a small inte-ger, i.e., m q < (n i d i + m i q i ) and m q < T , the calculation of the pseudoinverse Obviously, step 11 and step 12 require O (2T (n i d i + m i q i + m q)) ≈ O(T (n i d i + m i q i )).Hence, with the total L iterations, the computational complexity of Algorithm 2 becomes If we re-train the whole model based on Algorithm 1 and assume (n i d i + m i q i ) < T is assured, its computation complexity will become Obviously, increment of enhancement nodes can save much training time.

Increment of feature nodes
In this case, after n (i.e., a small integer) groups of feature nodes with d feature nodes are added, each sub-system in D&BLS is further expanded by the corresponding additional m (i.e., a small integer) groups of enhancement nodes with q enhancement nodes.
After the total n groups of the additional feature nodes M ex are generated from X i , the total m groups of the additional enhancement nodes E ex are accordingly generated from a portion of the additional feature nodes M ex .After denoting U ex M ex |E ex , the pseudoinverse of U m i n i + n U m i n i |U ex can be calculated by Eqs. ( 12) and ( 16).As such, the output weight vector can be obtained by . The incremental algorithm for this case is summarized in following Algorithm 3.
The computational complexity of Algorithm 3 can be analyzed in terms of L sub-systems.Obviously, the training of each sub-system in this incremental case mainly has three stages, namely, the generation of the additional feature nodes and the corresponding additional enhancement nodes and the calculation of the pseudoinverse.According to Algorithm 3, the computational burden from step 4 to step 7 about the generation of feature nodes takes O K T D n d .Step 8 requires O T n d .The computational complexity from step 9 to step 13 about the generation of enhancement nodes is O(T m q p i ).Then steps 14 and 15 require O(2T m q) ≈ O(T m q), step 16 requires O T n d + m q .Because the number of both the additional feature nodes and the additional enhancement nodes is generally a small integer, i.e., n d + m q < (n i d i + m i q i ) and (n i d i + m i q i ) < T , and the calculation of the pseudoinverse What is more, step 18 and step 19 require . In summary, with the total L iterations, the computational complexity of Algorithm 3 becomes If we re-train the whole structure based on Algorithm 1 and suppose (n i d i + m i q i ) < T is assured, its computation complexity will become Obviously, increment of feature nodes can save much training time.

Increment of input data
In this case, D&BLS is expanded by adding . Suppose that {X T + T , Y T + T } denotes the all input data after incremental operation.After the total n i groups of the additional feature nodes M ex with d i feature nodes are generated from X T , the total m i groups of the additional enhancement nodes E ex with q i enhancement nodes are generated from a portion of the additional feature nodes M ex .After denoting U ex U ex can be calculated by Eqs. ( 21) and (22).As a result, the output weight vector can be obtained by The incremental algorithm for this case is summarized in Algorithm 4.
The computational complexity of Algorithm 4 can also be analyzed in terms of L sub-systems.The computational complexity of each sub-system in this incremental case mainly contains three stages: the generation of the additional feature nodes and the additional enhancement nodes and the calculation of pseudoinverse.According to Algorithm 4, the computational complexity from step 4 to step 7 about the generation of feature nodes takes O K T Dn i d i .Step 8 requires O T n i d i .Then the computational burden from step 9 to step 13 about the generation of enhancement nodes is O T m i q i p i .Step 14 and 15 require O 2 T m i q i ≈ O T m i q i , then step 16 requires O T (n i d i + m i q i ) .Suppose What is more, steps 18 and 19 require O 2 T + T (n i d i + m i q i ) ≈ O T + T (n i d i + m i q i ) .Suppose T < T and K D < T are assure, with the total L iterations, the total computational complexity of Algorithm 4 can be described as If we re-train the whole structure based on Algorithm 1 and suppose (n i d i + m i q i ) < T + T is assured, its computation complexity will become Obviously, increment of input data can save training time.Finally, let us state how to set appropriate parameters in the above algorithms.As pointed out in the above, the total number of layers generally take a small integer (i.e., from 2 to 4).Besides, according to our extensive experiments, the weight matrix W and the bias matrix β involved in hidden nodes generally are drawn from standard normal distributions, respectively.Since the number p i of selected feature nodes should generally be much less than the total of feature nodes, it is set to be at most 75% of the total number of feature nodes.As for the regularization parameters λ and r involved in sparse autoencoder, they are simply set to be 0.001 and 1, respectively.The regularization parameter λ involved in ridge regression is simply set to be 2 −30 .In our experiments, the number of additional hidden nodes for each sub-system is simply taken to be 10 and the number of additional input data is simply set to be 3000.

Experimental results
In this section, to evaluate the performance of D&BLS, we organize four groups of experiments.The first group about classification is carried out on both five image datasets and 1 UCI [24] classification datasets, the second about regression on ten UCI regression datasets, the third about incremental learning on popular datasets MNIST [25], and the last about running time is compared between BLS and D&BLS.All the experiments are carried out under a computer that equips Intel-i3 3.40 GHz CPU, 4 GB memory.The experiments of LSSVM are carried out in MATLAB environment and the other experiments are carried out in Python environment.Before the experiments, all the datasets are divided into training and testing set.
The comparative methods include BLS [7], support vector machine (SVM) [26], least squared SVM (LSSVM) [27].Because the regularization parameter C and the kernel parameter γ of SVM play an essential role in performance, they need to be chosen appropriately for a fair comparison.In our experiments, we determine them by taking a grid search from {2 −24 , 2 −23 , . . ., 2 24 , 2 25 } for the parameters (C, γ ) in SVM, while the parameters (C, γ ) of LSSVM are decided using LS-SVMlab Toolbox.
In our experiments, the original codes of BLS in [7] (classification version) and in [16] (regression version) are taken.As for SVM, its classification version SVC and its regression version SVR are taken from Scikit-learn [29] package.
As for BLS and D&BLS, to make a fair comparison between them, the weight matrix W and the bias matrix β are drawn from standard normal distributions and the tansig functions are taken as the activation functions of enhancement nodes.All other parameters involved in BLS and D&BLS keep the same and their settings have been stated in the last section.

Classification
This group of experiments consists of two parts.The first part is to observe the performance of D&BLS and the comparative methods on image data, while the second part deals with the 1 UCI non-image classification datasets.Five popular image datasets are MNIST [25] and USPS [29] about handwritten digits, COIL20 [30] and COIL100 [31] about 3-D objects, and extended YaleB [32] about human faces.A UCI classification dataset is isolet [33] about spoken letter recognition.To obtain good performance, we use "one vs. one" strategy for LSSVM when faced with multi-class problems and thus omit its parameters setting here.
According to the accuracies obtained by BLS on the adopted datasets, we determine D&BLS by the trial and error strategy.This involves the number of lightweight BLS subsystems in D&BLS.To make our comparison fair, here we report average results (i.e., mean and standard deviations) for ten trials on each dataset.

MNIST
MNIST contains 70,000 handwritten digits of ten classes.Every digit is represented by an image with the size of 28 × 28 gray-scaled pixels.60,000 out of all the handwritten dig-its are partitioned into the training dataset and the remaining 10,000 handwritten digits into the testing dataset.Here a subset of MNIST is chosen for our experiment with the first 10,000 training images and all the 10,000 testing images.
In our experiments, we perform a grid search for the optimal parameters of BLS, including the size of feature node groups, the number of feature node groups and the number of enhancement nodes from [10,30] × [10,20] × [100, 5000] with the steps being 1, 1 and 100 due to the large searching range of enhancement nodes.The parameters setting can be checked in Table 1 and the classification accuracies of all the adopted methods are given in Table 2.

USPS
US Postal Service (USPS) handwritten digits recognition corpus is another well know dataset, which contains 9298 gray images of size 16 × 16 pixels, with ten class digit 0-9.We randomly select 7500 images for training and the rest 1798 images for testing.In our experiment, the search range for BLS is [10,30] × [10,20] × [100, 5000] with the steps being 1, 1 and 100.The parameters setting is shown in Table 1 and the classification accuracies are listed in Table 2.

COIL20
COIL20 contains 1440 gray images of 20 different 3-D objects with 128 × 128 pixels.The images of each object were taken 5 degrees apart as the object is rotated on a turntable and each object has 72 images.There are 36 images of each object for training and the remaining 720 images for testing.In our experiment, the search range for BLS is [10,30] × [10,20] × [100, 3000] with the steps being 1, 1 and 100.The parameters setting are given in Table 1 and the classification accuracies can be checked in Table 2.
Table 1 Parameters setting of the adopted methods on the classification datasets

123
parameters setting are listed in Table 1 and the classification accuracies are shown in Table 2.
As we can see, D&BLS exceeds the accuracies of all the adopted methods with much less enhancement nodes (i.e., only 2710 nodes), which is about 61.59% of the number of enhancement nodes in BLS.

Extended YaleB
Extended YaleB face consists of 2414 cropped images of 38 objects with the size of 32 × 32 pixels.The images contain large variations in terms of illumination conditions and expressions for each subject.There are 30 images of each person for training and the remaining 1274 images for testing.Following [16], BLS contains 60 × 30 feature nodes and 1 × 6000 enhancement nodes.The parameters setting of other models are shown in Table 1 and the classification accuracies of all the adopted methods are listed in Table 2.
From Tables 1 and 2, D&BLS matches the accuracy (97.92%) of BLS with only 232 enhancement nodes, which is about 3.87% of the number of enhancement nodes in BLS.D&BLS also exceeds the accuracy of SVM (89.32%) and LSSVM (89.95%).

UCI
To verify the effectiveness of D&BLS for non-image data, we choose the isolet dataset here, which contains 1560 spoken letter recognition samples of 26 classes with 617 features in each sample.We randomly choose 1092 samples for training and the rest 468 samples for testing.In our experiment, the search range for BLS is [10,30] × [10,20] × [100, 3000] with the steps being 1, 1 and 100.The parameter details are shown in Table 1 and the classification performance are given in Table 2.
D&BLS indeed gets the best result (95.49%) among all the adopted methods and only needs about 76% of the number of enhancement nodes in BLS, which actually demonstrates the effectiveness of D&BLS for non-image samples.

Regression
In this subsection, we attempt to verify that D&BLS can achieve smaller regression errors than BLS when they both have comparable number of enhancement nodes.We take ten UCI regression datasets, whose details are listed in Table 3 Following [16], the same parameters setting of BLS are taken by a grid search from [1,10] × [1,30] × [1,200] with step being 1 and the best experimental results about the rootmean-square errors (RMSE) from ten trials are taken to make a fair comparison.The detailed parameter settings and the experimental results by SVM, LSSVM, BLS and D&BLS on these datasets are given in Tables 4 and 5, respectively.Clearly, D&BLS performs comparably to BLS on the datasets Bodyfat and Weather Izmir and better than BLS with comparable number of enhancement nodes in the remaining eight datasets.It can be concluded that with almost the same sizes of the architectures, D&BLS attractively outperforms BLS in the sense of testing accuracy on these datasets, which actually hints that D&BLS may take smaller sizes of the architectures to achieve comparable regression performance to BLS.Compared to SVM and LSSVM, the D&BLS also achieves the best performance.

Incremental algorithms
In this section, we divide it into three cases to observe the average performances (with their standard deviations) of incremental learning algorithms of D&BLS for ten runs on the dataset MNIST: (1) increment of enhancement nodes (corresponding to Algorithm 2); (2) feature nodes (corresponding to Algorithm 3); and (3) input data (corresponding to Algorithm 4).
Our experiments begin with a three-layer D&BLS, whose whole structure is (696,275) consisting of (30 × 8, 10, 1 × 100)-(29 × 8, 10, 1 × 95)-(28 × 8, 10, 1 × 80) on the MNIST dataset where the first 10,000 samples from the total 60,000 training samples are taken as the initial training samples, and the same 10,000 testing samples are still kept.The parameter p i of the additional enhancement nodes is also set to be 10.Besides, the testing accuracy, training time and testing time of D&BLS with one-shot construction are also provided to highlight the advantage of the proposed three incremental algorithms.

Increment of enhancement nodes
Two groups of experiments are arranged to observe Algorithm 2 for this case.The first group considers the total 4  times of incremental operations.That is, it takes the addition of ten enhancement nodes in an incremental way.The second group keeps the same except 12 rather 10 enhancement nodes.Thus, the whole network with three sub-systems has 30 or 36 more enhancement nodes after enhancement nodes are increased in a step-by-step way.The results are given in Table 6.
As we can see, the testing accuracy continues to be improved as the additional enhancement nodes are inserted into D&BLS in an incremental way and finally reaches the performance comparable to that of D&BLS with the oneshot construction (see the last row in Table 6).Obviously, the training time of Algorithm 2 is much less than that of D&BLS with the one-shot construction, e.g., 1.8399 s vs. 5.5217 s.

Increment of feature nodes
Two groups of experiments are arranged to observe Algorithm 3 for this case.The first group considers the total 4 times of incremental operations.That is, it takes the addition of ten feature nodes and ten corresponding enhancement nodes in an incremental way.The second group keeps the same except 12 rather 10 feature nodes and enhancement nodes.All the results can be checked in Table 7.
It can be seen that the testing accuracy continues to go up as the additional feature nodes are inserted into D&BLS in an incremental way and finally reaches the performance comparable to that of D&BLS with the one-shot construction (see the last row in Table 7).Obviously, the training time of Algorithm 3 is much less than that of D&BLS with the oneshot construction, e.g., 3.3657 s vs. 6.4858s.

Increment of input data
Two groups of experiments are arranged to observe Algorithm 4 for this case.The first group considers the total 4 times of incremental operations.That is, it takes the addition of 3000 training samples each time in an incremental way from all the remaining training samples.The second group keeps the same except 8000 rather than 3000 training samples.The incremental results are shown in Table 8.
It can be observed that the testing accuracy has a slight decline (e.g., from 94.18 to 93.54%) when the additional input data are fed in the first time into D&BLS, while the testing accuracy begins to go steadily up as more input data are fed in an incremental way and finally reaches the performance comparable to that of D&BLS with the one-shot construction (see the last row in Table 8).Obviously, the training time of Algorithm 4 is much less than that of D&BLS with the one-shot construction, e.g., 3.6203 s vs. 10.4513s.

Running time comparison
In this subsection, we will show D&BLS's another attractive property-it only needs much less running time (i.e., training time and testing time) than BLS when they both share the same sizes (i.e., feature nodes and enhancement nodes).The COIL-100 and the USPS datasets are taken for experiments.For a fair comparison, here we report average results and standard deviations for ten trials on each dataset.D&BLS takes a three-layer structure and each sub-system keeps the same parameters setting, especially the number of selected feature nodes is 1.In other words, if a BLS has 30 × 10 feature nodes and 1 × 2100 enhancement nodes, each subsystem in the corresponding D&BLS will contain 10 × 10 feature nodes, 1 selected feature node, and 1 × 700 enhancement nodes.The obtained experimental results are given in Table 9.
According to Table 9, D&BLS indeed needs much less running time than BLS.What is more, D&BLS runs much more quickly than BLS with the increasing number of enhancement nodes.This tendency becomes considerably surprising in the case of 3600 enhancement nodes for the USPS dataset.That is, BLS occupies 27.2773 s while D&BLS occupies only 9.8075 s for the USPS dataset.

Conclusion
While the recently developed broad learning system BLS has been exhibiting its promising performance without the need of deep structure, it often requires a huge heap of hidden nodes, which may perhaps inevitably cause both overwhelming storage requirement and overfitting phenomenon.In this paper, D&BLS is proposed to synthesize both deep and broad learning so as to guarantee enhance performance, downsize the structure and reduce the running time of the general BLS.D&BLS consists of several lightweight BLS sub-systems.D&BLS has its structure novelty in its joint stacking way of feature augmentation and residuals boosting.The whole learning algorithm and three incremental algorithms of D&BLS are proposed and then verified by our experiments about image datasets and UCI datasets.
Since the parameter settings of each sub-system in D&BLS are given only by the trial and error strategy, a deeper D&BLS sometimes cannot assure better performance than a shallower D&BLS.As a result, how to determine the structure of D&BLS in a moderate time and how to design a deeper D&BLS with desirable ability for practical application scenarios are our research works in near future.

Fig. 1
Fig. 1 An illustrative structure of BLS while step 17 and step 18 require

COIL100COIL100
contains 7200 color images of 100 different 3-D objects with 128 × 128 pixels.The images of each object were taken 5 degrees apart as the object is rotated on a turntable and each object has 72 images.In this experiment, each image was resized into 32 × 32 and converted into gray image.We randomly select 5000 images for training and the rest 2200 images for testing.In our experiment, the search range for BLS is[10,30] × [10, 20] × [100, 5000] with the steps being 1, 1 and 100.The values are the best one among all the results

Table 3
Details of the UCI regression datasets

Table 4
Parameters setting of the adopted methods (1) on the UCI regression datasets and (2) on the UCI datasets for regression

Table 5
RMSEs of the adopted methods on the UCI datasets for regression

Table 6
Experimental results of the MNIST dataset using increment of enhancement nodes

Table 7
Experimental results of the MNIST dataset using increment of feature nodes