Next Article in Journal
Changes in the Soil Bacterial Community in a Chronosequence of Temperate Walnut-Based Intercropping Systems
Previous Article in Journal
GIS-Based Multi-Criteria Assessment and Seasonal Impact on Plantation Forest Landscape Visual Sensitivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Error Structure for Additive Biomass Equations on the Use of Multivariate Likelihood Function

Research Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Forests 2019, 10(4), 298; https://doi.org/10.3390/f10040298
Submission received: 27 January 2019 / Revised: 26 March 2019 / Accepted: 26 March 2019 / Published: 30 March 2019
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)

Abstract

:
Research Highlights: this study developed additive biomass equations respectively from nonlinear regression (NLR) on original data and linear regression (LR) on a log-transformed scale by nonlinear seemingly unrelated regression (NSUR). To choose appropriate regression form, the error structures (additive vs. multiplicative) of compatible biomass equations were determined on the use of the multivariate likelihood function which extended the method of likelihood analysis to the general occasion of a contemporaneously correlated set of equations. Background and Objectives: both NLR and LR could yield the expected predictions for allometric scaling relationship. In recent studies, there are vigorous debates on which regression (NLR or LR) should apply. The main aim of this paper is to analyze the error structure of a compatible system of biomass equations to choose more appropriate regression. Materials and Methods: based on biomass data of 270 trees for three tree species, additive biomass equations were developed respectively for NLR and LR by NSUR. Multivariate likelihood functions were computed to determine the error structure based on the multivariate probability density function. The anti-log correction factor which kept the additive property was obtained separately using the arithmetic and weighted average of basic correction factors from each equation to assess two model specifications on the comparably original scale. Results: the assumption of additive error structure was well favored for an additive system of three species based on the joint likelihood function. However, the error structure of each component equation calculated from the conditional likelihood function for compatible equations might be different. The performance of additive equations corrected by a weighted average of basic correction factor from each component equation performed better than that of the arithmetic average and held good property of compatibility after corrected. Conclusions: NLR provided a better fit for additive biomass equations of three tree species. Additive equations which confirmed the responding assumption of error structure performed better. The joint likelihood function on the use of the multivariate likelihood function could be used to analyze the error structure of the additive system which was a result of a tradeoff for each component equation. Based on the average of correction factors from each component equation to correct the bias of additive equations was feasible for the hold of additive property, which might lead to a poor correction effect for some component equation.

1. Introduction

Allometric research characterizes the scaling relationship between various response variables and different measures of body size, which has been dominant for many years in a variety of different areas, such as physiology, numerical ecology, and morphology [1,2,3]. Kittredge (1944) [4] described the biomass of tree components with tree dimension variables based on an allometric equation to quantify the tree biomass in the form of Y = a X b , where Y is tree component biomass, X is tree dimension variable and a, b respectively represents allometric coefficient and exponent. Up to date, thousands of biomass equations have been developed for various tree species and regions all over the world for the purpose of accurate quantification of forest biomass dealing with carbon reduction and climate change [5,6,7]. However, researchers witness a heated issue recently regarding fitting methods which concentrate largely on the topic, linear regression on log-transformed data (hereafter, LR) with a multiplicative error in arithmetic domain or nonlinear regression on original scale (hereafter, NLR) with an additive error.
For decades, LR was the most commonly adopted pattern in allometric research. The conventional practice is to fit a straight line from log-transformed data using ordinary least square and then to back-transform the resulting equation to yield the estimate on the arithmetic scale [8,9,10,11]. Nonetheless, the effectiveness and accuracy of applying LR have been subject to criticism mainly because of the following aspects: (1) Back-transformation from a straight line fitted to logarithm obtained the geometric means for prediction values instead of arithmetic means, which decreased the estimation on the original scale using direct back-transformation [12,13,14]. Although this bias from anti-logarithm could be modified by a certain form of correction factor [8,11,15], some research argued that using anti-log correction factor might cause overestimation [16,17]. (2) While log-transformation could stabilize the variance, it produced an insidious rotational distort for allometric equations which created a new distribution that differed in a fundamental way from the original scale [18,19]. (3) This nonlinear distort unduly emphasized on small values but compressed large-individual values which led to a poor fit for the end of the curve graphically [20,21,22]. (4) The artificial transformation might cause outliers undetected which made the data favorable [19,22,23]. Generally, the focus of controversy for allometric equations fitted by LR lied in the injudicious use of log-transformation [14,22,24].
NLR, directly fitted to the original data by iteration method for allometric equations, has been broadly used by more and more researchers because of convenient and user-friendly statistical software [25,26,27]. However, heteroscedasticity of arithmetic values fitted by NLR directly is of general occurrence which may fail to satisfy statistical assumptions [13,24]. Nonetheless, researches have shown that heteroscedasticity of observations does not necessarily invalidate the deterministic equation fitted by NLR [28] and even the failure of satisfying the constancy of variance, it performed better than LR which yielded more accurate estimates on the original scale [22,29]. It was worthily noted that a weight factor could address the problem of heteroscedasticity by generalized least square as well as the log-transformation of LR [26]. Nonetheless, the debate on which fitting method (NLR or LR) performs better and which error structure confirms the statistical assumption more appropriately has not subsided.
Xiao et al. [30] and Ballantyne [31] proposed the approach of likelihood analysis to determine the error structure (multiplicative vs. additive) for allometric equations so that the suitable fitting procedure (LR or NLR) could be adopted. Recently, the likelihood analysis has come to be applied in the area of forestry and ecology. Lai et al. [32] used the likelihood analysis to compare the allometry of coarse root biomass from LR and NLR for Castanopsis eyrei (Champ. ex Benth.) Tutch., Schima superba Gardn. et Champ., Pinus massonoana Lamb., and mixed species and concluded the empirical data supported a multiplicative error. Ma and Jiang [33] applied the likelihood analysis to determine the error structure of individual tree volume model for Larix gmelinii (Ruprecht) Kuzeneva. and Pinus sylvestris Linn. var. mongolica Litv. which supported the multiplicative error, but the comparison of model assessment indicated NLR performed better than LR. Dong et al. [34] adopted the likelihood analysis to determine the error structure of compatible or additive biomass equations for three conifer species in Northeast China, which favored the multiplicative error. However, the proposed approach of determination on error structure by Xiao et al. [30] and the following application including additive equations developed by Dong et al. [34] were all based on the one-dimension likelihood function which was considered only appropriate for a single allometric equation. For a compatible system of several equations, there were significant contemporaneous correlations when it was simultaneously estimated. Therefore, the analysis based on the one-dimension likelihood function seems to be unreasonable when it applies to determine the error structure for additive biomass equations.
An additive system of biomass equations ensures the logically equal relationship that the predictions for the components sum to the predictions from a total equation. To achieve the additivity property, there are different methods to develop compatible equations. At first, the total predictions could be obtained simply from the sum of components equations developed independently to ensure the additivity [35,36,37]. Up to date, simultaneous estimation for a system of equations widely known as seemingly unrelated regression (SUR) has been broadly used for a compatible system of biomass equations [26,34,37,38]. Back-transformation from a straight line (LR) that fits the logarithm to the original scale could introduce the systematic bias. To remove and reduce this bias, researchers have computed different forms of correction factors and compared the corrected effects. However, little information could be provided when the additive biomass equation was developed by LR. Dong et al. [34] used the correction factor separately from each equation to correct the bias of each component, which did not take the additivity property into account. To our knowledge, the correction factor for additive biomass equation has not been reported, which corrects the bias from anti-log transformation and ensures the additivity property at the same time.
Cinnamomum camphora (L.) Presl, Schima superba Gardn. et Champ. and Liquidambar formosana Hance are widely distributed in Southeastern China and are also the dominant broad-leaved tree species in Guangdong province. There are many differences in morphology and physiology between broad-leaved and conifer tree species. But researches on the biomass equations centered mostly on the conifer species while there are limited studies on broad-leaved tree species [6,27,34,39]. The purpose of our study is (1) to develop a compatible system of biomass equations between branch, foliage, stem wood, stem bark and total aboveground for three broad-leaved tree species separately based on NLR and LR by SUR, (2) to compute the multivariate likelihood function for determination on the error structure of additive biomass equations which extend the method of likelihood analysis to the general situation of a contemporaneously correlated set of equations, (3) to formulate the correction factor for a compatible system of biomass equations to correct the bias introduced by anti-logarithm transformation and ensure the additivity property at the same time, (4) to compare the fitting result of two procedures based on NLR and LR on the same arithmetic scale and evaluate the effect of assumption for different error structures on the result of model fitting.

