Parsimonious and Completely Compatible Taper ... - Semantic Scholar

6 downloads 0 Views 105KB Size Report
Jul 22, 2005 - Acknowledgments: The authors would like to thank Alexander Clark III, Wood Scientist, of the USDA Forest Service Southern Research Station ...
Parsimonious and Completely Compatible Taper, Total, and Merchantable Volume Models Lewis Jordan, Kenneth Berenhaut, Ray Souter, Richard F. Daniels Abstract: A parsimonious and compatible system of taper, total, and merchantable volume equations sharing a common vector of three parameter estimates was derived from Clark and Saucier’s (1990) volume ratio model. Clark and Saucier’s model expresses volume inside bark as a function of outside bark diameter and is itself not compatible with the derived taper equation. To ensure compatibility, models were fit to express volume outside bark as functions of diameter outside bark. By constraining the limits of integration of the taper function to go from the tip to groundline diameter, a completely compatible system of taper, total, and merchantable volume equations were developed. Systems of equations were fit using nonlinear seemingly unrelated regressions with nonlinear cross-equation constraints to account for contemporaneous correlations in the data. The system fit was found to account for at least 96% of the variation in diameter outside bark, total, and merchantable outside bark volume. FOR. SCI. 51(6):578 –584. Key Words: Compatible equations, merchantable volume, seemingly unrelated regression, taper.

T

HE DEVELOPMENT OF STEM PROFILE TAPER and mer-

chantable volume models are common in forestry related literature. Numerous taper and subsequently derived merchantable volume functions of various forms have been developed especially in the last 30 years, from simple taper models (Kozak et al. 1969, Ormerod 1973, Hilt 1980), segmented polynomial models (Bruce et al. 1968, Max and Burkhart 1976, Cao et al. 1980), and geometryoriented models (Parresol and Thomas 1996, Fang and Bailey 1999, Zhang et al. 2002). Newnham (1988) states that the two reasons for continuous study in this area are: 1) because no single theory has been developed that adequately explains the variation in stem form for all kinds of trees, and 2) as a method of estimating volume, a single taper equation can estimate both total and merchantable tree stem volume. If tree form can be accurately described, then volume for any merchantability limit can be accurately predicted. Merchantable volume equations allow for the prediction of volume to a specified upper stem diameter. For utilization purposes, it is desirable to merchandize trees into multiple products, hence the development of a taper function is essential (McTague and Bailey 1987). Clutter (1980) showed that, for a given merchantable volume equation, there is an intrinsically defined compatible taper function. This idea implies that integration of the taper function from the ground to total height of the tree

would return the appropriate volume, and subsequently the merchantable volume equation. Clutter (1980) states “The accuracy and precision of such taper functions will depend on the accuracy and precision of the merchantable volume equations from which they are derived.” Conversely, the same holds true for merchantable volume equations derived from taper functions. Precision and accuracy are a direct result of the model forms, be they empirical or geometry-oriented taper and volume functions, that best fit a unique data set. Burkhart (1977) developed a merchantable volume ratio equation based on upper stem diameter, which eliminated the need for separate volume equations for differing stem diameters. Cao and Burkhart (1980) developed a volume ratio model based on the distance from the tree top to a merchantable height. Knoebel et al. (1984) developed a technique to derive compatible taper and merchantable height functions from the Burkhart (1977) and Cao and Burkhart (1980) volume ratio models. However, McTague and Bailey (1987) showed that the integration of the Knoebel et al. (1984) taper function from the groundline to some merchantable height does not equal the value of merchantable volume from the Cao and Burkhart (1980) ratio model. Van Deusen et al. (1981), Clark and Saucier (1990), and Tasissa et al. (1997) presented exponential volume ratio models for the prediction of stem volume to some

Lewis Jordan, Research Coordinator, Warnell School of Forest Resources, The University of Georgia, Athens, GA 30602-2152—Phone: (706) 542-7587; Fax: (706) 542-0119; [email protected]. Kenneth Berenhaut, Assistant Professor of Mathematics, Wake Forest University, 1834 Wake Forest Road, Winston-Salem, NC 27106 —Phone: (336) 758-5922; [email protected]. Ray Souter, Principal Analyst, Rork Associates, P.O. Box 6412, Athens, GA 30604 —Phone: (706) 369-5341; [email protected]. Richard F. Daniels, Professor, Warnell School of Forest Resources, The University of Georgia, Athens, GA 30602-2152—Phone: (706) 542-7268; [email protected]. Acknowledgments: The authors would like to thank Alexander Clark III, Wood Scientist, of the USDA Forest Service Southern Research Station for the data used herein. We would also like to thank the editor, associate editor, and the three anonymous reviewers at Forest Science for their helpful and insightful comments. Personally, I would like to thank the co-authors of this article, for their enthusiasm, support, and especially their never ending patience. Manuscript received December 18, 2004, accepted July 22, 2005

