VOLUME AND TAPER EQUATIONS FOR COMMERCIAL STEMS OF Nothofagus obliqua AND N. alpina ABSTRACT: Timber volume of standing trees is essential information for management

Timber volume of standing trees is essential information for management decisions. The increasing need to optimize the potential capacity of forests maintaining their conservation, requires the quantification of the different potential possible timber products. The aim was to adjust taper equations to determine volumes of different timber products for commercial stems of Nothofagus alpina and N. obliqua. Trees of both species were randomly selected in harvesting areas of Lanin National Park (Argentina). Trees were felled and cut into commercial logs, measuring diameter with bark at different heights up to the beginning of the crown, and for each tree the diameter at breast height and total height. Five taper equations were selected and non-linear regression processes were employed for the fittings. We obtained the volume through the integration of the stem profile equation and the rotation in the space thereof through solid of revolution. The Bennet and Swindel (1972) model was selected for both Nothofagus species, obtaining similar equation parameters and differences were observed at the top of the stems of larger trees. For this the use of an integrated model is not recommended. With the obtained equations it is possible to: (i) estimate volume at different heights and for different commercial diameters, and (ii) predict the height at which both species reach to a certain diameter. The model presented some statistical limitations (e.g. multicollinearity), however, the fitting of the equation and the easy understanding of the outputs support it as a useful tool in a broad range of forest applications. EQUAÇÕES DE AFILAMENTO E VOLUME PARA TAMANHOS COMERCIAIS DE Nothofagus obliqua E N. alpina RESUMO: O volume de madeira das árvores em pé é uma informação essencial para as decisões de gestão. A crescente necessidade de otimizar a capacidade potencial das florestas mantendo sua conservação, requer a quantificação dos diferentes possíveis produtos de madeira possíveis. O objetivo era ajustar equações de afilamento para determinar volumes de diferentes produtos de madeira para hastes comerciais de Nothofagus alpina e N. obliqua. As árvores de ambas as espécies foram selecionadas ao acaso nas áreas de colheita do Parque Nacional Lanin (Argentina). As árvores foram derrubadas e cortadas em troncos comerciais, medindo o diâmetro com a casca em diferentes alturas até o início da coroa, e para cada árvore o diâmetro na altura do peito e altura total. Foram selecionadas cinco equações cônicas e foram empregados processos de regressão não-linear para ajustar. Obtivemos o volume através da integração da equação do perfil do caule e da rotação no seu espaço através do sólido da revolução. O modelo de Bennet e Swindel (1972) foi selecionado para as duas espécies de Nothofagus, obtendo parâmetros de equação semelhantes e observando-se diferenças no topo das hastes de árvores maiores. Portanto, o uso de um modelo integrado não é recomendado. Com as equações obtidas é possível: (i) estimar o volume em diferentes alturas e para diferentes diâmetros comerciais, e (ii) prever a altura em que ambas as espécies atingem um determinado diâmetro. O modelo apresentou algumas limitações estatísticas (por exemplo, multicolinearidade), no entanto, o ajuste da equação e a compreensão fácil das saídas o suportam como uma ferramenta útil em uma ampla gama de aplicações florestais.


