predict.pkmod | R Documentation |
Predict concentrations from a pkmod object - can be a user defined function
## S3 method for class 'pkmod' predict(object, inf, tms, return_times = FALSE, ...)
object |
An object with class pkmod. |
inf |
A matrix with columns "begin","end","inf_rate" indicating when infusions should be administered. Can be created by 'inf_manual' or 'inf_tci'. |
tms |
Times at which to calculate predicted concentrations. |
return_times |
Logical. Should prediction times be returned along with responses? Defaults to FALSE. |
... |
List or vector of values to be passed on to update.pkmod. |
Matrix of predicted concentrations associated with a pkmod object and and infusion schedule.
# dosing schedule dose <- inf_manual(inf_tms = c(0,0.5,4,4.5,10), inf_rate = c(100,0,80,0,0)) # pkmod object my_mod <- pkmod(pars_pk = c(CL = 15, V1 = 10, Q2 = 10, V2 = 20)) # predict at specific times predict(my_mod, inf = dose, tms = c(1.5,2.5,3)) # predict with an Emax PD function my_mod_pd <- pkmod(pars_pk = c(v1 = 8.995, v2 = 17.297, v3 = 120.963, cl = 1.382, q2 = 0.919, q3 = 0.609, ke0 = 1.289), pars_pd = c(c50 = 2.8, gamma = 1.47, gamma2 = 1.89, e0 = 93, emx = 93), pdfn = emax, pdinv = emax_inv, ecmpt = 4) predict(my_mod_pd, inf = dose, tms = c(1.5,2.5,3)) # predict with a subset of new PK-PD parameters predict(my_mod_pd, inf = dose, tms = c(1.5,2.5,3), pars_pk = c(ke0 = 0.8), pars_pd = c(c50 = 2, e0 = 100))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.