Plotting optimized output seems to used wrong times

[krimobe]

Plotting optimized output seems to used wrong times

[andrewhooker]

Hi.

It seems the ff.ode.compiled function is a bit outdated and is doing a strange filtering. I would change it to

ff.ODE.compiled <- function(model_switch,xt,parameters,poped.db){
  with(as.list(parameters),{
    A_ini <- c(A1 = DOSE*Favail, A2 = 0)
    times_xt <- drop(xt) #xt[,,drop=T] 
    times <- sort(times_xt)
    times <- c(0,times) ## add extra time for start of integration
    out <- ode(A_ini, times, func = "derivs", parms = c(CL,V,KA),
               #jacfunc = "jac", # not really needed, speed up is minimal if this is defined or not.
               dllname = "one_comp_oral_CL",
               initfunc = "initmod", nout = 1, outnames = "Sum")
    y = out[, "A2"]/(V)
    y=y[match(times_xt,out[,"time"])]
    y = cbind(y) # must be a column matrix
    return(list( y= y,poped.db=poped.db))
  })
} 

Which gives the plot after optimization of :