2. Material and Methods

2.1. Data Collection

Tree dimension variable and biomass data, including Cinnamomum camphora, Schima superba and Liquidambar formosana covering whole Guangdong province in Southeastern China with 90 trees for each species, were sampled in 2013 by Guangdong Forestry Survey and Planning Institute (Figure 1). The sample trees were classified by the diameter class of 2 cm, 4 cm, 6 cm, 8 cm, 12 cm, 16 cm, 20 cm, 26 cm, 32 cm and 38 cm (above 38 cm). Among them, 60 trees were evenly distributed following above 10 diameter classes with six trees for each class and the remaining 30 trees were chosen based on the actual distribution of diameter class and the number of trees from the 8th National Forest Inventory in Guangdong province.
The destructive sampling procedure was processed for the living sample trees avoiding severe defects. Before the tree was felled at the ground level height, the diameter at breast height (D, at 1.3 m aboveground) was measured. After felled, the living crown was evenly marked into three parts (top, middle, and bottom) and weighted separately, then the branches and leaves from each part summing to about 500–1000 g of fresh mass was randomly sampled and placed in a labeled bag for moisture content determination. The stem was also marked into three sections including 0–2/10, 2/10–5/10 and above 5/10 tree height and weighed separately. At each section of the stem, a 2–3 cm thick disk separately from the upper and lower part was cut and weighed, then taken to the laboratory for moisture content determination. All samples were dried at 85 °C to constant weight. The dry biomass of each component was calculated by multiplying the fresh weight of each component by the dry/fresh ratio of each component sample. The total foliage biomass was the sum of foliage dry biomass. The total stem wood biomass was the sum of all stem wood’s sections dry mass. The total stem bark biomass is the sum of all stem bark dry mass. The aboveground biomass was the sum of branch, foliage, stem wood, and stem bark dry biomass. The above procedure of moisture content determination was conducted by the laboratory center of College of Forestry and Landscape Architecture, South China Agriculture University according to the related technical regulations [40]. The data statistics were summarized for 90 sampling trees of each broad-leaved species in Table 1.

2.2. Model Specification and Estimation

To fit the allometric equation, either NLR on the arithmetic scale or LR on the logarithmic scale could yield the estimation values. The fundamentally substantial difference between these two approaches largely relies on the assumption of how error term manifests in the equation, which is known as the error structure (Xiao et al., 2011) [30]. NLR assumes the equation with the normally additive error on the arithmetic scale such that:
Y = a x b + ε , ε ~ N ( 0 , σ 2 )
In contrast, LR assumes that the error is normally distributed and additive on the logarithmic scale such that:
log Y = log a + b log x + ε , ε ~ N ( 0 , σ 2 )
which corresponds to log-normally distributed, multiplicative error on the arithmetic scale,
Y = a x b × exp ( ε ) , ε ~ N ( 0 , σ 2 )
To determine which model specification was the most appropriate for a compatible system of biomass equations for three broad-leaved tree species in this study, two model forms that correlated with additivity among four component biomass equations and total aboveground biomass equation were specified as follows with cross-equation constraints on the structural parameters:
(1) The first model specification assumes the error structure is additive (Equation (1)) and a compatible system of five biomass equations as follows:
W B R = a 1 · D b 1 + ε 1 W F L = a 2 · D b 2 + ε 2 W S W = a 3 · D b 3 + ε 3 W S B = a 4 · D b 4 + ε 4 W A B = a 1 · D b 1 + a 2 · D b 2 + a 3 · D b 3 + a 4 · D b 4 + ε 5
(2) The second model specification assumes the error structure is multiplicative on the arithmetic scale (Equation (3)) and logarithmic transformation was taken on both sides of equations (Equation (2)) such that
l o g W B R = log a 1 + b 1 log D + ε 1 l o g W F L = log a 2 + b 2 log D + ε 2 l o g W S W = log a 3 + b 3 log D + ε 3 l o g W S B = log a 4 + b 4 log D + ε 4 l o g W A B = log ( a 1 · D b 1 + a 2 · D b 2 + a 3 · D b 3 + a 4 · D b 4 ) + ε 5
where W B R ,   W F L ,   W S W ,   W S B and W A B represent branch biomass, foliage biomass, stem wood biomass, stem bark biomass and the aboveground biomass in kg, respectively, D is the diameter at breast height in cm, log denotes natural logarithm, a i and b i are regression coefficient for Equation (4), log a i is the intercept and b i is the regression coefficient for Equation (5), ε i and ε i are the equation error terms for NLR and LR additive model, respectively.
Above two model specifications for additive biomass equations were estimated using nonlinear seemingly unrelated regression generally known as NSUR. The logarithmic transformation tends to balance the heteroscedastic variance. For comparison, Equation (4) was fitted to data using weighted NSUR as demonstrated by Parresol (2001) [26] to stabilize the variance. The weight of each component equation was obtained by the weight function w = f ( D ) 1 , where f ( D ) was the prediction value for estimated equation [26,41,42].

2.3. Multivariate Likelihood Function to Analyze Error Structure

