Next Article in Journal
Smallholder Forestry in the Western Amazon: Outcomes from Forest Reforms and Emerging Policy Perspectives
Previous Article in Journal
Variation in the Growth Traits and Wood Properties of Chinese Fir from Six Provinces of Southern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developing an Improved Parameter Estimation Method for the Segmented Taper Equation through Combination of Constrained Two-Dimensional Optimum Seeking and Least Square Regression

1
Research Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China
2
Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China
3
Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Prague 6, Czech Republic
4
Penn State Hershey Cancer Institute, Penn State College of Medicine, Hershey, PA 17033, USA
5
College of Computer and Information Techniques, Xinyang Normal University, Xinyang 464000, Henan, China
6
Center for Statistical Genetics, Pennsylvania State University, Loc T3436, Mailcode CH69, 500 University Drive, Hershey, PA 17033, USA
*
Author to whom correspondence should be addressed.
Forests 2016, 7(9), 194; https://doi.org/10.3390/f7090194
Submission received: 21 July 2016 / Revised: 26 August 2016 / Accepted: 27 August 2016 / Published: 31 August 2016

Abstract

:
The segmented taper equation has great flexibility and is widely applied in exiting taper systems. The unconstrained least square regression (ULSR) was generally used to estimate parameters in previous applications of the segmented taper equations. The joint point parameters estimated with ULSR may fall outside the feasible region, which leads to the results of the segmented taper equation being uncertain and meaningless. In this study, a combined method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) was proposed as an improved method to estimate the parameters in the segmented taper equation. The CTOS & LSR was compared with ULSR for both individual tree-level equation and the population average-level equation using data from three tropical precious tree species (Castanopsis hystrix, Erythrophleum fordii, and Tectona grandis) in the southwest of China. The differences between CTOS & LSR and ULSR were found to be significant. The segmented taper equation estimated using CTOS & LSR resulted in not only increased prediction accuracy, but also guaranteed the parameter estimates in a more meaningful way. It is thus recommended that the combined method of constrained two-dimensional optimum seeking and least square regression should be a preferred choice for this application. The computation procedures required for this method is presented in the article.

1. Introduction

The mathematical function (or equation) describing the variation of the tree diameter at any point of the stem with the distance from the tree top is known as the stem taper equation [1]. Using this equation, one may calculate stem diameter at any arbitrary height and conversely, calculate tree height for any arbitrary stem diameter. Consequently, stem volume can be calculated for any log specification, and a volume equation can be developed for classified product dimensions [2]. In many countries, especially the United States and Canada, stem taper equations have been commonly used to replace volume tables and volume equations for estimating stem volume, merchantable volume ratio, 3-D reconstruction of stem and simulation, and optimization of bucking [3,4,5]. Therefore, developing a stem taper equation with high predictive precision has always been a matter of interest to forest managers and forest engineers as the taper equation provides technical support to forest management.
Numerous forms of the taper equation have been developed over the past century ranging from simple to complex. These taper equations can be divided into three categories: simple taper equation [1,6,7,8,9,10,11,12], segmented taper equation [13,14,15,16,17], and variable-exponent taper equation [3,18,19,20,21,22]. Simple taper equations are easy to fit, and can be used to calculate the stem volume applying the integration method. Relatively simple taper functions can effectively describe the general taper of the trees. However, they lack the ability to describe the entire stem accurately. They may provide reasonable estimates in a mid-portion of the stem, but usually are less accurate in estimating the profile in the butt or upper stem segments [3,13,18,23]. As compared to the simple taper equation, the segmented taper equation exhibits a higher precision of diameter estimation and also could be used to describe the entire stem profile more efficiently [13,24,25]. The variable-exponent taper equation also has a high estimation precision, but it cannot be used to estimate volume by integration without the aid of a computer program [3,17,26]. Overall, the segmented taper equation has greater flexibility than other taper equations and therefore has frequently been applied along with different modelling approaches (e.g., traditional least square regression, seemingly unrelated regression and nonlinear mixed-effects model) [14,17,24,25].
The Max and Burkhart [13] taper equation is one of the most typical segmented taper equations to describe tree stem taper [14,25,27]. This equation consists of three different quadratic functions joined at two joint points, and therefore is flexible enough to describe complex and irregular stem profiles. Many researchers have used this equation for the prediction of tree taper [15,17,24,25,28,29,30,31,32]. Two joint point (upper and lower) parameters (the feasible region of two parameters lies between 0 and 1) in the Max and Burkhart [13] taper equation were usually estimated by applying a traditional unconstrained least square regression (ULSR) method and subjective experience (SE) method [13,17,25,33]. However, in a practical application, it has been found that the estimates of the two joint point parameters obtained from ULSR sometimes fall outside the feasible region of the corresponding estimated parameters, which makes the results meaningless [34]. The study involved the application of a simplex method with constraint to estimate parameters of the Max and Burkhart [13] taper equation. However, this failed to converge and was found very slow computationally. In addition, the study also revealed that the estimates of the two joint point parameters obtained by ULSR varied with the starting values of the parameters. Therefore, the developed Max and Burkhart [13] taper equation based on ULSR is subject to a great uncertainty. The SE method has a great uncertainty in estimation of the join point parameters due to unpredictable man-induced factors and inevitable measurement errors [35]. These uncertainties directly influence the estimation precision and practicability of the taper equations.
The objective of this study was thus to propose a combined method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) as an improved method to estimate parameters of the segmented taper equation based on Hua’s optimum seeking theory and golden section search approach [36]. The data from three tropical precious tree species Castanopsis hystrix Miq. (C. hystrix), Erythrophleum fordii Oliv. (E. fordii), and Tectona grandis L.F. (T. grandis) in southwest China were used to evaluate the methods. Also, CTOS & LSR was compared with the commonly used ULSR approach by means of developing the Max and Burkhart [13] taper equation at both individual tree-levels and population average-levels.

2. Materials and Methods

2.1. Data

We used data from a total of 30 square-shaped permanent sample plots (PSPs) established in 2014 across three pure even-aged C. hystrix, E. fordii, and T. grandis plantations (each tree species with 10 PSPs) located in the Experimental Center of Tropical Forestry, Chinese Academy of Forestry Sciences (Figure 1). C. hystrix, E. fordii, and T. grandis are tropical precious tree species in southern China. The data were collected by the Research Institute of Forest Resources Information Techniques, Chinese Academy of Forestry. Stands were selected on the basis of the following criteria: (1) pure species; (2) fully stocked; (3) even-aged; and (4) minimum or no disturbance. Sample plots were avoided near major roads because road construction might have altered soil drainage and affected tree growth. Once a stand was located, a 600 m2 permanent sample plot was established. Four trees per sample plot, each representing dominant height, codominant height, average height, and suppressed height were selected. These trees were felled, and cross-sectional slices of stem (disks) were obtained from 0.33 m, 0.67 m, 1.3 m, and every 1.3 m interval thereafter until an outside bark diameter of 7.0 cm was reached. Disks were then taken every 20 cm between 7 cm and 4 cm outside bark diameter. Many tree characteristics such as outside bark diameter at breast height (D), total height (H), height to live crown base, and crown width, were measured for each sample tree. Summary statistics of D and H data are listed in Table 1. We applied two approaches for fitting and testing the Max and Burkhart [13] taper equations; (i) at individual tree-level using sub-dataset of nine dominant and codominant trees (three trees from each tree species) (Figure 2); (ii) at population average-level using full dataset of 120 trees (Table 1).

