Nonlinear mixed effects models [Software]

posted by Helmut Homepage – Vienna, Austria, 2018-06-14 15:03 (2114 d 21:16 ago) – Posting: # 18907
Views: 2,802

Hi libaiyi,

❝ I can't get the same results from winnonlin and nonmem.


Well, we are dealing with nonlinear mixed effects models with all the bells and whistles. The software tries to find the optimum in an n+1 dimensional parameter-space. NONMEM and Phoenix/NLME employ maximum likelihood, whereas WinNonlin least squares. See this post for details and the presentation given at the PAGE meeting 2005 about comparing PopPK software in a blinded manner. Also interesting: Bauer et al.1, Colby and Bair2. The latter deals with Phoenix.

❝ Is there any difference in algorithm in two softwares?


Well, in both there are more than one. It would be very difficult to get exactly the same results even if you use the same algorithm, initial estimates, parameter constraints, convergence criteria, max. number of iterations, :blahblah:
We can only hope that an algo converges to the global – and not a local – optimum.
That’s the reason why in BE only NCA is acceptable.

Or does it exist some mistakes in code?


Always possible (see this presentation). The software industry generally considers one defect / 2,000 lines of code ‘stable’. Oh dear! Package PowerTOST has currently ~10,200 LOC – I hope, we did better. Of course, higher standards are used in certain fields (see what I wrote about white-box validation in this post).


  1. Bauer RJ, Guzy S, Ng C. Survey of Population Analysis Methods and Software for Complex Pharmacokinetic and Pharmacodynamic Models with Examples. AAPS J. 2007;9(1):E60–83. doi:10.1208/aapsj0901007. [image] PMC Free Full text.
  2. Colby E, Bair E. Cross-Validation for Nonlinear Mixed Effects Models. J Pharmacokinet Pharmacodyn. 2013;40(2):243–52. doi:10.1007/s10928-013-9313-5. [image] PMC Free Full text. [image] supplementary material.

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
98 visitors (0 registered, 98 guests [including 5 identified bots]).
Forum time: 11:20 CET (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5