Comparison of different methods for fitting nonlinear mixed forest models and for making predictions

Nonlinear mixed models have become popular in forestry applications, and various methods have been proposed for fitting such models. However, it is difficult or even confusing to choose which method to use, and there is not much relevant information available, especially in the forestry context. The main objective of this study was to compare three commonly used methods for fitting nonlinear mixed models: the first-order, the first-order conditional expectation, and the adaptive Gaussian quadrature methods. Both the maximum likelihood and restricted maximum likelihood parameter estimation techniques were evaluated. Three types of data common in forestry were used for model fitting and model application. It was found that the first-order conditional expectation method provided more accurate and precise predictions for two models developed from data with more observations per subject. For one model developed on data with fewer observations per subject, the first-order method provided better model prediction...

[1]  Peter L. Bonate,et al.  Pharmacokinetic-Pharmacodynamic Modeling And Simulation , 2005 .

[2]  Machelle D. Wilson,et al.  Modeling response curves and testing treatment effects in repeated measures experiments: a multilevel nonlinear mixed-effects model approach , 2005 .

[3]  Ismo Strandén,et al.  Estimation of non-linear growth models by linearization: a simulation study using a Gompertz function , 2005, Genetics Selection Evolution.

[4]  Rafael Calama,et al.  Interregional nonlinear height-diameter model with random coefficients for stone pine in Spain , 2004 .

[5]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[6]  Hailemariam Temesgen,et al.  Analysis and comparison of nonlinear tree height prediction strategies for Douglas-fir forests , 2008 .

[7]  Robert L. Bailey,et al.  Nonlinear Mixed Effects Modeling for Slash Pine Dominant Height Growth Following Intensive Silvicultural Treatments , 2001 .

[8]  D. Bates,et al.  Approximations to the Log-Likelihood Function in the Nonlinear Mixed-Effects Model , 1995 .

[9]  LiYun,et al.  Evaluation of nonlinear equations for predicting diameter from tree height , 2012 .

[10]  Robert L. Bailey,et al.  Modeling and Prediction of Forest Growth Variables Based on Multilevel Nonlinear Mixed Models , 2001 .

[11]  Marie Davidian,et al.  Nonlinear models for repeated measurement data: An overview and update , 2003 .

[12]  Harold E. Burkhart,et al.  An Evaluation of Sampling Methods and Model Forms for Estimating Height-Diameter Relationships in Loblolly Pine Plantations , 1992, Forest Science.

[13]  Oliver Kuss,et al.  How to use SAS ® for Logistic Regression with Correlated Data , 2002 .

[14]  S. Huang,et al.  Assessing the goodness of fit of forest models estimated by nonlinear mixed-model methods , 2009 .

[15]  A. Robinson,et al.  Imputing missing height measures using a mixed-effects modeling strategy , 2004 .

[16]  Marie Davidian,et al.  Consequences of misspecifying assumptions in nonlinear mixed effects models , 2000 .

[17]  Guillermo Trincado,et al.  A multilevel individual tree basal area increment model for aspen in boreal mixedwood stands , 2009 .

[18]  D. Bates,et al.  Nonlinear mixed effects models for repeated measures data. , 1990, Biometrics.

[19]  Mo Zheng.  Fitting Linear and Nonlinear Growth Curve Models Using PROC NLMIXED , 2010 .

[20]  E. Vonesh,et al.  Non-linear models for the analysis of longitudinal data. , 1992, Statistics in medicine.

[21]  Russell D. Wolfinger,et al.  Two Taylor-series approximation methods for nonlinear mixed models , 1997 .

[22]  S. Meng,et al.  Improved calibration of nonlinear mixed-effects models demonstrated on a height growth function. , 2009 .

[23]  John F. Monahan,et al.  Monte Carlo Comparison of ANOVA, MIVQUE, REML, and ML Estimators of Variance Components , 1984 .

[24]  Robert E. Weiss,et al.  Modeling Longitudinal Data , 2005 .

[25]  Harold E. Burkhart,et al.  A Generalized Approach for Modeling and Localizing Stem Profile Curves , 2006, Forest Science.

[26]  Guillermo Trincado,et al.  Nonlinear mixed-effects modeling of variable-exponent taper equations for lodgepole pine in Alberta, Canada , 2009, European Journal of Forest Research.

[27]  Timothy G. Gregoire,et al.  Nonlinear Mixed-Effects Modeling of Cumulative Bole Volume With Spatially Correlated Within-Tree Data , 1996 .

[28]  E. Vonesh,et al.  Linear and Nonlinear Models for the Analysis of Repeated Measurements , 1996 .

[29]  J. Ware,et al.  Applied Longitudinal Analysis , 2004 .

[30]  L. Skovgaard.  NONLINEAR MODELS FOR REPEATED MEASUREMENT DATA. , 1996 .

[31]  S. Huang,et al.  Evaluation of population-averaged and subject-specific approaches for modeling the dominant or codominant height of lodgepole pine trees , 2009 .

[32]  Barry Gardiner,et al.  Variable-top merchantable volume equations for Scots pine (Pinus sylvestris) and Sitka spruce (Picea sitchensis (Bong.) Carr.) in Northern Britain , 2012 .

[33]  Lewis B. Sheiner,et al.  Evaluation of methods for estimating population pharmacokinetic parameters. I. Michaelis-menten model: Routine clinical pharmacokinetic data , 1980, Journal of Pharmacokinetics and Biopharmaceutics.