# NOT RUN {
mod <- mrgsolve::house(end=48, delta=0.2) %>% init(GUT=1000)
out <- mrgsim(mod)
plot(out)
plot(out, subset=time <=24)
plot(out, GUT+CP~.)
plot(out, CP+RESP~time, col="black", scales="same", lty=2)
# }
# NOT RUN {
plot(out, "CP RESP, GUT")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab