Abstract
An EM algorithm was used to analyse data arising from non-linear mixed-effects models. The fixed parameters were determined by maximum likelihood using simplex minimization, and the random effects were estimated using the EM algorithm after linearization with respect to the random effects. Applications to a simple linear model and population pharmacokinetics are described. The use of posterior parameter estimates to investigate covariate relationships is briefly described. The implementation of the estimation-maximization (EM) algorithm described here has proved in practice to be robust but slow. We intend to use a Newton-Raphson minimization routine in place of the simplex method to hasten convergence. The alternative linearization of the non-linear mixed effects model suggested by Lindstrom and Bates ( Biometrics 46 (1990) 673–687) is much more unstable than the usual linearization, especially during the initial iterations. In the case of indomethacin the two linearizations produced very similar results. The individual posterior parameter estimates provided by the program are very useful for the detection of covariate relationships in population pharmacokinetic studies. In addition, the posterior means can be used in the estimation of pharmacokinetic-pharmacodynamic relationships from sparse pharmacokinetic data where individual modelling is impossible.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.