RE: parent drug and metabolite
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
>
>
>