578

Forest Science 51(6) 2005

Copyright © 2005 by the Society of American Foresters

merchantable diameter limit. The exponential ratio models possess desirable properties such that, as diameter outside bark approaches infinity, the ratio goes to zero, thus predicted that volume goes to zero. Secondly, when diameter outside bark is equal to zero, the ratio becomes one, thus volume is equal to total stem volume. Clark and Saucier (1990) introduced a widely used and accepted merchantable tree volume exponential ratio equation for plantation-grown loblolly pine (Pinus taeda L.) in the southeastern United States. Clark and Saucier’s (1990) exponential ratio model was developed to predict merchantable volume inside bark as a function of diameter outside bark. In itself, Clark and Saucier’s (1990) model can not be differentiated to a compatible taper equation, because volume inside bark is calculated as a function of outside bark diameter. The purpose of this article is derivation of a compatible system of taper and merchantable volume models from the popular exponential ratio model presented by Clark and Saucier (1990). To ensure compatibility, models were fit to express volume outside bark as a function of diameter outside bark.

D ⫽ f共T兲 .

Differentiating Equation 5 with respect to T gives kf共T兲 ⫽ ⫺VT



f共t兲 dt ⫽ VT ⫺ VT exp共␤1 Dm␤ 2D␤ 3兲,

␤1 D␤ 3␤2 关 f共T兲兴␤ 2/ 2⫺1 2

䡠 exp共␤1 关 f共T兲兴␤ 2/ 2 D␤ 3兲

df共T兲 , dT

(6)

or, following some rearrangement, k ⫽ ⫺VT

df共T兲 ␤1 D␤ 3␤2 关 f共T兲兴␤ 2/ 2⫺2 䡠 exp共␤1 关 f共T兲兴␤ 2/ 2 D␤ 3兲 . 2 dT (7)

Equation 7 is now a separable differential equation involving T and f(T). The integration of Equation 7 gives T

␤1 D␤ 3␤2 k dT ⫽ ⫺VT 2



2 Dm

关 f共T兲兴 ␤ 2/ 2⫺2

0

䡠 exp共 ␤ 1 关 f共T兲兴 ␤ 2/ 2 D ␤ 3兲 df共T兲.

(8)

From Equation 8, using substitution, let x ⫽ f共T兲 ␤ 2/ 2 ,

x

2/ ␤ 2

dx ⫽

␤2 f共T兲␤ 2/ 2⫺1 d关 f共T兲兴, 2 (9)

2 2/ ␤ 2⫺1 df共T兲 ⫽ x dx. ␤2

⫽ f共T兲,

From Equation 9, D2m ⫽ f(T), thus x ⫽ D␤m2. Also, when f(T) ⫽ 0, x ⫽ 0, hence the integration constant equals zero. Substituting the results from Equation 9 into Equation 8 gives kT ⫽ ⫺VT

␤1 D␤ 3␤2 2





D m2

x1⫺4/ ␤ 2exp共␤1 xD␤ 3兲

0

␤1 D␤ 3␤2 2 ⫽ ⫺VT 2 ␤2



2 2/ ␤ 2⫺1 x dx ␤2 共10兲



D m2

x⫺2/ ␤ 2exp共␤1 xD␤ 3兲 dx.

0

(3) From Equation 10, let

It follows from Clutter (1980) that, T

(5)

0

(2)

Let, T ⫽ the distance from the top of the stem to the merchantability limit Dm, and M ⫽ the merchantable height (distance from the ground to the height of the merchantability limit Dm), so that H ⫽ T ⫹ M. By definition of the taper function, 2 m

f共t兲 dt ⫽ VT ⫺ VT exp共␤1 关 f共T兲兴␤ 2/ 2 D␤ 3兲 .

(1)

