Problems when implementing M3 for BLQs

2 messages 2 people Latest: Aug 17, 2022

Re: Problems when implementing M3 for BLQs

From: Leonid Gibiansky Date: August 17, 2022 technical
Hi Roeland, ERROR_TERM in Y and in DUM should differ: one should include EPS() while the other should use only SD. F_FLAG should be set to zero outside M3 block. What is DRUG.EQ.1, is it parent? there should be a separate M3 part for metabolite. It would be easier to debug if you would provide the full code of the ERROR block Thank you, Leonid
Quoted reply history
On 8/16/2022 6:53 AM, Roeland Wasmann wrote: > Hi all, > > I have a problem while trying to implement the m3 method in my model with a parent and metabolite. I could really use some help. > > About 50% of the parent data is BLQ while the metabolite only has a couple of BLQs. When modeling them separately there is no issue but when I have them both in one model, I get the following error: > > /Recompiling certain components. > starting wait > ending wait > Exiting lpreddo > Compiling FSUBS > FSUBS.f90:1322:6:/ > > / B000056=PHI(DUM)/ > > / 1/ > > /Error: Unclassifiable statement at (1) > Building NONMEM Executable > gfortran: error: fsubs.o: No such file or directory > No nonmem execution. > It seems like Fortran compilation by NONMEM's nmfe script failed. Cannot > start NONMEM. > > Go to the NM_run1 subdirectory and run psn.mod with NONMEM's nmfe script to diagnose the problem./ > > The NM_run1 folder does not provide any answers… The M3 part of the error block looks like this: > > Y = IPRED + ERROR_TERM > > *DUM = (LLOQ-IPRED) / ERROR_TERM > CUMD = PHI(DUM) > IF (ICALL.NE.4.AND.BLQ.GT.0.AND.DRUG.EQ.1) THEN > F_FLAG = 1 > Y = CUMD > MDVRES = 1 > ENDIF* > > When I put the “DUM” part within the PHI (i.e., “*CUMD = PHI(((LLOQ-IPRED) / ERROR_TERM))*” the error changes slightly: > > /Compiling FSUBS > FSUBS.f90:1324:6:/ > > / B000058=PHI(B000056)/ > > When I remove the “*.AND.DRUG.EQ.1*” from IF I get: > > /FSUBS.f90:1355:6:/ > > / B000061=PHI(DUM)/ > > Please let me know what info I should provide to better diagnose this problem challenge. > > Thanks in advance for helping out! > > Cheers Roeland > > ------------------------------------------ > > *Roeland Wasmann, PharmD, PhD* > > Division of Clinical Pharmacology > > Department of Medicine > > University of Cape Town > > K45 Old Main Building > > Groote Schuur Hospital > > Observatory, Cape Town > > 7925 South Africa > > phone: +27 21 650 4861 > > Disclaimer - University of Cape Town This email is subject to UCT policies and email disclaimer published on our website at http://www.uct.ac.za/main/email-disclaimer or obtainable from +27 21 650 9111. If this email is not related to the business of UCT, it is sent by the sender in an individual capacity. Please report security incidents or abuse via https://csirt.uct.ac.za/page/report-an-incident.php .
Make sure also that you do not have a data item or variable called PHI. It may interfere with PHI() as a function. Robert J. Bauer, Ph.D. Senior Director Pharmacometrics R&D ICON Early Phase 731 Arbor way, suite 100 Blue Bell, PA 19422 Office: (215) 616-6428 Mobile: (925) 286-0769 [email protected]<mailto:[email protected]> http://www.iconplc.com/
Quoted reply history
From: [email protected] <[email protected]> On Behalf Of Leonid Gibiansky Sent: Wednesday, August 17, 2022 1:48 PM To: Roeland Wasmann <[email protected]>; [email protected] Subject: [EXTERNAL] Re: [NMusers] Problems when implementing M3 for BLQs Hi Roeland, ERROR_TERM in Y and in DUM should differ: one should include EPS() while the other should use only SD. F_FLAG should be set to zero outside M3 block. What is DRUG.EQ.1, is it parent? there should be a separate M3 part for metabolite. It would be easier to debug if you would provide the full code of the ERROR block Thank you, Leonid On 8/16/2022 6:53 AM, Roeland Wasmann wrote: > Hi all, > > I have a problem while trying to implement the m3 method in my model > with a parent and metabolite. I could really use some help. > > About 50% of the parent data is BLQ while the metabolite only has a > couple of BLQs. When modeling them separately there is no issue but when > I have them both in one model, I get the following error: > > /Recompiling certain components. > starting wait > ending wait > Exiting lpreddo > Compiling FSUBS > FSUBS.f90:1322:6:/ > > / B000056=PHI(DUM)/ > > / 1/ > > /Error: Unclassifiable statement at (1) > Building NONMEM Executable > gfortran: error: fsubs.o: No such file or directory > No nonmem execution. > It seems like Fortran compilation by NONMEM's nmfe script failed. Cannot > start NONMEM. > Go to the NM_run1 subdirectory and run psn.mod with NONMEM's nmfe script > to diagnose the problem./ > > The NM_run1 folder does not provide any answers… The M3 part of the > error block looks like this: > > Y = IPRED + ERROR_TERM > > *DUM = (LLOQ-IPRED) / ERROR_TERM > CUMD = PHI(DUM) > IF (ICALL.NE.4.AND.BLQ.GT.0.AND.DRUG.EQ.1) THEN > F_FLAG = 1 > Y = CUMD > MDVRES = 1 > ENDIF* > > When I put the “DUM” part within the PHI (i.e., “*CUMD = > PHI(((LLOQ-IPRED) / ERROR_TERM))*” the error changes slightly: > > /Compiling FSUBS > FSUBS.f90:1324:6:/ > > / B000058=PHI(B000056)/ > > When I remove the “*.AND.DRUG.EQ.1*” from IF I get: > > /FSUBS.f90:1355:6:/ > > / B000061=PHI(DUM)/ > > Please let me know what info I should provide to better diagnose this > problem challenge. > > Thanks in advance for helping out! > > Cheers Roeland > > ------------------------------------------ > > *Roeland Wasmann, PharmD, PhD* > > Division of Clinical Pharmacology > > Department of Medicine > > University of Cape Town > > K45 Old Main Building > > Groote Schuur Hospital > > Observatory, Cape Town > > 7925 South Africa > > phone: +27 21 650 4861 > > Disclaimer - University of Cape Town This email is subject to UCT > policies and email disclaimer published on our website at > http://www.uct.ac.za/main/email-disclaimer > or obtainable from +27 21 650 > 9111. If this email is not related to the business of UCT, it is sent by > the sender in an individual capacity. Please report security incidents > or abuse via > https://csirt.uct.ac.za/page/report-an-incident.php.