PAGE. Abstracts of the Annual Meeting of the Population Approach Group in Europe.
PAGE 24 (2015) Abstr 3586 [www.page-meeting.org/?abstract=3586]
Click to open
Oral: Stuart Beal Methodology Session
Yasunori Aoki, Rikard Nordgren, and Andrew C. Hooker
Department of Pharmaceutical Biosciences, Uppsala University, Uppsala, Sweden
Objectives: As the importance of pharmacometric analysis increases, more and more complex models are introduced and computational stability starts to become a bottleneck. We can observe such instability, for example, in the computation of the covariance matrix with the error message "R MATRIX ALGORITHMICALLY NON-POSITIVESEMIDEFINITE” in NONMEM. In this work, we present a preconditioning method for nonlinear mixed effect models to increase the computational stability of the covariance matrix computation.
Methods: Roughly speaking, the method re-parameterizes the model with a linear combination of the original model parameters so that the R-matrix of the re-parameterized model becomes close to an identity matrix. This approach will reduce the chance of the R-matrix being non-positive semi-definite due to computational instability and gives a clear indication when the R-matrix is fundamentally non-positive semi-definite (e.g., if the model is not identifiable). Based on the re-parameterized model results, the parameters and covariance matrix in the original parameterization can be calculated.
Results: We have conducted a number of stochastic simulation and estimation experiments, using three published models [1,2,3] and NONMEM. We have simulated various datasets and then estimated both parameters and covariance matrices using the published parameterizations. In these studies there were 85/200, 30/50, and 69/200 of cases where the covariance step failed. However, after preconditioning the covariance matrices were successfully computed for all these cases. In addition, to illustrate the danger of computational instability, we have conducted a similar test using an unidentifiable model. Surprisingly, covariance matrices could be computed for 48/100 cases with reasonable relative standard errors (RSE). However after preconditioning, the RSE typically grew by a few orders of magnitude (e.g., 47% to 1243%) clearly indicating identifiability problems with the parameters.
Conclusions: Computational instability can potentially influence pharmacometric analyses and we propose a preconditioning method to reduce that instability and increase the chances of getting a covariance matrix in NONMEM (if the model parameters are identifiable). The method is automated and made available as a part of PsN [4,5]. Computational instability can also influence the parameter estimates and an investigation of this correlation using the preconditioning method is presented in a separate abstract .
Acknowledgement: This work was supported by the DDMoRe (www.ddmore.eu) project.