xt outside of design space when having 0 and using discrete_xt

[zparra]

It seems that when you specify a value of 0 as minimum of your design space, it gets replaced by 1e-05

poped.db$design_space$minxt (done with the vignette example)
obs_1 obs_2 obs_3 obs_4 obs_5
grp_1 1e-05 1e-05 1e-05 240 240
grp_2 1e-05 1e-05 1e-05 240 240

However if in your initial design you have a time 0, it does not get replaced

poped.db$design$xt
obs_1 obs_2 obs_3 obs_4 obs_5
grp_1 0 2 8 240 245
grp_2 0 2 8 240 245

Apparently, this is not a problem apparently for a continous space, but if instead you provide a discrete space, you get an error because the obs_1 is out of the design space:

poped.db.discrete <- create.poped.database(poped.db,discrete_xt = list(0:248))
Error in eval(substitute(expr), data, enclos = parent.frame()) :
xt value for group 1 (column 1) is not in the design space

[andrewhooker]

Hi Thanks for the example. You are right, this does not seem to be logical. I was trying to avoid situations where numerical problem occur (oral PK data at time zero should give a prediction of zero which causes numerical issues in the calculation of the FIM), but didn’t get it quite right. You can avoid this problem by specifying the ourzero option in create.poped.database(). Best regards, Andy

…

On 25 Aug 2017, at 10:10 , zparra ***@***.***> wrote: It seems that when you specify a value of 0 as minimum of your design space, it gets replaced by 1e-05 poped.db$design_space$minxt (done with the vignette example) obs_1 obs_2 obs_3 obs_4 obs_5 grp_1 1e-05 1e-05 1e-05 240 240 grp_2 1e-05 1e-05 1e-05 240 240 However if in your initial design you have a time 0, it does not get replaced poped.db$design$xt obs_1 obs_2 obs_3 obs_4 obs_5 grp_1 0 2 8 240 245 grp_2 0 2 8 240 245 Apparently, this is not a problem apparently for a continous space, but if instead you provide a discrete space, you get an error because the obs_1 is out of the design space: poped.db.discrete <- create.poped.database(poped.db,discrete_xt = list(0:248)) Error in eval(substitute(expr), data, enclos = parent.frame()) : xt value for group 1 (column 1) is not in the design space — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#13>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AFLmZE-_R6aJCv-bbrr8VwNDgDbr45USks5sboF4gaJpZM4PCWaB>.