|
51 | 51 | { |
52 | 52 | "cell_type": "code", |
53 | 53 | "execution_count": 2, |
54 | | - "metadata": {}, |
| 54 | + "metadata": { |
| 55 | + "collapsed": true |
| 56 | + }, |
55 | 57 | "outputs": [], |
56 | 58 | "source": [ |
57 | 59 | "init = State(S=89, I=1, R=0)\n", |
|
68 | 70 | { |
69 | 71 | "cell_type": "code", |
70 | 72 | "execution_count": 3, |
71 | | - "metadata": {}, |
| 73 | + "metadata": { |
| 74 | + "collapsed": true |
| 75 | + }, |
72 | 76 | "outputs": [], |
73 | 77 | "source": [ |
74 | 78 | "init /= sum(init)\n", |
|
177 | 181 | { |
178 | 182 | "cell_type": "code", |
179 | 183 | "execution_count": 7, |
180 | | - "metadata": {}, |
| 184 | + "metadata": { |
| 185 | + "collapsed": true |
| 186 | + }, |
181 | 187 | "outputs": [], |
182 | 188 | "source": [ |
183 | 189 | "state = update1(init, system)\n", |
|
224 | 230 | "cell_type": "code", |
225 | 231 | "execution_count": 9, |
226 | 232 | "metadata": { |
| 233 | + "collapsed": true, |
227 | 234 | "scrolled": true |
228 | 235 | }, |
229 | 236 | "outputs": [], |
|
243 | 250 | { |
244 | 251 | "cell_type": "code", |
245 | 252 | "execution_count": 10, |
246 | | - "metadata": {}, |
| 253 | + "metadata": { |
| 254 | + "collapsed": true |
| 255 | + }, |
247 | 256 | "outputs": [], |
248 | 257 | "source": [ |
249 | 258 | "# Solution goes here" |
|
306 | 315 | { |
307 | 316 | "cell_type": "code", |
308 | 317 | "execution_count": 12, |
309 | | - "metadata": {}, |
| 318 | + "metadata": { |
| 319 | + "collapsed": true |
| 320 | + }, |
310 | 321 | "outputs": [], |
311 | 322 | "source": [ |
312 | 323 | "tc = 3 # time between contacts in days \n", |
|
358 | 369 | { |
359 | 370 | "cell_type": "code", |
360 | 371 | "execution_count": 14, |
361 | | - "metadata": {}, |
| 372 | + "metadata": { |
| 373 | + "collapsed": true |
| 374 | + }, |
362 | 375 | "outputs": [], |
363 | 376 | "source": [ |
364 | 377 | "plot_results(system.S, system.I, system.R)\n", |
|
416 | 429 | { |
417 | 430 | "cell_type": "code", |
418 | 431 | "execution_count": 16, |
419 | | - "metadata": {}, |
| 432 | + "metadata": { |
| 433 | + "collapsed": true |
| 434 | + }, |
420 | 435 | "outputs": [], |
421 | 436 | "source": [ |
422 | 437 | "tc = 3 # time between contacts in days \n", |
|
425 | 440 | "beta = 1 / tc # contact rate in per day\n", |
426 | 441 | "gamma = 1 / tr # recovery rate in per day\n", |
427 | 442 | "\n", |
428 | | - "sir = make_system(beta, gamma)\n", |
| 443 | + "system = make_system(beta, gamma)\n", |
429 | 444 | "run_simulation(system, update1)\n", |
430 | 445 | "system.results.head()" |
431 | 446 | ] |
|
440 | 455 | { |
441 | 456 | "cell_type": "code", |
442 | 457 | "execution_count": 17, |
443 | | - "metadata": {}, |
| 458 | + "metadata": { |
| 459 | + "collapsed": true |
| 460 | + }, |
444 | 461 | "outputs": [], |
445 | 462 | "source": [ |
446 | 463 | "frame = system.results\n", |
|
457 | 474 | { |
458 | 475 | "cell_type": "code", |
459 | 476 | "execution_count": 18, |
460 | | - "metadata": {}, |
| 477 | + "metadata": { |
| 478 | + "collapsed": true |
| 479 | + }, |
461 | 480 | "outputs": [], |
462 | 481 | "source": [ |
463 | 482 | "# Solution goes here" |
|
506 | 525 | { |
507 | 526 | "cell_type": "code", |
508 | 527 | "execution_count": 20, |
509 | | - "metadata": {}, |
| 528 | + "metadata": { |
| 529 | + "collapsed": true |
| 530 | + }, |
510 | 531 | "outputs": [], |
511 | 532 | "source": [ |
512 | 533 | "system.beta = 0.333\n", |
|
553 | 574 | "cell_type": "code", |
554 | 575 | "execution_count": 22, |
555 | 576 | "metadata": { |
| 577 | + "collapsed": true, |
556 | 578 | "scrolled": true |
557 | 579 | }, |
558 | 580 | "outputs": [], |
|
575 | 597 | "cell_type": "code", |
576 | 598 | "execution_count": 24, |
577 | 599 | "metadata": { |
| 600 | + "collapsed": true, |
578 | 601 | "scrolled": true |
579 | 602 | }, |
580 | 603 | "outputs": [], |
|
596 | 619 | { |
597 | 620 | "cell_type": "code", |
598 | 621 | "execution_count": 26, |
599 | | - "metadata": {}, |
| 622 | + "metadata": { |
| 623 | + "collapsed": true |
| 624 | + }, |
600 | 625 | "outputs": [], |
601 | 626 | "source": [ |
602 | 627 | "# Solution goes here" |
|
646 | 671 | { |
647 | 672 | "cell_type": "code", |
648 | 673 | "execution_count": 28, |
649 | | - "metadata": {}, |
| 674 | + "metadata": { |
| 675 | + "collapsed": true |
| 676 | + }, |
650 | 677 | "outputs": [], |
651 | 678 | "source": [ |
652 | 679 | "tc = 3 # time between contacts in days \n", |
|
669 | 696 | { |
670 | 697 | "cell_type": "code", |
671 | 698 | "execution_count": 29, |
672 | | - "metadata": {}, |
| 699 | + "metadata": { |
| 700 | + "collapsed": true |
| 701 | + }, |
673 | 702 | "outputs": [], |
674 | 703 | "source": [ |
675 | 704 | "run_simulation(system, update1)\n", |
|
686 | 715 | { |
687 | 716 | "cell_type": "code", |
688 | 717 | "execution_count": 30, |
689 | | - "metadata": {}, |
| 718 | + "metadata": { |
| 719 | + "collapsed": true |
| 720 | + }, |
690 | 721 | "outputs": [], |
691 | 722 | "source": [ |
692 | 723 | "system2 = make_system(beta, gamma)\n", |
|
707 | 738 | { |
708 | 739 | "cell_type": "code", |
709 | 740 | "execution_count": 31, |
710 | | - "metadata": {}, |
| 741 | + "metadata": { |
| 742 | + "collapsed": true |
| 743 | + }, |
711 | 744 | "outputs": [], |
712 | 745 | "source": [ |
713 | 746 | "plot(system.results.S, '-', label='No immunization')\n", |
|
729 | 762 | { |
730 | 763 | "cell_type": "code", |
731 | 764 | "execution_count": 32, |
732 | | - "metadata": {}, |
| 765 | + "metadata": { |
| 766 | + "collapsed": true |
| 767 | + }, |
733 | 768 | "outputs": [], |
734 | 769 | "source": [ |
735 | 770 | "immunize_array = linspace(0, 1, 11)\n", |
|
800 | 835 | { |
801 | 836 | "cell_type": "code", |
802 | 837 | "execution_count": 35, |
803 | | - "metadata": {}, |
| 838 | + "metadata": { |
| 839 | + "collapsed": true |
| 840 | + }, |
804 | 841 | "outputs": [], |
805 | 842 | "source": [ |
806 | 843 | "plot(infected_sweep)\n", |
|
870 | 907 | { |
871 | 908 | "cell_type": "code", |
872 | 909 | "execution_count": 37, |
873 | | - "metadata": {}, |
| 910 | + "metadata": { |
| 911 | + "collapsed": true |
| 912 | + }, |
874 | 913 | "outputs": [], |
875 | 914 | "source": [ |
876 | 915 | "spending = linspace(0, 1200, 21)\n", |
|
918 | 957 | { |
919 | 958 | "cell_type": "code", |
920 | 959 | "execution_count": 39, |
921 | | - "metadata": {}, |
| 960 | + "metadata": { |
| 961 | + "collapsed": true |
| 962 | + }, |
922 | 963 | "outputs": [], |
923 | 964 | "source": [ |
924 | 965 | "percent_reduction = compute_factor(spending) * 100\n", |
|
982 | 1023 | { |
983 | 1024 | "cell_type": "code", |
984 | 1025 | "execution_count": 41, |
985 | | - "metadata": {}, |
| 1026 | + "metadata": { |
| 1027 | + "collapsed": true |
| 1028 | + }, |
986 | 1029 | "outputs": [], |
987 | 1030 | "source": [ |
988 | 1031 | "tc = 3 # time between contacts in days \n", |
|
1004 | 1047 | { |
1005 | 1048 | "cell_type": "code", |
1006 | 1049 | "execution_count": 42, |
1007 | | - "metadata": {}, |
| 1050 | + "metadata": { |
| 1051 | + "collapsed": true |
| 1052 | + }, |
1008 | 1053 | "outputs": [], |
1009 | 1054 | "source": [ |
1010 | 1055 | "spending_array = linspace(0, 1200, 13)\n", |
|
1076 | 1121 | { |
1077 | 1122 | "cell_type": "code", |
1078 | 1123 | "execution_count": 45, |
1079 | | - "metadata": {}, |
| 1124 | + "metadata": { |
| 1125 | + "collapsed": true |
| 1126 | + }, |
1080 | 1127 | "outputs": [], |
1081 | 1128 | "source": [ |
1082 | 1129 | "plot(infected_sweep)\n", |
|
1113 | 1160 | { |
1114 | 1161 | "cell_type": "code", |
1115 | 1162 | "execution_count": 46, |
1116 | | - "metadata": {}, |
| 1163 | + "metadata": { |
| 1164 | + "collapsed": true |
| 1165 | + }, |
1117 | 1166 | "outputs": [], |
1118 | 1167 | "source": [ |
1119 | 1168 | "num_students = 90\n", |
|
1136 | 1185 | { |
1137 | 1186 | "cell_type": "code", |
1138 | 1187 | "execution_count": 47, |
1139 | | - "metadata": {}, |
| 1188 | + "metadata": { |
| 1189 | + "collapsed": true |
| 1190 | + }, |
1140 | 1191 | "outputs": [], |
1141 | 1192 | "source": [ |
1142 | 1193 | "for doses in dose_array:\n", |
|
1216 | 1267 | { |
1217 | 1268 | "cell_type": "code", |
1218 | 1269 | "execution_count": 50, |
1219 | | - "metadata": {}, |
| 1270 | + "metadata": { |
| 1271 | + "collapsed": true |
| 1272 | + }, |
1220 | 1273 | "outputs": [], |
1221 | 1274 | "source": [ |
1222 | 1275 | "plot(infected_sweep)\n", |
|
0 commit comments