mod <- mrgsolve::house() %>% ev(amt = 100)
out <- mrgsim(mod)
out_df <- dplyr::mutate(out, time <= 72)
plot(out)
plot_sims(out, CP, RESP)
if (FALSE) {
plot_sims(out, .f = ~ CP + RESP)
plot_sims(out, .f = CP + RESP ~ time)
}
Run the code above in your browser using DataLab