Next Article in Journal
Effects of Soil Microorganisms on Carbon Sequestration under Different Mixed Modification Models in Pinus massoniana L. Plantation
Next Article in Special Issue
Genetic Parameters Estimated in the Early Growth of Dimorphandra mollis Benth. Progenies
Previous Article in Journal
Woody Plant Structural Diversity Changes across an Inverse Elevation-Dependent Warming Gradient in a Subtropical Mountain Forest
Previous Article in Special Issue
A Comparison of Probability Density Functions Fitted by Moments and Maximum Likelihood Estimation Methods Used for Diameter Distribution Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Different Modeling Approaches for Estimating Total Bole Volume of Hispaniolan Pine (Pinus occidentalis Swartz) in Different Ecological Zones

by
Santiago W. Bueno-López
1,*,
Luis R. Caraballo-Rojas
1 and
Juan G. Torres-Herrera
2
1
Vice-Rectory of Research and Innovation, Pontificia Universidad Católica Madre y Maestra (PUCMM), Santiago de los Caballeros 51000, Dominican Republic
2
Consorcio Forestal Serrano, San Jose de las Matas 51000, Dominican Republic
*
Author to whom correspondence should be addressed.
Forests 2024, 15(6), 1052; https://doi.org/10.3390/f15061052
Submission received: 30 April 2024 / Revised: 11 June 2024 / Accepted: 12 June 2024 / Published: 18 June 2024
(This article belongs to the Special Issue Forest Biometrics, Inventory, and Modelling of Growth and Yield)

Abstract

:
Pinus occidentalis (Swartz) is the primary timber species in the Dominican Republic (DR). Despite its economic importance, studies conducted on this species are scarce, making it difficult to estimate current inventory levels. This study aims to enhance the accuracy of estimating the total bole volume of P. occidentalis in different ecological zones (EZs) within La Sierra, evaluating and comparing two established volume equations—combined variable (CV) and Schumacher and Hall (S&H) across nine modeling variants. An indicator variables analysis determined the necessity of distinct equations for two EZs. Fitting included both linear and nonlinear models. Our comprehensive statistical analysis included goodness-of-fit metrics to evaluate each model variant’s performance rigorously. The second modeling variant (SH02) for the SH equation was most effective in the Dry Ecological Zone, showing superior performance in both the fitting and validation phases. Similarly, the third modeling variant (SH03) for the SH equation emerged as the best fit for the Combined Intermediate and Humid Ecological Zones, achieving the lowest overall ranking sum among tested variants. SH02 and SH03 provide reliable and precise volume estimations, allowing for the optimization of forestry management practices for P. occidentalis trees. The SH models outperformed the CV model variants’ consistency in parameter estimation. This tailored approach ensures more accurate volume predictions, which is crucial for sustainable management and conservation efforts.

1. Introduction

Accurate measurements of tree dimensions like diameter, height, stem form, and volume are not just academic exercises but crucial for estimating forest attributes. This research provides practical tools like individual-tree and species–specific volume equations, which are necessary to predict inventory levels and ensure long-term wood yields. Volume equations use diameter at breast height (DBH) and height (H) to calculate stem volume, providing reliable estimates of tree volumes by accounting for size and shape variations [1]. These practical tools are essential in forestry practices, empowering forest managers and timber industry professionals with accurate and reliable data.
Stem volume models estimate timber volume, biomass, and carbon sequestration potential for individual trees and larger areas. They are useful in forest management planning and growth simulations, driven by the need for the accurate estimation of tree bole volume. Because the same regression may not be equally suitable for predicting total volume estimates in different ecological conditions [1], specific details and techniques for developing statistical models for stem volume estimation may vary.
The history of excurrent tree bole volume estimation dates back several decades [2]. Various statistical techniques can be used to develop stem volume models, including linear regression, nonlinear regression, mixed-effects procedures, and machine learning algorithms [2]. The choice depends on the data set’s underlying assumptions and characteristics [2]. They are developed in a two-step process: model fitting and validation. The model is fitted to the data set using a statistical technique and then validated using an independent data set. Evaluation uses measures such as mean squared error, root mean squared error, bias, and the coefficient of determination.
The most common statistical procedure for developing volume tables are ordinary least squares (OLS) regression analyses, which relate bole volume to explanatory variables such as DBH, H, and sometimes stem form. When dealing with biological data, the constant variance assumption is often violated whenever a direct measure of stem content is used as the dependent variable in a regression equation [3]. If the assumption of constant variance is unmet, the equation must be weighted by a factor proportional to the standard deviation of the dependent variable.
Weighted least squares can be used when the ordinary least squares assumption of constant variance in the errors is violated. It will produce a new regression model which results in the dependent variable having constant variance [3]. In the weighted linear regression model, each observation is assigned a weight W i . The weighted sum of squared residuals, Q = i = 1 n W i × ε ^ i 2 is minimized.
Assumptions related to linear and nonlinear statistical models involve linearity, independence, homoscedasticity, normality, and zero mean of residuals [4]. Additionally, nonlinear models must have a correct expectation function [4]. These assumptions guide the development and application of these equations in forestry practices, and their violation may lead to biased predictions and incorrect inferences.
The scarcity of data on volume models for tropical and sub-tropical tree species negatively impacts the accuracy of estimating tree volume for both research and operational purposes [5]. Even when developed, stem volume models should be periodically updated to account for changes in forest structure, climate conditions, or management practices to ensure accurate predictions and relevant decision support. P. occidentalis is the main timber species in the DR, growing on approximately 302,500 hectares and comprising approximately 95% of all the timber harvested [6]. Despite its economic importance, estimating inventory levels and accounting for harvested volume is difficult because no standardized system is used for volume appraisal and inventory purposes.
The Schumacher and Hall [7] equation (SH) has been widely used in forestry to estimate tree bole volumes. Other researchers have shown that this equation and the Combined Variable equation developed by Bennett et al. [8] provide accurate results for estimating tree volumes, good performance in graphical analysis, and reliable estimates without bias [9].
The objectives of this study were as follows: (1) to fit two commonly used volume equations, the combined variable (CV) and the SH equations with nine different modeling alternatives, to estimate the total bole volume content of P. occidentalis trees growing in three ecological zones (EZ) within La Sierra, DR; (2) to evaluate goodness-of-fit statistics of these nine modeling alternatives in each EZ; (3) to conduct and indicator variable analysis to determine if separate equations were needed for each EZ; and (4), based on the ranking of performance in terms of accuracy and precision, recommend the bes alternative for estimating total bole volume in individual P. occidentalis trees.

2. Materials and Methods

The study area is in the north-central part of Cordillera Central, DR, covering approximately 1800 km2 (Figure 1). Even-aged natural stands of P. occidentalis are located within three ecological zones according to the Holdridge [10] classification: Subtropical Dry Forest (Dry Zone), Subtropical Humid Forest (Intermediate Zone), and Subtropical Very Humid Forest (Humid Zone). Average elevations above sea level are 500, 650, and 800 m, respectively [11]. The climate varies depending on the altitude and precipitation. The average annual temperature is between 12 °C and 24 °C. These forests usually develop in shallow, carbonate, lateritic, low-producing soils with rugged topography.