where, VDIB ⫽ merchantable inside bark volume (m3) to an outside bark top diameter (Dm) (cm); VT ⫽ ␾1D2␾2H␾3 ⫽ total inside bark stem volume, where D is dbh (cm), and H is total height (m); VR ⫽ exp(␤1D␤m2D␤3) ⫽ the volume ratio, where e is the base of the natural logarithm, and ␾ ⫽ (␾1 ␾2 ␾3)T and ␤ ⫽ (␤1 ␤2 ␤3)T are coefficients to be estimated. The objective is to express outside bark diameter (do) (cm) as a function of distance from the top of the stem t. Thus, the taper function to be derived is given as d 2o ⫽ f共t兲 .

T

0

The general form of the model is given as V DIB ⫽ VTVR ,



k



Derivation

k

3 into Equation 4 gives

z ⫽ ⫺␤1 xD␤ 3, (4)

0

where k ⫽ ␲/40000, for conventional metric measurement units. Substituting for D2m from the relationship in Equation

z , x⫽ ⫺␤1 D␤ 3

dz ⫽ ⫺␤1 D␤ 3 dx, dz ⫽ dx. ⫺␤1 D␤ 3

(11)

From Equation 11, x ⫽ D␤m2 thus z ⫽ ⫺␤1D␤m2D␤3, and when x ⫽ 0, z ⫽ 0, hence the integration constant equals zero. Forest Science 51(6) 2005

579

Table 1. Average, standard deviation and range of total height, dbh, and age of the 126 loblolly pine trees sampled for analysis.

Variable

Mean

Standard deviation

Min

Max

Total Height (m) dbh (cm) Age (Years)

20.65 23.81 22.63

2.54 4.54 1.77

15.0 15.49 20

25.52 36.57 25

Substitution of the results from Equation 11 into Equation 10 gives

␤1 D␤ 3 kT ⫽ ⫺VT ⫺␤1 D␤ 3





⫺ ␤ 1 D m2D ␤ 3

0

⫽ V T共⫺␤1 D␤ 3兲2/ ␤ 2

冕 冕

z⫺2/ ␤ 2 e⫺z dz, 共⫺␤1 D␤ 3兲⫺2/ ␤ 2

V C ⫽ VT ⫽

kH

冉 冊



冉 冊冊

␤2 ⫺ 2 ␤2 ⫺ 2 ⌫ 共⫺␤1 D␤ 3兲2/ ␤ 2 ⌫ ⫺␤1 Dm␤ 2D␤ 3, ␤2 ␤2

. 共15兲

In our case, we need to constrain Equation 15 to go from the tip to groundline diameter (Do). The compatible total volume equation now becomes a function of total height, dbh, and groundline diameter. Substituting Do for Dm in Equation 15, the compatible volume equation takes the form VC ⫽



⫺ ␤ 1 D m2D ␤ 3

simpler solution is to solve Equation 13 for VT. Solving Equation 13 for VT and substituting H for T gives

z⫺2/ ␤ 2 e⫺z dz,

kH

冉 冊



冉 冊冊

␤2 ⫺ 2 ␤2 ⫺ 2 ⌫ 共⫺␤1 D␤ 3兲2/ ␤ 2 ⌫ ⫺␤1 D␤ 3D␤o 2, ␤2 ␤2

0

⫽ V T共⫺␤1 D␤ 3兲2/ ␤ 2

z共 ␤ 2⫺2 兲 / ␤ 2⫺1 e⫺z dz,

0

⫽⌫



(16)

(12)

␤ ⫺ ␤ 1 D m2D ␤ 3





冉 冊冊

␤2 ⫺ 2 ␤2 ⫺ 2 V T共⫺␤1 D␤ 3兲2/ ␤ 2 ⌫ ⫺␤1 Dm␤ 2D␤ 3, . ␤2 ␤2

Equation 16 is a better total volume equation, or at least completely compatible with the taper model from ground to tip. Equation 16 is also equal to the total volume equation (VT) presented in Equation 1. Let, V C ⫽ VT ⫽ ␾1 D2 ␾ 2H␾3,

Rearrangement of Equation 12 gives kT ⌫共关 ␤ 2 ⫺ 2兴/ ␤ 2 兲V T共⫺␤1 D␤ 3兲2/ ␤ 2



冉 冊冊