2.2. Taper Equation

The Max and Burkhart [13] taper equation is:
y = β 1 ( x 1 ) + β 2 ( x 2 1 ) + β 3 ( α 1 x ) 2 I 1 + β 4 ( α 2 x ) 2 I 2 + ε
where y = d 2 / D 2 , x = h / H , h is the height above ground to measurement point (m), d is the outside bark diameter at height h above ground (cm), H is the total tree height (m), D is the outside bark diameter at breast height (cm), β 1 to β 4 , α 1 , α 2 are parameters, I 1 and I 2 are bole portion indicator variables, and ε is an error term with zero expectation. In Equation (1), α 1 and α 2 are upper and lower joint point parameters, respectively, while I 1 = 1 if α 1 x and zero, otherwise, and I 2 = 1 if α 2 x and zero, otherwise. The value of α 1 indicates relative height where shape of the bole changes from paraboloid of the middle bole portion to conical of the top bole portion, and that of α 2 indicates where shape changes from neiloid of the butt portion to paraboloid of the middle portion of the tree bole.
The methods for parameter estimation of Equation (1) evolved from simple traditional least square regression and seemingly unrelated regression to nonlinear mixed-effects model. An implicit assumption shared by these methods is that the joint point parameters α 1 and α 2 are estimated without any constraint, which is in conflict with the condition of the feasible region (i.e., 0 < α 1 < α 2 < 1 ). Therefore, we applied a constrained two-dimensional optimum seeking method (CTOS) based on the residual sum of squares minimum theory to solve this problem. The corresponding objective function for estimating the parameters is given by
θ ^ = ( β ^ 1 , β ^ 2 , β ^ 3 , β ^ 4 , α ^ 1 , α ^ 2 ) = arg  min 0 < α 1 < α 2 < 1 ( f ( θ ) = i = 1 n | y y ^ | 2 )
where, θ ^ = ( β ^ 1 , β ^ 2 , β ^ 3 , β ^ 4 , α ^ 1 , α ^ 2 ) is the estimates of equation parameter vector θ, f ( θ ) is the residual sum of squares; y and y ^ are observed and estimated values, respectively.

2.3. CTOS & LSR Method

To solve for θ ^ in Equation (2), the CTOS & LSR method based on the optimum seeking method was proposed in this study. The optimum seeking method was first proposed by American mathematician Jack Keifer in the 1950s when studying the single-factor optimum seeking methods such as the golden section method (also called the 0.618 method). The purpose of the optimum seeking method was to reduce the number of experiments and to rapidly find out the optimal experimental points by using mathematical principles in the production and scientific research [37]. The CTOS method belongs to multivariate optimization approach, in which the effects of two variables or factors on the objective function are considered, and an optimal point is sought under a given constraint condition. The general idea of the CTOS method was as follows (see Figure 3): (I) The optimal point was first sought on a midline (e.g., red line); (II) Suppose the optimum was taken at point P 1 , and a horizontal line passing through P 1 is drawn (green line); (III) Then the optimum was sought on this horizontal line (green line); (IV) Suppose the optimum occurs at point P 2 , and an optimum was found on the vertical line (yellow line) passing through P 2 ; (V) The optimum was set at point P 3 ; (VI) The above procedures were repeated several times. The shadow indicated that this section was removed from the feasible region. Therefore, the feasible region was gradually narrowed until the final optimal point was found. For CTOS & LSR, the joint point parameters α 1 and α 2 were first estimated iteratively in the CTOS method, and then other parameters such as β 1 to β 4 in Equation (1) were estimated given α ^ 1 and α ^ 2 applying LSR. All estimates of these parameters in Equation (1) in each iteration were finally used to compute the values of the objective function (residual sum of squares).
The golden section method which is the foundation of many optimization algorithms was applied in this study to improve the computational efficiency of CTOS & LSR [37]. Two straight lines in each iteration were drawn at 0.382 and 0.618 of the value ranges of variable 1 and variable 2, respectively. The intersections of these straight lines were called interior points indicated by G 1 , G 2 , G 3 , and G 4 , which formed a rectangle U (see Figure 4). The values of objective functions f ( G 1 ) , f ( G 2 ) , f ( G 3 ) and f ( G 4 ) were calculated. The interior point that had the minimum of the objective function was denoted as temp point. The rectangle U with the new interior points G 1 , G 2 , G 3 , and G 4 was retained according to the temp point. After each iteration, the area of the rectangle U was gradually decreased, until the desired precision was obtained (Figure 4). A four-step iterative algorithm given below was applied to implement the CTOS & LSR.
Step 1:
The feasible region of variable 1 was ( a 1 , b 1 ) and the feasible region of variable 2 was ( a 2 , b 2 ). The four corners of a rectangle U formed by the overlapped area were represented by P 1 , P 2 , P 3 , and P 4 , respectively. In rectangle U , four interior points G 1 , G 2 , G 3 , and G 4 were obtained from the intersections of the following four lines:
x 1 = a 1 + 0.382 ( b 1 a 1 ) ,
x 2 = a 1 + 0.618 ( b 1 a 1 ) ,
y 1 = a 2 + 0.382 ( b 2 a 2 ) ,
y 2 = a 2 + 0.618 ( b 2 a 2 ) ,
The distance from the point P 1 with coordinates ( a 1 , a 2 ) to the point P 4 with coordinates ( b 1 , b 2 ) was calculated by the following formula:
D = p 1 p 4 2 = ( a 1 a 2 ) 2 + ( b 1 b 2 ) 2 .
Step 2:
θ in Equation (1) was estimated at each interior point with known coordinates by LSR, the objective functions f ( G 1 ) , f ( G 2 ) , f ( G 3 ) and f ( G 4 ) were computed, and the interior point that had the minimum objective function value was denoted as temp point.
Step 3:
If f ( G 1 ) was the minimum, then the rectangle with P 1 G 4 as the secondary diagonal was reserved; if f ( G 2 ) was the minimum, the rectangle with G 3 P 2 as the secondary diagonal was reserved; if f ( G 3 ) was the minimum, the rectangle with P 3 G 2 as the secondary diagonal was reserved; if f ( G 4 ) was the minimum, the rectangle with G 1 P 4 as the secondary diagonal was reserved. The four corners of the new rectangle U were denoted as P 1 , P 2 , P 3 , and P 4 , respectively. U , the four interior points G 1 , G 2 , G 3 , and G 4 in the new rectangle also were determined by the golden section method. Computed distance would be D = P 1 P 4 2 .
Step 4:
Steps 2 and 3 were repeated until a desired precision was obtained, for example, | D | < 0.00000001 . The final α ^ 1 and α ^ 2 were obtained, and then β ^ 1 β ^ 4 were estimated by LSR.
The CTOS & LSR with the four-step iterative algorithm for estimating in Equation (2) was implemented in VB version 6.0, see Appendix. It should be noted that, to simplify the LSR calculation, the model parameters β 1 β 4 in each iteration after giving α ^ 1 and α ^ 2 were estimated by transforming Equation (1) into the following linear form:
y = β 1 x 1 + β 2 x 2 + β 3 x 3 + β 4 x 4 + ε
where, x 1 = h / H 1 , x 2 = ( h / H ) 2 1 , x 3 = ( α ^ 1 h / H ) 2 I 1 , x 4 = ( α ^ 2 h / H ) 2 I 2 , y = ( d / D ) 2 , and the corresponding indicator variables I 1 and I 1 are given by
I 1 = | ( α ^ 1 h / H ) | + ( α ^ 1 h / H ) 2 ( α ^ 1 h / H )
I 2 = | ( α ^ 2 h / H ) | + ( α ^ 2 h / H ) 2 ( α ^ 2 h / H )
Therefore, Equation (8) was a typical linear model, which was easily solved by LSR.

