202 Fortran Syntax of wings for NM7
Dear All,
I have exact same NONMEM code run in two command windows, NM7 and
wings for NM7. The NONMEM run is successful in the NM7 but cannot
start in wings for NM7. Below is the error message.
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
AN ERROR WAS FOUND ON LINE 25 AT THE APPROXIMATE POSITION NOTED:
MU_1 =THETA(12) LOG(THETA(1))
X
THE CHARACTERS IN ERROR ARE: LOG
202 FORTRAN SYNTAX IS INCORRECT OR INAPPROPRIATE IN THIS CONTEXT.
Error in NM-TRAN_Execution with ......
After several tests, I realized that it may be due to the setting of
wings for NM7, such as parameter size.
In my model, I have 12 THETAS, 10 ETAS and 3 SIGMA. If I delete
THETA(12) , everything is working.
We have NONMEM 7.1.0 and Wings for Nonmem version 703 installed
with their default setting.
Can anyone help me to handle this problem? Is this due to the size of
parameter? How to increase the parameter size of wings for NM7.
Here is my NONMEM code:
$PROB
$DATA
$INPUT
$SUBROUTINE ADVAN6 TOL=9
$MODEL NCOMP=5
COMP=(GI,DEFDOSE) ;COMPARTMENT 1 DEPOSIT IN GI TRACT
COMP=(CENTER) ;COMPARTMENT 2 FOR CENTER COMPARTMENT
COMP=(PERI) ;COMPARTMENT 3 FOR PERIPHERAL COMPARTMENT
COMP=(MET) ;COMPARTMENT 4 FOR METABOLITE CENTER COMPARTMENT
COMP=(FREE) ;COMPARTMENT 5 FOR FREE UNBOUND COMPARTMENT
$PK
IF (CNIT.EQ.1.OR.CNIT.EQ.3.OR.CNIT.EQ.-1) CSP=1
IF (CNIT.EQ.2.OR.CNIT.EQ.4) CSP=0
IF (ALB.NE.-1) NALB=ALB/3.4
IF (ALB.EQ.-1) NALB=1
MU_1 = LOG(THETA(1))
ALAG1 = EXP(MU_1+ETA(1))
MU_2 = LOG(THETA(2))
K12 = EXP(MU_2+ETA(2))
MU_3 = LOG(THETA(3))+0.75*LOG(BW/70)
Q = EXP(MU_3+ETA(3))
MU_4 = LOG(THETA(4))-0.25*LOG(BW/70)+THETA(12)*LOG(NALB)
K24 = EXP(MU_4+ETA(4))
MU_5 = LOG(THETA(5))-0.25*LOG(BW/70)
K40 = EXP(MU_5+ETA(5))
MU_6 = LOG(THETA(6))+CSP*LOG(THETA(11))
K41 = (EXP(MU_6+ETA(6)))
MU_7 = LOG(THETA(7)) + LOG(BW/70)
V2 = EXP(MU_7+ETA(7))
MU_8 = LOG(THETA(8)) + LOG(BW/70)
V3 = EXP(MU_8+ETA(8))
MU_9 = LOG(THETA(9)) + LOG(BW/70)
V4 = EXP(MU_9+ETA(9))
MU_10 = LOG(THETA(10))
FR = EXP(MU_10+ETA(10))
K23 = Q/V2
K32 = Q/V3
S2 = V2
S4 = V4
S5 = V2
$DES
DADT(1) = -K12*A(1)+K41*A(4)
DADT(2) = K12*A(1)+K32*A(3)-(K23+K24)*A(2)
DADT(3) = K23*A(2)-K32*A(3)
DADT(4) = K24*A(2)-(K41+K40)*A(4)
DADT(5) = FR*(K12*A(1)+K32*A(3)-(K23+K24)*A(2))
$ERROR (ONLY OBSERVATIONS)
$EST METHOD=ITS NITER=50 PRINT=5 NSIG=3 SIGL=9
$EST METHOD=IMP NITER=500 ISAMPLE=1000 PRINT=1 NOABORT CTYPE=3
CITER=10 CALPHA=0.05
$THETA
$OMEGA
$SIGMA
Thanks for your help,
--
Hong Li, Ph.D.
Postdoctoral Associate
Department of Pharmaceutical Science
University at Buffalo, SUNY
363 Hochstetter Hall
Buffalo, NY 14260
Tel#: (716)645-4827