R示例：SEIR模型代码

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

不再提供，谢谢

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：SEI ggplot2 Library Solve gplot SEIR SEIR SEIR

这也太棒了吧！我空有matlab的相关代码，但完全不知道如何使用。看了博主的，想用R试试。非常感谢博主！

ppeeii 发表于 2020-10-2 14:57
这也太棒了吧！我空有matlab的相关代码，但完全不知道如何使用。看了博主的，想用R试试。非常感谢博主！

只是示例。实际应用需要更复杂的模型，考虑人口流动、社会网络等模型配合。单就SEIR模型而言，采用ODE就不如用DDE，但是DDE的模型我始终没有调好，由于变量之间的延迟微分误差的积累，曲线画到某个时点后就会崩溃。后来忙，就放弃DDE了

您好，请问触发事件这里的代码是什么意思，主要起什么作用呢？希望您不吝赐教，谢谢！

目百魔苦院 发表于 2021-4-20 19:43
您好，请问触发事件这里的代码是什么意思，主要起什么作用呢？希望您不吝赐教，谢谢！

var = c("S", "S", "S"),
  time = c(10, 20, 90),
  value = c(1, 1, 1),
  method = c("mult", "mult", "mult")
类似一个矩阵，三列是三个事件，比如第一个事件是在第10次计算时，变量S的取值乘以value中的1，method中的mult是乘的意思。示例代码这里实际上未做任何时间触发。

llb_321 发表于 2021-4-21 07:33
var = c("S", "S", "S"),
  time = c(10, 20, 90),
  value = c(1, 1, 1),

感谢！