Solving Equation 13 for Dm and substituting the relationship (H ⫺ M) for T gives







冉 冊冊

k共H ⫺ M兲 ␤2 ⫺ 2 , ⌫共关␤2 ⫺ 2兴/␤2 兲VT共⫺␤1 D␤ 3兲2/ ␤ 2 ␤2

1 ⫺␤1 D␤ 3



冋冉

␾1 ⫽ k ⌫

␾2 ⫽ ⫺

D m ⫽ ⌫⫺1

(14)

Table 2. Parameter estimates, standard errors and P-values from the COMP (Equations 17, 18, 19) and ground line diameter (Equation 20) system fit.

Parameter

Estimate

Standard Error

P-value

␤1 ␤2 ␤3 ␣1 ␣2 ␣3

⫺0.9994 4.0419 ⫺3.9200 ⫺4.5683 1.3274 0.3167

0.0444 0.0076 0.0151 0.4948 0.0175 0.0301

0.0001 0.0001 0.0001 0.0001 0.0001 0.0001

Forest Science 51(6) 2005



␤3 ␤2



⫺1

,

and ␾3 ⫽ 1.

Substitution of VC into the taper model gives

1/ ␤ 2

⫹ ␧.

冊冉

␤2 ⫺ 2 ␤2 ⫺ 2 ⌫ ⫺␤1 D␤ 3D␤o 2, 共⫺␤1 兲2/ ␤ 2 ␤2 ␤2

1/ ␤ 2

A compatible total volume (m3) outside bark equation, denoted as VC, can be readily obtained by integrating the taper function (Equation 14). However, an equivalent and

580

(17)

where,

␤2 ⫺ 2 ⫽ ⌫ ⫺␤1 Dm␤ 2D␤ 3, . (13) ␤2

D m ⫽ ⌫⫺1

.





冢冉

冉 冊冣

k共H ⫺ M兲 ␤2 ⫺ 2 , ␤2 ␤2 ⫺ 2 ⌫ VC共⫺␤1 D␤ 3兲2/ ␤ 2 ␤2

1 ⫺␤1 D␤ 3





1/ ␤ 2

1/ ␤ 2

⫹ ␧.

(18)

The taper Equation 18 is now completely compatible with the total volume Equation 17. Subsequently, the volume ratio (VR) presented in Equation 1 is still applicable. Thus,

Table 3. Fit statistics of the COMP (Equations 17, 18, 19) and ground line diameter (Equation 20) system fit, where R2 is the coefficient of determination, RMSE is the root mean square error, MB is the mean bias, and MAB is the absolute mean bias

Equation

R2

RMSE

MB

MAB

Dm (cm) VC (m3) VDOB (m3) Do (cm)

0.9647 0.9703 0.9871 0.7972

1.43752 0.0366 0.0256 3.2796

⫺0.2545 ⫺0.0037 ⫺0.0017 ⫺0.1301

1.0946 0.0279 0.0172 2.4867

Table 4. Cross-equation correlation matrix of residuals for the COMP (Equations 17, 18, 19) and ground line diameter (Equation 20) system fit.*

Dm Vc VDOB Do

Dm

Vc

VDOB

Do

1

0.4355 1

0.3100 0.7879 1

0.2142 0.4149 0.4244 1

* All correlations were found to be significant at the 0.0001 level.

merchantable outside bark volume can now be defined as V DOB ⫽ VC exp共␤1 Dm␤ 2D␤ 3兲 ⫹ ␧,

(19)

where, VDOB ⫽ compatible merchantable outside bark volume, and all other variables previously defined. The total volume Equation 17, the taper Equation 18, and the merchantable outside bark volume Equation 19 (call them the COMP equations) are now completely compatible, sharing a common vector of parameter estimates, ␤ ⫽ (␤1 ␤2 ␤3)T, and specified predictor variables. Because the COMP equations rely on groundline diameter, a prediction equation was developed to relate groundline diameter (Do) to dbh and total height (H), further enhancing the utility of this compatible system. The equation takes the form D o ⫽ ␣1 ⫹ ␣2 dbh ⫹ ␣3 H ⫹ ␧.