2.1. Tree Data Sets

Three data sources were used for fitting total bole volume outside bark ( V S O B ) models, one for each ecological zone. Trees from either dominant, intermediate, or overtopped classes in each zone were selected for destructive sampling, given that their phytosanitary status would allow it. After recording DBH outside bark at 1.30 m from the ground, measured with a diameter tape (0.1 cm precision), trees were felled as close to the ground as possible, and H was measured with open-reel tape (0.1 cm precision). Besides DBH, other diameters were recorded at every meter along the bole from the base (10 cm above ground) up to the apex. Stem analysis data for model fitting purposes included 37 trees from the Subtropical Dry Forest, 48 from the Intermediate Zone, and 72 from the Humid Zone. In addition, an independent data set collected from each zone was available to validate the resulting best total V S O B model. These data included 85, 90, and 75 independently measured P. occidentalis trees for the Dry, Intermediate, and Humid Zones. Table 1 showcases descriptive statistics for both the fitting and validation data set.

2.2. Data Exploration

To visualize the form of the relation between the dependent variable V S O B and the chosen explanatory variables DBH and total tree height (H), we proceeded to merge corresponding data from all three EZs and plot the V S O B against DBH first, and later against the combined variable diameter square times total tree height, D2H.

2.3. Approaches to Individual Tree Volume Prediction

2.3.1. Indicator Variables Analysis

An indicator variable analysis was conducted to check if three different equations were necessary for each EZ. We searched for statistically significant differences in the intercept and slope parameters in a regression equation fitted to the data for all three EZs, employing a single combined effect variable (D2H) and indicator variables. A statistically significant test for the intercept and/or the slope for two zones would indicate that each would require different equations. The following model was fitted to indicator variables and the continuous variable D2H in the indicator variable analysis:
V i b = β 0 + δ 1 Z 1 + δ 2 Z 2 + β 1 D 2 H + δ 3 Z 1 D 2 H + δ 4 Z 2 D 2 H + ε
where
  • D ,   H ,   V i b ,   ε are as previously defined;
  • Z 1 , Z 2 are dichotomous variables;
  • β 0 ,   β 1 ,   δ 1 ,   δ 2 ,   δ 3 ,   δ 4 are the parameters to be estimated.

2.3.2. Total Bole Volume Model Fitting

Observed outside-bark bole volume (m3) computations were performed using Smalian’s formula [12] for each 1 m section from the base to the apex, except for the last portion where the cone formula was used. The outside-bark volume computed for each section was then summed up for each tree to calculate V S O B .
To estimate/predict V S O B of Pinus occidentalis individual trees, we choose two models.
The combined variable equation (CV):
V S O B = β 0 + β 1 D 2 H + ε
  • Schumacher and Hall’s [7] equation (SH):
V S O B = β 1 × D β 2 H β 3 × ε
where
  • V s o b = total stem volume content outside bark (m3);
  • D = normal diameter at 1.30 m from the ground outside bark (cm);
  • H = total tree height (m);
  • L n = natural logarithm;
  • ε = error term;
  • β 0 ,   β 1 ,   β 2 = coefficients to be estimated.
We fitted five variants to our model [2] and four variants to our model [3], including linear, nonlinear, and weighted regression estimation techniques. The ordinary and weighted least squares method fitted the five Model [2] variants. These variants are as follows:
  • (CV01) Original equation
  • Weighted linear regression using four different weights:
    • (CV02) Weight 1 = 1/fitted values from the original linear regression between the dependent variable “observed volume” (Vol) and the predictor normal diameter squared times total tree height (D2H);
    • (CV03) Weight 2 = 1/fitted value resulting from fitting the absolute values of original residuals against the fitted values of original combined variable regression;
    • (CV04) Weight 3 = 1/fitted value resulting from fitting squared values of original residuals against the fitted values of original combined variable regression;
    • (CV05) Weight 4 = 1/ D 2 c , where the variance of ε is assumed to be proportional to D 2 c [13].
Where
  • CV0i = variant identification code for model [2];
  • C = exponent to be assumed or estimated.
The four variants from model [3] were fitted by the ordinary least squares and nonlinear methods, using multiple regression techniques between the dependent variable, V S O B , and the predictor variables, DBH and total tree height (H).
  • (SH01) De-transformation of the logarithmic conversion ( L n V S O B = L n ( β 0 ) + β 1 L n D + β 2 L n H + L n ε ), solved by employing linear regression and correcting for bias. The correction is achieved by adding one-half of the estimated variance from the fitted regression before exponentiation [14]. The resulting expression is as follows:
    V S O B = e β ^ 0 + β ^ 1 × l n D + β ^ 2 × l n H + σ ^ 2 2
    where
    • V S O B = corrected estimate of the stem volume outside bark;
    • μ ^ = mean volume outside bark estimated in log scale;
    • σ ^ 2 2 = half-estimated variance in log scale.
  • (SH02) Nonlinear SH (model [2]) version.
  • (SH03) Nonlinear weighted version SH version assuming exponent c = 2;
  • (SH04) Nonlinear weighted SH version with modeled variance (exponent c), where the variance of ε is assumed to be proportional to D 2 c [13].

2.3.3. Statistical Analysis

Weighted linear and nonlinear least squares were used to maximize the efficiency of parameter estimation. All tests on the full model parameters were conducted at α = 0.05. Data analysis and model development procedures were performed using lm, nls, and nlme commands in RStudio [15] to obtain parameter estimates for V S O B .
Starting values of the coefficients for nonlinear and weighted nonlinear variants were obtained by applying ordinary least squares to the log-transformed data, ensuring faster iteration. Even though log-transformed and nonlinear models are not mathematically equivalent, the coefficients of the former estimated by multiple regression may serve as starting values for the algorithm that estimates the coefficients of the latter [13].
To estimate coefficient c for the weighting of the observations in variants CV05 and SH04, observations were divided into five DBH classes containing approximately the same number of observations. Then, we calculated the standard deviation of V S O B in each D class. Following a methodology from Picard et al. [13], we plotted the standard deviation of V S O B against the median DBH in each of the five classes in the log scale. The five points on the plot should be roughly aligned along a straight line to confirm that the power model was appropriate for modeling the residual variance. If that were the case, we would proceed to fit a linear regression of the log of the standard deviation of V S O B on the log of the median D for each class. The slope of such regression corresponds to exponent c. The standard deviation of the stem volume would be approximately proportional to D2C, and a weighting of the observations would be inversely proportional to it.
Statistical analyses performed on model [3] variant SH01 included initially working on the log-transformed data and fitting an ordinary least squares (OLS) multiple regression of ln(V) against ln(D) and ln(H), and then transforming to original units the coefficient estimates of this initial equation and correcting it for bias by adding one-half of the estimated variance from the fitted regression before exponentiation.

