From 66b579c87b3c840f8293153f81fd3d737a148fd6 Mon Sep 17 00:00:00 2001 From: GitHub Actions Date: Tue, 7 Nov 2023 09:11:40 +0000 Subject: [PATCH] differences for PR #55 --- compare-interventions.md | 2 +- ...imulating-transmission-rendered-plot-1.png | Bin 14897 -> 15072 bytes md5sum.txt | 8 ++++---- model-choices.md | 10 +++++----- modelling-interventions.md | 8 ++++---- simulating-transmission.md | 6 +++--- 6 files changed, 17 insertions(+), 17 deletions(-) diff --git a/compare-interventions.md b/compare-interventions.md index 49d9676e..0f7ea821 100644 --- a/compare-interventions.md +++ b/compare-interventions.md @@ -64,7 +64,7 @@ The diagram below describes the flow of individuals through the different compar -See `?epidemics::epidemic_vacamole` for detail on how to run the model. +See `?epidemics::model_vacamole_cpp` for detail on how to run the model. ## Comparing scenarios diff --git a/fig/simulating-transmission-rendered-plot-1.png b/fig/simulating-transmission-rendered-plot-1.png index d025a2355fcc941d1d60afc347e10414aa0c077e..50c38b6e3f43c96f1468e59c3caa66a088fb0773 100644 GIT binary patch literal 15072 zcmc(`cT`hdw>O#)Aas=8QBjaCU3x;d0Meza2m}z2E};b?BGN=e1O&Dq(m{|80hQi2 z(z{Y4RUp#K7oPWh&ikD2oH6eGCz=WKEBJBFJHNGg`b~aKtSN?)vJPnfgwtm8XB6Knp#>~+S=MWIyyIQ+_-u3=B-<|Zr{GGtE;Q0r>C#4f9K8} z0|Ns?Lqj7YBV%J@6B83tQ&TfDv%7cinwy(jSXfwET3T6ISzBA%*x1rIoa&~rhadCO{=#i_dtDBpfySw}2$B#WcJUl%; zy}Z18TsPHi>Rom=;-K}n3&ku*tod3moHz&$H%{V^(rAD;q~j+ ziHV7C-n>alN_zYDZE|vQN=iyn%=W@ctqR@R3PAF{Kv zKYsj}larI1o12%Hm!F?sP*6}U1kU%q_(`t{qlZ%s{2&CSg%EiJ9Bt!-^> z?d|Oy9UYyWon2jB-QC?iJw3g>y?uRsAPDyN_kaKXePCcMTo1dRwSXfwGTwGdO!sGG3e*Icr zUS3&QSzTRSTU%RSU*Fi+`2G9$=H}+s*4Fm+Hi1C+^XJdb&d%=c?%v+s{{H^K!NK9- z;nC62@$vD=$;s*IsmbikNMh(yKGre!f_4N5AF^Xm#np+eo#JP(M@f7EyFI?ypx}d>=Johby zwDuhNrKq+nzt*|yu~}ZTR&Day>yMCkUS7F)CuFj?cx|dK!P{)a^3Rb9uUhR)aWQ)N z_+ygAv~#A~1)iG;PsAal3_TFk^BfZ3VmcM(*%KiM$86r&00S8~lpf_LT!$a)436f( z3-z_sC8ljYvR^DXhjD1GON62zYedTb~;rI*!d@ zX>qHRATKfrAIfypFR&{W@0AC2RsQatY-`}NGR*X~%=~4!wP4=*&9Be$g%(&=Vp(&V-5N^a zF(aAwy{qMC0wln@@+KyE5*hDBfF)h(lbLPhx!8TXB1NHqR<{NTcoA9NmAAX!Rp_an z`+Mkxwst{5tuSlby+^-KU#lOi`5CHw+Q54XIc)ooNSod=83NNxkWdNnUO^(?0tG+4 z^+|vKM@v_!@%~2G)7|SlwfkC)KTe0atGJ(!1#&}TMhq(ys%|o#+n(hj(R~9e0nuUB zjhkQd^H1@=vwc^Jq}QlF`$P-^krw$pE7PAHQ`ar+Y~HB@9iEW#^m`*A<3S@yM*PMm zMThu=>GxIlw`23g1ngvI-eUQk%~{9GnFDV``Fqg0MQ9N<7J(wa(?|HD@`qe)it!lJ-7CNC=>&8%p0cU0JNMKJ^%UQ zYbWDlbrH8YxA>xZ->|6qxUFk;r%UF3`su4`T&smjSsAaj%~;EG;etrCy(-m2T8YRD zSkF77M%FkoLt(cDApLx02JUkDiDz2%@e!HF;gIm6d#stSxl#~8h!39P6lgy~InhXM zG?wi$yTkj<*woVQgrz0ap*x2ZtcRd%AmtYzDXb7$qE0*z(UiJscEmOb;131EASf+J zIcj~@0n$T_StLVHlK|{cFta8q6)G)$wnmf!<3NgtS`>35{|~G4|DQ6s#Exn?VLlV> z@4<@F^aU%liu27}!fV9bij7*<%t|6BO2Pw96g&=nR`zk0<}t8DZdJF_bt`+bO!)W_ zXQ&y<<_B{Mt676usy^kJSwx$W5TgDZ2pCEkVj=D4A%$7g6ars_Jtsf2*#;tgW|b5P zpy|ab2SeF>fVoEg?=xp?EY3lCC@?HN){`{<;m-FrU$Qe_l0m`fLC$g#fCwE#?*Ce& z14E`}LjU$X(XY`ERR3g4&zWC|?q>O?27j;le|FSopT>2L~BI+zQ#qsKY z`Ri{pl}Kydx_4|N|)JNcswyT5)V)LcMM-YWVe6?Idsc8 z?_Xjhtjl3lThmT*QM%LAGbAWsh(#PJ5Z7%q)Wz=8R`xl<(UEto^2(Fl<`QAv<5l#& z^x8!)FHe438N-}mYv#q;>irZUw~{&*hl6WY^HNRw+a)|U;T)C*BH1Loo^E`ANYHBBqQV`TUg6^Xk7T46{DQ3Nz0my z&+Ja^3K7WitB-5{91&SAiUip2Y-!@kI0%aY4`VZ!^m&WR#;+w>X9-cm51Q4s-;f_P zjeVuR$hXqaOC8qx{A7YPaQ>e5Y#dkO?60)bb`P_EvxC}k+&lmO#=1QcrnA6EYVl?o`VPlmRjY$5J zh_{rG!ry-VS0?^96s(DY1N;9yHUEyxO(6U~`Sp)Qm#9&lQW%!M5+6#8&XoVc!%K3g zc+?qY5}_1^f)v63gPecI-rumH=?RkjCv6}W^nU@&Z9l~Z{(nK&{|L-ip{=H5n8lVg zVE=Em*66_Ie8~TmIy=m$qb4TIA}NA;Rt&srZv*dPraeO_VutNvw9a^(X{W_S9z(2E z@|#R3Pm?E<<+GY#xg#Wo9FZE&E3Ge$VbPEUeVzVY1IR&$0@J@iiwJi$Q3AZE4`*A# zfLT1lR2|ZkLK=kEUqDEN-H3lxm&XXpf^S~HS|>zWz=ovgWdX?b z)-vxEiC$hOwT`dmqd@zcW!@DbLKr)eGtr-ckRLl9n^@MIc4o4xrXJWg$vSQ{x?~S2 zBD0ixA8@u1qWi>Dl0O{br~(Bleh6#V!4ccDdZ3_q1fs7VF3v(kca57M+TMngTW6es z1Xs(5872YNbSTj)F`?)1Wgi2e-ZB_1zAN|^9P{z-elcO_}c65gEjW zn*|j3)k7cxmR#8Ms}F2I;ymqS5-~QPe*l=>!-6Z4$P|2_M8HNpBt4!~X2C~f5cilP zCK_mSv;?-bNKmbU&5#sMxiB7u9unnLu4Qj$Rf#GzQ*FV|3hsiaF4hfgywJ2nH~^JqphW6(zQCkQ(jcS|jVtHDXjPMq zP@aS+OuC-{(PEv;q{sJENYEyZAIMU9{?&JoO4}GVJ|HD1M(!X^Ameag<_`~Etb`AM z^~^v)$aw5iASOt#RTI0)*a(&WzDC3*D#7rc z-!wpt2QS{Bh7JE`4N~5SP&7$40(jfo5F0mgKTUwYj_E8$H;5^@gq$&m4mEvMS!Vn1 z6k0<1*@j}*OO@z=(<8%T^f8RcJYM>ocR|l}e4Wi>9zuKu!720PSI!~STM40PoU9HA zGwl2XB9s@D@C{eSV{m!J_2Lz+sm_(2m|Iw#*;g_NL<&J zWsxOocqRPC_Y~6YAF7kI7T}p3ZAg8t8<^?jC!hzIe=vcIXvmjygjmx86#$uia+10% zk`pFqkZDJ$0OZgh+Ti7UiSndq9Vpb|S9PlgsH?WWM4E(M+*KsPq6aZBZ-7;zh>4v^eC>nXRgjmVL(GtN=@jSMb${80)nx=eF4;ZS1DChxc zx!vp_oMGn`(8G5^*)MTwSnce$5S_vVhaE*cIgV0XGpnr15#fPe2w_b_AeWOiNmlG{tA+uIaBAF5qiV5BHZc zcu&lBv5VgIT)^C6&N2=Dyb%h-|Mu<0oLmzuPqTZfh-|Dphwz${1+i34;6bAOjs8;7 zsgQB^)YIUmV@E4Wg+|vlN}}tlLiMb5u)Q4~9E3jcTsd-8HOIR1gx)tURpq0`HwtMl z5ovOee&)JZ+2Tv{7w}I%J4f*jS(v>G)7;tnaL8mo$Ttn=-rnVLL`*bwjsi>B^0Wsh z-P^m&$LFfQ#sGDqrmtf}>7XcIEArYSR(omug7Q1^WBU@Rjnl##2@hDq#GHy~V}>0s z)+qcY<~SEQZcS^&jFYfvs)rrk>Y3JGUK@jgM!aCKUSM-7VB0=EI~pOVjo@=^9tXNs zlUySZye{n?T)h~cJu8%Q#qMRi{Fq3Ql{bm;LW9?~>fsk>ts5`+6?p}pK) zBt4TRY%Aj7G8=UHr!SArp-mbZtL(;Qss+31)6qk^5m9Au`2hh z1P`VD((WsP%()xNWzE-vS-XVOlrA zsW;DrQeXbM2u`gE2k-&ObyPu~G0gQ_HGV0;cEn3zQerW<#n`2mh#R-Q^Y#%CFWO6i z@-(}jt7jX=bNZGN?JxHr=wAGJLRbUZ-aaxr9qJm_85$~f*b)rk?)*W*NN8|4j3kPp z0-Fm4pT?u@Z_bcS);@&c_$0oZS&hF1D!)kv5_?MmXuJ#m0{>(-JOlAsNrHY z)i9NN4V>`0gC#}?wUc`A*QbCVzbRXy z>m%kLEd1;!??W!Cl!^i{4k%jbL6z9S+GWBt)WBQl z+QsZWFES@NDRQ_Lw|(=$*?RBgz=AS!?7Q4=s&UG%mjh##h5I%7nSG?SL0*BGb`B3u^_osNIgmth|PV^(CO4WU!gstCU|=A z?$FNS-kqwd$=ihtid&%#jx2`e1{@{VY7lla$JhSF$knz~VsI^pg|VPHHV}w!DrT~` zR&>YM3?;%-N(;aU!1e2Uw~%WYRsyCUVBFVwMF)h1wdK5;m-UCB;P87YCRA&&8XP%7 zM|@yur8a%oQ7+Nq7%96gXa`dmV+P)T}#3dYtodC>mRY;9toOm%Mkg9cobo9H*j) z-2>xcCl@3QyMFi$(KsZT4itDsc`l}#_C@b&iP}C4n?82<-BvxImC6XPj>AsFYjNe~ zWTkqAI&cL)X^Uie#OqZcLbU~Eo&HHY7C3cWtCDfoH~Jrc;6^1kU9v>qmI?}DpY(x0 z>f;{vs&C9W;i&7_|0K7`j+W62H@m%=~Q5r&RocX z0D)-1le+w1rL~SQ@j^w_GzEvR{slcykb9VcxTEUDOj2c`Vop(JTB{sdACoA($t1gJ zC>8w5!`eausOV)sYcqYXOY@n=@{A!G->G4GIXlE|3lV#OOqn0RHZMUL?)Nbp~q_8SRkc7x3emD-TqeMK$*2b?RHDV9E41 zxb7_cs=_{=L33<(|6HW6%6+sJST{18DqQ(;;rvCc9GY*Qc{ky4PAJa*uq0# z9R2z7Ix5?m(sv*NfcJesBxnzHQgO(aeI?a(&l%!l#K$xCQ&E6b4{{v$&U?>jGND!> zo{>oRGJO+Cq<&}gnKnVEH02~D(tkYq{eX#!InA6?Ck~cmaq>7)puffDL7ak1SbA4& zv0jB0un~B^2m?uzkHjlBh8F7}8Z_CaY~FmSE)QktK=9Qm>4F7{9rU$Jl2CjkUl1Mf z&rmz<=ZXH%+n32>R#~NmC#INGM<+yd`9arGQ5jOyRyEm7s711hx(eBv43u z;MpkAlWexwR@+a7ZoJppU_k%&N@3^g^pYENDQveBNOjYnE=YlUf?TFmsxP>~D}wm< z`kxV5+N?Jz4~lElK=G+{RZyp9yL79_EhCH$epMJnSc#_{R}sBs2;zQMu`x*dm%r#k z^s$ApxcdZ1G&t$#f!3y_2-Y8*l+oWWlT7ibg3m~GgIE5rqXN0$1rsB)0tq)M6~I76 zYJxCsfn3)u^b#^*O$8Q{@%6sXA z?6DK9#etTIscE#)qq-?5I9eBNJqA6xC-V!z7~u&5@$aVsDL8mbzpVRi4N6lY`x;}+ z*?0Oa)WG>_Qzq1d!k$i}*`-}JJPW=$M5oFw;z>It$6-ShhQF#b;$s5bdz1St)iIpczd>T({k8etLY^#l7!FK%}zp35jmt1uwU`H!y%#Yo`DJ zGJ@}_fji$9f3gZ@AsGWqiA*SEZ*$aoteE~n662Jg4XT_mn#j+N7n!AfMxq-%vm_2S z_kZ=RvNM}WWb`&8K!iyNXac%yE0n)+>JH#H`xle%d|V*&Je!ljx-T4IXP0{kNXAOH z!)il3*uYV`D~+MO=rCByq>?x<-f%vUxll6uSRomipx$H%B#MoeHWC%(LO=69OxtEs z8h_-=zpt?;F#knm^R^SF97tuIA`Q^r zW4(Jntd-s|S1Fq!V+IDqJ4gE2ewfJUe#zzi^&St(^NUs>IEY zj^!wF9NR?14`Fcrio;7c$5>}P8f@SCc;7bW4p{j815`IgSPeUkSv&c|jM{2uDxcUa zpwOmF+d|f_u+=ZA1C{&S}-*=;t7Y$%r}i2loM^NGK|meHDX-PLqRe>BpnYJ8<|S=8F1Pst`E4+( z|Ga$=S*@_dEbI8%t$>@f95ZZD*pDT=BJ`*uG}9w+pua?`fh4GDVlM?zmnK7!rvbFt z)pb1vvwL#ae9JH@I^1L(keD>xUjpDB)JdZX?lu}Lr&TdwM~5CY3IuPuIM=0#0C``H z%z=^jGrUKn9I6*$R(Dd9&KDVmw(b_30}c%*)*eESwcv*0Fai&uj47A|cTJP@Al+ON z%TG?)iX1qb*9{!8$C_GP%*ublTW-Db85_wQ>Z7=0EkC#-8Okpv6Nc3iyN7zsPmTYE zyllgQpDr0~bOAIoaYTW=FIRnABWJcS{5wQ>QR|m(C6EywZspqWqX_l1s-gT1*UwMs zy*n?)irQL%(xGL1>WMl>{yb!s#Y?oOx=@rO%99W|@d0U{+-6cVT&uhkQ32-tD+PBF z1p`vv%=LW7eAFs`5Y`*=z+tT4(@YZ+CW|bo?$Au&r9R82GBw4 zMT%Nxx@V3sDnTls!o2;zsjT>MOWT%5AdJ*0i2B>aXz+vcx2{3alK!;pHxL_HHrTg7 ze?mqOY07wsn)l&tPpJwGK(UC$x;6{%)_s?w zS)35N7MF@|k7>R$y-JB%%_l{V0&7qqBBNS_^3a;S?(WaE{*@gz?8jCKKBTbg@1_Ky zYF?vb2!s`Tb{0TitDx*|7=PT>uqVFkFI;8R!vrO&cb7uLV!|*XS*Bko>H(S3-~kEkaE%BXc^^+`Plm1);uEcn8?|q8rrt93eVV)a*vX3ww?8eoxqZr3 zTxPLdb*sOm+oCf%Cs#}{u-PhSZ|Zo_GLgOUw@@9E!)9S;x9ZSUOvA(Z7wCmo*G8Iq zD!w}AOEc(PH@>kz8vHqP?EdOb2eBZQdpsXS0tv{fXAGJa{}VD6MolVYn)LrXg2u(Ra_J zOi2Gl8;5U~cC2Lw_71-Db?|7sNPmII8(#U`gxRqD*yH#nh3Jk%Uir)_&qv_Ae}A`| zzVJ7IX8u z_EIO2NCz?Gej~%xgqtVZBI4f1X?&;);%`X1I=c2P*pCwBoTN-mC|fVh6M05_;=Mw} zNMQ={SIFgkLun)^ff96+IJQ^;O|LU*k^vXf~KfY8Y2;mnLSPB zfv^Kpw3PJJ8WD@$%UyvgiyrA6-%}{OBIKn^>_V%5EjFrFBWH5o2?U?n~an>J`t5`3H{AZjsa;QiapNyw9Vj-=`vZC0QkXRHp2$ zKu!zx#C=lfITVToUo$%Pk(e)c!#-*!4Shz?_q76=*1sm7DpDz!TAg82P9s|#sMX7 zEm4vb^K@;5RI9x`Rx3dHMTs5p`R2AsZTV*!(7wm_d}0l;tsK349yQnOn}M-;d$WX# zA7AFseyfHFfXNwVO@2lSE3fq{K7_c^^I#bkunjJHxkL?u(N|rF&kpQL{0uP;t`&Az z9tYdEYUE%?Gjr!3vPEn!<(l9wC2CeAIBW78D!8NzdSKFOswj&OJhXftg=};uuy5ZZ zJC6!7%hg9nRgFs^m-r%{+#tG~q?m;tZ?dx<0+X(iWyLsf7v2#i!M(58F%1uMNEDDm zc`f{C)3X>VLzHMd`w}(1a+e0_DAj9q4_o|6-0C6}Uo-!9X1;x+PcB?RA zrr$ALKuC3NljAS;T{DZNOSvnaYGFR%V-M*;Z~qK69Z%6_tK;={l0p}x4U0jzb;LHQ z@#!U+bP6x)ZTP{OxYC<*n1&Ih)g~Me#}QOnR-mTo!iq8d?zm6ttvO4#*m*z9biI^r zCmu@p;kW`l5cG=3J4E(z7>9rfl1g;rR`P?1q0Hy;g2Qf}>JZS7QJVJA_|pk}4>+@d9;*teV5r3z%N1zH>;}H*cl!XS z#J|pmFk0xKA=)im_t=yU?j0@82u*5exaF7p=bRfWgBUF_& zSiqJlgX~CTiN1;s#$lMn5A1N4r67Ep4u(vn)_$aerdWx)R}$ykR#)y8B9^vfb+!=j z@sd|iLM11i&_e?JIda3L;grcsQkVx4!&An~j>?|9c>_TYXn$wMno9B{K?z63zSchx z!pF7*6}(BkC?XI{KJfW!O|S4Jz_AjUP{F)Q-1w*Q@5G2^E6e%6-9<_@FmllU2(x)i zO<5a1#lr`tT7GI|q(W>^TU*16qD?I75yDy-*?orXY;hDYwDp!pUzn^}@0|lHp1Jc@ zp#8xQ(PA>!V8yphDva8s7VLK~K=HCG9WAlv%D=~R7;|J`;RwI?aoSQm_d`Dpcj5L8 zF})mpg1V4{G3YuRClh&~hf=x3g*RY)VzeC3ONeduw22G^!;iBeQvnxie*8pkpc+yk zJ-MR2fNeF_$m|IyT)X{okWUpqYrNk16FWT>7{SmpQaXy6PSW93bT;O{PID5gHi8X~ zuF7~xqHD6-^>EXXZXX$!8OJ2Y`xRdFVd)kG?IG@5!$*8J@TCG^&gyepz->Q;Tg$&X zu=~OzE{<|psK3p^Mu{J*iU&A{)+F|g?AJ^sgc&hnxchJec5913qzim``_3~mT#AAx zJfNNc_zA54WIiypHJ5gp8P^K4kk5v zk^+hEznDg3r8Nxc>Uxu9?4dcMH1xHhqOWMD*L<*DRV+VQo z>Lf4Xdv3t>wmufR+e;8ZOtGSQ$L16&f_k`}rlljRLuGbDRq|QhFC-lhi$x z_6+Or==HmG2T(nacp|4w15&G$T>C8J7&V66&C$}=;^8ALdgxHY={D7<}xCX&I5oco6f=u1| zJtRnj%I;J8>YQh0Z$Row82jFvX&*&Oi2cwrnRE8v8`0l+9&5Uu@#E|;L4<7OJWrNd z*hd%3!8y#Q@+V>T05yjr{9!E?sB~f8;-hOD-g*Q4+~%-Pq~P`?&Dz-D7Hyj z3FQ?-MksSuWzII_jtB4Y9#h$~D<6A;#Ci#Dv4B!DeQ{ zxn;x5yY@smX;w3{qF))C=?Uqa9i{b~+NFrBA@(LVrHO>2YLt?(9;ijTeZcj=%>F4v zw=tG`ZRg!{T`uGX>u$gku(MkFE3XR|f|Gi4UVT<2>YVq5#BI9EQur@L;@l609$q0T z)MaU4_Dq|HmmWM&5**h+pmwb&##R+-SKR z_nMvXO|OlU+M+{uq29+Sk{6E??GUbGObmw?6}r~E5ms;mw%tQc&jrLsYfxDT;%~2q zQ79YU1P`|0ZvHCWofnZ(({Jes0Z%7ZHX*JeP59_*ydS^077bN?>ceZqei^`YW|CRh zV#-qA$WWA9rztayg3v2K&)hZxfW5>{QBm*bi&##4%#DDi`dtkY3;)ZB zet>;FVrj}p`gVpUU9FPO^*PJ1Xy8w{xUNlD3D-erJjPVEQt;06w{gBa)DdDx>188&a(>f1qa>?AP`qaNP`qy<#YDz>}rz zJhhgm0i#DG?A3E+)PgWlI;y`D@tWONlRq9dSkK$UW7e#}g#oO~oV$7V!~vgC%9T zm|;gHaeuEpmASBD2aKsuSwiF-Ra`r{h3)D>QMsM)T4GcBJF=8N3Li=NtJ<%`QuXt2 zu!0kHOB%m(IMxdig1CNe^}?QaSqv7Gp@_`tfQ2%OgJ}BN6sUy(i;qs*G&Yut_VW`l zqHcfc--rw0JAdBpB_+CtI#}kmyOh27*ajVzuE)o^O#(_IU03$lHti+MuFEq)>s%YX zBDF=vQn%@?OTVO`cxM)r{Y>?T`|TIywPUx^$L=cmg%00h8*|7DE9RP>CpPVFW^d?G zn^%J+bRv*0pK0c++T%V<_ya}yDWmbz6J@W~)8m44_p-5oy40j{c@sx@Uv)wh!t65w zJ~NU ztQ%Dud-iXuN>E{RJ&$i>q}s;iPSgNBat^bD;nqe0l^h3Zf`-X^J|Zsw^`JL1i=p4! z4{lQ^>^-6ZUzAl{`135l(XA9g?Qn-4WNkZn?{>&AA4|WcxLHYF%;k!dQ|X`Pc_mft z`82@`=~ZhjU+w*@Y zI#_W^SPU+ilIUm{c;wRDX&^MLe2S`KQBNAzw7SzrAwWDCmi}^R4QWQ*n&YQ~_a5O( za1fud5<2pVV?kZ{ZR(-tJd3%sV)wQGabox{UkCfn6EgwY^>^S&?0J99Tep14Nh>rd z%MCZR?rSZXHpPFf6U&{U6riU*;8-z^xkFd45`^k&o*#VA=pbfQ#rLJSm~Dp3huq; zvM)dH5ZbH_lM^iWf4->3eO#Zs@7p=tn|g$w!cEAy4p(YBAHb6BTJ<|JYJ$6lx3*vc zdhec`8i_U@F5ZsWzrCn%%{^nu*mu3$W}MYF=wLIFkvX8MFeODX@Lt@as2i_L$m#vM z%Bbbt*-q!DJxj$_)F=0neD4KqzK-ABsmU`}n^@0-ZS--aDcT=&%T+!av$f$7S!;GL z|E{8H8j)N67v7jqAb!^ZWF}ER*E@oK7;F>#g5y%_aZ+Rj9 zA3OR|Wp$T98T&7=G5sEE+^OC2oIlNWL%dOF8mFm(0ydL-wcpxPiaB{)nV-&X?)Tumf7|6ly*bsa_sEPkHm6BP{PJ=a@m%h8o0Yow zjQ3e)=C8$9Ou1wFz)TH~k=^XljG-dSy!0;lqF|E;VX{O=!$1e(0ZlQIqF@D|?y&9r zgUxaQ>V^v&XER#F88;GckN+R95^VCJq-<#p>J0{x=<}Th2_<@{!TIG6n)Pq;kr}ef z$uO+gvyhh;gj;2uJ$*U@ZMsIVDsF>PyzGP3+caj5a>UUCrLFX|Ox!ZaO5eZNS#vi03)IprR@vn8 zr^i&104u!u*b_TFldO-bDe7qtZA%_qeKa`l+^=BMmEqa&Lq_xUoP5$ZRh^er4rdxo zbJ%o+EFBda^aTnn1@6Hq%hr;K?0RiArx`WDb2>+`zb_i~!+N-a)$Y>!CRahmLFxV+KTBX@d+4JXS zBj14Bs#93V8f|*;LFo#dS|dT--q^~@oKUH5r@lJO*Y?bj+$&w)(ES5t@jG(;(-HP8 zDt&6A746m4lA}fxo5m37W+NzgPM-{+uEY$X{d-{Zb^>nGhNO=S(I%X8Vr`;-#5aBG zXAK1Zo?&hbO;m$DMfPPewjoOq<}fP98^3n*3RPJ)e=y6}kJWgB!wK z;vIN}zR zFZ@9^FLqzuI9vY~@XKF*Z+*Eg&TqYoj7nuHJp5fh@%M2qO1(UJ%T(fjBlA^JZAL2w4qONbIZk|6pi zTJ#_oJ!16sj(@qX`@Wz1x!y1Dal9YQam+gRUVH7ef9to_KCOw=(Nd+QVxuA=BBE8h zr=&+jL;?~K5${uw5+nwhf@p%^L`TCwnedgEn3#lwgp`yN0)dc`k&%;=L!nR#3JOX} zN-8QUYHDg48X8(!S{MvQPfyRlz;NNh1x7|jCMKqf7cVk1GqbR;u(GmVx^#(+jg6h1 zor8melarH+i;J6^`|{<>SFT*Sdi5#~4-YRd@3m{!`1ttv`T4J3zb+skASfs(BqVg> z#tmU%VG$7#QBl#GH*boGiQT$&OI%zW4u{{qeOp38LQ+!l&Ye3_Qc}{=(lRnKva+&r za&q$W@(KzH2m}I&MBcr7S5Z+>Nl8grSy@FzMO9Vx-o1NjYHIiI-&a>x*U-?=)YR0{ z($dz}*3r??)z#J0)6>`2H!v_TG&D3aGBP$cHZd_VH8n+{P-bRk=H})W78aJ4mR43) z4<0;t`0(MQM~|$ntsg&rY-3|%Yiny~XNN|k?d|Oy92^`S9iKdT;^gGy?Ck8~;^OM+ z>gML=?(Xj4;o<4&>E-3+?d|R3*wc(!C?IT{R0950s{ksf`WpBgF`|>LPJBb zSZr8W*wd#^pFMjP9v&VM5%K)_^T^1^sHmvu=;)Z3nAq6ZxVX6Z`1n8m_~XTk7YPXo ziHV6XU%pIAN=i;nPDx35_3Bk>YHC_qT6%i=>({R{GBPqVGqbX?004ON=1q2Xc1})C zZf@?|w{P?E^78ZZ3knJf3k!>iii(SiOG-*gOH0ej%F4^j-@SWRQBhG@S^56``>Lv{ z>gwv6nwr|$+7BN-)YaA1*Vi{RG&D9gHZ?UhH#fJmw0!*dv9-1J)2C0LKY#x6 z=;-M8@88G9#>U6TCnhE)Cnu+-rlzN-XJ%%8{P;0DJ3BWwH$OlB^XJcng@wh%#igYs z91gd({TXt*zg`e{XMZ-8z>04Bp`e%FnRjHy(a2GXaJ#E_r%zLya63`ivgKi1X<6cpt2Vj*8cp0_`Lq(q-X zL2>nWkRnAcf{WqlsWhvHPJ**w zyi@tb#iw4)%Y<`~HdO6+Ap25!%aVwOFBu84YLg+YWiLWvw+!gv?fU0}Yvh2s9LuJW zcia0Ui#Sf)h5ojNq>P;>HtogqSmrMcNmtO2^}@^l+XwCgR(?7=tESE>;zbMH6lx}a z>V!M_<5*Brci4(29FHDAu(A0nE{`RS8E)%;pg!gqe&+y6{< zGiW}p-$ibAGJunODy;C@ONHrFcftg7Qq5$>uc!{Y0xF!3YVynDvSDsgN%^4KT%}cq z>8(%jrNvcQbba!J%#4}MSKDoMRa-v$@dDJ=ey$3;eAO2&1Ux%Wk@ig{dUUJvAB8%otJF*#=OwO=&qJl{ClzTzTeS;yPv<`v)=8L&IzHO1#6*(*FTahBI_%o zUKx9vj&B6nh;g778-tMY%=g|l9wY=AjR%!!$Orl*HJ;2}^E;|8Zg(j&HWJB__fD!H z?Hy{aLd`Nu?REI@_?E9<(`8Q@Gw2uF%>3T!wQcym6q>M{jN z!pBeiKM8IVYd-&P$_ho$w8(TMs6!8Df(TP$nIXvTMG+U3e+$0E142A?c~-?Q7sDIt z?x)mRHHf&hIPV=B{@!I@-1e`nq{okM#Fm3Pzb>-fzuSD5T{{wZ zvrk)fjSZ>(ruGj=PJi07kdA* zWQf(%1ddIHoCAykMSs5Wd>>wL9x#rJ7=g=8d&5p=0I@nW_d=6zh2P}aOHk=iq$gJLdm6^RL zV$4$7$Kd(?V@&IWX@4L~qesW@1;?(dbwyUfiHr4(b$ThaDOQQK##4XP8XJ6_C^x>; za!Kgp2hPo|vbK`SVzIK4yP~zuOw{go<(n*QjC0CLWbcC!9-d7B2E!G2Jj3#Ggk&jf zQtn`NE<(PgLR2O<;K#>V$rzmE&qAk$sh?u$(1yd;R*N(8sc%FB<^?1gYMp|O z;{VDNk4V^<*JZ$;MZXC7`xb;2f>ckq1JrfV!wE4>Nti|<{TD~w^=SwPdauLSnAR1* zFGYxR3MG_1HsE8MbNvJ+;6Xz;sm}Ly9-u=|+Q{82Xh%a|a=2A03nE&U2Q5S~{+~$Z z3l;qGcs$T+I4cDj5tKgnK`;uneW3oAPySDexz~e;RzA8=RUz029WVKB&Txp}C0Ku+ z`zVP~*T=H}o`1ai$DI%fK}8Tjf6k$VwD@;&ppkAAA^7?>B3kM1cgi5+f@b4&(q{g}a;|Cc^6z&p|eRf2n&<+uOYQXT9E8`^5mITIB(gff9 zPi<*{)W9hjfi&(A45GfkParG;EfL&E_!{w7{`(2shS2hbB1Ho|Gc-K!VnyO}@0MRr zIgih80eQb4A-l?rI`hlxA8sEfNkrSpw}d4{I(Sr67HilTJ6enle|5c7>T1RI_h@vT%U`X{7tvq?TLN{X|vNm&yV zPwDsFRO`{C8)!(elazvQI=gGI>y)fJ*W9nrq?;GuBqwb6Iwf_zwdD`Rhlglip1(4} z7ixSIvobL#l(5HxwhA3-q{7SUm}Exrjjn$-Ax#GlT_A+tk=_OU~L)Q^!F+8P$X#zPb?uh67-5OYJ!Tk?Zo(g zA*^KCi4Q?-=@5zM5Aip|n3d;uLYz+nf3*h>se%;a_k|IXW_H-VRc0s{4?8RBtFyMa z#o)<-XfRBT1|%hJ+(No&3VL$mVJ9n)Tlg$WA{0*Q`}(&aRZwS+>ehYI+X0%b75_?Z&)x#7y~se6ZvfGM*=Dco#4BhV4Yk9xkG0kPnJnKh(| z1}o!?*{aloi$t&+2N#Gj8;Cum6?N`hKO@~$G{oucB9Pvq8q(*oB#(0y^uL1UM}{SJ z<{2~KYjRH@xTzhgIS9%n5|}3$&*+*YZ7NAr2^m3^yK;KscxE>`GZ@*ca%UI_j2V(C z^PFVLsxkfZ7VMeAa=--wbqL)i&qzgv(s6yiLch+7`vu+~UW@_M8LEJ?R&9O2j@FF@ z?<_<~$F||yhHUtqzruwRDc1#9dyj&71vm6wD&{Njw^q9nh8IW zOUi`HR$(`XTI>Wj@6ScUPHB-HlvV6Wb>z}5IVzze^&;Nvc>Y+$BPkU=Y4p2uBkzP{ z2J@3*zUn)NZ6ux?kp~CtO~3!sbpp7K_L~Io+?8M-l%YIiNx`dxYT`WqUhLO04SMAl zuQ6g&E?=^)#F+1sfQdV((;);Kf2LqKikM064zU~=a{~ABEE?@vwa7gmcT`0KVF$K*cEh`C@@Qaw z7`lF(Rb&Q?-X1mqEtEcI5==7DjGIS*ggnpJE;$9kwR)J7TjaN4&;FE&--etNBROy} zlp+SaJEZq~ zED5KSj_1%Hy#`|1-h>Ww-yA=`YP1CdN6jBYEk|GbX}C(OB8}0EL8d;xRRu$sy=AfC zbMiQE7dIeQ&aW&q(_}B7$J=VSn0sD;;A$LF=45*Q0`BVkEhGdz-@shNm;sscvDt}N zzY7woS?ToHaxyhA&|`|sV*ZQsf)Wfhv_-)dwF+#@pk;Bh+LzI1q2BvV9oImqT55z{ z?4L<%0TH3XWrOVJ4gns5*f>Ijm_sw1V?H6VK{?{v*J+VQL{;xM5GctKw^k~zV!ZVZ~SJHx#Fb&q5nyg*%MeEkIKF~=W~Mzo%*hL6D7rJFFkFwV4< zw3?*cIx$?oI4rIR@Ttkzw~UXFgB&7G)*pd9Ty~CoLO7YlaAdi?5l#-Y4(AY@{`MLp zX`5$E1K*mBC-3DvU&{Y(L^Y%qE<@E33RBeu)7hLT3?xpxdt6e#>A#!};dZtUtA_j|8xQAmz=pZy*FxHZ3{3)H!(HD5-voyjgPBS+J*%vN? zCnZz8G`MNeTf;03(0X@Y6l(m@?)0CPxtLRhnG_LDwZjBpp%xuxh)MvAB^BrkpZE#Q zRvU%s(s<@eUIzE<2ID;$@pW!JoLP`li%WY;MZM9fHQ!5Lzd!Dx%$s@7jB8T+Y0G-M zDFJvN!J2fj+I*G){L;qH3gWA$CdmSA@9!_Uyw8^nzOt|WYqz70c*s6UC0z!HPHVIy z&{tx5nn8;YCFaMI_mvOhNuLRjW3(++=Y>fRms__;Ep9zTE^gc5-(Q-+9a)MYHq9CL zU5-=lnO)}Im@@MeNVwIsmH;>Wy;zLFb?8d4r*TrA$a>4}2Vc2*;%0)L zSxsp0utzQvs-~>7Z#I9uR#7==U#C8d|E-o&93#5^qs0GiQXp&nzy__`L%VKAGt*~V z6BqL;x2sFm5`g|i*t;MF{FPw?o!XGHmi_wU15p)=)RTx$XEX)e2lJ82?VNvJMM|Ov z3}4jTQ(UX*_%Zi79G_<)Xb7Br2d3m$;M+tj|pxu_)2YL^9F9i55GQ_ z_m?MM9!0(w*GlN6$Y4y%`~RR{}^M! zMdrv6rAgPOlbt+rfn+=~$4}yRGpSQmmGX?~Jtt}lx-6M3(kl(bz=_1?K03_!sG4S| z1^Tp=ed<()dG)Gn69(`Z5d{P%(swKVm=16gf+YH zQKK3TG~+{>y%_Yzu5RVueT&F{@j1-{800&mg-_LhgIr!+l1JBUPegZGh)unG94I^h>E_RO|r z2*xWO7O3+?$b$KG`5_VU-_j#YET3?o{eaN!XlK@aMIDhO0lcq7Mp0a!-%2SPE1 zj;lV!6Q!idua~`j5DE`1t-^@pobGI2A4bFjWd`kGL5o^eVSUR&B*N3ZZ`g1z7S*O9 zIOBTpcs*Lr;1#MJS=`&9A5WV}H7&W)aJzlw{zRUl9@n+9D`RP4H(}A|p3pMA3T=AB zP{%Xw14J;Uzf`JK)w0lI_x#;}?W(ff6oU;8kz{CO;)I9-LohzH*yPAv4xHteWh13p z1wrXFMFNUf)A?Qg0Q{E(X)#t7E8i`-sz$u)vXHgBD&tq`16w&KF?w`xHv>{iX-!vb z>Rn(|@=dY@QT-q8dVJRi8~XyH|9!jKRb)R$HXJ7@Jr>Hf^|Tx@*rv}< z$v9=H<5HVAwLyh`^vmS}Ja}3BEn*2&y-S7Am&raMPeaK)uAMnf{%eH6jmpY(HN4+I zRA|s|Vuvv@U@6gf2G5RmFXtRyCtbAHcuAV^i1VScQC{_f_^nX`5;1Ulxt=i)MSmY?QatIq<3ga(H@aD4QDL;&5;6`7S$d;%lJEbN!^sH`VlNT{>z=tGrCzPxLh_|{_`mhCfuluqlN)VZ%hm? zNg!10QByynJg;xYbu%ZC9T<#!kx?`lT_790C0&sd+}_AeF^6EJ+}KL{ z5rd!0bfz*p8kxRSO;`|}1nC#tYCYGbt(J2>$Uo4r@%>G-NmU2J2hj4%ug{Rzk|p?q z1Sk8uhlvfZ=PVf^4XBEDP(>jeHB02KoRS2F1en>@Hx)261z*oQADQ~Ia{!J$MaMm| zLCnadSh`;!+?7&W4yp0s?L!!#p^4ptQnd|{uRvwgzL{VEh5J9$yn<^T!Bu2JH7_=Q z5b{dg@b=&k7QY{w7gzQO;i2kLKw5wCeYY-!^k?aiPAecQPc-kY&(Y(>`G`k=Of#`j zUTkLWfG)Y`!>J5>0yR>>tK2ZP9AOui&fTu5GW(y8yqkVP$&E;&A2 zrH3$b$#hiU7YoiG=h3k2Y}FfESbPC}I3z=kPhjd0^kDUr_wRW{kFU2ivS}unS~pzb zn_7Ri1O=_-tc!8KUYL^*jyIPR0QLAtE z=1o!#J#K8sh)h2yh2UEjK_7n>M|e>M?qVOj}Q?e#N)ZSn9PucOR*|1swjg3iI$wsVyQ zB=yhP7foo!>l0$i5%2lmzMJZ6=C}^Q2sGuVmTIWy#V_iR;gflH>Fk1p+{%DELE#uqJ8(j8;f;Vsj&sxAPg z;_pat;0Ew&LPa5c3uS?%ArfZ~P9TRCyLSbj)V>ev{Zcl;3Dt~o%iEa20h->x+NnE z;Ss#?q;|8!_@j3N{(zBy%6eKT#_lBWr=}rdMkYFmIwYttUpF-#q0p0g+~eiHdT~h8 z3n6aox+FB{oGnjZOj5t3{)CeVSKYptOUYKyuMR+EKbMe|O)MvVx=7L#BSo7G)!b~= zCG&K&ZDArnCFpTPbr!i>S>EOepV6GSMIESTN}EKI{n1L|=5!>@ zAyECQzp_K&6R!UNbcHjp+ za(Hd208xfA7fl%jPG;%7VF*yYIIE65qJVF&zk4^4^|sGG(xO&|X8e1|OM^V&wN4{G z0bwMT4z|P} zf;pce_z$4)qjROyh2p#~iF=(iDG3$OwRTq>Qn-I|&`o0S?Ojki94dS79++~~vNiPz z5fO9Wc`ty#V7G{x{xz58BJS$nwq3ds9dRB$UZyIOcY<#H%KCZGUEZ}6e3G`sZZ?*L zGWJe2UwIjD)D^eums+Tyhr%Z%bu`{`6fmCZF?PfIx>k_wwTy%ntX((Bw6HO>wXvxn z4HQcdI!ucGW}YZ%H5_l^XZwpmLN(sQNA70g2e;$LKNLE*Q&PH@f}6hW3y&lvt)@I* zI0k0D32Wrqi6yDs7kXb<(q-Q4o>$q8F%=e)aS&Pf)Esazw2LHuZd0MMdKPMAGdomP z$%^%!TaTXUV`>)DX4jg0KE>{_CRDuW1BCWmiQCh$PP9IBEu-*F%qi$Rdec>+N%t^t zdMww=x3GqnL7DHl9Zap$XW@)8CFqsT9A4_=gnJS(*!u6sQU3y_+c%}ZA+I^0+L@h>l)_SbN0ea z?}zlU2^SDt#N!+A&IK{k<77uNH9f($X2}4W!jz3Kxm{DahJa=DrCcWEWSf}~1@^X< zV5C!-TIOGLX+cFhCF4|yDGqy8O0a3HwJ%J{gv;;p3=^Av3W;z*zYs8E-#HBPlo5*y zUaOzd1&?oM5QE)4+NLTtI`)*r^$RlQWP0Ec`;#3Qj^)F{8*iMg1}BdV3dnKU;l~>& zF4P3YsTad@og1yN^h)lyMMMc~5M!4X)2!;)S4FJ6kS*Iovja z>L#F1??vGjS{$YIF)`}+tv-F#yHWm_fmhBf{g;D)(%8bg8^1SNLx8N!UIV@|VPZsW zH_L7+2R_^Rlt|7k#&rRLj5boYd|Yz&y(vx#iW%0-{zmt#A;aZ+g26|FWB?#rN_xg@ zM2E4{8&LD6lAb1#d(O@#ZQwvr^`L1@ZTyYcl*nu3)eyq%N7w!`x^(_RUL}OADLeMT zENQnOwDj$bZ8BW?%!M*e0=n2PIhqiI)^*WPs0SS+L*hl@4J1;5NBRxOg;z8UB5Hcy zjQcvS$LeWY-|pGqN*+m6uN_xl^$?V5de!x`x{VuuXG@FjLRCTI(%~U7Uhj0}K^l|i zV;*Vma_hKI5r-Nh12ae!P_wS)U-gBoG&OL<%2m&U38;n3wreY@JOoF&Qsp%o-4P%e z?cGtF0TEt$R@e9fh^kGweGxHJb^p7KH&ET@Ld|z|juM;`ZN3InH@|7EzQA0yvwRu) z{R(ARA4iS=5yn3OcP~!QA6(XlFfwdGaDD1+Li{wSD$Bg5bZPN~Tt-g#FA(-^3fA)5 zi&708)FiS+khM-&5RAs0>m~#JSE9Y`le-Y*Xwxs=4G-ujiN>X^Vv00L>Tmjfu>Ta6 z$tNpYImG4p@Pecz*=?w3h~;=>`2-;Yc3(`R0wTpRh2wo{eT@zNc_hFlV+JW(0E7}{ zV6o*?m)+UpON^h#CT=`IIDQTry{Q=h7AGoaQL4WMmhg7S?dmp8&j=grHJJP-W;}eu zK*@p5l1R>%&&sFZJ|&I?oB!d&he-euP}+0ykbAi%eWYItqwhlg2S;yRj+4DxenNa~ayqSD#1$eWfMFaIUMQr8^nv z`+#_RRC41G9scJLBV&~Sia}&!U_aBMSLd;WUi8kVXvsl zOKPJ(-#dpdqqHChMR{1gPG>Bj8aJ1EuT1(&iHjFnYLnNt7CRQHVCZDQ_~X?>f5Lm2 zh2|$k4{fwUxURCv@c#zT2uH>`ap@k)mryxZ<`2$mh?ah<)zD0}*GM}U`O(SiUSx8c zz;Eksb5xjM@zxJH?P08-xz`ntoy{#Qi8bDs(-Zlb_E%R0F{d=B_ic5HJn%@i^_HsO zds4Y98M@3O;WEIguE$3+JGagjQZ|mMb*VKR2d0<A$UmjKKYV|mQsfE|fR5=Xc?r0v`^1G$A?IxBZv z#KV2Lp`<80*_P^f-{Xzd4`~L>;-E@a;FnFMWUOtGxC(b*S&nZKvmP20KVSb1U^xs;zD1C0zvRnA~hrf*T~GKe$)5Y#WVa@mks)4d*JAz;#62m(C!k zZ4?uAslOD2^~DC23fKB+H}8E9JfSDD_#Vlw*GQh!#(~#qwPA^jP^Lzir8O$z;rFo; zRCm5dDSUkvV9o?oSNZ?Qe@=?dnb4+vEx5WzK!$okEwcZIJldMFMk2YnkP0z4Y)YC5 zyO?=e3dhyeaz1v%B1v6(WXwL#6t*$oNBAaJVLQ`|v}h;Z{SGz-aGq=KFa~3e2Wt5| z_|21xXzSli=^b=M62d8^$HMjNtE+K<>f>M0>PEJEJBk$evKmlDbrRv?QdvC?6s|m1 z*Smq{VO3I72HU^&H}CK^xZ|MTfH$KR6-{sCI5)y1n14*={gzzBAsi3DKJ+M({eAGy zQ~lra#pGzG@0U-l%Rd1=)h4#@-Lpw8zNgQKv1Jmn?B*-y?P=@US_yZK*`D;a^nops zlc4$7Cs|@Jrcl0O9^|H4>KDWO%V_H{RhYcPdv{h4T>W0-{gJ-3EKI_-1+b{-%_^E| zCYt~FJpw}k+wozq65AxcQd(KJxCAae|)c22)<0hZsM=_%vbeQ^Vg^IB_4){er+ zIRmx4qSgMJk<9qeD^1?i4$+J_JLhTZQM+N1Qel_&xIT=D6!?N{1~cfv_TK$NgbFpv ze=jd`ZN%i)taDgg4<>jA~swxY;qi$$7%=5oqMq<3+Wd z-k;|v6EI}Eb%ECNaUHc@2fC+gmb*JN+_Ix;Cu4)L@pc2v_^T~NIcldF|227O)1^5M z*1_-wl~>&WkH^%2TTx&ZG0?m6o!X_WQ6S?`BY3|6XByc?D2`sViF~mB5oEhS!MB7| zEVmvnLvkLFKA$I*e*anGNWrinUh8M4O^d^=B)N@NL-2*)#+q2=*XI7qVEg-bKy|O2 ze&?zn8M0f6pAT0m2VEO1NaA^IPQ-@c9E{3zkPe?R^!{Q9&iXOOt$h%PBIpy=TzS$t zELY#$XJEVCzD9(j6@YO(&+bQnZ;8-OXV#j0rYSqrp~=s+Zi2JV`Woy#ao0LfQgfZk zKb2UgUqqs3_02yNhiLU-R+Iq0NX^|0rH(qVvWG|MMKw;MFgoPuq` z+o?@1Hhh7h1HxJ1cNBF6_lu@h;hXQRlYA;e+r9Rl%>HIOhvG8vQtb%UCzc(LpwgAL zN41^NW&ON>qBb8%slEE)ZO#lUC2aVj4lQWhep2VK{b(*uyF_NWhOO-N4?X6ebm#>8 zfcZ-(sf12q(6!c)ov&J(Yd1!$-sZ!^tkf4ePge=3CHKI7vzI~b=cKWwIJYGGXN2^E z&pcMSVtVOLgh?pfnV2wcbWP0qOHhM+B9>0}4Xfu58|O96QtqCy`ytH6#!68UkS z;Yv$y=I`TYKRESEWMXU5%3pt+b25`%=nJ2!EEW9~pD8Vr!2#h6lGm+u953qFL{C

