NC State
BioResources
Gao, Y., Hua, J., Chen, G., Cai, L., Jia, N., and Zhu, L. (2018). "Prediction of fiber quality using refining parameters in medium-density fiberboard production via the support vector machine algorithm," BioRes. 13(4), 8184-8197.

Abstract

Fiber quality greatly influences the performance of medium-density fiberboard (MDF). To evaluate the fiber quality more accurately during refining, a novel quantitative parameter-property relationship model was developed based on the support vector machine (SVM) algorithm. Based on the mill production conditions, a total data set of 1173 experimental fiber quality data points under a wide range of five refining parameters was employed to train and verify the model. By comparing the effectiveness between the model using nonlinear SVM and the model based on multiple linear regression (MLR), the values of mean absolute error (MAE), mean relative error (MRE), root mean square error (RMSE), and Theil’s inequality coefficient (TIC) were reduced 92.19%, 92.36%, 87.29%, and 87.21%, respectively. The results showed that the performance of the predictive model developed using SVM was superior to the MLR model. Furthermore, the variations of the percentage of qualified fibers with each production parameter were predicted using the established model. The prediction model that resulted can be applied to predict the fiber quality during the refining process in an MDF production mill.


Download PDF

Full Article

Prediction of Fiber Quality Using Refining Parameters in Medium-density Fiberboard Production via the Support Vector Machine Algorithm

Yunbo Gao,a Jun Hua,a,* Guangwei Chen,a Liping Cai,b Na Jia,a and Liangkuan Zhu a

Fiber quality greatly influences the performance of medium-density fiberboard (MDF). To evaluate the fiber quality more accurately during refining, a novel quantitative parameter-property relationship model was developed based on the support vector machine (SVM) algorithm. Based on the mill production conditions, a total data set of 1173 experimental fiber quality data points under a wide range of five refining parameters was employed to train and verify the model. By comparing the effectiveness between the model using nonlinear SVM and the model based on multiple linear regression (MLR), the values of mean absolute error (MAE), mean relative error (MRE), root mean square error (RMSE), and Theil’s inequality coefficient (TIC) were reduced 92.19%, 92.36%, 87.29%, and 87.21%, respectively. The results showed that the performance of the predictive model developed using SVM was superior to the MLR model. Furthermore, the variations of the percentage of qualified fibers with each production parameter were predicted using the established model. The prediction model that resulted can be applied to predict the fiber quality during the refining process in an MDF production mill.

Keywords: Fiber quality; MDF; Refining; Predictive model; SVM

Contact information: a: College of Electromechanical Engineering, Northeast Forestry University, Harbin, 150040, China; b: Mechanical and Energy Engineering Department, University of North Texas, Denton, TX 76201, USA; *Corresponding author: huajun81@163.com

INTRODUCTION

Medium-density fiberboard (MDF) has been widely applied in the furniture and interior decoration markets due to its dimensional stability, workability, flatness, smooth appearance, good bond strength, and screw-holding ability (Hua et al. 2012). A crucial step in the fiberboard production is the refining process (Runkler et al. 2003). It is impractical and costly to evaluate which fiber qualities are affected by certain production parameters through experimental determination techniques during the refining process. Production parameters are mainly adjusted based on the experience of workers, which suffers from poor accuracy due to the lack of theoretical guidance on the model of fiber quality relating to production parameters. Therefore, it is essential to develop new models to predict fiber quality according to production parameters.

To investigate the influence of production parameters during refining on the fiber or the fiberboard quality, some studies have been carried out over the past decade. Chen and Hua (2009) developed a constraint relationship between the fiber productivity and fiber quality using third-order polynomial and linear regressions, and an optimization method was developed to adjust the fiber quality by modifying fiber productivity that depends on the feeding screw revolution speed and the opening percentage of the discharge valve. The relationship between the content of bark and the fiber quality was investigated by Jia et al. (2015). It was demonstrated that the fiber screening value first presented an increasing trend and then descended with an increase in bark content. Xing et al. (2006) investigated the influence of thermo-mechanical refining on the properties of MDF panels made from black spruce bark, and the results showed that the preheating retention time was an important factor for both the modulus of rupture and the modulus of elasticity. The steam pressure was an important factor for internal bond strength, modulus of rupture (MOR), and modulus of elasticity (MOE). The previously described studies revealed the relationship between the production parameters during refining and the quality of fiber or fiberboard based on linear or polynomial regression, which was shown to be disadvantaged by low accuracy and the inability to predict the quality online.