Xiao et al. [30] outlined the approach of likelihood analysis to facilitate the objective determination of the error structure based on the single one-dimension likelihood function. When applied to the additive biomass model with cross-equation correlation, it seems. Considering that, in this study, we computed the multivariate likelihood function including the joint likelihood function and the conditional likelihood function to respectively analyze model system and each component equation for the correlated error structure of additive biomass equations. Based on the joint probability density function, the joint likelihood function can be calculated by:
(1) For the p components system of NLR (Equation (4)), the joint likelihood function that the data are generated from a normal distribution with additive error is calculated as follows:
L N L R ( p ) = j = 1 n [ 1 ( 2 π ) p 2 | Σ N L R | 1 / 2 exp [ 1 2 ( w i j · ( Y i j a N L R X i j b N L R ) ) Σ N L R 1 ( w i j · ( Y i j a N L R X i j b N L R ) ) ] ]
(2) For the p components system of LR (Equation (5)), the joint likelihood function that the data are generated from a lognormal distribution with multiplicative error on the arithmetic scale:
L L R ( p ) = j = 1 n [ 1 y i · ( 2 π ) p 2 | Σ L R | 1 / 2 exp [ 1 2 ( log Y i j log ( a L R X i j b L R ) ) Σ L R 1 ( log Y i j log ( a L R X i j b L R ) ) ] ]
where L N L R ( p ) , L L R ( p ) are the joint likelihood function of NLR and LR for p component equations, respectively. Where n is the sample size. X i j ,   Y i j (i = 1, …, p; j = 1, …, n) are the jth value for predictor and response variable of the ith component equation, respectively. Σ N L R , Σ L R are the error variance-covariance matrix of NLR and LR, respectively. | Σ N L R | , | Σ L R | are the determinant of responding matrix respectively, w i j is the weight of the jth predicted value of the ith component equation. a N L R , b N L R , b N L R and b L R are the responding regression coefficient for NLR and LR, respectively.
According to the definition of conditional distribution for the multivariate probability density function [43], the conditional likelihood function for ith component equation can be defined as follows:
L N L R ( L N L R ( i ) | L N L R ( p i ) ) = L N L R ( p ) / L N L R ( p i )
L L R ( L L R ( i ) | L L R ( p i ) ) = L L R ( p ) / L L R ( p i )
where L N L R ( L N L R ( i ) | L N L R ( p i ) ) , L L R ( L L R ( i ) | L L R ( p i ) ) is the conditional likelihood function for ith component equation of NLR and LR respectively, L N L R ( i ) , L L R ( i ) is the likelihood function for ith component equation of NLR and LR respectively. L N L R ( p i ) , L L R ( p i ) is the value of the joint likelihood function calculated from (pi) components without the use of ith component equation for NLR and LR, respectively.
To compare different candidate models fitted to the same dataset statistically, Akaikes’s Information Criterion (AIC) can be used to evaluate the goodness-of-fit of a model by involving both the likelihood and a penalty for extra parameters. The lowest value for AIC identifies the candidate model conveying the most information about the relationship between predictor and response. AICc which is a second-order variant of AIC for small sample size is computed as
A I C c = 2 k 2 log L + 2 k ( k + 1 ) n k 1
where k is the number of parameters. L is the joint likelihood function for the model system (Equation (6) for NLR and Equation (7) for LR) or the conditional likelihood function for each component equation (Equation (8) for NLR and Equation (9) for LR). If A I C c n o r m A I C c l o g n < 2 , the assumption of additive error is favored and the result from Equation (4) should be processed. If A I C c n o r m A I C c l o g n > 2 , the assumption of multiplicative error is favored and the result from Equation (5) should be processed [44]. If | A I C c n o r m A I C c l o g n | 2 , neither of these two error structures is appropriate and model averaging is suggested. Besides the difference of AIC respectively from NLR and LR, evidence ratio (ER) (see Appendix A) was also taken to provide the evidence for the appropriate model selection [44].

2.4. Back-Transformed Correction Factor for Additive Equations

To obtain the arithmetic value of prediction, a correction factor (hereafter, CF) is commonly used to correct the systematic bias introduce by anti-log transformation from a straight line (Equation (2)) fitted to logarithmic data. For the additively log-transformed biomass equations (Equation (5)), not only the systematic bias should be corrected, but also the additivity property of the value of prediction from back-transformation need to be satisfied. Thus, based on the basic CF, we formulated the specific correction factor for a compatible system of biomass equations. The two basic CFs for ith component can be calculated as follows [8,15]:
C F i 1 = exp ( δ i i 2 / 2 )
C F i 2 = j = 1 n y i j / j = 1 n y ^ i j
where δ i i 2 is the ith diagonal element of the error variance-covariance matrix, y i j is the jth observed value for ith component, y ^ i j is the predicted value of the jth observed value for ith component. Then, the arithmetic and weighted average CF for the compatible system can be respectively obtained by
C F a t = 1 p i = 1 p C F i t           t = 1 ,   2
C F w t = 1 2 i = 1 p W i · C F i t           t = 1 ,   2
where C F a t , C F w t are the arithmetic and weighted average of the tth (t = 1, 2) basic correction factor from each component equation, respectively. C F i t is the tth (t = 1, 2) basic correction factor for ith component equation, W i is the proportion of the ith component biomass accounted for the total aboveground biomass.

2.5. Model Assessment

This study used the entire empirical data to fit additive biomass equations [45]. Model fitting and predicting was assessed by the statistics as follows.
Coefficient of determination
R 2 = 1 j = 1 n ( y j y ^ j ) 2 / j = 1 n ( y j y ¯ ) 2
Standard error of estimate
S E E = 1 n k j = 1 n ( y j y ^ j ) 2
Total relative error
T R E = j = 1 n ( y j y ^ j ) / j = 1 n y ^ j × 100 %
Average system error
A S E = 1 n j = 1 n [ ( y j y ^ j ) / y ^ j ] × 100 %
Relatively mean absolute error
R M A = 1 n j = 1 n | ( y j y ^ j ) / y ^ j | × 100 %
Mean prediction error
M P E = t α · ( S E E / y ¯ ) / n × 100 %
where y j , y ^ j are the jth observed value and the responding predicted value, y ¯ is the average of the observed value, k is the number of parameters, t α is the t value when the confidence level is α (usually taken by 95%).
To ensure that the estimated mean function captures dominant pattern in the arithmetic scale, the fitted model not only needs to be assessed by several statistics but also should be validated graphically, which was a critically important oversight by many researchers [23,46]. In this study, the additive biomass equations based on different assumptions of error structure (additive vs. multiplicative) were validated graphically.

3. Result

3.1. Error Structure for Each Component Equation and Additive System

Additive biomass equations were fitted to original and log-transformed data, respectively (Equation (4) and Equation (5)), to yield the parameter estimation value, then calculated the AICc, respectively, namely AICc-norm for Equation (4) and AICc-logn for Equation (5) based on the conditional likelihood function from the conditional probability density function. The difference between AICc value and parameter estimate were computed in Table 2. The AICc-norm for Schima superba and Liquidambar formosana was clearly lower than AICc-logn with a difference between −731.2 and −221.5 supporting the additive error for each component equation. But there existed different error structures for each component equation of Cinnamomum camphora. The AICc-norm for branch and foliage was larger than AICc-logn favoring the multiplicative error structure with a difference of 38.3 and 263.2 while for other component equation, the AICc-norm was lower than AICc-logn with a difference between −105.9 and −15.9 favoring the additive error structure. Components for two other tree species had a large ER (evidence ratio) more than 100 supporting the additive error as well, while the branch and foliage for Cinnamomum camphora had a smaller ER less than 0.01 supporting the multiplicative error as well.
The joint likelihood function was calculated for the whole model system based on the joint probability density function. The analysis of error structure for additive model system was shown in Table 3. The AICc-norm for three tree species was all lower than AICc-logn with a difference of −7.0, −810.2 and −846.7 and got relatively large ER as well, supporting the additive error, which meant that the approach of NLR for additive biomass equations was appropriate for three broad-leaved tree species, especially for Schima superba and Liquidambar formosana in this study.

3.2. Assessment of Anti-Log Correction Factor for Additive System

