Hamiltonian Analysis of Subcritical Stochastic Epidemic Dynamics

<div><i>Transient dynamics of probability density</i> in the subcritical SIS system, displayed as <svg height="15.6876pt" id="M188" style="vertical-align:-3.4294pt" version="1.1" viewbox="-0.0498162 -12.2582 98.3347 15.6876" width="98.3347pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M542 242C542 369 469 429 379 468C370 472 372 496 376 513L416 681L400 691L344 648L219 51C167 69 113 128 113 217C113 324 166 387 266 428L252 459C124 418 23 336 23 202C23 95 90 13 205 -8L182 -132C165 -222 191 -254 209 -261L215 -258L270 -7C415 -4 542 98 542 242ZM469 232C469 126 403 33 279 39L355 409C399 394 469 331 469 232Z" id="g113-243"></path></g><g transform="matrix(.013,0,0,-0.013,7.346,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z" id="g113-41"></path></g><g transform="matrix(.013,0,0,-0.013,11.844,0)"><path d="M536 404C536 423 520 448 491 448C445 448 398 404 308 283L286 338C255 416 242 448 217 448C182 448 138 402 92 341L111 321C149 368 169 378 178 378C188 378 198 363 214 320L254 212C185 117 138 65 107 65C98 65 85 69 82 75C79 82 73 86 65 86C44 86 23 60 23 39C23 7 44 -12 71 -12C119 -12 168 33 265 177L306 61C321 17 347 -12 373 -12C413 -12 465 33 507 96L491 119C459 84 432 60 413 60C395 60 378 92 358 148L321 250C341 279 369 310 389 332C417 363 439 382 456 382C466 382 475 376 481 368C486 361 492 358 496 358C513 358 536 381 536 404Z" id="g113-121"></path></g><g transform="matrix(.013,0,0,-0.013,19.112,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path></g><g transform="matrix(.013,0,0,-0.013,24.255,0)"><path d="M324 430H196L233 583L223 592L145 529L120 430H54L29 396L31 388H111L56 126C33 15 54 -12 77 -12C137 -12 214 57 250 95L233 119C208 92 155 59 138 59C126 59 120 70 131 125L186 390L298 394L324 430Z" id="g113-117"></path></g><g transform="matrix(.013,0,0,-0.013,28.688,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z" id="g113-42"></path></g><g transform="matrix(.013,0,0,-0.013,36.818,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path></g><g transform="matrix(.013,0,0,-0.013,48.08,0)"><path d="M645 631C614 643 545 666 457 666C215 666 23 519 23 283C23 90 158 -16 337 -16C412 -16 489 2 522 10C543 39 590 127 606 167L580 181C519 89 459 18 348 18C201 18 122 136 122 287C122 464 244 632 435 632C544 632 602 595 608 472L639 475C643 526 645 581 645 631Z" id="g113-68"></path></g><g transform="matrix(.013,0,0,-0.013,56.765,0)"><path d="M391 364C391 409 353 448 295 448C249 448 198 426 152 393C65 331 23 225 23 139C23 14 96 -12 146 -12C198 -12 280 9 367 101L351 124C300 78 242 48 194 48C129 48 109 107 109 162V191C208 213 391 266 391 364ZM313 350C313 305 268 261 113 223C132 334 187 381 217 398C227 404 244 405 261 405C290 405 313 385 313 350Z" id="g113-102"></path></g><g transform="matrix(.0091,0,0,-0.0091,62.147,-5.741)"><path d="M556 236V289H56V236H556Z" id="g54-33"></path></g><g transform="matrix(.0091,0,0,-0.0091,67.707,-5.741)"><path d="M870 650H622L613 619C694 615 708 606 708 559C708 531 702 477 686 396L636 136H633L281 650H127L120 619C180 616 198 611 212 592C227 572 232 558 220 503L166 254C148 169 139 131 126 89C112 43 85 35 33 32L24 0H280L286 32C208 37 189 44 189 93C189 124 198 179 212 250L267 526H271L622 -8H655L732 396C747 477 762 532 770 556C785 600 798 616 862 619L870 650Z" id="g50-79"></path></g><g transform="matrix(.0091,0,0,-0.0091,75.38,-5.741)"><path d="M528 97L513 125C484 95 440 68 431 68C425 68 421 76 427 107C449 227 476 342 505 451H495L418 424L385 278C338 191 220 62 157 62C145 62 133 76 145 130L200 368C215 434 210 451 183 451C157 451 92 418 24 352L37 324C78 356 106 377 115 377C123 377 122 364 116 340L64 115C58 89 56 68 56 52C56 2 78 -12 103 -12C127 -12 159 -4 190 18C247 59 309 113 366 187H368L351 108C328 1 348 -12 368 -12C403 -12 472 33 528 97Z" id="g50-118"></path></g><g transform="matrix(.0091,0,0,-0.0091,80.404,-5.741)"><path d="M309 -142C211 -61 151 99 151 271S210 601 309 683L288 710C141 611 75 451 75 272V271C75 90 141 -72 288 -169L309 -142Z" id="g50-41"></path></g><g transform="matrix(.0091,0,0,-0.0091,83.661,-5.741)"><path d="M545 403C545 425 527 451 498 451C448 451 400 404 312 286L290 341C262 412 246 451 219 451C182 451 138 404 92 345L112 323C152 368 169 378 181 378C192 378 201 359 216 321L257 214C184 115 142 69 113 69C102 69 93 73 88 79S78 89 68 89C47 89 24 61 24 40C24 8 49 -12 75 -12C125 -12 175 31 271 175L312 64C330 15 357 -12 382 -12C421 -12 475 35 515 98L496 121C466 88 439 62 420 62C402 62 384 92 363 147L325 247C349 280 374 309 397 333C422 361 446 381 462 381C471 381 481 374 488 366C493 360 499 357 505 357C521 357 545 380 545 403Z" id="g50-121"></path></g><g transform="matrix(.0091,0,0,-0.0091,88.84,-5.741)"><path d="M98 134C72 134 46 117 46 90C46 73 55 65 60 64C95 55 124 32 124 -4C124 -42 95 -68 44 -89L57 -123C122 -104 194 -60 194 22C194 94 136 134 98 134Z" id="g50-45"></path></g><g transform="matrix(.0091,0,0,-0.0091,90.996,-5.741)"><path d="M329 433H203L239 587L230 596L147 534L123 433H57L30 395L34 388H115L61 129C37 16 59 -12 85 -12C147 -12 222 58 260 98L241 125C212 95 160 62 144 62C132 62 127 71 138 126L192 386L305 394L329 433Z" id="g50-117"></path></g><g transform="matrix(.0091,0,0,-0.0091,94.154,-5.741)"><path d="M283 271C283 451 220 610 70 710L48 683C146 603 207 443 207 271S149 -59 48 -142L70 -169C220 -67 283 90 283 271Z" id="g50-42"></path></g></svg> versus <svg height="6.1673pt" id="M189" style="vertical-align:-0.2063904pt" version="1.1" viewbox="-0.0498162 -5.96091 7.39387 6.1673" width="7.39387pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M536 404C536 423 520 448 491 448C445 448 398 404 308 283L286 338C255 416 242 448 217 448C182 448 138 402 92 341L111 321C149 368 169 378 178 378C188 378 198 363 214 320L254 212C185 117 138 65 107 65C98 65 85 69 82 75C79 82 73 86 65 86C44 86 23 60 23 39C23 7 44 -12 71 -12C119 -12 168 33 265 177L306 61C321 17 347 -12 373 -12C413 -12 465 33 507 96L491 119C459 84 432 60 413 60C395 60 378 92 358 148L321 250C341 279 369 310 389 332C417 363 439 382 456 382C466 382 475 376 481 368C486 361 492 358 496 358C513 358 536 381 536 404Z" id="g113-121"></path></g></svg> using the same data points as in Figure <a href="../fig4/">4</a>, with <svg height="8.8423pt" id="M190" style="vertical-align:-0.2064009pt" version="1.1" viewbox="-0.0498162 -8.6359 44.761 8.8423" width="44.761pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M822 650H589L583 622C660 617 677 607 674 561C672 534 664 481 647 390L600 137H596L273 650H126L120 622C176 620 194 615 207 594C221 571 225 557 214 504L161 257C141 166 129 112 121 85C108 42 83 30 29 28L23 0H260L266 28C193 33 173 42 176 89C178 122 186 172 202 255L256 527H259L583 -8H612L690 390C708 481 720 535 728 558C744 603 756 619 816 622L822 650Z" id="g113-79"></path></g><g transform="matrix(.013,0,0,-0.013,14.542,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path></g><g transform="matrix(.013,0,0,-0.013,25.805,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path></g><g transform="matrix(.013,0,0,-0.013,32.045,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path></g><g transform="matrix(.013,0,0,-0.013,38.285,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path></g></svg>. Each curve is normalized to total probability one. The quasistationary distribution (<a href="https://static.hindawi.com/articles/cmmm/volume-2017/4253167/figures/#EEq9">23</a>) is plotted in gray.</div>

Computational and Mathematical Methods in Medicine

Hamiltonian Analysis of Subcritical Stochastic Epidemic Dynamics

Figure 5