To address these drawbacks, some researchers developed models between the production parameters during refining and the quality of fiber or fiberboard with intelligent algorithms. To determine the process set points leading to minimum production cost for the target quality, Gerstorfer et al. (2001) established a Takagi-Sugeno-Fuzzy model for the refining process based on the experts’ knowledge as well as the data collected. The neuro-fuzzy modeling methods were used by Runkler et al. (2003) to model the wood chip refiner process for fiberboard production to provide online predictions of the two important quality indices (flexural strength and water uptake). The results showed that the model accuracies achieved were approximately ± 5 N/mm2 for flexural strength and approximately ± 10%/24 h for water uptake. However, the fuzzy rules were determined on the dependence of experts’ technological experience, which only leads to a certain extent of improvement on the prediction accuracy. Because the artificial neural network (ANN) can model highly nonlinear systems without using complex deduction rules or large data (Huang and Lu 2016), it was used as a predictive method to determine the moisture resistance of particle- and fiberboards under cyclic testing conditions by Esteban et al. (2010). However, the training procedure for ANN models is not only time consuming but it is also possible to get trapped in local minima (Hong et al. 2013).

With the advantages of a simple structure, good generalization ability, nonlinear modeling properties (Wang et al. 2009; Zhou et al. 2016; Sun et al. 2016), and avoidance of the problems of over-fitting, local extremum, and dimension disaster (Shi et al. 2010; Zhao et al. 2014; Sun et al. 2016), the support vector machine (SVM) has become a promising classification and regression algorithm. The SVM can be used for classification of data and text, system modeling and prediction, pattern recognition, anomaly detection, and time series prediction (Jiao et al. 2016) in many fields (Mokhtarzad et al. 2017; Roushangar and Ghasempour 2017; Huang et al. 2018). Among these fields, Zhang et al. (2016) employed the SVM method to establish parameters-properties models in the papermaking field. Although SVM has been employed in many fields because of its advantages of nonlinear relationship expression, it has not been used for the modeling of fiber quality during the refining process in MDF production.

This study is aimed at increasing the percentage of qualified fibers (QF) by adjusting five parameters, i.e., the conveyer screw revolution speed (SR), accumulated chip height (CH), opening ratio of the discharge valve (OV), content of log bark (CB), and content of Chinese poplar (CP) during refining. The QF were evaluated by fiber size, i.e., the screen mesh grades of the fibers. Firstly, a large amount of fiber quality data was collected from the MDF production mill under a wide range of SR, CH, OV, CB, and CP to uncover the relationship between the refining parameters and fiber quality. Secondly, SVM was utilized to construct the nonlinear predictive model of fiber quality during the refining. The experimental values and predicted outputs of the model were compared, and the accuracy of the model was established. Thirdly, the results from the SVM-based model and the multiple linear regression (MLR)-based model were compared, showing that the SVM model predicted fiber quality more effectively and accurately than that from the MLR.

EXPERIMENTAL

Materials

Data collection

The data were collected in a production line at a MDF mill in Northern China and two major species, namely Chinese poplar (Populus lasiocarpa Oliv.) and Chinese larch (Larix potaninii Batalin), were used in the production. In the production line, the 4000 kW refiner (model: 50-ICP; Andritz Group, Graz, Austria) utilized had double 1372-mm diameter disks running at a rotational velocity of 1,500 r/min.

Figure 1 shows the major components of the refiner. A feeding screw (2) transfers the wood chips in the hopper (1) to the pre-heater (3). The preheating retention time was determined by CH. Chips were discharged to the refiner (6) using a conveyer screw (4) after the steam-softening. Through the discharge pipe (5), the refined fibers were unloaded under the steam pressure in the refiner. The opening ratio of the valve installed on the discharge pipe (5) was used to adjust the amount of unloaded fibers.

Fig. 1. The major components of the refiner: 1) Hopper, 2) Feeding screw, 3) Pre-heater, 4) Conveyer screw, 5) Discharge pipe, and 6) Refiner

In the chip washing and steaming processes, the moisture content of the wood chips increased. However, during the chip transportation, the feeding screw squeezed the moisture from the chips in the feeding pipe and brought the final moisture content of the chips to 50%.