2.4. ULSR Method

The method ULSR based on the minimization of the residual sum of squares has been commonly used to estimate parameters in forest growth and yield models. Equation (1) was difficult to solve directly using ULSR due to the unknown parameters α 1 and α 2 in both indicator variables I 1 and I 2 . Based on Equations (9) and (10), we transformed Equation (1) to the following form:
y = β 1 ( x 1 ) + β 2 ( x 2 1 ) + β 3 ( α 1 x ) 2 | ( α 1 x ) | + ( α 1 x ) 2 ( α 1 x ) + β 4 ( α 2 x ) 2 | ( α 2 x ) | + ( α 2 x ) 2 ( α 2 x ) + ε
This transformed Equation (11) is a general nonlinear model and could be directly estimated by ULSR.

2.5. Taper Equation Evaluation

With application of Equation (1) the individual tree-level equation using sub-dataset from nine isolated dominant or codominant trees (Figure 2) and the population average-level equation using the full dataset of 120 trees (Table 1) were developed using both ULSR and CTOS & LSR. The predictive ability of Equation (1) at the individual tree-level obtained from the two methods was assessed using the coefficient of determination ( R 2 ) and the residual sum of squares (RSS) statistics that are accounted for in Equations (12) and (13). The predictive ability of Equation (1) at the population average-level was assessed using R 2 , RSS, the mean of prediction errors ( e ¯ ) (Equation(14)), the variance of prediction errors ( σ 2 ) (Equation (15)), and the root mean square error ( δ ) (Equation (16)) by a defined cross-validation [38,39].
R ¯ 2 = 1 ( y i y ^ i ) 2 ( y i y ¯ ) 2
RSS = ( y i y ^ i ) 2
e ¯ = ( y i y ^ i ) / N
σ 2 = ( y i y ^ i ) 2 / ( N 1 )
δ = e ¯ 2 + σ 2
where y i and y ^ i are the observed and predicted dependent for the i th observation; y ¯ is the mean of the observed dependent, N is the total number of observations; e ¯ and σ 2 are the mean and variance of prediction errors, respectively; and δ is the root mean square error that combines the mean bias and the variation of the residuals and was used as the primary criterion for taper equation evaluations [40].
In this study, the cross-validation meant that the full dataset was randomly divided into 10 sub-datasets, each having 12 trees and Equation (1) was fitted 10 times. Each time, one sub-dataset (i.e., a total of 12 trees) was removed from the full dataset and the remaining nine sub-datasets (i.e., a total of 108 trees) together were used to fit the equation using both methods of ULSR and CTOS & LSR. The equations obtained were then employed to make tree taper predictions for the removed sub-dataset and the differences between the estimated and observed values were calculated. After all the sub-datasets yielded the predicted values of tree taper, the statistics ( R 2 , RSS , e ¯ , σ 2 and δ ) of the differences between the estimated and observed values were calculated. For details of the cross-validation approach, readers can refer to the studies by Vanclay [38] and Arlot and Celisse [39]. Moreover, all calculations were carried out using ForStat 2.2 version [41].

3. Results

3.1. Taper Equation at Individual Tree-Level

Table 2 shows the fitting results including parameter estimates and statistics for nine individual trees using both ULSR and CTOS & LSR. The fitting of Equation (1) using both ULSR and CTOS & LSR for each tree reached convergence with meaningful parameter estimates. The fitting of Equation (1) using the ULSR did not reach convergence with one-time starting values. Thus, different starting values were needed to make the equation convergence with the global minimum. However, with CTOS & LSR, Equation (1) easily converged for all nine individual trees. The differences of some parameter estimates of Equation (1) between ULSR and CTOS & LSR for all trees were significant at a risk level of α = 0.05. For example, for E. fordii, the estimate of the joint point parameter α 1 = 0.2431 from ULSR was almost four times larger than that ( α 1 = 0.0603 ) from CTOS & LSR, and the estimate of joint point parameter α 2 = 0.2818 from ULSR was almost one third of the corresponding parameter estimate ( α 1 = 0.8665 ) from CTOS & LSR.
A t-test showed that the mean biases of all trees (Table 2) for the two methods were not significantly different from zero (p > 0.05). The coefficients of determination R 2 of each tree for the two methods were very high. For many trees, such as C. hystrix No. 1 and No. 2, E. fordii No. 3, T. grandis No. 2 and No. 3, the two methods had identical fitting precisions. However, for several trees, such as E. fordii No. 2 and T. grandis No. 1, the RSS and R 2 from the CTOS & LSR method were slightly better than those of the corresponding trees obtained from ULSR. Especially for E. fordii No. 2, the value of RSS from CTOS & LSR was 40% smaller than that from ULSR. These results suggested that the fitting accuracy of Equation (1) at the individual tree-level obtained from CTOS & LSR was slightly higher than that of the corresponding equation obtained from ULSR for those trees whose taper curve has more variation.
For testing CTOS & LSR, the residuals of Equation (1) with the whole feasible ranges ( 0 < a 1 < a 2 < 1 ) of the two joint point parameters α 1 and α 2 by step size 0.01 and satisfying the constrained conditions α 1 < α 2 were calculated for the nine individual trees. The residual distribution of each tree is displayed in Figure 5. The red points in Figure 5 denote the corresponding final solutions of Equation (1) obtained by CTOS & LSR for each individual tree. The points with the minimum R S S of Equation (1) with the constrained 0 < a 1 < a 2 < 1 were almost identical to the corresponding red points that were determined by CTOS & LSR for the nine individual trees. This indicated that the CTOS & LSR proposed in this study could minimize the RSS of Equation (1) more efficiently.

3.2. Taper Equation at Population Average-Level

