Re: ADVAN2 and ADVAN6: biased results

From: Nick Holford - n.holford@auckland.ac.nz Subject: Re: [NMusers] ADVAN2 and ADVAN6: biased results Date: 2/14/2004 2:15 AM Hi, When I tried your ADVAN2 code it terminated by running out of function evaluations. You can get it to converge and run the $COV step if you allow more function evalutions and tweak your initial estimates a bit. I can't see any problem when I look at graphs of your fit. THETA: K KA V ETA: PPVK PPVKA PPVV ERR: SD advan2.lst 101.072 eval=122 sig=+3.4 sub=12 obs=132 CCIL=NNNN NV1.1 PIV1.1 THETA = 0.0809 2.82 0.481 ETASD = 0.0181384 2.49199 0.070214 ETAR12 = 0.241 ETAR23 = 0.159 -0.213 ERRSD = 0.618061 THETA:se% = 11.4 25.0 4.8 OMEGA:se% = 81.8 85.3 37.9 SIGMA:se% = 26.0 $PROBLEM FIRST ADVAN2 -- POPULATION DATA $INPUT ID TIME CP=DV DOSE=AMT MDV $DATA A2IN2.dat IGNORE # $EST MAXEVAL=9990 SIG=3 PRINT=0 POSTHOC $COV $THETA (.001,0.0806,5) ;K (.001,2.81,10) ;KA (.001,0.481,20) ;V $OMEGA BLOCK(3) 0.00032 ;PPVK 0.0104 6.2 ;PPVKA 0.00021 -0.0369 0.00492 ;PPVV $SIGMA 0.382 ;SD $SUBROUTINES ADVAN2 $PK CALLFL=1 K=THETA(1)+ETA(1) KA=THETA(2)+ETA(2) V=THETA(3)+ETA(3) S2=V $ERROR CONC=F Y=CONC+ERR(1) $TABLE ID TIME DOSE Y ONEHEADER NOPRINT FILE=advan2.fit $SCAT (DV PRED) VS TIME BY ID -- Nick Holford, Dept Pharmacology & Clinical Pharmacology University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand email:n.holford@auckland.ac.nz tel:+64(9)373-7599x86730 fax:373-7556 http://www.health.auckland.ac.nz/pharmacology/staff/nholford/ _______________________________________________________