In the pre-heating process, the steaming pressures usually slightly fluctuate from 0.766 MPa to 0.990 MPa (corresponding to a steaming temperature that ranged from 168.6 °C to 179.5 °C). The gap between the two refining disks was pre-set to 0.1 mm. The sensors installed on the production line measured five refining parameters, including SR, CH, OV, CB, and CP, hourly under usual production conditions.

Generally, large fiber size would produce panels with poor board appearance, while smaller-sized fibers could cause the reduction in panel strength (Shi et al. 2006). Technically, the good fiber shape for MDF requires a moderate ratio of length/width of the fiber (Chen 2012). Based on the mill’s practices, the fibers with a size between 20 to 120 screen mesh were considered to be qualified in this study.

In the fiber processing, 10 g of fiber were collected and weighed with a balance off the production line to determine the fiber sizes for each measurement. The percentage of the QF in the total amount = (The weight of qualified fiber (g) / 10 g) × 100%.

It was hypothesized that the response variable (QF) during the refining could be estimated by correlating the predictor variables (SR, CH, OV, CB, and CP). A large sample size, 1,173 measurements for each variable, was analyzed in this study. Two models, namely, MLR and SVM, were developed for predicting the fiber quality. The accuracies of the two models were validated and compared using the data collected in the MDF production line.

Methods

Support vector machine algorithm

The SVM, proposed by Vapnik (1999), is a relatively new and promising classification and regression algorithm based on the statistical learning theory and structural risk minimization principle. Based on this principle, SVM possesses an optimum network structure that is beneficial to reduce the global error of the model (Xiao et al. 2014).

The basic principle is as follows (Chu et al. 2017). Training data are presented in the form{xiyi}n i=1, where xi values are the input data, yi values are the corresponding output data, and nis the number of training data points. The SVM is used to search an optimal regression function, which can estimate all of the training data (Drucker et al. 1997). The regression function can be expressed as,

 (1)

where wRn denotes the weight vector, φ(x) denotes the nonlinear mapping function, and b denotes the bias. As mentioned previously, SVM is based on risk minimization, while w and bare estimated by minimizing the regularized risk function as below,

 (2)

where 1/2‖w‖2 is the flatness of the function, C is the penalty factor, which is the correlation between empirical error and flatness of the model (Yan and Shi 2010), ε is a prescribed parameter, and Lε(yif(xi)) is the ε-insensitive loss function that can be defined as:

 (3)

By introducing the slack variables ξ and ξ*, Eq. 2 can be written as:

 (4)

 (5)

The dual objective optimization problem can be solved using the Lagrange multipliers (Ma et al. 2003). Finally, the regression function is obtained as the following equation,

 (6)

where aand a* i are the Lagrangian operators, nsv is the number of support vectors, and K(xiyi) is the kernel function. It is important to select the kernel function and its parameters because the generalization performance of SVM depends on the type of kernel function, its parameters, and several internal parameters of SVM (Zhao et al. 2016).

The Gaussian Radial Basis Function (RBF) is mostly used for the kernel function due to its properties of good generalization and nonlinear forecast as well as its characteristic of few parameters that need to be adjusted (Bishop 1995; Keerthi and Lin 2003). Therefore, in this study, the Gaussian RBF was selected as the kernel function using the following formula,

 (7)

where σ is the width of RBF.

Therefore, there are two variables that need to be selected in the SVM model, which are the constant “C” and the width of the Gaussian RBF kernel “σ”. In this study, the optimization of these parameters was performed by a systematic grid search of the parameters using the cross-validation on the training set.

The SVM predictive model for fiber quality

In this paper, the software used for analysis was Matlab (MathWorks, R2010a, Natick, MA, USA). The diagram of the SVM predictive model is illustrated in Fig. 2. The details of the operation process for the SVM model are discussed below.

Step 1: Data preprocessing

To ensure the training stability of the SVM and avoid the bad influence caused by discrepancy of quantitative dimension, the data from experiments were normalized with the following mapping function,

 (8)

where xis the normalized data, x is the original data, and xmax and xmin denote the maximum and minimum raw input values, respectively. The original data were normalized to the range of 0 to 1.

Step 2: Cross validation (CV) to select the best parameters of regression, C, and σ