Log-transformed equation predicts the logarithm of the response variable. To obtain the unbiased value in the original scale, the anti-log correction factor is necessary. The arithmetic and weighted average ( C F a t and C F w t ) of basic correction factors from each equation as well as responding evaluation statistics for total aboveground biomass was listed in Table 4. The arithmetic average of basic CF from each component equation was represented by C F a t (t = 1, 2) and the weighted average of basic CF from each component equation was represented by C F w t (t = 1, 2). CF0 represented the model was not corrected.
The uncorrected model performed worse with a lower R2, larger standard error of estimate (SEE, hereafter) and mean prediction error (MPE, hereafter) than that of NLR for three tree species. Thus, NLR model could yield relatively better prediction compared to the uncorrected LR model. After applying the CF for Cinnamomum camphora, R2 of LR model improved by 0.018 to 0.025, SEE decreased by 6.07 kg to 8.66 kg and total relative error (TRE, hereafter), Relatively mean absolute error (RMA, hereafter) and MPE dropped in varying degrees. Importantly, R2 and SEE might not be the best with C F w 2 , but the remaining statistics including TRE, average system error (ASE, hereafter), RMA, and MPE were even better than NLR. LR model for Schima superba obtained worse fitting and predicting accuracy when it was corrected with R2 decreasing by 0.006 to 0.041 and SEE increasing by 1.73 kg to 10.91 kg, but among the different four CFs, the correction effect for C F w 2 was relatively better than that of other correction factors, and TRE, ASE and RMA statistics was better than that of NLR model. The different assessment statistics of corrected LR model got dropped and increased to different degrees for Liquidambar formosana. Using C F w 2 increased the R2 by 0.003, dropped the SEE by 1.43 kg compared with LR0, and the ASE and RMA statistics was better than NLR model, but in terms of R2 and SEE statistics, it was slightly worse than NLR model.
C F w t (t = 1, 2) corrected better for Schima superba and Liquidambar formosana especially the CF based on secondly basic correction factor (Equation (12)), that is C F w 2 . Although C F a t could yield higher R2, TRE and ASE got relatively larger and reached −3.82%, −1.35% and −8.42%, −6.06%, respectively. Generally speaking, the value of C F w t was larger than that of C F a t . The approach of weighted average apparently reduced TRE and ASE for additive biomass equations. Taking Liquidambar formosana as an example, based on two basic CFs using the approach of weighted average, TRE reduced by 8.08% and 8.12% while ASE reduced 7.99% and 8.02%. As far as all the evaluation statistics, C F w 2 corrected best for aboveground biomass of additive equations for three broad-leaved tree species. When CFw2 for additive system was used, it did not perform better than NLR model except Cinnamomum camphora which is slightly better than NLR, but the difference between NLR and corrected LR model was small for three tree species total aboveground biomass.

3.3. Comparison of Model Fitting and Error Structure

Major violation indicated the inappropriateness of the model and potential invalidity of the result. To assess the correction effect of C F w 2 for each component, the result of evaluation statistics was listed in Table 5. LR0 represented the uncorrected model and LRw2 represented the corrected model by C F w 2 . For Schima superba and Liquidambar formosana, the NLR model of each component including the total aboveground biomass obtained better estimation supporting the additive error which was consistent with the determination (Table 2 and Table 3). For Cinnamomum camphora, the NLR model for stem wood and stem bark yielded relatively better prediction with higher R2 and small SEE, TRE, RMA and MPE, favoring the additive error which was consistent with the determination (Table 2). Nonetheless, for the branch, R2 for LRw2 model was slightly higher than NLR, but it got the worst TRE and ASE, while NLR model performed better for foliage component. In addition, the LR model for total aboveground biomass corrected by C F w 2 got a relatively better fitting result (Table 4) favoring the multiplicative error which was slightly inconsistent with the determination (Table 2). It could be seen that the error structure of component equation for Cinnamomum camphora was different and this would be discussed later in detail.
Observed values of biomass components together with the fitted curve were plotted against diameter at the breast for three tree species, respectively (Figure 2, Figure 3 and Figure 4). All models showed a good fit to the small untransformed observations. There was no visually apparent difference between LR0 and LRw2 model for three tree species. Except for Schima superba branch and foliage as well as Liquidambar formosana foliage, the fitted curve from three model curves showed a close path following the path of data. Nonetheless, most NLR models estimated slightly larger than LR0 model especially for a large diameter and the larger the diameter was, the clearer this pattern exhibited. The mean function from NLR model could capture a relatively dominant pattern, which followed the path of the data, especially for the larger individual. The additive error structure for a compatible system was fairly more appropriate graphically to formulate the model specification and fit on the original scale which was consistent with the determination in Table 3.

4. Discussion

The one-dimension likelihood function was derived from univariate normal distribution and solved the estimate issue of a single function called maximum likelihood estimates (MLE), which was used later to determine the error structure of allometric equation by Xiao et al. (2011) [30] and Ballantyne (2013) [31]. In this study, each equation of a compatible system of additive equations was estimated simultaneously to ensure the significantly contemporaneous correlations by NSUR. However, the one-dimension likelihood function only considered the single variate but ignored the correlations of multiple variates in the additive equations, which might be inappropriate to determine the error structure of additive equations. In contrast, the multivariate likelihood function took the relationships of multiple variates into account and reflected the multivariate error distribution for additive equations more accurately.
This study computed the joint and conditional multivariate likelihood function for additive biomass equation of three broad-leaved tree species respectively based on the joint and conditional probability density function to analyze the error structure (additive vs. multiplicative) of each component equation and model system. The model satisfying the responding error structure fitted better and major violation indicated the inappropriateness of the model and potential invalidity of the result. The NLR model for Schima superba and Liquidambar formosana indeed yielded better estimation than uncorrected and corrected LR model statistically and graphically, which verified our determination on additive error structure properly in this study. However, for Cinnamomum camphora, the corrected model of total aboveground obtained more accurate estimated value than NLR, especially the total aboveground model corrected by C F w 2 which had six evaluation statistics relatively better than that of NLR while for the foliage component the NLR model performed better. It indicated that the error of total equation might be the additive, but the error of components was not necessarily the same. This is mainly due to the different error structures for each component equation determined by the conditional likelihood function (see Table 2), but either NLR or LR could be taken to estimate additive biomass equations. Nonetheless, to hold the property of compatibility for each component NSUR compromised the error among component equations [26,37]. The error structure for additive system based on the joint likelihood function was the result of a tradeoff for component equations and might cause the inconsistence of error structure between additive system and each component equation, leading to the determination on error structure and model assessment for aboveground and foliage component was inconsistent just like Cinnamomum camphora in this study.
Likelihood analysis based on AIC provided a method for analyzing the error structure to determine more appropriate regression (NLR or LR) especially for a compatible system of biomass equations [30]. Nonetheless, using AIC as a direct indicator to compare candidate regression equations (NLR or LR) has been criticized by some researchers [46]. Through the graphical validation of NLR and LR equation, Packard (2013) thought that AIC was not a sufficient way to choose alternative statistical models between NLR and LR regression. In addition, some researchers thought the individual AIC, AICc, or BIC values were not interpretable in absolute terms as they contain arbitrary constants and are much affected by sample size [46]. Evidence ratio rescaled these information criteria and was good evidence to compare candidate models, which overcame the shortcomings of direct comparison from AIC. It is noteworthy that the larger the difference between AIC, the larger the evidence ratio. Evidence ratio might be more appropriate to compare the model which had a close AIC that could not be directly differentiated in absolute terms. Moreover, as Packard (2014) [23] said, the good fit must capture the dominant pattern in the untransformed data, Figure 2, Figure 3 and Figure 4 in this study clearly indicate that the lower the AIC of an equation, the better the capture of pattern in all range of data. When considering the candidate equation, statistic test might not be enough to assess the appropriateness of fit and several criteria, as well as graphical validation, were quite necessary.
Both log-transformation for LR model and weighted estimation for NLR model can stabilize the heteroscedasticity and make the constant of variance. However, the log-transformation was thought to create a newly logarithmic scale to estimate parameters [10,14]. Thus, the value of back-transformed prediction would not reflect the real relationship and relied largely on which variance of each response value changed on the arithmetic scale [19,21]. In our study, it is noted that the uncorrected model LR0 substantially underestimated the predicted value, especially for large observed value (TRE much larger than zero and the curve was lower than others). When the log-transformed observed value did not fall on the real linear curve, it is easy to understand from log function curve that the back-transformed linear model put much weight on the predicted value for a small individual and compressed the predicted value for a large individual. This nonlinear transformation caused an accurate estimation for the small value and a poor estimation for the large estimation on an arithmetic scale.
To obtain the accurate predicted value closer to the arithmetic scale, a correction factor is necessary to correct the systemic bias introduced by log-transformation [8,10,13], but the compatible property of each component value summing to the total value is also needed for the additive system. Based on the two basic CFs, the arithmetic and weighted average from each component equation were computed in this study. The first basic CF (Equation (11)) has been the most used CF from the log-normal function but only satisfying the assumptions strictly, can it yield perfect correction effect [8]. Because of overcompensating the bias with the standard error of estimate, it might cause an overestimate [16,17]. The second basic CF is independent on the model distribution and corrects the bias from the observed and predicted value, which might cause the value to be lower than 1.0 [11,15]. Using second basic CF to formulate the systemic correct factor for additive biomass equations performed better in this study, which was consistent with the result proposed by Snowdon in 1991 [15]. In addition, weighted average applied the proportion of each component accounting for the total as a weight to calculate the CF, which considered the relationship among components. Thus, it could obtain a better correction effect compared with that of the arithmetic average in this study. But because the CF for the additive system was calculated from the average of each component, it might lead to a poor correction effect for some specific components. For example, the fitting effect got worse when it was corrected for Schima superba branch and the total aboveground.
NLR has become a commonly used approach with a feature of inexpensive, user-friendly software in allometric studies. So, does it imply NLR perform better definitely than LR model and does it mean the conventionally log-transformed model is unnecessary [13,19,24]? This debate on NLR (additive error) or LR (multiplicative error) model which one is better suited for allomeric research has never subsided. Because of the unbalanced weight put on the predicted value on the original scale for LR model and better fitting of a large value for NLR model, it is suggested that the LR model might be appropriate for small individuals, such as young forests, while NLR model might be appropriate for large individuals, such as mature forests. However, to choose a better model, both statistical analysis and graphical validation for the real empirical data are needed. This research provided a statistical analysis of the determination on the error structure for additive biomass equations. For a compatible model system, especially when the error structures of component equations were analyzed differently, describing the error structure accurately and improving the fitting accuracy could be an interesting research area in the future.

