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