Starting with version 0.10.0 (released October, 2019), mrgsolve can
implement dosing records calling for a zero-order infusion at steady
state. To do this, set rate
to a non-zero value,
amt
to zero, and ss
to 1.
For a 10 mg/hr infusion at steady state
library(mrgsolve)
library(dplyr)
<- house(end = 72, delta = 0.1) %>% update(outvars = "CP")
mod
<- param(mod, CL = 2)
mod
<- ev(amt = 0, rate = 100, ss = 1, cmt = 2) dose
The dosing event looks like
dose
. Events:
. time amt rate cmt evid ss
. 1 0 0 100 2 1 1
With CL
set to 2, we expect the system to start out at
50
mrgsim(mod,dose,recsort = 3) %>% plot()
rate
to -1 in the input dataR_CMT
in [main]
rate
to -2 in the input dataD_CMT
in [main]
mrgsolve: mrgsolve.github.io | metrum research group: metrumrg.com