5. Conclusions

In this study, we developed the multivariate likelihood function to analyze the error structure for additive biomass equations of three broad-leaved tree species, which extended the likelihood function proposed by Xiao [30] and Ballantyne [31] to the general occasion of a contemporaneously correlated set of equations. To compare NLR and LR on the original scale, the correction factors specific for additive equations were developed by the arithmetic and weighted average of two basic corrections from each component equation to hold additive property. The main conclusion could be that: (1) the multivariate likelihood function could be used to analyze the error structure of additive biomass equations and the result of model assessment confirmed our determination. The conditional likelihood function could be used for component equations. The joint likelihood function could be used for additive system. The determination on error structure was a result of tradeoff for additive biomass equations and the error for total equation might be the additive but the error for components was not necessarily the same, (2) the correction factors developed in this study could yield a good effect of correction especially for the approach of weighted average based on the second CF (Equation (12)) which could be used for additive equations to hold the compatible property after corrected, (3) the additive equations confirming the responding error structure got more fitting accuracy while violating the responding assumption caused the accuracy loss. In this study, NLR got relatively better goodness-of-fit for the additive biomass equations of three broad-leaved tree species.

Author Contributions

H.L. provided the data and idea; L.C. analyzed the data and wrote the paper.

Funding

This research was funded by National Natural Science Foundation of China, grant number 31770676.

Acknowledgments

The study would not have been possible without the corporation of Guangdong Forestry Survey and Planning Institute and South China Agriculture University. The authors thank Yuancai Lei and my colleagues in Chinese Academy of Forestry for their suggestions. We also appreciate the valuable comments and the constructive suggestions from editors and anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Appendix A.1. Multivariate Normal Distribution

The multivariate normal density is a generalization of the univariate normal density to P ≥ 2 dimensions. Recall that the univariate normal distribution, with mean μ and variance σ2, has the probability density function
f ( x ) = 1 2 π σ 2 e x p { 1 2 ( ( x μ ) / σ ) 2 }
Assume that the P × 1 vectors X1, X2, …, Xn represent a random sample from a multivariate normal population with mean vector μ and covariance matrix Σ. The multivariate normal distribution with mean µ and covariance Σ > 0 has the density function (Equation (A2)) and can be written X ~ Np (μ, Σ).
f ( x ) = 1 ( 2 π ) p / 2 | Σ | 1 / 2 e x p { 1 2 ( x μ ) Σ 1 ( x μ ) }
The joint density function of all the observations is the product of the marginal normal densities:
{ J o i n t   d e n s i t y   o f X 1 , X 2 ,   ,   X n } = j = 1 n { 1 ( 2 π ) p / 2 | Σ | 1 / 2 e x p { 1 2 ( x j μ ) Σ 1 ( x j μ ) } }
when the numerical values of the observations become available, they may be substituted for the Xj in Equation (A3). The resulting expression now considered as a function of μ and Σ for the fixed set of observations is called the likelihood function.

Appendix A.2. Multivariate Log-Normal Distribution and Correction Factor

It is a well-known fact that if a random variable X is normally distributed with mean μ and variance σ2, then the moment generating function of X is
m ( θ ) = E ( e θ X ) = 1 2 π σ e θ X e 1 2 ( x μ σ ) 2 d x = e μ θ + 1 2 σ 2 θ 2
For a proof see a relative statistical book, such as Hogg and Craig (1995).
Let Y follow the lognormal distribution, so W = ln Y is normally distributed with mean of μ w and variance σ w 2 . Then the moment generating function of W is
m ( θ ) = E ( e θ W ) = e μ w θ + 1 2 σ w 2 θ 2
Since W = ln Y implies Y = e w , it follows that
E ( Y ) = E ( e W ) = m ( 1 ) = e μ w + 1 2 σ w 2
According to Baskerville (1972), the correct factor for back-transformation of univariate log-transformed equation defined as
C F = e x p ( S E E 2 / 2 )
where SEE represented the standard error of estimate of the regression which was the unbiased estimate of σ w 2 . Equation (A7) was calculated from the univariate log-transformed equation, but for the multivariate log-transformed equations, multivariate log-normal distribution should be considered.
Let X = [ X 1 , X 2 , , X p ] be a p-component random vector having normal distribution with mean μ and covariance matric Σ . Now we use the transformation Y i = exp ( X i ) and define Y = [ Y 1 , Y 2 , , Y p ] . The density of Y is multivariate log-normal distribution and has the following form:
f Y = 1 ( 2 π ) p / 2 | Σ | 1 / 2 Y 1 e x p { 1 2 ( l n Y μ ) Σ 1 ( l n Y μ ) }
where l n Y = [ l n Y 1 ,   l n Y 2 , , l n Y p ] is p-component column vector.
Mean of Yi and covariance matrix D i j respectively is
ν i = exp ( μ i + 1 2 σ i j 2 )
D i j = exp { ( μ i + μ j ) + 1 2 ( σ i i + σ j j ) } · { exp ( σ i j ) 1 }
where σ i j is the ijth element covariance of X.

Appendix A.3. Multivariate Conditional Distribution

Assume X(1) represent r × 1 random vector and X(2) represent (p − r) × 1 random vector. If p random vector X = ( X ( 1 ) X ( 2 ) ) , the conditional density of X(1) given that X(2) is defined by
f ( X ( 1 ) | X ( 2 ) ) = { c o n d i t i o n a l   d e n s i t y   X ( 1 ) o f g i v e n   t h a t   X ( 2 ) } = f ( X ( 1 ) , X ( 2 ) ) / f ( X ( 2 ) )
where f ( X ( 2 ) ) is the density function of X(2).