The parameter estimates of Equation (1) obtained from both ULSR and CTOS & LSR based on the full dataset (Table 1) were listed in Table 3. Equation (1) for each tree species reached convergence with meaningful parameter estimates using either method. As with Equation (1) at the individual tree-level, different starting values were required for ULSR to make Equation (1) converge. The parameter estimates of Equation (1) with ULSR and CTOS & LSR for each species were not identical. Especially for C. hystrix, the value of β 3 estimated by ULSR was 7% larger than that estimated by CTOS & LSR.
Based on cross-validation, the assessment results of five prediction statistics [Equations (12) and (16)] of Equation (1) using both ULSR and CTOS & LSR for the three tree species are listed in Table 4. The mean prediction errors e ¯ of Equation (1) for each tree species using either method were very small. Equation (1) for each tree species using both ULSR and CTOS & LSR slightly over-predicted the d 2 / D 2 . The coefficient of determination R 2 and RSS of Equation (1) for each tree species from ULSR and CTOS & LSR were correspondingly identical. The prediction accuracies of Equation (1) for each tree species from CTOS & LSR were much higher than those of Equation (1) when the d 2 / D 2 was estimated using from ULSR. For example, the values of σ 2 and δ of Equation (1) from CTOS & LSR for each tree species were smaller than those of Equation (1) from ULSR. Especially for C. hystrix, the variance of prediction errors σ 2 from CTOS & LSR was 8.34% smaller than those from ULSR, and the value of root mean square error δ from CTOS & LSR was 5% smaller than those from ULSR. This suggested that the prediction ability of Equation (1) using CTOS & LSR was more powerful than that of Equation (1) using ULSR.
The residuals of Equation (1) for each tree species using both ULSR and CTOS & LSR were plotted against the observed values (Figure 6). The results showed that Equation (1) for each tree species using either ULSR or CTOS & LSR did not lead to any trend of heteroscedasticity, which further indicated that Equation (1) was a better segmented taper equation for the prediction of tree taper.

4. Discussion

Equation (1) is able to describe the shape of a tree. The bottom trunk section, which includes the butt, is modeled as a frustum of a neiloidal solid, the middle of the trunk is assumed to take the shape of a frustum of a parabolodial solid, and the top portion is assumed to be conoid [13]. In this system, three additive equations were grafted together by incorporating two joint points between the three segments. The joint point parameters α 1 and α 2 with the constrained condition 0 < α 1 < α 2 < 1 in Equation (1) were previously commonly estimated by ULSR and SE, which made the results of the equation uncertain. We proposed the CTOS & LSR method to estimate the parameters of Equation (1). It was based on the constrained two-dimensional optimum seeking algorithm with Hua Luogeng’s theory [36]. The CTOS & LSR method was evaluated using the stem analysis data sets of three tropical precious tree species C. hystrix, E. fordii, and T. grandis in the southwest of China. Our results (Table 2 and Table 4) show that this method is more reliable than ULSR for estimating the parameters in the segmented taper equation.
The literature show that ULSR is a widely used method for parameter estimation [15,24,28,29,30,31,32]. This method was also applied in the CTOS & LSR. For each iterative calculation in CTOS & LSR, the CTOS was first used to estimate parameters α 1 and α 2 , after the two joint point parameters were determined, other remaining parameters in Equation (1) were estimated by LSR without constraint. The parameter estimates of Equation (1) for the two methods were obtained by minimizing the similar objective function (residual sum of squares). However, for ULSR, the estimates of parameters may fall out of the feasible range and make the results meaningless. For example, for the individual tree of E. fordii (Figure 7), the fitting accuracies of Equation (1) obtained from ULSR ( R S S = 0.0460 , R 2 = 0.9791 ), and CTOS & LSR ( R S S = 0.0391 , R 2 = 0.982 ) were almost similar and very high (Figure 7). However, the estimate of α 2 from ULSR is α ^ 2 = 1.3130 (the other estimate was α ^ 1 = 0.1915 ), which was larger than one and out of the feasible range. For CTOS & LSR, the estimates of α 1 and α 2 were 0.2692 and 0.6943; they were within the feasible range. Figure 7 also shows that the predicted taper curve from ULSR almost tends to zero when the relative height h / H is just 0.92, which conflicts with the assumption that the relative outside bark diameter d / D should be equal to zero at the relative height h / H = 1 . This situation always occurred in our preliminary analyses.
The overall prediction accuracies of Equation (1) at individual tree-levels and population average-levels from CTOS & LSR were much larger than those from ULSR (Table 2 and Table 4). This is due to the fact that the optimum solutions of Equation (1) obtained from CTOS & LSR were closer to actual values than these obtained from ULSR (Figure 5). In addition, the ULSR approach does not always converge for parameter estimation and this problem becomes more serious while fitting the taper equation at the individual tree-level. This is because the convergence of ULSR based on the iteration optimization algorithm, i.e., Gauss-Newton, depends largely on the given initial values of the parameters. In contrast, CTOS & LSR does not have a divergence problem. We also tried using other existing constrained optimization method e.g., simplex method and constrained Newton Raphson method, to estimate the parameters in Equation (1), but these methods failed to converge for some individual trees or tree species. Especially for the simplex method, the computational speed was also very slow. Besides forestry, CTOS & LSR could be also applied in agricultural engineering, industry or other fields requiring optimal selection of parameters. It is to be noted that the CTOS & LSR method only works on unimodal models. For multi-modal functions, there is no guarantee that the solutions of CTOS & LSR are the global optimum. We will focus on this limitation of the method in our future work.
Equation (1) with the estimates of parameters obtained from CTOS & LSR (Table 3) could be used to demonstrate differences of tree taper among three species C. hystrix, E. fordii, and T. grandis. The prediction accuracy with Equation (1) at the population average-level for each tree species is high (Figure 8). In addition, there are apparent differences among the tree taper of the three tropical precious tree species. For example, when the range of h / H was from 0.00 to 0.03, the values of d / D for T. grandis were much larger than those for the other two species (C. hystrix and E. fordii), which have almost the same d / D . When the range of h / H was from 0.03 to 0.7, the values of d / D for C. hystrix was the largest, following by E. fordii and T. grandis. When h / H is larger than 0.8, the three species have the same tree taper. Equation (1) also could be used to describe differences in tree taper between individual trees.
Measuring tree diameter and height may be subject to errors, even though stand and tree variables are commonly assumed to be measured without error [42,43]. The measurement errors made by field crew or faulty instrument or both might be substantial [42]. In all existing taper equations including Equation (1), it is assumed that (i) the dependent variable is a random variable; and (ii) the independent variables are fixed and measured without errors. It is well known that the violation of the second assumption may lead to biased estimation of parameters and (or) of the standard errors of the parameters, and consequently to misleading the hypothesis test [44,45]. If the predictor variables in Equation (1) are considered to have measurement errors, a new parameter estimation approach needs to be developed. We are in the process of developing such a method to deal with this problem.

5. Conclusions

In this study, the CTOS & LSR was proposed as an improved method to estimate the parameters in the Equation (1). This method was compared with ULSR for both individual tree-level equation and the population average-level equation using data from three tree species including C. hystrix, E. fordii, and T. grandis in the southwest of China. The differences between CTOS & LSR and ULSR were found to be significant. The segmented taper equation estimated using CTOS & LSR resulted in not only increased prediction accuracy, but also guaranteed the parameter estimates in a more meaningful way. It is thus recommended that the CTOS & LSR should be a preferred choice for this application.

Acknowledgments

