Dear All:
why do I get following error message with the following code? Any opinions regarding this problem and code will be highly appreciated. Thank you very much for your consideration and help. I am modeling parent drug and active metabolite at the same time with this code.
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
CL1=THETA(1)*EXP(ETA(1))
X
THE CHARACTERS IN ERROR ARE:
202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
4
$PROBLEM Parent drug & active metabolite ;Units: Time=mins, Concentration=uM(umols/L), Amount=umols
;Age=yrs, Weight=kg
$DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL
COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
$PK
CL1=THETA(1)*EXP(ETA(1))
V1=THETA(2)*EXP(ETA(2))
K20=THETA(3)*EXP(ETA(3))
Q=THETA(4)*EXP(ETA(4))
V2=THETA(5)*EXP(ETA(5))
V3=THETA(6)*EXP(ETA(6))
S1=V1
S2=V2
K12=CL1/V1
$DES
DADT(1)=-CL1*A(1)/V1
DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
DADT(3)=Q/V2*A(2)-Q/V3*A(3)
$ERROR
W=F
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
P1=0
IF(CMT.EQ.1)P1=1
P2=0
IF(CMT.EQ.2)P2=1
Y1=W*(1+ERR(1))
Y2=W*(1+ERR(2))
Y=P1*Y1+P2*Y2
$THETA
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
$OMEGA
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
$SIGMA
(0,,0.1000E+07)
(0,,0.1000E+07)
$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
$COVARIANCE
$TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
ID Time DV Anum AMT Gender WT Age CMT
1 0 0 1 3371.43 F 110.23 44 1
1 1 187 1 0 F 110.23 44 1
1 1 14.24 1 0 F 110.23 44 2
1 5 85.55 1 0 F 110.23 44 1
1 5 28.49 1 0 F 110.23 44 2
1 15 12.1 1 0 F 110.23 44 1
1 15 13.43 1 0 F 110.23 44 2
1 30 1.03 1 0 F 110.23 44 1
1 30 3.25 1 0 F 110.23 44 2
1 45 1.19 1 0 F 110.23 44 2
1 60 0.83 1 0 F 110.23 44 2
1 120 0.4 1 0 F 110.23 44 2
1 180 0.35 1 0 F 110.23 44 2
1 240 0.13 1 0 F 110.23 44 2
Zheng
parent drug and metabolite
Sam,
There is no problem with using CL1 in this code. This is not a reserved word and CL is not a required variable name when using $DES.
I fixed the more obvious errors so that the following code will run with NM-TRAN. However, the equations for parent and metabolites dont make sense to me. It looks like Amounts are being multiplied by volumes when they should be divided in order to get concentration.
Nick
$PROBLEM Parent drug & active metabolite
;Units: Time=mins,
;Concentration=uM(umols/L), Amount=umols
;Age=yrs, Weight=kg
$DATA AW2.CSV IGNORE=#
$INPUT ID TIME DV ANUM AMT SEX=DROP WT AGE CMT
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL
COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
$PK
CL1=THETA(1)*EXP(ETA(1))
V1=THETA(2)*EXP(ETA(2))
K20=THETA(3)*EXP(ETA(3))
Q=THETA(4)*EXP(ETA(4))
V2=THETA(5)*EXP(ETA(5))
V3=THETA(6)*EXP(ETA(6))
S1=V1
S2=V2
K12=CL1/V1
$DES
DADT(1)=-CL1*A(1)/V1
DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
DADT(3)=Q/V2*A(2)-Q/V3*A(3)
$ERROR
W=F
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
P1=0
IF(CMT.EQ.1)P1=1
P2=0
IF(CMT.EQ.2)P2=1
Y1=W*(1+ERR(1))
Y2=W*(1+ERR(2))
Y=P1*Y1+P2*Y2
$THETA
(0,,100000) ; CL1
(0,,100000) ; V1
(0,,100000) ; K20
(0,,100000) ; Q
(0,,100000) ; V2
(0,,100000) ; V3
$OMEGA
0.1 ; PPV_CL1
0.1 ; PPV_V1
0.1 ; PPV_K20
0.1 ; PPV_Q
0.1 ; PPV_V2
0.1 ; PPV_V3
$SIGMA
0.1 ; RUV_P
0.1 ; RUV_M
;$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
;$COVARIANCE
$TABLE ID TIME DV ANUM AMT WT AGE CMT EVID MDV IPRED IWRES
NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
Sam Liao wrote:
>
> Hi Zheng:
> It is complaining the CL1. You need to replace all CL1 with CL. But
> you have other problems in your $INPUT statement and the data set too.
>
> Sam
> > Dear All:
> >
> > why do I get following error message with the following code? Any
> > opinions regarding this problem and code will be highly appreciated.
> > Thank you very much for your consideration and help. I am modeling
> > parent drug and active metabolite at the same time with this code.
> >
> >
> > AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
> >
> > AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
> > CL1=THETA(1)*EXP(ETA(1))
> > X
> > THE CHARACTERS IN ERROR ARE:
> > 202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
> > 4
> >
> >
> > $PROBLEM Parent drug & active metabolite ;Units: Time=mins,
> > Concentration=uM(umols/L), Amount=umols
> > ;Age=yrs, Weight=kg
> > $DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
> >
> > $SUBROUTINE ADVAN6 TRANS1 TOL=5
> > $MODEL
> > COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
> >
> > $PK
> > CL1=THETA(1)*EXP(ETA(1))
> > V1=THETA(2)*EXP(ETA(2))
> > K20=THETA(3)*EXP(ETA(3))
> > Q=THETA(4)*EXP(ETA(4))
> > V2=THETA(5)*EXP(ETA(5))
> > V3=THETA(6)*EXP(ETA(6))
> > S1=V1
> > S2=V2
> > K12=CL1/V1
> >
> > $DES
> > DADT(1)=-CL1*A(1)/V1
> > DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> > DADT(3)=Q/V2*A(2)-Q/V3*A(3)
> >
> > $ERROR
> > W=F
> > IPRED=F
> > IRES=DV-IPRED
> > IWRES=IRES/W
> > P1=0
> > IF(CMT.EQ.1)P1=1
> > P2=0
> > IF(CMT.EQ.2)P2=1
> > Y1=W*(1+ERR(1))
> > Y2=W*(1+ERR(2))
> > Y=P1*Y1+P2*Y2
> >
> > $THETA
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > $OMEGA
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > $SIGMA
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> >
> > $ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
> >
> > $COVARIANCE
> >
> > $TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
> > NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
> >
> > ID Time DV Anum AMT Gender WT Age CMT
> > 1 0 0 1 3371.43 F 110.23 44 1
> > 1 1 187 1 0 F 110.23 44 1
> > 1 1 14.24 1 0 F 110.23 44 2
> > 1 5 85.55 1 0 F 110.23 44 1
> > 1 5 28.49 1 0 F 110.23 44 2
> > 1 15 12.1 1 0 F 110.23 44 1
> > 1 15 13.43 1 0 F 110.23 44 2
> > 1 30 1.03 1 0 F 110.23 44 1
> > 1 30 3.25 1 0 F 110.23 44 2
> > 1 45 1.19 1 0 F 110.23 44 2
> > 1 60 0.83 1 0 F 110.23 44 2
> > 1 120 0.4 1 0 F 110.23 44 2
> > 1 180 0.35 1 0 F 110.23 44 2
> > 1 240 0.13 1 0 F 110.23 44 2
> >
> > Zheng
> >
> >
> >
> >
> >
--
Nick Holford, Dept Pharmacology & Clinical Pharmacology
University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand
email:n.holford
http://www.health.auckland.ac.nz/pharmacology/staff/nholford/
I cannot see any problems with the equations, looks correct to me:
K=CL1/V1
DADT1=-K*A(1)
S1=V1
etc.
Leonid
Nick Holford wrote:
> Sam,
>
> There is no problem with using CL1 in this code. This is not a reserved word and CL is not a required variable name when using $DES.
>
> I fixed the more obvious errors so that the following code will run with NM-TRAN. However, the equations for parent and metabolites dont make sense to me. It looks like Amounts are being multiplied by volumes when they should be divided in order to get concentration.
>
> Nick
>
> $PROBLEM Parent drug & active metabolite
> ;Units: Time=mins,
> ;Concentration=uM(umols/L), Amount=umols
> ;Age=yrs, Weight=kg
> $DATA AW2.CSV IGNORE=#
> $INPUT ID TIME DV ANUM AMT SEX=DROP WT AGE CMT
>
> $SUBROUTINE ADVAN6 TRANS1 TOL=5
> $MODEL
> COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
>
> $PK
> CL1=THETA(1)*EXP(ETA(1))
> V1=THETA(2)*EXP(ETA(2))
> K20=THETA(3)*EXP(ETA(3))
> Q=THETA(4)*EXP(ETA(4))
> V2=THETA(5)*EXP(ETA(5))
> V3=THETA(6)*EXP(ETA(6))
> S1=V1
> S2=V2
> K12=CL1/V1
>
> $DES
> DADT(1)=-CL1*A(1)/V1
> DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> DADT(3)=Q/V2*A(2)-Q/V3*A(3)
>
> $ERROR
> W=F
> IPRED=F
> IRES=DV-IPRED
> IWRES=IRES/W
> P1=0
> IF(CMT.EQ.1)P1=1
> P2=0
> IF(CMT.EQ.2)P2=1
> Y1=W*(1+ERR(1))
> Y2=W*(1+ERR(2))
> Y=P1*Y1+P2*Y2
>
> $THETA
> (0,,100000) ; CL1
> (0,,100000) ; V1
> (0,,100000) ; K20
> (0,,100000) ; Q
> (0,,100000) ; V2
> (0,,100000) ; V3
> $OMEGA
> 0.1 ; PPV_CL1
> 0.1 ; PPV_V1
> 0.1 ; PPV_K20
> 0.1 ; PPV_Q
> 0.1 ; PPV_V2
> 0.1 ; PPV_V3
> $SIGMA
> 0.1 ; RUV_P
> 0.1 ; RUV_M
>
> ;$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
> ;$COVARIANCE
> $TABLE ID TIME DV ANUM AMT WT AGE CMT EVID MDV IPRED IWRES
> NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
>
>
> Sam Liao wrote:
>> Hi Zheng:
>> It is complaining the CL1. You need to replace all CL1 with CL. But
>> you have other problems in your $INPUT statement and the data set too.
>>
>> Sam
>>> Dear All:
>>>
>>> why do I get following error message with the following code? Any
>>> opinions regarding this problem and code will be highly appreciated.
>>> Thank you very much for your consideration and help. I am modeling
>>> parent drug and active metabolite at the same time with this code.
>>>
>>>
>>> AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
>>>
>>> AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
>>> CL1=THETA(1)*EXP(ETA(1))
>>> X
>>> THE CHARACTERS IN ERROR ARE:
>>> 202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
>>> 4
>>>
>>>
>>> $PROBLEM Parent drug & active metabolite ;Units: Time=mins,
>>> Concentration=uM(umols/L), Amount=umols
>>> ;Age=yrs, Weight=kg
>>> $DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
>>>
>>> $SUBROUTINE ADVAN6 TRANS1 TOL=5
>>> $MODEL
>>> COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
>>>
>>> $PK
>>> CL1=THETA(1)*EXP(ETA(1))
>>> V1=THETA(2)*EXP(ETA(2))
>>> K20=THETA(3)*EXP(ETA(3))
>>> Q=THETA(4)*EXP(ETA(4))
>>> V2=THETA(5)*EXP(ETA(5))
>>> V3=THETA(6)*EXP(ETA(6))
>>> S1=V1
>>> S2=V2
>>> K12=CL1/V1
>>>
>>> $DES
>>> DADT(1)=-CL1*A(1)/V1
>>> DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
>>> DADT(3)=Q/V2*A(2)-Q/V3*A(3)
>>>
>>> $ERROR
>>> W=F
>>> IPRED=F
>>> IRES=DV-IPRED
>>> IWRES=IRES/W
>>> P1=0
>>> IF(CMT.EQ.1)P1=1
>>> P2=0
>>> IF(CMT.EQ.2)P2=1
>>> Y1=W*(1+ERR(1))
>>> Y2=W*(1+ERR(2))
>>> Y=P1*Y1+P2*Y2
>>>
>>> $THETA
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> $OMEGA
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>> $SIGMA
>>> (0,,0.1000E+07)
>>> (0,,0.1000E+07)
>>>
>>> $ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
>>>
>>> $COVARIANCE
>>>
>>> $TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
>>> NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
>>>
>>> ID Time DV Anum AMT Gender WT Age CMT
>>> 1 0 0 1 3371.43 F 110.23 44 1
>>> 1 1 187 1 0 F 110.23 44 1
>>> 1 1 14.24 1 0 F 110.23 44 2
>>> 1 5 85.55 1 0 F 110.23 44 1
>>> 1 5 28.49 1 0 F 110.23 44 2
>>> 1 15 12.1 1 0 F 110.23 44 1
>>> 1 15 13.43 1 0 F 110.23 44 2
>>> 1 30 1.03 1 0 F 110.23 44 1
>>> 1 30 3.25 1 0 F 110.23 44 2
>>> 1 45 1.19 1 0 F 110.23 44 2
>>> 1 60 0.83 1 0 F 110.23 44 2
>>> 1 120 0.4 1 0 F 110.23 44 2
>>> 1 180 0.35 1 0 F 110.23 44 2
>>> 1 240 0.13 1 0 F 110.23 44 2
>>>
>>> Zheng
>>>
>>>
>>>
>>>
>>>
>
> --
> Nick Holford, Dept Pharmacology & Clinical Pharmacology
> University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand
> email:n.holford
> http://www.health.auckland.ac.nz/pharmacology/staff/nholford/
>
>
Hi Zheng:
It is complaining the CL1. You need to replace all CL1 with CL. But you have other problems in your $INPUT statement and the data set too.
Sam
> Dear All:
>
> why do I get following error message with the following code? Any opinions regarding this problem and code will be highly appreciated. Thank you very much for your consideration and help. I am modeling parent drug and active metabolite at the same time with this code.
>
> AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
>
> AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
> CL1=THETA(1)*EXP(ETA(1))
> X
> THE CHARACTERS IN ERROR ARE:
> 202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
> 4
>
> $PROBLEM Parent drug & active metabolite ;Units: Time=mins, Concentration=uM(umols/L), Amount=umols
>
> ;Age=yrs, Weight=kg
> $DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
>
> $SUBROUTINE ADVAN6 TRANS1 TOL=5
> $MODEL
> COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
>
> $PK
> CL1=THETA(1)*EXP(ETA(1))
> V1=THETA(2)*EXP(ETA(2))
> K20=THETA(3)*EXP(ETA(3))
> Q=THETA(4)*EXP(ETA(4))
> V2=THETA(5)*EXP(ETA(5))
> V3=THETA(6)*EXP(ETA(6))
> S1=V1
> S2=V2
> K12=CL1/V1
>
> $DES
> DADT(1)=-CL1*A(1)/V1
> DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> DADT(3)=Q/V2*A(2)-Q/V3*A(3)
>
> $ERROR
> W=F
> IPRED=F
> IRES=DV-IPRED
> IWRES=IRES/W
> P1=0
> IF(CMT.EQ.1)P1=1
> P2=0
> IF(CMT.EQ.2)P2=1
> Y1=W*(1+ERR(1))
> Y2=W*(1+ERR(2))
> Y=P1*Y1+P2*Y2
>
> $THETA
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> $OMEGA
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> $SIGMA
> (0,,0.1000E+07)
> (0,,0.1000E+07)
>
> $ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
>
> $COVARIANCE
>
> $TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
> NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
>
> ID Time DV Anum AMT Gender WT Age CMT
> 1 0 0 1 3371.43 F 110.23 44 1
> 1 1 187 1 0 F 110.23 44 1
> 1 1 14.24 1 0 F 110.23 44 2
> 1 5 85.55 1 0 F 110.23 44 1
> 1 5 28.49 1 0 F 110.23 44 2
> 1 15 12.1 1 0 F 110.23 44 1
> 1 15 13.43 1 0 F 110.23 44 2
> 1 30 1.03 1 0 F 110.23 44 1
> 1 30 3.25 1 0 F 110.23 44 2
> 1 45 1.19 1 0 F 110.23 44 2
> 1 60 0.83 1 0 F 110.23 44 2
> 1 120 0.4 1 0 F 110.23 44 2
> 1 180 0.35 1 0 F 110.23 44 2
> 1 240 0.13 1 0 F 110.23 44 2
>
> Zheng
begin:vcard
fn:Sam Liao
n:Liao;Sam
org:Pharmax Research Inc.
adr;dom:;;199 Pierce Street, Suite 817;Somerset;NJ;08873
email;internet:[EMAIL PROTECTED]
tel;work:2019882043
tel;fax:720-2946783
version:2.1
end:vcard
I agree with Sam,
Getting rid of TRANS1 would help - as I undertand it, it isnt needed
with ADVAN6 and that's why the model is looking for CL.
Regards,
David
David Foster, PhD
Lecturer
School of Pharmacy and Medical Sciences
Room P4-08
City East Campus
University of South Australia
Adelaide SA 5000
CRICOS Provider Number: 00121B
Phone: 61 8 8302 2055
Fax: 61 8 8302 2389
Email: [EMAIL PROTECTED]
Note: This message (including any attachments) is confidential and is
intended only for the individual(s) named. If you have received it by
mistake, please let the sender know by e-mail reply and delete it from
your system.
Quoted reply history
> -----Original Message-----
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Zheng Lu
> Sent: Monday, 5 March 2007 4:18 AM
> To: [email protected]
> Subject: [NMusers] parent drug and metabolite
>
> Dear All:
>
> why do I get following error message with the following code?
> Any opinions regarding this problem and code will be highly
> appreciated.
> Thank you very much for your consideration and help. I am
> modeling parent drug and active metabolite at the same time
> with this code.
>
>
> AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
>
> AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
> CL1=THETA(1)*EXP(ETA(1))
> X
> THE CHARACTERS IN ERROR ARE:
> 202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
> 4
>
>
> $PROBLEM Parent drug & active metabolite ;Units: Time=mins,
> Concentration=uM(umols/L), Amount=umols
> ;Age=yrs,
> Weight=kg $DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT
> SEX WT AGE CMT EVID MDV
>
> $SUBROUTINE ADVAN6 TRANS1 TOL=5
> $MODEL
> COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
>
> $PK
> CL1=THETA(1)*EXP(ETA(1))
> V1=THETA(2)*EXP(ETA(2))
> K20=THETA(3)*EXP(ETA(3))
> Q=THETA(4)*EXP(ETA(4))
> V2=THETA(5)*EXP(ETA(5))
> V3=THETA(6)*EXP(ETA(6))
> S1=V1
> S2=V2
> K12=CL1/V1
>
> $DES
> DADT(1)=-CL1*A(1)/V1
> DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> DADT(3)=Q/V2*A(2)-Q/V3*A(3)
>
> $ERROR
> W=F
> IPRED=F
> IRES=DV-IPRED
> IWRES=IRES/W
> P1=0
> IF(CMT.EQ.1)P1=1
> P2=0
> IF(CMT.EQ.2)P2=1
> Y1=W*(1+ERR(1))
> Y2=W*(1+ERR(2))
> Y=P1*Y1+P2*Y2
>
> $THETA
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> $OMEGA
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> (0,,0.1000E+07)
> $SIGMA
> (0,,0.1000E+07)
> (0,,0.1000E+07)
>
> $ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
>
> $COVARIANCE
>
> $TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED
> IWRES NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
>
> ID Time DV Anum AMT Gender WT Age CMT
> 1 0 0 1 3371.43 F 110.23 44 1
> 1 1 187 1 0 F 110.23 44 1
> 1 1 14.24 1 0 F 110.23 44 2
> 1 5 85.55 1 0 F 110.23 44 1
> 1 5 28.49 1 0 F 110.23 44 2
> 1 15 12.1 1 0 F 110.23 44 1
> 1 15 13.43 1 0 F 110.23 44 2
> 1 30 1.03 1 0 F 110.23 44 1
> 1 30 3.25 1 0 F 110.23 44 2
> 1 45 1.19 1 0 F 110.23 44 2
> 1 60 0.83 1 0 F 110.23 44 2
> 1 120 0.4 1 0 F 110.23 44 2
> 1 180 0.35 1 0 F 110.23 44 2
> 1 240 0.13 1 0 F 110.23 44 2
>
> Zheng
>
>
>
Sam,
There is no problem with using CL1 in this code. This is not a reserved word
and CL is not a required variable name when using $DES.
I fixed the more obvious errors so that the following code will run with
NM-TRAN. However, the equations for parent and metabolites dont make sense to
me. It looks like Amounts are being multiplied by volumes when they should be
divided in order to get concentration.
Nick
$PROBLEM Parent drug & active metabolite
;Units: Time=mins,
;Concentration=uM(umols/L), Amount=umols
;Age=yrs, Weight=kg
$DATA AW2.CSV IGNORE=#
$INPUT ID TIME DV ANUM AMT SEX=DROP WT AGE CMT
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL
COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
$PK
CL1=THETA(1)*EXP(ETA(1))
V1=THETA(2)*EXP(ETA(2))
K20=THETA(3)*EXP(ETA(3))
Q=THETA(4)*EXP(ETA(4))
V2=THETA(5)*EXP(ETA(5))
V3=THETA(6)*EXP(ETA(6))
S1=V1
S2=V2
K12=CL1/V1
$DES
DADT(1)=-CL1*A(1)/V1
DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
DADT(3)=Q/V2*A(2)-Q/V3*A(3)
$ERROR
W=F
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
P1=0
IF(CMT.EQ.1)P1=1
P2=0
IF(CMT.EQ.2)P2=1
Y1=W*(1+ERR(1))
Y2=W*(1+ERR(2))
Y=P1*Y1+P2*Y2
$THETA
(0,,100000) ; CL1
(0,,100000) ; V1
(0,,100000) ; K20
(0,,100000) ; Q
(0,,100000) ; V2
(0,,100000) ; V3
$OMEGA
0.1 ; PPV_CL1
0.1 ; PPV_V1
0.1 ; PPV_K20
0.1 ; PPV_Q
0.1 ; PPV_V2
0.1 ; PPV_V3
$SIGMA
0.1 ; RUV_P
0.1 ; RUV_M
;$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
;$COVARIANCE
$TABLE ID TIME DV ANUM AMT WT AGE CMT EVID MDV IPRED IWRES
NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
Sam Liao wrote:
>
> Hi Zheng:
> It is complaining the CL1. You need to replace all CL1 with CL. But
> you have other problems in your $INPUT statement and the data set too.
>
> Sam
> > Dear All:
> >
> > why do I get following error message with the following code? Any
> > opinions regarding this problem and code will be highly appreciated.
> > Thank you very much for your consideration and help. I am modeling
> > parent drug and active metabolite at the same time with this code.
> >
> >
> > AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
> >
> > AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
> > CL1=THETA(1)*EXP(ETA(1))
> > X
> > THE CHARACTERS IN ERROR ARE:
> > 202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
> > 4
> >
> >
> > $PROBLEM Parent drug & active metabolite ;Units: Time=mins,
> > Concentration=uM(umols/L), Amount=umols
> > ;Age=yrs, Weight=kg
> > $DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
> >
> > $SUBROUTINE ADVAN6 TRANS1 TOL=5
> > $MODEL
> > COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
> >
> > $PK
> > CL1=THETA(1)*EXP(ETA(1))
> > V1=THETA(2)*EXP(ETA(2))
> > K20=THETA(3)*EXP(ETA(3))
> > Q=THETA(4)*EXP(ETA(4))
> > V2=THETA(5)*EXP(ETA(5))
> > V3=THETA(6)*EXP(ETA(6))
> > S1=V1
> > S2=V2
> > K12=CL1/V1
> >
> > $DES
> > DADT(1)=-CL1*A(1)/V1
> > DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> > DADT(3)=Q/V2*A(2)-Q/V3*A(3)
> >
> > $ERROR
> > W=F
> > IPRED=F
> > IRES=DV-IPRED
> > IWRES=IRES/W
> > P1=0
> > IF(CMT.EQ.1)P1=1
> > P2=0
> > IF(CMT.EQ.2)P2=1
> > Y1=W*(1+ERR(1))
> > Y2=W*(1+ERR(2))
> > Y=P1*Y1+P2*Y2
> >
> > $THETA
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > $OMEGA
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> > $SIGMA
> > (0,,0.1000E+07)
> > (0,,0.1000E+07)
> >
> > $ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
> >
> > $COVARIANCE
> >
> > $TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
> > NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
> >
> > ID Time DV Anum AMT Gender WT Age CMT
> > 1 0 0 1 3371.43 F 110.23 44 1
> > 1 1 187 1 0 F 110.23 44 1
> > 1 1 14.24 1 0 F 110.23 44 2
> > 1 5 85.55 1 0 F 110.23 44 1
> > 1 5 28.49 1 0 F 110.23 44 2
> > 1 15 12.1 1 0 F 110.23 44 1
> > 1 15 13.43 1 0 F 110.23 44 2
> > 1 30 1.03 1 0 F 110.23 44 1
> > 1 30 3.25 1 0 F 110.23 44 2
> > 1 45 1.19 1 0 F 110.23 44 2
> > 1 60 0.83 1 0 F 110.23 44 2
> > 1 120 0.4 1 0 F 110.23 44 2
> > 1 180 0.35 1 0 F 110.23 44 2
> > 1 240 0.13 1 0 F 110.23 44 2
> >
> > Zheng
> >
> >
> >
> >
> >
--
Nick Holford, Dept Pharmacology & Clinical Pharmacology
University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand
email:[EMAIL PROTECTED] tel:+64(9)373-7599x86730 fax:373-7556
http://www.health.auckland.ac.nz/pharmacology/staff/nholford/
I cannot see any problems with the equations, looks correct to me:
K=CL1/V1
DADT1=-K*A(1)
S1=V1
etc.
Leonid
Nick Holford wrote:
> Sam,
>
> There is no problem with using CL1 in this code. This is not a reserved word
> and CL is not a required variable name when using $DES.
>
> I fixed the more obvious errors so that the following code will run with
> NM-TRAN. However, the equations for parent and metabolites dont make sense to
> me. It looks like Amounts are being multiplied by volumes when they should be
> divided in order to get concentration.
>
> Nick
>
> $PROBLEM Parent drug & active metabolite ;Units: Time=mins, ;Concentration=uM(umols/L), Amount=umols
>
> ;Age=yrs, Weight=kg
>
> $DATA AW2.CSV IGNORE=# $INPUT ID TIME DV ANUM AMT SEX=DROP WT AGE CMT
>
> $SUBROUTINE ADVAN6 TRANS1 TOL=5
> $MODEL
> COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
>
> $PK
> CL1=THETA(1)*EXP(ETA(1))
> V1=THETA(2)*EXP(ETA(2))
> K20=THETA(3)*EXP(ETA(3))
> Q=THETA(4)*EXP(ETA(4))
> V2=THETA(5)*EXP(ETA(5))
> V3=THETA(6)*EXP(ETA(6))
> S1=V1
> S2=V2
> K12=CL1/V1
>
> $DES
> DADT(1)=-CL1*A(1)/V1
> DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> DADT(3)=Q/V2*A(2)-Q/V3*A(3)
>
> $ERROR
> W=F
> IPRED=F
> IRES=DV-IPRED
> IWRES=IRES/W
> P1=0
> IF(CMT.EQ.1)P1=1
> P2=0
> IF(CMT.EQ.2)P2=1
> Y1=W*(1+ERR(1))
> Y2=W*(1+ERR(2))
> Y=P1*Y1+P2*Y2
>
> $THETA
> (0,,100000) ; CL1
> (0,,100000) ; V1
> (0,,100000) ; K20
> (0,,100000) ; Q
> (0,,100000) ; V2
> (0,,100000) ; V3
> $OMEGA
> 0.1 ; PPV_CL1
> 0.1 ; PPV_V1
> 0.1 ; PPV_K20
> 0.1 ; PPV_Q
> 0.1 ; PPV_V2
> 0.1 ; PPV_V3
> $SIGMA
> 0.1 ; RUV_P
> 0.1 ; RUV_M
>
> ;$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
> ;$COVARIANCE
> $TABLE ID TIME DV ANUM AMT WT AGE CMT EVID MDV IPRED IWRES
> NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
>
> Sam Liao wrote:
>
> > Hi Zheng:
> > It is complaining the CL1. You need to replace all CL1 with CL. But
> > you have other problems in your $INPUT statement and the data set too.
> >
> > Sam
> >
> > > Dear All:
> > >
> > > why do I get following error message with the following code? Any
> > > opinions regarding this problem and code will be highly appreciated.
> > > Thank you very much for your consideration and help. I am modeling
> > > parent drug and active metabolite at the same time with this code.
> > >
> > > AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
> > >
> > > AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
> > > CL1=THETA(1)*EXP(ETA(1))
> > > X
> > > THE CHARACTERS IN ERROR ARE:
> > > 202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
> > > 4
> > >
> > > $PROBLEM Parent drug & active metabolite ;Units: Time=mins,
> > > Concentration=uM(umols/L), Amount=umols
> > > ;Age=yrs, Weight=kg
> > > $DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
> > >
> > > $SUBROUTINE ADVAN6 TRANS1 TOL=5
> > > $MODEL
> > > COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
> > >
> > > $PK
> > > CL1=THETA(1)*EXP(ETA(1))
> > > V1=THETA(2)*EXP(ETA(2))
> > > K20=THETA(3)*EXP(ETA(3))
> > > Q=THETA(4)*EXP(ETA(4))
> > > V2=THETA(5)*EXP(ETA(5))
> > > V3=THETA(6)*EXP(ETA(6))
> > > S1=V1
> > > S2=V2
> > > K12=CL1/V1
> > >
> > > $DES
> > > DADT(1)=-CL1*A(1)/V1
> > > DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
> > > DADT(3)=Q/V2*A(2)-Q/V3*A(3)
> > >
> > > $ERROR
> > > W=F
> > > IPRED=F
> > > IRES=DV-IPRED
> > > IWRES=IRES/W
> > > P1=0
> > > IF(CMT.EQ.1)P1=1
> > > P2=0
> > > IF(CMT.EQ.2)P2=1
> > > Y1=W*(1+ERR(1))
> > > Y2=W*(1+ERR(2))
> > > Y=P1*Y1+P2*Y2
> > >
> > > $THETA
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > $OMEGA
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > > $SIGMA
> > > (0,,0.1000E+07)
> > > (0,,0.1000E+07)
> > >
> > > $ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
> > >
> > > $COVARIANCE
> > >
> > > $TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
> > > NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
> > >
> > > ID Time DV Anum AMT Gender WT Age CMT
> > > 1 0 0 1 3371.43 F 110.23 44 1
> > > 1 1 187 1 0 F 110.23 44 1
> > > 1 1 14.24 1 0 F 110.23 44 2
> > > 1 5 85.55 1 0 F 110.23 44 1
> > > 1 5 28.49 1 0 F 110.23 44 2
> > > 1 15 12.1 1 0 F 110.23 44 1
> > > 1 15 13.43 1 0 F 110.23 44 2
> > > 1 30 1.03 1 0 F 110.23 44 1
> > > 1 30 3.25 1 0 F 110.23 44 2
> > > 1 45 1.19 1 0 F 110.23 44 2
> > > 1 60 0.83 1 0 F 110.23 44 2
> > > 1 120 0.4 1 0 F 110.23 44 2
> > > 1 180 0.35 1 0 F 110.23 44 2
> > > 1 240 0.13 1 0 F 110.23 44 2
> > >
> > > Zheng
>
> --
> Nick Holford, Dept Pharmacology & Clinical Pharmacology
> University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand
> email:[EMAIL PROTECTED] tel:+64(9)373-7599x86730 fax:373-7556
> http://www.health.auckland.ac.nz/pharmacology/staff/nholford/
Zheng,
When you edit your control stream, make sure there is no format issue that
NONMEM can not recognize.
Alan
Quoted reply history
-----Original Message-----
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Zheng Lu
Sent: Sunday, March 04, 2007 12:48 PM
To: [email protected]
Subject: [NMusers] parent drug and metabolite
Dear All:
why do I get following error message with the following code? Any
opinions regarding this problem and code will be highly appreciated.
Thank you very much for your consideration and help. I am modeling
parent drug and active metabolite at the same time with this code.
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
CL1=THETA(1)*EXP(ETA(1))
X
THE CHARACTERS IN ERROR ARE:
202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
4
$PROBLEM Parent drug & active metabolite ;Units: Time=mins,
Concentration=uM(umols/L), Amount=umols
;Age=yrs, Weight=kg
$DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL
COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
$PK
CL1=THETA(1)*EXP(ETA(1))
V1=THETA(2)*EXP(ETA(2))
K20=THETA(3)*EXP(ETA(3))
Q=THETA(4)*EXP(ETA(4))
V2=THETA(5)*EXP(ETA(5))
V3=THETA(6)*EXP(ETA(6))
S1=V1
S2=V2
K12=CL1/V1
$DES
DADT(1)=-CL1*A(1)/V1
DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
DADT(3)=Q/V2*A(2)-Q/V3*A(3)
$ERROR
W=F
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
P1=0
IF(CMT.EQ.1)P1=1
P2=0
IF(CMT.EQ.2)P2=1
Y1=W*(1+ERR(1))
Y2=W*(1+ERR(2))
Y=P1*Y1+P2*Y2
$THETA
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
$OMEGA
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
$SIGMA
(0,,0.1000E+07)
(0,,0.1000E+07)
$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
$COVARIANCE
$TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
ID Time DV Anum AMT Gender WT Age CMT
1 0 0 1 3371.43 F 110.23 44 1
1 1 187 1 0 F 110.23 44 1
1 1 14.24 1 0 F 110.23 44 2
1 5 85.55 1 0 F 110.23 44 1
1 5 28.49 1 0 F 110.23 44 2
1 15 12.1 1 0 F 110.23 44 1
1 15 13.43 1 0 F 110.23 44 2
1 30 1.03 1 0 F 110.23 44 1
1 30 3.25 1 0 F 110.23 44 2
1 45 1.19 1 0 F 110.23 44 2
1 60 0.83 1 0 F 110.23 44 2
1 120 0.4 1 0 F 110.23 44 2
1 180 0.35 1 0 F 110.23 44 2
1 240 0.13 1 0 F 110.23 44 2
Zheng
Dear Zheng,
Could it be strange characters? Fortan is very picky about characters
in its control steam, a property that trips me up every time I use the
<TAB> key/character in my control stream. Please go through you code
very carefully, deleting any whitespace or manually replacing it with
spaces. Many of us use programming-specific text editors for this
reason over something like WordPad or Notepad. There are tools to allow
you to view the ASCII codes of every character in a file, which
sometimes reveals odd linefeed/carriage return issues.
----
Mike Dodds, PhD
PK/PD, Pre-clinical Development
ZymoGenetics, Inc.
1201 Eastlake Avenue
Seattle, WA, 98102, USA
Tel: 206-442-6737
Email: [EMAIL PROTECTED]
Quoted reply history
-----Original Message-----
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Xiao, Alan
Sent: Sunday, March 04, 2007 9:07 PM
To: Zheng Lu; [email protected]
Subject: RE: [NMusers] parent drug and metabolite
Zheng,
When you edit your control stream, make sure there is no format issue
that NONMEM can not recognize.
Alan
-----Original Message-----
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Zheng Lu
Sent: Sunday, March 04, 2007 12:48 PM
To: [email protected]
Subject: [NMusers] parent drug and metabolite
Dear All:
why do I get following error message with the following code? Any
opinions regarding this problem and code will be highly appreciated.
Thank you very much for your consideration and help. I am modeling
parent drug and active metabolite at the same time with this code.
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
AN ERROR WAS FOUND ON LINE 13 AT THE APPROXIMATE POSITION NOTED:
CL1=THETA(1)*EXP(ETA(1))
X
THE CHARACTERS IN ERROR ARE:
202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
4
$PROBLEM Parent drug & active metabolite ;Units: Time=mins,
Concentration=uM(umols/L), Amount=umols
;Age=yrs, Weight=kg
$DATA AW2.CSV IGNORE=C $INPUT ID TIME DV ANUM AMT SEX WT AGE CMT EVID
MDV
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL
COMP=(PARENT) COMP=(METAB1) COMP=(METAB2)
$PK
CL1=THETA(1)*EXP(ETA(1))
V1=THETA(2)*EXP(ETA(2))
K20=THETA(3)*EXP(ETA(3))
Q=THETA(4)*EXP(ETA(4))
V2=THETA(5)*EXP(ETA(5))
V3=THETA(6)*EXP(ETA(6))
S1=V1
S2=V2
K12=CL1/V1
$DES
DADT(1)=-CL1*A(1)/V1
DADT(2)=CL1*A(1)/V1+Q/V3*A(3)-Q/V2*A(2)-K20*A(2)
DADT(3)=Q/V2*A(2)-Q/V3*A(3)
$ERROR
W=F
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
P1=0
IF(CMT.EQ.1)P1=1
P2=0
IF(CMT.EQ.2)P2=1
Y1=W*(1+ERR(1))
Y2=W*(1+ERR(2))
Y=P1*Y1+P2*Y2
$THETA
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
$OMEGA
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
(0,,0.1000E+07)
$SIGMA
(0,,0.1000E+07)
(0,,0.1000E+07)
$ESTIMATION MAXEVAL=9999 PRINT=5 POSTHOC
$COVARIANCE
$TABLE ID TIME DV ANUM AMT SEX WT AGE CMT EVID MDV IPRED IWRES
NOPRINT ONEHEADER FILE=OUTAMIWR.FIT
ID Time DV Anum AMT Gender WT Age CMT
1 0 0 1 3371.43 F 110.23 44 1
1 1 187 1 0 F 110.23 44 1
1 1 14.24 1 0 F 110.23 44 2
1 5 85.55 1 0 F 110.23 44 1
1 5 28.49 1 0 F 110.23 44 2
1 15 12.1 1 0 F 110.23 44 1
1 15 13.43 1 0 F 110.23 44 2
1 30 1.03 1 0 F 110.23 44 1
1 30 3.25 1 0 F 110.23 44 2
1 45 1.19 1 0 F 110.23 44 2
1 60 0.83 1 0 F 110.23 44 2
1 120 0.4 1 0 F 110.23 44 2
1 180 0.35 1 0 F 110.23 44 2
1 240 0.13 1 0 F 110.23 44 2
Zheng