![asreml-r users variance-covariance matrix asreml-r users variance-covariance matrix](https://media.springernature.com/lw685/springer-static/image/art%3A10.1007%2Fs13253-016-0267-0/MediaObjects/13253_2016_267_Fig1_HTML.gif)
In addition, an automated preconditioning routine is made available as a part of the software package Perl-speaks-NONMEM (PsN) version 4.4 and above ( 3).Ĭomputational instability related to the R-matrix can influence pharmacometric analyses in two ways. Our hypothesis is that preconditioning will reduce the failed variance-covariance matrix computations when the R-matrix is fundamentally positive definite, and correctly indicate a fundamentally singular R-matrix in case of the existence of non-identifiable model parameters. To test this preconditioning method, we have conducted numerical experiments using published non-linear mixed effects models (and data) from applications in pharmacometrics. This approach should reduce the influence of computational issues and reduce the chance of the R-matrix being non-positive definite and also give an indication of when the R-matrix is fundamentally singular. Inspired by this technique, we reparameterise the model with a linear combination of the model parameters so that the Hessian of the −2ln(likelihood) (R-matrix) of the reparameterised model becomes close to an identity matrix. Preconditioning is a widely used technique to increase the computational stability for numerically solving large sparse systems of linear equations ( 1, 2). In this paper, we propose a preconditioning method for non-linear mixed effects models to increase the computational stability of the variance-covariance matrix. The pharmacometric analysis is usually based on the assumption that the parameter and uncertainty estimates of the non-linear model are estimated correctly by a numerical method however, these estimates will unavoidably be influenced by the numerical stability of the computational algorithm when the numerical method is based on finite precision computations. As a result, pharmacometric analysis based on non-linear mixed effects models, also known as the population approach, has become an essential step in drug development.
![asreml-r users variance-covariance matrix asreml-r users variance-covariance matrix](https://media.springernature.com/original/springer-static/image/chp%3A10.1007%2F978-3-319-55177-7_3/MediaObjects/325391_1_En_3_Figc_HTML.gif)
Non-linear mixed effects models have been shown to be an effective tool for the analysis of clinical trial data.