INTRODUCTION
Lanin National Park (Argentina) has a long tradition in forest management of mixed forests of Nothofagus alpina Poepp. et Endl. (sensu N. nervosa Phill (Dim. et Mil.)) commonly named raulí and N. obliqua (Mirb.)Oerst commonly named roble pellín, and recently, in-situ and ex-situ plantation tests begun.These species are regionally considered as the best potential alternative species for commercial plantations.To date, exotic species (e.g.Pinus contorta Douglas ex Lauden), pino ponderosa (P.ponderosa Douglas ex C. Lawson) y pino oregón (Pseudotsuga menziesii Mirb.Franco) are employed in plantations, however, native Nothofagus species are prefered due to the small commercial benefits (e.g.growth and timber quality) and the negative impacts over the native flora and fauna of the exotic trees.
The approach to forest management has shifted towards a multifunctional vision of the forest considering: (i) integration of biodiversity conservation in forest management process, (ii) characterization and maintenance of the multifunctionality of forests, (iii) the effective participation of stakeholders, and (iv) the use of new technologies (GÓMEZ GARCÍA et al. 2016).In these sense, the development of new tools for forest management planning allowed to improve the current silvicultural proposals.
The hardest task to achieve during forest inventories is to estimate the individual volume of trees and classified it into the different timber products for local markets (MARTÍNEZ PASTUR et al. 2000).The increasing needs to optimize the timber potential of the forests, requires the quantification of the different possible timber products (ANDRADE, 2014).Therefore, it is critical to use biometric tools to promote appropriate silvicultural planning, and encourage sustainable use of forests.
Usually, forest yield is calculated using volume values provided by the total forest standing trees.However, volume is very difficult to measure directly (even more so when trees are standing) for this, volume is generally estimated from one or more variables that are positively correlated, such as height and diameter at breast height (GREGOIRE et al. 1986, WIANT et al. 1992).Stand volume can be obtained from individual tree models, resulting from the addition of all individual trees.However other models provide predictions at stand level directly (HARI, 1996).
Total tree volume can be estimated from volume equations (CHAUCHARD et al. 2005) which are useful and accessible; nevertheless, these functions only allow to estimate the total values for the tree without classifying the data in potential timber products that can be obtained along the stem (PRODAN et al. 1997, MENÉNDEZ-MIGUÉLEZ et al. 2014).One of the purposes of the managers is to perform a more efficient use of the forests, considering different alternatives for the timber products obtaining through different harvesting proposals (JUNIOR et al, 2012, CHAVES et al. 2014).In this sense, the information about different wood products is an advantage during forest planning resulting in a more efficient use during the production chain, as well as maximizing benefits and reducing the impact caused by market fluctuations of certain products (DOS REIS et al. 2016).To overcome this issue, stem profile equations have been proposed, which allows to estimate the diameter of a tree at any height and vice-versa (CLUTTER et al. 1983), and from mathematical conversions of these equations, total volume of each product can be estimated.
A taper equation is defined as the diameter reduction rate of the tree stem from the base toward the apex, which determines a form of thereof (CRECHI et al. 2009, SOUZA et al. 2016).The outputs generate more precise volume estimations, and are more efficient their use to build growth and yield models and sawmill simulators (GEZAN et al. 2009).There are several types of equations (KOZAK, 2004;SÁNCHEZ GONZÁLEZ et al. 2007;POMPA GARCÍA;SOLÍS MORENO, 2008;GEZAN et al. 2009;ANDRADE, 2014), where the most commonly employed independent variables were: diameter at breast height (DBH), tree height, crown ratio, and pairs of diameter/height along the stem (HANN et al. 1987).These variables are easy and are standard measurements during forest inventories.To construct taper equations, it was required data from individual trees, which can derived from models of longitudinal structure or multiple measurements (diameter and height data pairs) of each individual.
To date, stem profile models for Nothofagus are developed both in Argentina and Chile.In Argentina, a taper model for N. pumilio was constructed, where site index was added as an independent variable (CELLINI et al. 2002, MARTÍNEZ PASTUR et al. 1997, CELLINI et al. 2012).In Chile, spline functions were developed for N. pumilio (TRINCADO; VIDAL, 1999) and taper models for saplings of N. obliqua, N. alpina and N. dombeyi (GEZAN et al. 2009), but it has been recommended to use in DBH < 50 cm for N. alpina and DBH < 60 cm for N. obliqua.Regional changes in taper equations for N. obliqua were also studied.These differences in growth patterns (GEZAN et al. 2009, ATTIS BELTRÁN et al. 2015) supporting the idea to develop specific stem profile functions for Nothofagus alpina and N. obliqua in Lanin National Park.For this, the aim was to adjust taper equations to determine volumes of different timber products for commercial stems of N. alpina and N. obliqua.

Study area
The distribution range of Nothofagus obliqua is located between 33º and 41.5º SL, while N. alpina covers a slightly narrower latitudinal range between 35º and 40.5ºSL (DONOSO et al. 1993;SALAS;GARCÍA, 2006).The sampling was performed within the basin of the lakes Lácar and Nonthué (40° 07 ' SL -71° 27› WL to 40° 07' SL -71° 39› WL) covering 1,500 hectares where forest management units of the Reserve of Lanín National Park are located in the province of Neuquén (Argentina).The forests are dominated by pure and mixed Nothofagus forest types located between 600 and 1,100 m a.s.l.(CHAUCHARD;SBRANCIA, 2003).Climate is temperate humid with average annual temperatures of 9°C and average rainfall of 1,800 mm .yr -1 .The region was shaped by glaciers and later covered by volcanic ash (MERCER, 1976), and from this process emerged the dominant soils developed with the Holocene volcanic ash, and belonging to the order Andisols.They present a sequence of horizons O-A-Bw-C or sequences of these horizons because of successive ash deposits which resulted buried soils (FERRER et al. 1991, FRUGONI et al. 2005).

