seirah_solve: SEIRAH solve function

In sistm/SEIRcovid19: Interface Monolix for multilevel modeling of epidemic dynamics

Description

SEIRAH solve function

Usage

seirah_solve(
  init,
  t,
  par,
  pred = TRUE,
  ns1all = NULL,
  ns2all = NULL,
  typecov = "scale"
)

Examples

# Set parameter values
b<-1
r<-1
alpha<-1
De<-5.2
Di<-2.3
Dq<-2
Dh<-30
popSize<-10000000
dailyMove<-0

init<-c(9999467, 346, 80, 0, 80, 27)
t<-seq(0,365)
par <- c(1.75, 0.41, alpha, De, Di, Dq, Dh, popSize, dailyMove,30,30,0,3000000)

# Solve system using lsoda
sol <- seirah_solve(init, t, par)
plot(sol)

# Plot solution
plot(sol$time,sol$I,ylim=c(0,100))

sistm/SEIRcovid19 documentation built on Aug. 9, 2021, 7:37 a.m.