We are grateful to Shouzheng Tang, Yuancai Lei, Yuanchang Lu and other researchers in the Chinese Academy of Forestry for their valuable suggestions on this study. We also thank the Experimental Center of Tropical Forestry, Chinese Academy of Forestry Sciences for providing friendly help during the data collection. Financial support for this study was provided by the National High-tech R & D Program of China (863 Program) (No. 2012AA102002), the Basic Scientific Research Business of Central Public Research Institutes (No. IFRIT2013), the Forestry Public Welfare Scientific Research Project of China (No. 201404417) and the National Natural Science Foundations of China (Nos. 31300534, 31570628, 31470641).

Author Contributions

L.F. and L.P. conceived the study. L. F. L.P., Y.M. amd X.S. performed the analysis and wrote the initial draft of the paper. All authors contributed to interpreting results and improvement of the paper.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Appendix

A VB program for estimating the Max and Burkhart ([13] 1976) taper Equation (1) using CTOS & LSR with the four-step iterative algorithm is shown here.
Function Do_Function_Two_dimensional_optimization_selection(Independ_Col() As Long, dePendent_Col() As Long, Independ_Name() As String, dePendent_Name() As String, a_ValueRegion() As Double, b_ValueRegion() As Double, startPoint As Long, EndPoint As Long)
Dim Independ_value() As Double, dePendent_value() As Double, n As Long, i As Long, j As Long, xyData() As Double, Para_value() As Double, Mat_result() As Variant
     n = EndPoint − startPoint + 1
     ReDim Independ_value(n, UBound(Independ_Col))
     ReDim dePendent_value(n, UBound(dePendent_Col))
   For i = 1 To n
     For j = 1 To UBound(Independ_Col)
       Independ_value(i, j) = frmData.DataTab1.TextMatrix(i + startPoint − 1, Independ_Col(j))
     Next j
     For j = 1 To UBound(dePendent_Col)
       If frmData.DataTab1.TextMatrix(i + startPoint − 1, dePendent_Col(j)) = "" Then
         frmData.DataTab1.TextMatrix(i + startPoint − 1, dePendent_Col(j)) = "00"
       End If
         dePendent_value(i, j) = frmData.DataTab1.TextMatrix(i + startPoint − 1, dePendent_Col(j))
     Next j
   Next i
     ReDim xyData(1 To n, UBound(Independ_Col) + UBound(dePendent_Col))
     For i = 1 To n
   For j = 1 To UBound(Independ_Col)
     xyData(i, j) = Independ_value(i, j)
   Next j
     Next i
     For i = 1 To n
   For j = 1 To UBound(dePendent_Col)
     xyData(i, j + UBound(Independ_Col)) = dePendent_value(i, j)
   Next j
     Next i
     ReturnVal = two_dimensional_optimization_selection_Tang(xyData(), a_ValueRegion(), b_ValueRegion(), Para_value())
     ReDim Mat_result(1, 7)
     Mat_result(0, 1) = "a1" : Mat_result(0, 2) = "a2" : Mat_result(0, 3) = "b1"
     Mat_result(0, 4) = "b2" : Mat_result(0, 5) = "b3" : Mat_result(0, 6) = "b4"
     Mat_result(0, 7) = "E2" : Mat_result(1, 1) = ReturnVal(1) : Mat_result(1, 2) = ReturnVal(2)
     Mat_result(1, 3) = Para_value(1) : Mat_result(1, 4) = Para_value(2)
     Mat_result(1, 5) = Para_value(3) : Mat_result(1, 6) = Para_value(4)
     Mat_result(1, 7) = ReturnVal(0)
     With frmOutput.Text1
     Call TextPrint(Mat_result, PrintFormat0, , , 0)
     End With
     Call Do_WriteOutputTxt
 Exit Function
	
Function two_dimensional_optimization_selection_Tang(xyData() As Double, a_ValueRegion() As Double, b_ValueRegion() As Double, Para_value() As Double)
Dim i As Integer, j As Integer, Sa1 As Double, Sb1 As Double, Sa2 As Double, Sb2 As Double, x1 As Double, x2 As Double, Y1 As Double, y2 As Double, Nx1 As Double, Nx2 As Double, Ny1 As Double, Ny2 As Double, r1 As Double, r2 As Double, R3 As Double, R4 As Double, temp As Double, dis As Double, Best1 As Double, Best2 As Double, ER As Double, reVal() As Double
     Sx1 = a_ValueRegion(1)
     Sx2 = a_ValueRegion(2) 'the feasible range for the first variable
     Sy1 = b_ValueRegion(1)
     Sy2 = b_ValueRegion(2) 'the feasible range for the second variable
     x1 = Sx1
     x2 = Sx2
     Y1 = Sy1
     y2 = Sy2
dis = ((x1 − x2)^2 + (Y1 − y2)^2)^0.5
While (dis >0.0000000001)
      Nx1 = x1 + (x2 − x1) × 0.382
      Nx2 = x1 + (x2 − x1) × 0.618
      Ny1 = Y1 + (y2 − Y1) × 0.382
      Ny2 = Y1 + (y2 − Y1) × 0.618
      r1 = ErrorC(Nx1, Ny2, xyData, Para_value)
      r2 = ErrorC(Nx2, Ny2, xyData, Para_value)
      R3 = ErrorC(Nx1, Ny1, xyData, Para_value)
      R4 = ErrorC(Nx2, Ny1, xyData, Para_value)
      If (r1 <r2) Then
        temp = r1
      Else
        temp = r2
      End If
      If (R3 <temp) Then
        temp = R3
      End If
      If (R4 <temp) Then
        temp = R4
      End If
      If (temp = r1) Then
        Y1 = Ny1
        x2 = Nx2
        Best1 = Nx1
        Best2 = Ny2
      ElseIf (temp = r2) Then
        x1 = Nx1
        Y1 = Ny1
        Best1 = Nx2
        Best2 = Ny2
      ElseIf (temp = R3) Then
        x2 = Nx2
        y2 = Ny2
        Best1 = Nx1
        Best2 = Ny1
       Else
        x1 = Nx1
        y2 = Ny2
        Best1 = Nx2
        Best2 = Ny1
      End If
      dis = ((x1 − x2)^2 + (Y1 − y2)^2)^0.5
      ER = ErrorC(Best1, Best2, xyData, Para_value)
          ReDim reVal(2)
          reVal(0) = ER
          reVal(1) = Best1
          reVal(2) = Best2
          two_dimensional_optimization_selection_Tang = reVal
          Exit Function
 End Function
		
