search

JohanWBorg / pharmpy

A library and toolkit for pharmacometrics
https://pharmpy.github.io
GNU General Public License v3.0
0 stars 0 forks source link

IF(ETA(1).GT.10) EXIT 1 2 does no seem to be converted? #43

Open JohanWBorg opened 1 year ago

JohanWBorg commented 1 year ago

IF(ETA(1).GT.10) EXIT 1 2 IF(ETA(2).GT.10) EXIT 1 3 IF(ETA(1).LT.-10) EXIT 1 4 IF(ETA(2).LT.-10) EXIT 1 5

not converted whatsoever.

JohanWBorg commented 1 year ago

https://github.com/nlmixrdevelopment/nlmixr/issues/252 Has information regarding this feature. Specifically, it states "Large ETAs and Omegas are handled different from NONMEM. In SAEM, you cannot exit and reset the system.

In nlmixr FOCEi, there are two different types of automatic reset if all of your ETAs or some of your ETAs are becoming too large which are controlled by resetThetaP and resetEtaP in foceiControl. You can bound Omega values in FOCEi with diagOmegaBoundUpper and diagOmegaBoundLower. Other than that, there is currently no way to control etas getting too high by an "exit" statement."

JohanWBorg commented 1 year ago

These statements are not stored in pharmpy in any place I have found. They are only stored in the model code from what I have found and hence I have a hard time translating them

JohanWBorg commented 1 year ago

Per talks with supervisor, this will not be implemented