First, the fitness function (mean square error) was determined based on 3-CV, ε was pre-set as 10-4, and the range of C and σ were defined and meshed. Second, the fitness function was recalculated through updating C and σ within the grid. Finally, the C and σ that generate the minimum average mean square error (MSE) of the three models were selected as the best parameters.

Fig. 2. Diagram of the support vector machine predictive model

Step 3: Building the SVM model for fiber quality

The SVM model for fiber quality was trained and established based on the best parameters obtained in step 2, which can be used to investigate the relationships between the parameters and the property.

Step 4: SVM regression prediction

The built SVM model was used to predict the training set and test set, and then the fiber quality under prediction was compared with the experimental data.

Step 5: Evaluation of predictive performance

The predicted performance was evaluated in terms of mean absolute error (MAE), mean relative error (MRE), root mean square error (RMSE), and Theil’s inequality coefficient (TIC). They are defined according to the following formulas, respectively,

 (9)

 (10)

 (11)

 (12)

where yi are the actual outputs (experimental qualified fibers),  the outputs of models (predicted qualified fibers), and n is the number of compounds in the analyzed data set.

RESULTS AND DISCUSSION

Results of the SVM Algorithm

The inputs for the fiber quality model using the SVM algorithm were SR, CH, OV, CB, and CP, and the output was QF. Out of a total set consisting of 1173 groups of data from the experiment, 887 groups of data were used to build the model as a training set, while 286 groups of data were used to validate the model established as a test set. In selecting the SVM parameters, K-CV can effectively avoid the occurrence of overlearning and under-learning. In this paper, the optimization of the SVM parameters was carried out by 3-CV. Figure 3 illustrates the SVM parameter selection process and results (contour map and 3D view). As shown in Fig. 3a, contour lines represent the C and σ corresponded MSE based on the K-CV method. The most optimum condition where the MSE of the 3-CV method equaled 0.00025373 was selected. The final optimization results were C = 1 and σ = 5.6569.

Fig. 3. SVM parameter selection process and results: (a) contour map, and (b) 3D view

To estimate the accuracy of the fiber quality model, the experimental data were compared with the predicted outputs as shown in Fig. 4. In Fig. 4a, the scatter plots of the predictive versusthe experimental outputs for the qualified fibers were observed. If the model exactly matched the actual values, all data points would be on the main diagonal. It was revealed that the data points were indeed close to the main diagonal, indicating good model accuracy. As described in Fig. 4b, the relative deviations of the predicted qualified fibers values in comparison with experimental values indicated that the new model established by the SVM algorithm clearly possessed promising prediction properties.

Fig. 4. Predicted versus experimental qualified fibers using the SVM algorithm: (a) predicted qualified fibers, and (b) relative deviation

As described in Fig. 5, the percent of value shows 97.02% in the range of 0% to 5%, 2.39% in 5% to 10%, 0.43% in 10% to 15%, and 0.17% in 15% to 20%. The associated estimation error for the majority of data points (99.41%) lies in the range of 0% to 10%, which demonstrates a good predictive accuracy for the fiber quality.

Fig. 5. Percent of value for relative deviation range under the SVM model

The Comparison of Results for SVM and MLR

The performance of SVM was compared with that of MLR based on the training set and test set. Figure 6 shows the predicted qualified fibers versus the experimental ones using the MLR algorithm. In Fig. 6a, the predicted qualified fibers were not in good agreement with the corresponding experimental fiber quality, and fewer data points of qualified fibers were close to the diagonal. As shown in Fig. 6b, approximately 11.2% of the relative deviation of predicted qualified fibers was above 20%, and the maximum value even reached up to 51.2%. This demonstrated that the fiber quality prediction was not a simple linear problem, and the fiber quality model established by the MLR linear algorithm had certain limitations.

Fig. 6. Predicted versus experimental qualified fibers using the MLR algorithm: (a) predicted qualified fibers, and (b) relative deviation

The errors of the MLR algorithm were compared with that of the SVM algorithm, and the detailed results are listed in Table 1. As can be seen in Table 1, on the basis of MAE, the SVM achieved a decrease of 92.19% compared to the MLR. Similar results were obtained from the other error analyses, such as a MRE with a decrease of 92.36%, RMSE with a decrease of 87.29%, and TIC with a decrease of 87.21%. It was demonstrated that the SVM algorithm performed better than the MLR in predicting fiber quality.

