Dear NonMEM users,
I am constructing a transduction model to capture disease onset time-delay
using NonMEM. Because
I need to use too mamy compartments (means too much ODE) to handle this delay,
I want to use DoLoop
to save a little room. But I can not find a way to use it in $DES block. The
block and error message are
listed as below,
$DES
I=4
DOWHILE(I.LT.22) ;(I also tried Do I=5,22, same as DOWHILE)
im1=I-1
DADT(I) = KT*(A(im1)-A(I))
I=I+1
ENDDO
Error message:
DOWHILE(I.LT.22)
X
THE CHARACTERS IN ERROR ARE: DO
438 THIS WORD IS NOT APPROPRIATE FOR THIS BLOCK OF ABBREVIATED CODE.
I am very appreciated if anyone can tell me how to incorporate DoLoop into $DES
block.
Thanks a lot!
Best regards,
Liu, Dongyang, PhD, Postdoc Fellow
Department of Pharmaceutical sciences,
State University of New York at Buffalo.
Tel(o):01-716-645-4840,
Cell: 01-716-908-6644,
About DoLoop in NONMEM
Dear Liu,
In nonmem you can not define differential equations through a loop; each
has to be explicitly written out.
However, I notice that you have the same rate of transit among all of
your (transit) compartments.
Therefore, an analytical solution would be the most efficient way of
solving your problem.
Please see the reference below.
Kind regards
Jakob
Analytical solution:
Implementation of a transit compartment model for describing drug
absorption in pharmacokinetic studies.
Savic RM, Jonker DM, Kerbusch T, Karlsson MO.
J Pharmacokinet Pharmacodyn. 2007 Oct;34(5):711-26. Epub 2007 Jul 26.
You may be familiar with this already?
Transit compartments versus gamma distribution function to model signal
transduction processes in pharmacodynamics.
Sun YN, Jusko WJ.
J Pharm Sci. 1998 Jun;87(6):732-7.
Quoted reply history
________________________________________
From: [email protected] [mailto:[email protected]]
On Behalf Of Liu Dongyang
Sent: 25 March 2011 04:17
To: nmusers
Subject: [NMusers] About DoLoop in NONMEM
Dear NonMEM users,
I am constructing a transduction model to capture disease onset
time-delay using NonMEM. Because
I need to use too mamy compartments (means too much ODE) to handle this
delay, I want to use DoLoop
to save a little room. But I can not find a way to use it in $DES
block. The block and error message are
listed as below,
$DES
I=4
DOWHILE(I.LT.22) ;(I also tried Do I=5,22, same as DOWHILE)
im1=I-1
DADT(I) = KT*(A(im1)-A(I))
I=I+1
ENDDO
Error message:
DOWHILE(I.LT.22)
X
THE CHARACTERS IN ERROR ARE: DO
438 THIS WORD IS NOT APPROPRIATE FOR THIS BLOCK OF ABBREVIATED CODE.
I am very appreciated if anyone can tell me how to incorporate DoLoop
into $DES block.
Thanks a lot!
Best regards,
Liu, Dongyang, PhD, Postdoc Fellow
Department of Pharmaceutical sciences,
State University of New York at Buffalo.
Tel(o):01-716-645-4840,
Cell: 01-716-908-6644,