NONMEM Users Network Archive

Hosted by Cognigen

[NMusers] Use of ACCEPT in $DATA

From: Dennis Fisher <fisher_at_plessthan.com>
Date: Thu, 24 Aug 2017 15:16:23 -0700

NONMEM 7.4.1

Colleagues

I am trying to use the ACCEPT option in $DATA in order to select a =
subset of records (to evaluate the impact of the # of samples/subject on =
confidence intervals).

I used the following code:
        ACCEPT=(TIME=0, TIME=1, TIME=2, TIME=4, TIME=6, =
TIME=24)

NMTRAN then creates a dataset but — to my surprise — =
TIME=6 is not in the dataset (all the others are).

I am copying the first few rows of the input dataset so that you can see =
what is being provided to NMTRAN:

=
ID,AGE,MONTHS,SEX,WT,AMT,RATE,TIME,EVID,MDV,REPLICATE,IPRED,CWRES,DV,PRED,=
RES,WRES
1101,12,144,1,30.054,210.38,841.51,0,1,1,1,0,0,0,0,0,0
1101,12,144,1,30.054,0,0,1,0,0,1,187.42,0,179.28,199.26,-19.979,0
1101,12,144,1,30.054,0,0,2,0,0,1,180.92,0,187.92,194.09,-6.1659,0
1101,12,144,1,30.054,0,0,4,0,0,1,169.84,0,177.66,184.37,-6.712,0
1101,12,144,1,30.054,0,0,6,0,0,1,160.61,0,153.43,175.39,-21.96,0

The underlined / boldfaced value (6) in the final row is the problem.

I assume that NMTRAN is reading that value as something other than 6.0 =
(e.g., 6.000001) and thereby omitting it.

I have reviewed NMHELP to see if there is some other way to accomplish =
this. Ideally, there would be something like:
        TIME.GT.5.9.AND.TIME.LT.6.1
but that does not appear to be supported.

The alternative is to modify the dataset to include many possible =
MDV/EVID columns. However, it would be more elegant to do this in the =
control stream.
Or, if there is some way to find out the exact value that NMTRAN sees, I =
could specify that value.

Any help would be appreciated.

Dennis

Dennis Fisher MD
P < (The "P Less Than" Company)
Phone / Fax: 1-866-PLessThan (1-866-753-7784)
www.PLessThan.com <http://www.plessthan.com/>






Received on Thu Aug 24 2017 - 18:16:23 EDT

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request_at_iconplc.com. Once subscribed, you may contribute to the discussion by emailing: nmusers_at_globomaxnm.com.