Table 1. Comparison of Errors by the MLR and SVM Algorithms

Application of the Predictive Model of Fiber Quality

According to actual producing conditions, the ranges of production parameters selected in the prediction model are shown in Table 2. Because the established predictive model of fiber quality was effective and accurate, the effect of variation of single parameters on fiber quality was analyzed.

Table 2. The Range of Production Refining Parameters

One parameter was selected as a variable and the other four parameters remained unchanged in each prediction. A parameter was selected as a variable with its range provided in Table 2. However, the values of parameters selected as constants were SR = 50 r/min, CH = 5 m, OV = 45%, CB = 12.55%, and CP = 24%. By inputting those values into the developed fiber quality model, the variations of each production parameter with the fiber quality were obtained as shown in Fig. 7.

Fig. 7. Variations of the percentage of qualified fibers with each production parameter: (a) SR, (b) CH, (c) OV, (d) CB, and (e) CP

As shown in Fig. 7a, when the SR was less than 62 r/min, the percentage of qualified fibers increased with the increase of SR. However, when it was more than 62 r/min, the percentage of qualified fibers decreased with the increase of SR. The increase in SR led to the increase in the amount of chips and resulted in enhanced friction among fibers between refining disks. The increased friction greatly reduced the excess crushing and cutting of fibers due to its direct contact with the disks and improved the fiber quality. However, too much feed caused by the excessive conveyer screw revolution speed led to the increase of the gap between disks. This reduced the force between the wood chips, resulting in insufficient defibrination. Consequently, more coarse fibers were produced, resulting in a reduction in the amount of qualified fibers.

As shown in Fig. 7b, when the CH was less than 5 m, the percentage of qualified fibers increased with the increase of CH; however, when it was more than 5 m, the percentage of qualified fibers decreased with the increase of CH. The CH had a positive correlation with steaming time, which is a vital parameter deciding the degree of softening of wood chips. The increase in CH extended the steaming time of the chips and weakened the inter-fiber connection. As a result, the percentage of qualified fibers increased due to better fiber shape and less mechanical damage of fibers during refining. Nevertheless, excess steaming time destroyed the fiber itself due to the increase of the proportion of rupture of long molecular chains of the cellulose and hemicellulose in fiber cells, and the intensification of the loss of intercellular layer material, which was likely to generate more tiny fibers and resulted in a decreased percentage of qualified fibers.

As shown in Fig. 7c, the percentage of qualified fibers generally decreased with the increase of OV. The increase in OV caused a gradual increase of the pressure difference between the inside and outside of grinding chamber, where the chips defibered, and the aggravating of fluidity outward of fibers. This led to the decrease of the percentage of qualified fibers due to the increase in coarse fibers caused by the insufficient grinding, which resulted from less refining time. However, there existed an optimum point at 18.2%, where the fluidity outward and degree of grind of fibers achieved a balance, leading to the optimum fiber size.

As shown in Fig. 7d, when the CB was less than approximately 12%, the percentage of qualified fibers increased with the increase of CB. However, when the CB was more than approximately 12%, it decreased with the increase of CB. The fiber quality was improved by the longer bast fibers in the bark. However, if the bark content was too high, excess skinny fiber in the bark resulted in the reduction of the percentage of qualified fibers.

As shown in Fig. 7e, when the CP was less than approximately 24%, the percentage of qualified fibers increased with the increase of CP. However, when the CP was more than approximately 24%, it decreased with the increase of CP. Chinese poplar is a broad-leaved wood whose fibers are shorter. Therefore, the proportion of coarse fibers decreased and the proportion of finer fibers slightly increased when the CP increased, improving the ratio of qualified fibers. Nevertheless, too much CP resulted in the exceeding proportion of finer fibers, which reduced the percentage of qualified fibers.

CONCLUSIONS

1. Using 1173 groups of mill production data, an effective and powerful model based on the SVM algorithm was developed to predict the fiber quality during the refining process in the MDF production. The model was constructed with five predictor variables including the conveyer screw revolution speed (SR), accumulated chip height (CH), opening ratio of the discharge valve (OV), content of log bark (CB), and content of Chinese poplar (CP), which enable the evaluation of fiber quality more comprehensive and more suitable to be used in the mill conditions.

