From 184be47a3379efc30d913de337fa10cbc90f8007 Mon Sep 17 00:00:00 2001 From: noah Date: Mon, 1 May 2023 15:03:31 -0500 Subject: [PATCH] Model v1, with call to python --- Analysis/SERIDS.pdf | Bin 0 -> 11569 bytes Analysis/ode.R | 50 ++++++++++++++++++----------------------- Analysis/prep_model.py | 7 ++++-- 3 files changed, 27 insertions(+), 30 deletions(-) create mode 100644 Analysis/SERIDS.pdf diff --git a/Analysis/SERIDS.pdf b/Analysis/SERIDS.pdf new file mode 100644 index 0000000000000000000000000000000000000000..d8cd1e262d15ae27005846c4cfcae75f87396c1a GIT binary patch literal 11569 zcmZ{K1z42b^EQe|DBZn)gw(P~r_|CY-OUoxuqYkUEgcflNOyNECEYEJq%;!pE%?6g zum0a&7t8g`GiPSboOAY>`?-hqtvHwk#L9t2>$B>!<}>UwuH%3P0Dz zI|AM~LG7HN@L_KqY)zbuArAL%J@_ffF_6Ab{U;$jB@ z0L2WUhE}#_0N`6gGl(PL6f|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)1$uJT(wof9C>xpSlZ&ENikro z9WZ;f^+sRV^ZM8};O6e;&`pEtCI%>k^?Q~6Zf|MrXc^dgRULIVDC>5!@ca5!Ry*Z4 zP;WIoXF+jkW?&5$LBauL$ghs~r9FkWK<^?!a`8rrH~A{=P}YD{mq#6GV&dsqr`!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>;ToyEZo6d)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=>RvX4IP37Ln*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&9Y4@O6hDLd+eBDo>Y86U4Uy>(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!gUUzzgW3A)aIP<%gvW{KKa-$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*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%-7g6O`hSNWO*T<=&oRRoAYnvk^;5iLhSV;*u;<8I3E-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#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;xt^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%RXcNqIBKsbm8pE@b87nr2#$Bx5E$+&8)@j^EDkMG3Mf4-B{vux?hc_sx9TfRO=yINh2Cd>Hd4E4zFn28&a58ldFx|u`B4c zDKya&=UI=2Hmkj<^1)MW3b`TvxG7pF#Y>D?O;=u3F%a9cH1pQen;bcta-Cw^26y8C^$Had{e#AYynV=_W|qHp1AT08?kSM=xStXB~?W-kVR>&{rB~=Pt@T7>SWX zB5nQ8c}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}9n=~6@Cbd6$ z4Br(N;-0tmgVfjrfl@CzRxPKE$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)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{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;hHgx6+HSz&$SfzE%i|LFMKf^5`}J1M&P-PFQ#N_gqIdpA;|3c_Ut>nj4)RKxHBsb9ch+@bi7;K+jD z@X~^$#=O<9e=Vjah(xRsb}BBSWDQePCD%{|DUpSjmsNq^IUYah>v~%JZIRRV+ds-$ zwJI7LyrQ%@?Q%ei7xkOpo`OXi3oj`cJavGVxn~1n<953e&Qh_+<3t&pJ}8 zO=UW@ekJOT4_Emdq0vXfxid1htWt=TAopai>zaIY`&7~qC!BR(|MYit_akwfa0As* zDfLd0wDQ+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#SZzhLc+h{(*C!jd!|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`uXMyb&oRm!5zk^bf z*Eo(x$mkQzHrp89JcbiN4Dt_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;33P90-)&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_nNZH-F-)UT242kP=qUf)h>|U20@?7#|x-?usXD-NxJDc+=TMnJm_65n$_Kzq% z>uw{PHs^xITwHo{1aN81nuRwrVeTE3Z)hkrOCbXF23_2-srv?AIpu5NYuV zWADD$;rvCd!n4Yqt@}dWBN2LO>a>(z&h4AGXEsNk9DDjyXN$dGWEf1x5%~{l~~8^gLKVu!R$C zw=z8b!3rFzq&nXPz0}V~uzH8%NoU8HrzIFf)g*))<>zv2zYFzEo`>Lt@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+|-{0fBpSaF9jIf~p<-8T15fZ~w%7RHc!qd)~Xn-2Vo6KCh-<@!H4@%xzm;lxwpKd30r z5(Ex6b2q<@#tMlYu{2kDEQQ`I6~)R3LdBA4`-DKvC9Yy8t{9bz!p`vBJ4-C8oI1q+ zkz)Do6JqfRdpkVRcxjVj;^iin!@CMB7W``bA4G5oZe33fjG*G>tXK3K9 zY;^~Phng86KYIQX0nrNKyIp0aD8VK8anfV=nJ}v)+)Kjd1`NQjYxr3XYe%?CkFD3F z1|hsn4BI-soYf&F?BfxHQ}?xR`Vwn6^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($(-ctz)&|A6KXZ@7;{;Qes8EP3^Sp|y7Z&UMyCR8UF-l3gg z%Y!X52S)yBZ#S9RYUOww#7#I$SchSl-6uZ>jv$v{cT7(d;fQo92HY^z5e5_pN$~9})t8KFFlWg!kLyza-Bf9Yx&!ams#Wc;$4(ugaqOMzv*! zZjx(~si1UnGv6iuXi{ckZW27P^!3$*$He$V?bi=qm&adQBwF}cuD{zGv1HoS6$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;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}gRtK}35@$3R6wDHQ$eWS!bm(-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 zG$+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_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=E8o{zY_pbdVSJCfyV9 z!tdnze@boMf67p_zZyU#Y~a0?7B=t>(EGCWUsviDHX=5T7Jm!SfZ`UWrVx17=Y3nJ zHh_Z@Xl-HR~*c$=b@$2?_w(n%UYwEaBb^A1piGSMIw1)&cua>;GMd z|7n4JKmtHzCnM-XN9?^24lxuAt>GcT(ZNGr{^&?StlmSQ7RH8Xz`yeXxX79O52#)!qM4%K7dNLwot%*yYWMA9HE8|(1*-`K