2.3.4. Evaluation Criteria

Model Validation and Goodness of Fit Statistics

The goodness-of-fit statistics used to determine how well the regression functions fitted the sample data were (1) root mean square error (RMSE); (2) Bias; (3) the sum of squared relative residuals (SSRR); (4) the residual variance estimator (RVE), and (5) Akaike Information Criteria (AIC), used only in the fitting phase to compare different models and balance the goodness of fit with model complexity.
The “validation” statistics used to determine how well the regression functions performed on the independent data representing the population were (1) RMSE, (2) Bias, (3) SSRR, and (4) RVE. The best model for total bole volume in each of the three zones was selected based on the ranking of these evaluation criteria. We also considered the significance of parameter estimates [13,16]. The computational formulas for the goodness-of-fit statistics are as follows:
A I C = 2   l n l ( θ ^ ) + 2 q
R M S E = i = 1 n V s o b V ^ s o b 2 n p
B i a s = i = 1 n V s o b V ^ s o b n
S S R R = i = 1 n V s o b V ^ s o b V s o b 2
R V E = n n p 1 1 R 2 S Y 2
where
  • l : is the model’s likelihood;
  • p: is the number of free parameters estimated;
  • V s o b : is the observed stem wood volume outside bark;
  • V ^ s o b : is the estimated stem wood volume outside bark;
  • n : is the total number of observations;
  • S Y 2 : is the empirical variance of the response variable.

Ranking of Models

To rank the nine modeling variants tested in each EZ, fit and validation goodness-of-fit statistics values were ranked. Rank valued No. 1 was linked to the best value for each of these statistics, rank No. 2 the second best, and so on. The overall rank for each model variant was determined by summing the ranks for the various goodness-of-fit statistics for the total volume and then choosing the lowest sum as the best model variant.

Residual and Quantile-Quantile Plot Graphs

Scatterplots were constructed to check regression assumptions for the best-ranked model variants in the Dry Zone and the Combined Intermediate and Humid Zones to check that the hypotheses assumed for the residuals were satisfied.
The hypothesis that the residuals were independent has already been satisfied due to the sampling plan adopted. The constant variance hypothesis of the residuals was visually checked by plotting the cluster of points for the residuals ε i = V o l S O B V o l ^ S O B in function to the predicted values. The hypothesis that the residuals are normally distributed was visually inspected with the quantile–quantile graphs, plotting the residuals’ empirical quantiles against the theoretical quantiles of the standard normal distribution. To further assess how well predictions from model variants align with the actual data, we plotted observed-versus-predicted stem volume values and observed volumes versus predictor variables.

3. Results

3.1. Data Exploration

On average, sampled trees were smaller in terms of DBH, H, and V S O B in the Dry, and largest in the Humid. DBH ranged from 11.50 to 42.00 cm in the Dry Zone, 16.50 to 53.50 cm in the Intermediate Zone, and 21.50 to 46.50 cm in the Humid Zone. Following the same zone order, H from 10.30 to 24.00 m, 14.00 to 25.10 m, and 14.50 to 35.00 m. Volume outside bark ranged from 0.06 to 1.10 m3, 0.32 to 1.30 m3, and 0.20 to 2.73 m3, respectively (Table 1).
The plotted points in the left panel of Figure 2 show that the relationship between DBH and volume is not linear. Volume variance increases as DBH increases. The relationship between V S O B and the combined effect variable D2H, shown on the right side of Figure 2, shows that the relationship between these two variables is linear; but, as before, the volume variance increases with D2H.

3.2. Indicator Variables Analysis (IVA)

The results of the Indicator variable analysis show that the Humid and Dry Zones have different intercepts and the same slope. Humid and Intermediate Zones have the same intercept and slope, and Dry and Intermediate Zones have the same intercept but different slopes (Table 2). Therefore, the Humid and Intermediate Zones can have the same equations to estimate the V S O B of individual P. occidentalis trees, while the Dry Zone requires its own equations. Tree stem form has been observed to change in different zones with variation attributable to environmental factors [1]. Based on these indicator variable analysis results, we combined the observed data from the Humid and Intermediate Zones and fitted the tested models to the combined data.
Merged fitting data from the Humid Zone and Intermediate Zone totaled 121 observations. DBH ranged from 16.50 to 53.50 cm, H ranged from 14.00 to 35.00 m, and V o l S o b ranged from 0.20 to 2.76 m3. Respective averages and sample standard deviations (within parenthesis) were 31.77 (±7.05) cm, 22.67 (±4.73) m and 0.84 (±0.48) m3. After merging data from both zones, the validation data set had 165 observations with DBH ranging from 10.60 to 54.20 cm, H ranging from 9.40 to 27.80 m, and V o l S o b from 0.11 to 2.22 m3. Respective averages and sample standard deviations were 28.53 (±8.15) cm, 20.27 (±3.78) m, and 0.59 (±0.36) m3.

3.3. Total Bole Volume Model Fitting in the Dry Ecological Zone and the Combined Intermediate and Humid Ecological Zone

