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)
if (FALSE) {
plot(out, "CP RESP, GUT")
}
Run the code above in your browser using DataLab