Appendix A.4. Toward Akaike Weights and Evidence Ratios

ER (evidence ratio) is imperative to rescale these criteria such as AIC. Usually, this is done by subtracting to AIC of each model AIC of the model with the minimum one:
Δ A I C i = A I C i A I C m i n
where A I C i represent AIC calculated from the ith (I = 1, 2, …, R) model, A I C m i n is the minimum AIC of R models.
It is convenient to normalize the model likelihoods such that they sum to 1 and treat them as probabilities. Hence, we use:
ω i = exp ( Δ A I C i / 2 ) i = 1 R exp ( Δ A I C i / 2 )
The ω i is called Akaike weights and the ratios ω i / ω j is called evidence ratio which is identical to the original likelihood ratios.
In this study, NLR and LR were compared using evidence ratio and ER was calculated by:
E R = ω N L R ω L R

References

  1. Bolte, A.; Rahmann, T.; Kuhr, M.; Pogoda, P.; Murach, D.; Gadow, K.V. Relationships between tree dimension and coarse root biomass in mixed stands of European beech (Fagus sylvatica L.) and Norway spruce (Picea abies [L.] Karst.). Plant Soil 2004, 264, 1–11. [Google Scholar] [CrossRef]
  2. West, G.B.; Brown, J.H.; Enquist, B.J. A general model for the origin of allometric scaling laws in biology. Science 1997, 276, 122–126. [Google Scholar] [CrossRef] [PubMed]
  3. Huxley, J.S. Problems of Relative Growth; Dial Press: New York, NY, USA, 1931; pp. 1–31. [Google Scholar]
  4. Kittredge, J. Estimation of the amount of foliage of trees and stands. J. For. 1944, 42, 905–912. [Google Scholar]
  5. Jenkins, J.C.; Chojnacky, D.C.; Heath, L.S.; Birdsey, R.A. National-scale biomass estimators for United States tree species. For. Sci. 2003, 49, 12–35. [Google Scholar]
  6. Ter-Mikaelian, M.T.; Korzukhin, M.D. Biomass equations for sixty-five north American tree species. For. Ecol. Manag. 1997, 97, 1–24. [Google Scholar] [CrossRef]
  7. Zianis, D.; Xanthopoulos, G.; Kalabokidis, K.; Kazakis, G.; Ghosn, D.; Roussou, O. Allometric equations for aboveground biomass estimation by size class for pinus brutia ten. trees growing in north and south Aegean islands, Greece. Eur. J. For. Res. 2011, 130, 145–160. [Google Scholar] [CrossRef]
  8. Baskerville, G.L. Use of Logarithmic Regression in the Estimation of Plant Biomass. Can. J. For. Res. 1972, 2, 49–53. [Google Scholar] [CrossRef]
  9. Brown, J.H.; Gillooly, J.F.; Allen, A.P.; Savage, V.M.; West, G.B. Toward a metabolic theory of ecology. Ecology 2004, 85, 1771–1789. [Google Scholar] [CrossRef]
  10. Smith, R.J. Allometric scaling in comparative biology: Problems of concept and method. Am. J. Physiol. 1984, 246, 152–160. [Google Scholar] [CrossRef] [PubMed]
  11. Smith, R.J. Logarithmic transformation bias in allometry. Am. J. Phys. Anthropol. 1993, 90, 215–228. [Google Scholar] [CrossRef]
  12. Beauchamp, J.J.; Olson, J.S. Corrections for bias in regression estimates after logarithmic transformation. Ecology 1973, 54, 1403–1407. [Google Scholar] [CrossRef]
  13. Gingerich, P.D. Arithmetic or geometric normality of biological variation: An empirical test of theory. J. Theor. Biol. 2000, 204, 201–221. [Google Scholar] [CrossRef] [PubMed]
  14. Zar, J.H. Calculation and miscalculation of the allometric equation as a model in biological data. Bioscience 1968, 18, 1118–1120. [Google Scholar] [CrossRef]
  15. Snowdon, P. A ratio estimator for bias correction in logarithmic regressions. Can. J. For. Res. 1991, 21, 720–724. [Google Scholar] [CrossRef]
  16. Madgwick, H.A.I.; Satoo, T. On estimating the aboveground weights of tree stands. Ecology 1975, 56, 1446–1450. [Google Scholar] [CrossRef]
  17. Zianis, D.; Mencuccini, M. Aboveground biomass relationships for beech (fagus moesiaca cz.) trees in vermio mountain, northern Greece, and generalized equations for fagus sp. Ann. For. Sci. 2003, 60, 439–448. [Google Scholar] [CrossRef]
  18. Glass, N.R. Discussion of calculation of power function with special reference to respiratory metabolism in fish. J. Fish. Res. Board Can. 1969, 26, 2643–2650. [Google Scholar] [CrossRef]
  19. Packard, G.C.; Boardman, T.J. Model selection and logarithmic transformation in allometric analysis. Physiol. Biochem. Zool. 2008, 81, 496–507. [Google Scholar] [CrossRef]
  20. Jansson, M. A comparison of detransformed logarithmic regressions and power function regressions. Geogr. Ann. Ser. A Phys. Geogr. 1985, 67, 61–70. [Google Scholar] [CrossRef]
  21. Packard, G.C. On the use of logarithmic transformations in allometric analyses. J. Theor. Biol. 2009, 257, 515–518. [Google Scholar] [CrossRef]
  22. Packard, G.C.; Birchard, G.F.; Boardman, T.J. Fitting statistical models in bivariate allometry. Biol. Rev. 2011, 86, 549–563. [Google Scholar] [CrossRef]
  23. Packard, G.C. Multiplicative by nature: logarithmic transformation in allometry. J. Exp. Zool. 2014, 322, 202–207. [Google Scholar] [CrossRef]
  24. Kerkhoff, A.J.; Enquist, B.J. Multiplicative by nature: Why logarithmic transformation is necessary in allometry. J. Theor. Biol. 2009, 257, 519–521. [Google Scholar] [CrossRef]
  25. Li, H.; Zhao, P. Improving the accuracy of tree-level aboveground biomass equations with height classification at a large regional scale. For. Ecol. Manag. 2013, 289, 153–163. [Google Scholar] [CrossRef]
  26. Parresol, B.R. Additivity of nonlinear biomass equations. Can. J. For. Res. 2001, 31, 865–878. [Google Scholar] [CrossRef]
  27. Wang, C. Biomass allometric equations for 10 co-occurring tree species in Chinese temperate forests. For. Ecol. Manag. 2006, 222, 9–16. [Google Scholar] [CrossRef]
  28. Finney, D.J. Was this in your statistics textbook? v. transformation of data. Exp. Agric. 1989, 25, 165–175. [Google Scholar] [CrossRef]
  29. Asselman, N.E.M. Fitting and interpretation of sediment rating curves. J. Hydrol. 2000, 234, 228–248. [Google Scholar] [CrossRef]
  30. Xiao, X.; White, E.P.; Hooten, M.B.; Durham, S.L. On the use of log-transformation vs. nonlinear regression for analyzing biological power laws. Ecology 2011, 92, 1887–1894. [Google Scholar] [CrossRef] [Green Version]
  31. Ballantyne, F.T. Evaluating model fit to determine if logarithmic transformations are necessary in allometry: A comment on the exchange between packard (2009) and kerkhoff and enquist (2009). J. Theor. Biol. 2013, 317, 418–421. [Google Scholar] [CrossRef]
  32. Lai, J.; Yang, B.; Lin, D.; Kerkhoff, A.J.; Ma, K. The allometry of coarse root biomass: Log-transformed linear regression or nonlinear regression? PLoS ONE 2013, 8, e77007. [Google Scholar] [CrossRef]
  33. Ma, Y.; Jiang, L. Error structure and variance function of allometric model. Sci. Silva Sin. 2018, 54, 90–97. [Google Scholar]
  34. Dong, L.; Zhang, L.; Li, F. A compatible system of biomass equations for three conifer species in northeast, china. For. Ecol. Manag. 2014, 329, 306–317. [Google Scholar] [CrossRef]
  35. Cunia, T.; Briggs, R.D. Forcing additivity of biomass tables: Some empirical results. Can. J. For. Res. 1984, 14, 376–384. [Google Scholar] [CrossRef]
  36. Kozak, A. Methods for ensuring additivity of biomass components by regression analysis. For. Chron. 1970, 46, 402–405. [Google Scholar] [CrossRef]
  37. Parresol, B.R. Assessing tree and stand biomass: A review with examples and critical comparisons. For. Sci. 1999, 45, 573–593. [Google Scholar]
  38. Bi, H.; Turner, J.; Lambert, M.J. Additive biomass equations for native eucalypt forest trees of temperate Australia. Trees 2004, 18, 467–479. [Google Scholar] [CrossRef]
  39. Zou, W.; Zeng, W.; Zhang, L.; Zeng, M. Modeling crown biomass for four pine species in China. Forests 2015, 6, 433–449. [Google Scholar] [CrossRef]
  40. State Forestry Administration of PR China. Technical Regulation on Sample Collections for Tree Biomass Modeling. In LY/T 2259-2014; Chinese Standard Press: Beijing, China, 2014. [Google Scholar]
  41. Zeng, W. Comparison of different weight function in weighted regression. For. Res. Manag. 2013, 5, 55–61. [Google Scholar]
  42. Zeng, W. Research on weighting regression and modeling. Sci. Silva Sin. 1999, 35, 5–11. [Google Scholar]
  43. Richard, A.J.; Dean, W.W. Applied Multivariate Statistical Analysis, 6th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2007; pp. 149–173. [Google Scholar]
  44. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2002; pp. 70–79. [Google Scholar]
  45. Kozak, A.; Kozak, R. Does cross validation provide additional information in the evaluation of regression models. Can. J. For. Res. 2003, 33, 976–987. [Google Scholar] [CrossRef]
  46. Packard, G.C. Is logarithmic transformation necessary in allometry? Biol. J. Linn. Soc. 2013, 109, 476–486. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The locations of trees for three broad-leaved tree species.