The SH02 variant, fitted by nonlinear least square (NLS) procedures, was ranked number one for estimating V S O B  P. occidentalis in the Dry Zone, while in the Combined Intermediate and Humid Zone, the best variant was SH03, which was fitted by weighted nonlinear least square (WNLS) procedures. Based on the ranks of the goodness-of-fit statistics criteria applied to the data fitting and validation phases, these results indicate the superiority of the SH model variants in this study. The fitting and validation statistic values, along with rank sums and the overall ranks for each modeling variant in each EZ, are presented in Table 3 and Table 4.
Considering both the fitting and validation phases, SH02 in the Dry Zone was ranked No. 1 in terms of RMSE (0.0146) and RVE (0.0002) and No. 2 in AIC criteria (−199.86) in the fitting phase (Table 3). Likewise, it was ranked No. 1 in terms of RMSE (0.0702), Bias (−0.0099), and RVE (0.0030) in the validation phase. Percent bias calculated as V S O B V ^ S O B / V S O B × 100 in the fitting phase indicated that SH02 tends to slightly overestimate V S O B by approximately −0.059%. Percent RMSE was calculated at 4.77% for SH02, suggesting that predictions deviate by about 4.77% from the mean of the actual values. We consider these values to be adequate for estimating V S O B for the individual trees of P. occidentalis in this EZ.
SH02 had an overall sum rank of 31, being the lowest ranked of all volume variants tested in the Dry Zone. It was followed by SH04, the fourth variant of model (3) with a sum rank of 33 and fitted by WNLS procedure, where exponent “c” in the variance model ( V a r ε = k D B H c 2 ) was a parameter that needed to be estimated.
SH03 in the Combined Intermediate and Humid Zone was ranked No. 1, achieving the lowest sum in the ranking (23) among the nine variants tested in this combined EZ (Table 4). For weights, the conditional standard deviation of V S O B derived from DBH was proportional to DBH4. In the fitting phase, SH03 first regarded AIC criteria (−311.66). Similarly, it was ranked No. 1 in terms of RMSE (0.0943), Bias (−0.0614), and RVE (0.0059) in the validation phase. Percent bias calculated as V S O B V ^ S O B / V S O B × 100 in the fitting phase indicated that SH03 tends to underestimate V S O B by approximately 0.025% on average. Percentage RMSE was calculated at 8.96%, suggesting predictions that deviate by about 8.96% from the mean of the actual values for estimating V S O B for the individual trees of P. occidentalis in the Combined Intermediate and Humid Zone. Regarding the sum rank, SH03 was followed by modeling variant SH01, the unweighted nonlinear version of the SH function, with a sum rank of 35.
The parameterization of SH02 and SH03, respectively, are as follows:
V S O B D E Z = 0.00005815 × D 1.7802 H 1.0787
V S O B C I H E Z = 0.00005669 × D 1.7019 H 1.0786
Although different fitting procedures estimated them, the intercept and coefficient corresponding to the Dry Zone and Combined Intermediate and Humid Zone tree height are very similar. The intercepts differed by 2.54%, and the tree height (H) coefficients differed by 0.009%. The coefficients for D differed by 4.49%.
In checking regression assumptions, our plots show no curvature. However, there is one point beyond the cutoff limit of two standard deviations in the top left panel corresponding to the Dry Zone. We assumed modeling variants SH02 and SH03 do not seriously contradict the constant variance assumption. Even though the quantile–quantile plots of the residuals appear to have a slight structure in both zones, most of the points are aligned along a straight line and remain relatively consistent across all the levels of the predictor variables, indicating that the constant variance assumption is sustained.
All nine variants tested from models (2) and (3) in the Dry and Combined Intermediate and Humid ecological zones resulted in parameter estimates that are statistically significant and logically consistent for outside-bark volume estimation (Table 5 and Table 6). The estimates for intercept, interpreted as the average value of the total stem volume accumulated by a tree until it reaches breast height (when DBH is zero and total height is 1.3 m) were positive and significantly different from zero. All other parameters were consistent in terms of sign and magnitude.
In Figure 3, scatter plots of standardized residuals plotted against fitted values (panels a and b) and quantile–quantile plots (panels c and d) are shown as tools to check the constant variance and normal distribution hypotheses of the residuals in the Dry Zone and Combined Intermediate and Humid Zone using modeling variants SH02 and SH03. The standardized residuals plotted against fitted values show randomly clustered values around the 0 line for outer-bark volumes in the Dry Zone (a) and the Combined Intermediate and Humid Zone (b). The points around the line in the quantile–quantile plots remain relatively constant across all levels of the predictor variables, confirming the constant variance assumption and closely following the observed values, suggesting that these models capture the variability in the data well.
In Figure 4, the top panel (a) shows predictions carried out by modeling variant SH02 against the volume predicted by the SH02 variant. Likewise, the bottom panel (b) shows the predictions carried out by modeling variant SH03 against the predicted volume in the Combined Intermediate and Humid Zone. These plots support evidence of the agreement between observations and predictions. SH02 and SH03 efficiently estimated the stem volume outside bark for P. occidentalis in both ecological zones and effectively captured the variability in the data.

4. Discussion

Alvarado-Segura et al. [17] used a nonlinear version of the SH model to estimate the stem volume of Pinus patula (Schl. et Cham.) in Hidalgo, Mexico, and reported an RMSE of 0.0914. As a measure of the average magnitude of the residual error, RMSE lower values indicate better goodness-of-fit of a regression model [18]. Our RMSE statistic value is much smaller and better for our SH02 modeling variant (0.0146), while the corresponding value for SH03 (0.0943) is close to that reported by these authors. The SH model has been found to perform well in predicting stem volume for many species and genera in diverse environments [19]. Castillo-López et al. [20] and Valerio-Hernández et al. [21] used this model in studying volume equations for pine species in Mexico and Nicaragua. Based on the evaluation criteria, this study’s nonlinear version (SH02) and the nonlinear weighted variant (SH03) of the S&H model performed best.
All nine variants tested from models (2) and (3) in both EZs resulted in parameter estimates that are statistically significant and logically consistent for outside-bark volume estimation (Table 5 and Table 6). The observed relationship between the dependent variable stem volume outside bark and the predictors is unlikely to have occurred by chance [18]. The estimates for intercept, interpreted as the average value of the total stem volume accumulated by a tree until it reaches breast height (when DBH is zero and total height is 1.3 m) were positive and significantly different from zero. All other parameters were consistent in terms of sign and magnitude.
The best-ranked Combined Variable model variant was CV02, which achieved an overall ranking of 3 in the Dry Zone and 4 in the Combined Intermediate and Humid Zone. Values for intercepting the coefficients of variants in model (2) resemble, in a certain way, the shape and form of a geometric solid, with higher numbers representing better and more cylindrical form, and thus, greater volume [12]. On average, model (2) variant intercepts are lower (105.77%) in the Dry Zone than their counterparts in the Combined Intermediate and Humid Zone, indicating that the stem volume accumulated by individual P. occidentalis trees below breast height is lower in the former zone and that trees are smaller. A t-test assuming unequal variances indicates that intercepts from the two zones are statistically and significantly different (p = 0.0026). The slope coefficients for D2H in the same model (2) variants are also statistically and significantly different (p = 0.000126).
In mathematics, the volume of a circular base solid with base diameter D and height H is expressed as follows: V = β D 2 H [22]. If the diameter is measured in cm, β = π 40,000 for a cylinder, π 80,000 for a paraboloid, π 120,000 for a cone, and π 160,000 for a neiloid. This β is equivalent to β 1 , the coefficient in the second term of the CV. The averages of the five β 1 values from our results are equivalent to 0.0000347 in the Dry Zone, and 0.0000321 in the Combined Intermediate and Humid Zone, differing by about 12.4% and 20.47% off the perfect paraboloid, respectively. Similar values for β 1 using the CV equation in fitting data from 150 P. patula trees were found by Alvarado-Segura et al. [17] in Mexico. Therefore, tree shapes in our study approximate a paraboloid solid, which is described mathematically by the expression π/80,000 = 0.0000393 (diameters are in cm units). This contrasts with the results of Sharma [22] who found that the shape of the trees is not a solid described by a cylinder, paraboloid, cone, or neiloid while assessing the tree volume of twenty-five species in the natural stands of Canada and northeastern United States, including balsan poplar, eastern white cedar, Engelman Spruce and European larch. These are temperate forest trees, and the results may not be comparable, but it has been assumed that all trees should approximate one of these mathematical shapes, regardless of species.
Differences in the parameters of the CV models between EZs may result from differences in the height-dbh relationship, tree form, tree taper, or a combination of these factors [1], although all were logically consistent. These results agree with the findings reported by Sharma [23], who employed the CV equation to compare goodness-of-fit statistics, logical consistency, and the predictive accuracy of several models in Pinus resinosa Sol. ex Aiton. Exploratory analyses of variance not reported in the study showed that the form and quotient coefficients of P. occidentalis trees in these zones were statistically significant. That, combined with higher values for the β 0 coefficient in the Humid Zone, may indicate that the form is better there and, therefore, better suited for wood production.

