Función do({}) en R no me permite nombrar a ode para poder añadir las variables dentro de mapings de ggplot

Question

Good afternoon.

I am implementing SIR model in R, I need to vary beta and gamma for it.

library(deSolve)

par(mar = rep(2, 4))
N = 1000
vi <- c(S =