Figure 1. The locations of trees for three broad-leaved tree species.
Forests 10 00298 g001
Figure 2. The curve fitted from nonlinear regression (NLR), uncorrected linear regression (LR0) and corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 ((LRw2) for Cinnamomum camphora (Cinnamomum camphora (L.) Presl). The scattered points were the observations from data. (AD) represented branch, foliage, stem wood, and stem bark, respectively. The solid line represented the nonlinear regression model (NLR), the dashed line represented the uncorrected linear regression model (LR0) and the dotted line represented the corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 (LRw2).
Figure 2. The curve fitted from nonlinear regression (NLR), uncorrected linear regression (LR0) and corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 ((LRw2) for Cinnamomum camphora (Cinnamomum camphora (L.) Presl). The scattered points were the observations from data. (AD) represented branch, foliage, stem wood, and stem bark, respectively. The solid line represented the nonlinear regression model (NLR), the dashed line represented the uncorrected linear regression model (LR0) and the dotted line represented the corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 (LRw2).
Forests 10 00298 g002
Figure 3. The curve fitted from nonlinear regression (NLR), uncorrected linear regression (LR0) and corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 ((LRw2) for Schima superba (Schima superba Gardn. et Champ.). The scattered points were the observations from data. (AD) represented branch, foliage, stem wood, and stem bark, respectively. The solid line represented the nonlinear regression model (NLR), the dashed line represented the uncorrected linear regression model (LR0) and the dotted line represented the corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 (LRw2).
Figure 3. The curve fitted from nonlinear regression (NLR), uncorrected linear regression (LR0) and corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 ((LRw2) for Schima superba (Schima superba Gardn. et Champ.). The scattered points were the observations from data. (AD) represented branch, foliage, stem wood, and stem bark, respectively. The solid line represented the nonlinear regression model (NLR), the dashed line represented the uncorrected linear regression model (LR0) and the dotted line represented the corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 (LRw2).
Forests 10 00298 g003
Figure 4. The curve fitted from nonlinear regression (NLR), uncorrected linear regression (LR0) and corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 ((LRw2) for Liquidambar formosana (Liquidambar formosana Hance). The scattered points were the observations from data. (AD) represented branch, foliage, stem wood, and stem bark, respectively. The solid line represented the nonlinear regression model (NLR), the dashed line represented the uncorrected linear regression model (LR0) and the dotted line represented the corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 (LRw2).
Figure 4. The curve fitted from nonlinear regression (NLR), uncorrected linear regression (LR0) and corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 ((LRw2) for Liquidambar formosana (Liquidambar formosana Hance). The scattered points were the observations from data. (AD) represented branch, foliage, stem wood, and stem bark, respectively. The solid line represented the nonlinear regression model (NLR), the dashed line represented the uncorrected linear regression model (LR0) and the dotted line represented the corrected linear regression based on the weighted average of secondly basic correction factor from each component equation that is CFw2 (LRw2).
Forests 10 00298 g004
Table 1. The descriptive statistics for 90 sampling trees of each broad-leaved tree species.
Table 1. The descriptive statistics for 90 sampling trees of each broad-leaved tree species.
VariableCinnamomum camphora (L.) PreslSchima superba Gardn. et Champ.Liquidambar formosana Hance.
RangeMeanSDRangeMeanSDRangeMeanSD
Diameter/cm1.9–41.014.510.51.7–51.514.410.91.8–43.514.410.6
Branch/kg0.1–547.044.8105.10.1–371.036.365.90.1–468.730.062.4
Foliage/kg0.1–89.85.812.50.1–45.75.78.70.1–46.44.69.6
Stem wood/kg0.2–519.762.1100.10.3–570.769.2110.50.2–569.981.5127.8
Stem bark/kg0.1–75.510.516.10.1–89.012.619.50.1–121.613.120.5
Aboveground/kg0.4–1016.9123.1221.70.6–897.4123.8192.20.3–955.5129.2206.1
Note: SD was standard deviation.
Table 2. Results of parameter estimate and likelihood analysis based on conditional.
Table 2. Results of parameter estimate and likelihood analysis based on conditional.
SpeciesComponent Δ A I C c ERLRNLR
ab A I C c - logn ab A I C c - norm
Cinnamomum camphoraBranch38.3<<10−20.01231 (0.00218)2.74507 (0.06317)357.80.01213 (0.00197)2.75077 (0.05655)396.1
Foliage263.2<<10−20.01429 (0.00356)2.04953 (0.09561)228.30.01360 (0.00246)2.09921 (0.07636)491.5
Stem wood−55.5>>1020.07100 (0.00741)2.26782 (0.03982) 350.80.08086 (0.02089)2.26358 (0.07766)295.3
Stem bark−15.9>>1020.02016 (0.00307) 2.0898 (0.05842)216.50.03274 (0.01061)1.99888 (0.09682)200.6
Aboveground−105.6>>102396.1290.5
Schima superbaBranch−495.5>>1020.03682 (0.00596)2.37522 (0.06300)340.00.05756 (0.01954)2.18189 (0.10335)−155.5
Foliage−263.2>>1020.07793 (0.01242)1.44971 (0.07602)276.30.04990 (0.01759)1.66296 (0.11131)13.1
Stem wood−530.9>>1020.08755 (0.00892)2.24571 (0.04176)317.80.14024 (0.02612)2.13536 (0.05512)−213.1
Stem bark−311.7>>1020.02301 (0.00295)2.13066 (0.05245)253.00.04092 (0.01291)1.9895 (0.09415)−58.7
Aboveground−617.1>>102364.8−252.3
Liquidambar formosanaBranch−595.5>>1020.01603 (0.00312)2.51909 (0.07487)357.30.01702 (0.00551)2.51466 (0.10100)−238.2
Foliage−221.5>>1020.00385 (0.00117)2.26453 (0.12098)167.80.00475 (0.00219)2.33859 (0.13845)−53.7
Stem wood−670.7>>1020.08077 (0.00991)2.33577 (0.04822)347.80.09299 (0.01878)2.30736 (0.06021)−322.9
Stem bark−401.6>>1020.02048 (0.00310)2.19051 (0.05960)238.50.02329 (0.00480)2.17340 (0.06284)−163.1
Aboveground−731.2>>102379.6−351.6
Note: Values in parentheses are the standard error of the mean. Δ AIC c = AIC c - norm AIC c - logn . AIC c - norm was AIC value calculated from the nonlinear regression model and AIC c - logn was AIC value calculated from the linear regression model. NLR and LR represents the nonlinear and linear regression respectively. ER represents the evidence ratio. The symbol “<<” and “>>” denotes far less than and far greater than.
Table 3. Results of likelihood analysis based on the joint likelihood function for additive system of three tree species.
Table 3. Results of likelihood analysis based on the joint likelihood function for additive system of three tree species.
Species Δ AIC c ERLRNLR
logL LR AIC c - logn logL NLR AIC c - norm
Cinnamomum camphora−7.033−961.81941.3−958.31934.3
Schima superba−810.2>>102−1018.02053.8−612.91243.6
Liquidambar formosana−846.7>>102−952.41922.6−529.11075.9
Table 4. Evaluation statistics of aboveground biomass applying different correction factors for three tree species.
Table 4. Evaluation statistics of aboveground biomass applying different correction factors for three tree species.
SpeciesModelBasic CFAdditive Model CFCF ValueR2SEETREASERMAMPE
Cinnamomum camphoraNLR---0.89773.141.36−5.5217.9210.22
LRCF0-1.000000.88078.919.624.3819.3511.02
CF1 CF a 1 1.139750.90570.25−3.82−8.4218.179.81
CF W 1 1.075440.89872.841.93−2.9417.8310.18
CF2 CF a 2 1.111110.90371.08−1.35−6.0617.829.93
CF W 2 1.098240.90171.63−0.19−4.9617.7610.01
Schima superbaNLR---0.90360.320.92−6.6521.718.54
LRCF0-1.000000.89064.152.615.7221.619.08
CF1 CF a 1 1.126390.84975.06−8.90−6.1519.8610.63
CF W 1 1.065020.87568.42−3.65−0.7420.059.69
CF2 CF a 2 1.082850.86870.11−5.24−2.3719.889.93
CF W 2 1.032480.88465.88−0.622.3920.689.33
Liquidambar formosanaNLR---0.93353.850.23−2.7719.797.31
LRCF0-1.000000.92955.426.545.3021.057.52
CF1 CF a 1 1.176470.91560.36−9.44−10.520.328.19
CF W 1 1.080100.93254.18−1.36−2.5119.677.36
CF2 CF a 2 1.163120.91959.06−8.40−9.4720.128.02
CF W 2 1.068420.93253.99−0.28−1.4519.717.33
Note: CF is correction factor. SEE is standard error of estimate (see Equation (16)). TRE is total relative error (see Equation (17)). ASE is average system error (see Equation (18)). RMA is relatively mean absolute error (see Equation (19)). MPE is mean prediction error (see Equation (20)).
Table 5. Evaluation statistics for each component equation from nonlinear model (NLR), uncorrected linear model (LR0) and linear model corrected by CFw2 (LRw2) for three tree species.
Table 5. Evaluation statistics for each component equation from nonlinear model (NLR), uncorrected linear model (LR0) and linear model corrected by CFw2 (LRw2) for three tree species.
Tree SpeciesComponentModelR2SEETREASERMAMPE
Cinnamomum camphoraBranchNLR0.78251.651.63−2.3941.7019.18
LR00.78051.872.09−2.5641.5819.26
LRw20.79949.59−7.05−11.2740.6518.42
FoliageNLR0.5808.17−3.93−5.7846.8524.91
LR00.5548.427.430.8147.7325.68
LRw20.5728.25−2.18−8.2146.3125.16
Stem woodNLR0.87535.621.74−3.1022.4710.07
LR00.85438.4914.269.2224.6510.88
LRw20.87335.884.04−0.5522.5710.14
Stem barkNLR0.8486.311.01−10.5430.0310.52
LR00.8067.1222.1616.0335.5511.88
LRw20.8366.5411.235.6530.3210.91
Schima superbaBranchNLR0.63639.957.21−5.2536.4219.31
LR00.53345.25−12.14−6.6634.2421.87
LRw20.50846.44−14.9−9.6033.9322.45
FoliageNLR0.6555.165.160.5051.0715.84
LR00.5605.8231.528.7756.9817.86
LRw20.5745.7327.385.3555.9517.59
Stem woodNLR0.90833.73−2.07−7.7925.348.55
LR00.90534.278.7613.0327.098.68
LRw20.90634.005.349.4725.698.61
Stem barkNLR0.8138.49−1.04−10.2828.6011.79
LR00.8048.6810.6713.0830.5412.05
LRw20.8058.677.199.5229.1212.04
Liquidambar formosanaBranchNLR0.60939.220.87−2.3737.8922.94
LR00.60739.315.522.5838.8422.99
LRw20.60839.24−1.23−3.9937.6822.95
FoliageNLR0.5976.14−0.397.0368.2723.45
LR00.4687.0556.9756.2594.1326.93
LRw20.4946.8846.9246.2588.2526.28
Stem woodNLR0.93233.40.05−3.4121.597.19
LR00.93133.814.843.8323.207.28
LRw20.93133.86−1.87−2.8221.647.29
Stem barkNLR0.9096.200.17−1.5327.978.28
LR00.9026.467.697.4630.348.63
LRw20.9106.190.790.5728.468.27

Share and Cite

MDPI and ACS Style

Cao, L.; Li, H. Analysis of Error Structure for Additive Biomass Equations on the Use of Multivariate Likelihood Function. Forests 2019, 10, 298. https://doi.org/10.3390/f10040298

AMA Style

Cao L, Li H. Analysis of Error Structure for Additive Biomass Equations on the Use of Multivariate Likelihood Function. Forests. 2019; 10(4):298. https://doi.org/10.3390/f10040298

Chicago/Turabian Style

Cao, Lei, and Haikui Li. 2019. "Analysis of Error Structure for Additive Biomass Equations on the Use of Multivariate Likelihood Function" Forests 10, no. 4: 298. https://doi.org/10.3390/f10040298

APA Style

Cao, L., & Li, H. (2019). Analysis of Error Structure for Additive Biomass Equations on the Use of Multivariate Likelihood Function. Forests, 10(4), 298. https://doi.org/10.3390/f10040298

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