5. Conclusions

Nine model variants, five from the CV and four from the SH, were evaluated in their capacity to estimate the stem volume outside bark (SVOB) of individual P. occidentalis trees growing in natural stands within three EZs in La Sierra, D.R. All the variants produced results consistent with those of other pine species. All modeling variants gave good results on the fitting data set, but the performance was poorer when used in the validation data set. Ercanli et al. [24] and Sahin [25] encountered the same situation while modeling the tree volume of conifer species in Turkey.
An indicator variable analysis indicated that only one equation was required to estimate the total stem volume outside bark in the Intermediate and Humid Zones. Therefore, observations from these two EZs were combined to evaluate the different modeling variants further.
The CV model variants were logically consistent in parameter estimation. They complied with regression assumptions but were inferior to the chosen SHVE model variants in terms of ranking based on goodness-of-fit statistics and predictive ability. Modeling variants SH02 and SH03, strategies for modeling V S O B individual P. occidentalis trees, were efficient and selected as the volume equations for the Dry Zone and Combined Intermediate and Humid Zone, respectively.
The maximum bias among all modeling variants in the Dry Zone was 0.20% and 5.54% in the fitting and validation phases, respectively. Likewise, in the Combined Intermediate and Humid Zone, bias was 0.007% and 8.84%, respectively. According to Sharma [22], volume equations resulting in a bias larger than 10% are very imprecise, leading to poor performance and unreliable predictions. It should not be recommended for forest management decision-making. The minimum observed variability explained by all variants in both zones was 94% for outside-bark volume estimates.
Developing forest management strategies requires accurate tree volume estimates, employing tools such as tree volume models. These volume models are simpler to use if the entire stem volume is of interest [26]. The assumptions underlying the regression methods employed in the study for constructing these nine volume model variants were validated. The results showed that these modeling variants performed well in both the fitting and validation phases, and therefore, we propose their use as useful tools for predicting the total stem volume outside bark in these EZs within La Sierra, D.R. The equations presented here provide more scientifically accurate and consistent predictions of individual P. occidentalis trees. They are a valuable tool for supporting present and future management decisions.

Author Contributions

Conceptualization, S.W.B.-L. and J.G.T.-H.; methodology, S.W.B.-L. and L.R.C.-R.; validation, S.W.B.-L. and L.R.C.-R.; formal analysis, S.W.B.-L.; investigation, S.W.B.-L.; data curation, S.W.B.-L. and J.G.T.-H.; writing—original draft preparation, S.W.B.-L., L.R.C.-R. and J.G.T.-H.; writing—review and editing, S.W.B.-L.; project administration, S.W.B.-L.; funding acquisition, S.W.B.-L. and L.R.C.-R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Higher Education, Science and Technology of the Dominican Republic, National Fund for Innovation and Scientific and Technological Development grant number FONDOCyT 2018-2019-2B4-088. The APC was funded by Pontificia Universidad Católica Madre y Maestra (No. 121-2019-PUCMM).

Data Availability Statement

The authors will make the raw data supporting this article’s conclusions available upon request.

Acknowledgments