Data
Trees of N. alpina (n=84) and N. obliqua (n=98) were randomly selected in the harvesting areas and close unharvested areas in order to cover all the site qualities gradient.The sampling was designed to cover all the diameter and height distribution, including only trees with no visible damage or widespread diseases.Trees were cut and the following variables were measured in each tree: DBH in cm and total height in m (H).After that, commercial stem (from the stump cut to the bottom part of the crown) was divided into timber logs (2.2 to 5.0 m length) according to the market demands.In each cut, diameter and height data pairs (d-h) were taken.To detect possible anomalies in the data scatter plot between relative diameters and height (d-h data pairs related to H) were examined.
Their expressions correspond to linear and nonlinear equations that can be grouped into three different types according to their mathematical construction, where d: diameter (m) at each h (m); DBH: diameter at 1.3 m; h: height (m) at each d; H: total height (m).b 1 to b 7 : function parameters.

Fitting and self-evaluation
Nonlinear regression process was applied for fi ttings, where parameters were obtained by the Marquardt method (1963), which it is an iterative method for solving nonlinear problems in parameter estimation by least squares.To evaluate the prediction of selected equations, the adjustment and the analysis of residual were analyzed.Statistical evaluation was performed by the coeffi cient of determination (R 2 ), the adjusted coeffi cient of determination by the number of parameters (R 2 aj ), the standard error of estimate (SEE), and the mean absolute error (MAE).Root mean square error (RMSE) and bias were used to analyze the residuals.The percentage values of each goodness-of-fi t statistic were obtained dividing by the average of the observed values and multiplying by 100 (RMSE% and bias %) ([6] and [7]).Errors for each equation was the difference between observed and predicted value.Further, statistical errors were graphed through the predicted dependent variable (KOZAK, 2004).The residual serial correlation was examined graphically and using the Durbin-Watson statistic (Durbin;Watson, 1951), where, n is the number of observations and, y o and y e are the experimental and predicted response values, respectively. [6] [7]

Estimated volume from the stem profi le function
Through the integration of the stem profi le function and the rotation in the space thereof, can be obtained through the revolution volume estimates [8] (CHAVES et al. 2014).Through this function, volume can be estimated between two determined heights (h i and h j ).
To assess the accuracy of the model to predict volume, fi rst the volume of each tree logs cut up to the height of the beginning of crown was calculated using the Smalian formula (AVERY; BURKHART, 1994).The volume between the base and the beginning of crown was estimated by the sum of the individual volumes of each section (observed volume), and the volume of each sample was estimated with the function [8].
[8] Subsequently errors were calculated as the difference between the observed and predicted data by function [8], which were used for analysis of residual data.The following statistical data were calculated: mean absolute error (MAE), root mean square error (RMSE and RMSE %) and bias.Besides this, graphs were made with: (i) errors estimates of stem volume according to DBH classes, and (ii) errors estimates of log volume depending on the average diameters.

