Deterministic SIR

1 collaborator

Olivier Restif (Author)

WHAT IS IT?

(a general understanding of what the model is trying to show or explain)

HOW IT WORKS

(what rules the agents use to create the overall behavior of the model)

HOW TO USE IT

(how to use the model, including a description of each of the items in the Interface tab)

THINGS TO NOTICE

(suggested things for the user to notice while running the model)

THINGS TO TRY

(suggested things for the user to try to do (move sliders, switches, etc.) with the model)

EXTENDING THE MODEL

(suggested things to add or change in the Code tab to make the model more complicated, detailed, accurate, etc.)

NETLOGO FEATURES

(interesting or unusual features of NetLogo that the model uses, particularly in the Code tab; or where workarounds were needed for missing features)

RELATED MODELS

(models in the NetLogo Models Library and elsewhere which are of related interest)

CREDITS AND REFERENCES

(a reference to the model's URL on the web if it has one, as well as any other necessary credits, citations, and links)

Comments and Questions

Please start the discussion about this model! (You'll first need to log in.)

Click to Run Model

;; System dynamics model globals
globals [
  R
  I
  S
  dt
]

;; Initializes the system dynamics model.
;; Call this in your model's SETUP procedure.

to setup
  ca
  set dt 0.001
  ;; initialize stock values
  set R 0
  set I I_0
  set S 1 - I_0
  reset-ticks
end 


;; ===========================================================================================
;; Step through the system dynamics model by performing next iteration of Euler's method.
;; Call this in your model's GO procedure.

to go
  if  ticks >= t_max  [ stop ]
  do-SIR-step
  do-plot
  tick-advance dt
end 

;; ===========================================================================================

to do-SIR-step
  let dS ( bd * (1 - vac) - transmission * S * I - bd * S )
  let dI ( transmission * S * I - recovery * I - bd * I )
  let dR ( bd * vac + recovery * I - bd * R )
  set S ( S + dS * dt )
  set I ( I + dI * dt )
  set R ( R + dR * dt )
end 


;; =========================================================================================
;; Plot the current state of the system dynamics model's stocks

to do-plot
  if plot-pen-exists? "R" [
    set-current-plot-pen "R"
    plotxy ticks R
  ]
  if plot-pen-exists? "I" [
    set-current-plot-pen "I"
    plotxy ticks I
  ]
  if plot-pen-exists? "S" [
    set-current-plot-pen "S"
    plotxy ticks S
  ]
end

There is only one version of this model, created almost 8 years ago by Olivier Restif.

Attached files

File	Type	Description	Last updated
Deterministic SIR.png	preview	Preview for 'Deterministic SIR'	almost 8 years ago, by Olivier Restif	Download

This model does not have any ancestors.

This model does not have any descendants.

NetLogo