Analysis Seemingly Unrelated Regression Development of compatible taper and volume functions involves parameter estimation of the taper function, which is integrated to provide volume. This approach will minimize the error associated with diameter estimation, but does not ensure minimal error in volume estimation. A set of nonlinear equations that has contemporaneous cross-equation error correlation is known as nonlinear seemingly unrelated regression (NSUR) system. A system of this type, which contains statistical dependencies, is best fit using nonlinear joint-generalized least squares regression. At first look, the taper, total, and merchantable volume equations seem unrelated, but the equations are related through the correlation in the errors. A set of nonlinear regression functions are specified such that 1) each equation can contain its own independent variables, 2) each equation can use its own weight function (if needed), and 3) conformity (e.g., a volume function derived from a taper function) is ensured by setting constraints on the regression coefficients (i.e., parameter sharing). The structural equations for the system of nonlinear models can be specified as: y1 ⫽ f1 共X1 , ␤1 兲 ⫹ ␧1 y2 ⫽ f2 共X2 , ␤2 兲 ⫹ ␧2 ⯗ yk ⫽ fk 共Xk , ␤k 兲 ⫹ ␧k ,

(21)

(20)

Data The suitability of the derived equations were examined with data from 126 plantation-grown loblolly pine trees representing 42 stands from the Southern Atlantic Coastal Plain of the United States. The stands were located on land owned by forest products companies, and included only stands with similar silvicultural history: 1) site preparation with no herbaceous weed control; 2) no fertilization at planting except phosphorous on phosphorous-deficient sites; 3) stand density of at least 617 trees per hectare at the time of sampling. Outside bark diameter was measured with a diameter tape at the groundline, 0.30-, 0.61-, and 1.37-m height. The trees were then felled and measured for stump height and stump diameter, total height, and height to a 17.7-, 15.2-, 10.1-, 7.6-, and 5.0-cm outside bark diameter. Cross-sectional disks 2.54 cm thick were then cut for use in another study at 1.37 m and then at 1.5-m intervals to 5.0-cm diameter outside bark. Disk diameter outside bark was then measured with a diameter tape at the USDA Southern Research Station wood properties lab. Volumes (m3) were determined with Smalian’s equation applied to successive diameter pairs, except the topmost bolt was assumed to be a cone and its volume was calculated accordingly. Tree size characteristics are summarized in Table 1.

where yi is a vector containing the dependent variable from the ith equation, Xi is a matrix containing the independent variables from the ith equation, ␤i is the parameter vector for the ith equation, and ␧i is the random error vector for the ith equation. When the stochastic properties of the error vectors are specified, along with the coefficient restrictions, the structural equations become a statistical model for efficient parameter estimates and reliable prediction intervals. The use of NSUR offers the best unbiased estimator for ␤i, which has a lower variance than the estimator of the ordinary nonlinear least squares estimator of ␤i because it takes into account contemporaneous correlation in the different equations. The availability of econometric software, such as SAS/ETS (SAS Institute, Inc., 1999, v8doc.sas.com/ sashtml/) makes complicated statistical procedures like NSUR easily implemented. It would be unrealistic to expect that the equation errors would be uncorrelated (Borders 1989, Parresol 1999, 2001). In an attempt to simultaneously minimize the error associated with these equations, the models in this article were fit as NSUR using the SAS/ETS Model Procedure (SAS Institute, Inc.).

Comparison Criteria The comparison of the models was based on four statistical indices: coefficient of determination (R2), root mean square error (RMSE), mean bias (MB), and mean absolute Forest Science 51(6) 2005

581

bias (MAB) (Loague and Green 1991, Mayer and Butler 1993). These criteria are given as 共Y ⫺ Yˆ 兲 冘 R ⫽ 1⫺ 冘 共Y ⫺ Y៮ 兲 , 共Y ⫺ Yˆ 兲 冘 RMSE ⫽ 冑 , 2

n i⫽1

i

n i⫽1

i

n i⫽1

2

i

2

(22)

Results

2

i

i

n⫺p



(23)

n

MB ⫽

1 共Y ⫺ Yˆi 兲, n i⫽1 i

(24)

MAB ⫽

1 Y ⫺ Yˆi , n i⫽1 i



(25)

冘冏 n

models should be refit to obtain more appropriate parameter estimates. The models fit in this study were evaluated using the goodness-of-fit statistics and graphical analyses described above.

