seir模型数学建模python_R示例：SEIR模型代码

说明：本例对常见的SEIR模型作了改变，由于DDE模型始终存在缺陷，故仍采用ODE建模。

实现环境：R 3.6.3+deSolve+ggplot2。

#-------------------------------------

#基于微分方程组建模

#dS/dt=-b_1*S*I/N-b_2*S*E/N

#dE/dt=b_1*S*I/N+b_2*S*E/N-k*E

#dI/dt=k*E-z*I

#dH/dt=z*I-d*H-m*H

#dR/dt=m*H

#dD/dt=d*H

#变量

#S易感人群

#E暴露潜伏

#I发病

#H入院

#R治愈

#D病故

#参数

#b_1,b_2接触感染概率

#k潜伏到发病的转化率

#z发病入院比例

#m治愈率

#d感染期死亡率

#--------------------------------------

library(deSolve)

library(ggplot2)

#模型

seir

with(as.list(c(state, pars)), {

dS

dE

dI

dH

dR

dD

dN

list(c(dS, dE, dI, dH, dR, dD, dN))

})

}

#参数

N