From 564780fa41383bd7f04ba58a8da63e5ec3094314 Mon Sep 17 00:00:00 2001 From: Vijay Ivaturi Date: Mon, 11 Dec 2017 23:23:38 -0500 Subject: [PATCH 1/4] parallel-first order absorption --- absorption/parallel-firstorder.Rmd | 63 ++++++++++++++++++ absorption/parallel-firstorder.md | 57 ++++++++++++++++ .../unnamed-chunk-5-1.png | Bin 0 -> 39123 bytes topics.Rproj | 13 ++++ 4 files changed, 133 insertions(+) create mode 100644 absorption/parallel-firstorder.Rmd create mode 100644 absorption/parallel-firstorder.md create mode 100644 absorption/parallel-firstorder_files/figure-markdown_github/unnamed-chunk-5-1.png create mode 100644 topics.Rproj diff --git a/absorption/parallel-firstorder.Rmd b/absorption/parallel-firstorder.Rmd new file mode 100644 index 0000000..a8c9c34 --- /dev/null +++ b/absorption/parallel-firstorder.Rmd @@ -0,0 +1,63 @@ +--- +title: "Parallel-firstorder absorption" +output: + github_document +--- + +```{r,echo = FALSE, message = FALSE} +knitr::opts_chunk$set(comment='.', + message = FALSE, warning = FALSE) +``` + +```{r} +library(mrgsolve) +library(tidyverse) +``` + +## Two parallel first order absorption models + +```{r} +code <- ' +$PARAM CL = 5, VC = 50, KA1 = 0.8, KA2=0.6 +GUT1BIOAV = 0.5, GUT1LAGT=0, GUT2BIOAV = 1, GUT2LAGT=5 +BIOAV = 1, LAGT = 0 + +$CMT GUT1 GUT2 CENT + +$MAIN +F_GUT1 = GUT1BIOAV; +ALAG_GUT1 = GUT1LAGT; + +F_GUT2 = 1- GUT1BIOAV; +ALAG_GUT2 = GUT2LAGT; + +F_CENT = BIOAV; +ALAG_CENT = LAGT; + +$ODE + +dxdt_GUT1 = -KA1*GUT1; +dxdt_GUT2 = -KA2*GUT2; +dxdt_CENT = KA1*GUT1 + KA2*GUT2 - (CL/VC)*CENT; + +$TABLE +double CP = CENT/VC; + +$CAPTURE CP +' + +mod <- mcode_cache("parabs", code) +``` + +```{r} +mod <- update(mod, end=72) +``` + +### two-parallel first-order events + +```{r} +ev1 <- ev(ID=1, amt = 10, cmt=1) + ev(ID=1, amt = 10, cmt=2) +out <- mod %>% Req(GUT1, GUT2, CP) %>% mrgsim(ev1) +plot(out) +``` + diff --git a/absorption/parallel-firstorder.md b/absorption/parallel-firstorder.md new file mode 100644 index 0000000..089214b --- /dev/null +++ b/absorption/parallel-firstorder.md @@ -0,0 +1,57 @@ +Parallel-firstorder absorption +================ + +``` r +library(mrgsolve) +library(tidyverse) +``` + +Two parallel first order absorption models +------------------------------------------ + +``` r +code <- ' +$PARAM CL = 5, VC = 50, KA1 = 0.8, KA2=0.6 +GUT1BIOAV = 0.5, GUT1LAGT=0, GUT2BIOAV = 1, GUT2LAGT=5 +BIOAV = 1, LAGT = 0 + +$CMT GUT1 GUT2 CENT + +$MAIN +F_GUT1 = GUT1BIOAV; +ALAG_GUT1 = GUT1LAGT; + +F_GUT2 = 1- GUT1BIOAV; +ALAG_GUT2 = GUT2LAGT; + +F_CENT = BIOAV; +ALAG_CENT = LAGT; + +$ODE + +dxdt_GUT1 = -KA1*GUT1; +dxdt_GUT2 = -KA2*GUT2; +dxdt_CENT = KA1*GUT1 + KA2*GUT2 - (CL/VC)*CENT; + +$TABLE +double CP = CENT/VC; + +$CAPTURE CP +' + +mod <- mcode_cache("parabs", code) +``` + +``` r +mod <- update(mod, end=72) +``` + +### two-parallel first-order events + +``` r +ev1 <- ev(ID=1, amt = 10, cmt=1) + ev(ID=1, amt = 10, cmt=2) +out <- mod %>% Req(GUT1, GUT2, CP) %>% mrgsim(ev1) +plot(out) +``` + +![](parallel-firstorder_files/figure-markdown_github/unnamed-chunk-5-1.png) diff --git a/absorption/parallel-firstorder_files/figure-markdown_github/unnamed-chunk-5-1.png b/absorption/parallel-firstorder_files/figure-markdown_github/unnamed-chunk-5-1.png new file mode 100644 index 0000000000000000000000000000000000000000..59a91fa19211cf3885baa0023e88494e676e3158 GIT binary patch literal 39123 zcmeFZWmJ@18#c@g-5t^(B`HXW#L(T{B`pZ3bPOGm(mixHC?y>NN(oZZASvDXUUP@f z{r-FZeSf~SuH_nM?b+Af=RWf|&P{}xiYzt;IR*j(0=B%Ilm-F<7>0lVT1G;?~kMojtD+$y~|aVR`cQB$sp?O_ynb8Im4Y7ZQx46V4Rzz zWf5(-3{u;AU$mR{;#rzqOYNf7zRyY4LALjPv&{kWR&$wNay$`;B)y7*lxOu#SZufm zTU>A$j1(nTSDlK@xKs(uuZrp3Xt*TV?QJ?cwQ!$P`ALR~nAU~HqW)J2W_T>nM$-9h2J@e=uX2Qs3TzH(cFmNR zCG=%>x%GbHGdMnefz%n%Vwiq3;CkbznSH#uvBQ@>?CmLad)?+Kckb(&^ZfU=&(=wB zG6(8yVM<~Uac@w#f3A%9`R3|ywi?&7;u}B3bX%tCz zuLCq-4^s8-U7^+bDCO!=(W-MC%{#kS>U|@mJ~=)i{N(}u?b!9`$SkQ(!#-mu>mlO5 zWPzY8k<8eZ2H~La5LFUAiZso&CZ*Wy6h0PeNc+~+-ko^=<%?<1RntfV>&r0FK$1H0 z+n5TH?^Z|UAx4&jVt(ubU@@jgv;?z3qb5A^=kYE)AbQV8Vde7#$~wZA2iG$SucS?c z1%BMMMLp9f%8*`o-p~p&dw9mGWTZJ=!H)AW^zrbM8)W0pxBaHIM%(NH>K|7hWwq#z zT$x_&-twHY-P#u{`>*RBeR3Tz_}W^h)U)?W?~S`Qxg*LN$>DffUqm0(iObQtF0PDU z@?hF#^{KJ2{%cMgnTizTYx^smM=H@$6{kU`VGZn`nSELA`agg39rFF=JL21A432*> zoVC`z!s$%&#ZOif@eCIdwo0H?`nX<$G;yt!?MTCsaDGl?B_y38NS!Uq&VJYh5(y zp5)RescvZ4;c33mBrL6+6vqv#FLHx9f^TlArb?aaG zkXQ4wCa3Rgg_n?peo>0C)eM_+G@w8}SUinA&ERRpB{bxe&Ah?tY!^G)yWw1oX+TJ( zga!;eb|M=cLR~}_JPDqwoA08Em-E3~p1jyQ@3`ji@VlGb{@MGcVv_^fDDWbyW9i`S z4^hZNIf6~JKM0oNqpfCk{~7@Sf*>y?uH^yR%RN^`)s=KHg1`juPa;D%TM}I&g z??fcVt0t2VN0UcM&`S*)4T|u; zrvwQGU6nup;r)4>lMsm1bY1d!{`VCEhy~V!@_)_6h=4q){j&t?f73${Fj%VoGY#;p@qn>nhUo}LY6U4bRjc8!6S12N;+0zY9x?)%M~ z6SFos%BFr-g>Eesbq|+ov2Rx+m?ub0edsDH3JVKQcTFGCRmASzo*k^d=PKISvDaQa zK6X9xy6off6aPK`sd(P0HM)r_4p1h76a=EN-rwF0lt`4!YTDY;miyLWG>$)ztddE}naJwni;B)0_5#qGa@~*6#*KUGAn^)h;#^LwYcz)2~Aj@vz6HP75 z-ya?4ZNC>O_qx7$cyoP{^Pu#5t}#EFLSQ1D&msI~B?wm(&7QCwuq_4>z+57xmTWoD zq>?(zw9B>1rw>eVpC20IE=|oP)B-6jP(PdOHGBP&2UJ_k( zDYMD;EE4Kn8D!2W0K{0HthWCE%;w;@>>!owAVG*hL5l#cB9L=w;2|)-xGw+PJ4k z9uxS6JzzqA+_yBSzxwPm77kU#?s$EV2bm|!b-A>$W?#|>IAxY~9wAX@`}D4fPFES1 ztH^v%*3seAA8fnXZ!PtIxMx+(-s7_wp=>$q#+*73mVd&=rnr!Nl&*%{;}72)J`}t| zZm)#|K2n7){NAI+QKDoDAGjO*u-~%F>9t$=Dcq*b{b13UP?3LjcxY%c!bo7rr=In3 zYg3NruSW`siuu$cUUGfSj`JVOl9HLpBPR2tF{e02dngP(IVt;49lB(CyclET8tGz2 zgF;sEI?2X*w;dXYkn&WZm8MUgJUKlc<|;RLn*sONLS~8CBv-j}s@sk+?v@_za)dsht=c+P@IHUBNBG%3u-~gZ*qB z#5I&SgSE!c^y}u8Ak(-C+;z4(3qK{V8+HLHiB)8KN89gmB;{^nKd)30jj3 zbjiVUiekm}U8Myo>+y#@o}r zQj%r;HRSU-F7o|*3gnyfz@K$rnykN(ubS{3-mtxG;9%IQ2-kwy}pa_PHOtucD&XAdoTfB{#-w|_3VeC6~Ulda{Z|2{4*|!i=S^|$rYz- z#6P?o$q~z;_FTu8FK|)8Vu6j+O3D9`De!>=a4wAq=wZN7{dojs5im;Lu&WdN?<*)2 z8IUyxq?_))KOTYG5T#-w(*KT$)+_)x0L(aNwEzD24bq3ir=kOyY$}T%%F4=Ye^!3d zi@a%cKQNhV_ocWyPO7U7W0&9YC5dJdU}lb5 zItB_N-e(_Vems_q8S603Pz2%+cf8&9MB(YyV_Z$Z``XW<=k+;#21MVl6%}H}E|#i`?(c3gsM&fzH2$~o^t%WOs;VjCl-i-+RbE6h zEnIA-e$g~GR(6l-$Pyf>(Vy$^7js-}h`2ue?ci2IbcIRc0|`RKE-GL;p8u#*bsKto z@?f6RV_I@ILs=jc>Ij9kjh_nx#}}u-;$br-cD_L$4X!C^Mu9*YA95 zB4xG}hq)c${V?F9=`4!Y8`W7ntgF6TRXl-+D~wl*^#ng>yNmG68`mxO6s%L&I^H8B zAEb8UdhQ|?4GoPF$@w2|uV`QEEm>e-CMn9p+S{{bq7)-bUo}6U^&shMKviot28%iU zlyO_wJf0u?tQyf1df2f6dq}|I&M+}d52=~uYPTcTxh?xyUY(VKeff18O66hhy^t2H97M}?)+t+ zi-e6WkGqd0Q%Q@5WQ(ox@3?O$i{JF7sCL}Oq`!S0CVGEM@lN-90!rzUr)i0ppBgGt z6<5c>Yum13Q{?@Ael$7N<)s5ZsT31;+(w~0jWcR{uVXWz?OXNJ5&smFLBHiO{4+m~ z9h#HC^Y>f=51!^%ivBIiXzE!Y9r%iS9j-`ow ze{W=5YY<6}i}dnf$7ms_V*T-l*=nA$HCOD>z z9=U1W3$AJibu4i;oG=MrgNfZ3gT?na;1mh(?4q()f2H zBE*$%=OsUBs(fdEN3(>xR6Q>>`qeAp$Ez8I2VX1uvD(X3Y9K!B*Cy8$%_n~9ydd?J zHyzDUVa^%TPl*4-CDK9jPV9m;y{yfy{d6nKAbTo$$<(XeZL6%2LL2#XPgJbfTyk{k z{Cza3iEr5D%3Q6lLoH?w?J_zzhuu(qTIG`M?rgvI?es&`*h|q`LlE}!VlJzK;wQ;p zbhl%5jmv`R6}}%=zRH|xB|XZK!xedvf$^dtKVBD$pmlQeF{4$l!XcFaR&j(ttj0x_ zrr}@T7mNBOBZxqB3%%SGQQ!JYn^J*h zid+jC&wb(+>Drh@m~cUB&YL7WP!27v2uSBv=%Yg=tG`#bN}Cnh<6Gn-N zr`OF+lJ`|e*ueJ?#6U|@_w7j&->!9U!6#2XKI=GJ?-#{GOX-}a;Vcv_&MMedZ*p1R z&?=8YcARc0w5u)7@(J|%+V9v3gQibJj)2r4`8`vFg;nznSwll?+e#bXxV11EEJB6) zBkq>k1tfhX7xWvZBUaT=fT(VE(x4UB|{!V5yn9;H)>O@Hi(CmJeHh3X4fc+^6`14C`w+)#X&B?;12{W9Oblr z@vS>`n1KxIE?-;E_rCAbjl5c7s*b04y2TLI@$=?#SgYXz;q%#Q4^B$dDuk9F=_Bh} z+0`OoBZrCo0-|~q(I5`(Y)4bS*5t+8WzFkvarLDz$@v*xltztHUX0G9YzV=Xk?Tg1 zws_sKQn$3J4@JZ4n>x>-OQogbGzo5OA;^{MS2w-6P{!34XdXcbqYzxBlSM-%fkjf` zX7e3if$zVk?>G-@N*u**_bAhkCpxx|ykBco)~lDJAqPRg@H~U;)888uV7bY|p#pHUIS|2Y$5!lY!~ZU#Nv>sz)myh)&N$8a%3?S=rl2 z<@L?jz)#0Iw8fa1US%tOg2ym<0zFn#RkL_hn_2oO%7XzfoD{tG#8MUPO(ns2E_}$n z7ei|ITlP)EMo@;8Kw(3JR(!}P837bbE_qp7b=_Nq+^g8_+j9G&#i@yFpLUo5#pdMv z+_=k<=egAPA(}v%k3N*ed7n|IMyrsoc9`!G1Egp>ElWHtkPO8$Sf4pO-}vNMcz?tn z=8LX$OEBQSM?J0B*5Oe6WuGnvQOL1>_jL8KJ1;q5nJ(#h0CXDd_O@WZcMakh^#=A~ z@b@wucQN+o7(YJ^f$OcDfy!*pZ*P(g1H0;%KmXYGNW8dYlMw!hPN+r1g_|_R`pgv9F`?b;g5v!e&8DD(2}fr(&2dPz1$5Bb&(6Muw;&ZeL70^REFSN z%+`%1qr)eadW|tCISXMP$sENFp~h16i+Yg9{UTECMq)-0>N5#3EL9iNzx*p=D16sK zQ1e;2Fi=6WXwI1`?fcbaFCwn7Thpq@gdtg;?b;DsO_#&3r_#2B$8!OHKQ5KSb0i~V zb6GiLZ|Od9*3lt_{k_3%-pT8>&|lr5Pow9(L~z2)@5|2b^^>{0j8+@G#u0aqQZL)1 ze=&4-w$xyQ#~GYwL+{jlhmizOVqSD?4$3+Ct22{GB$X%lO|CIV%1m4euYR!*__suKV(da zjq7bxo0IzuLRetfxy!9K=(lJ4VyqXMLcXW70;jRlOS-y0*GMaN{oJ`ld^XRTeuzxe z%SUER(YzXXOd!khE0UOPtADpzJvqsTuA6{hjn@zFxGt?CuBafjBj=EnK$uf2%2yFf zq1IvJbL?lOW@M4=Mel>?2)qL zr{=S5JtvUJVJ9NKqFI%>lo7AWDLp8gYDOfgvrVu2>-*1GY5B^qHv{L%A6gV0&ssPJRV;lvJGegF-}U;ggfzU_RI#5P znM-vQn9Gj=Mdu!mbj*JCI;%TnFRQ31`0#(v}~Va`O3zPjSXJETaVu__yPQo{&Cz>Hah80RGYtF$v_^kne82;Q)wW* zkjQA$alOOVay-N~WraErO_5^DGG>AoP7Y+H>XPpR#s5OM+hoiJVyF`1<8xVa*u8T^ zeQg1FZ?f7{?4AB1on)v-rF?sS*)UEXsBrFf7hrZ(z>?f zdtq#WxpjMW(s&1?&v@aPz-GS8t|rvA9za`CJDmmpbF@$q@B|AC*mOPsAO>7&761TT zA^_`BdjX1S2WX>(fG8jv2|Pic%cnqTt~^{>UY`Hxne|vOI^ovc#gZFvqrcQ{(aciOj^<6q@8nk|x1)v{AolFwRBKj-!| zZLT50OS=RCiITDN**idr@2*Gv?B&g|$_OhdD3C&N!Ll;j zr8`NA;Hd7N0NjFPvozBHr%|m1vWrm<35P*R>+AL_>+50Ddx`JGs-$7azJ3pP7iEpS zwlA9pgF9P zVkS>3&`9Zp{e_XldH{ad&RegCnQFlAc@AYh?3@xLPyo4{3ey~S#4RO9*2OzK<)cf zj-T7x+v8>3wXU<|qqyannUn)V4Tfq42I(90diqwFC2FEqy$;o%llPm}V&bfvv83{E zZw@+gJX}MQcB7L4D9bnRbze(<{FZsl-c)v$M0IbwhE39kBo{ zk$eq7@w@_!8Ue7T6`t(p;%KCmu$}e61a0IT%Yc#0=ij%S7usa0J%ihCj};e$&(Dv% zj&(wlYL5r~bz-K$op_IVq8v4(!ro)!X#^<%LkV z%+2MM?ZpKTf5$lruL-~j2%UU=2FI*-2yAZ$?Plw&c5t0oHWFf=uSGq(loU9U^Z?#n zcAW1G=^QCw$eBpXOs}G(iKAO#661YdtFalzsVzoj zjYL<=lLQ9#9Gy90f*8d|6<=!KeZ#EdA<&WWu*ywLoH=#uaau`HV8b9VgfJ_o7qGrx zsl(YYH`AYa|J0%vZD@4#HH$(1K}E^&YwP)sqKtLpAf^{Z{t8+~J1=a1XJaY8ODD2R z{kAsRY^-xhPO>iX3t*eHkt6*ntWr@utCq?7kD%B>eUrMVPbS6{bICe!c z#XkIOuuHgpd;{waC5b#VqVD5X?W5+3_A4FCW0xW;f;#tPbyrL&B=UO-039Ms_CM4 zY40K6fK-C5ncByVP74M}0bw7%6#3Mdb%TdUC>ZeKvJg~vp1IZz(?5L3@%a!X61$aR z>@&~yX^5T*aCQ+kUBhF~oq>NaVV0Z*N_Ga*(MDsM$gDb}56!T}Jy}Gv<@C^ltW`RQEKp~kcC8h9l3~i~t_<=S zx4g|2^}Y760MEoL(A*nY#VET!P!=70d~-?Hi+FY#@<_H@62!Nqr1uhQnHCd?lAe|w z=CSeyjU}>H86Xmo;bb1)x~aK&E_tNm;0MMAfVAYhr&B#%)63;4e(ZQ0Qca3=u zIdR|0ZGQ|E+h+=sre(v+e5|{?`aGzg?CCch7Fe_p@e2zH$eg?!ckVb$L+>RyLYFRH zI2gn|i+f>6w4C>)2K|uXaZFQQv5%`cK!(({m6^;Wz3>QBgS@$;%kw6^#(2<4LytuQ z^WZ_j_G6eeQ1VW35;Q?6NIjH(FHKTPJ&Fn}Jd*JJXY-qP?^hljwjqR%=drXfw=^YK z?&z?v`SQ!t4~-mZ`?Q#X?#zqd2f|qKp3(xD1&ERbhTA=oiFOUfO1Qnpu)tiNmSECN{ znvlCSG5z~ftvt78Dess|>zQg^L<m)XAtKfrd8FYCFb2;hu9Iz-W32Y#zz~=ZTF8F;jNfGUGIZ0I2Rpl*XBgY7b_i%LE0ut zyZ0^{I2go9588k2k_Bc&X>v^&2*@EBI*!B;O5YFN7fP%NA9I$6PrHToXC1Ughek$J zKM!%#rMV*%I$4kipV<}Z$z&7}Cg~EQTk)yu=PuO3+i9R;72N+U1x*9OjXy|N#=@kI zj#gM_j&}{AzWV07};*g5qUTJ{C0}&JlLM0$g%hBwx-{zpcMtTt33bOuk z9WMi*l1O0FQ#|i!x)K&9b97shF}(WaM^QDkt8Y22a) z)M+fwzIV2ft?9g%my@xC7;fS*vi&N)vyy6a^C9OAAiFTI(1{%$fi3-UgusR(^$+3V z|A@w}R|A!ap_GR=dUSV73=IzEGAbsIC36^-R8_@Te@mA}M8T@jmTgd}DmzQcOe1dz zWr{o=%uRKpW?0pOWf_Kul02uyLA^uM%A z1N_2m57yQuY`l>$Ovz_YL?!G2r?DSb(wA*cT^`Uq`_w{$NkkZ*H%G%F%nPVS?kpbr zzYzuSH(`)GGFh}kK~)u?DP@Y|_4VoYTP05ii;K3Xh=I4R3X12xq1QCIsS@#BexJ}% z$+i~=BuBi=KuPtcwAw!Sof?@(Q#nUt>9h4V-&46QZM%cfcgV(m{xo8J@+7y?@RNpu zLS$Jc6CN9;9$@Knf27t|vWz^aREYc+?kC#CS+EL zBa%J@9u@zxU0V{hatxNhaps*uf)lq4qde{ut|D+JGL7bKiLf{t5!6G-;Y{TQ740jme^*Z?^fB;Vi3Je-elVE*p4Eo{~70}ZBP7T*-7H=y-8%iBU1m- zH8d?qEm1LG08x(MSD-9jI|6MTozTzcH$XKAV-@g{xP^C%oqQ9}1nBhBx*(v$It5D8 z$+D!zMjqUx8x>VWqxS4%IoYrP3Sie9o|7k zPPBU2t}|=k;QKZP11R%ih`I{vl)8B+GA-{i2A^8(BKNJu2dtQ{*818KqJ4n?@2u=& zV0hGm#6yNR^U*Owh<|c9F4amu3&e#4n5vs4QyU~ul@y*{L9yW};9q_`2E9I7BRjQg zS8bdH1!#%ubm%*A&E#PMJi)6+DsGJ~F0i#6PU5Fj#X1iw83Rfe%;6JHT0H+7b-{D9ef@Qw)R7vN98pIsb zD1kX3sNA{g{o(tR0at7fPTR$zH?P~C-Y)a=w+uvZZ<3lycnP!wmD+tT+7bG}Oht<-@dP(Bz)tG^zDW87;EK3pqYwg>50I!a>Wcu#i zJ6?w$?0kB@B`21RVn5~OJe4D647{2GOX@;do^j}W$X}FVpp6C(aW`{*KGkT%y;Rm@ z@WWgHfDLQVLX%7W7yHxa-si7YL91R79Hn8nG&&Eb$A}kHB4t05CWCj9nAcW%hJ>sL zm~pheM$G>T=@na~oQ@%;nSzw`{AN%XKZ&P@0wn%X@ zYCVoDT7Z_TEugJ(>KH8G5qiK@-A2&FX6xM(imj!k^)dn-)r%9rnAccPh$og6@6?kM z5N!#Ql5|Z)qPuP*+b$s}m2=zk_ zwy7M;7Fd+4Zoe+G|Lp@bhfaEj9%OJ^_KyJ#lcj#Rg+6s(Z_R$R`uB>qx-6YGA#wFZ zc&o(XpEt`wqsf4hyJd)3k=uTf<=t84ss&?U>&w!QBg|^{WJQeJIBzfL2FVd!RY6;;e^qlxo|x%@+5@{?`)hTH3d?E2Bev)w6t_bsBMqwsKshdqH&_LMD0<~XB0vvCA9(Q zu``L^_)^Qe@-n%4SNnqAevH&#fw*gfRej(j?*{#fxV%Lc0<<-r=sY4J=^Wql*Y5+5%iwgZ&h%_V&8V=w31x3ya!kOyO*U{sMhSk4GIQrY*kFA&;A#(J)cJ z)s85t-|qTsXsuGsm);MeJ8prHgT>F=K>ypDypobc`~oO7%^g8fiy8}>9%E;Z@%8!d z@_T&|<{txgmXMyY*bro*+*FGN?+D!v-h{Z9(zwW4Ke7XijwEPp2?%}!#-=z)zp#MN z5|Fml0C6;0JdDvJNe!aRiO>t?0|mdrtC2?_QOiJZ5d!K&Th2=~c9Fj2;970Z2OvN=0 z@Xao~mwO*GO5gx}caUJDhQ{Td)P=5c8op6_i;Z{w`oY9&ym;a8!jjMyNCUxWv$SSG z5ZF!Ddwg6q3h7lXJUKmZIx;U7sc;X$8ibor)SN3T6`qHqyRpESpmEq>el%V!$cjr1 zh_IMR7^Dos?g0;t9UESb{g!xo9o!N66A1_y!Ti+W|BOZ&fo@uW7moNN`$}o!;*jXm zf_q`X$V9^~OzT#J^)f*AB94n$f=>G~2#O8R*mKY@O>r0q#DQa*N-~C11LdcUsCI%w zArg;(W&Ax%4dU3@Q%XAqg<6B6wLpFLa@gUk4!rOvmy3x!iG)`t-iZ$fg6UeBJtM}gN=}ZB_as`LV$ch0Z|~C0i_RaaUcW{oDL>0Le^4}=q4LOfgJIb z!H;4Yh_C+)OneFxM^A_gbY8|p3zmiib1O!Vfm<*+FKay3%)H@)4}<>qc218Yghq7+ z736?mM!Ztw1}n|HD6PB+);yjKTy&(L0If9{%=2dv%YmQ%Swy^w#Hz(p!gi37PqYw_ zWiy4Z1Pc76FygRWa3rZH#ZwpiZo;wdry<*cHqPP5>4zz?y9mJBX$KI;|J^$VOFtmX z2u;FlNCBq@3Nt$_6+P!LCEfu(5auK(^KIQ3u3o{Cur78656bRjioqrW<}$r1=O&tdAXL_2eIsu9A2O(f9Xbk{f|aH zMm{Zx?gqv&2!imX>G{42!lIM`fGF}w7Z*Kg27vP*@%I;`AxQ9}&Ct>kjS+PW<58pj zvoCg7eA9Bq&7n9*`EvLRAAk@jt4#}#>jKgoXIvLDtW<#KnL7Deb`tgje0t3-( zFe(OuI=tXU;`*bqXsAKf@M4vM*Tp`MaBLMSV+L}}vqk*v>WgHKtbEEgm~{@Q?Fu_P zFH;yD1<+p)2H*-VCgtT2E^xRYBP&u2E9oZC8#-0DySuv$goa{NY7eOB^;Uk#d~c0Yza@1#0iUqq?dnxS}FqWBd7!It8k*f z5S7L!P9~OG8R+uGHul;sJe{$OeXgC9zkxR0e8>`P*i;0iC$VkXa`r(&qW4jagfYct zgWpHVzQD2Y4uq;ffUUJ1dgjfqntTBo-DCe^54H;uwgft^qV3Tan%^*6y?i;=@@}^T zXkLetpbE;$Ns%E%P?^h@2S5QslUSyY$E1QT<&oC^1Ut*xoTLqhXzwpT2$}bJL`lT) zz>+@t#vmR5&tZZqY?+vuzXKN7y4b zffoUt59*9&YZ}pIMiWFhDKK7RkoamGv}UH}?3r}!+$4+R`jpYc#__9uI#1rg@XobZ zF%^Jq7jOhvE^4nWy0X$zb>fj{BZ5IXRYvlfMM~L535>$fDFd_79fSb@7fC?GpB8JRn zfkx@WxnaOs-C?q7!H?hu8JxI90Y?SZ>0$a1;$7zCwJ3s@L=^ZS9Z^zFBROHFQ()RM zrzObu2sUZwIgdOv(&H5skarpG+rrEo;WJ3Zn0m-@$iq2Y)_ zoF3f!gj~M-&{Zwj|2G)Jxq!16!=)Q})Uwu-enN-w_DL+FNN%O(L@p(strEZwTQZ4z zF9+gw2XN6Yvw&&|J{lfmHB;S|q{OwrnEbbW)-A4ToHH)Z~HUtThVX(q%8 zg`1s*Bijd4c;eu70zp3fPtd?43#Yb(I6`|%tMGxUnnF2)|6AjAg-$^qsL@WR+ci@s zVoeZU=ceNYV=n$3azP@PEcB`~2$7okD$oJvV=0E&57wZ@Qs^(%6NS1Y^Z1P^sjwj{ z2gdaXX1Kk15&lVM#4$=9rqN7LDIhNPljw(JV`|TS=9@H&BO8!0Ww=0^G4LZJg2o>Q zKpN0yd~_p|^^OzYhvf?Ug5+0Yfq<=>sQNW9w&VWH_5Eq`>N=h>Pz2Xf`7x0j0jYgJ zu@D<5dYI2xI|2p}%b zz9)IIhX(s@gw>(FcBoABSCHk05_l38Nc?Bj^h2~^l|Oz2ej=Ui3FQ?eMoJ99Yil90 z$*R7Xp!6tnAdCmDZYGD*gQ7ru0e`YGkRz5Y6w^e{^%c>3O>tTaaJF&e_k-uo5+IH* zCoA)Qv=&i<_(*^8LE=b%k|_vjn=_RKNgFmyaEd>mvrzC1qrQ%)c+qowu{gxswG)|( zXv`64IbLu?H*QDlzrOqxCNFoi&|z!*3Jz3v!V$rNkD6*PW{o4Ck)UN|25W77f@kH@ zs!S?bz-e5Qll+?w1`?@tY@qz>`(}XU;nOK$n5zbaQ+9>fyFjRFpdVlFQzRJs1edfV zvxu}2KUZdu2lvOD=*R&vMdNlOk_mD9rFj7F*|U1G6oN=4tnh$?L7KC zkR+s+J0Qd%oYg+unw^og@|aIcX;JarmkOzrz1~-NQF~anCQ)!+G4rx-Ls*I7`DpM zuDVYm#Fo4dMAC&lc{ohle_s7CWWoPI2HeXDe~V=AkQ88e-vHdcdG~tY{kc*q6YD$JnT4u*O{b19ZB`VtpKPe z#|Tx0#S@s6^We3xw#3cZZn3CTmaLf(K%|hLh*$X_F{xw%q?HyC3@AyGyS4Gbm#Y&WN5>#6z)dtL=^o>e5@+Y?b6 zW;&NQfG{vhq_Svi@&r$-Ugiv`_f3oVTso)oJCeWK8Vi8mGi2rBDqd)H*VNG=;V^DV zb_owRzUh(O0u73agE;*HvCgS5#Mjdi$+uiom3>kHz=)XyA0@#laYJ zj*oygCVfixr+|z^`vPar-vhPCZK09>&_VrPv zFivn9CUo-czV_5{&m&QexIYUyP<@ZMGp1?Ww9Jjd`;)r+NC#c;ZK**e=cx;jeCE6T zo+~9$72S#TzkPfMK&)lQgG;lK2`MSEdn&J%+C>iR6Rp3diUS_ssUgQbF)^`3Z07cm zS5Y>$?x3r9U5CvNn{6{S=~1Q6ktAz6oxVlQP=J-&yl1Ua4v+0!RXzRFs({MdLlABT1 zn9%B**~`@1s0X+B-nf-oDdM^t{QI{6uwH5!87aZ9#eV(zwR1tkss6b4hk@TwJR_fD zS*kD=K}Sh6jaar5(BSosWw_1rg!e}UTHF2gF7L}PadC{9ZgX}G9!INjSAcCa@8s6F z(f~TCKAg<)L_|d7rRO(SK+O7Si5YA2Bd^l)jYjXe=tkf2OR~HXek*?7JG~CnDk5)2 z4}LZ$z%QOYdEzSUqz6ipeiXBAKu8!j1D-2pQx9FOnQus8S@Pd5+J1r7Q&kp35Lz0s zS~8a}NbKoTsu}4x6Sy4Q#E3V?yW*N~g&JTS1jP8fSkde7OL4YjDc~=XY);<6uN4k! z;G1KW?0!+t#RkKnv@-x(QxQIrn);RsTykd3Gg=2a5x0S!R_#g{f0^^CmNc^tq1aY0 znHqcFHj4!ji7ZkKOH4jrYlO$_++s-kOf&Q54f*A7{-t8gNcyE%Ah-FFniaN+MxRcP z4SH@Oe<$eCjURhLEQH5Dj#xrRv~Ndh&&#xwCCBQ1{CP#DP?Cqd z3Xaq$VNT0Ow=%6nM6f7W;|Yvi!>MYCnRHANuxr;&0e691Xv;7~<8w)9SOc$;P1ZY| z9%JqY>H_FBO{N;ea*+cMq5*aQ()r@WtSQ{H4oxBIrC%Jy>@Xu8W@_ICP;P1}@&~|o zx^%r*;Aa$VrLH2I}M63n`MIeD?LIQ~(CsO#?Ly0F#pDGzuslDLMQVdPx^^4~y4M7Bki6Vr33hZMT zo}@_Y){i~ofGHAG2wgCgaU}TxoB#=_L0)dCIUV`*QV0(}FxZsa2Nfcl4<+T^ z{W&!7hOIxK{~&}oKUul}daOW@!S6_KBATGPVN3SunT%tELZfq+EHWDgGAIcEp0usI z08*IR^a;KQCRqJ3OaMKC4YL6nSv#qMPDeDx;oS2gglOXuJ-86Ku6ebUA4|ZcI}Ols ztHhD;);}S@6zqV$QWtb$&da7f9bN?j0DmgGCFnxtk}rl}*#ROoJ`6~R)0`0n(>i}WaME1^ata~mQAVlh{CIpkmb1~oX-%G`j-V2a zN5xNcyQ%{_Ed$BkiKP$7~T=FRn~x8fJXLrp8ton zw+@T4?Y=;nVL-YYq@+O_MRI7QyE~*r1VKV#K#&%YMq=m=Y3W8n8VNy4kQ7i_gmcgP ze&6pq=l^qE^N$YC%>6ub$KGqNwRR_{tT1~#W|1VdPGdQsiXf+i3FA7#8Lv!Wd2!Jl z3imSEiH4!WBBzd;o-MWV#nadnTE&rh5;9%2P!qsd2gSLvLfuRrDP3LNo%(5weEX4< zNOSk!*p&WDo~$&?OrZ8JynGl}T9v4!6#;{i5E0qJPkE4&B$}Z%JuC%|#2QN;_qmZM z7lzSMHNPKO-ZYbobq<90E!xuNX=vQAe|GuoY1b>UK^J*^`k@&KTZt8A25i_&xLBWL z5R&NEgNdh^FY9ODUx| zPQEv;F9JTmz3-l>`*`GA%n;wgc({=<>nPd%@)0=FIw|~ifEz;*sfK9qH=hlBub2K3 z%p2}VW<_g>MnYGRRS|>w=Re_`(JAIGK1+?u9PjYRNPbKV9>5Xu0638f>IF*ZcL)jq zb(8Z*%dkw1kew zZULO{VGM-kvzWp{8a7D(ryoPDIgxU2@U+c|b74M!IJuqqVhlU^VBv2K!mn|UvH3e{ z7l_KFLg(^zV2$t*_x@P`VmMtHnc&5BAbx1MSdVV_0^|ViJzl(U0`-8t#23pd85s55 zVn@$r;?j?vX=${H^>}nQxqdb9Qm7FRx%F8m4u7(Kq$*+XrgP91aAD|Qk-NSNM(q%6 zmxY&>;pGCj7H_H^8&Qi6i11bYno{K}#rwyg<7~k}WCL}a12)K*{$?Tuu`PTVubYK8 zg!e{C+(3>*{%8jlT*HHu*wW=gMTwrDMdnJ}V*M^HR_`7u>+M{@3!{z)u$l5a< zo1x@_Pzhg9_>MH+5Bj@jxdVJN{^8xwrf9yLtRvZd8C8f;??P+o^6=2Bb+&{r z8?{5dfx6#Z{wQ~QEn)Nu)VgvBe(wizAhQp|uy`Lqhz)D7P2qo4-EBzTea}ikTmSGr z00b*#>mQw_*>+btrA1u6sCkY2))LO9K71)?g2B01TW|0sef%?vn9k$+CV6A;4&!_6e%>p@`4~17B$IRG);b* z|2K8#FCDrkaOS87LAIjr7$^!HXL!3*1vBq|gK%^`>``1u$ER9GM%7O$44A|b{Z5WW$t=babN%}g!EKMAF zAF0f^3T73V7CYC84K-EVR=>sTQ47WiKwnVv;X5($AK{=dP8avCa>jO@Uj#^3`38JJ zl7Z{KNqq!J^YY)fgZNd%-`dyr-VdUJs~+VDhj849~n@=a?HT+svwhLjx$MRcG%lK zGis#NT(^opv2kBhhin_i-X^63ji;QL_HW(N%7~fqhYigIQcJ#dv3Lc8Y(g@9j3hU* z4oGDY!CkSwLMYcjsqrRa11NPWhymyHwwLGbUbiWFcz&J=f}N_;qTR&&ZQ_48N8{!9 z+?9*I+ig&Lkg{vwaWbTe1CqO{d!6qX26soBU}DJFUb_Vy zOh!An3-vF!1x+O}d40+Yp3`B4D#3j^LwgB;mgiz^dz*OGC@0|0nf`5_Vf5)cTS@Fl z*3I;zS8s`2wFZKBIgwvXU*D5KaD)<@GvJ4ttSC@!1!H$HJ~V|~EBm-_Z|zqBOQYv6 zDry3y?$?zozCZ4Iwc!!Yye;l4uYDeZ`)v?%`x-F~>my=bb46gx{O40KXAY>4OK(tv zr4GeNGCxvgv`R(TV$2bU(0+sCkwY{g-+i$ewZ0VF-U6f`5|7)#Dj5||l z5qIc^$IgWD0IaIBq?D8*UWI|HUm#4l~9GALsD$V;OE`A53S)84?wMrj*;=_ z{k4#AfDx5&^iV*kQ`dQZn@Za@;1#1~l#NDPv3lkUVtTf(zd4+h`3}JM9(Jn*jB$lJ z?htmyfTQj^Gj$G{UU=L3hUN-Dpe5yXu&A%E(n~4xO=vYFwaJaq*q;d0F|8^5{XHG& z>U)NDGVc)5>!q%HmVXqCxg~QwNI3j@RjQ2Qt;#E+;?L+GjBO*Sd~*iRSFs;U2xX_% z+OP_Ldlw;sN9o^o%CmS6H~5V$JWRFPM1DC%LF-cjKP=NIYTC(2a1=r5D)!;-WKa^{ z9uk4WauE{`rs6iR8-qVR{w|8l6opNLyz^d#$<1^qF1>Zy_d{;wZFML(V#0;B7%NoRL3fyVKPyYzy9i~h^v_ZKwOe6eyou}< z4Fkl%IZ_I$-b`wL*syw1W($k%|g;W`>*xPikC zx53wP=KWio^)&cAPKi~w%uacEd1+&#K~^Zro=iPM-Yk>Dpln=iQ*lF2d*e^vEkGp4 zNKjj*=A8kA*a&4|{~ZW52#{GK-I2kI1O^(e6I((DwNicbeu>nv>;Zgx0j!|h1N9&60E4|eif8Zb~Hu#`WrmCmIS0F zw&2uWeUS&W7sMan7$p9RAm#xZ%0XljU_B7_%I|vB&w^QHezcz5TSNKq6~t>u13ns( z%0$3XM}o6`B`NPuIJj^^xWE_QoWQp)fOMbsz{LlV;Tl0b8g4i$kELl@6{$w&K@v~& zaUQF{+u|b~bSl24zkfW_@ran#q491+gK)`8IjB$2oMckPVd9oRSOD)2ABvAY323tD zY#OnqP-cn-HGzi;@T&$aWrn6X^LkoLK(egy`AX)NgBi{YbOfNRjej2rh6!*1es`56 z5&)kW(J8!nbGEyhxE`J*sEkhCrVPLYnQEsmRL$u4989<{X7&jHCbSVyosm4!&7=;3 z?OsRc7O;&>H|kzPo(Ime#;L{knf^1h)CctVb*ulm3Jurs`~LGM>qOv7)KJnO(jQmR z2C9Ho9Zyc71wNb+eE3r$#}Y);D}r8&r||w#cO^t>`*)U<@UwvEgIfvGY;qe5(hibx zjDbcW2@`Rko_I^TF$oH1h)^CB%^gTcI^{FwAmD1P{?cPbK#ory!Wci_9bbG`OZIRw zE-ix0Ej5$-AVCAN43fw8tKD*NiWM=UB_}5f0}eMbVDg4}AYXsOn#y&CW5pWUX7>|_ z!eoH#PlmPHq#+8#f$hkE57!~0$PWc^E$y<%Z_=20MURz7t0)J8{WeM7>5SDtCzV<; zP;PYU8N!h`I5>J3p1-~ip|;HODwbF6_0wBT|z-0rRQz7P@)LZSYX-i(WA)m zG#uMp52oUa`sXTogb507pVM&iAO&L~_LshN+k3<8PkO?q--{X$#j8L90Jax#Ul@po z?=`NP`EDjBWoA}XiSxl8TzzMiy@ZsfUMPi+EO!L7TD+X-XlZ#leD2A`_o&2wOWR4I z(j$Ax<*l)kq1HXyi^8%pt~U+gmR=ytxnNh=i2*iUT0n&-o6=#^k5g8;x#!hdEB2Ul zeRjlX;re@`W{F|71=BbwbxzyEY^{VH@Ic8>NOC>Q~cavv3Ub$Kx z(=TSP(<-OPV|bnzk+?25#%5g1MG#7`XM6>=hLb$sCTrV0#}5}-$MtruuKdP14%b7# z!n612u6tU8h@$UtV1X5A!@WUGnVOoq?A<4P5FHbu;^M2*Jta_P25>CV^ckoZ2!9v2 zw(8g^dS6mg!~b1?`5tTkk3+|m2b#-X6^}CtfA?KJ3#_&IOvucA;GtWk<@G_f2S@Ko zQEk$j?~%)MsBGtyJ-X!M$B%FFjwJaVCm4ZQpl)32W{T;|ALpah?)s+-g1s*u{RYXc zcFLFEfBtyB1(?!EDPXInR5byjJIa*r!xRlyVbPVxuisy}?XMbz$R>St#{7vE4zW^|&w6vwi*B%lqNcVyae^z2Hu@-?c(i#OHzj6d2@Za7})fnFUD*%#5-Da_k-$c2re(O zX(jOeybjkC>bSl3o>ZVrP)p%R(R?3Dh9tDRI3-iNA(f06muCmrEJ;6PQ) z+v<^a!r@)s$yS^L!q*649F7qQfWypAJ+dDrtHcHoQ=#wgtx4CW<)VWZF6ecu%gnKU znat0DYnl$wf|ml5$@*u@VwLlq*PDEvpq@e!SDwPT!-;eXEazam9h%XNXH9|G-9;Z( z3#AP_9Wx@4UAACh4A<<_d$;ivQbivY1BQq9)uCY@W|nFOhe)83a#nLoCd3Yw5G4B-^XkQRHC$V(?Z#=k)0K)&8n zaluc@#c} zaHPK!_o+5UBkX&J$LLnC`PQV>2jpD6Am1?QB=>x}3y8$OkWdp47B zDJ)rpYCBR@QIMAO(7X=M#pVD7co<9=LwGe6%NBkQ;I``-h*4f2)+PR9Q7sLZXfk#1 ztvn=H7Yj{3iXnH|WCfU5)Iz)@GWKvF#2&)@s}fJTTSMrRLmwkH(mQ~^q;>Y@=voi7 ztYU=OUTm46*ZVU)2U6!{e1YuahuDuFXxPwk$o72tA{5U~kF`B21H>U4xeo1L-+=t@ z3OK~lh#*p%0zd{+|F{;!qJjCgc@1(#aC=ki$MuFW=zocEN&r(H7gf}Uhy1tj43VS+ z!}tn{`mUs+&sPyNLaDeoX0B1GMBkd+7y#gEP-(j{5Zky{YXXD-MasK%WH=l=NP}0v z^8@mj&%syhjxBt12UnRS2zFwII`VlgF818`>LUy{O+R`Rwjg`|$9MLXA+z`%+_bDG zm7WgBEalL}d`V(jmzVS!g>u7XK)37{A5yJPn88?B+nh+il!qv1 zu8RiG%;LPAU0sokoRI- zGSw&Y>+rD3@aSm0dB-fY55{8b&=vhMbF#mw!KQrvwI&fiO2m^ECsva|*xfq;Bde7G z!AuJdhk|bJX#g6fiHh6^X-V+mbkwNjMoPqN&0N!5J$=CzP}sbd zi@3VNzYPbY7Z#=!xl^~Qjs^&ZPl+aYkdXXIz?u4&00HXezy&6|>eo6R%}o&k^$_ff5$@9X&tQU-UH%D z4vhwN#GtQ7#=l^Z(dx<0w7xxJ`FHMA!<>aGIFd`BJuo$@tE+PYA$ZGWu(Qk4A&S|Cb3R47h*uWgorv47?sla2G!g&ZVW9Ioo>g z_F(?S_1^Amxl>;RaUbOHGl{0pqNQ;WeM)wQ%?@3iV19a#bN|c23)_>8Y)|LZl}Nj6 zQIyN*sRAR=17#@Fakf@esOd|YHaKr8`)ZK&oQIlobAU=yb{AMr{}PJ9O7Gc2^fNp~ z^p`eb9l`DykQr8*6mn%LC@7Sfb|8_6tMfz6p`oG1VZtBZ7yVZD(np6g9X0D>Y)^#j z^?I&s@(z;`C(Oc$p9a~nU?@*CWjO~>O}NW3k5X?>yZ}(ZE|D3TZ8Tg_5fd1_82-iO zaRep$lP{)js3eZQS={-n&3y5`I~)Tc)N)0pMW6!s=~~^z^udVSTfc>Q_6mXLch*hN zAU{w988imej2)+D6;nd6bAk=J=^aJHw|kl{_vtyHQ>|W+K8^HoW>X(F0EK>kK7Z?i z9a)U9q42*xZ%sA*#Ou=}%P+tVAjUbrIf3l=euGd@BL5>sL)Zj`23? z!76;R?|3NgBb_zmU;4ZZIM1R4kpEH?l|Pcad-Xf}slfNtp}H2*_NJ(aN+h6(uf^;o z<{~mY`3C-vmQxP!q+~FYFpV;41Rx6rvcDFov@ll=VL8yN^gP5{QSSAsu1I%{u6y|^ zHgmYB-PU^l3Qc5?57a94eD3le^9BAf-%2qy#@+$7pIsH@XZ3O%NDKFFp~Rag=rH!i zKO`3^FxV@z`t}-xGe->s`j(U^aYX^{SA}$Hsm}Q2JXkIl%xI;;_?(rc{dHwUXo0sy z!Jpl@tDhv4a&d}WWiP*t<3`jb|M5wzrP4Vjc76B&&KLnCwKAeKLBwCuc3TE~=h2~z zR&n6m!5lksP1)hDE10ra5T%7Nvs!tx7AWt{ifbCdbtUY0*Z+L17Not|Xow>Byf^6< zm^&x9uj+{+p;e)c{>v?oh(}<_GKivQN3W+!1tN}BEnOuaw`W?0st-E718@Rsff!+B zDA|0a1c9DDd`=F6|2^bj100t04!NxuxFx(4`!Hae5qeq|dDMF+X`bixwiX7lsm$+M zbio>IC@(Qq1IatpS(jNyZYJMP(BYms=x`)EDC{BVA#vtYe(DiFrD7nlyp~DX&9E4< z3A#?&jtVaKOXMbF7fhJ4O{s@J(qe-=L2x?a%f8M)a7sk`1uXc#0|+xdYB8!Ig*_$Z z*MToUhf85NC4{gWO(Q5EZ}5E%`7yS>=F)@H;$z^t#X`(oBhWT@Y>~D^m?ChgEQMM?O$>?5fkJywCo}F zkTfP4pkuAA*Ech170ZXLeixyX7+ESu*WROBX*t?EKLWJ{Q4KVlZFh{PuC%?A=U$YQ3swlC#eeerNxO z2E{(x7u7FkV9;Y#P?gzjUJZkqfVFn~3WBsx!8)?R%WSJBE7r7GxdjAT9ahG6Dh;Bv zxPv|vlU{ZO_o2UjW!P(XVBm~kh1&c*f5XjmJQRWMjLk%#f(8Y)8vO#uCCR&{Bw6!g zO(lJl?@5ro6M(IYM3}&n*DU#P0F;ghlY@6LL+(*ZtRqZ9BYM}f5vhZSs4;4fGPk|g zhyC=x^PL#@L>tOv)IJa?b~b_sFaqX-C>)O4Td)zLM5lEJEkB}}5rp3Kc)cuQR-ZUK zDhUmvB`;uLWu{btnB%^6j|?Lb1q3BJ#~=fu9-FAw8*zv^^P{DEfCewZ)W8V{19IdW4Jz+g}wIw+I^n+g>goFSR%#Te4;9^yGW*~{ta`#yxE~z}ne|NUa<1@)bdG;$z_*xBn zP|xE2*}6~y;Ak{V`~hS7&y_~nK!7Z)Rw63#6AKHsgmNx~a6iUiK|Y|BoE?Ti1Mwbj*H6LhdMgT-(o z-)FV-#0VaH$w&!fX0jr?lT<00?$km4iOZRp^X?A>+@FZ3{EqK3#(4Yj-nR~^wMgZ( z8Cr(9;w%F&P&lFOzlPvMY^Y^HdM7D*`0+V8)7i)NB0YD8#+g7S%I*}6ybxU5b%)bfux)-^$+>Pd9wW@;zj=K zs+ps&P;iY<%)kETAQu`nJNlfItSj0L-Q7Ukk|x62jNz>fOuMt%AFGq~$0+`$r3FM8 z`~ghRhUctUxEz2Ti^A*>3eYDwjt>S9CM=P?>*)w8^xac_wcr5HO+NpXtw8c-%CC`% zsGFjvw1t<0Fb~5-6b1@Q{3h6uszekMI?bzNSb3E16U?YCi`{%0T^X*$K zxu?(}mM;&$O85()U@Tf|2n7KI5)B(A(-Oc0&L9d%(NS&atF!%FL(72V2d;}`de7|b zSwXe)zN4jZ=gglC#J~uZ)|39awLe*Kz0y#lSLa_Hr{op{NtF{6nt>J_CB|3lzXXxI zEd7gV^@o$^E8O?qd+9hiaROkB__idi@aigMn90d=wJUHj(L(8T)_C~*N&UIiM~}eL zg`KJ`hTqctmGvX~-LiTHKr9Y0{sKF~KC;s4>f-!iHR=@O6g>+I4Ny|XnM>D zmZL}7p@@*6?*!ov{Vx#GOx!*eVnd&#e|5Ay+fpbdEP zKBgKk3&LuR=#N{ecaMOVe+C->tuZh#%AjHZ90#`gZGNIFuiFCt95jWCEl>VE-VSkq zcN$h8bKDsHTpvwF)B#MiM?w}(gMk95QT5^sF|V?#mkXe>UqcKy9I<*jI(e_Y0o9~3 zil<}D+4Pu{6{qoC?F)~jjX3Q;DIJH&NHDK!@$s@6_Mr92fC zJSv*!ghwu(p%UU;TwF7E+7o(yG}b>^3Is*|X<+&gVcCqb)-yi_hU{>ECY`2d%iA9Y z>2VWOJ6+uX0m2SxE)|fV$ASJ0@e~rVY&Wi4Ul#7005$s90aiS zC&gmfmb6OTOc*^MWvTD@4Y~}=7kB}q9-@X3{L%iG%tY~u5PNRyPJ0SU%4xt9Q5SL; z!+U?rS)57oE3>{{7$kEAF)=amgAp+W1^K&A{7hOrB(z;RdCg3cEBsS6A@Q*^H=uj&Hl9h2hLWaVvIj~7ww-}l zGj|>WZTbbktSL8=8N`obmM6Q<3JF*-!c(vZ-cwN$V~Qz#4{?fn$vZT!388|dpMc|1 zB@%zIYKX1Qm~P_9jkf}bFDiemn`Q&Pb7}9hc_L<#SKHtI&-}6fnW0oaaLcHn@Rzf? z_qXcYI7ei(`&^!pF$+3u@E)KW7W%UGJBY0tpHMBh{SNUbYelQoE>IbJ^m)wto~0=^ zr@TFgj5v@?Xx82{SZbM~+6|r>>4FO7`M}&l*zKD{%R%GU`haIoBLKsV>H))px7|?~ z_3i7{OIm96o&od48?ZRicnwH$g2V43sZ!uIHZLAOR1T{bMKT430W;8Y$-RJ&z-LkY zYm-|scU3cUiNLiDn3ZB)Av#)$WSuitilN$VA6KcV!&9za2lb6t6l#Y?+fl8;$FKcZ*qEY;@_81a$N zkf>C)&b< zN$cb{Qjx70r2uKc-izWylGCG*#7K%NrmHzh`Jk3A0Z$HP4-jPEe@SF2L&%X}lB9*4 zVe_O?j4oW`2d0)n7ZV`l#Y6FAP-SjU!a@jC!Ec_FY@crYA2K=6S? zN1|mFqk3=bj1$Sn$JXM+$ob*50)f;*=6tMS3lSgawj#o#>mp!MU@&aD-j+6s`zHF!AMf9}JI zhFj8~-7l&YjeJgpWyE6DNaoNecoIfAzLv`;LBN~M&_$83#!vZQn-~52x%lOoQUNVA z3~3gj0PFL-qSXnQxSqTikW6{wZGLMmU=sIZrGT{`fTPt38>6NK%dA ztzb0$U;fmh4omgorpw#2F<4vF?*7It3`-SDbzWd!)tEW@sltaS$FM~lgHm$$jZ)dd zVFUP;%-XA^#s)LrdRSywhI+ zYO~P!&3NL_Xy&86!_h3?ACf=Xe)!Z_HC{$%IV|AdB$4~4ioZDJ$Dwnx03UAlAmzPWkE%~%1KxvQxOV7oa7 zB2-s%Vj8f4m($Krt_df_i=YZe<-dkJ8i65#%*#feOS(lMLCdcP>2qCCto0dA=-cGx z@>+?C8E&%oR)wwRUhTJ+a`?1oUj}l8)Z<@S_=mZs8$I!I&#xHS9?3r@$yc5O=y zg`tjhi(Q|GS*ClugNwgBS=}+s=M5>0u-WW-imOUYLZuw4IT>Q|(IJzosNyNxxdYEz zXo-T~6H19egb8fGiiPM6LK6*NfTSHnzfD25739n1>S-|Fbl6iCM%ApSW8Vdmop;N4 ziPB#!x$kf1og?g>PkJ)TA?Zw-sckS>iL1n>smYm9GPq_9t`v=bj2v_pD~(}%M`o1& zF>DeMCBA-RZK5c1DZPWe|FwxQPVy|x_EeaiiM5=XEhzrrJ+za`F-<(PbNCh~E-?WO z{H&GEILQ5wCajolWkx*6i%>(X%S5On)tk0`2}n9CA8A77^DVKXQentp)4iCXaMhcB zLQVFcALO@WAT!JXx@emonH0eSRM??RR)x6PG!y(pz6idCDLd`cTvT_dWdRlA}mJvBT2=Po-C^{$iw^nPA|LXO4j5*FDEO_ zrMxlM5|Hy++55>wt7=OK!9h|y$}m5o0~r)*?ovSwgpsF$2u$bhSQ`;0rX0439->Wq zajXlzergU$ZdG^8D9Tt%r;h?stg|hv>IR*EIJFEV(JgMtF6>fwMgnSytSH<1Q-np3 zU>IzgB?@HQ$_SVZqm=?$^7lg8n}Vz5$opC^IMnbK1i7$*Dt}KW(^{y8s%r*OWJB4; z=nSS+LCXq*8Mc}!j{k{x$E*9_#ImszAu~2YkzJg!C0NMXV+BWIIxtuwr~z1LWmrf{ zdI8wP&&t;x0xZ!{&1_aLdH?NYCEr zBhC#tCqfbJtqka;1j)e^R&#^0Gh-l>s+&ol_v&O8VG8TFt`;T~;@)&mcBQcg|Cu$K zaJho0?+vzJrU#?LdL*vyas-FhI+nn++?hkC640>w6uYACY7KRh!50s3x3amp)?f-(vQf)+1Zr8sN$j#u0(-Q${YmZgl~J;Xfm- zScH)OY4vi3g!OkhHB2VJ%ZT_t^8@y%NA!v))kQyi!QQON1sD9yMafC$o_^cnlazaZ#M+?csTCoy=N?Nh;WF3(_|Oc zwR#{dfhq~6%*Yt544ZhV9>el7h?bwz?;Tu+gx{h21)k&SmMhXvG3AElQyF?&3x7`@?_cn&yu!MXL_jJ;VT#g(wux0(Bh!9V0I#O~oSMmk1|8eFAvks{`ct+w1NF z4X-TGQ`~Uq-gf-_X9O>umgo6HU+t3_%0;2kzJ3RXL!}@yJ-QqZ2bE=8zHTcs>SKWo z%wP4k0BkE__m&REA)w_2LpaltOjl3-%3TaERZg)C6Rokf7pjRhbQY{uykT~^*B{Ag zCJxBp_a|)P?U<~<4zc=P%r5a;00zH16=Pe~EpCmy+o#pGZ4zi(<-Mm7yu+~2W^=me zCy3P`l8AQ*O>~+8f)bXMc9Hm^o&fhgpCS2l*8BuSF;<#OcO4i{-Cu1d86`qAjVz3+ zK)U%BX@9ym$~>qC^6OA97gBie0jpa;6kCsJl@xwoB5b!z1Su|p-pBKKtXH38ofDgd z6GD^M^sUCG;3oYmA|W7_$wT@T-hA7jpZ^=rY($t+=kG9&7jH8ozHfc8UTBl=+SSQ? zk~Qj%+<)^h@_Eia49wC4go%sO+R61nm*yX|ONB0d6|8_m=`ZVoNx#9itN6*xz?$eB z?A0`k7dfyzmepy5$#zI8k0_~>%0$eTF?5Sc>t(m8o2>Qe`rR%n`w9z+s^+utu9Z2{ zxbgTgVlXwCT?9b-ajmQEc z;h@R@JY#$Y3<*W1U>s-y4l-|YJ9ktdawAF^^UjUtGpLl@iNEgP%!FJkI>R_-;5>GG zVSTuCS$%%FkhU3(Zfy+RJ;=8&+GxTKJujGb$d15hxUl585GwT@2Qi<_4ho>m-jCJc zqi-WOOpY+ajc#&^4Ft_Z@Y8(B+CMVhi)7YBfAhm&~vH&SK-8ohT zaIp1=j2$*D zUE|gGshFups*E8@opl&5UOX$pJW*}n!o7rV5*98SUgo^7$?80L=E)3tJHAYN6pPvW0fW z$dwlxJV{t+z7rRY^+}Zl4wb>a08yOy14q$Y++WN$*bT@VWgj|rD9Dccfv_@CA)3Pa zTbr*BAO;<;g)^*=QrkoOnSQ*-7=PGs&?l6%`sj49aNZ*)0MU`3lW)bU-To5L{se{L76Tv|BV$2W_~PiS0Yd;e8*7DT$h;<*`+JF(%Zp^Mq0%uih>s66%&}zyr`H6y&UCFa~h72J+PDDAaL4OjxSU zI0fPr6S-E1u@_>TfQm1i;PWW=QX+gG2Iw>Hz^u#mc)h1TtP=XadHpFbMqYJoN+R5)b?SbM-Y$H9>C5|8@_hN8f-j@r!uDC7 zuXaID{bZ3xoR>?U)9Rvu*P_4$#o3_|GHpZY59e_|ZENPPv3cj$#9gQUuxwL8$eHR) z&dmtsX<#x{8f3_p_0nCRRh%Xm>X<~f$BjSAf2z;29KqeFn2wd70LLNa$IvS0X1dI8 z@rrOiyO0)rHJgJQXUjg3`#K{xZ)t?DlaBubU`Yujc zD~k7j_;Mf3=ERUEdyESiX!~8Ei@*-8S!2dvKLfSXUk;;3&0J`=d_sk+9#1-`1St*< z#?f}(bedN-!j+&bR#}Kgn0iltDeMLo3^>Q`mj}Lb#ydMqK_6eCcf!Hbgmz=|7}vr8 zU7{D%i6;SO1Nb#W?3&qqz`I?$?uoj&sCx9)2&Yd9q#Aq$BZZP=olvaEroFi@_1W)X zw$<6sT^)zP4s+D+)f{OmT8{X%dtb}H_6x2xah;0t~f?MViDjUT1&{31Yeb451@B6zys@(IF> z8$lXe?`CR+&?&Nr`={rWE+TR^&hjlZV)hv(z6LI|L3gS@q5hJlVScG2U?AW81nsCVwyW!HLP1XrixnT_@^3_C58w`N^PyUmNpz7Q4$Pq}E}4v@0q1 z<1Q3GImC>ObL5Gbz*D}A$=ecj1NX#aGH*^foI6QYEtawD$@kM#HAmzVNBvs(7lm<& zf*-baeBjDFX;E>j_`ZO&ArcV9?~#*}q}xbC&rxk`z&<}LX`T9B%!>=v!p405DNu$u&z6M&7+r~BoXvzbz5gp{ z9MW2goAGx3lI*3Ej74@K`QqJ==eF*@h%ZK6%fxx^_nVlk@TZP zxAC+7%b(9*`g=+JJ(<2j9K%J#Jh$Vkw8<--MGY;gr0&O8jr)y`SNo{GiJZ}GAD>e* zSOkP$t08u;nK|{0n=n$g4E)s1hWW8Ko|!7l9fVI`x#Sj*hnxL)V(;TUJXB*@HYGQK z41QM-(wxj~xL%lQF%@)2|8rlaxVOf`ns%#h=n@cILUiMK6^XT8e+2A*SvI@SKX(uQ zNZswtaSZNRKid5)n@f&S5p`!YU%yS}154h}&?goVCagOb*rK74 zIX#w_(Ze>#xxFe6$X7ws%P2}EDhwGnZ(%zRf|zK zK;nuwCIT(FYndg;Z#J+C^8hV-uHu-NLwNn3hrjUaje0Vz;v*!j?LAJT9kL_M+g<#} z#p>^z$cIB8v5q0lIk$JRUO73OkEcaa1Y19-Nun#L7D(gFN-qGL(%rewYS05#xag@G z(kJrhw~q1M@4B6Dr#d*3eJV`q0g-Y8P+{;L^Rs#D&qPfd8c{avK-+e`{q*Bqbe>3mbh?L4)~di8Q&kEW5nBzSpn>VxH9 zf=*qU(2d=-O6jEJ7#*{+!06p9-sQ}^*kmS1#ubzPU-gUy%q!R7h6xd2+plfsQ-_~& zYrp8NhJ?stB}In%lg&cBouh1X(9TTz^!`*VrjizJ?rjCyCE#8-Zl;!MDOVP%)fSK zSR&JBD9XR2L0|Oe{BX{7HjF!wwrXUTDX;s5Gldh6mlUvNvl}R{h$0DMm+n6@X%_R+ zV{pgvUa(KK2GqGUQFmrkdo9ia{|{2vN4=uA{)KFK>wTl66~=WTX8I}gU#fvlqq67b zZPJ4L^2M02lNl7t+@Orw+2v^+EeLB;Co z;-W8!vDVC9!UMQi1TgR=^uHCZ2lC(zd32heo7CmhMYPYFnH)L$+NJ^Lo)S!U>?QGC zydtwR$@S$~hWZ*}i};mX` zP!4LCF~HT2QNgh$gM~;E;@eF?BdolVrx_;j7y%0*ex8g*vuEf~qcW8P(}_f=#Ww@; zHY2n+W+$oI;RD+w5Jq_(ev?ers;OeJ5)TH&$TiUUyB~D!daC&z`yHGM^$A?7m$T_A zO$>W9e9$R-_OF0lypfDC%^IAHv?}ZHCx$(wD}$Y zA*IM9{<-Ua@15j%2!KY}IfMxR^Y4G&V*uIz|F3@+oFvdYQ=)wj^uMA8F8)XQQu8v_ z)=VlYD)C^R^F)+i|Ni}Zvf16%FqzZH+77!Fkm~OV2?-ey&0XY3_>ut6gDM;D^P^h( z5t4xtnb(+9z;YkVH-<*MsV`o8UuK(Xf=TIC`yxsGu)jZAd+B|)Z%jr`?z?r#axMd0 ziNMLeiBq3ZlnWgDt?d(5CP3Nm*^W8O0nfofp|>{hbWX_r8MvIXlnS&8C6=9mJwGwZ zsYhB%>g#9p&x6~dQNEAm`dC!F>VK&k5Qgt2euIi`$Ymb3emnU3HNfEP1&!w9FkbYF zICX&WS%b*EibJzmRxMkCI$tSnK!MVnjMuTJn$NO>4s?(k7)=)xcy)47nxG->nsaD!lw4)kY^w&pwp<`GcPF(CGG&5cut z78*+>0Oe~tf&w=tXF2CHmW0H_*C2&lNHiiV{q#ZshJ3D(y;sXM4f-|oVtDm?`Bnu1 zrD4@2=%Q*r5JMgR>c#d%Zf?lO^+nJ~vj#LV1mhYc-xR>LR>*aDgI}*njbEbI0V2X+ zD!t$0ekUCAq$J!+&e((9a>X2isMhPte8@_GgC@x zhx-EveHdvJs_Ps48C4G>>OX$Xa!dy11N+#GKsO=^f)@5+n#ZlJmYVdBU8zz~i}63N zu2TRPwRA)p!T&d%m3hH>NGz8$rv2alJ(mPCyvmK4BL2TuB4F=O5YKmK%`b2M_kU)Z zjD=6SO8|E}|Y{*VO$ z?g0jONcnfO_-D%*29Ko3W3LXNl>closBH)hJQ7m0EaP|oemv?MxJUs?hCJ>3duB>T zMh0tMC#XK`yrl;$TnO&QKy7Mj>ITS)D-6D^uiNaebarM0UyfXRyN1Dmr#28wge>TB6PW9j^6;S`&`2cwY*dmS@Ta?Tn z!od$aGx89~c0BjK;DbVLSdww(`w^+28-F}{iDr^E&{>WF7!9;^9!UYCufjcz90}0y z3y}9&LH9zgll=u!+6PaM$gP2(3V565k6G*wRGKsld&k~E-7pa$^w}%LtGR?j<{+9g z`r?^ZWo#o)a#)dCdR{EGkfQFuU_*81jmaK3u^K%#RpyR90wq#C(2gqU(kQgYq?&m| zHabo+xBEZ4{|BHk4Io5t8NO(X+fZT%DB$c}Sp+tT^?;w@*tFyk-Ej#O^*;ip@lM05 zb>Ft$odW87y((Y39$#2eHCVQbQ-!LhzYQ&X6?A1wES@dDzC}&kt|zE58dTD+DCkGd zzsPnzWS%^l)F%QY<;pW)Z+hW#eRJM9_dR?IlU?=)Xl2y^xI7n8FAm$!odHj&A@Juw zZS>KXazj3INw(~i~~`@yVr=X z*$Fp#m;W~g4urvyw4kpN+la8^M3=V`I24$qrB}Rvrl1o`&1yjcM%RH47XHJO0I>C3 z0L!)9`cJbpcOQ8akT@8BH8IOGVi5Dhij(*i!OK(5MCoS?ddd_Y5;91|fOasugIm%c zl|XNdI=rPsG;9_=NR#Bf*9lv#98~Gs*SZJh@{sy6o+`e7AH5%pw7C!6 zj3B;;9LST!r|yy?lne#GqP_!nDC;TFc1OytrCEP%r9Zd;0)D~ z&!HJmK-&Fn1XZ5wXl9BM15tRmr!}tW^w5Wz@zGJ?cDiL%%GOcD!z@P$1+5ehRpauE zC+~i}+@9#h8HajgJ0vo3D`rFQtp2z~qwaK64|Eoe z&Y5{<03~fH(os<@PSs53nF5gE|I^#Ge?ytZab}m2qaByBtXy)x$E>B@ZC7#`_pCia z(bYuKxV52{(N5>2wxbo<5L%OaA}mS_5}BbImE6^;*18sQS*`s%PH9emKy&7YIq!Mj zcb@n8KJWMY`99D0`FygO`TCy58L7M)FcmqJ4*2bf7yNXT}C$B6{YGpeJBY7yWRIzk_7{-LGeEq}|(MvY%5^k{{GklkIC{o_sT3L{H2^ zHopu=BhXW#F0=dXz$ z+fPXoQ`ygTN`D`n=$Rx0Bt=gADs)V>4gUVN!iIC&Z{aP8dbwl)(@*(ZtXV9Bzmqgv z93svi#LKoEK}i(wzOdr%@&WvKN-jX`j z<}clrzNsLmg3+qTxqSJ_@YUKsHGkH!BK_y9WhwfAcctkXMoZ;+E)Zle7_U$!=w!SX zn-4=k0Hm(3&nx$1)RBkg>Qc-uww#<>in@BzRH@wB;m%U+^9Ei|#~b&Or=@#Qumpav zz4R&W{&IaWw%YPR@#Pxbd{!;u$q9ML9*R$m_p5+${LC0!=8jE#$RwNw1=Qv*Ikj`G z%*+7AwFaInIOfglM{|=ew&^jTzgXlI1Xf7kxKnOfw!6XQm{dQ(CWh&8VJU8{sG#zY zLJVxA08VkiXb=aUJQ)zASfgw&yena67Y@+$$` zjnc+Y>v-)Z6gM^?F}h#D-05Sn4&nIqSnp|Q8D?VlS_DXwOi++H(o-c1nb!2ZjTH15 zM~lFW&apx>fCMFp@Ub2R`JKl~j97ia4N}c~{#3wsZ%WZ9okvAJr?GAApFU7p-bT=l z2P0?};`#OdZf5Nj$^7*!vg`J7Wr+GMCiWA}Posz{%d3An_sR>hGpdqJ zE|D2#je#{JbrSqwpp9o0NK<2JG94Uy*|NY6mkLe2uWW{vMo7V@7Lt7hhE1VlM5VBF zRXxm_T3eq{r%5 znPw3}L6x#yYs%}@RmtIg7v~_<5+x8CY_7@me)hm{-zz-#2dAAbjF38ma3>7 eG3!Tc(!^p^)TZdIf@pOSeC%yptc$F69seK1&(E>| literal 0 HcmV?d00001 diff --git a/topics.Rproj b/topics.Rproj new file mode 100644 index 0000000..8e3c2eb --- /dev/null +++ b/topics.Rproj @@ -0,0 +1,13 @@ +Version: 1.0 + +RestoreWorkspace: Default +SaveWorkspace: Default +AlwaysSaveHistory: Default + +EnableCodeIndexing: Yes +UseSpacesForTab: Yes +NumSpacesForTab: 2 +Encoding: UTF-8 + +RnwWeave: Sweave +LaTeX: pdfLaTeX From 2c905a0a543b92e0f769c4115414729b01a28ada Mon Sep 17 00:00:00 2001 From: Vijay Ivaturi Date: Mon, 11 Dec 2017 23:24:04 -0500 Subject: [PATCH 2/4] mixed zero and first order absorption --- absorption/mixedzeroandfirst.Rmd | 60 ++++++++++++++++++ absorption/mixedzeroandfirst.md | 55 ++++++++++++++++ .../unnamed-chunk-5-1.png | Bin 0 -> 37027 bytes 3 files changed, 115 insertions(+) create mode 100644 absorption/mixedzeroandfirst.Rmd create mode 100644 absorption/mixedzeroandfirst.md create mode 100644 absorption/mixedzeroandfirst_files/figure-markdown_github/unnamed-chunk-5-1.png diff --git a/absorption/mixedzeroandfirst.Rmd b/absorption/mixedzeroandfirst.Rmd new file mode 100644 index 0000000..ae8ccf9 --- /dev/null +++ b/absorption/mixedzeroandfirst.Rmd @@ -0,0 +1,60 @@ +--- +title: "Mixed first order abs and zero order infusion" +output: github_document +--- + +```{r,echo = FALSE, message = FALSE} +knitr::opts_chunk$set(comment='.', + message = FALSE, warning = FALSE) +``` + +```{r} +library(mrgsolve) +library(tidyverse) +``` + +## Mixed first order abs into gut and zero order infusion + +```{r} +code <- ' +$PARAM CL = 5, VC = 50, KA = 0.5, GUTBIOAV = 0.5, GUTLAGT=0 +BIOAV = 1, LAGT = 5, MODE = 0, DUR2 = 4, RAT2 = 10 + +$CMT GUT CENT + +$MAIN +F_GUT = GUTBIOAV; +ALAG_GUT = GUTLAGT; + +F_CENT = 1-F_GUT; +ALAG_CENT = LAGT; + +if(MODE==1) R_CENT = RAT2; +if(MODE==2) D_CENT = DUR2; + +$ODE + +dxdt_GUT = -KA*GUT; +dxdt_CENT = KA*GUT - (CL/VC)*CENT; + +$TABLE +double CP = CENT/VC; + +$CAPTURE CP GUT CENT +' + +mod <- mcode_cache("mixed", code) +``` + +```{r} +mod <- update(mod, end=72) +``` + +### mixed zero and first-order events +```{r} +ev2 <- ev(ID=1,amt = 10, time=0, cmt=1) + ev(ID=1,time = 0, amt = 10, cmt=2, rate=-2, MODE=2) +out <- mod %>% Req(GUT, CENT, CP) %>% mrgsim(ev2) +plot(out) + +``` + diff --git a/absorption/mixedzeroandfirst.md b/absorption/mixedzeroandfirst.md new file mode 100644 index 0000000..5865239 --- /dev/null +++ b/absorption/mixedzeroandfirst.md @@ -0,0 +1,55 @@ +Mixed first order abs and zero order infusion +================ + +``` r +library(mrgsolve) +library(tidyverse) +``` + +Mixed first order abs into gut and zero order infusion +------------------------------------------------------ + +``` r +code <- ' +$PARAM CL = 5, VC = 50, KA = 0.5, GUTBIOAV = 0.5, GUTLAGT=0 +BIOAV = 1, LAGT = 5, MODE = 0, DUR2 = 4, RAT2 = 10 + +$CMT GUT CENT + +$MAIN +F_GUT = GUTBIOAV; +ALAG_GUT = GUTLAGT; + +F_CENT = 1-F_GUT; +ALAG_CENT = LAGT; + +if(MODE==1) R_CENT = RAT2; +if(MODE==2) D_CENT = DUR2; + +$ODE + +dxdt_GUT = -KA*GUT; +dxdt_CENT = KA*GUT - (CL/VC)*CENT; + +$TABLE +double CP = CENT/VC; + +$CAPTURE CP GUT CENT +' + +mod <- mcode_cache("mixed", code) +``` + +``` r +mod <- update(mod, end=72) +``` + +### mixed zero and first-order events + +``` r +ev2 <- ev(ID=1,amt = 10, time=0, cmt=1) + ev(ID=1,time = 0, amt = 10, cmt=2, rate=-2, MODE=2) +out <- mod %>% Req(GUT, CENT, CP) %>% mrgsim(ev2) +plot(out) +``` + +![](mixedzeroandfirst_files/figure-markdown_github/unnamed-chunk-5-1.png) diff --git a/absorption/mixedzeroandfirst_files/figure-markdown_github/unnamed-chunk-5-1.png b/absorption/mixedzeroandfirst_files/figure-markdown_github/unnamed-chunk-5-1.png new file mode 100644 index 0000000000000000000000000000000000000000..e77d844629626959e2c4787940f5917f17d33cc4 GIT binary patch literal 37027 zcmeFZWn7fo7dFfcj0_;%NDG2AiiFhAB^}a8gVNm$-6aT04&4GuNh2VL(kb03-TmC- zIsbFc`{Dike0luM2iVTu_qx~KtFCozB2<)QurVnxQBY8@`f`W!^Eg_*ICm}(r;$&}OZDWpt!XFc(X<$XENhH$fZ26F(U)_c$#HC0Sn-M!v zizmCf!GCs^*DhI(cZz<)>~y#P{r7(Ur;Rt+%2Fym+*@fxovn3nTILhnX_0EqcbJS* zqqNK-)#pJf8&L(@9|WI%ux_ZDGvD<&%Gk^F-u+^^2if>iqMe*T1R_cO$WHoX$vG@8 zT$nXJI1E90FZ-n`6{|t<3ox%TmV1rfyu?>;qhDjQw^`+N(p1E>jx;8p*NU*h<3Rg` znm^n?UY7|pqFcJurzRV9+^|kgulhyPi4JEGzZh$a6}F|3aEPY*>+lr3xSw6;()X{p zE60rxim$gzHdb>?(p+Z|RStW`w^_FBRv|~W9fboQvP6RSjPrn zhAzFa^{u5`m`+Yhz!$w2KQotFQGmZEsb4mYUec14S6&lI=eHird&IWH9zMXg-mBg| zSz40VmEP{wS@%r$@K6w~Euuj$^`P77%2qw|aD8>_S?UjOPs!`cCQsQ@U#BdAolT#O zqu^wA^y~bTq#)wXpm6_eX|dDwrQ50*c|26|*A%OA{v2;}ar8Ks4QvZG1{Cg-@{Q6c zkZfNDs3GnoKSw#iD|GIasz%2sPqBa5+P+Zj>L>Nd@(JNB4e)Ql{)_=(e$N%wjH#%N zif_gYyEjiVX;mDAbMFaJIT3lJQMLu?gY`C{!{^l>`r2E*Cf)Y*7zJI__gAx+g^2`` ze4@CHEh8B*KPU~+HzgGHW8(vhJ|d?j_!Tr@_(bkB!SM-*-ZN52@idY06QSAO<)plg zl%Wvc_)SyvQ?-ILsab*QM#QVTrz|h@)hEi>aB4#z{&;)^F=)R2W>lrW$;PKzyF{MR z@Us8H=wkc&$tmmg>%0a3m6r!~PTjgc8$Z41_-&)@?5;^+dvBR!f9OM3L>JYO
H zTxq}L?;qAHjtzu#>^X3x%TgehuP?O7m0~2zj)RWFs@a;Ee3@^)HTU`U`u6$u`?ec^ z69j)`EVnFjIM6s7eLH=9M4>zHP2x#bXStdDrc9php$#!LiPKl!MOEV`hfMdOrsArT z0mpS(aKkszzq;SnB6{%qieo9Z z8YP($9?<>Jj%=V8eGbBZ6g>55x}7RP)(2~0^z8R(>*W&!hxSo4`&CV%%AB_9co%g=TmDA%BK;bC=ekLUmL5e#Fay%5! zzn?)U4Ad2Z0;c`v^NJP3jHq3bHvM-91*`&d#Qy)RC69txR9-(s_V4m=6e1OU)c;PZ zC4hrpv)b|8`_Jkq5Tc3y$OMO%qaol^kq=@5{xb@Q^z0vbfFJ&!NdLFAG)GtIxGeqp z^=mb!wwRchYSZpmgs*zlo?e|@{?yde&1IRc(U0CWS0{yR(g`Ym2l{ zIJ^!xTPMj)+`4$-f&n2|q{D3isGBQ^$_3ZgKBc!AJ?%|_4REYpHEGt z*zw}O+{i7KiK9Ngn!n4`CXodW;T$6hWVG>9caK<`6WG<&mB+l7sNwGR%DN*Io5y+W zMJm4|{edI#@ZIg1sBV)7ca?aZ*U4)|<5UsfhVf^bSN-Et-xHbY+6zsLXuLM)`otSO zk1TXs{e&G`b#Bvl>z8r1=9&e|_3C1^%CyDU!p`RXTl?PcqoTpu^nerXr7%av4()u` zHk_~UtaKt6gRt7IU-Zt0g^evQPUMpL<4$UMQ+|F3wZ{Ue_vg=1C4rmck-3 z;NiW~a!eVL!rrGmp~WUfS-z(_4QF%S#S&pSNvWyw6|J|=zi%(}@y+MGsc23f9vQK2 zgM(bYGo}mzzgtA38^KV}nq&YAJu{2N4}#4=ZKkRtE>E^63T|6`uhdjkRntr(F1hVy z_$dV4%7(6^I9p^p!f*>?sRZ9f`WjWf0{QN@qr05VxL*uksx%vb;Q-8*nl&37I)K3ewZl zc};)z_ZLJ)M)G!RoF5D;a9VSiXdUtwkGmU|IbW=%jHC-V^ZdvbeL8m8dt>}#9oS`Y zT+!>w+QE&j(uxA(j?f`h87e6o&$Vq>L0K6aV+F=9QJ+8*XlE7*(ct-0X)M|Wa&Xf| zc)9Uyg#T4z+|3ceNyhBHvd{rqUm9N=j^J;u@&_d+H<#NL6GZ%mZ5>f@9&_ux?^tWP zWJ^xoy;IP7V&8IMK_TFj_brjsyDw8XO-cAfV~3AZZELhR@zkK}jrca%E#Lzur2)m- z77gJ+>$6ZAr+*3m&f7@2y*#Ef+<5A`rLlV1yL+a<*{tE*?(C7@+biF9`8}$M|A6Hd zoyJGybF4F3;7dc#E^S2KSJwa$fK-{`ir* zp6MMc9Yc}d=^>E3HdXWHY#WO+b91(luXx72cKnG#Y(=kJO-pi0%1EAE;>OipYhh^U zz2`?ipGFYhuB@1`bKQ8ZXUK@&p68wh9yZS3#M-sqoZhU!kx^uh4lrQBPv2z4+IYd& zhCdqTWK&9QvGoSy6E8oS79q^G?>1%eOiuk|T^69h08DQhX?l9MWWxW3rAZ|e zXPB)_>^(6QG{7v7pZ}k)VHWN53;}-Ye}Tv&4Txy=KWK8%|BaW^7%&hvJl%A%T<7y6 zSZrz-38!J1MK}03^2-w=ro6Xr-`)iJ3TDrj96NSJ5_h(bX2jRpPOrUh{{H688|#b1 zRi(C)lLXTbyu7^63tCxRliZ`1>Wh`+jQnqENbV*uZf-UUv%)391n3SP5k26waS(s` zaO)w_AWFT+-L-SW!M6vfD@E?%9hH>sGmcxXY_f3$ifPhOBW0dH0F zd=)Wwy~fu{rzM^dZ_$^U4;gS5H@aYq4Qdf&%FTo?~}Outm{PaXnwq} zz+x~V<aP(Xyhwxq&?>HnJ?jN=2z0hPii+JBORSAoD}!gv+5|LqiVvz~+H z=%bZOUB`;`)C~-%4E$$n3YPrd0g)z|J#F*7X{@F@OnN|fv-xbUu*$TXLPq2Eli7EI zt@#QXS9UgyFKN2-y!LZ1#xyluHjKj_k4(%qx{s7rG%HT~}9kqYV}Jh0DvYwXK!|cDl7T0x^^VPi%CEF?|ld zcjr{(G7AeEgzQ@Rl0A1}R&2eA4L&YxY013h+PSyH-?6BsqeE`1-I^JSO+)>$h~UC` zy+ye&%k4R%tIB-s$*lXl+CXM=k@wlId?nUcL{7P0aZ^rUEdNrtN|#OF@|&X=0edwu zF|pAdFKLi@*ZLC1`jGHQfh0B-Qy?q{!>4}#0fo@9N5jdu9^vr2D$7hI%NrjkvCDTk z*()m|G!w?QM=4Ps8(ED(Bc6l-RhH@Uxbv6O+F^l5VzkC+B6>|p4LCx9yA9hPi>~4c zPCnoWaWk<-`MI8wp`TJ3+yiB_ow1HH5=qI+W1eh{zx$T4c7Iz>+&J5j5rc=1Z=`F2 ziLu&VK<9%838MzsN=OS_I&NE_J-fa&{oOP__~>7=LqSWTM}cTOb}ah!fqh!OsHH`$ zuRr@GErKi{qmiY!x=g&vFk6szK?&??qO`M?Lp3@&y3!A7@X*L~YH!qzO~&0ByYqW- zcRA2~yz0(c)Z~QgikH}nyD^`(+j=4Y{H$74+Pjg(O?@KQPb_ZkK=AS=2`hT}@W(TE zA!+@W7sKwG6MS>ZoXr9+Eg;uE8?~hyxG(j|xNyKz!=D>%r9G@;4)?A6gZKt29lrs8^>p z>Mq}hHF>UQ=7%jSMCCM})wrA{=}oQF-8udqz$xgm7oE7}@(Y5SVhCM+-9FyKB^iEe zURs=;ouH%iG@%A@vwpbP!?b#hSesAT8rSh%6ERB!pU&4Yy!Hs7cBb$l#dGZGW z7b%4!A97^|V>22zYTVa`YW%&@FVA+BEMo*dVwZMzV(Iz8L-aEQU2W9u>cvXwMOBETGnz1ruvK{O?96}o9cjkHuO)fp7Zxp) z?xV(?cV-P&_MED%hO^c=9&e<;ph5!vyDdTMop1D)x3i>HmHPxe$5HAsYL1&6iUIlu ziwa2@V>%W?oCEv`lT+jzu6+afaiVvj7FYxAD@XlM?Iz2*K;E;#trI;S2+^4)OrL+l z0sLF)eEaj=)BxU->mIt%a$05Y&Bp!o#G3VYlv@{;v5#6(uGK9H@0U}%&7To+M#!EL z0^U&6IbY9g+j4K$Amo9Mc$?px+ll0JZRT9>uBJ3u;z2m^5-pVfeOqP;#U-GR&33(t z*SS8ubi(fQF+6)e_KLErRH-8Yq{CP=cWsR%)ooz@*^k#(VOEmwvbmu3aFoW`Fl=E4 z1nF)`#B!19Psq;PGtS?#znt5wXV+|1U>(De${;rJNZgm*;|=Fw@22`P($gs&?%QBA zYcgj$#VWRm^FzVZbo)DTZ>kfh-AA)QGCm7!BFEDT-&Dm6^|lq-zKY@8?tI*?3449@ z#YQZ08F*|e%D@GsquIh$%mmB6E?TjEDEMF;9jn^La2lRO2CjrZ3T=E_uBACa(74BO z+4b1=vO6lic4zW2m~T(;Cl=eQm|xwf^G$LKY32>!-lM#Sa;VafQ?7mbrU!mUrM#La z1k@2%b`5II>`IjWB^DUA>L{YHwuflE5jD_jq2~@TGx_E8v-JngFwdd6m-7{QeDilN zUW5n)T%K;N;d^ai6Ji#5J(2T>d;R*XcC+KbUVmY!oqZ#8DPg+e*Phs~HHK>x5H4PP zRA@<;FyVKsU$zsKq{eqUFMh33GgrMyDD2jEy~^!8l1qTKMWH>nG_&gX_|f_*8nZhv zx2h*8DM_-qfi!%50&=`Qf2T|R%@Pb$MT#gcnrT0Ke@jZ>Y$(r{Yfl$Mt*0M4KX*Xf ze$(OnY7pK0J5JiK@V@Vuok2W*c4KMDLocoe#LHMvet6pe4td#$)5&EY4?;N(E!ntD(3mScx+Gx&Y1 z()ul2UNGHum-?bIiilJyPm_fP^e&P?`4igk(Q$_xn$N>DL5+9waUz9&-A1s1=EHgT zMFOReezj|Ec*2iG)4&>1IbUWhuo6QXr>_uhb$-5i>+TWHr)*Pyi9*FoiL|Ec8A@Af z-hr;r!(aCl#r*+QG5eiRWF3}PiC>iEI#;AqJDla_myYjDY#NIS2G89(&^fAj-F5b^ zug|j;qSB{VEE}P>1i^dJCf_)G|57-NI?zOKICFe~_|#f%GJKp7cN63xHqWUfY)Bou zdDOqnITkpq?0zDkj6_Tex-RkTRNj4s-8r4( z0L!%3iMeE#uLtw zg7Eh?wlpoZxKQjs7$t@a38?(S&(rE#DbH0o?uPQ$t?v=<`Ym@aP!)+u(oQ$W(JiY* z&9141O{`QPNmvJVCW)P~3Y<~3@kpDY#?upb@0BG@ZK_0IN8h^eTxq-I;+8|_Q}D-z>nFO0S(ZPug5OeraD3aHttMy5IG*wH z4khm}P6gh-lUrEqmuL|H*WJ3_TRj#|vTHuka%R%j-y&U4ms0N+$o!VUJG_)St3G@B z_DK~AJqD*?WMs9D&DeXlw(p>Iys4qp)gLJ%Ql!y^tIj8+#DY^U&JUUl0}8zllItEn zt}On=xZOydB~>oqyZNzkH?l!ami4vEoIsWTUR$X*D!#0uVz$_g7o_(qH{L!P$gnyL z$aP=gf=*%Q_5p8{ORc&v!I^3>1-Im(AQN0220RCr02#{`jye9jx0SmS8WV%0^BjfG z9@H|F;aZWQnb5TM}xOA6rb5sR|P~w^4O`+u z)plOh$uJ{0TQ2aiDf7Z;Fvl!RB~Tf$Qr;65nD9j1 zEL^SoL8ogx%Pr^XU67`1gqFw$iOr9g&I9yxZ*Kj=F9hhU2k2Wkh!b3yEfkk@I`8-v zb%-L>_Jv|^sJ1`%JjGH(NE8kdG#GPPclDK*n4=VZeD5m zTMX88$kLJtvPyM9Z=%%H(t5r8y)VXZ;ItCP#?m%n1|;QGZ0Yd;BvIt;?z*;a?F=I< zYC5|66-MRuzSweTuG!lffIUXb_4t^=Q;`*f@71a1n)ySb>fvEE0JF(c%J^Dg3)PIH z@+^xR9UD^>3wS=_?XhcFHbCDAs#)F=d+^_S=N^!<+`~lTBmdS)K|l^C$efG)?+^+I zIfe>p`KJm+7U5z5;y{%&bnoBV6b?s@1%Q!A&p)5v$a)N?*9rik{uymUj==&r|0{|c z9*}IHn^%MWTY4eK5N)LY6-5w;R5d|Qrv4Ea_<;);d#190^tCIiMO znD-(76$J-?k#`W$y#Gpvf_5Jm!zMZh|5p@rkU7!d__)^PK{ri=sQKmZ&jTCDdUo05 zPc2k*bUp}0EB{wYbV3zrPR80ki12n=?%!Wk0i(fa`2(`#!%un924rxclXA3O{om-Vw_^o& z5da_(!qKonzf*m2v_5ivu)O&di!%@S2iH_gXfMqehYr2NFNy0NgUJY?ohfZ+X%W)1GA={VoOOSy`C1i0=zBUu z9(ZbUvPk%RKyXjU^nR{PT*1tY{?r#|CcuO!x^G%5KlS!1BC0L=D7@A_SiTT$`f*mE z!pb-_H_oQ&H=hKDojkfEV@3Szd@L_#Zq8_KW~Soh<&`rOV%u6Zbld142SxR;Ryzr5+Mi!spxb7-qMA(beQK z@1_a*dU$jc z<0II@BL=ELeL95_5Rb;<)5}E45(MF4U^NSbae;|NIKS>TOJ3V3v_GGIRpc=92)Wo! z_m899m6^sMT#+p49(L*oc?H6b!2Cdqmju6S?wD^}MPrY=k5?H2 ze`5UJhTJwdrVTmj#_Sx%FRC0Ncs00l(L9p+i5X978W8ya2hJc6SkrEB! zImAQ2y!ckjffb*MRe4as91m#?_J!D(5gddcCBjf4Ji!aVSK&cyc$E<-YzW=}>8Bqb zE!cxXnA*p3+^T%PH zPxQW~&f|3j^}*q}7#fk@(l)7<4D}0u3wR3yT!BPe7?=W#^9k1W2z*PTVn=|hoYt00 zK7rRJyh5Sc@R}C*N;sk8ze^sV@kMw7FCz2wBVFaZFdAvo8?)_|005<<6J|!xEzn}2 z+0fI{WQOQ?1y#bQ(5UAZL{gmw(HMk@^qB2Kv;RDzEHMM$D!SwZ=Gz@P`vp0UPTOy)KC6biLseI&do z?wnDnrN@{hs#NlKPdie2-k*POEh@(>?7mzKJO@aC!Lv)Qg+TWZeuJk_ zktTPXDYMgf*vgC~hAia(ixy z-4RCPC%ZjSQCOJi6**p}llwS3-ErEsZU_idF9w9q^kMqg*x1!?CO*PFeJVeGjNRJa zjv$_OD6Xk_2RL5Eci37JA6filc3}eK`pTu1@W}AvBuU1`v<>*zSC&$KP5Dx1n`z`P zS8DXULx1O9SJ%32ze$`3IC&=X@}W`U3{wTe9Yj{hqAc(%a%GGLd3c`${;{njhDw`X zPo8h)B{DNIhLiw(wRi~}?j+^ieF&*OY{J0ZV6esiYNyU+&T}2oS5aJ?eX!g;0`!Ig zlM~3V5tin>5lnW{ALLmYBcG0Rjk7mOreTi?FdAL=BTl3mtsTbanBUCa@f5?y$0IQc zBAK)@>su|%5MNd5D}deyw*k6VY*2*)C!^_Yp*?<<25d3#Y*-VqrkxEBjmP?9WBCjL zgi8DY$f)~#J>7X@zXNxI`f{)JP9(!~?ZYS#@^nOHqRH^9_qwUvDG$l8W3)jmUE%&gcBy5{EjZUBnO z2m7cSHYsnQr6mnC-84!~#k|xtG$et=CJZA3iK@BMLH9rDwF|IKFrX>k++m^`rUB?X9;m6gHCD$`>ByW1&n9?;c#Qgmg| z1*&nUSz0X4Ng8B&htVT|pbU&q&kTPw)K01xI3bX?Uq-@!CyN3kjC6>OzqX9g0XXl& z`6-6-M{9i#Yeh2LrmRyTTa_I+4{IHsKPKOmnD6*pYA<@8p{Yaxxx@rVJRH2GUa}%+ z!YkI8k)Xt?_EuWD3{7rUMEc|li9hx&PuK%A<(vnEUv7$oZtL5k;?D4%&-Gy3@46_3k|L1|1iF3kznl=BB{&qkb+^!oImfRqb7F@KX-)6SUw9-WGVFhl zMe2Sq>@WWfEG#>b-MT70{dS}-ko3$l9RlAn5NoT1UwRT?`f@tuL%UWm$pu~UfJw4A zxb=SmrsE^L9&j)B4wlbu$CIB%QJ>LNnjb*Xyb^*Squ6gB4`v&060** zaMA!6CUc!;_{Ez8@mMbdS;Qx48Y2CD4x*18VQ4^%MKS(I@&Ns?8HO@G$9r*G22zOO zI?OqL7A#&Dhj-yrf@Dbd873wBCrB;5pnVU~5E0_`XYU8YPGXrVUv>`LTG7iBR=EgX znrt2-_6gQKh8j#zkP#4{G1rRA968#+GqN?|EU2 z^T*)^@Ds-4D*I2+@kNALrnl#MP2fRK2_^(Gm>L|@T_bZrV9F=9B!{2Lfqtt0H%t!?bBC#M=T$Rfn1}KG6ZAY?rH?xX%Dsphz30pq?&2* z(&6cM$n{9Q>YD$a_FEwPF!f=2}ktaGH z5B-0Wf}PtHwiPkNxw`^HN?fo4R5;Gf4vsAulX4Y*mpT0i=f+_|CM?prlNgWx>@X!X zSTA@~SjN8E%$<@#<{j+<7+j`Jai)r=9Oy{1DAf-kaG%2itSNz9f8Q8|?g82o47R)Q zmaxsWr_c4r_=6A#56-@7J}Y!z8mBDSqJXff=BWxxXOuoI8@~l%0BuUXV!0rKoEd?G zN6@hG97vxDZHNJnG!a3G`sqEeF&`e#-Yc@fJWF5*7?8wMF$A%*u;c(;Y7TsxYF=)1 z9^;Ln$}o(cFK+^tNS`FzRnB;vZ?z5f;Q?DnyX*c@f(0SS^{(r_T2j?}SY+3Ul!1Y* zExRwkPq@lxby}78I{eOLy}!_|JNKO_tI6nt2X!0epvSP(EG^!9xCSgp$KE5V|9%yf zrEO3P71NzG;4CI(U>qx0wP|E2os@&VP5XU*$cmuRZ70XMUqI&J;AczaAW5LDxij-6 zu|O%KgbHBRWF$=7{s6l4lEBZl{@Pc4$)iB4e#fW{?1!$WJ87ByAdJ=RWt4Mg!AN9* zR_4k>g1A16Vjz_%DJ?B-X)(Omb)wOe^dWdRd3Y0_v6F^UEY!Xs$Xw3X&W6x`Dk{Xc z%diw8s{+fljCndM2U9Y|;EhX_!B!3wG~Y{7eVhwF-5d#Px!M&><*_D+p%zZv{jyFi zudF;#mNTL69Ft6n^V-zZRHZBQC`V3_U2Z5n)EOx0O=$60;I(DRvtg-V7Oj1p-A8dt zglJ`SZJgdTzQ_1g2RS)8fsn+5N`(HrskSwqI@C`NMVh&l|3>+<1V zmWHU}0Tumf3ysfyd$rXNeZx_|pmk?BekVzv?C@%e5q0(XTSf(|U!Tm8ZQAPSpBsb5 zMjCU+u^C+gPh^^vYvH@UK6j7SYitPX5kGzUl(NcE`*YXX~ZT2n}s**_+FgT!31l3(#L|XX+CsqmGuMIERZ=bCcHS4ffjs5!a1D zn!3Rs_tMocnMafq74rc~Lsq}g(TnGI*GDyVeUkE^(gi{}8e8uKlzC_Gt{Ej9su7*w zc#Pzg4ni#EK>V)`gwta8jWB~F1iot}Z$F8oXX0^~3fRTtoLdSStnBtT<0JG{DZJ|mO6Tgv&5Po!~HLZ>m zd;))EeH~pxQjkfRg4*Xh*qCchJbr zOYx;|-P^bYFFuI#&AKG=q9?`B2~GHo6ey;?jD0tf$52CBq!p17&u<>$(C#eaO;G!r zV2Gj9qy=Pf0VRM=N&umWtRJY2K*2V6mG4pP>3SYdcLpW?Fotld!#=%5;pEb7@eyD* zXwXo99u*b!tj=e?#n*bJCm|_4-NxM6vG@nj;wsYvyovsPt{Re(5E=kPYozCbccHxy z<*%Y?>RM6F5Cz5ZbNie->}1_J7W&LDA2FDPslyy;1f=plZWj(a!!!{wQnD3H~q3`;d!h3 zK}AHN%Psv0TVmI&m@D z)@nEV-1djnoH@Jx65^cQ#4qB*r(4&m^^YI96>ZiA<-d$S@2>^cA7HNa+(_-xS^`TiTgkx(@9J<)7C@`&@ieW|5I!Zm+K z%bImP@IvSAh~>3T@o%vY*5$?nineRSH?u|cPkeOlcmsXcTb_%7%psA)gQ*nr{iN;p z^4bmnSv=gDYiqTpwFHxSfkTiz5#5Dt39B7&glN1Mr3SwW;&OSs zS?#`O*xSR!P1A`xg~i#la8|nW5@*i^Kp2te;YCk7RyWvSrrgGw7L*U>OHEdFVEjqgQ;V}>*6C5C&)|Nvkq;MP`k0{HO`RM)% zv)X>Lrkdsvp+IQCJjWg}d)he4cTYMBZT&G4f+5QWYWlyAK?nHt*9`_z>UM~I6`qQL z3fN|K`Aaqn(-)M-GaRx`_lFbVJN%6JK2rd^Y(>@!SPY6dFI4*lEdiEQoPP7c&N~u! z@K^v04-CX6#yMQWRkNnf50&skey?yVg7ip7noZjM2@)VF?|Aji@UI{b9D5-lMQxb( zC6S5{KRXGy1WO)Xx1|5Y$$Zd<_zAsd#PaRwM$Q?(w)V-(Oj@0_8Atw z)+GWp2UqoCAnsr99w1xBmoXgf6r7)A8kQUdA@a;%-@RXQ#1Qbsd3_jY)D>4&#seft z866#+>Q{g(tap6qc%DgWiTfb~^|c2<<$Gn#goaN>z$=Cm14;riRuC4( z=(DXVcTaPUyQhhhn7x?HdiDV2i7xa>_P_rAQ3hb9SHu~@-n@z#&o!<2nlxa2HfwFD zNtwYr#8=p)=rCG4G%g_j(yd?kOYgF@!{|(Ww%Ydl`_(E0bLEkw!}*=^mX3~9Ysq`lb)1f51uJHA zq6V{OO+6o==HGH<@jEtC^H~v;(JhvwuHB)TU&4?({YjGYuf&*yApE%fyW)>9Mm7Wq zCRkIoJC@4&e1CCDv~V$eO)6Fx!v?f^B|E3C^nmYjn+?oW9hN})kHTna(kXSn+q~uQ z=mbAgxW77sbj}|l@>zAr&hXa)BYhKHc#t$wpPC$novo9qf9Rry^;|Mli|HkZmfQw*}R zgzr+aLMqxYDptF0H>3|AS_~k$hxbXWO&pv6Y(jx&mtD4Y7c2d;O2qpPpoMJssI7GJ zj2F)pkJCJ3q%V<+4e`fEAe%;j1^+LM#KZCyP0dDzy-@j@7y^|m2*NEWus**%Q{m=y z+26`Hy{ZR_c`=;7;~^m9HrLxPZ)Td((IZ1-Nib101>VA&!stJI`-<8XiW5szmn#T> zXwl~jf14LwKv6kgOMOTsT$lO2#;PoM=!NNGfIJrSE}cx|Y(hpTP$<}hqf-3+Dnq*V zKuScSf)x@_LMcZ(WzFW95iu-f=HQ1>M(iw`0Hd8l9u&I!q&|NenC!M9gCp#V=VoBo z$v1Ih1d9woJf~xVlIdbQ1y%wV?r6d%cilqZB@)Q~E|7q8Fa~KwW%XcBx>w>vA7Ety z?{%ceL_^C>k=-YkLHek+5b!?|HQ+WzX(Q5KfyM+DD%>$7y>Yp@JHU~o0(*uq-dN;T zOj4?2BUz@@;9&?&@l~aJ=3Fa%Jmu%*jqy7wIw_474)V0lF-A?ul>+YzqCs1Sj&5d*^j$OXx2W`eOvzNQEB$?`O1N|+&T1~QHHgLo9CHwc6pU6H z?6j+!Cs^NOaj>7nu)k>zu4?KbcZMAhtgkBc)K*D=o&`{l{6)_+VW5@ZlA=J>;ElR? z^A=5{`2}DkQVsEOA6)wY$fAAYN0-oRwA`~pO|yYXK(vKuDcMJ~*-&VwEQ|#)#Fn{G z05}tlb23Q42$dpjGCx92&V&VPWLYGK>rg<(*z<|SVXMy1&VPr;Szto zEe9TRy{r^qk>tU^O329$bF@hW=%6ryq5~b6sc(V= zn@X4=3eFZUFz4&Pa~$wIl`?XxXB_q~17Fbl8MDpiPK315E<6KQwe++@_G(xkIWC8Q zUk3`aA^uhX?Ksk|q9K@kKG9#>He(su23T%W*fM`Nr3$Aq2xHWTf4I9i?zIKTs0@&}e_zyNoK#sIWTq|jD6o=%N5&S3up z)yR?+eM>;fR7h4(wcnP{1UiCId6M%Nj37}fmM$4|XP5*`oCK;ZNjtvJJ%hd69wC4r zk%Nr!N5Fj|qIi=&P?q`tZ7nf^c1t0k)r>9n3<69U6iZHo0_7n3@x^);>#VtTbHwW6 znB*WZ)Tv=?VsD%7IiSM>Y@SSr7n(53r}CuP=a%WayLm538T=xp{!%K@Y-6 zW2geN?TY=To zuOMF`rB-0=&J?&&qoJcC=jz(Bqjmjs+SlINdT4i1D{-1YbK;=0hC$xCobaTslr1nz z*lEOCPl}+N^JC#HWUW>bxUDm+%k}glw1FTc2e`CjQdp>^{v=!hoxUvTs+Uvbtz&SZ zo?+9MFFMPXo?Vz4%+1g$V4Jxv*S3RxeNqSz%f|Z1h)4oo0#x*2WSRlg?EpQPgXQsK zy4sn_nl)rs+%nTM^ke-J=5#O}>IE;))HnkLIN&{0-*rqgDI-G>xQdXqnH@%*gG850 z?Wd}(M{C}^L~@iR#s^3=S|$SX0p+VID#jE9-(!6kcNEl8QWE*xg^O>bu%RW*?X~)> zU5O2MA%Ab}wWznnItGL%M$5K1dxmKK+vF_-?32=)aD87z@`FU# zd;2!nS~sBc_e^TpC%8Ps9IPGf)cHl==DBo6@*Dex7q92*!XzrCyMIYNZ}hvd1Y&Se?gV8 zIcdkt)H`dr?YD$%d_OgGyEHt>9J2x$Yk3Yj> z2d^Y%_Ra^hB{n|~h)nO*3vmG#aCUBg&ows+5gmW7uYZ4dcnDC*uixtPz8$|XJid3l z`6%ffVl1?eAHg*AwolSq_pnEjujVrEN%shmhVpN6lc-v9@*=L}oE-X{RVGHMm~QO( z5vB7xCW>F=B!l~P`aN4$q#M7|dc{d39zi=LVKI03op=D5rR*lG7gt#c=zB4i>(
)P|#{O#XaPl=(YZ zA-7l9yUVjCTS-nk6G&#GS&cJ}-@GR-0&cxbHM&>Kt3GQz42mSScsKcpb>zj?L9!(G>?zLRS#y;O-?i=gxEqD<{$_d;MDns-# zrC{+^#2h8F>E!5W_Gx3kP5gBzmo)oW@#@815_{6J%}(QMe{In0C-W%NGCHC1K!Z+l za{8&NrHI*X8`~vP8cntJ{tQL;5QF);87C}c%f$t_EMa{lu9#)xELO_KH)C52H@Qvi z6PI+xV0ty?+xgkTgKqd+hNKS-71`Wk!hTgQYX8Gmc23FI%O)v3Kp9&u%m!R{Fd(zUBM02&iX zBE5|vuo(MZ=Yq-a%;UkNU$ilnt{u>ls^CsCzi#@Eg!1)t+%jLb0rnQr;Ic`N`|I%W zEDqp?+GC`JUI3A|>USW+f$?oB&~&N5uVletcHR#2L>M;u*w9n}m{#1Js7rJnw^CUceE4cchMPm& z#BLZI`DN6wG4l6=iYFW%SD`j7jP(nBpT*_y!b9Nik{@V^=wPhAJDL2>NK(9@Ux;5h zoJd>v=@m!K$Vc${F)0vr$cbs@dZ{aZ}AXyq%AXaw18C-39lZ1Qyt}Gzdy@xXr zO3j$*BlbX~GlP)&{Anj+VQ5=_S+x6p@g@zwdr_?mxe^eujs#?ktk7^{c;lo z^w!;m!MeysxO@*dVDitLxi&At#uq9 z3UF&R+=SQ%DhsrKu%NO(6q&F{lpWmv#8eU}x+Z@p+$UkXEUYrNc)p1atsA@=d?mXM zZyiTY4lD$MC8o%qUR#cinQ*FA@A`pF0c*f14@d2N@3pq&t&#@EDvU&zgg+x)KOt+8oxB>T3-za>T9>v413rmFg9n%-!xr8eJ6Ei&&C} z&-YNkTA^_X9aP#y{(43%aSMcKGYYiSoQ|mZ;7ISDW>O`vjD;>}2e|tQkfjEv?OU6T z+}(dqiUFaVsW!CXu^bk&{=VH?mTx)q0rf&czL?Bx2o`xt-%qJ9A;YD^rxBl?hL>!G z_m-6Ja&H|y2aS|2{OY<*o4*dVYuw>Ad1D3e5*eA8MgWuvHXw4vModIB2*lxCmdq!3 z!dmT(@?~VvfoS*7vO|t0@=h#-XTmV$DGmNwy*|!|=B6^lDyTA$1R6qwfmTNUN*AiWof77}_H#(9`whCmAv%dl3Y4HjB@<|3l&7Pz5Ucn( z&^>2u7D2&G0mLeqy-#cq(~nJj+z7&DLFS~_WjdcX`UUs$X#B4eSXfx9 z9bR-(I^WOF*Qq=S%y`H*q!sxd5&!nZC*>m}B~`b~7vZC6(7igbz;V()6#1$_B8Qfxg%3RDT?`Hx7yLYnlz>yRgrt|L8yo`O zf<%}`N6_Nop3$MNtlgBk1E*DrSn^5U3;J(iykbD@xO70!?Mux-N^i`f&?L3-ycZ>XK=5v8F2qij>NTC8hILoq=71e`Oz4|?)9zYO48PNd0OqsBm36xcOHxQ zVHQjK3c%&EDf{(@)`H+Kc*WVpJC4^tGpI7u?ReoK3d&!93jo$cojxN-9*J=TD^DPP zt(pRxVMDZ75!8o`zH3(ZCP{!K1r4IhXZajr7`Xj1AFU}J(4r}dmx%j58(&FPC;eRu#q;R+)uwvQdUw4;TuwB}y4cKci;COMvi@vCbtY zE(iB#2=Zy;cT+yIEJzN@wbRUjWh!*mQ65em4<{c(mVr_t0zcX%2VThHw>f zZ9O}toTh%s)n}57CeM8MxU+!ren$LN8kE!Rr@Y9=NxC1GwJ70>*!g zc=zs|R?YT%uluAzTJ3_;78^yEXl`2g&DSAw1D>VF%C$?K537iih*fYN0(?8_U{%W; zsL<0mD|^(T@k3wX@b0MeV}>?xL)_DrE~U1&0*ESiiCGGMJYl9Xfcvq)3Dp*R_01>8 z>p-o$eaSqsJ1B8G`3Pp*oRD|0c@He)YCbOR%g`imP~f2rTNSRxAyJ7S&mzHtT7$|< zN``Cac&7US!bwfJKwq-ExB!5bfVa@Bv*VO#+s+6Y=LDL)h=+~^Bt`9_=B07bT}Gyt zE(&C67f8ULN%40`o-X-`yx*uZ4SzH+C_{p)EIuzpZ-}ePmX0Vp|CWng*8^GRQlfGsg_nl+Q#SdF#} zmeH9C+rNrmd?mW^7P_S#HK=f92fFtQQW+Sul~zJj#KJQw;-{0R_ zpKBHifvfhl_qoqFj^lLfYE><~rw+vscu17|n*i?%BmKh4XiP6wmRFR}%lU@}(SSg) z35)cpSnzhQbzlT}%_E{;PqyLvk1V0ZBGOf5_CQ5ZhXA}VDnyk?WdI+iBDYN z*ojY85HF}y;RGp^a&B;)Zjq}Y$N}zYr+mpUDtY8N`-9JIaohC9)^cG# zeS0;Pb`*ukWkSAtEnU**Kws+4S$)caDRUq}z3Io!fX;t(JYzB$PQYNbHk5ocYg;qk z=6SfW8Y|(pMx|~I>7!?KUqHE++zt}8#i`NGGakWc&u-G&gXC;qF~La?A|&cGamxw; zl;YI~$BJuUobPs5%zHv6Ks8t`jK-_6P1$#I7lTR53j1FhNcYC&o|JAU0~*~dwckwSm(7x!hVW9U&r(ql+;Je054J| zX!IUB@FJZ)kCJfcjD?UH%io@TlBoC0c?m&IQJf|pY+sE7t<4|@0Dom#d#x%IhR~)S z2jFUQ`uje9#n;ynmo>pgq5sGa}@bQ!k;yjw+Dk-0bE?PSHZrGU+pl+IJ- z=w3}9>?hTj%m-Xg!Gq;NBy>!>P*KRlAum+dMV`g9~LJ_^rhyv^0 z!Sznmz2}>Gd+Y@FiPbYs15aPTe}sVoyc+6amJ}kyJA8(Mb|(WZ@^XMpHzxf$7?dKo zL$PQ`$L<#b9913$G1e&^%K|CCh5~o)Gcr-+H3KkD@%a<8VXV2TtwB>;JfNg-ur#

