What does nonmem do after iteration stops?
From: "Eleveld, DJ" d.j.eleveld@anest.umcg.nl
Subject: [NMusers] What does nonmem do after iteration stops?
Date: Fri, June 3, 2005 6:28 am
Hi Everyone,
I've noticed something curious:
I use the PRINT=1 option in the $ESTIMATION to see the progress of the model fitting. With a fitting that
I am doing (PD mdoel fitting involving about 3000 data points) the iterations stop after about 20-30 minutes
and then nothing is printed on the screen for several hours. Apparently nonmem is calculating something
after convergence that takes a very long time. The problem does have 2 tables but I cant really imagine that
these would require so much computing time.
What does nonmem do after finishing iterations that could possibly require so much computation?
Is there a way for me to avoid this computation and still see the results of the estimation?
If its useful my control file is at the end of this message.
Thanks very much,
Doug Eleveld
------------------------------------
$PROB Potentiation fitting
$DATA potpd_.prn
$INPUT ID TIME CPLA DV MDV AMT RATE
$SUBROUTINES ADVAN9 TOL=6
$ABBREVIATED COMRES=1
$MODEL NCOMPARTMENTS=2
NPARAMETERS=6
COMP(POTENT NOOFF)
COMP(EFFECT NOOFF NODOSE)
$PK
CALLFL=0
KEO=THETA(1)*EXP(ETA(1))
EC50=THETA(2)*EXP(ETA(2))
GAMM=THETA(3)*EXP(ETA(3))
POTR=THETA(4)*EXP(ETA(4))
POTK=THETA(5)*EXP(ETA(5))
SCAL=THETA(6)+ETA(6)
IF (NEWIND.EQ.0) COM(1)=POTR
F1=ABS(COM(1))
$DES
DADT(1)=-POTK*A(1) ; Decay potentiation
DADT(2)=(CPLA-A(2))*KEO ; Effect compartment conc
$ERROR
DPOT=A(1) ; The degree of potentiation
CEFF=A(2)
DPD1=CEFF**GAMM ; Degree of NMB
NMB=DPD1/(DPD1+(EC50**GAMM))
Y=SCAL*(1+DPOT)*(1-NMB)+ERR(1) ; Twitch prediction
COM(1)=POTR*(1-NMB)
$THETA (0,0.13)(0,1460)(0,3.45)
(0,0.0008)(0,0.146)(0,101)
$OMEGA 0.1 0.1 0.1 0.1 0.1 10
$SIGMA 10
$ESTIMATION MAX=9999 SIG=6 PRINT=1
$TABLE TIME KEO EC50 GAMM POTR POTK SCAL
NOHEADER NOPRINT FILE=potent2.txt
$TABLE TIME CEFF DPOT NMB
NOHEADER NOPRINT FILE=potent3.txt