Re: NM VI and NMV gave different results

----- Original Message ---- From: Leonid Gibiansky <[EMAIL PROTECTED]> To: Jian Xu <[EMAIL PROTECTED]> Cc: [email protected] Sent: Wednesday, April 9, 2008 1:17:57 PM Subject: Re: [NMusers] NM VI and NMV gave different results Hi Jian I did extensive testing of NM V and NM VI on all the datasets and models that I was able to find, and there was no difference in the results (OF and the fit; I have not looked on ETA p-values) for the FOCEI (for runs that converged). Unstable models sometimes stopped at different places, but the fit was similar or identical. If you had stable model in NM V, try to re-run it in NM V again: it could be that you changed the dataset somehow. For NM VI, try to use SLOW option at the estimation step, it should reproduce NM V results. It would help if you post NM V and NM6 output (error messages if any) Leonid -------------------------------------- Leonid Gibiansky, Ph.D. President, QuantPharm LLC web: www.quantpharm.com e-mail: LGibiansky at quantpharm.com tel: (301) 767 5566 Jian Xu wrote: > Dear all, > > I did a Pop PK a year ago in NM V, and got a stable model and reasonable > estimation. Recently, I tried to rerun the same code in NM VI (with > latest patches from Metrum) . The results are huge different. Some > p-values for Eta from NM VI are even significant. I really want to know > what the real difference is between NM VI and NM V. I appreciate any > help or reference. > > Thanks in advance. > > Jian > > Below is the control stream for your consideration: > > $PROB 2CM > > $DATA DATA.CSV IGNORE=# > > $INPUT ID AMT TIME CON=DV MDV FST > > $SUBROUTINES ADVAN4 TRANS4 > > $PK > > V2=THETA(1)*EXP(ETA(1)) ; VOLUME OF CENTRAL > V3=THETA(2)*EXP(ETA(2)) ; VOLUME OF PERIPHERAL > CL=THETA(3)*EXP(THETA(7)*ETA(1)) ; CLEARANCE > Q =THETA(4)*EXP(ETA(3)) ; DISTRIBUTION CLEARANCE > IF (FST.EQ.1) THEN > KA=THETA(5)*EXP(ETA(4)) ; ABSORPTION RATE CONSTANT FOR FASTED > ELSE > KA=THETA(6)*EXP(ETA(5)) ; ABSORPTION RATE CONSTANT FOR FED > END IF > > S2=V2 > > $ERROR > > Y=F*(1+EPS(1)) > > $THETA > (0,150) ;V2_F > (0,200) ;V3_F > (0,6) ;CL_F > (0,50) ;Q_F > (0,0.5) ;KA > (0,0.2) ;KA > (0,1) ;CL_V Correlation > > $OMEGA > 0.01 0.01 0.01 0.01 0.01 > $SIGMA > 0.04 > > $ESTIMATION SIG=3 MAX=5000 PRINT=10 METHOD=1 INTER NOABORT > $COVARIANCE > $TABLE ID TIME NOPRINT FILE=BASE > > > &&&&&&&&&& > &GO GATORS! & > &&&&&&&&&& > > __________________________________________________ > Do You Yahoo!? > Tired of spam? Yahoo! Mail has the best spam protection around > http://mail.yahoo.com > __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com