library(Pkdesign)
paramList                      <-  list()
model                          <-  getExamplePkModel() 
paramList$pkParameters         <-  getParameters( model )
paramList$kineticsPlot         <-  plotMMKinetics( model , doseRange = c( 0 , 10 ) , concentrationRange = c(0,5) )

settings                       <-  list( NsimTimePoints = 100 ,  nSimSchemes = 1000 , nSubjectsForTimePoints = 3 , weightAUC = 50 , weightCmax = 50  )
paramList$settings             <-  settings 
paramList$optimalTimePoints    <-  "{ 1, 2 , 3.5 , 8}" 
paramList$optimalScheme        <-  getExampleSetOfSchemes()[ ,, 1]

Pharmacokinetic model

Input parameters

knitr::kable( paramList$pkParameters )
parameter value coeffVariation explanation
F 1.00 0.0 bioavailability
volumePlasma 10.00 0.2 volume of the central compartment (plasma)
Cld 15.00 0.2 distribution parameter between central and peripheral Compartment expressed in volume per time unit
volumeTissue 15.00 0.2 volume of the peripheral compartment (tissue)
VmaxAbsorption 5.00 0.2 maximum absorption rate ( absolute rate is rate per dose x dose)
kappaMMAbsorption 2.50 0.2 Michaelis-Menten constant for absorption
KaConstant NA 0.2 constant absorption rate per unit of dose (overrides Michealis-Menten kinetics)
VmaxClearance 30.00 0.2 maximum clearance rate (absolute rate is rate per concentration x
kappaMMClearance 0.25 0.2 concentration)
ClConstant NA 0.2 constant clearance rate ( overrides Michealis-Menten kinetics )

Model kinetics

plot( paramList$kineticsPlot )

Settings

setttings          <-  paramList$settings 
settingsTable      <-  data.frame(  Setting = names( settings ) , Value = unlist( settings ) )
rownames( settingsTable ) <-  NULL
knitr::kable( settingsTable )
Setting Value
NsimTimePoints 100
nSimSchemes 1000
nSubjectsForTimePoints 3
weightAUC 50
weightCmax 50

Optimal design

Time points

paramList$optimalTimePoints
## [1] "{ 1, 2 , 3.5 , 8}"

Scheme

paramList$optimalScheme
##      timePoint1 timePoint2 timePoint3 timePoint4
## rat1       TRUE      FALSE       TRUE       TRUE
## rat2      FALSE       TRUE      FALSE       TRUE