NONMEM questions
From: Ping Ji ping.ji@bms.com
Subject: [NMusers] NONMEM questions
Date: Fri, 18 Mar 2005 11:53:14 -0500
Dear users,
I am new to this community. I am trying to fit multiple dose data from
healthy subjects (TYPE=0) and patients (TYPE=1). I use
two-compartment with zero- and first-order absorption and absorption
lag time in patients. Shown below is the control file. My problem is
that output is always "minimization terminated". Can anyone help me
out to see what is the problem?
> $PROB MAD&POC
> $INPUT C ID AMT RATE TIME DV EVID TYPE
> $DATA full-data-set.csv IGNORE=C
> $SUBROUTINE ADVAN4 TRANS4
> $PK
>
> IF (TYPE.EQ.0) TVCL=THETA(1)
> IF (TYPE.EQ.1) TVCL=THETA(2)
> CL=TVCL*EXP(ETA(1))
>
> D1=THETA(3)
>
> IF(TYPE.EQ.0) TVV2=THETA(4)
> IF (TYPE.EQ.1) TVV2=THETA(5)
> V2=TVV2*EXP(ETA(2))
>
> Q=THETA(6)
>
> V3=THETA(7)
>
> KA=THETA(8)*EXP(ETA(3))
>
> IF(TYPE.EQ.1) ALAG1=THETA(9)
> S2=V2/1000
> S3=V3/1000
>
> $THETA (0,10) (0, 10) (0,1) (0,20) (0,20) (0,4) (0,160) (0,0.3) (0,0.1)
> $OMEGA BLOCK(2)
> 0.1 0.05 0.2
> $OMEGA 0.5
> $ERROR
> IPRE=F
> Y1=F*(1+EPS(1))+EPS(2)
> Y2=F*(1+EPS(1))+EPS(3)
> Q1=0
> Q2=0
> IF (TYPE.EQ.0) Q1=1
> IF (TYPE.EQ.1) Q2=1
> Y=Q1*Y1+Q2*Y2
>
>
> $SIGMA 0.4 10 10
> $ESTIMATION NOABORT SIG=6 MAX=9999 PRINT=10
> $TABLE ID TIME IPRE DV PRED EVID TYPE