2. The predicted outputs of the SVM model were compared with the experimental data. The results showed that the predicted values agreed with the experimental ones well, and the estimation relative deviation for the majority of data points (99.41%) were in the range of 0% to 10%, indicating that the new model established by the SVM algorithm had a good predictive accuracy.

3. The performance of the SVM was compared with that of MLR based on the training set and test set. The results demonstrated that the MAE, MRE, RMSE, and TIC values of the total data set of the SVM were decreased by 92.19%, 92.36%, 87.29%, and 87.21%, respectively. Compared to the model established by MLR, the model established by SVM possessed better properties in accuracy and validity.

4. Additionally, the variations of the percentage of qualified fiber with each parameter were presented using the developed predictive model, and the reasons of variations were described. It was demonstrated that the model established by the SVM algorithm can be used as a promising guideline for a fiberboard mill manager to predict fiber quality.

ACKNOWLEDGMENTS

This research was funded by the Graduate Student Independent Innovation Fund Project of Central Universities (Grand No. 2572017AB17) and the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20130062110005).

REFERENCES CITED

Bishop, C. M. (1995). Neural Networks for Pattern Recognition, Oxford University Press, New York, NY.

Chen, G. (2012). Model Analysis and Experimental Study on Mechanism of Thermomechanical Refining Plate Separating Fiber, Ph.D. Dissertation, Northeast Forestry University, Harbin, China.

Chen, G., and Hua, J. (2009). “Effect mechanism of fiber quantity and quality on refiner energy consumption,” Journal of Northeast Forestry University 37(6), 37-40.

Chu, F., Dai, B. W., Dai, W., Jia, R. D., Ma, X. P., and Wang, F. L. (2017). “Rapid modeling method for performance prediction of centrifugal compressor based on model migration and SVM,” IEEE Access 5(99), 21488-21496. DOI: 10.1109/Access.2017.2753378

Drucker, H., Burges, C. J. C., Kaufman, L., Smola, A., and Vapnik, V. (1997). “Support vector regression machines,” Advances in Neural Information Processing Systems 28(7), 779-784.

Esteban, L. G., Fernandez, F. G., De Palacios, P., and Rodrigo, B. G. (2010). “Use of artificial neural networks as a predictive method to determine moisture resistance of particle and fiber boards under cyclic testing conditions (Une-En 321),” Wood and Fiber Science 42(3), 335-345.

Gerstorfer, E., Runkler, T. A., Schlang, M., Junnemann, E., and Hollatz, J. (2001). “Intelligent control of the refining process for improved fiber quality,” in: 2001 European Control Conference, Porto, Portugal, pp. 751-756.

Hong, W.-C., Dong, Y., Zhang, W. Y., Chen, L.-Y., and Panigrahi, B. K. (2013). “Cyclic electric load forecasting by seasonal SVR with chaotic genetic algorithm,” International Journal of Electrical Power & Energy Systems 44(1), 604-614. DOI: 10.1016/j.ijepes.2012.08.010

Hua, J., Chen, G., Xu, D., and Shi, S. Q. (2012). “Impact of thermomechanical refining conditions on fiber quality and energy consumption by mill trial,” BioResources 7(2), 1919-1930. DOI: 10.15376/biores.7.2.1919-1930

Huang, H.-X., and Lu, S. (2016). “Neural modeling of parison extrusion in extrusion blow molding,” Journal of Reinforced Plastics and Composites 24(10), 1025-1034. DOI: 10.1177/0731684405048201

Huang, S., Cai, N., Pacheco, P. P., Narrandes, S., Wang, Y., and Xu, W. (2018). “Applications of support vector machine (SVM) learning in cancer genomics,” Cancer Genomics Proteomics 15(1), 41-51. DOI: 10.21873/cgp.20063

Jia, N., Liu., B., Hua, J., and Lin, X. (2015). “Effects of bark proportion on defibrator energy consumption and fiber quality,” China Wood Industry 29(3), 35-38.

Jiao, G., Guo, T., and Ding, Y. (2016). “A new hybrid forecasting approach applied to hydrological data: A case study on precipitation in northwestern China,” Water 8(9), 367-381. DOI: 10.3390/w8090367

Keerthi, S. S., and Lin, C. J. (2003). “Asymptotic behaviors of support vector machines with Gaussian kernel,” Neural Computation 15(7), 1667-1689. DOI: 10.1162/089976603321891855