We recognize the effort and contribution of the acquisitions department of PUCMM and the administrative staff of the Vice-Rectory of Research and Innovation at PUCMM. We are also grateful to José Reyes and Roque Herrera, who assisted us in field data collection.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kangas, A.; Pitkänen, T.P.; Mehtätalo, L.; Heikkinen, J. Mixed linear and non-linear tree volume models with regional parameters to main tree species in Finland. For. Int. J. For. Res. 2023, 96, 188–206. [Google Scholar] [CrossRef]
  2. Li, R.; Weiskittel, A.R. Comparison of model forms for estimating stem taper and volume in the primary conifer species of the North American Acadian Region. Ann. For. Sci. 2010, 67, 302. [Google Scholar] [CrossRef]
  3. Clutter, J.L.; Fortson, J.C.; Pienaar, L.V.; Brister, G.H.; Bailey, R.L. Timber Management: A Quantitative Approach, 1st ed.; Krieger Publishing Company: Malabar, FL, USA, 1983; 333p. [Google Scholar]
  4. Bates, D.M.; Watts, D.G. Nonlinear Regression Analysis and Its Applications, 1st ed.; Wiley India Pvt. Ltd.: New Delhi, India; John Wiley & Sons, Inc.: New York, NY, USA, 2014; 365p. [Google Scholar]
  5. Vibrans, A.C.; Moser, P.; Oliveira, L.Z.; Maçaneiro, J.P. Generic and specific stem volume models for three subtropical forest types in southern Brazil. Ann. For. Sci. 2016, 72, 865–874. [Google Scholar] [CrossRef]
  6. MMARN—Ministerio de Medio Ambiente y Recursos Naturales Inventario Nacional Forestal de la República Dominicana. Programa Regional de Reducción de Emisiones de la Deforestación y Degradación de Bosques en Centroamérica y Republica Dominicana (REDD III); MMARN: Santo Domingo, Dominican Republic, 2021; 292p. [Google Scholar]
  7. Schumacher, F.X.; Hall, F.S. Logarithmic expression of timber-tree volume. J. Agric. Res. 1933, 47, 719–734. [Google Scholar]
  8. Bennett, F.A.; McGee, C.E.; Clutter, J.L. Yield of Old-Field Slash Pine Plantations; No. 107; U.S. Department of Agriculture Forest Service: Washington, DC, USA, 1959.
  9. Azevedo, G.B.; Tomiazzi, H.V.; Azevedo, S.; Pereira, L.; Teodoro, R.; Pereira de Souza, T.; Silva, T.; Philipe, B.; Guerra, S. Multi-volume modeling of Eucalyptus trees using regression and artificial neural networks. PLoS ONE 2020, 15, e0238703. [Google Scholar] [CrossRef] [PubMed]
  10. Holdridge, L. Ecología Basada en Zonas de Vida, 1st ed.; Instituto Interamericano de Cooperación para la Agricultura: San José, Costa Rica, 1987; 304p. [Google Scholar]
  11. Bueno-López, S.W. Understanding Growth and Yield of Pinus occidentalis, Sw. in La Sierra, Dominican Republic. Doctor Dissertation, State University of New York, College of Environmental Science and Forestry, Syracuse, NY, USA, 2009; 256p. [Google Scholar]
  12. Burkhart, H.E.; Tomé, M. Modeling Forest Trees and Stands, 1st ed.; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; 457p. [Google Scholar] [CrossRef]
  13. Picard, N.; Saint-André, L.; Henry, M. Manual for Building Tree Volume and Biomass Allometric Equations: From Field Measurement to Prediction, 1st ed.; Food and Agricultural Organization of the United Nations, Rome, and Centre de Coopération Internationale en Recherche Agronomique pour le Développement: Montpellier, France, 2012; 215p. [Google Scholar]
  14. Flewelling, J.W.; Pienaar, L.V. Multiplicative regression with lognormal errors. For. Sci. 1981, 27, 281–289. [Google Scholar]
  15. RStudio Team. RStudio: Integrated Development for R; RStudio, PBC: Boston, MA, USA, 2020. Available online: http://www.rstudio.com/ (accessed on 28 April 2024).
  16. Hu, N. Development of a New Variable-Form Taper Equation to Investigate Differences in Stem Form Following Release in Eastern White Pine (Pinus strobus L.). Master’s Thesis, State University of New York, College of Environmental Science and Forestry, Syracuse, NY, USA, 2003; 97p. [Google Scholar]
  17. Alvarado-Segura, A.A.; Zamudio-Sánchez, F.J.; De La Cruz-De La Cruz, K.I. A Procedure For Choosing Tree-Stem Volume Equations Previously Fitted in a Forest. J. Sustain. For. 2020, 39, 595–607. [Google Scholar] [CrossRef]
  18. Schabenberger, O.; Pierce, F.J. Contemporary Statistical Models for the Plant and Soil Sciences; CRC Press: Boca Raton, FL, USA, 2002; p. 738. [Google Scholar]
  19. Abreu, J.C.; Soares, C.P.B.; Leite, H.G.; Binoti, D.H.B.; Silva, G.F. Alternatives to estimate the volume of individual trees in forest formations in the state of Minas Gerais-Brazil. CERNE 2020, 26, 393–402. [Google Scholar] [CrossRef]
  20. Castillo-López, A.; Quiñonez-Barraza, G.; Diéguez-Aranda, U.; Corral-Rivas, J.J. Compatible Taper and Volume Systems Based on Volume Ratio Models for Four Pine Species in Oaxaca Mexico. Forests 2021, 12, 145. [Google Scholar] [CrossRef]
  21. Valerio-HernánDry Zone, L.A.; Campos-Vanegas, W.A.; Cruz-Tórrez, L.E.; Pena-Ortiz, J.A.; Vargas-Larreta, B. Improving Volume and Biomass Equations for Pinus oocarpa in Nicaragua. Forests 2024, 15, 309. [Google Scholar] [CrossRef]
  22. Sharma, M. Total and Merchantable Volume Equations for 25 Commercial Tree Species Grown in Canada and the Northeastern United States. Forests 2021, 12, 1270. [Google Scholar] [CrossRef]
  23. Sharma, M. Increasing Volumetric Prediction Accuracy: An Essential Prerequisite for End-Product Forecasting in Red Pine. Forests 2020, 11, 1050. [Google Scholar] [CrossRef]
  24. Ercanli, I.; Senyurt, M.; Bolat, F. A major challenge to machine learning models: Compatible predictions with biological realism in forestry: A case study of individual tree volume. In Proceedings of the 3rd International Conference on Environment and Forest Conservation (ICEFC), Kastamonu, Turkey, 21–23 February 2022; p. 39. [Google Scholar]
  25. Sahin, A. Analyzing regression models and multi-layer artificial neural network models for estimating taper and tree volume in Crimean pine forests. iForest 2024, 17, 36–44. [Google Scholar] [CrossRef]
  26. Sharma, M. Inside and outside bark volume models for jack pine (Pinus banksiana) and black spruce (Picea mariana) plantations in Ontario, Canada. For. Chron. 2019, 95, 50–57. [Google Scholar] [CrossRef]