where, Yi, Yˆi, and Y¯i are the actual, predicted, and average values of the dependent variable, n is the total number of observations used to fit the model, and p is the number of parameters in the model. The models evaluated in this article were also examined visually by plotting the fitted taper and volume ratio equations. To truly assess performance of the models fit in this study, validation of the models with an independent data set would be the most desirable approach, because quality of fit does not always ensure the quality of prediction (Huang et al. 2003, Kozak and Kozak 2003). The commonly used validation methods of data splitting and cross-validation, as shown by Kozak and Kozak (2003), do not provide any additional information on model performance compared to the statistics obtained from models fit to the entire data set. Models validated with an independent data set prove that either the data are from the same population and will perform as per se validation utilizing data splitting, or the data are from a different population entirely, in which case the

Initial parameter estimates for the COMP (Equations 17, 18, 19) models can be obtained by simultaneously fitting the outside bark total and merchantable volume equations (Equation 1). The COMP and groundline diameter prediction models were fit to the complete data set. Convergence was acquired after two iterations. All parameters in the model were found to be significant at the 0.0001 level (Table 2). As seen in Table 2, all parameter estimates are logical and ensure compatibility between equations. Fit statistics from the criterion described above are presented in Table 3. All of the COMP models were found to account for at least 96% of the variation in diameter outside bark, total, and merchantable volume. The mean bias values for all of the models were found to be negative, indicating that the models are slightly over predicting, but are generally small, suggesting that the models are overall unbiased. In addition to negative mean bias values for the COMP fit, positive correlation was found to exist among the error components (Table 4). This is intuitively reasonable because, for example, if diameter outside bark is overestimated, we would expect volume to be overestimated. The taper function defined by Equation 18 was plotted with various dbh, total tree height, and corresponding predicted groundline diameter values (Figure 1). When height above ground is equal to zero, outside bark diameter is equal to the predicted groundline diameter. Similarly, when height above ground is equal to the total height of the tree, outside bark diameter is equal to zero. The volume ratio, defined as

Figure 1. Outside bark taper curve from the COMP system fit with predicted groundline diameter at varying dbh (cm) and total height (m) values.

582

Forest Science 51(6) 2005

Figure 2. Volume ratios from the COMP system fit with predicted groundline diameter at varying dbh (cm) and total height (m) values.

cumulative volume divided by total volume, is plotted in Figure 2 with various dbh, total tree height, and corresponding predicted groundline diameter values. When height above ground is equal to total height and zero, the volume ratios are equal to one and zero, respectively.

Discussion A logical, consistent, and parsimonious system of compatible taper, total, and merchantable volume outside bark equations were developed from Clark and Saucier’s (1990) volume ratio models. The equations described herein were fit simultaneously to account for cross-equation correlation among the error components. We apply this taper model to the tree stem from the top down and constrain it to predict the same total stem volume when integrated as a direct volume prediction method for the total stem. This constraint was imposed by defining the limits of integration of the taper and total volume equations. The imposed restriction ensures diameter outside bark is equal to groundline diameter and zero when height above ground is equal to zero and total tree height, respectively. Similarly, the merchantable volume equation will yield volume ratio values of zero and one when height above ground is equal to zero and total height, respectively, ensuring compatibility between the merchantable and total volume equations. The models presented herein were found to account for a significant amount of the variation in taper, total, and merchantable volume. Here, we reparameterize our total volume equation to share the same parameters as the volume ratio. McTague and Bailey (1987) state that there is no reason why the parameters of total volume and merchantable volume equations should be estimated separately, stating that total volume may be regarded as a special case of merchantable volume when the upper stem diameter is zero.

Even though McTague and Bailey (1987) simultaneously estimated the parameters of a total and merchantable volume equation, the parameter estimates for the total and merchantable parts of the model differed. The simultaneous system of equations presented in this study is unique in that the taper, total, and merchantable volume equations are composed of a single common vector of three parameter estimates, ensuring parsimony and complete compatibility. As pointed out by one reviewer, the data used in this study were collected from trees in stands with similar growing conditions, ages, and size characteristics, from within a narrow geographical range, severely restricting the use of the model. Thus, caution should be used when extrapolating beyond the natural range of the data on which the model is based. However, the original purpose of this study was derivation of a completely compatible system of equations. The models presented in this article provide alternatives from which data collected in other regions, or for a larger more inclusive data set could be based.

