Re: NONMEM questions
From: Ping Ji ping.ji@bms.com
Subject: Re: [NMusers] NONMEM questions
Date: Mon, 21 Mar 2005 09:49:04 -0500
Dear all,
Thanks a lot for the quick reply. I truly appreciate it. To start with
MAD or POC data only, and use the following control file for both
dataset, I can get both minimization successfully.
> $PROB MAD
> $INPUT C ID AMT RATE TIME DV EVID
> $DATA full-data-set-MAD.csv IGNORE=C
> $SUBROUTINE ADVAN4 TRANS4
> $PK
>
> CL=THETA(1)*EXP(ETA(1))
> D1=THETA(2)
>
> V2=THETA(3)*EXP(ETA(2))
>
> Q=THETA(4)
>
> V3=THETA(5)
>
> KA=THETA(6)*EXP(ETA(3))
>
>
>
> S2=V2/1000
> S3=V3/1000
>
> $THETA (0,20) (0,2) (0,100) (0,4) (0,40) (0,1.5)
> $OMEGA 0.1 0.1 0.1
>
>
> $ERROR
>
> IPRE=F
> Y=F*(1+EPS(1))+EPS(2)
>
> $SIGMA 0.4 0.5
> $ESTIMATION NOABORT SIG=2 MAX=9999 PRINT=10 METHOD=0
>
> $TABLE ID TIME IPRED DV CL KA D1 V2 V3 Q
For the two dataset combined, the following control file would work.
> $PROB MAD&POC
> $INPUT C ID AMT RATE TIME DV EVID TYPE
> $DATA full-data-set-1.csv IGNORE=C
> $SUBROUTINE ADVAN4 TRANS4
> $PK
>
> IF (TYPE.EQ.0) THEN
> CL=THETA(1)*EXP(ETA(1))
> ELSE
> CL=THETA(2)*EXP(ETA(2))
> ENDIF
>
> D1=THETA(3)
>
> V2=THETA(4)*EXP(ETA(3))
>
> Q=THETA(5)
>
> V3=THETA(6)
>
> IF (TYPE.EQ.0) THEN
> KA=THETA(7)*EXP(ETA(4))
> ELSE
> KA=THETA(8)*EXP(ETA(4))
> ENDIF
>
>
>
> S2=V2/1000
> S3=V3/1000
>
> $THETA (0,20) ;CL1 1
> $THETA (0, 25) ;CL2 2
> $THETA (0,5) ;D1 3
> $THETA (0,200) ;V2 4
> $THETA (0,20) ;Q 5
> $THETA (0,400) ;V3 6
> $THETA (0,1) ;KA1 7
> $THETA (0,0.75) ;KA2 8
>
> $OMEGA 0.3 0.3 0.7 0.7
>
>
> $ERROR
>
> IPRE=F
>
> Y1=F*EXP(EPS(1))+EPS(2)
> Y2=F*EXP(EPS(3))+EPS(4)
> Q1=0
> Q2=0
> IF (TYPE.EQ.0) Q1=1
> IF (TYPE.EQ.1) Q2=1
> Y=Q1*Y1+Q2*Y2
>
>
> $SIGMA 0.4 0.5 1 1
> $ESTIMATION NOABORT NSIG=2 MAX=9999 PRINT=10 METHOD=0
> $TABLE ID TIME IPRE DV
This is initial step for the fitting. Both would require optimization.
I run into difficulties when I move METHOD=0 to 1, can anyone give me
some suggestions on this? Thanks a lot.
Ping Ji, PHD
Clinical Discovery
Pharmaceutical Research Institute
Bristol-Myers Squibb
Princeton, NJ, 08543
_______________________________________________________