RESULTS
Sample was composed by 84 trees of Nothofagus alpina and 98 trees of N. obliqua used for the fi tting of the equations (Table 1).Minimum and maximum diameters were upper and under the diameter of timber trees.The data distribution along diameter and height gradients for both species cover the natural range in the studied area (Table 2).Finally, the distribution of the data pairs of d-h of the sample trees not presented anomalies in the scatter plots (Figure 1).The statistics obtained for the proposed equations showed good fi ttings and low bias values (Table 3).However, among the fi ve tested equations, [1], [3] and [4] were those which obtained the least bias, in addition to presenting signifi cance in all of its parameters (Table 4).
Bias gradients were analyzed along the stem (Figure 2) for both species through the root mean square percentage error (RMSE %) and average bias    Solving the derivative, the volume function was defi ned as [9]: The values of the parameters for the volume equation (b1 to b4) are shown for both species in Table 4.With these functions [10] which estimate the stem volume up to the beginning of crown for N. alpina and N. obliqua [10], we obtained an acceptable statistical precision and bias, not exceeding 10% of RMSE, and where MAE was lower than 0.15 m 3 per individual in average (Table 5).If we compared the values of RMSE% at different DBH classes, they not exceed 16%, while bias% remained within the range of -7% to 2% (Figure 5).With this tool it is possible to know the volume up to each diameter according to a specifi c commercial purpose and calculate it, using the inverse function determined for the stem profi le, cleared in  The RMSE% of N. obliqua were lower than 13.6% and below 14.5% for N. alpina for the three equations, while the bias remained between -8.8 and 5.5% for N. obliqua and between -8.9 and 7.4% for N. alpina.Equation [1] had the lowest RMSE% values to relative heights between the beginning and the middle of the stem.The bias maintained between -3.8% and 5.0% along the stem.The RMSE% of N. alpina was lower than 14.0% for the three equations and the bias remained between -8.9% and 7.4%.Equation [1] presented, as for N. obliqua, the lowest RMSE% values to relative heights between the beginning and the middle of the stem, and the bias was maintained between -8.4% and 5.0%.Equation [3] and [4] presented the same trend.
Based on the bias analyses and the general accuracy, the equation [1] was selected to estimate the stem profi le of N. alpina and N. obliqua.Residuals plots (Figure [10] [11] [12] [13] [14] the function [1] height (h) to which a tree reaches a certain diameter (d) [15].Thus the log volumes can be estimated according to different trade lengths and commercial diameters using function [10], where Z is Beside this, during the estimation of log volumes, RMSE% values were close to 45% in logs of 0.20 cm diameter, but from larger diameters RMSE% remained below 10-20% in the basal portion of the stem for both species.Bias% estimated for the log volumes had similar behavior, being greatest in small diameters and in larger diameters for N. alpina (-1.7% and 5.0%) and for N. obliqua (0.3% and 3.0%) (Figure 6).

DISCUSSION
The data obtained for both species fulfi ll the ranges of the studied variables (Table 1).Adequate distribution of values for height and DBH classes were found by using equivalent data employed in other research for other Nothofagus species (CELLINI et al. 2002;GEZÁN et al. 2009).Similar data thresholds were found for both species at ranges of total height and DBH, as well as the data distribution classifi ed by height and DBH classes (Table 2).In addition it was also observed similarity pattern distributions between scatter plots of relative heights in relation to relative diameters for the studied species (Figure 1).
The selection of the best equation for each species was based on multi-criteria, as was stated before.This was convenient, due to all models had higher coeffi cient of determinations (Table 3), however along the nonlinear adjustment the selection decisions should not only be made based on R 2 values, although this measurement gives an intuitive idea of how the variability explains a  Bennett and Swindel (1972).
particular fi tting.Thus, we given a greater importance in the decision to the statistical error, the signifi cance of the equation parameters and the residual dispersion analyses.Therefore, three adjusted equations were chosen out of the fi ve initial candidates ([1], [3] and [4]) due to their lowest error values, greater accuracy, and especially because N. obliqua and N. alpina presented good signifi cance for all the studied parameters.Particularly, model [5]can be successfully adapted to local conditions and different tree species by changing one or more of the three components, to produce models with better predictive ability (KOZAK, 2004).
The adequate prediction of the stem basal section is very important because it is where the greatest amount of wood quality is located.RMSE% and bias% behavior depending on the relative height allowed to evaluate this premise and it was found that the behavior was similar between the two species.In this context, equation [1] was superior for N. alpina as well as N. obliqua among the 3 selected equations.This equation meets the condition that at the height of DBH the estimated diameter is equal to the DBH and at total height estimated diameter equals to zero.Due to we not used constrained models, an error were detected: parameter "b1" is not equal to 1. Anyway parameter "b1" for both species is 1.025 and 1.021, less than 1 cm in a tree of 100 cm (close to 1% error in average).
In equation [1] for N. alpina RMSE% values increased in the fi rst quarter of the relative height without reaching unacceptable values.Probably the loss of precision in the basal portion of the stem was due to N. alpina database had two individuals with DBH greater than 1.20 m, and both trees presenting the same outlayers magnifying the error values in this portion for these large diameters.It is common for old trees with large diameters, that they tend to have a markedly neiloide form in the basal portion of the stem (ROJO et al. 2005), but these trees are less frequent in the timber forests.However, to use the equations for bigger trees (up to 1 m DBH) and if subsequently it is necessary determine the volume, an alternative is to consider the shape of a cylinder of diameter equal to DBH between the beginning of the stem and the DBH.Using this constrain, the estimation of the basal part of the stem will avoid the basal overestimation in the large old trees.The small differences in the coefficients of the adjusted equations are reflected in stem profiles of Figure 4, which presented slightly different shapes, however, larger trees have different stem profiles as was discussed before.Nothofagus alpina showed a constant diameter decrease along the stem profile, while N. obliqua decrease is moderate at the beginning and then decrease more quickly at the end (Figure 4).However, among the studied trees, these differences may occurred within the crown, where the equations are not completely valid.For example, the beginning of the crown height was in average at 12 m, and the differences in the stem profile are more evident only in the smaller trees (Figure 4).In future studies it must be considered adding the crown height variable for the equation adjustment as well as for the model selection.
There are many reasons for the discussed observed behaviours from those expected in the equations, many of them can be explained due to differences between expected and observed patterns.In this sense, native forests presented greater heterogeneity than plantations, due to natural irregularities along the stem cross section.Therefore wild forest rarely provides theoretical perfection for the equations, as was observed in the plantation trees.Irregularities in the shape are usually seen where tress has lost their branches by natural lopping, or where they still persist.Other reasons for irregularities in the shape of the stem can be due to fire damages, natural diseases or insect attacks.Most of the trees in the nature are eccentric in shape, so they are wider in one direction than in another (WEST, 2009).Also, this can be a response to a wind exposure or a result from the heterogeneous spatial arrangement among the trees (ROBERTSON, 1991) and their intraspecific competition.
As limitations and weaknesses we can mention that equation [1] contain also several terms with the same regressed variable resulting in multicollinearity problems.In addition, autocorrelation is not very high, but maybe can be improved using another methodology (e.g.Rojo et al. 2005).Since the database contains multiple observations for each tree (i.e., hierarchical data), it is reasonable to expect autocorrelation within the residuals for each individual.The hierarchical structure may lead to a correlation between residuals of observations below DBH in the same way as residuals of observations above DBH.Multicollinearity and autocorrelation are some of the most important problems that come when models are applied to forest data (ROJO et al. 2005) Volume function [10] presented good outputs when stem volume and logs of different diameters were tested.The estimated value using the Smalian equation is not the real volume and includes an unknown error, but when the ratio among the two diameters (which is the most common in our samples) is higher than 0.82, the error is less than 1% (WEST 2009).When we estimate stem volume, the values of the RMSE% tends to decrease in trees of greater DBH, which is desirable from a practical point of view and also from management, due to less mistakes were produced in the estimation of large trees.Usually, these trees (40-60 cm DBH) are of highest value (CELLINI et al. 2002), and the bias% tends of the equations were similar and presented lower values across the different lower DBH classes (Figure 5).
When volumes of different logs were estimated, RMSE% and bias% presented similarly outputs for full stems, but decrease towards the larger diameters (Figure 6).When the estimation considered each individual log, the error outputs were higher for diameters between 0.15 m and 0.20 m when were presented as percentage, but when the value was presented in cubic meters they are not significant.Nonetheless, the use of the model [10] for N. alpina and N. obliqua is recommended for trees with DBH between 0.25 m and 1.00 m, and medium log diameters between the same values.
Any natural dynamic process or silvicultural practice affecting the future shape of the trees, and has effect on tree growth modifying the stem profile (Souza et al. 2016).However, the shape equations showed enough flexibility to estimate the stem profiles of the sampled trees, which were not exempt to the influence of natural or anthropogenic phenomena.It is important to recognize that the stem diameter estimation at different heights, as well as the volume estimation, presented some bias when the models are applied at the beginning of the stem or at the crown.But the timber material is outside crowns.The other wood material is commonly used for lumber.
Thus it has been possible to simplify this inconvenience if we adjust stem profile functions especially for broadleaved, and not use those equations created for conifers.This results in a disadvantage, because estimates cannot be made within the crown height, although the volume can be obtained by difference with conventional models (CHAUCHARD et al., 2005), as it is the volume conformed by the products of lesser timber value.
To summarize the objective of this work, it was possible to obtain the stem profile equations and to estimate the volume for different timber products, including the most profitable timber products of the harvested trees.Some literature suggests that to use equations designed for conifers can be applied in broadleaf species, and they propose a transformation consisting in subtracting to each measured diameter the measure at the beginning of crown, however, this proposal do not solve the mentioned disadvantage (TRINCADO, 1996;TRINCADO;SANDOVAL, 2002).
Adjusted equations are useful management tools, which also allow knowing the growth and development of the stems for both species according to their dendrometric magnitudes.The biometric knowledge of the species under this study can be used for planning and the design of appropriate silvicultural proposals that allows to achieve sustained use of Nothofagus forests.Usually these trees were harvested for timber wood for sawmills and fire woods, however, one option is to increase the uses, allowing a significant increase in the incomes from a forest and improved forest management (CAMPOS et al. 2014).The equations developed attend to maintain a structural richness of the forest, allowing to estimate forest production and based on estimations, to determine extractability, so that the extraction does not exceed the growth.In addition, the classification of the product becomes an important tool to evaluate the role of forest and wood products, especially in light of in a climate change context, due to fixation and carbon sequestration (RODRÍGUEZ et al. 2015).

CONCLUDING REMARKS
Stem profile functions have been developed for Nothofagus alpina and N. obliqua in order to be applied in the basin of lakes Lacar and Nonthué.The equation chosen for each species was Bennett and Swindel (1972), and presented similar parameters between the two species, but with differences at the top of the stem in larger trees.For this, their use in these cases were not recommended.Volume functions were also obtained derived from the stem profile equations, being possible to estimate log volume at different heights and for different commercial diameters.Finally, a model to predict the height at which N. alpina or N. obliqua trees reaches at certain diameter was obtained.With these functions it is possible to determine volumetric existences in standing trees of N. obliqua and N. alpina growing in the Lanín National Park, which can be used for a better forest management of these forests.The model presented some statistical limitations (e.g.multicollinearity), however, the fitting of the equation and the easy understanding of the outputs support it as a useful tool in a broad range of forest applications.

FIGURE 1
FIGURE 1 Dispersion of relative height according to relativediameter, of the trees used to fi t the stem profi le model.

3
) and Durbin-Watson statistic indicated that residual autocorrelation of the equation [1] was low for N. alpina (DW= 1.28) and for N. obliqua (DW= 1.02).Figure4showed the stem profi le of two trees of N. alpina and two trees of N. obliqua, for heights of 35 m and DBH of 45 cm, and heights of 20 m and DBH of 30 cm.With the obtained equations of stem profi le [1] it is possible to obtainfi nd the volume function, producing therewith a solid of revolution of each tree stem, where, V ij : volume between the height ijd: diameter (m) at height h, DBH: diameter at 1.3 m, h: height (m) at diameter d, H: total height (m).b 1 to b 7 : equation parameters.

