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
Question about error message regarding invalid MDV field
From: "Serge Guzy" GUZY@xoma.com
Subject: RE: [NMusers] Question about error message regarding invalid MDV field
Date: Tue, 6 Sep 2005 14:57:39 -0700
Is it possible that you forgot to have the first column blank.
It seems like the program is reading DV values for MDV.
Serge Guzy
President POP_PHARM
From: "Bachman, William (MYD)" bachmanw@iconus.com
Subject: RE: [NMusers] Question about error message regarding invalid
Date: Wed, 7 Sep 2005 10:04:20 -0400
As Serge suggests, I think you have a blank first column in your .csv file
that doesn't show in your email snippet of the data file because:
1. "0.1552822" should be item 9 not 10
2. if I run your control stream with a data file I make from the snippet, it
runs without error.
so you either need to remove the blank first column or add a dummy item to
the $INPUT line to shift the reading of data items, e.g. (where C is the
dummy item):
$INPUT C ID TIME AMT RATE ADDL II EVID CMT DV MDV
_______________________________________________________