×

Maximum likelihood estimation of nonlinear mixed-effects models with crossed random effects by combining first-order conditional linearization and sequential quadratic programming. (English) Zbl 1419.92012

Summary: Nonlinear mixed-effects (NLME) models have become popular in various disciplines over the past several decades. However, the existing methods for parameter estimation implemented in standard statistical packages such as SAS and R/S-Plus are generally limited to single- or multi-level NLME models that only allow nested random effects and are unable to cope with crossed random effects within the framework of NLME modeling. In this study, we propose a general formulation of NLME models that can accommodate both nested and crossed random effects, and then develop a computational algorithm for parameter estimation based on normal assumptions. The maximum likelihood estimation is carried out using the first-order conditional expansion (FOCE) for NLME model linearization and sequential quadratic programming (SQP) for computational optimization while ensuring positive-definiteness of the estimated variance-covariance matrices of both random effects and error terms. The FOCE-SQP algorithm is evaluated using the height and diameter data measured on trees from Korean larch (L. olgensis var. Changpaiensis) experimental plots as well as simulation studies. We show that the FOCE-SQP method converges fast with high accuracy. Applications of the general formulation of NLME models are illustrated with an analysis of the Korean larch data.

MSC:

92C80 Plant biology
62P10 Applications of statistics to biology and medical sciences; meta analysis
90C20 Quadratic programming

Software:

SASmixed; MEMSS; S-PLUS
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Sheiner, L. B. and Beal, S. L., Evaluation of methods for estimating population pharmacokinetic parameters. I. michaelis-menton model routine clinical pharmacokinetic data, J. Pharmacokinet. Biopharm.8 (1980) 553-571.
[2] Davidian, M. and Giltinan, D. M., Nonlinear Models for Repeated Measurement Data (Chapman & Hall, New York, 1995).
[3] Vonesh, E. F. and Chinchilli, V. M., Linear and Nonlinear Models for the Analysis of Repeated Measurements, (Marcel Dekker, New York, NY, 1997). · Zbl 0893.62077
[4] Pinheiro, J. C. and Bates, D. M., Mixed-effects models in S and S-PLUS (Springer, New York, 2000). · Zbl 0953.62065
[5] Davidian, M. and Giltiman, D. M., Nonlinear models for repeated measurement data: An overview and update, J. Agric. Biol. Environ. Stat.8 (2003) 387-419.
[6] Hall, D. B. and Clutter, M., Multivariate multilevel nonlinear mixed effects models for timber yield predictions, Biometrics60 (2004) 16-24. · Zbl 1130.62391
[7] Panhard, X. and Samson, A., Extension of the SAEM algorithm for nonlinear mixed models with 2 levels of random effects, Biostatistics10(1) (2009) 121-135. · Zbl 1437.62573
[8] Fu, L., Wang, M., Lei, Y. and Tang, S., Parameter estimation of two-level nonlinear mixed effects models using first order conditional linearization and the EM algorithm, Comput. Statist. Data Anal.69 (2014) 173-183. · Zbl 1471.62067
[9] Laird, N. M. and Ware, J. H., Random-effects models for longitudinal data, Biometrics38 (1982) 963-974. · Zbl 0512.62107
[10] Littell, R. C., Milliken, G. A., Stroup, W. W., Wolfinger, R. D. and Schabenberber, O., SAS for Mixed Models, 2nd edn. (SAS Institute Inc, Cary, NC, 2006).
[11] Lee, Y., Nelder, J. A. and Pawitan, Y., Generalized Linear Models with Random Effects: Unified Analysis Via H-Likelihood (Chapman & Hall/crcBoca Raton Fl, 2006). · Zbl 1110.62092
[12] Raudenbush, S. W., A crossed random effects model for unbalanced data with applications in cross sectional and longitudinal research, J. Educ. Stat.18(4) (1993) 321-349.
[13] Rasbash, J. and Goldstein, H., Efficient analysis of mixed hierarchical and cross classified random structures using a multilevel model, J. Educ. Behav. Stat.19(4) (1994) 337-350.
[14] Clayton, D. and Rasbash, J., Estimation in large crossed random-effect models by data augmentation, J. Roy. Statist. Soc.162(4) (1999) 425-436.
[15] Goldstein, H., Multilevel Statistical Models, 4th edn. (Wiley, Hoboken, N. J, 2011). · Zbl 1274.62006
[16] Sun, S. and Pan, W., A methodological review of statistical methods for handling multilevel non-nested longitudinal data in educational research, Int. J. Res. Meth. Educ.37(3) (2014) 285-308.
[17] Pinheiro, J. C. and Bates, D. M., Approximations to the log-likelihood function in the nonlinear mixed effects model, J. Comput. Graph. Statist.4 (1995) 12-35.
[18] Wolfinger, R. D. and Lin, X., Two taylor-series approximation methods for nonlinear mixed models, Comput. Statist. Data Anal.25 (1997) 465-490. · Zbl 0900.65409
[19] Lindstrom, M. J. and Bates, D. M., Nonlinear mixed effects models for repeated measures data, Biometrics46 (1990) 673-687.
[20] Fletcher, R., The sequential quadratic programming method, in Nonlinear Optimization, Lecture Notes in Mathematics, Vol. 1989 (Springer, Berlin, 2010), pp. 165-214. · Zbl 1192.90002
[21] Henderson, C. R., Sire evaluation and genetic trends, J. Anim. Sci.1973 (1973) 10-41.
[22] Goldstein, H., Multilevel covariance component models, Biometrika74(2) (1987) 430-431. · Zbl 0619.62064
[23] Harville, D. A., Maximum likelihood approaches to variance component estimation and to related problems, J. Am. Stat. Assoc.72 (1977) 320-340. · Zbl 0373.62040
[24] Wolfinger, R., Tobias, R. and Sall, J., Computing gaussian likelihoods and their derivatives for general linear mixed models, SIAM. J. Sci. Comput.15(6) (1994) 1294-1310. · Zbl 0805.62069
[25] Lindstrom, M. J. and Bates, D. M., Newton-Raphson and EM algorithms for linear mixed-effects models for repeated-measures data, J. Am. Stat. Assoc.83 (1988) 1014-1022. · Zbl 0671.65119
[26] Hestenes, M. R., Optimization Theory (Wiley, New York, 1975). · Zbl 0327.90015
[27] Tang, S. Z., Lang, K. J. and Li, H. K., Statistics and Computation of Biomathematical Models (ForStat Course) (Science Press, Beijing, China (in Chinese), 2008).
[28] Yang, Y. and Huang, S., Comparison of different methods for fitting nonlinear mixed forest models and for making predictions, Can. J. For. Res.41(8) (2011) 1671-1686.
[29] Swallow, W. H. and Monahan, J. F., Monte Carlo comparison of ANOVA, MIVQUE, REML, and ML estimators of variance components, Technometrics26 (1984) 47-57. · Zbl 0548.62051
[30] Yang, Y., Huang, S., Trincado, G. and Meng, S. X.. Nonlinear mixed-effects modelling of variable-exponent taper equations for lodgepole pine in Alberta, Canada, Eur. J. For. Res.128 (2009) 415-429.
[31] Pinheiro, J. C. and Bates, D. M., Unconstrained parametrizations for variance-covariance matrices, Stat. Comput.6 (1996) 289-296.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.