FIGURE 3
FIGURE 3Bennett and Swindel (1972) [1] residuals for different stem height.(%)for the three selected equations ([1],[3] and [4]).The RMSE% of N. obliqua were lower than 13.6% and below 14.5% for N. alpina for the three equations, while the bias remained between -8.8 and 5.5% for N. obliqua and between -8.9 and 7.4% for N. alpina.Equation [1] had the lowest RMSE% values to relative heights between the beginning and the middle of the stem.The bias maintained between -3.8% and 5.0% along the stem.The RMSE% of N. alpina was lower than 14.0% for the three equations and the bias remained between -8.9% and 7.4%.Equation [1] presented, as for N. obliqua, the lowest RMSE% values to relative heights between the beginning and the middle of the stem, and the bias was maintained between -8.4% and 5.0%.Equation[3] and [4]  presented the same trend.Based on the bias analyses and the general accuracy, the equation [1] was selected to estimate the stem profi le of N. alpina and N. obliqua.Residuals plots (Figure

TABLE 1
Descriptive statistics of the trees that constitute the sample for the study.Standard deviation was presented between brackets.

TABLE 2
Distribution of sample trees by DBH and height classes, where (a) N. obliqua and N. alpina

TABLE 3
Statistics of the stem profi le equation adjustment.

TABLE 5
Statistics for error analyses in the volume function[9].RMSE: root mean square mean error; RMSE%: root mean square percentage error; MAE: mean absolute error.