Function ErrorC(Best1 As Double, Best2 As Double, xyData() As Double, Para_value() As Double)
Dim ReturnVal() As Double, i As Integer, j As Integer, iniParValues() As Double, X() As Double, Y() As Double, n As Integer, x1() As Double, x2() As Double,x3() As Double, x4() As Double, EstVal() As Variant, ResVal() As Variant, Sum_Res As Double, Rel() As Double
        n = UBound(xyData, 1)
        ReDim x1(n): ReDim x2(n): ReDim x3(n): ReDim x4(n)
        For i = 1 To n
      x1(i) = xyData(i, 1) − 1
    x2(i) = xyData(i, 1)^2 − 1
    If Best1 >= xyData(i, 1) Then
       x3(i) = (Best1 − xyData(i, 1))^2/(Best1 − x1(i))
    Else
       x3(i) = 0
    End If
    If Best2 >= xyData(i, 1) Then
       x4(i) = (Best2 − xyData(i, 1))^2
    Else
       x4(i) = 0
    End If
       Next i
       ReDim X(n, 4)
       For i = 1 To n
    X(i, 1) = x1(i)
    X(i, 2) = x2(i)
    X(i, 3) = x3(i)
    X(i, 4) = x4(i)
        Next i
        ReDim Y(n)
        For i = 1 To n
          Y(i) = xyData(i, 2)
        Next i
        Rel = LinReg(X, Y)
     ReDim Para_value(4)
     For i = 1 To 4
        Para_value(i) = Rel(i, 1)
     Next i
     ReDim EstVal(n)
     ReDim ResVal(n)
     Sum_Res = 0
     For i = 1 To n
        EstVal(i) = 0
      For j = 1 To 4
        EstVal(i) = EstVal(i) + Para_value(j) × X(i, j)
      Next j
        ResVal(i) = (Y(i) − EstVal(i))^2
        Sum_Res = Sum_Res + ResVal(i)
     Next i
        ErrorC = Sum_Res
        Exit Function
   End Function
   

References

  1. Demaerschalk, J.P. Converting volume equations to compatible taper equations. For. Sci. 1972, 18, 241–245. [Google Scholar]
  2. Duan, A.; Zhang, S.; Zhang, X.; Zhang, J. Development of a stem taper equation and modelling the effect of stand density on taper for Chinese fir plantations in Southern China. PeerJ 2016, 4, e1929. [Google Scholar] [CrossRef] [PubMed]
  3. Kozak, A. A variable-exponent taper equation. Can. J. For. Res. 1988, 18, 1363–1368. [Google Scholar] [CrossRef]
  4. Fang, Z.; Borders, B.E.; Bailey, R.L. Compatible volume taper models for loblolly and slash pine based on system with segmented-stem form factors. For. Sci. 2000, 46, 1–12. [Google Scholar]
  5. Jordan, L.; Berenhaut, K.; Souter, R.; Daniels, R.F. Parsimonious and completely compatible taper, total and merchantable volume models. For. Sci. 2005, 51, 578–584. [Google Scholar]
  6. Kozak, A.; Munro, D.O.; Smith, J.H.G. Taper functions and their application in forest inventory. For. Chron. 1969, 45, 278–283. [Google Scholar] [CrossRef]
  7. Ormerod, D.W. A simple bole model. For. Chron. 1973, 49, 136–138. [Google Scholar] [CrossRef]
  8. Amidon, E.L. A general taper functional form to predict bole volume for five mixed-conifer species in Califomia. For. Sci. 1984, 30, 166–171. [Google Scholar]
  9. Biging, G.S. Taper equations for second-growth mixed conifers in northern California. For. Sci. 1984, 30, 1103–1117. [Google Scholar]
  10. Thomas, C.E.; Parresol, B. Simple, flexible, trigonometric taper equations. Can. J. For. Res. 1991, 21, 1132–1137. [Google Scholar] [CrossRef]
  11. Sharma, M.; Oderwald, R.G. Dimensionally compatible volume and taper equations. Can. J. For. Res. 2001, 31, 797–803. [Google Scholar] [CrossRef]
  12. Zakrzewski, W.T.; MacFarlane, D.W. Regional stem profile model for cross-border comparisons of harvested red pine (Pinus resinosa Ait.) in Ontario and Michigan. For. Sci. 2006, 52, 468–475. [Google Scholar]
  13. Max, T.A.; Burkhart, H.E. Segmented polynomial regression applied to taper equations. For. Sci. 1976, 22, 283–289. [Google Scholar]
  14. Cao, Q.V.; Burkhart, H.E.; Max, T.A. Evaluations of two methods for cubic volume prediction of loblolly pine to any merchantable limit. For. Sci. 1980, 26, 71–80. [Google Scholar]
  15. Valenti, M.; Cao, Q. Use of crown ratio to improve loblolly pine taper equations. Can. J. For. Res. 1986, 16, 1141–1145. [Google Scholar] [CrossRef]
  16. Parresol, B.; Hotvedt, J.; Cao, Q. A volume and taper prediction system for bald cypress. Can. J. For. Res. 1987, 17, 250–259. [Google Scholar] [CrossRef]
  17. Jiang, L.; Liu, R. Segmented taper equations with crown ratio and stand density for Dahurian Larch (Larix gmelinii) in Northeastern China. J. For. Res. 2011, 22, 347–352. [Google Scholar] [CrossRef]
  18. Newnham, R.M. Variable-form taper equation for four Alberta tree species. Can. J. For. Res. 1992, 22, 210–223. [Google Scholar] [CrossRef]
  19. Muhairwe, C.K. Tree form and taper variation over time for interior lodgepole pine. Can. J. For. Res. 1994, 24, 1904–1913. [Google Scholar] [CrossRef]
  20. Valentine, H.T.; Gregoire, T.G. A switching model of bole taper. Can. J. For. Res. 2001, 31, 1400–1409. [Google Scholar] [CrossRef]
  21. Lee, W.K.; Seo, J.H.; Son, Y.M.; Lee, K.H.; von Gadow, K. Modeling stem profiles for Pinus densiflora in Korea. For. Ecol. Manag. 2003, 172, 69–77. [Google Scholar] [CrossRef]
  22. Zeng, W.; Liao, Z. Research of taper equation. Sci. Silvae Sin. 1997, 33, 127–132, (In Chinese with English Summary). [Google Scholar]
  23. Demaerschalk, J.P.; Kozak, A. The whole-bole system: A conditioned dual equation system for precise prediction of tree profiles. Can. J. For. Res. 1997, 7, 488–497. [Google Scholar] [CrossRef]
  24. Cao, Q.V.; Wang, J. Calibrating fixed- and mixed-effects taper equations. For. Ecol. Manag. 2011, 262, 671–673. [Google Scholar] [CrossRef]
  25. Sabatia, C.O.; Burkhart, H.E. On the use of upper stem diameters to localize a segmented taper equation to new trees. For. Sci. 2015, 61, 411–423. [Google Scholar] [CrossRef]
  26. Rojo, A.; Perales, X.; Sanchez-Rodriguez, F.; Alvarez-Gonzalez, J.G.; von Gadow, K. Stem tape functions for maritime pine (Pinus pinaster Ait.) in Galicia (Northwestern Spain). Eur. J. For. Res. 2005, 124, 177–186. [Google Scholar] [CrossRef]
  27. Jiang, L.; Brooks, J.R.; Wang, J. Compatible taper and volume equations for yellow-poplar in West Virginia. For. Ecol. Manag. 2005, 213, 399–409. [Google Scholar] [CrossRef]
  28. Figueiredo-Filho, A.; Borders, B.E.; Hitch, K.L. Taper equations for Pinus taeda plantations in Southern Brazil. For. Ecol. Manag. 1996, 83, 39–46. [Google Scholar] [CrossRef]
  29. Muhairwe, C.K. Taper equations for Eucalyptus pilularis and Eucalyptus grandis for the north coast in New South Wales, Australia. For. Ecol. Manag. 1999, 113, 251–269. [Google Scholar] [CrossRef]
  30. Sharma, M.; Burkhart, H.E. Selecting a level of conditioning for the segmented polynomial taper equation. For. Sci. 2003, 49, 324–330. [Google Scholar]
  31. Coble, D.W.; Hilpp, K. Compatible cubic-foot stem volume and upper-stem diameter equations for semi-intensive plantation grown loblolly pine trees in East Texas. South J. Appl. For. 2006, 30, 132–141. [Google Scholar]
  32. Trincardo, G.; Burkhart, H.E. A generalized approach for modeling and localizing stem profile curves. For. Sci. 2006, 52, 670–682. [Google Scholar]
  33. Brooks, J.R.; Jiang, L.; Ozcelik, R. Compatible stem volume and taper equations for Brutian Pine, Taurus Cedar, and Taurus Fir in Turkey. For. Ecol. Manag. 2008, 256, 147–151. [Google Scholar] [CrossRef]
  34. Pang, L. Stand Harvest Analysis and Information Management System for Target Tree. Ph.D. Thesis, Chinese Academy of Forestry, Beijing, China, 2015. [Google Scholar]
  35. Tang, S.; Li, Y.; Wang, Y. Simultaneous equations, errors-in-variable models, and model integration in systems ecology. Ecol. Model. 2001, 142, 285–294. [Google Scholar] [CrossRef]
  36. Hua, L. Optimum Seeking Methods; Science Press: Beijing, China, 1967. [Google Scholar]
  37. Schwefel, P.H.P. Evolution and Optimum Seeking: The Sixth Generation; John Wiley & Sons, Inc.: New York, NY, USA, 1993. [Google Scholar]
  38. Vanclay, J.K. Modelling Forest Growth and Yield: Applications to Mixed Tropical Forests; CAB International: Wallingford, UK, 1994. [Google Scholar]
  39. Arlot, S.; Celisse, A. A survey of cross-validation procedures for model selection. Statist. Surv. 2010, 4, 40–79. [Google Scholar] [CrossRef]
  40. Yang, Y.; Huang, S. Comparison of different methods for fitting nonlinear mixed forest models and for making predictions. Can. J. For. Res. 2011, 41, 1671–1686. [Google Scholar] [CrossRef]
  41. Tang, S.Z.; Lang, K.J.; Li, H.K. Statistics and Computation of Biomathematical Models (ForStat Course); Science Press: Beijing, China, 2008. (In Chinese) [Google Scholar]
  42. Omule, A.Y. Personal bias in forest measurement. For. Chron. 1980, 56, 222–224. [Google Scholar] [CrossRef]
  43. Gertner, G.Z. The sensitivity of measurement error in stand volume estimation. Can. J. For. Res. 1990, 20, 800–804. [Google Scholar] [CrossRef]
  44. Fuller, W.A. Measurement Error Models; John Wiley and Sons: New York, NY, USA, 1987. [Google Scholar]
  45. Rencher, A.C.; Schaalje, G.B. Linear Models in Statistics, 2nd ed.; John Wiley and Sons: New York, NY, USA, 2008. [Google Scholar]
