An interesting undocumented feature of PROTECT. If you have EXP(0), e.g., KA = TVKA *EXP(0)
The fortran translation in this case apparently is: KA=TVKA*PEXP(*)B000092 Which then throws a syntax error. Doesn't happen without $ABBR PROTECT Not sure if that counts as a bug, but does prohibit using PROTECT in automated methods to run models with/without BSV, or at least makes it a little more difficult. Mark Sale M.D. Vice President Integrated Drug Development mark.s...@certara.com Remote-Forestville CA Office Hours 9 AM - 5 PM Eastern Time +1 302-516-1684 www.certara.com -----Original Message----- From: owner-nmus...@globomaxnm.com <owner-nmus...@globomaxnm.com> On Behalf Of Leonid Gibiansky Sent: Friday, January 26, 2024 7:42 AM To: Jeroen Elassaiss-Schaap (PD-value B.V.) <jer...@pd-value.com>; Elashkar, Omar I. <omar.elash...@ufl.edu>; nmusers@globomaxnm.com Subject: Re: [NMusers] NONMEM Error on SAEM but not FOCEI CAUTION: This email originated from outside of Certara. Do not click links or open attachments unless you recognize the sender and know the content is safe. This part is easy to solve: if you do it in log space, Q will always be positive, no PROTECT is needed (and I am not sure it is relevant in the case when Q is going to be negative). On 1/26/2024 10:27 AM, Jeroen Elassaiss-Schaap (PD-value B.V.) wrote: > Dear Omar, > > > What you describe is not that strange. It can happen because SAEM does > random sampling whereas FOCE-I uses a directed search. Therefore > extreme values can occur with SAEM. Please refer to the section I60 of > user guide for v7+ "Stable Routines for Estimation Methods and > Automated Protection Against Floating Point Exceptions" for more > detail, but you could start with including > > $ABBR PROTECT > > which probably is adequate in your case. > > > Hope this helps, > > Jeroen > > > https://nam11.safelinks.protection.outlook.com/?url=http%3A%2F%2Fpd-va > lue.com%2F&data=05%7C02%7CMark.sale%40Certara.com%7Cecb21da5ca6e4a50d8 > b308dc1e865d9a%7C7287abd30220456e98514352bae208c9%7C1%7C0%7C6384188096 > 75762447%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiL > CJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=4uyLVZH1Zx4EWEQws7xn4 > %2FKR%2BeFIro5ul10iI6LMUM4%3D&reserved=0 > jer...@pd-value.com > @PD_value > +31 6 23118438 > -- More value out of your data! > > On 26-01-2024 16:01, Elashkar, Omar I. wrote: >> Hello NMusers, >> >> I have created a model that fit ok with FOCE+I. When I try to fit >> with SAEM, I get the following error: >> >> """ >> ERROR IN TRANS4 ROUTINE: Q IS ZERO >> Elapsed estimation time in seconds: 0.00 >> INDIVIDUAL OBJECTIVE FUNCTION VALUES ARE ALL ZERO. PROBLEM ENDED """ >> >> The data is very small, but the model seems acceptable from bootstrap >> on FOCE+I. Any ideas are appreciated! >> >> Thank you, >> Omar >> >> >> >> *Omar Elashkar, *B.Pharm, MSc >> >> Graduate Assistant >> >> University of Florida | College of Pharmacy >> >> >> > This message (including any attachments) may contain confidential, proprietary, privileged and/or private information. The information is intended to be for the use of the individual or entity designated above. If you are not the intended recipient of this message, please notify the sender immediately, and delete the message and any attachments. Any disclosure, reproduction, distribution or other use of this message or any attachments by an individual or entity other than the intended recipient is prohibited.