Literature Cited BORDERS, B.E. 1989. Systems of equations in forest stand modeling. For. Sci. 35:548 –556. BRUCE, D., R.O. CURTIS, AND C. VANCOEVERING. 1968. Development of system and taper volume tables for red alder. For. Sci. 14:339 –350. BURKHART, H.E. 1977. Cubic foot volume of loblolly pine to any merchantable top diameter. South. J. Appl. For. 1:7–9. CAO, Q.V., AND H.E. BURKHART. 1980. Cubic-foot volume of loblolly pine to any height limit. South. J. Appl. For. 4:166 –168. CAO, Q.V., H.E. BURKHART,

AND

T.A. MAX. 1980. Evaluation of Forest Science 51(6) 2005

583

two methods for cubic-volume prediction of loblolly pine to any merchantable limit. For. Sci. 26:71– 80.

MAX, T.A., AND H.E. BURKHART. 1976. Segmented polynomial regression applied to taper equations. For. Sci. 22:283–289.

CLARK, A., III, AND J.R. SAUCIER. 1990. Tables for estimating total-tree weights, stem weights, and volumes of planted and natural southern pines in the southeast. Macon, GA Georgia Forestry Commission. Georgia For. Res. Paper 79. 23 p.

MAYER, D.G., AND D.G. BUTLER. 1993. Statistical validation. Ecol. Model. 68:21–32.

CLUTTER, J.L. 1980. Development of taper functions from variable-top merchantable volume equations. For. Sci. 26:117–120.

MCTAGUE, J.P., AND R.L. BAILEY. 1987. Simultaneous total and merchantable volume equations and a compatible taper function for loblolly pine. Can. J. For. Res. 17:87–92. NEWNHAM, R.M. 1988. A variable-form taper function. Petawawa Nat. For. Inst. Information Rep. PI-X-83, 33 p.

FANG, Z., AND R.L. BAILEY. 1999. Compatible volume and taper models with coefficients for tropical species on Hainan Island in southern China. For. Sci. 45:85–100.

ORMEROD, D.W. 1973. A simple bole model. For. Chron. 49:136 –138.

HILT, D.E. 1980. Taper-based system for estimating stem volume of upland oaks. USDA For. Serv. Res. Pap. NE-458. 12 p.

PARRESOL, B.R. 1999. Assessing tree and stand biomass: A review with examples and critical comparisons. For. Sci. 45:573–593.

HUANG, S., Y. YANG, AND Y. WANG. 2003. A critical look at procedures for validation growth and yield models. P. 271–294 in Modeling forest systems. Amaro, A, D.D. Reed, and P. Soares, (Eds.). CABI Publishing, Wallingford, UK.

PARRESOL, B.R. 2001. Additivity of nonlinear biomass equations. Can. J. For. Res. 31:865– 878.

KNOEBEL, B.R., H.E. BURKHART, AND D.E. BECK. 1984. Stem volume and taper function for yellow-poplar in the southern Appalachians. South. J. Appl. For. 8:185–188. KOZAK, A., D.O. MUNRO, AND J.H.G. SMITH. 1969. Taper functions and their application in forest inventory. For. Chron. 45:278 –283. KOZAK, A., AND R. KOZAK. 2003. Does cross validation provide additional information in evaluation of regression models? Can. J. For. Res. 33:976 –987. LOAGUE, K., AND R.E. GREEN. 1991. Statistical and graphical methods for evaluating solute transport models: Overview and application. J. Contam. Hydrol. 7:51–73.

584

Forest Science 51(6) 2005

PARRESOL, B.R., AND C.C. THOMAS. 1996. A simultaneous density-integral system for estimating stem profile and biomass: Slash pine and Willow Oak. Can. J. For. Res. 26:773–781. TASISSA, G., H.E. BURKHART, AND R.L. AMATEIS. 1997. Volume and taper equations for thinned and unthinned loblolly pine trees in cutover, site-prepared plantations. South. J. Appl. For. 21:146 –152. VAN DEUSEN, P.C., A.D. SULLIVAN, AND T.G. MATNEY. 1981. A prediction system for cubic foot volume of loblolly pine applicable through much of its range. South. J. Appl. For. 21:186 –189. ZHANG, Y., B.E. BORDERS, AND R.L. BAILEY. 2002. Derivation, fitting and implication of a compatible stem taper-volumeweight system for intensively managed, fast growing loblolly pine. For. Sci. 48:595– 607.