Question about error message regarding invalid MDV field
From: zexun zhou zexunzho@usc.edu
Subject: [NMusers] Question about error message regarding invalid MDV field
Date: Tue, 06 Sep 2005 14:06:00 -0700
NONMEM users,
I have a following control file for a two compartment PK model with i.v. infusions for a period of 1 hour
(3 hours for some sujects) every 24 hours for 5 consecutive doses. and the plasma concentrations are
measured immediately prior to the last dose and after the last dose.
$PROB 2 CMPT POPULATION PK MODEL
$INPUT ID TIME AMT RATE ADDL II EVID CMT DV MDV
$DATA TEST.CSV IGNORE=C
$SUBROUTINES ADVAN3
$PK
K=THETA(1)*EXP(ETA(1))
V=THETA(2)*EXP(ETA(2))
K12=THETA(3)*EXP(ETA(3))
K21=THETA(4)*EXP(ETA(4))
S2=V
$THETA
(0, 500) ;[KE]
(0, 200, 2300) ;[V]
(0, 0.065, 10) ;[K12]
(0, 0.056, 10) ;[K21]
$ERROR
IPRED=F
IRES=DV-IPRED
IWRES=IRES/F
Y=F* (1+ EPS(1)) +EPS(2)
$OMEGA DIAG(4)
0.16 ;[KE]
0.16 ;[V]
0.49 ;[K21]
0.49 ;[K12]
$SIGMA
0.13 ;[P]
0.02 ;[A]
$EST SIGDIGITS=4 MAXEVAL=9999 PRINT=20 METHOD=1 POSTHOC
$COV
$TABLE ID TIME DV IPRED IWRES NOPRINT ONEHEADER FILE=xxx.TAB
$TABLE ID K V K12 K21 NOPRINT ONEHEADER FILE=xxx1.TAB
and the data file looks likes this
C ID TIME AMT RATE ADDL II EVID CMT DV MDV
1 0 1309 436 4 24 1 2 . 1
1 95.9 0 0 . . 0 2 0.1552822 0
1 98 0 0 . . 0 2 0.124313 0
1 100 0 0 . . 0 2 6.468202 0
1 102 0 0 . . 0 2 4.24409 0
1 108 0 0 . . 0 2 1.330367 0
1 120 0 0 . . 0 2 0.1491756 0
2 0 1309 1309 4 24 1 2 . 1
2 95.9 0 0 . . 0 2 0.1474309 0
2 98 0 0 . . 0 2 5.74108 0
etc.
running the control file gives me these error messages
RECORD 2, DATA ITEM 10: 0.1552822
INVALID MDV FIELD.
etc.
I altered the code where I input every dosing event instead of using the ADDL and II data.
$PROB POPULATION PK MODEL
$INPUT ID TIME AMT RATE EVID CMT DV
$DATA TEST1.CSV IGNORE=C
$SUBROUTINES ADVAN4
$PK
K=THETA(1)*EXP(ETA(1))
V=THETA(2)*EXP(ETA(2))
K12=THETA(3)*EXP(ETA(3))
K21=THETA(4)*EXP(ETA(4))
S2=V
$THETA
(0, 500) ;[KE]
(0, 200, 2300) ;[V]
(0, 0.065, 10) ;[K12]
(0, 0.056, 10) ;[K21]
$ERROR
IPRED=F
IRES=DV-IPRED
IWRES=IRES/F
Y=F* (1+ EPS(1)) +EPS(2)
$OMEGA DIAG(4)
0.16 ;[KE]
0.16 ;[V]
0.49 ;[K21]
0.49 ;[K12]
$SIGMA
0.13 ;[P]
0.02 ;[A]
$EST SIGDIGITS=4 MAXEVAL=9999 PRINT=20 METHOD=1 POSTHOC
$COV
$TABLE ID TIME DV IPRED IWRES NOPRINT ONEHEADER FILE=xxx.TAB
$TABLE ID K V K12 K21 NOPRINT ONEHEADER FILE=xxx1.TAB
and the data file was as follows
C ID TIME AMT RATE EVID CMT DV
1 0 1309 436 1 2 .
1 24 1309 436 1 2 .
1 48 1309 436 1 2 .
1 72 1309 436 1 2 .
1 95.9 0 0 0 2 0.155
1 96 1309 436 1 2 .
1 98 0 0 0 2 0.124
1 100 0 0 0 2 6.468
1 102 0 0 0 2 4.244
1 108 0 0 0 2 1.33
1 120 0 0 0 2 0.149
2 0 1309 1309 1 2 .
2 24 1309 1309 1 2 .
2 48 1309 1309 1 2 .
2 72 1309 1309 1 2 .
etc.
and running the program, NONMEM gave me the following error message,
ROCORD 1, DATA ITEM 7:2
THE DV DATA ITEM IS POSITVE, BUT THE MDV DATA ITEM IS 1
etc.
What exactly is MDV? and I thought I didn't need to secify MDV for the latter set of control and data files.
What is exactly wrong with my code?
Thanks for your time and any help would be greatly appreciated.
Susan