Ma, J., Theiler, J., and Perkins, S. (2003). “Accurate on-line support vector regression,” Neural Computation 15(11), 2683-2703. DOI: 10.1162/089976603322385117

Mokhtarzad, M., Eskandari, F., Vanjani, N. J., and Arabasadi, A. (2017). “Drought forecasting by ANN, ANFIS, and SVM and comparison of the models,” Environmental Earth Sciences76(21), 729-738. DOI: 10.1007/s12665-017-7064-0

Roushangar, K., and Ghasempour, R. (2017). “Prediction of non-cohesive sediment transport in circular channels in deposition and limit of deposition states using SVM,” Water Science and Technology Water Supply 17(2), 537-551. DOI: 10.2166/ws.2016.153

Runkler, T. A., Gerstorfer, E., Schlang, M., Jünnemann, E., and Hollatz, J. (2003). “Modelling and optimisation of a refining process for fibre board production,” Control Engineering Practice 11(11), 1229-1241. DOI: 10.1016/S0967-0661(02)00233-2

Shi, F., Wang, X. C., Yu, L., and Li, Y. (2010). 30 Case Analysis of MATLAB Neural Network, Beihang University Press, Beijing, China.

Shi, J. L., Zhang, S. Y., and Riedl, B. (2006). “Multivariate modeling of MDF panel properties in relation to wood fiber characteristics,” Holzforschung 60(3), 285-293. DOI: 10.1515/HF.2006.046

Sun, Z., Wang, C., Niu, X., and Song, Y. (2016). “A response surface approach for reliability analysis of 2.5D C/SiC composites turbine blade,” Composites Part B: Engineering 85, 277-285. DOI: 10.1016/j.compositesb.2015.09.025

Vapnik, V. N. (1999). “An overview of statistical learning theory,” IEEE Transactions on Neural Networks 10(5), 988-999. DOI: 10.1109/72.788640

Wang, D. C., Wang, M. H., and Qiao, X. J. (2009). “Support vector machines regression and modeling of greenhouse environment,” Computers and Electronics in Agriculture 66(1), 46-52. DOI: 10.1016/j.compag.2008.12.004.

Xiao, C., Hao, K., and Ding, Y. (2014). “The bi-directional prediction of carbon fiber production using a combination of improved particle swarm optimization and support vector machine,” Materials (Basel) 8(1), 117-136. DOI: 10.3390/ma8010117

Xing, C., Deng, J., and Zhang, S. Y. (2006). “Effect of thermo-mechanical refining on properties of MDF made from black spruce bark,” Wood Science and Technology 41(4), 329-338. DOI: 10.1007/s00226-006-0108-3

Yan, K., and Shi, C. (2010). “Prediction of elastic modulus of normal and high strength concrete by support vector machine,” Construction and Building Materials 24(8), 1479-1485. DOI: 10.1016/j.conbuildmat.2010.01.006

Zhang, X. Y., Li, J. G., Zhang, Y. Z., Cai, W., and Liu, H. B. (2016). “Data-driven method for pulp properties estimation in stock preparation,” BioResources 11(2), 4947-4963. DOI: 10.15376/biores.11.2.4947-4963

Zhao, Y., Zhang, X., Deng, L., and Zhang, S. (2016). “Prediction of viscosity of imidazolium-based ionic liquids using MLR and SVM algorithms,” Computers & Chemical Engineering 92(9), 37-42. DOI: 10.1016/j.compchemeng.2016.04.035

Zhao, Y., Zhao, J., Huang, Y., Zhou, Q., Zhang, X., and Zhang, S. (2014). “Toxicity of ionic liquids: Database and prediction via quantitative structure-activity relationship method,”Journal of Hazardous Materials 278 (8), 320-329. DOI: 10.1016/j.jhazmat.2014.06.018

Zhou, Z., Yin, J. X., Zhou, S. Y., Zhou, H. K., and Zhang, Y. (2016). “Detection of knot defects on coniferous wood surface using near infrared spectroscopy and chemometrics, ” BioResources 11(4), 9533-9546. DOI: 10.15376/biores.11.4.9533-9546

Article submitted: July 1, 2018; Peer review completed: August 28, 2018; Revised version received and accepted: September 8, 2018; Published: September 13, 2018.

DOI: 10.15376/biores.13.4.8184-8197