AQC1?Zn=lC?)5AV3^IagIx={HN> zBq38v1qSvXLM%V0NrEin!Moj@^7+RT;N_I=rUq0{1Ve44GW9SVpoEbD2y|$A(IE7LKCaw;tS%cLRYHY~lq1mWE#2t~g(=B8NeTxhMy z#;iLHlL-Y2yjQyTZTnaV)dke-jz;}&62I1%EQX!{GS(*Gjo3PWuBHVrq3fkx`WUx) zIYcrVi?djMtg5?b{2Z&fs1A|U%a>EcUjSfewA2g4xd0%f1VlmvV!j{udBiL28`{w9 z5p<~Fr`k2x9$9^|gC#vdY1+(F zcWV0zN4T=fcTy-alJh1jKIiGM!+?Z}nz>*MMuVFUZR>?n^O7ZO~(T7b;Ol$yUK{_1)&l?u+)wA0G z&%OH?p;3_Y9uZ{Vbj&Uw;^cZ-Cj&^aH6|aQ$-Y@-W23jqlNa+kf^%nQl)YH{ z=|1DP{c-ASIPw8Uj%9z5RIHMsq5vv814CEfj&f&dd3in%l-a%D*O*VC7ZNffvO;6F z{-CkkHF$1O%;zj}97Mtvonn2{tjao-F4(I(T>PrRR~3dX5KJ?>whe!MLLG9wd6 zw%Psee)hBeU{?svYrtl?mqLkZ7=}&wmGZHPYNJcG*aIZ)jf8|mLsmAV@xm8qG-3FD zfApIH^;f_F`vs_^2xjY05hIYUVV$WvP{xXsz67ajh}a3}^|w!$4z zP*Ct%jecZxb@m72jyQ0MoYeA0&K z6_)rZpo&is1>^zL>jQ|RW5=$!KuWvjI8gH`twWC>aHi@hqG!F^x{Sx>Nauph#ts!7 zO`P1W`B?o_YEJjy>TuL{`ARCp^y0Y2zp5(V!JqrKx~67k{%o!>WDz5;s7Tf}YS@0h zFLv8LDSlplt`Xy76ST}0bMBp2_+_c=kwT>$VZ>#6VPu5Pw#JggmoHz^IX(Ppwgzv7 z)r)GpPj>P?D=%X1vWlRBVx*uoKy)v_I6=`D1y-0Lq7*MHKQG(k#lV?M6oWlxO?ANe zFZ?3EN_Gjb2YiEQo7yv%RDGk0?!zE>dZxc( zo2NuZ(5n|e_2E5wF_n1`)CQVci;hvTq%sC}qi;C@E?pcJkmIV^KB>a%wY6i0woWS` zb);`-qw`-RUW2+aRXLm1)(?ZUnb(XujDlaSK%pv*@DlgXX9|~E*vjJcLe%j)xkB;C z4_zgnkJRI-R@rH|t>uErxMUXuvBMactf%oN%%2D7I$%9BBQqIoxZSCZC9u;7uJ^S7 zN6ow;`Ja+845CewR|1a#t(Hult= zWRH&KP(25F@x!u~(3tv;QCnZ&1s6Pbo}Lk_|G6hg+g4m+Can6l+XQe<-@~xU+N7{p zJtU_B%@X->^V@A6TScXOf0dtkiLp=5l#OQN{)7#^{q|ez$?opKF7s^xVC63(8DXR0 zKZlIY@gKaxj07~@S?8xG|pY%RM--*hs8@m3lbvOd3xA7w~xu4OU(?UYZ_XYFz|swY<)cE#^BXmsb`iPBZ&1rug@OwnyjE&t*Kj!UZzH)VP}q}|6ZD=j|DxL8HLQp zN?gA>vz3{)kXcod0L&a($oUHM9>a9!fMBQ+m}uZWDkW;|`^aNshSy^wqvAZ;uIQT@ zGQ;mo;BAuNK}#p?kRuYuI!d$+cT+ASaQSxi9vyu4b_hxa=iuHk8S zW_cTuAMPp^8%@iczddP)PMKoF-`2_eb1%6qR?9bDq*!8nZjr$EeZ1Am=&B@&YfChk zONoSEM&RX@Oe?pr7MF*fgpokm>mHDwN-@tFZGfRE3^B-yy4 zOia0{pnfb)@lc-SFu6GlBp^LAHh~0U_4Lx0~*E zx{1{72D~c>WDR}ltQ)1gJm5g_Yznuk-i&fQyL|{qXAm<2^hTk2;hc1fmuMwYb<6b%k(g&B)igL`90pWoh4b+QplZ;VZ-`J!^<-l`E1A)?GIsj)&0~WHKaHGrS6C^cpjwMB(nwJNHdbp$vh%CIZ}1sM zAd^yxy16-n-S@f>P-{VqiXOfxemU}w6yv*hzTa!j(STx`x8(4vcHp5GZrKpjiYeO= zh&Rv9>X%n8-&!WHeJIk&WYN)NO{1WoC{@GNFka`JtrdyTd0Cz3cV2~gqV+O#e2Q^C zgR2;@3sHW$cFq)Fg2dp{$|)Vq#;!9UX(s9WwUbZO}*Tf&e#*luMtukq00%m>U!B zDuL{8?~uMx77wfGyPy?84jRM8a#jCtCCt-Nhf{`TBZ|}|<`)mg!YU@LyS7MMmaHDE z&&i_>y}_=RWZMtBEKh)zqa=j4z;S3N@CnBUACie1v*u=!R!>X(A%g`Wu{Cd$TUQVI zQ>)|LzMo-074Zevk~BB!hPey2t6O6LdPG63kfXXe2&f9m0SpFnG_8dAa|t}*4LD#F zSdQ@l<%~sH@}CW)> z+_#dVuPY_cobeh_P|--J-os+H9amtBZJL(98v1^R$+4H<@fd)Al(R)t=+-=>_?ssJ zOitV%{gX zB)L*WQDfokz|`jNtD3`TyE+)I_c?b`!8mrfs`0<(@C8tnMsCE)$y(1O#@7t#kYi?v z#M07Ct5sEhFR%kgFu-#}oY?hNyv(ZCG(6(@(S}xwfO&00`Prk{FCgGa5gCY(LdOP$ zrnPeEw#wAKSNVOg%p<~2 z(y7_OsA*_6f0CD~x4>DGO6mi?2Ag=B*2#gHbxIVMLU|5d4<_jGT7|KBR>Z362lplj zm5OMXKGYbr&uy7 zi@O7LBd{7P2)5RIx?{ut)on>3JURDrGz@Cg|bM(Rvq=`Z)Kw9_;&#!kj)0Mr2=WM*8KjL>n#*CinQ_PcKHH5XAv@%6jD5UH*>GiG4(8M}|P@`-CBn>1ej z?7l^^C%TE^P9qaJABmQ0`e1In==Pcnz7%WA)u6A2QBZFLy?`HxUwFtG{>gO;LzdP9 zmkLR=RMBIXw1`n5YJ^~*OW1OO--Ao)Ea^!=vS%F)gGMkeOK8FiFFWcqYJm{PngQ-i((mMNH;9`QyzjrFM&dDJ8P zE=Qt8+un!;Yk+qUIm{)qESdAg96~sJ{viAa?9_cy0v{HC@V<2{DQ^mrqeLm=Yan%7Vz;i zHx%M{oM^ADXnE>?A_YqHftm117ogNVtpX#j| zaNxt#^pJou&gvFoE^2JzF?2`E;_Z;N-`oV5k@@f?1t{27(ex2dpB*d%l=mK~I*^15posK(i+Ifw98bjch_mUv zU$N4jbZw+=XwpHQDDO!xejK+43wv%`ZO_}s8}kd-bi7w%<>|nKk%NnkU7w$gjG6e^7E!iAf9-<(dO@w1bPUDfRKDR z99^Ewr(D`-pHs#f)t3xU&2SQf9xjHVS|7ky#~X%kQAhyxcdJjoVUjZMhmbZjP#!?q z+5f$7kc!T-`C$4zYuC9)8ghJ)ocXPG=F!$Jcu!_)iOK4!&9N@HXG);@6I6lFN7w=c zKGcSF%L1y}H2-cVX6Nw~L#mYl>Ph{{K|^L~ZBrv*L%4sYU}>vnqEGxK?XN&#?ckKt zxvtdrYnIGOQ?2J9I!ZRCT33c;-P6;?lFa{9Uf?%=8Y#6YrJYeU=jbUMa?o zR6)4!gExYzq(v6KK=%To$bdb$fUYW~fxmwr@&SEt5R)}(FEg(NfAMCdr@MEINUD?T z6hF*N*gj62=SnR2fzm)jp?(2!g2)>e8#lKx2PLW#!$g%kahwrPeOPE@l6ZdKPWQ=h zG8&i5J^}l)DuyfoGpB*Wz0W&p(g;1mfC)Dj-|A zLhD}0tspAe?jt)VaT02|+(bG!7|+Pq@aOTi{aLJx;g28tjEbDu$fRoE239auAkj(4 zS7aG8m;qk(NL3f%nYm6~5lf0`C0NF=s0fV~M-}+7?2LG;(DX8&sQ0At$}BSXNXM9I zq!J)eVNS!_l^@%+rwAW3$vD!3>HZi)9;=jcKCA#2&Q6Rsc3h-~vo}60MAzdsDb@>! zX&;ss^%C$it+eQn>jeL9SmywlB5PRRH*M;!2v$re4HLkg*uM*K!#}+}D{zq<>FSv* z(5=7jNg-8#)V%a|fu~ojcXUEcm&jBsSc&)IS$?r)zP~@EgIB4E1(dDa#eR3T0mzLfz58(v zG6%0yJ^<8>_4m5YKh!{Lp|qtX-EHR^!I96~@%c=aj^goZBa3I(I6~|chU*G zXj;CnGUx$L;bM)*Zb4SSM9)$6VE4D#)^lQLAS=2&mQ}jaZ#p?gARwfsLwM~>Fz%ol z4#iY`g^Emcb-hNL1OUK|_{TJPpJ=?_I63!yWgsZ7sBVb_B0WU=mD&|LB(7U=a-!NxynaW^4=?>5x;%IHj;Vr0xpIr7FVSa)j{q(qBaWQ42tt(@m7O! z`4768UDVy*zE~}i$r4VQiIQ5q)7!z!SP8T=Y*+Sl?D<&^*kLs zjNcrsthbOH{Wb~Zj8RS{RJX@v1&Ssb^R2pH;HGJJwkAX%J&{;}a=fEaiIk*tn+ycp zXHWt~Iau)`hviaRh|o(yEb&#D&S+9~xxQ?eYK>n=DW5RL-nF(WWOt$nOAJWDrxTXv z3qHT^)eC4W!W)S;$Q>R&YyT8Mcdja`dfV^(Ee-*Y;Xg%f9(`DPJA=-feLYDE;I}zpE$05Yoqyl$oA3^OK1NiDW9&c|L6LPGq-GZFam8%GFst6PX zO?Fl~Yf(J0cnJ!q0D%JzgIKq}8h#kG8|Kw#b&Wp1zdN|pqz!xHJzz1Z$Cv}E_0(UZ zO0_C0JT67Rtr;V?LMgUGnnq7236H#BdJrnZn8!Gc$I=ZFN?&`JDiIulyrUm)QmHrY zc)3|^f_pOM<8WT4bcFAFo`hf}67^r+YgSnYlZQuYSI0;tIE{g)Ez7m|Y?vN`ef+RO zkuTAQ5=Hn=Yj1WHE^0nWfak`5ekB0AYLjfjK1;8GbDjGi%3A)3b=zy_{U2YDuZ7O& z13pGAq97nlQ&bue8vAPXO3k0nn46Trm6Y!7H)zYt3V+;=*N7chA(lATgOe7msz7(S zbYcCyAFH(Gvr0$9Heh56;J|~cB^1Z>671uE&9}T(@;+a4qrW(ozh6ce|EfeZQmqN;1BN^Yyne8| z^f?30(%6@%xT_GaL%j)|tB-BrR+?}B4#5$!K3tXe`Dze*L$MbtK`^NQ&fdlPQY$40 z1f!+OLs5xACGqbqbbe`tqb4jb*~^echP}|_A2<24j~$HuDIcgUAdIL&2Wf}7UTm7q zd4$!vGTjKB{eFlX#5h(kJs~(H9x(wvj$VDxA!vuTsHi9R#>d_sXWz1QnHb2YzvXrqJDS3>U;_=j!(!VxK!U@c-sqdl!4cPC{#j5N_QUX zcOm2@__nEoNxyAD0&P(e?gX7kTq?ROjrHhvG93G=v#gy)B*dS&jC3&z;xtu5JE&i% zh-nOFv$A)x!IvNaXu1`0U-@QbXF6n~Me~@Uscrf1?Fb>A6D7^3X@mxnt>O1t4|R#! zKXN+R3rhPUPnAT_Ix9%06^x{<#{7=TT*A*gT>yCfDYhDqaPow?!41^N^F!CaaWJv(uQ3zswT#W{` zLL+3_+j81L0afk%k5??}7*uxx_hD_qsp*6JQiYGXC&+)CO6>5F-q^g-!2t7n$TrS{ zV-;wM{re;ah_|v&2*(qNj^8Pd`YIk>|2hfYr$mM=r0O+aAj!{FX{21bYGu4B!;9P5 zV^=1TbVoThSv^DnENoI#y0gN0<_guv!y(^R^{Y9B@gZNl*3k9vk_dV*z4~I|x|cR& zZCqlXgWmW`#9e5d>T?;c>jy!0JWQ*km2Rpr)~;oNQT+Rc`!wIXn9`mJaT1};7&ns4 zVcj(kUhP66o*f^1KX4spZaMYwJ$iH%W%bFvCc$wAaQo!g)H>++Cc*{piILuu!5}F4 zkRDKf;`!aF2l@+Df#0vH^(U=g&IV;I9}NcVGhrO}nM>?QQA#E^rGm?>sSHVjjP7U( zu*T2-zF?ZsJqmawMdp1xEdX-{G9@V7XTAb{;(D<-Ln4euyiuhT6)UG*T}^Q?3LkC{ z9@KoiyFWJvyt?gR#>5fQ>Q8Mm=kn<+~>aRuTK!wKF25EGk& zMP#G~jNzXaEMv@cqtf>^Q$`23EGNy`9}vF)m={VEhZXnidj8XT>SH1!?v>Yd79K0u z`9Hi**OoKhBPWie@wzU`(aG#!BmZ)n=% zAKBqImo;SgIJ*^U{;gNE-HM!?eexp=V6V9Q7D(va5zXB4f~MV}(Cej|+;3vpB+bOi zVV>2eX{L$^0^)mtHINhLHqJH-PqS9l1(ur_e>Co}0vn*uOYzzV#ESonO7qg|s56LZ z7r8UKMp=0NtFFabRyIL4hxLI}!36^JGA@YUO~-*jK5^Ht2)w?d6rR~%tEm@%hqj>O z{DPZBwmYb+RR9jQtQ`gvvMX`Y-H(5y6W3qYDcyXc&&TK{)L;D8LG0(91^8EoUg*Wo zX1c(}O7l@~wEZMUW{(R;hwzg&$j}PhE!#KP%rn!rRpJ;= zD5PpaKI&{a>L0l=?3y>spo;^hqJxT5X}7r8!{$t47PM-^pf6j z{*p3QZ0Wq7pVW#Wk0c)MJNX?+Hs0M&Eh$m*XwMbTfd0sxYs^D%>In9EDTKIa zEzy(bgw$H49C+oY13HhTf-T3LiiQjNAPDf5{Yn=J(+zczX{@HIIyPTNIFzvI- z;0fTqkb*Rzl=I>dW0s-{XoMJmKwgx1Y=CL@L6p^R+TNhT_!f>okxylGJHi8wZx<7K*?9;1R@P}ZmggWyLmDsK(s0+5KXKh| z{n+a3n!2Zx`TV7Xh92V&&lSB;_S6_*T9?U0g!5bG({uOE&KH5Pr`p_{;%lN1_FUZi z)294Y@S9hd1B32x#9c|Kc0fRTK#&`PBQ;IGWTWGcQtyl^RL2)Ge zk1&u!HlgBY z6TQBJAeTvochY&hNO;wD6y(Q^V>K$UG!WaF#ly1l7Ekk4AZ_ai%84WEJl4Y7-tmxb z8;8^H8ehL=4L1(Ze_VXjopMJjFq6*kgSz}jr}=8nt&W?oOuI$73=l#TpQyM(mRN34 zsL}J`kvArHoWa>kr)F=P_{ld9f9h^gUhfxGq{DoQKkJ2dt3ETLldx+}DiA_AkT;0F zHKIOk*T*tXg)HV^iFn*q$AE%IgK&LgI6b=KQu=4|bwaJhG72?DuL#5&mo;oDYY?M_ za>T)p&h>L8{gJhH__g~c5eFm@{p zJ1;)o*;IpkG{(yO_;cR--jj1!X+X*-YoW+gz!A?x)dG689u5Y%L-`8qn>W#0*zbjj zPFj}C9$P#J*S3LWy32ZMPCG(&N?`>R71iVR`m`35PB``6p@0@o%s`!%|L1z^qth;K zX@Q=f&xiDT_VMsWddNO<-*;DgVtc3N{k_?-R#3n)T<0xH$-K(jH?tNosh@9$ao`1lOzI)!rM`{qTbJ8xCv>iEf; z3(y>EQS{CT##BN#r~NjXN*FG0Jn*vkm8d>Lq+1Ps(coq!%idWW8q<4$cIO2s5KF^$ zLK7?888WOIo0;FZSEb& zqR)HpwaW^#a9AW>87bT(>33>r$lzU( z?MvVdMV?b~Ch|u;914m8zq4S=w|6}wEg3(I?X+=wlVwqQQdw_A8d*orCDoc{{If-U0Bk*T> zXADCf^eh(f>FHO)@zjSAbzVb8-_6$Sc4s3<>3?UrKfmpV|8AHuCK|qTH~nF>PZJ;6 z0jdHD8MmvQsolEg#>YTe^GfZ*@1C8Me1|J0WP{;BtObb+(1pdGFLd1sxF>4C(mVt2 zU4KSx_FV>!D+QeU-}U>BW<`gEx*-qGM*nuV-gWJvkqi0;Y1dW*;lR%Pt^hWvI_Q?*_#l!#fm^9?!wSg{IaVR7V-X)H3woo~F1 z%*^RS)&K?fmGH8k`x(%pb+apT8w;w+mW>KVml-vG5 zt7h_dJ(+e%|5!4Z?N;Y+*nK{Q8`tfKx?E^n zIb3t)a71xF+`A*5AY`BM5#l!h)vYmJ?`?KDrv1Z@xRY5D{gNfK4NM z4|F-|(>3w6V9x@{A-=FRpdu43`TB02m#ce6rO#S8X0(BV$uIAgh9r^$2Ck2LXl&9Z z2F|jsF3S#!PmO86#@c6T&YC1m_$TlciZHp^<8t?7OQm31;j^k5dYvy;k+>{oCS@rV z4&Pj}Att?4gFoMw@4N&{gP(02-;`ld^QSNULfpIR&m`1L1$}^<%~(hknpAF+Mpt1P~>pSVB(s78e_8RBy2fM{3#gXjO@b*oKQzmjfDt zey@E*vH9pwc!Y3C$%c^?6NC_aR+xx+hE2GH@*C8tul*G+7EuOsbd|sCZ29 zN<|ow31U=yU^TAsxFqFGb`+^edY7YnqxGa*)JY*q!neC|tmQ)I6R+7Z%P;COxdUv} zTZPi6Yhsw$5YpDll!I}8*TN}k?H3glXWJ2K@%@H5i{1U*GUbkxf^IUu>3V!b!gdBeh87|CnsXbX70OHL(82pDX0E+8nxPP^;Fj6PaHT3ZTM7%<9`nw zY1%)Sq!nGgb+P)S|L|%!K83E8)B1buF?LnO@j|w>+2gaXW8Ju?+@4vFGKoyG{Lgx5 zKVz(RnI*k_Ykg8`Q)A;u@7XI61Kk;-pf&L4>%}(v4Jkx?scH889MjQ^o~xI zzD3IG2}`qT!w+T@-;6 z$WKP5rg_%-kH2ht{2(DCllJwrUXl}?Rn*anb?aIz?uq|2`jcT%Ue3+f_4yi8nmCJu zq?bC8V;*Abdk~DgQHzJ!4wFd6R0XlH2HV(wZ>&jMt|TLFd9`za$J)AG^ggD!Op{;P zCa$9%s}CEAPK*KlK$HEZL}xsUgoK3XicD@Sj37GHA?IXkeq+(zOE~sa;0w}KNSBo zh4Ax1frY~G7!#XxYN5H}aJ&#biP?ipcCs!n!PJN~)?4v1A(#^gTR1p0#pVb~h$_lIP9S?)Y|9Z&rl4aS_yzJO}%QYG>Ev zbP-x!I+bqEhk00LDyBY%>2Zi|<>vx&b(R9MLB}+E2R(YbLN7ZGu9|2|=#l@NQ4RI+ z$I*_(1`{Q%8DlaIo^k`F#>%AXT}SVt)7`l+lj|*ANjqX2wGE)=0HinO`AF1_ei-PG zKS{9(k~8vL%h%Mr^W{mTqf%JD)82zSZw~z3UOzQ%%jYaUzaS(h)Acl$)S;$$Ir&ch zKr3RvfY-}rM?{T6U0cpIsRXgr0UK? z@pgkO=xu;0bbZmpl|WN#voSHx`pCF(XZ&W`$2zsdYS^j$wBk*zdQ3v$ygm51}ZGxnt>hK?@x)Y2Nay0HecN?cFzce|Wz zNvs8HTRi^ce^}9WmyKwf%x5Lfr_!xRQQr9wLkSuf{_eQ6SI!nC2f7fL&6*$eOhEsH zd5_g-$jr6c)aTEjQR@{CDwV@x2Q`;|J!=z{^}>zM1)IRjZ|WThvRHqb-zO6`oNl<5 zRJ?1LVJ3a9or2se`fohS&`q8uYDlt}XRQhDLbz)FV9#Jr0(9|ElGGQGYk;i9VY|+z zc^g|MJxVtreQkkh&?WKw#fywT?6D50ucHTc^+W5;%T=lOX4GaiD+)}K*$(5&$LqN! zVwW8=Ul6ND%QTy0Ee9;UXH+LaKiZumL-u}caoddq;La zmP=`oqY;L)LhPQQVc&ap5#A9Vn25+E+)sj1{MLPuSxi@;0)@a_0((Krk37thkk40(eA%y6vh#8IzDMt3(uJe6%j9|swtGC5@p}m*#0WjnzH9M zFa{7)W|L-w5W*pZ;igdSSDH!p zRBRpk=|TZu(k59l{Z=kI4kK*7fXAKgMT3DDOw z3|^%Wz<3f-FoF@l2Sktaz(O7yhw#DQUnFMIdn5fllK-1?n;aB~lLtue2a0-tZ^1e6 z*I~_UikUA^mpnyv@~q^NwbFxG~0t&LjU>?Tvaos%4Z#`;61F|Y3J>fBruPft%?O$lppc?AW@ zXo?4-fnV=qBYVWDka9_bil&OZU3~HTVt394S)k`zs`KU2*C3rIF(E;{$^J*_T8O$y zr5zAHj4wF|Y?S5Dt;T+UDY>T<*KGy6Rw8LT0|Fym_Eyzo?Vu3a2EcWCKtA`_{;+pZ z;2ek~22hK*c}Qg+1gG&?V1jPj(SRpZ)Mij{bJ}1#v$DV3llH(?XUd`$bu6-7El2#M zZ2RLZ{5MeGD?AkkN)Kt`zQRDVg?!?AybamO))H0X)k?2bB*8>09*t-?XpR$RbRKrZ zhh|T~YX)2JyR;LqpZ^v(mWe^n#f1ZGg(LGj#o+woNOr`UP}Pj9ZZUueCb?Iz+*BW5 zl;rqry^W3DmV7JZ%q-x5OV0PxtI2c40U0iXcp(H(PCTV(48qlgkAMqw-LHaID%Sxx z1zqK#v9i`@J42srOna?%<>GXYmZ&1XeL%A+v^R*HKjQ~z zFHS}D5TU&QR88|V3gj9$vTy7_BZ~nYA5rS~C*6I?)MHOxt+W02f#MLv=#J>(m}BRP z|Gm5L7vL#Yg^tPH|JTo;;3+f$MsySYS=`_&LX;>MN-e7+$ z-vL-bS`v7f|H;l%uoPkEQIEL6`^5$8N;2e6{tgop(-~kl*TgrydxyUWpaKtO3S1@E zgxq%Y-quKX9RYo=EjE6_t2^L|pYiwgsh?6%&E-$XWD3QmHil!Cn)rQi*(3yQsNIe> zC(&%~I1MGS-N)|Pi`k%|5p?WKasVVAWRvdxXiBu?r&*~MCiPYXfaH4*W$*io99ez+ zbWbDC959EM_Uv=o4u&|a7TSbJiho8Zx=6t^CEV&oQn`O>Nk5}nPv%v`KOo%y=Tp@ehv^60GF>n6o*o@XFU*& z=00H7*vZG{v#kO~gAJ!q%T&^6Mcp1yzqT674W#Q|_ob2mZKI{J{V6}i(mq1=F9O?D z#)A9nVKijNX}Y^PUj$}!d5qh>HSadgDYI#H0GEOS3E3Q7^gT6_PEyyjKZM?b#339K zhT%~a!$~dHC|eXE=uG2a7k{!pZp1X!;%d7J8p!Ye@z^RW!RE3ub}V-Rxtm47S1ZA4 znP6ehHUsGDXJ?+F>LG&~<96kc5fez3xNi%u-K>C#=*jONwo}zcb6SO&G?Yzd*m9}1KEoiszx1N8HMh0?Fx8%?tb|W zpi~F34$NSG+%_DV$dr8=GQtuWD&lkY?t|Ed*y;S89|^C(o0C_rcPpC~%1g()%3{6D z9wfVQ`jZBBY8&JX+Eh#a&pWH)%7=PAxZPNxHsaK*F%~n4Q3PgQM<-=(o&&AiXGja; zJhB2nG`KuN9oCb~FXtL@CZc?B8bQE|=J>U}YFakR)zXp~P$aF`hxyjJlPpFug?{aU zDG-;QZbLJ7YfTf8QrU=o#RMTNZ9&eL6xnn;ji=Aq-NhDW>DhC24;gk+%9;#C1o-n+}_V(2u?#X#?wK2Eqd( zjl_#fLH$bYJSt$spL4T6A419Yzwrc(ut%O(*0DSLulM;Z^U%gBfK*UMz|7XbGlB1< zs0D-z4IZ{pm9P5u)3hT>Y^TdZm{rmdVdCizY%TlCj_%9XGn7kD0V;_!@QwHTDIMT! zY=gvuwsp5a#F)7U>;!2CWrM6nMG~F|&nt~uNri)}UfPv>8WI;0`rh$#s|L4r8`ERX zI<|y;h`nh@e^wxIXs+5wziYX*ao{^!eJypm>J`lcnT#!bP~%r7@&Fv5Y^w(9yw&kXJ^|FfgUl>!|SZu0Kn@PeHkFK-!n@X#9tGz zIXLJ8GgSNU08P|?CrxHv`Rv)&AZv)am+DZsKB7J;$9H8a&r|wgYkqH3LrlwA?jYJV zn(Q*Y*w8qoy&hT)9S&kI0G-=Ik3=c>iaNDq3NSpnf>4Xum>=LxgQfw}%&#wdLBt<5 zUT>|*ueO1CR%Ot1anpYXX=TG&TX-x=u-Db?&Wg6Lg}51Y4w^Q8yO$R?Ty zg;UH(mFl9}wM%us8avvkAjrrf>O7+b^_5*}Q^mwq(> z@U&maKn7Ld{twT8JiosVBqP2l)DV@{v=ck~zh)YMhJl?%(1`vu9sjZX;rD^rC%F50 z2~3{<@%Ny>RmJ()#sA;k0Un5s+zbBSKH+0@;MP3@QpT*TtoM87+C2FQ2?-s61%KTA z0*>6{@rQv&(C@UR;Rc?3lu+-zB@S3J!Ju)GWIgiGawYH)T`~jB}A@9dAL2sJ2V0N#B2ftTu`zrT2gCDX2j`M1P;wl=gPXO?k9}$8>!$9IC7l?rHeu{(wFQJzOOs79@f4_ME zh6IM;0#KE&(eNIUjDZaV4$Q{P;pXiu|BI1lfQh^@|KYu-C-C#W)mdQS(u%)1*`3Rn z5x?4(*3{5OJ<4)l=mI-m?u^`H{}5YRUx3JV!?4$5X+tZsn3qgOcC;tPw5)3M$HmWHMAeI(-GyfxZakE(d>kft)p+D&WLF3ct96NU;|OGPiU!l*~ar zcWnh~L1V`1Mo(*sv-$Y>dA(0wwgHNK^JXms>GH=`5GU@{w5t!p(tbx4_4pE0#FoRL}pO0fQO4yQ7-bkEEEs*Qy-;WmJE!zLs9`KDJ@?~hx`;h Date: Mon, 11 Dec 2017 23:24:28 -0500 Subject: [PATCH 3/4] weibull not solving! --- absorption/weibull.Rmd | 65 ++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 65 insertions(+) create mode 100644 absorption/weibull.Rmd diff --git a/absorption/weibull.Rmd b/absorption/weibull.Rmd new file mode 100644 index 0000000..e543b79 --- /dev/null +++ b/absorption/weibull.Rmd @@ -0,0 +1,65 @@ +--- +title: "Weibull absorption" +output: github_document +--- + +```{r,echo = FALSE, message = FALSE} +knitr::opts_chunk$set(comment='.', + message = FALSE, warning = FALSE) +``` + +```{r} +library(mrgsolve) +library(tidyverse) +``` + +## Weibull-type absorption + +```{r} +code <- ' +$PARAM CL = 5, VC = 50, KA = 0.4, GUTBIOAV = 1, GUTLAGT=0, GAMA=4 +BIOAV = 1, LAGT = 0, MODE = 0, DUR2 = 2, RAT2 = 10 + +$CMT GUT CENT + +$MAIN +// Weibull parameters +double GAMA1 = GAMA; + +// Weibull function +double WB = 1 - exp(pow((-KA*TIME),GAMA1)); + +F_GUT = GUTBIOAV; +ALAG_GUT = GUTLAGT; + +F_CENT = BIOAV; +ALAG_CENT = LAGT; + +if(MODE==1) R_CENT = RAT2; +if(MODE==2) D_CENT = DUR2; + +$ODE + +dxdt_GUT = -WB*GUT; +dxdt_CENT = WB*GUT - (CL/VC)*CENT; + +$TABLE +double CP = CENT/VC; + +$CAPTURE CP +' + +mod <- mcode_cache("weibull", code) +``` + +```{r} +mod <- update(mod, end=72) +``` + +### First order events + +```{r} +ev3 <- ev(ID=1,amt = 10, cmt=1) +out<- mod %>% Req(GUT, CENT, CP) %>% mrgsim(ev3) +plot(out) +``` From 6577e039b6bcbda1ac5081475d521778b137f152 Mon Sep 17 00:00:00 2001 From: Vijay Ivaturi Date: Mon, 11 Dec 2017 23:26:54 -0500 Subject: [PATCH 4/4] update absorption readme --- absorption/README.md | 3 +++ absorption/mixedzeroandfirst.Rmd | 1 - absorption/weibull.Rmd | 6 +----- 3 files changed, 4 insertions(+), 6 deletions(-) diff --git a/absorption/README.md b/absorption/README.md index 42cce05..f8cfe16 100644 --- a/absorption/README.md +++ b/absorption/README.md @@ -1,4 +1,7 @@ # Absorption models [Sequential zero-/first-order absorption](sequential.md) +[Parallel-firstorder absorption](parallel-firstorder.md) +[Mixed first order abs and zero order infusion](mixedzeroandfirst.md) +[Weibull absorption](#) diff --git a/absorption/mixedzeroandfirst.Rmd b/absorption/mixedzeroandfirst.Rmd index ae8ccf9..31a12aa 100644 --- a/absorption/mixedzeroandfirst.Rmd +++ b/absorption/mixedzeroandfirst.Rmd @@ -55,6 +55,5 @@ mod <- update(mod, end=72) ev2 <- ev(ID=1,amt = 10, time=0, cmt=1) + ev(ID=1,time = 0, amt = 10, cmt=2, rate=-2, MODE=2) out <- mod %>% Req(GUT, CENT, CP) %>% mrgsim(ev2) plot(out) - ``` diff --git a/absorption/weibull.Rmd b/absorption/weibull.Rmd index e543b79..7c5a397 100644 --- a/absorption/weibull.Rmd +++ b/absorption/weibull.Rmd @@ -18,7 +18,7 @@ library(tidyverse) ```{r} code <- ' $PARAM CL = 5, VC = 50, KA = 0.4, GUTBIOAV = 1, GUTLAGT=0, GAMA=4 -BIOAV = 1, LAGT = 0, MODE = 0, DUR2 = 2, RAT2 = 10 +BIOAV = 1, LAGT = 0 $CMT GUT CENT @@ -35,11 +35,7 @@ ALAG_GUT = GUTLAGT; F_CENT = BIOAV; ALAG_CENT = LAGT; -if(MODE==1) R_CENT = RAT2; -if(MODE==2) D_CENT = DUR2; - $ODE - dxdt_GUT = -WB*GUT; dxdt_CENT = WB*GUT - (CL/VC)*CENT;