Figure 1. Location of study area: Experimental Center of Tropical Forestry, Chinese Academy of Forestry Sciences and spatial distribution of 30 sample plots.
Figure 1. Location of study area: Experimental Center of Tropical Forestry, Chinese Academy of Forestry Sciences and spatial distribution of 30 sample plots.
Forests 07 00194 g001
Figure 2. Observed height (h) above the ground-outside bark diameter (d) curves for nine dominant or codominant trees.
Figure 2. Observed height (h) above the ground-outside bark diameter (d) curves for nine dominant or codominant trees.
Forests 07 00194 g002
Figure 3. The general idea of the constrained two-dimensional optimum seeking method, a1 and b1 are the lower and upper limits for the feasible region of the first parameter, respectively; a2 and b2 are the lower and upper limits for the feasible region of the second parameter, respectively, P1, P2, and P3 are corresponding optimal points for different iterative steps.
Figure 3. The general idea of the constrained two-dimensional optimum seeking method, a1 and b1 are the lower and upper limits for the feasible region of the first parameter, respectively; a2 and b2 are the lower and upper limits for the feasible region of the second parameter, respectively, P1, P2, and P3 are corresponding optimal points for different iterative steps.
Forests 07 00194 g003
Figure 4. Diagram of the combination of constrained two-dimensional optimum seeking method and least square regression (CTOS & LSR), P 1 , P 2 , P 3 , and P 4 are four corners of the rectangle U determined by the feasible region of the estimated two-dimensional variables, G 1 , G 2 , G 3 , and G 4 are interior points of U determined by the golden section method, P 1 , P 2 , P 3 , and P 4 are four corners of the rectangle U obtained by one iteration calculation of CTOS & LSR; G 1 , G 2 , G 3 , and G 4 are interior points of U determined by the golden section method.
Figure 4. Diagram of the combination of constrained two-dimensional optimum seeking method and least square regression (CTOS & LSR), P 1 , P 2 , P 3 , and P 4 are four corners of the rectangle U determined by the feasible region of the estimated two-dimensional variables, G 1 , G 2 , G 3 , and G 4 are interior points of U determined by the golden section method, P 1 , P 2 , P 3 , and P 4 are four corners of the rectangle U obtained by one iteration calculation of CTOS & LSR; G 1 , G 2 , G 3 , and G 4 are interior points of U determined by the golden section method.
Forests 07 00194 g004
Figure 5. The residual distribution of Equation (1) with the whole feasible range ( 0 < α 1 < α 2 < 1 ) of joint points α 1 and α 2 by step size 0.01 and satisfying the constrained conditions α 1 < α 2 calculated by the combination of the constrained two-dimensional optimum seeking method and the least square regression (CTOS & LSR) for the nine individual trees. The red points denote the corresponding final solutions of Equation (1) obtained by CTOS & LSR for each individual tree.
Figure 5. The residual distribution of Equation (1) with the whole feasible range ( 0 < α 1 < α 2 < 1 ) of joint points α 1 and α 2 by step size 0.01 and satisfying the constrained conditions α 1 < α 2 calculated by the combination of the constrained two-dimensional optimum seeking method and the least square regression (CTOS & LSR) for the nine individual trees. The red points denote the corresponding final solutions of Equation (1) obtained by CTOS & LSR for each individual tree.
Forests 07 00194 g005
Figure 6. Residuals of the Equation (1) for the three tree species C. hystrix, E. fordii. and T. grandi using both the combination method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) and the unconstrained least square regression (ULSR) plotted against the observed values based on cross-validation.
Figure 6. Residuals of the Equation (1) for the three tree species C. hystrix, E. fordii. and T. grandi using both the combination method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) and the unconstrained least square regression (ULSR) plotted against the observed values based on cross-validation.
Forests 07 00194 g006
Figure 7. Observed relative outside bark diameters at different relative height location (black line) for an individual tree of E. fordii. The predicted taper curves from the Equation (1) using both the combination method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) (red line) and the unconstrained least square regression (ULSR) (blue line) are shown.
Figure 7. Observed relative outside bark diameters at different relative height location (black line) for an individual tree of E. fordii. The predicted taper curves from the Equation (1) using both the combination method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) (red line) and the unconstrained least square regression (ULSR) (blue line) are shown.
Forests 07 00194 g007
Figure 8. The predicted taper curves from Equation (1) using the combined method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) for the three tree species C. hystrix (red line), E. fordii (black line), and T. grandi (blue line) at the population average level.
Figure 8. The predicted taper curves from Equation (1) using the combined method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) for the three tree species C. hystrix (red line), E. fordii (black line), and T. grandi (blue line) at the population average level.
Forests 07 00194 g008
Table 1. Summary statistics of full dataset.
Table 1. Summary statistics of full dataset.
Tree SpeciesVariableNo.MeanMaxMinSDCV
C. hystrixD (cm)4022.1237.558.707.400.33
H (m)4021.4330.5010.105.210.24
E. fordiiD (cm)4021.1529.5010.054.520.21
H (m)4018.6422.4011.802.710.15
T. grandiD (cm)4025.9337.9512.705.490.21
H(m)4021.1627.409.303.520.17
D, diameter at breast height; H, total tree height; SD, standard deviation; CV, coefficient of variation.
Table 2. Parameter estimates and statistics of Equation (1) fitted by a combined method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) and the unconstrained least square regression (ULSR) for nine individual trees (three trees from each species of C. hystrix, E. fordii, and T. grandis), α 1 α 2 , joint point parameters; β 1 β 4 , model parameters; RSS , residual sum of square; R 2 , coefficient of determination.
Table 2. Parameter estimates and statistics of Equation (1) fitted by a combined method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) and the unconstrained least square regression (ULSR) for nine individual trees (three trees from each species of C. hystrix, E. fordii, and T. grandis), α 1 α 2 , joint point parameters; β 1 β 4 , model parameters; RSS , residual sum of square; R 2 , coefficient of determination.
Tree Species No.Method α 1 α 2 β 1 β 2 β 3 β 4 RSS R 2
C. hystrix(1)CTOS & LSR0.05790.3937−2.87961.3293139.5926−3.3680.02300.991
ULSR0.05850.3937−2.87911.329−3.3658129.3540.02300.991
(2)CTOS & LSR0.03520.5093−3.21361.6851182.6908−1.75460.00900.996
ULSR0.03560.5092−3.2131.6847173.213−1.75450.00900.996
(3)CTOS & LSR0.22920.3783−2.11260.853825.3279−7.45830.00740.997
ULSR0.28120.3407−2.17120.891232.6352−22.270.00660.997
E. fordii(1)CTOS & LSR0.15950.7361−4.72642.389427.0427−2.79650.02040.989
ULSR0.16490.7675−5.91433.060921.8090−3.44460.02060.989
(2)CTOS & LSR0.24310.2818−2.17160.959193.5096−53.88560.02990.985
ULSR0.06030.86652.0566−1.421292.61922.33320.04990.9774
(3)CTOS & LSR0.03400.3631−2.64551.3478886.3202−1.70840.00290.999
ULSR0.03420.3695−2.65541.3541844.7019−1.6670.00290.999
T. grandi(1)CTOS & LSR0.39960.5837−4.40122.34427.9341−5.25710.01790.990
ULSR0.49160.5549−4.09752.16519.0758−9.54580.01900.990
(2)CTOS & LSR0.13850.6910−2.89391.418933.1479−1.25890.00260.999
ULSR0.14030.8530−8.48214.500027.9491−4.25670.00260.999
(3)CTOS & LSR0.03250.5609−3.49951.7276558.3056−2.29660.01340.995
ULSR0.03280.5611−3.50111.7286532.4962−2.29680.01340.995
Table 3. Parameter estimates of Equation (1) obtained from both ULSR and CTOS & LSR based on the full data-set from three species (C. hystrix, E. fordii, and T. grandi), α 1 α 2 , joint point parameters; β 1 β 4 , model parameters.
Table 3. Parameter estimates of Equation (1) obtained from both ULSR and CTOS & LSR based on the full data-set from three species (C. hystrix, E. fordii, and T. grandi), α 1 α 2 , joint point parameters; β 1 β 4 , model parameters.
Tree SpeciesMethod α 1 α 2 β 1 β 2 β 3 β 4
C. hystrixULSR0.04810.6479−3.62511.7049158.3459−2.1074
CTOS & LSR0.04740.6479−3.6241.7043169.9877−2.1063
E. fordiiULSR0.04440.7190−2.59921.1148257.2861−0.7817
CTOS & LSR0.04400.6910−2.42561.0154274.0969−0.6923
T. grandisULSR0.09030.6993−2.59231.128251.1553−0.9210
CTOS & LSR0.08900.6910−2.51931.086152.6617−0.8752
Table 4. The statistics of Equation (1) fitted by the combination method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) and the unconstrained least square regression (ULSR) for three trees species (C. hystrix, E. fordii, and T. grandis) based on cross-validation approach.
Table 4. The statistics of Equation (1) fitted by the combination method of constrained two-dimensional optimum seeking and least square regression (CTOS & LSR) and the unconstrained least square regression (ULSR) for three trees species (C. hystrix, E. fordii, and T. grandis) based on cross-validation approach.
Tree SpeciesMethod RSS R 2 e ¯ σ 2 δ
C. hystrixULSR3.4040.959−0.25182.28981.5340
CTOS & LSR3.4040.959−0.18142.09891.4601
E. fordiiULSR4.4010.953−0.15682.18781.4874
CTOS & LSR4.4010.953−0.16112.17531.4837
T. grandisULSR3.4790.962−0.22113.73281.9446
CTOS & LSR3.4790.962−0.18703.54521.8921
RSS , residual sum of square; R 2 , coefficient of determination; e ¯ , mean of prediction errors; σ 2 , variance of prediction errors; δ , root mean square error.