e}*id*>$}#Klg~Id4vd>a~qiEgw{$6PpU6;^%e@;fjy!C(}rH0%Mo=lSu^v4a*>zS9-_asM2l=(oWTEouuOtf02Fa%e7hU*|N_LaT&LWD0P-3F3FxALj5(^c}nW`35#mT;JZf+FkXWu9+?Fx8EiV~S%*k=tH? zRA;o;{_(e-mXTuP%{&PU3?0`FZcTYL{&;~6zQnp@zzdaZX~LzhsV#rrWOwyIp9QN0}LK{YQfLXy_RZzM7AT9E{f6lsoHhH(ermyIPkj_SKgKYd-%TAz8< z;<|y4T>czRgG(GQRg9=<&2*$keQYxA3F-?ofkGPgKZ!?0?&;@sPl#G1)f$@FbF3l;mJK9CvNNI zd3W{St&)U3N|p0Imw%eB{}TVJYxtjW|DOo<8LD0=&>m+S4Q{fV=oT+qrJEFY?g*44 zgs+imkFKOflUr5<|H11LQ#}0hu99CAmv8GSGU0}HIYAu5v75LC9q<@OVtiRqir=|1Z{RxHgtc)X~51XZiFH##5P+&Rst2(oG%hXG=lRM>W3!3t7gjv&dj z|78G`W-lDZ60{EMM_1e(^2d>tr=@?gb{NyQWk-$asGFOnzDASHxrY zk{ERoV?x0%v>-4qbdri7a8s5qe7KS_ztweogR;8AInxx*OM49#=%oNv*6tXl$ zhU*B5`aD`jNrm$B23ZvF2eE-;m7Nv0bT3p_z6z}g|6b}XpSCO6O0a72vkk`eT$#4* z^5G6t^gdrl7ewv*x&2&^YNZdX#+E!E+8 zR8t;a^X5!=2O1bKrF=LOYy|AJoE=3CE1XD7tLJ&rdp2GmK8YUnI>`tY7~CY$2HSsk0|9;ihv1C+=1VK3X1+m3!cX6b)ReW9%I-eE{y*U3R~G;P diff --git a/md5sum.txt b/md5sum.txt index f8f192bf..24528de1 100644 --- a/md5sum.txt +++ b/md5sum.txt @@ -4,10 +4,10 @@ "config.yaml" "ba8182fefb4b3a156058872157fb2d7d" "site/built/config.yaml" "2023-10-31" "index.md" "d45b8843b661f7f6d628b6b57447b6c5" "site/built/index.md" "2023-10-31" "links.md" "8184cf4149eafbf03ce8da8ff0778c14" "site/built/links.md" "2023-10-31" -"episodes/simulating-transmission.Rmd" "9c915400856003c134ca46edd1c52d87" "site/built/simulating-transmission.md" "2023-10-31" -"episodes/model-choices.Rmd" "113c11e0dccc812744f685df84733fdf" "site/built/model-choices.md" "2023-10-31" -"episodes/modelling-interventions.Rmd" "01eff78117bf3f53015154ca9b961651" "site/built/modelling-interventions.md" "2023-10-31" -"episodes/compare-interventions.Rmd" "d594dc0466f735d14ccbcc7b2a9d56ef" "site/built/compare-interventions.md" "2023-10-31" +"episodes/simulating-transmission.Rmd" "ffefab0380c8477686b43d957ab45f36" "site/built/simulating-transmission.md" "2023-11-07" +"episodes/model-choices.Rmd" "2bce847d6462de90d1cb479698dba709" "site/built/model-choices.md" "2023-11-07" +"episodes/modelling-interventions.Rmd" "942c73e444c4ec7080a1727491ff4f0c" "site/built/modelling-interventions.md" "2023-11-07" +"episodes/compare-interventions.Rmd" "b22d68654bc38c33ac8e4916e1be4c5f" "site/built/compare-interventions.md" "2023-11-07" "instructors/instructor-notes.md" "ca3834a1b0f9e70c4702aa7a367a6bb5" "site/built/instructor-notes.md" "2023-10-31" "learners/reference.md" "451e0fd931938d38cf92a313192c1cdb" "site/built/reference.md" "2023-10-31" "learners/setup.md" "b3c6bfa13fd687f926bb1a3e772a2516" "site/built/setup.md" "2023-10-31" diff --git a/model-choices.md b/model-choices.md index 67bf2a5b..de85691c 100644 --- a/model-choices.md +++ b/model-choices.md @@ -62,7 +62,7 @@ Finally, interventions such as vaccination may be of interest. A model may or ma ### Available models The R package `epidemics` contains functions to run existing models. -For details on the models that are available, see the package [vignettes](https://epiverse-trace.github.io/epidemics/articles). To learn how to run the models in R, read the documentation using `?epidemics::epidemic_ebola`. Remember to use the 'Check model equation' questions to help your understanding of an existing model. +For details on the models that are available, see the package [vignettes](https://epiverse-trace.github.io/epidemics/articles). To learn how to run the models in R, read the documentation using `?epidemics::model_ebola_cpp`. Remember to use the 'Check model equation' questions to help your understanding of an existing model. ::::::::::::::::::::::::::::::::::::: checklist ### Check model equations @@ -85,9 +85,9 @@ You have been asked to explore the variation in numbers of infected individuals Which of the following models would be an appropriate choice for this task: -+ `epidemic_default` ++ `model_default_cpp()` -+ `epidemic_ebola` ++ `model_ebola_cpp()` ::::::::::::::::: hint @@ -118,7 +118,7 @@ Consider the following questions: + What is the outcome of interest? **Number of infections** + Will any interventions be modelled? **No** -#### `epidemic_default` +#### `model_default_cpp()` A deterministic SEIR model with age specific direct transmission. @@ -128,7 +128,7 @@ A deterministic SEIR model with age specific direct transmission. The model is capable of predicting an Ebola type outbreak, but as the model is deterministic, we are not able to explore stochastic variation in the early stages of the outbreak. -#### `epidemic_ebola` +#### `model_ebola_cpp()` A stochastic SEIHFR (Susceptible, Exposed, Infectious, Hospitalised, Funeral, Removed) model that was developed specifically for infection with Ebola. diff --git a/modelling-interventions.md b/modelling-interventions.md index 3808a76a..178fb072 100644 --- a/modelling-interventions.md +++ b/modelling-interventions.md @@ -54,7 +54,7 @@ In mathematical modelling, we must make assumptions about how NPIs will affect t We want to investigate the effect of school closures on reducing the number of individuals infectious with COVID-19 through time. We assume that a school closure will reduce the frequency of contacts within and between different age groups. -Using an SEIR model (`epidemic_default()` in the R package `{epidemics}`) we set $R_0 = 2.7$, preinfectious period $= 4$ and the infectious_period $= 5.5$ (parameters adapted from [Davies et al. (2020)](https://doi.org/10.1016/S2468-2667(20)30133-X)). We load a contact matrix with age bins 0-18, 18-65, 65 years and older using `{socialmixr}` and assume that one in every 1 million in each age group is infectious at the start of the epidemic. +Using an SEIR model (`model_default_cpp()` in the R package `{epidemics}`) we set $R_0 = 2.7$, preinfectious period $= 4$ and the infectious_period $= 5.5$ (parameters adapted from [Davies et al. (2020)](https://doi.org/10.1016/S2468-2667(20)30133-X)). We load a contact matrix with age bins 0-18, 18-65, 65 years and older using `{socialmixr}` and assume that one in every 1 million in each age group is infectious at the start of the epidemic. We will assume that school closures will reduce the contacts between school aged children (aged 0-15) by 0.5, and will cause a small reduction (0.01) in the contacts between adults (aged 15 and over). @@ -114,7 +114,7 @@ To run the model with an intervention we set ` intervention = list(contacts = cl ```r -output_school <- epidemic_default_cpp( +output_school <- model_default_cpp( population = uk_population, infection = covid, intervention = list(contacts = close_schools), @@ -150,7 +150,7 @@ To implement this intervention on the parameter $\beta$, we specify `interventio ```r -output_masks <- epidemic_default_cpp( +output_masks <- model_default_cpp( population = uk_population, infection = covid, intervention = list(beta = mask_mandate), @@ -200,7 +200,7 @@ We pass our vaccination object using `vaccination = vaccinate`: ```r -output_vaccinate <- epidemic_default_cpp( +output_vaccinate <- model_default_cpp( population = uk_population, infection = covid, vaccination = vaccinate, diff --git a/simulating-transmission.md b/simulating-transmission.md index b3f824d2..93ff7980 100644 --- a/simulating-transmission.md +++ b/simulating-transmission.md @@ -64,7 +64,7 @@ pak::pak("epiverse-trace/epidemics") ## Model structures To generate predictions of infectious disease trajectories, we must first select a mathematical model to use. -There is a library of models to choose from in `epidemics`. Models are prefixed with epidemic and suffixed by the infection name. In this tutorial, we will use the default epidemic model, `epidemic_default()` which is described in the next section. +There is a library of models to choose from in `epidemics`. Models are prefixed with `model` and suffixed by the name of infection (e.g. ebola) or a different identifier (e.g. default), and whether the model has a R or C++ code base. In this tutorial, we will use the default epidemic model, `model_default_cpp()` which is described in the next section. ::::::::::::::::::::::::::::::::::::: callout @@ -313,7 +313,7 @@ No we are ready to run our model. We will specify `time_end=600` to run the mode ```r -output <- epidemic_default_cpp( +output <- model_default_cpp( population = uk_population, infection = influenza, time_end = 600 @@ -402,7 +402,7 @@ output_samples <- Map( ) # run an epidemic model using `epidemic()` - output <- epidemic_default_cpp( + output <- model_default_cpp( population = uk_population, infection = influenza, time_end = 600, increment = 1.0