Figure 1. Study area within the La Sierra region in the north-central portion of Cordillera Central, Dominican Republic.
Figure 1. Study area within the La Sierra region in the north-central portion of Cordillera Central, Dominican Republic.
Forests 15 01052 g001
Figure 2. Scatterplots of the relation between V S O B and DBH (a) and D2H (b). On the left, the relationship is nonlinear and presents heteroscedasticity. On the right, the relationship is linear, but heteroscedasticity persists.
Figure 2. Scatterplots of the relation between V S O B and DBH (a) and D2H (b). On the left, the relationship is nonlinear and presents heteroscedasticity. On the right, the relationship is linear, but heteroscedasticity persists.
Forests 15 01052 g002
Figure 3. Residuals versus fitted values in the Dry Zone (a) and Combined Intermediate and Humid Zone (b) are shown on the left. Sample versus theoretical quantiles in the Dry Zone (c) and Combined Intermediate and Humid Zone (d) are depicted on the right.
Figure 3. Residuals versus fitted values in the Dry Zone (a) and Combined Intermediate and Humid Zone (b) are shown on the left. Sample versus theoretical quantiles in the Dry Zone (c) and Combined Intermediate and Humid Zone (d) are depicted on the right.
Forests 15 01052 g003
Figure 4. Predicted outside-bark volumes plotted against corresponding observed volumes in the Dry Ecological Zone (a) and Combined Intermediate and Humid Zone (b) are illustrated. SH02 and SH03 efficiently estimated the stem volume outside bark for P. occidentalis in both ecological zones and effectively captured the variability in the data.
Figure 4. Predicted outside-bark volumes plotted against corresponding observed volumes in the Dry Ecological Zone (a) and Combined Intermediate and Humid Zone (b) are illustrated. SH02 and SH03 efficiently estimated the stem volume outside bark for P. occidentalis in both ecological zones and effectively captured the variability in the data.
Forests 15 01052 g004
Table 1. Descriptive statistics of tree measurements in each ecological zone for the fitting and validation data sets used in the study.
Table 1. Descriptive statistics of tree measurements in each ecological zone for the fitting and validation data sets used in the study.
VariableEcological ZonenMeanStd DevMinimumMaximum
Fitting Data Set
Diameter (cm)Dry Zone3722.136.5711.5042.00
Intermediate Zone4832.638.3416.5053.50
Humid Zone7231.185.1821.5046.50
Height (m)Dry Zone3716.492.9210.3024.00
Intermediate Zone4819.782.7814.0025.10
Humid Zone7224.654.7614.5035.00
Volume (m3)Dry Zone370.340.220.061.10
Intermediate Zone480.660.240.321.30
Humid Zone720.980.570.202.76
Validation Data Set
Diameter (cm)Dry Zone8521.327.958.0042.10
Intermediate Zone9027.268.2011.0054.20
Humid Zone7530.067.8710.6050.10
Height (m)Dry Zone8516.334.467.3026.10
Intermediate Zone9019.964.0610.1027.80
Humid Zone7520.653.409.4027.40
Volume (m3)Dry Zone850.350.250.101.23
Intermediate Zone900.550.350.122.22
Humid Zone750.640.360.111.81
n: number of observations; Std Dev: Standard Deviation.
Table 2. Statistical test results for the intercept and slope coefficients from the indicator variable analysis of the dependent variable stem volume outside bark in the three ecological zones.
Table 2. Statistical test results for the intercept and slope coefficients from the indicator variable analysis of the dependent variable stem volume outside bark in the three ecological zones.
ZoneInterceptSlope
Humid versus DryDifferent: p-value = 0.0277Same: p-value = 0.1414
Humid versus IntermediateSame: p-value = 0.4851Same: p-value = 0.974
Dry versus IntermediateSame: p-value = 0.104Different: p-value = 0.0294
Table 3. Goodness-of-fit statistics and rank (ranking value in parenthesis below corresponding statistic) were obtained in the fitting and validation stage of five variants of the combined variable equation and four of the Schumacher and Hall [7] model to estimate the stem volume outside bark of P. occidentalis trees in the Dry Zone within La Sierra, Dominican Republic.
Table 3. Goodness-of-fit statistics and rank (ranking value in parenthesis below corresponding statistic) were obtained in the fitting and validation stage of five variants of the combined variable equation and four of the Schumacher and Hall [7] model to estimate the stem volume outside bark of P. occidentalis trees in the Dry Zone within La Sierra, Dominican Republic.
Fit StatisticsValidation StatisticsRanking
ModelVariant CodeRMSE (Rank)BIAS (Rank)SSRR (Rank)RVE (Rank)AIC (Rank)RMSE (Rank)BIAS (Rank)SSRR (Rank)RVE (Rank)Sum RankOverall Rank
Model (2): Effect Variable D2HCV011.61E−023.87E−191.06E−012.82E−04−1.95E+027.92E−02−1.32E−025.98E+003.52E−03466
(5)(1)(9)(9)(8)(4)(5)(1)(4)
CV021.62E−029.96E−199.66E−022.80E−04−2.07E+028.11E−02−1.28E−026.41E+003.63E−03383
(6)(2)(6)(5)(4)(5)(3)(2)(5)
CV031.62E−02−2.41E−049.67E−022.80E−04−2.06E+028.21E−02−1.37E−026.42E+003.71E−03537
(7)(6)(7)(6)(5)(7)(6)(3)(6)
CV041.63E−02−5.93E−049.57E−022.81E−04−2.12E+028.31E−02−1.38E−026.59E+003.77E−03588
(9)(8)(4)(8)(1)(8)(7)(5)(8)
CV051.63E−02−6.83E−049.62E−022.81E−04−2.11E+028.32E−02−1.42E−026.51E+003.77E−03619
(8)(9)(5)(7)(2)(9)(8)(4)(9)
Model (3): Effect Variables D, HSH011.48E−022.70E−049.55E−022.49E−04−1.08E+027.36E−02−1.30E−027.20E+003.19E−03414
(4)(7)(1)(4)(9)(3)(4)(6)(3)
SH021.46E−02−1.81E−049.85E−022.39E−04−2.00E+027.02E−02−9.89E−037.28E+002.97E−03311
(1)(5)(8)(1)(6)(1)(1)(7)(1)
SH031.47E−02−6.23E−069.56E−022.43E−04−2.11E+028.14E−02−1.89E−027.33E+003.71E−03445
(3)(3)(2)(3)(3)(6)(9)(8)(7)
SH041.47E−02−2.54E−059.56E−022.42E−04−1.97E+027.15E−02−1.02E−027.33E+003.06E−03332
(2)(4)(3)(2)(7)(2)(2)(9)(2)
CV01: Combined Variable function first variant; CV02: Combined Variable function second variant; CV03: Combined Variable function third variant; CV04: Combined Variable function fourth variant; CV05: Combined Variable function fifth variant; SH01: De-transformed S&H Model with Bias Corrected variant; SH02: S&H Non-linear Model variant; SH03: S&H Weighted Non-linear Model Assuming Exponent C variant; SH04: S&H Weighted Non-linear Model Modeling Exponent C variant; H: Total Tree Height.
Table 4. Goodness-of-fit statistics and respective ranking obtained in the fitting and validation stage of the five variants of the Combined Variable Equation and four variants of the Schumacher and Hall model [7] to estimate the stem volume outside bark of P. occidentalis trees in the Combined Intermediate and Humid Zone within La Sierra, Dominican Republic.
Table 4. Goodness-of-fit statistics and respective ranking obtained in the fitting and validation stage of the five variants of the Combined Variable Equation and four variants of the Schumacher and Hall model [7] to estimate the stem volume outside bark of P. occidentalis trees in the Combined Intermediate and Humid Zone within La Sierra, Dominican Republic.
Fit StatisticsValidation StatisticsRanking
ModelVariant CodeRMSE (Rank)BIAS (Rank)SSRR (Rank)RVE (Rank)AIC (Rank)RMSE (Rank)BIAS (Rank)SSRR (Rank)RVE (Rank)Sum RankOverall Rank
Model (2): Effect Variable D2HCV017.85E−026.11E−181.06E+006.31E−03−2.67E+021.01E−01−7.29E−022.22E+006.70E−03507
(5)(2)(9)(9)(8)(4)(8)(1)(4)
CV027.87E−02−4.01E−181.01E+006.19E−03−2.99E+021.03E−01−7.06E−022.39E+006.84E−03414
(6)(1)(7)(6)(5)(5)(4)(2)(5)
CV037.89E−02−1.29E−031.01E+006.17E−03−3.01E+021.06E−01−7.24E−022.56E+007.17E−03496
(7)(7)(6)(5)(3)(6)(5)(3)(7)
CV048.08E−02−5.31E−031.00E+006.19E−03−3.11E+021.14E−01−7.27E−023.23E+007.81E−03598
(9)(8)(5)(7)(2)(9)(6)(5)(8)
CV058.07E−02−7.43E−031.02E+006.21E−03−3.00E+021.14E−01−7.43E−023.10E+007.85E−03679
(8)(9)(8)(8)(4)(8)(9)(4)(9)
Model (3): Effect Variables D, HSH017.59E−02−5.22E−049.56E−015.94E−03−2.39E+029.49E−02−6.17E−024.14E+006.02E−03352
(4)(5)(1)(4)(9)(2)(2)(6)(2)
SH027.55E−026.86E−049.57E−015.86E−03−2.74E+029.86E−02−6.53E−024.33E+006.41E−03363
(1)(6)(2)(3)(7)(3)(3)(8)(3)
SH037.56E−022.08E−049.65E−015.83E−03−3.12E+029.43E−02−6.14E−024.22E+005.97E−03231
(2)(4)(4)(2)(1)(1)(1)(7)(1)
SH047.56E−02−9.23E−059.65E−015.81E−03−2.98E+021.08E−01−7.27E−024.75E+007.43E−03455
(3)(3)(3)(1)(6)(7)(7)(9)(6)
CV01: Combined Variable function first variant; CV02: Combined Variable function second variant; CV03: Combined Variable function third variant; CV04: Combined Variable function fourth variant; CV05: Combined Variable function fifth variant; SH01: De-transformed S&H Model with Bias Corrected variant; SH02: S&H Nonlinear Model variant; SH03: S&H Weighted Nonlinear Model Assuming Exponent C variant; SH04: S&H Weighted Nonlinear Model Modeling Exponent C variant; H: Total Tree Height.
Table 5. Parameter estimates, corresponding confidence intervals, residual standard error, and the adjusted coefficient of determination obtained in fitting the five variants of the Combined Variable Equation and four variants of the Schumacher and Hall model [7] to estimate the stem volume outside bark of P. occidentalis trees in the Dry Zone, La Sierra, Dominican Republic.
Table 5. Parameter estimates, corresponding confidence intervals, residual standard error, and the adjusted coefficient of determination obtained in fitting the five variants of the Combined Variable Equation and four variants of the Schumacher and Hall model [7] to estimate the stem volume outside bark of P. occidentalis trees in the Dry Zone, La Sierra, Dominican Republic.
CV (Model (2))S&H (Model (3))
ParametersStatisticsCV01CV02CV03CV04CV05SH01SH02SH03SH04
Residual Est. Error1.66E−022.79E−021.25E+006.07E+016.87E−051.48E−021.52E−023.16E−051.47E−02
Adjusted R29.92E−019.93E−019.92E−019.92E−017.60E−019.93E−019.93E−019.93E−019.93E−01
B0Estimate 1.59E−021.35E−021.34E−021.25E−021.29E−026.14E−055.81E−055.88E−055.84E−05
Lower Bound 95% CI5.54E−036.80E−036.36E−037.69E−037.66E−034.69E−054.29E−054.48E−055.84E−05
Upper Bound 95% CI2.63E−022.02E−022.05E−021.73E−021.81E−028.04E−057.85E−057.71E−055.84E−05
Pr (>|t|) Bo3.66E−032.43E−044.72E−046.80E−061.49E−055.59E−391.02E−071.00E−089.99E−09
B1Estimate 3.44E−053.46E−053.47E−053.48E−053.48E−051.82E+001.78E+001.82E+001.81E+00
Lower Bound 95% CI3.33E−053.37E−053.36E−053.38E−053.38E−051.73E+001.67E+001.72E+001.81E+00
Upper Bound 95% CI3.54E−053.56E−053.57E−053.59E−053.58E−053.63E+001.89E+001.91E+001.81E+00
Pr (>|t|) B11.44E−381.51E−391.11E−387.22E−395.88E−401.52E−301.07E−272.04E−302.97E−30
B2Estimate 1.02E+001.08E+001.04E+001.04E+00
Lower Bound 95% CI 8.79E−019.64E−018.95E−011.04E+00
Upper Bound 95% CI 2.04E+001.19E+001.18E+001.04E+00
Pr (>|t|) B2 3.88E−169.02E−201.75E−161.35E−16
CEstimate 1.74E+00 2.00E+001.90E+00
Table 6. Parameter estimates, corresponding confidence intervals, residual standard error, and adjusted coefficient of determination obtained in fitting the five variants of the Combined Variable Equation and four variants of the Schumacher and Hall model to [7] estimate stem volume outside bark of P. occidentalis trees in the Combined Intermediate and Humid Zone, La Sierra, Dominican Republic.
Table 6. Parameter estimates, corresponding confidence intervals, residual standard error, and adjusted coefficient of determination obtained in fitting the five variants of the Combined Variable Equation and four variants of the Schumacher and Hall model to [7] estimate stem volume outside bark of P. occidentalis trees in the Combined Intermediate and Humid Zone, La Sierra, Dominican Republic.
CV (Model (2))S&H (Model (3))
ParametersStatisticsCV01CV02CV03CV04CV05SH01SH02SH03SH04
Residual Est. Error5.82E−024.81E−024.61E−023.27E−023.63E−026.13E−055.86E−055.67E−055.57E−05
Adjusted R23.04E−022.55E−022.34E−021.65E−021.70E−024.63E−054.33E−054.28E−055.57E−05
B0Estimate 8.60E−027.07E−026.88E−024.89E−025.55E−028.15E−057.92E−057.50E−055.58E−05
Lower Bound 95% CI6.27E−054.80E−051.03E−041.13E−042.98E−043.84E−961.36E−091.33E−101.07E−10
Upper Bound 95% CI3.13E−053.17E−053.19E−053.26E−053.25E−051.82E+001.79E+001.78E+001.79E+00
Pr (>|t|) Bo3.04E−053.07E−053.08E−053.15E−053.14E−051.74E+001.71E+001.70E+001.79E+00
B1Estimate 3.23E−053.28E−053.30E−053.37E−053.36E−053.65E+009.72E−011.86E+001.79E+00
Lower Bound 95% CI5.65E−959.39E−921.66E−885.40E−902.43E−911.67E−754.29E−765.83E−753.44E−75
Upper Bound 95% CI 1.01E+001.06E+001.08E+001.08E+00
Pr (>|t|) B1 9.28E−011.87E+009.89E−011.08E+00
B2Estimate 2.02E+001.14E+001.17E+001.08E+00
Lower Bound 95% CI 5.59E−462.36E−482.71E−473.74E−47
Upper Bound 95% CI 2.03E+00 2.00E+002.20E+00
Pr (>|t|) B27.91E−028.08E−021.17E+001.29E+016.45E−057.59E−027.55E−026.80E−057.56E−02
CEstimate 9.73E−019.69E−019.65E−019.67E−019.68E−019.73E−019.73E−019.74E−019.74E−01
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bueno-López, S.W.; Caraballo-Rojas, L.R.; Torres-Herrera, J.G. Evaluation of Different Modeling Approaches for Estimating Total Bole Volume of Hispaniolan Pine (Pinus occidentalis Swartz) in Different Ecological Zones. Forests 2024, 15, 1052. https://doi.org/10.3390/f15061052

AMA Style

Bueno-López SW, Caraballo-Rojas LR, Torres-Herrera JG. Evaluation of Different Modeling Approaches for Estimating Total Bole Volume of Hispaniolan Pine (Pinus occidentalis Swartz) in Different Ecological Zones. Forests. 2024; 15(6):1052. https://doi.org/10.3390/f15061052

Chicago/Turabian Style

Bueno-López, Santiago W., Luis R. Caraballo-Rojas, and Juan G. Torres-Herrera. 2024. "Evaluation of Different Modeling Approaches for Estimating Total Bole Volume of Hispaniolan Pine (Pinus occidentalis Swartz) in Different Ecological Zones" Forests 15, no. 6: 1052. https://doi.org/10.3390/f15061052

APA Style

Bueno-López, S. W., Caraballo-Rojas, L. R., & Torres-Herrera, J. G. (2024). Evaluation of Different Modeling Approaches for Estimating Total Bole Volume of Hispaniolan Pine (Pinus occidentalis Swartz) in Different Ecological Zones. Forests, 15(6), 1052. https://doi.org/10.3390/f15061052

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