Share and Cite

MDPI and ACS Style

Pang, L.; Ma, Y.; Sharma, R.P.; Rice, S.; Song, X.; Fu, L. Developing an Improved Parameter Estimation Method for the Segmented Taper Equation through Combination of Constrained Two-Dimensional Optimum Seeking and Least Square Regression. Forests 2016, 7, 194. https://doi.org/10.3390/f7090194

AMA Style

Pang L, Ma Y, Sharma RP, Rice S, Song X, Fu L. Developing an Improved Parameter Estimation Method for the Segmented Taper Equation through Combination of Constrained Two-Dimensional Optimum Seeking and Least Square Regression. Forests. 2016; 7(9):194. https://doi.org/10.3390/f7090194

Chicago/Turabian Style

Pang, Lifeng, Yongpeng Ma, Ram P. Sharma, Shawn Rice, Xinyu Song, and Liyong Fu. 2016. "Developing an Improved Parameter Estimation Method for the Segmented Taper Equation through Combination of Constrained Two-Dimensional Optimum Seeking and Least Square Regression" Forests 7, no. 9: 194. https://doi.org/10.3390/f7090194

APA Style

Pang, L., Ma, Y., Sharma, R. P., Rice, S., Song, X., & Fu, L. (2016). Developing an Improved Parameter Estimation Method for the Segmented Taper Equation through Combination of Constrained Two-Dimensional Optimum Seeking and Least Square Regression. Forests, 7(9), 194. https://doi.org/10.3390/f7090194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop