Almost universally, forest inventory and monitoring databases are incomplete, ranging from missing data for only a few records and a few variables, common for small land areas, to missing data for many observations and many variables, common for large land areas. For a wide variety of applications, nearest neighbor (NN) imputation methods have been developed to fill in observations of variables that are missing on some records (Y-variables), using related variables that are available for all records (X-variables). This review attempts to summarize the advantages and weaknesses of NN imputation methods and to give an overview of the NN approaches that have most commonly been used. It also discusses some of the challenges of NN imputation methods. The inclusion of NN imputation methods into standard software packages and the use of consistent notation may improve further development of NN imputation methods. Using X-variables from different data sources provides promising results, but raises the issue of spatial and temporal registration errors. Quantitative measures of the contribution of individual X-variables to the accuracy of imputing the Y-variables are needed. In addition, further research is warranted to verify statistical properties, modify methods to improve statistical properties, and provide variance estimators.
1 March 2008 in Canadian Journal of Forest Research
Using an extensive Douglas-fir data set from southwest Oregon, we examined the (1) performance and suitability of selected prediction strategies, (2) contribution of relative position and stand-density measures in improving tree height (h) prediction values, and (3) effect of different subsampling designs to fill in missing h values in a new stand using a regional nonlinear model. Nonlinear mixed-effects models (NMEM) substantially improved the accuracy and precision of height prediction over the conventional nonlinear fixed-effects model (NFEM) that assumes the observations are independent, particularly when a few trees are subsampled for height. The predictive performance of a correction factor on a NFEM with relative position and stand-density measures was comparable to that of a NMEM when four or more trees were subsampled for height. When two or more heights were randomly subsampled, the NMEM efficiently explained the differences in the height–diameter relationship because of the variations in relative position of trees and stand density without having to incorporate them into the model. When only one height was subsampled, selecting the largest diameter tree in the stand would result in a lower predicted root mean square error (RMSE) than randomly selecting the height, regardless of the model form or fitting strategy used.
Bamboo forests play an important role in achieving the objectives of the United Nations program on Reducing Emission from Deforestation and Forest Degradation. We developed and validated a modeling system that simultaneously estimate aboveground biomass and its components for a common bamboo species (Bambusa procera A. Chev. and A. Camus) in tropical forests. Eighty-three bamboo culms were destructively sampled from seventeen 100 m2 sample plots located in different parts of the Central Highlands in Viet Nam to obtain total plant aboveground biomass (AGB) and its components. We examined the performance of weighted nonlinear models fit by maximum likelihood and weighted nonlinear seemingly unrelated regression fit by generalized least squares for predicting bamboo biomass. The simultaneous estimation of AGB and its components produced higher reliability than the models of components and total developed separately. With a large number of bamboo species, it may not be feasible to develop species- specific biomass models, hence genus-specific allometric models may be considered.
There are few allometric equations available for dipterocarp forests, despite the fact that this forest type covers extensive areas in tropical Southeast Asia. This study aims to develop a set of equations to estimate tree aboveground biomass (AGB) in dipterocarp forests in Vietnam and to validate and compare their predictive performance with allometric equations used for dipterocarps in Indonesia and pantropical areas. Diameter at breast height (DBH), total tree height (H), and wood density (WD) were used as input variables of the nonlinear weighted least square models. Akaike information criterion (AIC) and residual plots were used to select the best models; while percent bias, root mean square percentage error, and mean absolute percent error were used to compare their performance to published models. For mixed-species, the best equation was . When applied to a random independent validation dataset, the predicted values from the generic equations and the dipterocarp equations in Indonesia overestimated the AGB for different sites, indicating the need for region-specific equations. At the genus level, the selected equations were and for two genera, Dipterocarpus and Shorea, respectively, in Vietnam. Compared to the mixed-species equations, the genus-specific equations improved the accuracy of the AGB estimates. Additionally, the genus-specific equations showed no significant differences in predictive performance in different regions (e.g., Indonesia, Vietnam) of Southeast Asia.