From 9c2f92b91cbc41f9b3fa7eecdce40cde5f3664f3 Mon Sep 17 00:00:00 2001
From: noah <noah@theschricks.com>
Date: Mon, 1 May 2023 22:18:47 -0500
Subject: [PATCH] Infection count edge weighting

---
 Analysis/SEIRDS.pdf    | Bin 0 -> 11602 bytes
 Analysis/SERIDS.pdf    | Bin 11569 -> 0 bytes
 Analysis/ode.R         |  29 ++++--
 Analysis/prep_model.py | 230 +++++++++++++++++++++++------------------
 4 files changed, 154 insertions(+), 105 deletions(-)
 create mode 100644 Analysis/SEIRDS.pdf
 delete mode 100644 Analysis/SERIDS.pdf

diff --git a/Analysis/SEIRDS.pdf b/Analysis/SEIRDS.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..07d8447dd446eed5fd6f8b4a8dd08b79c16f4f58
GIT binary patch
literal 11602
zcmZ{K2UOF~wl#=|f`CZxQ3R0|dJ{xCf^?*}0HFo~B%y~Qy$VWEklsZ)NN<Al(53ex
zz4!it`oH(S`@WC0Sc{y@%sI39WoGZ}!>A%9!_LLQi_7RU@3Y|3>ochDh|2}w1VBx!
zaYaOMfs&435ZnrCCk29o0nAdu+??DzoP3;I+}vC|e0<Mvfr?PG|L=PmR&WS<mKs3T
z5oB*^WeNi*BjEN3IC`3jBh(CG3U<7CmzRTwgA*P5oeO#sdWQSr1uodm?4}Cte^mnj
zYPi~i0YFI*90Y+{0DvkW3os18gC3v;`1b^X|1;rc+|6i5upRvGaG)gA4vu~a0|@?}
z4palfpa@4(^m@GiwcgEuf8clxwzGg+0tEQ^xdA{KD+nCD1E3594UQDp6lw;(sT&4I
zV`hu%^0_Wip>*?`ga7if)<!CCJAwEEkJpDmT?LMG)sEoVJpg0Zvl<V^wh~q^v8lRY
zYyPde;0QBWEThD&U$^2%pWOAhy>fK&BVF`r=cF*s>tZW38M$||RugxmAG@S_eYuqf
zHH>ZBqv$gKcD27;NZ;UoDxiMlb+zqYB5^V(K!3O4^dNc4UF7O~X{4ce(d&GGuA%7q
zaxRlS=?Z8xFQy=rE>axV_Ox@hzinTbzD~68XjH#i?AxY?VS+$ogIIl)x_;6~^^3KJ
zr1|SiuVS)`Bz5hi^qK6US{N?#MRn4qOD%Hj+hr9z-pj+a*}kT=w9xtU5wg*(eEHi9
zPG3FY%JS}H+&RX$r#Gz|yGUI{**BPL62o`fndFIfL6l=-^;p$65pQ#$l!{#%ce<CQ
zboT7M-==XIzY4w9!tAjuvtWICyOej#Nn|9SmvUpkOwa2m@An^<dq!J!+@r^s36~bY
z1}a-IV@U7ii$~gr9>XU-rHY_oeuK0973YNSUvX%bVWQhu#>4!y>-&>i{Qf!wqr=U;
zS~z`dIHUN+2o5LO9}Vk1hKr)Wg57jf#p(ixkKg7P;S~WQos4MS`r5YAix9c!j@$IP
zleQZ}Z;kH(K5+RKD@>(oO(Wmt@{JXdV8qG77bX~uc-#p|$*b}_+`FeH=sZ*nKnfRR
zCmxJ0JYEkK7#B=Nwg-i9KnaDNZ2CD8!<c^(#+V-4gbH_A*|3nxM284dg8BP?>PKC@
zNXk2uKsFoLKi{CTC9Da~zqUjkcUe~y%)TcQGMH}bcNj&v$UkL}Ka^Z5%3GH+aHnHX
zz-N#LCkAm-uUHqzU5_TCLdG{*bTfc$bBDBnMw=+%xq6X}ftkSjh3y=3;Ftx-ZkGgl
zg>cE$(WHsk>MhqVL!9%kb)KemGr`ph<YwHoY^uLiSjet-6`v7$Bze9ouRB+UCMJqn
zSh`3LW4&Wdf0Uvm6HOKP9TbN`$O)v`$Ko>dkRZM{K~fSNrO%Eoiq)+G4I*SSlmi7u
zjx0kFx3UCBmaUiTM!h(64BwPaThDBt9EtSk;o(UySi1Erf6Iu?N0or)EQj#l?cHCE
z=Iu!@KOD=?GfLx|)+D6LJ6zW8mU#XR1!d4x)tgARX;=uHiOS6+1uv}Zl;J7j+t&2e
zGAEElL{#(8xn_R*g@Wy#PnYQBH_Ld%CsjU4XXQjujhG!^OKJ{T&_>c`RURrH(OZF~
z#~jK;LdjWdYicE8@bSRuc0WFmm-9iBLPh1{y6CDVN*9il-XVBA<HmvoN{^PU6wMfL
zfwY=K75Py<^ZlvWv6|80>8m=(?%44ok&&8>B75F-8GTpr=SwJ9adrr09s=q_*^r;K
zKXwsYdZ)t1^4(Q~t5#im$k_#^)8y`8;mOcINQncl=<Wl#O|pvMP69mHQoGi_MQGV9
z=~y@*KI4fFCu>`bi*(?nE7(x4Uu+s?k{;69M7$84zI>9A>3_38;5(d5(c{V&I{aS+
z;}JXsuTvBEepP+uH{ZkTo3}Xs=0+u)QOQs<Wqs5Kv6##+OM4L($>aIr-Qn<n?g1p$
zBA6(&OyRzc20t53cr~f4;lOaM1IIMcC}LzYYq2gliO=r!`umOe&AK*4o2j3IU?i!m
zGr76P$z`OjZs-=$bt9|7Q)&FC@EAi^fPcDb#eBQMc)URG5BFC#;UFz+I*wJVY2t#-
zu+soS?fwq!LM7I<(ufwirI5lrBmIHS8xTY$mIm3aY<caI{%nFZqBiGd6Qenn+6rS!
z0-SQMdzY{heg>%U-(7R!nka&AF2g8G(q619Z^)=lx+N87YHCK)%67p{QlxcjyljWk
zN<gWQp?typw3|nwxPHFGW{>8|45=qtL+Gzb?aA>`ZhRk6p46s?t2wkR8G0CDG!4p%
z9qxmPh)!kG=aw&TAmV7rBFv}4c!1V6(;{)ECsk0V(@hGmOpC*gRclH9^Cf|rvN}G<
zg|QZPeopTlB7P#3$vXl3g&%&`(mYtW(v-V-^rS-dQCF6KN%~fptYL1x>%RWR%VgyV
z43lL(rrw>oiM67}Uja_8XZhx8=dY|>kyBgB>o(l=twW^aWoZx^MP=DI>CH9RSGMP0
z1<8NXCUvux!KMY{?3~l_Ej<dx?87sc!({g7f>3Jt$2)h*Q#*^>LpC6vL89udgvm9T
zf1cA|A6M%NOZ%Ar+!1+#5zR2a(}Hbf>p-*J(w>spgh}HXZM34Dr8IKbXDp0m;p3b-
z>$(PEQWAsRZR<qIjej1o+LYh_9s6+(oNQeW&8+3UWdkvj$x@w}<(!p+)rPNVIF9*s
zI$XkvGGY*GI-io7lbr{~S2lG^A=(bS>*l{dDB-ZN*h)eE?${Yk<<~ihv$UQZabkt@
z15;pEq~~dbXPdxwqnZM1XKVjK%Yk4y|3>vzN%Qyquhb4p*qhPlCyYZTFycJThE0@R
zxQKMNBP-Ae|K40<$eHHFaLB@!$k^dDo7Z&D_bXLXXt_qU#?T+pHBB1VS9<$Luwnw#
z;&eVOCL}mxW?%UA!)K1XC8UmK7R4NXg6`dVFI|p4r(<1nlIaIgUC=H~b>D5^<UK*E
z(p}a25Sd0{@i8{YhkgaOX&OgWxHD7~bDgA>>>ewI7-M0JH#el$QWK9g@OS&S>IuU_
z8=x{R6pr3mi&~gb<v+P`=}TTiR4_G6>7G;3;CRj0-{oH+<h`h2B_2mDn+>tO(6n&f
zJr4J$kR<4fzF;9P8E-X=jHI;~;dk+lv6aP|n%z#oX}YuH6W|5y^m956yu+7ggkpn|
z%9t4hJqqj_x?s68zM$i-4(Mv3<Hg{1zb`0mvvVBoO@ViWfqxybXLnWFMjT5;J(tV^
zysDogg#3CCA8*7cra9q8x?}oDFdEZ+EM4zuypL+myEIHhhOfh7S#-oof9Ay7>7pS;
zZzV*4Z-?JkCgCPJ`oW`4H#UakM$8Htb_#mS`(d}-Nbn;VXv3JxRp{vDwg*lr+|_eA
zapzS@##FQ%6Hg+3Rx7&%3i%X`(o``{`lyZ)q`wpxYI>QAy|J-d?=AexPYp!!C7pR7
zy<k!J5pD9f!k86pDy>_z$vaPbeN^ou(!^~te2Fs2Vi(wq<(5BO7of*SHZzx`7W@)c
zrO89+v9uI3i4r)FDD1rMzLoxfMjES64u?UW{iN0UNK36x=dIRFK{bC1YQ(3?x2Nvv
z2!7JnQG8k&YKAT2!yj_zLu!P~prvfY4;;23ow12bbQ1k$=29PYl42^>$V&YPx=GTq
zmWKxdjNZ>UakEu<Im6(bnsbgqv6W>h`yEF&Wl5ypl(i^4M+vjl%H;8s-njK7Hl+b^
zTg*qbFDy;mFvGX_%c~{JD<_R=m+-G@&CDt?{r=7-EAj=zQ4U&#FIxw-Tckvq-n8GA
zN#zQ=^-V<<y$|Nr^Gt(kQDzJUe^tS-wDvpI@$0uiwfoOme2aM@K@#^(?g=a1X7o|T
zxSu}DI@BMehKDR>D%Z%}C*Z6i2zzTtBF#S(qyKcvB-<H6T62;=EcJp-JXsAZtdY4}
zy5Evx)mSPQyOKt{Fh>&Fj^4(2+MjKtk~LJyFtR&^n38T3eLhlTyDihi7WOvncU587
z>3D0L)y8OOjaA`0*>?>(5NWJ|4V?w{Ka_VL&Ip>wDP3TuN0LfoeSZ1Wl!zkR^qqqc
z;z<86SWMnu^&>7#_nqozBlmK&J!oRGU9`&YS-FVSg}O8{b4vFcGY<IssNucZ4ZMCl
zMG)_eE-EZ-rPIgu1}VwiWb6F9KcE*1B+Xt0OzoKIJH*mh?^I+LhEl{Mnt_&U5S8JF
zU5(6(eyZkJ1}MaLCNg@hOs1>uX>5*K!WTdIj7iMWa_(@8b8AI3x4v-=UMedX7z0!?
z2K%TQ-%Cf=6~8Ve&-P*%SAOH;{=>W-g~EH_M&@X#{&1dxz9uxjk9R-ItbH@aVdD$?
zE=%#7&OH^pRc|obPN0S<GRLsN<c68<O(eaWpq`EAt&Sx$PxQz@<gteJ$dERII|P!$
zTJ<)Ntr)&J`q>K-=h9v?-6=`)62|98li-X@AxQ`)*x2nXm*KQyO0)DPT($ff59#Q?
za0N`w_v*z^wJOSx{LaZrv+InRPbxfLhzpUxHS4#9m^X9n+@}%%q+*fy2~*a%-4^2a
z5`HxTO_N+`eG>AUHO54@D5R#RJltmKa~>sOF~qa##|ux0?ZbQ`fhMbmcbJe93$7Sb
z{DKw-YM{w(soxlun|!F3T6TmHdo_(>D_zw6VseyNka-ICM#GMkIepM#WJNIdp?>2<
zHu{0E^G;(P#G8t?0W;?jvaG97dBqm$$}hb0uz6riNF^Yv%4>WN!Cf7hLX)@r#9-e@
zF_Zr?>kAJ}=M5z#%{e*4EFg$94JE&5q^^fI8D05$_OYliE;K2sDcIBz-4}Ja<X|^!
z$lH*nl^mv}%Oo9mw(9yU95`K@mC0*+Os}9|1zEHs<YmwFd2?5xwb(^QS6A{~?E$1X
z?3FIYtr+A!)=$)Lo*<eU>r_`0on?S{rdYG4Mx~2g!Ay|$IRC42+(20#N2Zm9fxKki
z%eXL}!s=+IE{@4e+Tf%?;srIeNGFIkHg%+*CJ0DLSo>aAH#l{WSVc$ozVe%+Bu-s4
zDZ9BblK^2K-huN-TIbq_{iMc&q~ASZD#H!xILUp(m6}=3nE1isFiP4wL7BZ|BMY3R
zBe2jNIv#&4Wp=?kp&Jz@AP<O|Wtg^XQoEJh$D+ZsfN6Nc0P#LrUGDsd893uDvluJF
zO;FYUB=LMnbdIutw<3LyvcaQQu_z<I-}XahA4?F^f;b3$r6v=s28%<ipU>iTZ3Sm&
zI0NmO%K5f~f0mJn?M$rQ>v`okpzE#lmQ>e&#7rdewY3#KyM_YcFz8iij226cO{(tj
zU&MS8%n1iz`HNX&SVO;TT9{F-J%)75c=@-WBsj@Yn!^SD4yzr99{4E@Q0K7(ZM~a$
zZ=E{b(*#E{Nqyia>ejt6ytYlXG%fy`m_i}GKjrY&t9uh77<L{;Yb}0%TX>$KTA0=z
zRsDX)&v8S%kusP@Pvouln*yR<mY3ud*h|b*MKZz^qFP6N{Y7PWsz<?TYUGDOUp}as
zNup062_cIEEbOJv{iF}_PAr*K(Z`=71@_m}h}KlZtmB4*V_zaxmr>dfPoUUW^?ywv
zmx?tFjiN0IAw@93eN<H<(lBjOKJ+6%UHwJg=p4MMYH|8sT7{caAyj0W(fUuLO*e}x
zLL7TidRkbPsfgYGROYsWZ@$uq;1JbX6ub<4cuUyi93QhkH85m{C|&Md_4{@7j&K;H
zVi$I<<e8Vl_-LerEI4)nogiLgd<3tPxS3FqM3gX_idn9~0>{tbKetaRjh91_MzkcR
zYR{J9(`aHDX`!%r12TBFA%s?!xpWVe2xQD4SbwTxjRuGFw=h*=Zp(fmD=l<IWlgNv
zWP>g7Z`>5~o-GB9;G$t-qk7aSR35FyjPvJ^ln|2T_!+~XIq@)$ruacz-%>lGT8cXa
zmok7rhGT%VGk7jXrAIZkXOcvdIZ(BadQHB7bZH~*+Fj>dBO2z;+F6CBYx3^L{kKy>
zL=E{JXfBYZG_7w3ReT!+#I44CQgU=IPh0miiL8eGf0H8tTfY(`TWSPD^kw7G>|9pH
z@I_P$y8@@3n<P$!u=wv|rB#-`>)F)9LCIL$g)D&(Zc()A-IvKxnCM&CK{)HK&J%x<
zRxq%fJP07LL%dZwPvSrkE()INvhqvN7H{}drxN%$+xiJtMlucfqHl6^#&)$;MMyJA
zos|+o8v3t^ttv{62QFXKDq3q;VyoA{-fL3PV}5p>i1E9Kf^ii<8(BhW(Zo_1B~=Am
zV%1DIbv%W^A<HNTJ6wN-Np`<nVBklYUgj?rfuEEFU}E&OX2me{6H(L2RU~>ZA!<Wu
z@_kVs%K_8E=*koJ{3PJTk`9+sW-MDEct($1!x@u5TIrSYLDs2v28ueB%p$jn=*tf?
z1?tpzW{91iiAco@+YU!<r&NLBa;q8viOphKb?Ve3{ejeQl`tN4v{u>fOHjZ|EC11+
zu-J~0tJgp8M9CN7)~SA`PRt?AtJ>>c65`r@a#~f9L7)>e&D`@Etx9?iT6+=+w+?16
zUa5HKjVAeFYi0_B`1Z&t*o-VZ7i~^%&fUCH=`N{Bn$;1Q@5CE_w*6`P#3`A-KY^9(
zAUI`In#vOU&!pt)c}WGE4n>25&&5PWx;6gDL9(uZH>=^lh2imQXpi;_-PLgXF{Gya
zTgO%t2clRW@(-bn%ml)c`Cpl91`>$hy|gA2!9uGIc_@XJVV}}Xc_CAU!PA2K6CYAW
zBg?1$qdpW*Ts9_5t)|6S5e#(w^(DdKm$WU)JxkzdXZ2hPx^WM_D3kMR9UI0#k}r=w
zmlC7Un#f=k6|pHaL;&dpjTNz=v2j(^+fQ0WS?cb_C4GzzO&QIcu={>1Z-Bb>l|@d>
z&$2C`sDeGKwQb1aWC*zV%Ylr;FSH0bmV%nt^D=pqf<930q5|+OTLdQR#KTVl9OkL5
z!Bf;jE4vN@pLA#XA{9I=N&P-JEHa`?OB0;lg|&Q4BFp>SsG+<4g5`;0>Ki$wR^SFr
zCR(&`#fx0s3*sr~7ufxw#73<-Cq(&Ml|RiG9A1$^e=FjocdJNH^f2+`EU*lBSe0f`
zMps9npf)udl+9F*@6Lk;?v^umL~TaCtFlM-9C&eJO~<##=i7<>_mLuu);^Z7rwgo1
z9I1+|Gagwu&na>b$v?=(;^#%uxp5cv^|&e*wc)r*`LB-6=h3_b_9(~_e@cDht=#In
znh~pctNDNziW)KzNbNK{h<Lhg{s+4Dclh_eSrNSp$B)o&S?J~^-$gSn*?!~aL;Vz0
zejyR~z@gMogE1Ba$|(NtH)cwUIPdZ$JUK3QQaH81ydx;fm-Zy#yCeus73)W0%RVu8
zQy1!uqjA$C>196A!TCmlFGfoptzs@dH$}Vx^+Zu0J5l*xFjGNfg()r#x3OgRGu;N*
z>nUV2IgZ`vhyI?JQOU6WHCZ7O*+<X1DXDQgER9L_(k(B34K!(&ovAXf?9F`%>ZerD
zZvUW^?@3~7LQ)i1{hqbG7d#RTY{P3cRs0TItX;@owD@}RUVWSk_Dd#xbg5BMt_u4t
zb@YSgRz}O6<M@PR<F*fBq0Sr)os~|@57&PU(6&Tk@uF3tr_AdnyB!`A?bdJ^E{`sh
z`sWyX-+|UmS&QrdufnP@Y$p#OE>*95Lh;Tq%KfkC+u-W8^gj&iVKfjWJ7>fDbpNVe
zoDdR8V|jA?#WT*yqtlPdmSb;^?6C4H9@ip<Vn^W~?RJk$({l!XKc6yfO<RV5rWt_=
zAAA=q{A)wLt@NpjTzSP>FALYuq}H1$H;PYMn1Q^kFM{uxH!H*Rh<9z6$0sbd=Fv{p
z3@hKk_BU*0dMfmmWW9b%?Bmic)!)qtQgct=gz#u0JA#h_`1qPs+44Q_l{PZpmhH#j
z9cZd-AQzVnEyo$cexjLs-1b=8Lh$43cI!0a%3ak<wW@H&e!^5@qdEKFuZ_%oGX1^+
zR$>Dy)mKZF#Q^W>ABAc$)hlEwiD8d#t`g!rhmkj`;mU2!{SOs5$tq|zs#<k~@Y5SP
zhx%87;VVJ!UdjccnAi{<2-2^A@A~gNzvi6~Lf<IHZbgQjtWlrd3wvu$@{_zn0w;a>
zenLcm&ghr0LiChXzIF6*{!$B;7<6Cw;+EM7h?kQ{q0eIMB@-uZIGT-y!a1vxK021o
zN)54J1-;1ykr4b<a^fHR4_l&5@t)kiFENBA&+Rye{U{`j3yD^Xc2u$HAX*am6E3{%
zd6=#s4~xx>By}KnP9&`f)3%V`Ya%Dh-Wr1Ga-~~_9cTFVA3lg-c%k%Y>Y-l0(p1}c
zj?<c9#6L<4W~KFjR&QeqbAzIo`~#cWViGd#p)^k6Wgoub`8cVR?1n@nhts2+Ax3BH
zVkMa<d<5f2a6#DUzV2)%#wq3u#V+#p{JW@RWdJRRWY$lZs;X3zzyrfq7(gHE=C$FE
z;C8cIf+xTdy+f)h`Hdn_Ue6!n2VyF=V9EGA3ws-*BAagWm7|Sn(0xs}drQZ9P(u7?
z$DiW3NSrrGQ%@QUNSd!vo3fm^Vye75L1@UFAYH|+A&hJDrl~CDR)gQx%>HQpqTqL>
zaq@~6?UJ+K@oHOc<kolQpt}g(>S^}&;?^3{H50K=dDIlgKdkLB<JKB!w6TYRCg|Q-
z2|>zz3jg=aR1a<oQyO&LUD8I1JTchFSIa#<Vjmvc5$80!xnb``RN9OWEm|#RW%<O+
zr;u~vQm7ux%8JKL4*99!ky4(#Bte|KAW^VAjx$qvn%*9O=3gz6t##3U1O=kOpH>=8
zh$2HBfQ*<UD!f5B2;QyCCvf9ecM%vwY=^BRUy1*)qWMS6bvuzwqS6##d!Sk*9)*7(
z^A}8&F6j4^q3@D17!IqtkTR28#Ra^7_}cCR=gpWv7b*kDcVe2y=&}e^;e%m$+5LQx
zQKp9iH3Q`*cYq4Ex?Ep`iT(&jy7_|Ay+|Je6=UM0*+-6OIs4#7IO3$+L@>6k@jO^H
zsn1qR$kO?X*3ZFfJu%X&j<`eh$>BPYwT9jD9j}tz<%O;p(SnOT=L_O2Q6uEc8wm1r
zZ*zem)Z^l0j}VybMak&Zz}q?(K7JTzb@|F#?0MM);h>XynXb{T&8y}zJtO2A*5<R#
zb9&?JA+eQO8O-B`Z7<MeUGE2HC&u0U&9kM=lXke&zWT&qT*bDFr<MEJ&Jk@dp{?@F
zRQ<wf;^u1U1l{%AOzL$M%KKAqXJOe$p?h4fE2{`{s>TDQx4XJ9y<|E~SvYgVF|()X
zaaQJh8eK6z7*}V*a~XadM;PZJxZPY^Ic<o{)Fq!^FH#>f%nN(u1?<$CiFv`cSErJ@
zQTV3Sye^;F2r1wmxRSo15kWZhZf0hmo$bxFQp0N-7^-N5j-Dr&(Ayp5b$U5bD5h~r
zn3EHIMVL8225ZRyR*xBZ9Td_}Wc=Fgukbosr@vkRHe4Siyz|;W*&{C$TGzUmk!)G+
zE70bJo%!<5#&7SBpL*5rI744uYmWywShc|xExfLd3cb!|JUh;=X^Co_w!Ij<Vt*NF
zIXoazXqcZ7y43B)sjoH?E*wF6HaHm`-ES*6Jlq)9j!9xqTD(HuHmVSt-$#t9E3>v@
z)Yt9SE^v)RYHHSFC@Y5&j=!QV%q-pTz)V6g{8o?!l3`l>=RxJ`h1N%e$m9av^S~Dx
z$K&QRvIH_q!M#OmL(a7}CW9NtdCdHY^Ol;{;BcDvoy1Td^$AnnoDz2evNwOp8-4%A
z6=~Nd0(F&YyOTL4Rz`nSj@mNv7s)n1cI|6dr!QVqR>aEVulc=7;^UXAO}|{N1!}!b
z$BL{=|LEB(B9tAh|0MD#F9B&-a%EWBYcparvBAGJ8(($Mq6SKIdsh)>O{HBvMk-sX
z#+qlf^ijm7wk>+tf}0&p2on;~L=8=$-dP@jC-JY4$a9=&@{#&WEMNqBmHMyX#Q3`K
zk0rZ96A{x?|E>iQHFfZH1(j<{U|71rRGDxC&G7hO-o|nEJF!{StET>sP--l#=Uln&
z`^e_v2KDfI-g?>a<ZCIuvtvw|p_BdOoL0W}(2<iL${MDj&W45y7KOB$*B`|}UKT2&
zTu*~OWp7(b;X+?rpOPonxnAr|%y?e5wH4N%UkY5Gp2YROm@yQX^*W``2>Y3hHmm>o
zzJG1&8&~>|Yy6LG{l~2aDp;Ap0D6EMms$;AfOfN?2(&BB1pq37&8$EYP#1t6C;B};
zPM(`zga0tLZ;a`GjA*ofZmJH31N4As2OpgQ>;gBy1<IgJZGeafP=*`8b+gECTO26u
z0+&^X-x&BJA~(ySQ$nDQ>h>Vh-w_~8O;!TH_s{Tu8|gQpn|yk}SCWzvAQ;#Tz<)E!
z;2&%L=G8xY6u<Gp?QW6)rNB;Brr?`Jfl6q*9r)KY=jIdQ|9_h1H#PmsG*66m)=(cK
z3+$~Es8f0WASAla#!~H;9Dbc#6bB2}y$AA*kr<Epr8MlNRHL%)ax>3*r%Of^KMwK7
zR4rb=O(`|xU{69F`^qeza;DaG8<j9af9z?KT-Q5gv5g$A#)!P|#s2<U`=0mxRBe*E
z<sU)ep%%tCnD=)uupk(-_NAo~WG6DWKHc&d4TF3lJ|VBGCIFnCvkZ!bHDO#QM_1`Q
z=E6|<7`D83GNw;S-XSQe^_j6kx#QzEUd0|x`JTr*Z@D!;!1mbmJdQqr35E8#{X{6c
zU_JBToyEzg3*kZ!>%*dqDan`bo>XzyGnx9rUj@8ug)SCTdMNQ@Fu$MV_ytopgRq;v
zXXbv(j{d@JR<6TnX8@?(-Amj>_@k#=m*e4+>>x@eef5Cw55%E!y}gRdp+CAg2?W9!
z=D9{cQR>jlIbxXiUg(gx+ESK0em*KFB*<yYZd=bh<YjP?@c!5-*3i#ui&BdIr1?n8
z^AkQsQ5Lv->@3jMiWI|O9aE*|G7x19#0V*=M-U(#=&C&o)mrA+ELig<C-oIDlWm>>
zSGLcR3A`kYJbKmYr?!_sF8dz@-@~~i4F&AX-Lj-w+I!~hPw)Rq4m_d9xkAv95nVi~
z4+2I@OLb+v2E2|f*J~#Z!u})?@bV-%{)DjkDr2lFHhBDms4rIgR^}tBm!5PO%q&EJ
zd+uheecMR9R#~i5RA~RFkBR5KZ?BMAm9MYb$)0#cn953VS3Xqj*@(Z;y~PW_kba4^
z_L6!9i@5Rjw??8$48mEgTjI<?w;hbRX5W%TVoCwXNxaKbD5Bry8$ZptgGeFb_2svx
zfcu!v@(SW`_;eZHp2Sq+BuMgpwu9UK@`0a!`vX$7KsZ3Mou<TJJCb@43tN&%g(mtg
zqojn~Ez+0kDQY=1Ujo)r%~$kRnBi2_82FMEDXNwa&q&>H7&z}`HR%=+L9k1Wl?Crc
z_>nht>~IUxbG)r+jIqb^c(CQC+#EMcq>U@)lTGCPfJitkl#__;yNpDX%aeOVkMbfP
zX56xRIUY%uL0=Uj7nR+vM{-A6noH$%G+rk?ssDUhY^r7|fBF|y9F@eJ7ejA`n5}US
zi4<jQ(op?9D;RJ#Vami%4jHYKgg)Ap(iNfu!86i3zNw~RToDvruW6V_OK+7@mQoK!
z2>DaChJmt*wZ94v)2xT!|8OiafpC4pb${Ug&>j1E1k#sb9&Z0l&`4iIDTuY@&8)8d
zQ+q9Y;!4Wehtm(X!w7%u?{w{gPh>m@J@1M|BtK#%4!c+2Cnv!IVo!bFCnKm(so*01
z<l$kEQj2|>cV&R~D-kIXWg5nMMpZfUC;ZQRo;<RBwD~mg@zmowz=sYQGeSoZo!8S6
zc1+2SW}Y^vrm$$U+-BY((PSNBhO_9ax?}}u@oF4s=zslKsGaH}<si44`_QmME;zqM
zRX0;39j>t&Fi^(*W=uv%LrBuK)G>4qdr6jr?MKkO<eX{<^EtBxv#EyP7iF+|)evGx
zW=&^J5MPof3%gIvDzgk~>09=UG$bTZzDzzzKD^U`^cig`bwAe1>H+r|=nQf8{0;jX
z<u~<fOh5U5vVJN2`76gYXZNT4(D+Z8q3K+nA>`2DP(`kD?#$pztB+QGHjCC9{epSB
za01(iE}E^+?}Y-FCx`v`P`nJhT9#XbG<~0k2Zr|sI|sfD2&W%vBou6>CVr`tJoh`w
zAb8WJD4S9$xt!5s8i3pUA*v_G?0_o^&WE@}FjgN{`&RR~wYsUg`MYhp9Ut5~W4nBE
zNqbgzrhZ5G&Z9fnck=I4;O_>oHIp<O1Y8CbhKGflgz-}#C`Gw7c?dXr1nl_gZ6x>|
z1a!;?j2%%Vrk&i1JX5xAWB2(>1r)gQxTv`UQMBFIziLo{?0sZPWD=owLaARr9Q^R9
z83K;b9borm$78QyOV**({-I}Ex>j*!X$wuXo3cc<1(*<9s#^}4?&sxyVjua|dNi)T
zr$1Abx9rr5w6m~ltC%o#9!#>PH)&|9AGz}NP4{h2MS%wT8@c+6hj;IUw&I%FnkPC-
zJa-k!cM|?}Fg(_zv*9s%N=}}o5U&v0mG0WQf3i=t9A33n<*2`IkZE_nLaL(9PR-fC
zIl|rd9E%W1m?Zn5%`ZFG;>)sluZ7{T;p0Mr!V!qYk^Ip!Xv;)zO^bV9;5_k?`X+2f
zt5;yG;-_<G?=;kZ=9yQSSL8MSH8rL#=A`#8@7$MdjWIq*@r)E~@t@*^Z&7chn?xVk
z%Ne&HMb_-&FXJsc`@ik;?%PZ~{vOsG>yM<5xJQUrc+ZOhM6CLtzrnlOxZ`ascH!f?
z$6b%Nr9Mf$N^?smO(#xENH=P8X^U=?GtD>AYBdZQ3?gnpw3#lXFC@%olb^9oR@sX*
zy+{?DJA1q!?HbC#Hmv=*{^aE055~`ovJd-%xq{h)w}nza-|0CR@Xjn($gYHJacmJM
z;XrYq@3%XF%gRrxm|~ua$~!4eCy9KTPg0B8(`Y2qz)Q+L%O1|GSHe@ANIrTd%2X$c
z6&Ll<K|_ux>)ttcInO?KU}<-0NNGmb=eSLMen#Afd|=73(2|-TykosX-YrrkhTQQB
zvFac8nFNH7<&qqeIHAT>V|$8SPxp}f4Te6tO$PmXDSAcUjNzHfoo8?QzGUXPX{{G>
z=!(|2LQy{De&tp+Ud&}2Y(S4H4s!$LD9!%(wBke0b-hs15a}NR{FRWLiJ)=jQNc>9
zt*kGK^a?JF<O=x;tnxF^Uv~8;F(|DVBjot|GH_9MDx_DlTXSxBePXhLmRXCV7f3%r
zS5gPVF_rCzh<hHlQL6&dQmqIrQYcZ-rPLoY0IQUjuOSighjQUsithR5n_1%{HLf*<
z8@r=}n%Qr(bdT){jki&54QD-f(gQ+}ACW@mbo?>`M_#D2>Z(`ixn--(o^zh7K0#PW
zR0Ime9?uTWx6jK!>=}P=>Y8!!8k}l{<UzA5;^=M#J_)4rN;ug(Y!(l6e^(SyhO{^_
zu*E4)&mdmYpBDNVGr+hRu}=4b%B8HLEi-$;y=(iJeo#~sroZGecJF&|8ga7HujJ1!
z(ESsutBZrCbqa3VpMynV`M_G}4y12vf7zoqarJmZ%ys(kq8P7-Ji^e-V554`J{Q*L
znX-U$PA`0i_RjM_n!sw{bAv6)Bms~4kn@;9t-8Xb8Z!IMf%c22xpmbM=(Jqrqq>fo
z(-qCL{N)OmUGJn0j2v09m^4rqavpz7?a6ytY_J9peX%Agykp$t``mQ0_^BB9`oObx
z<<TkVhhLA>KEs4;8x!&=P|Wbe<>J-whfmo*#8|KO&$}1LoEuy&xH(Cm^Ihp5TnrW_
z809pKUJq=6x8b7)ms*sC=7w>#s$Q^J_I2Cqnvc%JJXPDp+pJ<AJy%XC4tBS)!<3Uf
zv5x)r&;JkA7Wx+r#r@kMDs6{u#<a3Sw}{@*)_=d!va%DmgIWC(zy(TKnVW;rZJ{@P
zpLzgZKA^3Y9Rdar;s-h)pl~n*Y!1ie;^hUJLD8M8FqjnpXyFKQLUUY@DdL9enp!!U
zB5cheU>5)oZUr$zhipNnj!-)Q(8LjZGY~x!WNM1;G{ps)S)o(GtY82jhaD7d1~vx(
z9e)>Tie~;02<SiJmIylwkR!qt0z$w6K&S=O4s3(Y2Kv1+?1s4;+_akh)f&IK{y#+g
zPmk<xB>+@Mn81Ix%-#gi5<>}Oi!KSR4u0eEXN(*e;sl0UnSyYE|KI|+LG(XcYF!+`
z=D3^yE-qZofBgXX`S^JF0Oo*y#kkRp!#6(wyMM*F_|ZD(f5!wlxzW1ezhYdR=#>8z
zL$?_J-xwcS1pV)EoZMXK7UBPn3H^5roq_kib8vBSbN)vya7U081nh|Jdj<g2t=z!q
u{Q`im?98EPb#{XXPz?%&18(|}e`5oKgB;<%k>TRv;pD+(WRzBx!Tlc;)#P*l

literal 0
HcmV?d00001

diff --git a/Analysis/SERIDS.pdf b/Analysis/SERIDS.pdf
deleted file mode 100644
index d8cd1e262d15ae27005846c4cfcae75f87396c1a..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 11569
zcmZ{K1z42b^EQe|DBZn)gw(P~r_|CY-OUoxuqYkUEgcflNOyNECEYEJq%;!pE%?6g
zum0a&7t8g`GiPSboOAY>`?-hqtvHwk#L9t2>$B>!<}>UwuH%3P0<ZyWjozUN2%rJQ
z93X~J3tJm;Lns8mAkN3m#{P<plMMvo;CRK(^AZgxZ)@`Zy{BRUwStdQ0!TU-+L>Dz
zI|AM~LG7HN@L_KqY)zbuArAL<Iapt@vcY8!zHIMBFVTdA&>%J@_ffF_6Ab{U;$jB@
z0L2WUhE}#_0N`6gGl(PL6<k3H@UH>f|1;oT?_S#hVgvn49VljN1BG8Y0(c&V1C=0-
zwoVSlaC;p8W$#|$?{vsPY|NnM0B$aJUH}kmVFiWH00_2%Cr2D&Y-<9!kJ}Ln&zUuv
zb4s)BKr~CO{lc(%3T^#Qd&F>0>z89_FGW2+sW@khbciToyyVDcmRVW;pb9Z{T3NDm
zl<L#6J5W(Iqaqy-m*hgI=^3WnKf3PmYQ9<fkbZX>)1$uJT(wof9C>xpSlZ&ENikro
z9WZ;f^+sRV^ZM8};O6e;&`pEtCI%>k^?Q~6Zf|MrXc^dgRULIVDC>5!@ca5!Ry*Z4
zP;WIoXF+jkW?&5$<Xkj=!bfT8rne6wYY>L<q)-1uyS(V(=I-8bIs2V2oVt4F?^RB6
zlcueaHt)LJaXiEVaM3>BauL$ghs~r9FkWK<^?!a`8rrH~A{=P}YD{mq#6GV&d<NSd
z&NBJ3PpI?W@|o?=xGA25t$K;N{cA@n;yg{(?aqUm4J3{MB+L1<`82|jB0`;k%Cqao
ziE*8$WZxrdx`Q{-+na3b_7_J#q;{}6xrS|N;*7xmAN;XuwN}6EZm8OD>sqr`!k$>&
zA3td+T=?9y`8!YgZ4=O2v(Tz?p7=vtJ2PRJr@%B5@QWsE*;NI5!q&{1ned2Q(^j{b
zd(|&NcY&fV)5(UkVdL7zV<3V%GwpPbS6RP=^8qCKzjD;t)xnWwmchP0MfNYrMbhh>
zyhuQ%{@4kba$sr9y8g4mgg38O76zRJ)yD!m%{7p$+yunHBkwHv(|74;&>EVAwr2>v
z@Po5_iafuP1>({Vp`-;?bOaSK#3+s-mN+2YdYmK_g&*m-jU*Vm9ZYvFjb18i6r9gQ
zZq@XnDC8x4+4dJkvg4q()RB?&($ii|UXHE^z)q-;6YqOk*N=Z1$?~q_x8D)QGRF1%
zexz}vOr`LxtI*A;zKOF)&)EV7@16_SQdeKlz57D$j*hb}^CZr&-hOs$TGH=~>LbFD
zTSi{|g{dcKS5UH}*Xen!2m9N@XU-{VJD*mNG}laRUy!Ru)}Zo`^pYfu*$dUOrc7RI
z#NvY&UiO>Bo<8QbN4gzo0J=GnRX+9LE!isD(JFGRSXcXQgTy*$=@(Bp8Gqm*3sWym
z=o$-|s7p>;ToyEZ<F)|a4Dc^h<Pyll9|`6-$f`cqna={|)jlbm^xdkN)gN~oIJW|4
z6fd=nhl~}`S7;5?35rjL=it}da|%jwyy{|g-b*4SC;`ug#1NLpaUbSHn$_+WuTePX
zj*YMGAAYVS^w=MjIbEcb(NfeCP9JdT0eZFI$}%~uy>o6d)56`otNMoiOg+o0yIpps
z|MP{>V0`n3#iewAFG809?w%zq?eXzq-T9g0^omRYOYxfhh>I=H7iwzju*uV^t(L_Z
zB4@}!+3LHv@H*#5{70M38ojQ!HpOl4SM}Cso9#~sI<+JYGZJXDf`v>ZT;AiJC!4rm
z*ie^OQFE8%_Ka1LTDq6;f6%B)hP+bH>rz{4A`Dx+a&q=>RvX4IP37L<vMEm}+aL9u
zVMtJm{lNTtfzhToUU;ReKV;!Pi6UL+1A3oj>n*w`8Ti-Mk@mViUpRmd8?U4i51$Eo
z9JPE&J*S>&KHXdt=igr;zMq3p#72bn$+9go^Zgvima@?-(K|HQxtVpsEAZ)+jhJkg
zLaVJo=8<#<)y{iI`%4pZg{DMqRe{{>q4b#>8%~yzESj)sqDVxhXAi5;#LzWYMw866
z9)N$*QL!SO>yzrTrq@s%dd^m&9fqk_NFRW?sIAE45vo?>mXlKPF;z9JnIvSMvP+ND
z5F&9<r<G<7!6ikE_}2WDe4v7u9(W(aK*N=9%N!{uo0Qz}W|n_T1X0eK^#PZ48Q4!P
z)@;%QT^d%rOTr!huo7MKDKndqN}CwZ)GMuEeW4F4Y3)CZT$H}#0z91w8_p9JnkjE6
z2T9InmRA!q6f+D5YRr=hyI)3YX!M2`{@~3Tcrg_~z2Cb~X2;pDmLeBi+~fk3ggSIm
znmjicd$o#QMp}#E!TeeH+tQN=F}2vhc0Cy8-N}nT0dSaon(3A_NaAPK8T@+&q1?_~
zq&A&}->Y4@O6hDLd+eBDo>Y86U4U<drJMIInHCfsxttS~xn^fxIP<sZ-xen;O=rU%
z_VEG4I^nJl9LeZhoR&{#t7qNBZRK6AJ9;{&@ky<0w@_Gzm{KT8mbo}%RNGawPs-|q
z^_3%cxK^cASNF4-VD7_iQ}L-)OiGa>y>(4v2itM_-S&XifThzpE1M|%jqm|}rf%yx
zd!67o{Q9+bE@A6ZMd4*rY4wIw-0;aa1>H~nJXVXFZ|V}URc$!u>AO0*sxOcD$YhxM
z`2hxOHkFK}EJlA8j~pGZY&U=o_1^~D1C4`4+I|{RT-}~Rn{v%2vIn_}jqTcp(v>vy
zmdeU6YL<>Dh6TzilGA!mXIRd<=2EYjh8E!GD1Pv-b2NjvRgf9yX=Ig@wt}A9ijpyD
zy_6`^;e)?%hQ-0|>5>fp9l1XOo5gxYgtlkw7nvZ5_(2{QY;Bk8<{6cd@8{!76B|=0
z%O?$qFLKNeg;JX1vOLbXajA!d(j@abe3xx!gUUzzgW3<H4yVF&TzP}+HzfG*JX?`y
zU~P{QJ4y>A)aIP<%gvW{KKa<JiA}_0VtBU3QO|t(=)H1S2o|i`H`MmGSSn%HSH2#s
zz(PbeKx<{hn*DNRE<rg5Oy3qMD1THh|0rYa`-L`{5i`u7>-$CT2PU!jzJQiC4b<}#
z8cNSr2-P2F??k!MG$o__iUv;8M}N$1!@L5Z@DK%2xSmiNT)M3{b_2FNu-2NmH#-k&
zak@0@0J%R4**`g5x~cj+Z{5Q)DXOphK|euq<0<>PEaM2FbD3tj_--b(INED(s9}Qg
z#xqnn%ty?Cv9;jFnMTK-(l%51g4Z|=Z1EJV=z)i>D9J+<_#5D`Cv^2<Nc$Mi-2kr{
z86z51gO)E?kX{(N5JB`w$%a>*66bBWkD?!igjf$Y2H1p!F>*1PNx6+2gac6{l2}#b
zm_bpmQ@DyyoT*_voLL`*urM%l4V-@EctS)ZF|WvxO}b(rK1IEEIhY!3`P=TLC6UCW
zBA!(clMcOf&r2QE(Q?UecfDAud%Oz@7OZR`DOu+EoL!?c%A;S5FU#cW%-7<a<cXsp
z)G}OcswFWY#&|F#2{N}BojAV;qtZktvsg_7JS53cEOS3!`^~QO2pXvN_aPRjwu!%D
z6{$pfKwXCuj!g;C3jQK2SKllKYE;}^Eu-n9C{7B_=);nvsUYZv*$m|{ejiG_VPguT
zo@OAi3M$m0FNQIk>g6O`hSNWO*T<=&oRRoAYnvk^;5iLhSV;*u;<8I3<t{@F;;)ED
z^2bfIS<!^CiWAK2wfy7YH*reBK8o4I5@Ifo6ySMl_=ZpU)Grhhu|IivSjCj5i_Q*r
zS~{v+iBmI|X5!^LYV{!HSUCec&H%+A7Wn)flI40xUg@s7%6&|#-Xcqau%mHx;g`gq
z80XaJaT2X@u=WAoPqW5~()uVqA(9Z&e^{NJP+>E-Rkb^`jG{szrF=Q?!~Dl|OSVHo
z_l#1&3CH#lfeMCCz@s{N=+Y8AL{Zf!tZMexNp&mBvvHhZnKcZ5+)NmiI2Ou;+V994
z1zxP~GJ^cnu(jAJiv7PJr!!+Tt_PK~`=~#UPa*`5YTzMJt)eYT8<e{m_b~NrkTcL?
zP$n{q;h+>#myU@2dh(Vy?@Eu?NJsKJJ|h37zhW%D1h>$L4-=GLb0Z|JIiFFlxTWNA
ze}HmTD^(n(mV|1bb3R=#`m?X0`3~J}p}kB`G31_SF))Tj)7ZuRe36&9yQeB|ul9ag
zg~Pq9j#Wv(9G_FdKA*OuLv)8a_wsWEszuk7YOo<5t%Rc8ixR9f6TDKs_nkjb%l(xC
zVorMym5__6vluXv#5FW$iN_r!@?^I0OHp3w&k|P}NchIs#T{pV>7r>?Md}U?XT|T7
zJevR1FH*?BKz!&7@s69EU>h}tT1kQPg%`dKwlN%|W_=f55O=|Ov(ST2_5)Qh@^Qep
z05>XGRfu!+5tby5Y~k0wjI9t&V1~M!TM{u?Q2`HKp+mKQL_S`l^ZoN_X&CL|x~}b{
zJS=TEoeY%Xn+!T9+P0S;OY3~S4G0N%O-&Z=e!M-iib`!0Y8I+l>wc_3mt`QPv!Nvs
z>vLvU+C6%UOO_oDQk%y7;x<y}kH_~V9*aKf2NNdw?#3{={Z#$DV^0yA4X=b^8C8ix
zaY(-_`4twk<a4Jlj&*ys+?d#Yu;Ky+5=>t^4vfM0^&l!7yIuS&YLgs$OQw)$bh|jq
zIQ1!(P1(|JP!q<;GOV`DKLjOQ+b9E(56zz4L5oiE0$Hq%z!kTWizVu9QqaSYNiDU{
zzWSEf+U}r?>g~6$SA8w>1Vq=Sw70$UMN`65OM}T4-6U6%g**G_PLuQDh3);3U~IZ%
z^gdzem%RXcN<IBy#_?Q5+8qCgb}ggy1;JDl`dimB@lcJz{j+ZiWg$<*BAuT?sRe=c
z<q=i3<I$vtW^Sx*N37246+{{m!F@dO6DYE~-&`UD$|=+|c9WV{Ol>qIBKsbm8pE<S
z?gox>@b87nr2#$Bx5E$+&8)@j^EDkMG3Mf4-B{vux?h<HMjqKv;FgbSSJK0y`XYPy
zvdfM)kCuD;s(P>c_sx9TfRO=yINh2Cd>Hd4E4zFn28&<VB*M{bM7;%m4%i-q^nO(^
zm~&7ovQo7gwzZi?G8R0%;t`7i@mIsPVsEjGURG)~Jjk5oaO=O-3$7+K*iRx=i?U`{
z$c0ykF%IMBs_BJe#`{~W<lU(QFO&S$6m9VWU(*y7aawAn9bifsX>a58ldFx|u`B4c
zDKya&=UI=2Hmkj<^1)MW3b`TvxG7pF#Y>D?O;=u3F%a9cH1pQen;bcta-Cw<vfxiC
z>^26y8C^$Had{e#AYynV=_W|qHp1AT08?kSM=xStXB~?W-kVR>&{rB~=Pt@T7>SWX
z<UdJ!LSbuAho_xZZD=7IhB`;$!8A#vC9$qZ#urgm^{H`NV9C2!-N+hiA@Un~p_Cx;
zeZhGsy>B5nQ8c}SDjhl(%I&koY7wBO=*~`&sbFNWp|x;YrJ*BY6`IlF(8<{F=ApcE
zi}q~7B{iG$)OzOeQB>}P@f!qG6><}E^wt?imPwd17kMYM>`0o;aSOSJ`Lq2ob<J@a
z6mSS0ZblPH%eYYasMVCH^wB_-R?2=G)-m?(mwGHcL0oPK)SsO(?ni%XrV7I!KT3Yd
z^4!xFec@&^j_@xY$3FDp0@QY~$-tQ^1rHwC;RyW5i5O8`sb-6q_ff8RYwOip@A__s
zbgQUZ(`5O^9sP`G?YF}=@Noi%bUdl-gX}}Zy2j?JvR?_(m=GCGe=>}9n=~6@Cbd6$
z4Br(N;-0tmgVfjrfl@CzR<VM`8)0%Uaj00MP&7278uP{RSaIg*9u+J;#v3qwDASfl
zLo!Q0S7j#RXRB3ptB7cg(#^i8rtZqYc*MlyHkC<Ry$dt4&UP>xPKE$u#Ko;PZqp31
zA&kccSHO1jE_OUrFipZ&DK?`N9PF~%hdazmucbi^MlNXhS!@MHe)&?5c%nStu=8Cp
z6L&ZC#fBLnImv(x<)g9J@PdQn`+^Ep$144HUJq^`kH&Fi$)*Qq6N6-wo~Vtfa?Anm
zTr?r7qub|dqrEFS+RQyWwx}0!hWj)hRNCt#Xe7#ULwQK#S98fi-ynHcHWDjk{=8-(
z27wt7R5QovSoE0$ci`W-i&cjG>UOg88TL$Is2wSR9ml6LwXi0$Gg^T*&8;m<w@>)6
zVxpHauJD>*oUWK9o;+4wOOIBJG`A+P*nxjCK1S1~#t2VaDJUzdkaRsL^>at}7$RhF
z5}grOvEosQLn@a=s;+J*nXjx`)0macQdK|z+~8=twMJXt_KwPWP`qkgq8Q_DCWXXb
zrH3l$ID?x|u)o{|^Z6NF&6d=k;9lcWE5HgiA{<UHzliCayqD4Z%%>d}%rN`JW|i+H
z_&i%T#>fP7#+3;?6&|4t`lyOe)185WDr-3{^-3`D293WgC$wJDhC*p}1)GIAIzQX2
zJKwA{jpDgk{u~se;mgkH-ek@d1BQf@hSqFh9=+>k7wlDZeWm*FAJ}pD0WWK?UM^%b
zA{Yffi{=XTC0LSCMJbFZk6N6$cz?=UN%0ut^$bO7vBO1^W89!Kqo%rB;6;;kVUPM{
z0TH%XcfM^hL8H;ciCZfkJTxL~CMMxud3oI<;xn}8bc`37qi^E>VBIusi)ukt0+)ye
z=|Wwn7=ofqA_&>%GQS+BUCraJOjwr@PUb?roZw$ObxU6BV(*0@QKmzdO6Siw(ys>i
z-BR{89n@-SLkeTQD=g)S9t;<P{EBcWKg~M2vnj{~{w?jK7IxG4juFi2qw|$$lrcmR
zcAqUK2leh=;q2v*QU5iBe4_S(kf_yl39AkGkSP^P_{LejYS8&Qrvrk0DnKk`1c!)W
zKjtbH^X)m^J5s{6tc3#EGvS$cj#f7?TdFjIu3tx|by{i_4BCri2NXDF|7MyDv}Y9U
z>hHgx6+HSz&$SfzE%i|LFMKf^5`}J1<D?X}=X0tMUJx$$ttHsM0HdgY=;gz@d~PAi
zHFL%e*3A-AIp9=}%>M&P-PFQ#N_gqIdpA;|3c_Ut>nj4)RKxHBsb9ch+@bi7;K+jD
z@X~^$#=O<9e=Vjah(xRsb}BBSWDQePCD%{|DUpSjmsNq^IUYah>v~%JZIRRV+ds-$
zwJI7LyrQ%@?<!%_v|srwJ$s^e5F_!Z495F2yj+bs9gU1l(&`;g)Bk$gP-!%AIQaV5
zKnzp$GyZ3f{GFl;Q$T?bJY6{PQ_aLuX+-g;^Tryv3uOgHT@S!44L|F^O%6T@FmgFE
zfi7-KpY)^;Je9DoDA-JvqYc6?DZNn|mP1CitPY*`%krJ}VOh8oK3&lcQxNCRRlOa{
z*{M>Q%ei7xkO<Pyq?rV|Oua4h>po`OXi3oj`cJavGVxn~1n<953e&Qh_+<3t&pJ}8
zO=UW@ekJOT4_Emdq0vXfxid1htWt=TAopai>zaIY`&7~qC!BR(|MYit_akwfa0As*
zDf<x7t}(u|eR{!+w%qSkzp5CKzRG9HqN)9X+t4h_*)^5)w?ReG&a1=xO9iAZoq7-&
z{b+lOwTYk&u^AQ4Q^7dW5?FIp5Yf%j8ExD-4pU%R@AfhHvt=0TKHUlqO;?*?y*la3
z<zj*iriO0n&g=7%5$h3Qk%{_i5eA_$Dx^o9%azk#7MpU`8Zrm&a&^?8KY&B35<xzQ
zZ2tQAK{Y8pqY=D#*~u)X>Ld0wDQ+U`hRhR`MnBz#8q9l|f)f)JLHa|r8U#k)5#(0X
zIjK&_&KA;GxRf);7grIGrMR^z4iR*b|KP=G2d_D*Ic(fR>XQxr;lIHtqP6a`L<>ww
z(ND8nHW+1%3l5ez*<>>#&Tv4hyNhDMJ&h2o>3-Vn9ttyMeg?8bk=sy@qdc*zq}unX
z7S3oUDvk<uF~FOI!z(90-_f&2ZzG=e5!N9+3;5(TI>#SZzhL<hp-`N@yeADdRWN1z
zA_aS$Ae}9b1ThOmFa(?jBn#tIXB`<W3o<tzgpp92=pHD$y?v-%SNMML&nJ_rZWEO|
zAfXLB{fhngd06H?1Ibvp9oy7Gz9=$XQdjpc!~TYKFn)X=!weQ2W=0Ry6b#1Gyhl({
zIZ|V&Rc>c+h{(*C!jd!|<Ya?#gw^>BD38oZXQA+Fr+o}axSlh;M~o#PmTM4+WXaW0
z2{nX(NTOYl;|a-%K~a6v#|VnXA3@qi+~k$4JNl}0S>dGq2_VN*Khh8KOdd5FA+LX!
ztgZNE6j{|0JzDWsG-vM`)ll}x$-1V%T`@xcSizNwp21Fx1b?GCEnRqoy7*%9l5be`
zU9f8EGY2eX1H4r5=qE!z-|}}+wPd=K&m)|%<73SY@chB|g}7bndXn{|+-nT?Y46t+
z-7vBp`uXMy<l~UcjHaJqL?jnWlzgsoH9?m4FVHHv<PlU%eVL&Z^V>b&oRm!5zk^bf
z*Eo(x$mkQzHrp89J<r0@4$3MoaEr_8@g!T-LQ(oQr&MZ`6!sC_-hcN+N1dfs^%}Wn
zS`4C(4<l+AbsGc8Rfrv%T5(2fRu2Zpsur1JSCRZ`eZR>cbiN4Dt_EsF5S70DHM7dh
zh(^N}vYBa~kKfVWo4Vrwy`!Vi*`rsVS|}=EC*Bw6cG!D0E&a}ZS~ZR!Yh_>*@SdlA
z0*=g+Bu0HZkJ{QboKWJ)$n!S$$k2N-cza#brMmC)wu&|_t$ReZS0y&kqgO9BfD*Cr
z1u2Vx$?h%BBL!G;qsjL6@04fLnKDV$8;WwwF{9y_!|#W5=ta7FDvNhsDy+q38$PL(
z2tf{`&SfCc3bykNOktg~;dCzfA%sKNW6xkYULt@0G)7RbS0D*#9AeD41UN`_+kI>l
z#~_jrHpDwdvF8XNapW}u19jxK<2W&O6_dhj_JEJ!+pQwFehykSNMOdVSjYC!jlv2s
zONuT`6q_mnCdC;TK}7$IFp%%1<@l^RW6Ylgjm42;33<E9p+rA>P90-)&CeK|nR=~e
z6aGXqRD88FNzfjIrZ$y}fAx16(iLTG?(a+XIbBB7#VFJ<=Br*CLq||%h8qg9vDAEY
zmYC?H2zmi3%W*X1xIc%0@3!^a>f}vJky&7AXgZ1KUD(>-U0;A!B+%=+?l}MW5<B;-
zrNxVN=X7AYN0Q*z<<Z8WzR)#0yO(g4b-DH9nctmdi{Plw{p7;W8?V~27(Wo6-zqPA
ztp>_nNZH-F-)UT242kP=qUf)h>|U20@?7#|x-?usXD-NxJDc+=TMnJm_65n$_Kzq%
z>uw{P<hONrNN$|?6r0GXN`68!oq0ucj()s8a`1MB?YgXP+}L&Jt#u6?t%SVE3!Z*=
z#kBnB8~7@Ec>Hs^xITwHo{1aN81nuRwrVeTE3Z)hkrOCbXF23_2-srv?AIpu5NYuV
zWADD$;rvCd!n4Yqt<wYxYu0*uX*YAaz*ZN%<8-f!tFu6~yGKZIdv|ZqtS#M@PLDIC
zuO=n$uF9K_BO_aG*Ytn?r1bi|Si9_gv)hv(>@}d<x^l-0(n;^Wa1{=T&?hlc@V?9H
zTyCQCnBM=H@at5hW_shz-RZvm@5|O(+TFWP)bA0Ng?)u9$;%nXg)dTWc*jdEgzUz=
zESi?yuX|ZH6{Rn({oL_{<?T=J5ZjNedUmIUz7n#Hq&X)V4e^}YeB^21yINT2+0oyR
z`QwL$eNb9ixebRSkS~a80IbvGo$5So*hX5o*s6V6t9l(c)WtOUJEQYxNo(KDGgx8R
z2YtHkOt6Kc@ZI?hPH}-U*`!QS(G5;&+GTRe%;i{0{fR8}uPW^~_GF9E=&*Ldt(f+X
zY_*rU?X7m<VlLxCe@q&a5ZUpRM~Z_wdXd(CNS0P`oRbXf4JOwrF8|M{uzJJEcV;?G
zb%m|Ra3f3Q>WBN2LO>a>(z&h4AGXEsNk9DDjyXN$dGWEf1x5%~{l~~8^gLKVu!R$C
zw=z8b!3rFzq&nXPz0<F0>}V~uzH8%NoU8HrzIFf)g*))<>zv2zYFzEo`>Lt@<akmd
zRi}qwsy8LAT==vcjl?#F7UbuVE8s7;3#a4t>X*}3Z@o)ye6v8TZJikUhrO+rEh{t2
z`QmP@dy!kTp64GJN0V`PshYadZ&PKecwzo4qg!5iINI*k=QzS1yJc5B&(WqoJoRVV
zd{Qp+J?s_7n&Q=-S-RI%_Zs2y@6BCd0ncBHhb_N`ds^MEW}WA~?>4*;4A8%y!MW(a
zXx+b9>3h=o4=MaNEB%Lz2FhBPI0Cc*_aw9uKo?F{+d9EXXAl4=4>7SY6t#5*XtTlZ
zak9O-|I__9-uj+*{+nlp)8od<5GX(!2q);_9w5$8T{IvV&RGKl1b|?60O;Px1M3Zx
zaE3}ML+|-{0fBpSaF<rL4$5|h#t#ypqmra3fb(zlfAQq^qI-YZKxr{CQA0<F34rTf
zOZRW){r>9jIf~p<-8T15fZ~w%7RHc!qd)~Xn-2Vo6KCh-<@!H4@%xzm;lxwpKd30r
z5(Ex6b2q<@#tMlYu{2kDEQQ`I6~)R3LdBA4`-DKvC9Yy8t{9bz!p`vBJ4-C8oI1q+
zkz)Do6JqfRdpkVRcxjVj;^iin!@<wXloy_s>CMB7W``bA4G5oZe33fjG*G>tXK3K9
zY;^~Phng86KYIQX0nrNKyIp0aD8VK8anfV=nJ}v)+)Kjd1`NQjYxr3XYe%?CkFD3F
z1|hsn4BI-soYf&F?BfxHQ}?xR`Vwn6<cHW~hNv|^va7{79y4pZpC>^uc~98=1c(P5
zhgKm!*B9v5!g+C8!lDd_3Aa!#>)BiAjD4Zf0k2`U>*d7m3S0;b(F?4bj>;xhEXL6c
z>>pXyToArX_38f{18RI5#vQ=y9{Q%mit{2bh?q`CIUqa+H*{rqSbi(C`x_eucR0-|
zXeNnRlWfHS!F2dW6VJt(xPqF0hKHAj&6>r!g<--=_vUl-#rt?YKd%E~amvfib9K)o
zbcC{8NcHShptS`)g6{64w~e=fg9bo^kct*33=eyX26sKBu2p6;rox2<9oa|emM?+Z
zhgazsUSj(0!%8iShZ}~@Cow$HP>`4&w~eW5u6S4P(5l_9=bPz)mz2oAohS&fZ(g(n
z0b?b^2Xf^Aa`DyLy|_V0NumL-FVmARF*|=}&(_BW&s}1V#A`gxp(1(hNrAw?hz&q>
zGhrGz^mqo7MEo@vI+~Q2dOiH{w}*P2-0xn37hVy@lH%-jIEq7i$v0Y$IRFR}uMu}%
zlWrs8wmqq7!>&WX{EqlogdymOy#eU^N4!su!~ukO-ql}-Vm}re(C0sO`hv~j%VkFd
z^)dd=!Gp}|Ghpy!;gJ#>Mw<7_pJ?A+WBK{_V&N+VLIGmEWEK7zpGemck;JgylEtFX
ziit`+#(%y0MJb=GAYdoMbX$9y0ZP(<fG$@1MbRAR3f~o(h7C2hU8@Y+3aQfI4G&6$
zA7OjnPj((k){nJqadwFASO<P@IupKQYoH1H<Y9YbVe@5%vSAZ+fJLL6U!Y=B6@J3W
zer)l2?h{2eWqpWLR9>$(-ctz)&|A6KXZ@7;{;Qes8EP3^Sp|y7Z&UMyCR8UF-l3gg
z%Y!X52S<mt5unU`#HpdIVD&GbN65D;x3N!ouJE7wW*CQoB8a@?$msAZA6F7rl8#64
z`V+&#40Fphiuk6;c0<s+9m<TXKs9J?SZ+9ONc<5VzBH>)yBZ#S9RY<PrY_a*T6Xky
z>UOww#7#I$SchSl-6uZ>jv$v{cT7(d;fQo92HY^z5<e+XMnjg2=n*iFY@MvL$P1ja
zAcZcwOz*k?4QT;!fj4BdEwqYKrZ2c&`n;gBrrM`ZrCy|N2E_D%O)wn<H072=ZRpae
zmg!p+zc6YrK4JKYr^YnF0A<usbj}S@=TJFS(J4wS)yQxbx0gElilf&j6<pk<sFkCV
z1y%VGFjmE`It%7i;T3bKbO=30+K|L!?haZNTTu*QxMnb8Fjny^cmq+cpKzK0?`ZDu
zpo_iAMH*4E$f>e5_pN$~9})t8KFFlWg!kLyza-Bf9Yx&!ams#Wc;$4(ugaqOMzv*!
zZjx(~si1UnGv6iuXi{ckZW27P^!3$*$He$V?bi=qm&adQBwF}cuD{zG<taRZVpvBE
zkR7B%^9F7$O!;vRa?o(7n;(pmjigMCO&yQ-j}?sZWu2*fF4@mWEoc_I_B+qUQ0<YI
z{8A~ll|5t}fYuokHI#313d)6YI^8<aHk>v1HoS6$xhlH)yY9PQoT6Sa-@dpdziPfx
ze#-Zh>M7FG;-|IfM*+;8c%8Zdw*jT$Vc|w$TtrU9g6wLqFxZB;Z8%#jMY-&`HBHA1
z90u`>``P7REn2(IKIf|BmIW1pNZA7i$-g0OHVy`|j1VXgh=x85C6&V&k4fsZf<$PI
zvG}q)V`*eg*Cf{H*0!$Psl76{woSEJH23HUFv2xgHXk=WDJ)K6nXZAI&*>cNEY}xq
zy&v|lF|+BZoj3k4p7xH?sI|Ri`nRudmTzyypyAkP8)&qA>gZ`G49(cuH1&fhzYA~i
zd%n%nso8eTJ@=VkgoL@W$+DpXSuU`X%M+5V@cNy42c2Er9GmC0;<e2-N*{DTM7a50
zBVu}Brb)*1_~m^yE7%elHq)EZqb|iLowhPNmpOlF+ciJj*yT17xQe@>yzjWIKFmE^
zJNco1c*)j(`K4Et*QYzKJJLs5j~2W)y}!QhX^Zpm5Xt_6Br++2`El^0M7toBos>cE
z`KQJc^sQ%GAN)U#c#rI7Ty%tW#`}9vMxbIoD@FApGQ?HH8g2D%FzEXjk5o$COg%t-
zD4rxPo#~o|pM{(GIZMCCxhJ+q%DC7_9i|sD9)#QF)MLDswf1>6kMN3lq25lQT_}TR
z<%)Vu!X=cIc}gRt<?<4zn>K}35@$3R6wDHQ$eWS!bm(-<JEvMUug>a#^#C^w*%sM0
z`mi6k_2xxAT^zlj%zOEzG=ZemG^MCxl{NyEXK8s?c~dzp3eV){)6ZWD(lrYrCPXFL
zt4LwzqF%FCzdB(Lto&9PQkgxFlCZDCMT>^R2@#tOt!V7#m>r(*?h>!iV^3a-SB^QM
z<L0}NN^?kKvo)xnJ(eG!KlV6j)$`G6*B#aVqFq*F0L|Izzfv73$SHJH-z{a;5^RCl
z4*FF4Ra;nkF;uZK1KsOcO?BTysf{LQmY;d<YKP*7NOW^^)mi1w2hB0e@YGox<QB+N
z$~w~$$`;Eq$t>G$+O%B84XVfKd(1^wLCU^mSPiRvQ(KwZonNRWXHaJy22##bR5Uvx
z8%y>@B=9HfHN7=dSF8;!ldX`|BG#GJg}kk<-tlmfIg<)kmv<{R-OrtyZggoZ-8-5Y
zSIbjX*SfGNH8>n}ZM_<LniUY@k?6sDO~D1`KKB~DYN(ga`dam)({si1hfffq$6&<Z
zAWJd}q}Z-7+v(V#qkUl6*=u|eW>sjLSDQfbIPgUvh1ciH{j*MyK)1-Uh$;`WOI>T^
z>a1+s9i1iK$+$7v{fJ!(Arj}R+Mb-eHMfDo3(9dpHAkHd=ULQ_;7r_fh0XL7p*W|V
zA3xT|jhkg%4^x5#9gBfYwm+>#W>2=<hf{xC><POpo!yi_8zPL*Gtu2^ShxG?*yj0V
z4f&dq?<x6H{y?(81|Yxg0dX3)`)bH_+_-vkX<8$J-TqkbP1MS+;<W9OR2@}wU*oTB
zwX5Q-T1T7V1x-gnkJ|OLvF4EL<O@<yj@xqG9e|+Fjv(JpgLYqj<Arj1VaVO7XVW&-
zFT-xXA@LKMdC49+4|<@m-lg-6^i)hzUbisQozC^Q_1O=t&Nu9A`23u|bxv=_OF!%9
zx6a&+9Y7ADGpD!e#HFTs2~CP#j^9~!t?wEWKj6NqKP*3F5>E8o{zY_pbdVSJCfyV9
z!tdnze@boMf67p_zZyU#Y~a0?7B=t>(EGCWUsviDHX=5T7Jm!SfZ`UWrVx17=Y3nJ
zHh_Z@Xl-HR<Oty90@^#-LLpWVQz#m|lND%U3vX0)bhH2f%^VEh!)shaW2gI4*Vw|r
z*vZ<|3gQd^LM^OJ;38{7cw3DP0BGa@xmSdbG&DAbH<F?OO)TK994#CHKvo-Ds0qXr
z0Cac=(imR(TUi<Yo7~*V#>~*c$=b@$2?_w(n%UYwEaBb^A1piGSMIw1)&cua>;GMd
z|7n4JKmtHzCnM-XN9?^24lxuAt>GcT(ZNGr{^&?StlmSQ7RH8Xz`yeXxX<Xncg#9F
zKupot03Z+=+dm%w7boW{PJk)k9~nEmFZligu=z&@;)3I#|CDj_^8SyEhm9MKCH_;7
z2gCz!to<(;`~TsC9R%+H{!cwNcHaM%f!KIK|Lp?=0&)LaEKmnS3oD2Nysa4kRJL%1
wz~>79O52#)!qM4%K7dNLwot%*yYWMA9HE8|(1*-`K<wOXXtcBvZ@_5(4^8;LmH+?%

diff --git a/Analysis/ode.R b/Analysis/ode.R
index 1092930..75ac6b1 100644
--- a/Analysis/ode.R
+++ b/Analysis/ode.R
@@ -1,6 +1,12 @@
 ## Set Working Directory to file directory - RStudio approach
 setwd(dirname(rstudioapi::getActiveDocumentContext()$path))
 
+weighted <- 'False'
+#conda_install(envname = "r-reticulate", packages="networkx")
+#conda_install(envname = "r-reticulate", packages="matplotlib")
+#conda_install(envname = "r-reticulate", packages="pydot")
+#conda_install(envname = "r-reticulate", packages="pygraphviz")
+
 seirds.f <- function(t, y, k) {
   S <- y[1]
   E <- y[2]
@@ -18,7 +24,7 @@ seirds.f <- function(t, y, k) {
 
   # Saying infec rate of S in contact with E same as contact with I
   dS <- epsilon - (beta*E + beta*I)*S + waning*R - gamma_d*S
-  dE <- beta*S*E + beta*S*I - (delta+gamma_d)*E
+  dE <- (beta*E + beta*I)*S - (delta+gamma_d)*E
   dI <- delta*E - (1+gamma_d)*I
   dR <- (1-mu)*I - (waning+gamma_d)*R
   dD <- mu*I
@@ -27,11 +33,22 @@ seirds.f <- function(t, y, k) {
 }
 
 library(reticulate)
-#conda_install(envname = "r-reticulate", packages="networkx")
-#conda_install(envname = "r-reticulate", packages="matplotlib")
-#conda_install(envname = "r-reticulate", packages="pydot")
-#conda_install(envname = "r-reticulate", packages="pygraphviz")
 source_python('prep_model.py')
+model_data <- prep_seirds(weighted)
+
+S <- unlist(model_data)[1]
+E <- unlist(model_data)[2]
+I <- unlist(model_data)[3]
+R <- unlist(model_data)[4]
+D <- unlist(model_data)[5]
+
+beta <- unlist(model_data)[6]
+delta <- unlist(model_data)[7]
+gamma_r <- unlist(model_data)[8]
+gamma_d <- unlist(model_data)[9]
+mu <- unlist(model_data)[10]
+epsilon <- unlist(model_data)[11]
+omega <- unlist(model_data)[12]
 
 # Obtained from prep_model.py
 seirds.params <- c(beta,     # beta
@@ -79,6 +96,6 @@ plot.seirds <- function(sol, method){
 }
 
 plot.seirds(seirds.ode.sol, "ODE45")
-ggsave("SERIDS.pdf")
+ggsave("SEIRDS.pdf")
 # Sanity check: Make sure sums to ~1.0
 sum(tail(seirds.ode.sol$y,1))
\ No newline at end of file
diff --git a/Analysis/prep_model.py b/Analysis/prep_model.py
index 5f30a7c..05194f9 100644
--- a/Analysis/prep_model.py
+++ b/Analysis/prep_model.py
@@ -11,104 +11,136 @@ print(os.getcwd())
 #os.chdir(os.path.dirname(sys.argv[0]))
 #print(os.getcwd())
 
-# AGraph preserves attributes, networkx Graph does not.
-# Many of the desired functions are in networkx.
-# So import AGraph to keep attributes, then convert to Networkx.
-A = nx.drawing.nx_agraph.to_agraph(nx.drawing.nx_pydot.read_dot("./1_mo_color_DOTFILE.dot"))
-A.layout('dot')
-#A.draw('tree.png')
-A.remove_node('\\n') # Remove "newline" node from newline end of dot file
-G=nx.DiGraph(A)
-
-color_map = []
-color_d = {}
-node_pos = {} # used for drawing/graphing
-
-# Compartments
-S = 0
-I_R = 0
-I_D = 0
-E = 0
-R = 0
-D = 0
-
-ep_tmp = 0 # counter for epsilon
-
-for node in A:
-    color = A.get_node(node).attr.to_dict()['fillcolor']
-    str_pos = A.get_node(node).attr.to_dict()['pos']
-    coords = str_pos.split(',')
-    x = coords[0] # layout for draw function
-    y = coords[1]
-    node_pos[node] = float(x), float(y)
-    if color is None or color == '':
-        color_map.append("white")
-        color_d[node] = color
-        in_edges = list(G.in_edges(node))
-        tmp_S = 1
-        for source in in_edges:
-            tmp_S = 1
-            # If previous node was infected, then we are recovered
+def prep_seirds(weighted):
+  print("Prepping the SEIRDS model, using trivial weighting=", weighted)
+  # AGraph preserves attributes, networkx Graph does not.
+  # Many of the desired functions are in networkx.
+  # So import AGraph to keep attributes, then convert to Networkx.
+  A = nx.drawing.nx_agraph.to_agraph(nx.drawing.nx_pydot.read_dot("./1_mo_color_DOTFILE.dot"))
+  A.layout('dot')
+  #A.draw('tree.png')
+  A.remove_node('\\n') # Remove "newline" node from newline end of dot file
+  G=nx.DiGraph(A)
+  
+  color_map = []
+  color_d = {}
+  node_pos = {} # used for drawing/graphing
+  
+  # Compartments
+  S = 0
+  I_R = 0
+  I_D = 0
+  E = 0
+  R = 0
+  D = 0
+  
+  ep_tmp = 0 # counter for epsilon
+  inf_ct = 0 # infection rate counter
+  
+  for node in A:
+      color = A.get_node(node).attr.to_dict()['fillcolor']
+      str_pos = A.get_node(node).attr.to_dict()['pos']
+      coords = str_pos.split(',')
+      x = coords[0] # layout for draw function
+      y = coords[1]
+      node_pos[node] = float(x), float(y)
+      if color is None or color == '':
+          color_map.append("white")
+          color_d[node] = color
+          in_edges = list(G.in_edges(node))
+          out_edges = list(G.out_edges(node))
+          
+          
+          tmp_S = 1
+          tmp_inf = 0
+          for source in in_edges:
+              tmp_S = 1
+              # If previous node was infected, then we are recovered
+              if (color_d[source[0]] == 'red'):
+                  R = R + 1
+                  tmp_S = 0
+                  break # No need to check the other nodes
+          for source in out_edges:
             if (color_d[source[0]] == 'red'):
-                R = R + 1
-                tmp_S = 0
-                break # No need to check the other nodes
-        S = S + tmp_S
-            #G[source[0]][node]['weight'] = 3 
-    elif color == 'yellow':
-        color_map.append(color)
-        color_d[node] = color
-        in_edges = list(G.in_edges(node))
-        tmp_E = 1
-        for source in in_edges:
-            tmp_E = 1
-            # If previous node was infected, then we are recovered
-            if (color_d[source[0]] == 'red'):
-                R = R + 1
-                tmp_E = 0
-                break # No need to check the other nodes
-        E = E + tmp_E
-    else:
-        color_map.append(color)
-        color_d[node] = color
-        # Check if node dies
-        out_edges = list(G.out_edges(node))
-        if not out_edges:
-            D = D + 1
-            I_D = I_D + 1
-        else:
-            I_R = I_R + 1
-        # Check if imported
-        in_edges = list(G.in_edges(node))
-        if not in_edges:
-            ep_tmp = ep_tmp + 1
-
-# Params
-beta = (I_R+I_D)/len(A) # rate of infec (I/total?)
-#delta = E/len(A) # symptom appearance rate (E/total?)
-delta = 1 # incubation period
-gamma_r = R/len(A) # recov rate (R/total?)
-gamma_d = D/len(A) # death rate (D/total?)
-mu = D/(I_R+I_D) # fatality ratio (D/I)
-epsilon = ep_tmp/len(A) # infected import rate
-omega = 1 # waning immunity rate
-
-
-print("Model Compartments:")
-print("S:", str(S))
-print("E:", str(E))
-print("I_R:", str(I_R))
-print("I_D:", str(I_D))
-print("R:", str(R))
-print("D:", str(D))
-print("\n")
-
-print("Model Parameters:")
-print("beta:", str(beta))
-print("delta:", str(delta))
-print("gamma_r:", str(gamma_r))
-print("gamma_d:", str(gamma_d))
-print("mu:", str(mu))
-print("epsilon:", str(epsilon))
-print("omega:", str(omega))
+             if(weighted == 'False' or not out_edges):
+                tmp_inf = 1
+            else:
+                inf_ct = inf_ct + 1/len(out_edges) # trivial weighting
+          inf_ct = inf_ct + tmp_inf
+          S = S + tmp_S
+              #G[source[0]][node]['weight'] = 3 
+      elif color == 'yellow':
+          color_map.append(color)
+          color_d[node] = color
+          in_edges = list(G.in_edges(node))
+          tmp_E = 1
+          tmp_R = 0
+          tmp_inf = 0
+          for source in in_edges:
+              # If previous node was infected, then we are recovered
+              if (color_d[source[0]] == 'red'):
+                  tmp_R = 1
+                  tmp_E = 0
+              if (color_d[source[0]] == '' or color_d[source[0]] == 'white'):
+                  if(weighted == 'False' or not in_edges):
+                    tmp_inf = 1 # add 1 for the inf counter
+                  else:
+                    inf_ct = inf_ct + 1/len(in_edges) # trivial weighting
+          E = E + tmp_E
+          R = R + tmp_R
+          inf_ct = inf_ct + tmp_inf
+      else:
+          color_map.append(color)
+          color_d[node] = color
+          # Check if node dies
+          out_edges = list(G.out_edges(node))
+          if not out_edges:
+              D = D + 1
+              I_D = I_D + 1
+          else:
+              I_R = I_R + 1
+          # Check if imported
+          in_edges = list(G.in_edges(node))
+          if not in_edges:
+              ep_tmp = ep_tmp + 1
+          else:
+            tmp_inf = 0
+            for source in in_edges:
+              if (color_d[source[0]] == '' or color_d[source[0]] == 'white'):
+                if(weighted == 'False' or not in_edges):
+                  tmp_inf = 1
+                else:
+                  inf_ct = inf_ct + 1/len(in_edges) # trivial weightin
+            inf_ct = inf_ct + tmp_inf
+  
+  # Params
+  beta = (inf_ct)/len(A) # rate of infec
+  delta = 1 # incubation period
+  gamma_r = R/len(A) # recov rate
+  gamma_d = D/len(A) # death rate
+  mu = D/(I_R+I_D) # fatality ratio
+  epsilon = ep_tmp/len(A) # infected import rate
+  omega = 1 # waning immunity rate
+  
+  
+  print("Model Compartments:")
+  print("S:", str(S))
+  print("E:", str(E))
+  print("I_R:", str(I_R))
+  print("I_D:", str(I_D))
+  print("R:", str(R))
+  print("D:", str(D))
+  print("\n")
+  
+  print("Model Parameters:")
+  print("infect counter:", str(inf_ct))
+  print("beta:", str(beta))
+  print("delta:", str(delta))
+  print("gamma_r:", str(gamma_r))
+  print("gamma_d:", str(gamma_d))
+  print("mu:", str(mu))
+  print("epsilon:", str(epsilon))
+  print("omega:", str(omega))
+  
+  return (S, E, I_R+I_D, R, D, beta, delta, gamma_r, gamma_d, mu, epsilon, omega)