Re: Use of ACCEPT in $DATA

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