NONMEM Users Network Archive

Hosted by Cognigen

Re: [NMusers] Use of ACCEPT in $DATA

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

Colleagues

My apologies — my question reflected a (complicated) error on my =
part. NMTRAN / NONMEM were doing the correct thing.
Of note, Leonid proposed using .EQN — that does not appear to be =
necessary.

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/>




> On Aug 24, 2017, at 3:40 PM, Leonid Gibiansky =
<lgibiansky_at_quantpharm.com> wrote:
>
> from the manual:
>
> "When the IGNORE option is used to filter records from the input file, =
the .EQ., =, .NE., and /= symbols perform literal string =
comparisons. To provide a numerical equality comparison, use .EQN. for =
numerical equals, and .NEN. for numerical not equals."
>
> May be there is a space there or something, try .EQN.
>
> Leonid
>
>
> On 8/24/2017 6:16 PM, Dennis Fisher wrote:
>> 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,PRE=
D,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/> =
<http://www.plessthan.com/ <http://www.plessthan.com/>>



Received on Thu Aug 24 2017 - 22:04:15 EDT

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