From a7b82c703de00b0e41b6f023ff9af361e08758ad Mon Sep 17 00:00:00 2001 From: Tobias Gerstenberg Date: Sun, 9 Jun 2024 14:55:18 -0700 Subject: [PATCH] code updates --- .gitignore | 1 + code/R/experiments/Experiment.Rmd | 2 +- .../optimal_intervention_model.Rmd | 5 ++--- data/aggregate/fixedint_pos.csv | 4 ++-- figures/plots/bars.pdf | Bin 13487 -> 13526 bytes figures/plots/scatter.pdf | Bin 13371 -> 13386 bytes 6 files changed, 6 insertions(+), 6 deletions(-) diff --git a/.gitignore b/.gitignore index dde3895..5938a89 100644 --- a/.gitignore +++ b/.gitignore @@ -1,2 +1,3 @@ .DS_Store *.pyc +.Rproj.user diff --git a/code/R/experiments/Experiment.Rmd b/code/R/experiments/Experiment.Rmd index 69ef7c2..c44c61b 100644 --- a/code/R/experiments/Experiment.Rmd +++ b/code/R/experiments/Experiment.Rmd @@ -498,7 +498,7 @@ df.exp5_summary = df.exp5 %>% normal_intervention_percentage = mean(normal_intervention == "normal") * 100, nopreference_intervention_percentage = mean(nopreference_intervention == "no preference") * 100 ) %>% - mutate(Outcome = "negative") %>% + mutate(Outcome = "positive") %>% mutate(Experiment = "fixedint") write_csv(df.exp5_summary, "../../../data/aggregate/fixedint_pos.csv") diff --git a/code/R/optimal_intervention_model/optimal_intervention_model.Rmd b/code/R/optimal_intervention_model/optimal_intervention_model.Rmd index c708a03..4f5f7f6 100644 --- a/code/R/optimal_intervention_model/optimal_intervention_model.Rmd +++ b/code/R/optimal_intervention_model/optimal_intervention_model.Rmd @@ -86,7 +86,6 @@ colnames(df.data) ### Model structure ```{r} - fun.success = function(p_abnormal, p_normal, causal_structure, outcome){ if (causal_structure == "conjunctive"){ p = p_abnormal * p_normal @@ -192,10 +191,10 @@ fun.fit_temperature = function(df_data, df_prediction, temperature){ pull(loss) } -fit.temperature = optim(par = 10, +fit.temperature = optim(par = 10, fn = fun.fit_temperature, method = "L-BFGS-B", - lower = 0, + lower = 0, upper = 100, df_data = df.intervention, df_prediction = df.choice) diff --git a/data/aggregate/fixedint_pos.csv b/data/aggregate/fixedint_pos.csv index 827a1dc..d473690 100644 --- a/data/aggregate/fixedint_pos.csv +++ b/data/aggregate/fixedint_pos.csv @@ -1,3 +1,3 @@ CausalStructure,abnormal_explanation_percentage,normal_explanation_percentage,nopreference_explanation_percentage,abnormal_intervention_percentage,normal_intervention_percentage,nopreference_intervention_percentage,Outcome,Experiment -Conjunctive,53.535353535353536,9.090909090909092,37.37373737373738,80.8080808080808,11.11111111111111,8.080808080808081,negative,fixedint -Disjunctive,27.27272727272727,18.181818181818183,54.54545454545454,65.65656565656566,12.121212121212121,22.22222222222222,negative,fixedint +Conjunctive,53.535353535353536,9.090909090909092,37.37373737373738,80.8080808080808,11.11111111111111,8.080808080808081,positive,fixedint +Disjunctive,27.27272727272727,18.181818181818183,54.54545454545454,65.65656565656566,12.121212121212121,22.22222222222222,positive,fixedint diff --git a/figures/plots/bars.pdf b/figures/plots/bars.pdf index e03091eedc750d4df98b7b4938b22460a92e2d58..c0c248a1bf4e6b4987888f4ee8c020d28db549fa 100644 GIT binary patch delta 9197 zcmZvBbwHF)*R~Q;Qj3VBbT>=)3QIQv(#_H#T{j{f0!xR&(kdw^-3S5#OG$SK(jgu1 z`ty08@B5zT+kfWXnK{=vb7tnuoNFrWf7mAqF$wXC@__}x{6cKFJZi4iau%L8w5)Ps z{Ji{N-ha0x8F8WuaPud{RHBIrCX3&vzdw;|djB!-Nd}GaapLQ0EE~?sbr@ssD0B9m z)bBw{4Id4k?ZL9srKDKvL(;w#g&Tu1f+p#`&b{<`VDV;ewGBAi-Md{inc^AvbPPx@ zw6z&|%d9u=@Oiz^yFRNEZpz)e#rJKiy?Y14v%9;!VwhgMyEWtcynB1L7bz{93^e04 zNp}K3+Z?9DW!v0{^rD-z#IJ;=`8xJ)?l3Pi^y{%R(S3;VSTI-97w#)T^U0guy_Qt2 ztL8oB{wr*zmpfNX`yM7|?i29g>B1T_AGM=FhlV|)nThc;i}`T)DU&Jd_65bBOiA$B z{8@GFe6BRawS-icuwP+uX1|p%!?E0ybpW^vn)`)m*mACC_oTMb@nHGOj`E_qTMThgN{wXm3_N&jsqU7Ya6rBpXsKXNSy&G{Yy&9GY0UzBbHm5GZkt2 zJHA_N!|bc?i>HXE#CNRXFhv3Z@$ZpS2~oqVqqKy$axTGWhhz5v!Yj=;Q9zvZ>vJHl ziUjDm3*x-h@+Oj|dI&^}By>JKqq|@`9x$=b<8`k~`;=gDR;AN85yal)(=>nPBiNCc zuX>;nO0QfO=6H{6%ek$y^f^(fIWdA#wJ9|gqqM$oKC4rTM^kR=(_Z#Ch{RK>jG0#D zo2HyS>`|OF*%q=b@10}!7QTFjX*v)?ua$&B+6Ak?seSPEIYoD6e5`e;2=_G|XTB$< ze-gbWAFKCk%nHjk>=y2zr$uJm8ORaB35f(mNoTYFZr(A8_Yf{Fq2$$)4*yX9DD1fr z!p8!>t;y$=P-H>PJ#mCcb`UZxDs|x=WOM8#|JckCXI$HunLm2uJZRSt<^#a5Zuhoq z(2mBoTBl#3F=b=D^6%7UCy?2-$u;+t3DD;EvcfHL9_T#I(4GPLY0;V%ps={ucy+Q5 zH{LjucnNrEzz^yK&DQ7_Ubt*W3*d3D{j8Bp0@Z<6J^c>Y@asaea8H9?_Yr)I8+@Lx zXlT%Sb#-tTk_J2xxb<`mBbhotk?I$QP~8X7PCN?K9T>?#?W_yAGW0h)FsS zD78T*5K3V_!0`4%YM%2t0B^pu@`jH_q@Jn|MrPtu=l_M$Huoy;X{xBLRGC-RvcLGU z5o`Qiye~*ZT>YebdYV{v!0@o9OIxvv0qmgm;il2{&2*AX7tl2Lu9i_k0=rU6NSu>t zX(7KkGtIj^`+WdBtgZZlp6qbQ*91Ic_Pb`fMMvRaBIxUZ?cqadCb9WgDYB@W(lp3rc%k)9y}=fP*(9@)PbQ z?(aHvNei}207?$7T@E&Kcdt=PVx_Q==u#MnY`Oo~d;NLLuZ__hxV^aU!d#W0p|Fmw zVQcvXXZnHqzYVfRehO_8g^_QE73oUk*dR5Lq7r(A6_ibywPUX$CO@!Yo<0!CLk&$= zI6du=2qVkt`h?y#&By9sqhopQmAC6sIRG__|H{tJw>%r)3{$+A!zGu)vSB$aa^C!F_y7d#iwb4bM(D;8c&>4)71rL9QMx>gl8PS zSu1ZYm#=VtR)rdc?^3tALXF}Fq{8}OA0^TUVINa}2)fozWOUP$hzZwxTT2WfowRx4 zEJ&ssh^l5<0$HNES8;Ei8;IFBub>t1Dtnly>G=<@piK7 zBf>x_wyLcNiXG z>~p{880dWycg+lh(=#!fryo0^>AMpiTX(Op4UW&@z|NP5l6RiS&=|wO%=s^mXB_Q< z#iERD>%2r_xUhudMM_4)y^SJ?HzRk4>vPO(>%xn=d-YSvT|>A9=Fm*bDzL&Ei5S$Tl95U2T*bTToSGIb@Jo|tjiv3ggJ&UCs&-C*t_?Ta(Y2t(f4qcBQn<&DhFf?-_ z*l(ATLA!vAcjJX3F;^raO%9I4*}dLS1~(dL1=C&Nb#)Ramn;R8Y)&|4bwer}P+98_H94M=iM#)Ng{m3UBmn$;@Vr_xkmIP8)Q zcr92nSq3r0OgA@b0H4IN@lkA+Zp>ytr*^5JB*`x>t_e)9 z6S-k9Yue($fo(J7g@$5GE!4NDNV{$u2X)TW825|!wSg4Ie`O2|ZFq*F&KAz}BXtkJ)nYrGi3qvlvIzy5)N~X# zEOB%CxLFZ?RuNOQ?k2LbZhFvZeq9m+4f*AT^> z8gRun6g1z^ffMotw;{0xy`#SRH6#BC6ldl>R$HRCy)YcxRxr3Vt%fHqKs;Gir8KNl z`m+Td2QMcOFWJ#EKd+d-3d(u+_^4`WR#aRyep_78^6roBw|{g~HYdqVy1mZTi#x<@ z>ZasZ#7Wge%@-9Pj&MV-xY=ZEcG$7}5JJ-Zm6C63D4uGkg1*vw6FZV2QYwlfjoB~~ zV3ucHEcrm@HgKu+qE=00kqQ3U{jHja3eh1Dx+a$Gws+qk^f+0FZXpOEoIri;2-}nb zPlxH2ie$U7Gb&g$57SzylN_L_wur97!jRPD>BsKZmA=Fhj1H6_;Fbbwrs8Nzl%E$=(C+SGm;QxTtKLRiScg*{rp$_(Vxd|b8_7V_ zO=|doW%|f2yZhG>_= zahhZfso7_})Oj^Uz;g5J)lzZOcy@4e@bpvmEHb-xe?KXPo-oHyniprd=}(Q z`-CjiVJseWvcE(&x`&|vLQn1=t_xoMW-%H1?;w*VglM8b;kA{VW*Fa2TR6-Sd{7fs zPW&o-9$pqHgqS-bmo~dUL)Jr$=*wf$?Ob#2?41njrCIX|GmDw=Y}dY=*jNtlnP2V* z3|8~iqmh^dfUABbj4|tI@w;@43ziJdCVlB$$l^fB1>x9xq4D(A#N&&Pk@JLRHN4BIyY*>X-%Uv;?6LD^ z;P~gHM%?h2YDzM?6BachEqSLl+k8b9SSX8+THDc($bx!(#kmt9r95Z#VfDP|g(h1+?Pmz!VqE zcv>h86py+&(UMXCJO=3uk@T#oA%2<#mBGO`vk7MVrZ=Jfub+wy=PuVbRlK`fjm(uY z8A>qi+H7^bmJa+qHrp25qJp=%pJ4hbSKrR|(e<($NIG`3^r8^mr8Y==I8x7=F^(Ts z$ADbLSN-b?$DW&O+pjoyR8#HH<`*v^WyCj1fV6_wJIAdg=k4XCrGQfTYyzx{)^;BY zAxk^%J2yXlAwQMJNV#Vl1{E8o&%QeHt#MBRl&*3e#n3(P!8*r|bSOuyR1HCg=-8BS zZqnjeQ^y6Sp#^SkIDJ)TuzZ;%}ln^6`2{jqu#G5iOnohIom zF5a#RE$^-h?`{dGht@-x`$Uv!aivYiFhInJxjf#FMyl8mQFo2GZ~SV5_wi48?ezPD z!q1m?SPwhAP<^{@+om++Y2FUsm>pbSb!Q)AV4W}}2>*Pagul9&!(tGB#lL_mMl7=z za(ZdbCqj|-??rNM2mSg1&meQxAri&rFI$6txN9$5QGcYdPh5w#a)nSL%|GJDvX4_p z7RYB&p#CWXCmRx9L#4=1HzZ{BOT=14d-dVhoNu>kJ)8Lh*Rnirs4uN1Q#6+zY0G}V zEHdbe^*>jf#+wb>=ZpbfW;@XCV+EgOu}3@IsH6nE{CII)4CwiVuQbkQ&+4*^8e;QA zVRKD1W_q}~aG5dMEGW4&Q{W;qZ4iMTg3&5HuN2Ki6gG_qDhp#36?V7=?k+-$mafZ< z9%h~}Cud3TIr@=gX`<`jrIG&jslF#yYHdLtqnr{j@ZsWG@A*G!jDDnSdKNUFcQ~N# ztP7tETNQyn8N<{hQt3%8`e$6~F7xPtlW_sPQLh0F;z&yS#z%G6xjY6!!)Dh#Jp}>E zmFrCW>gzU#Jf(+zFRQEc&11Eqm@*Gl8_34t^gM;(`mFC~Tk?98Yr%*(!g^LBajVon?Y_)j;c}M6&_$0f2C=^c9%gH0|yf)9Cxfb4(UNeoZIC^`fOQ-{+LZdjR)e^JBb-bo};P_>A1dU~EfJ&UC1F&Yr=8t$rAI z6yz2A>;Vb#ls11Z)@<-6SxZ01qsd*<>-=XXWj%OCb_V4r1+iv-6;NJ@6-3$=+yx1>K~0g0;c3G!Mhd^%iej(yLz`C7QeCD|!B?RcvJw{GR)XnE_#nYQr~?tl=cr z6jZJ}E8nBE8A=mZhY=Y29yQHw^5fneEe?Zt@@`EUP4t~^GTBf2j1&O|d98W>~aM$UFU{VSxT7NP^GWUod-z=M> ziTU1Vvyr*?9KGOPe!JGw3XS@7oWoHn_*uJF(RguWa`xyIA1}4v8|B+&3Xhfqfe1^7 zn_5$q<~(4k2^t80B=MrVusQ-nqu0gj{s~oKAP=W^j|=;=oW6?2yrs^|B-40|{qnu& z^AoPJz}dg~YENrV3y>^w%6U;g$f zjO}oHJ^p|)vv9cTbv(_Fe@a+P#Uy}WthruW^Cea0AILpCar;v)61-cLh(~qY_EmXF zr^sKu`eJ?RWuX9V-}ObvhjR?Y&@77d`|vPm0-K~F=*Jsq<|{_xVgpqYMNm8LTF5p; zpZ``u5u}__CkW}PNV833EaT|W1}8Ee1+WvsF$^|7#un3Y$l~^AOt?5K!wT|bhW}gXrng-u>aw;2pDC*mShZlZ&H?kWj>M< zP3-qhQ?v|ve=#$O2A9~|9>!EsKm8nnj{r+aHKeF} z^PrO=pb#qSE$jYQj9J69QAg}Hhud6C1+t~h;racJ^Iumc|A9kYJYBHh<|=&|j}kGi z-Vyj`ErKYM(clKRwzPeHf^zK-=AN%zq!cmcyhYPpqPXw;NwO5S58xmXfBhGmxdIQ( zE_uaR4}F*@$BG3vkQ!h@1xEtL4$<&QMqGH_f`208fg=OD5+@*#xK?P2H$^m)zm!{) zIKm`@>7v8;%a%I@U5Qa%hG<#PLF+(6hrLagbU|4|78RnX6fG|}$&QCS#1q3>sUIkY z2E|ss?Ol?E=Mih*lwwoh zr)vCC6-|^oCdc?^|DycABo;El6z)@mY_!FUqiy{cFAJZ2#W;1P@{|^*g81BF;mAUf z?^9MNV_F-zqK9@RrW1B*SqN4x-0hG!5ACPBVplgE${aX7ZPWpUyq9<)hEJ_aN8U5U z6EqaxL%>-5m%pJ_4E29Ha8F5dI2~o}_^KdbBG7##sYq`+lU)zbsoeg$`i|-na^g*= zobRq-YZu`!5qjHD&+-!UragOQtBQr}^OEXR@-$^QJ)SY?*;=#k;FyKR5O5cYYKm_k z!DZa(+tCO-rg+lmc^qGqs$#jhv0?#WU6WpOALL(wUcrZMa@e zkuZgz64p7e`sc5?rngh(HPuF$Z+g2_;KMm{G#+Cx)bZ(TSe=EyTB=o z*9G~fIgMjUzUZYUE4@~1#2VG0kN3Qag?03cjdL&tTmw}&n;~4kdJEVZ5QzNdhP+S* z6ut-!ipIu4;h%abWaK5Z?hPn=m!1}F7@7Ith91{fx<;spnT@;~cZ~J`mQQB16@{4~ z5u4@Z=*5#Di>@QeMG#)4X7OYir15=oS)Jh-((&fY+Q&a^wpviUUE$JYqc%p)$RxZn zXdZ9+AmQrLcB;_FV?um4mzkI$=Vdfgbz~GYhv z5r zJ1PT`#5cn<-RMo5g<}<1$x{qnit1tm=7+C zs?V4zmJZYegZo%zef&`g4ZUmMFI9ANIj*ziX>_VmTf}1ea|S z^AfHG{Yn6#tg2BW=Pw8LC7Zp4sbl=|1}u4{FvY zmrlOMdH{tavTs&iDRc>o9B;2xP`sVZW2~T&c_{cydZU@Qi==Z zkS;TBNI=mJtEVsi?WIQTsq7RVlamJ}T~@ZnBg*Yr&~_^sz{9xkb2o7qU`ryLcI?$6 z0ud19Z||p>b1RZq(h>zjU_=$gKHxN@ z0Cy%0qPmjpEyp`S2?#)PfGa@|#=P>&kU5MQ>_@=}_UC^;8T;-x`h&X@X+01w!Fe1q z&RZ$eXM7|%*>~e%7_>Ry9SK4drx|aQ7$nK!g^hnnLK<7Iq%Q4#!mmYus1Agf$ECAX z%+i6YdR8_Hnxn5}o#bj-;Ye4NR^>nUZ30Ax)F1&uw(`p+Kt$0z4i!-)rWja!h!FKd zFe+j>7`>x^5s|9i5DH)VZCPbTOVm#XPFp~* zUe1WBBk_Mu^TWSM%8w8qnp87+iqx@o@+%YF@}Gx;@6wXw)Y*H6X*jKtgviMLm+k9d z0B>yMW5-`q8{l|(YK7T8(;$$x*a#1t>XwZeJ= zP0GT=x>zAK2;n$-a7lm}GQPcMrYnU}=w1QCxO@k-a%^j*W|&RNQcM4NxbaDJI1|(5 ztNW{wuW=ggwngg=ul-YHRieWBr9Y!a>+4yJ`J6bbU8GY6RyCZyL`1<+jase#>%^Yu zL7pDDL=@w=m(Kxlr~>Rq0?{$VKg>3T_S)tEPgFHn4fWU~?G&Mq>UCX}nTT4hAbP}+ zhVub7KD$T@n|4@4vvJKc*a7Oac?K)%_QoExpC|KVXfIjsm{a@}m}2Y6$o(n=ldmA` z1bdKG_BaTp+}~-cieA>pRq6Uv6}iPS z5ek+|+rbf$k=z2Y;m;@EzWjD`I|w=@0vgZf8|R)t$mR`vE8CdDxFwa%T=5avGQU^- z9kJ)BHtnZW!jYstE|ew+6E~wHA?$)MWUCz8&JN-anbRooWb_ZoLvC_loM(zKV^=>h z;CQ|@oqW}hO4Sxn10TCO%TA@@;dRwNdcxLy7_Uro=;s#1$COU#jzWBU2@spF>Keez z^7V1g&6zU1ndmv=Vs$JKLD7cHOk|h)JnEhdn;ZM|s^En^6_Jyl*d#&KqIw?1te{5< zqAYLbene&?1qcjf$!6j#UZC*^o=F3~0e2hkHxeB0wT;>|N0I}=$L5`jL_J-!45{H% zWR#HIgit}kVyOZ|Eo@z_dnJ1@eBWdWJP%S|B+SsKeqGniwORwoqYVFsL}nuV?@<~|J!I$QITYG QFcb$YjLXC%uK~sVAEa)KVE_OC delta 9160 zcmZvCbwE^G_clW}IMUJ#0@58y!%)(VG|~bhNW&p!=xz|CBovVDl17k5LPA)olmbS3xb*uS~lB+9wH4*TWVkmJ%7V3?YH zdwFrzJzS0U=!bRO`=RR9j+}iD7nd6Xua5fLw_A9&w})Su<|=N^E7HGt-1;6*=Hv}Y zeupi_9ef8WgoB!IxP?zfJ)aN~^{@G+Sm>X5FQz0Tcx5|h5h_a?&HHPS8#h*n6~n?O zi@y;8wH^#*msB2{GyW${xp(+4tbeXUPTZ(hZnS5)n2)hz`-m|1$ zsc&vZI6Ksi-p?gFW|5rg|M5gvh*mb9_K4+-7NG*LURmwl)8Rp6wC%|Ptd|7~xm2QG z7INpSB$Kl{6hk-$;ucQPtBcXjXwUMVBwzEM;a7*eRolh5MFj1}y$ab~S+21?-$|HT zlK>!+>H%$R1SgRFX{|Hf-F;8xhGunZL;zx^=;3)XEAdrQE^LBx&utUjwx@q?MEZ`+ zKxPY*->DD`a^SRYD5@OAOV(&EyX<@2gC1gmU+PM7byq}Du zw~qt4-qT&gaPPuzo#`ypsiXqb(37R&=UK4AaoqZ@%C6*T?w%>mRMRy7_o+l^IyHxn zA$5!v{3X^ru?Pta3GCgi+`KU!2Oze-ng>1+ zn&P2692P-~)DO8U^U;`fhlxb@}<6 z1>3uodt!ADzXq6a_wOp^t}NQ98-Zy_VKc06Xi~v16*QNgIyO%a15=zzT~3be1M;mk z3t=^D1BzK$H1}$weKK9o7B7<{-jf#r$ zY)JATo;YXkj!mSp%6umGE8bra#LBMM3t7Oy12R^P`y2PPL5<0REu)>!&aOMVes$I6 zF{9>5ip8~A!AAEI3V_#g^k>pO*FJdU;ynAH{J-byUmj7`mMxnS8uf`;u#plyrxYM@ zCsZeS|FfNoWY0Z!rOrV)s0|eRu*;F;>u@16#z@a+FYmgq&HWB#{~H=*Ex!7wqcrcn zKW;XR&yE#{_75?kE-HK*l+I``-CgC9NaWj174vs^JzP}3g0rieBl0J{)vrexe|h{l zhRM}Dg)|;yv$dv2iEm^Qx_FZ8?@{(Kt!@-IYoav;1*=BoOHBj~;fzY-z>_z#vb`N? z>9+SDQ{T8X-~A^dC@bEpY5G+epBA-S(**q5)0{to`3u!x8S`=RjuFnHIcFAg)(+wo z8QRr(UiBN=2|ll;JL%DLbu?;)>shFT^$p9Ey3U0yBRfsbaMT*vv$>HO#VGA^=ecQL zs;FTYIGTjQ_1vU&J;{trP`M~=^1foXNS%@yfu@&`H_)gn3mV2AMbqOq8`=>7?G^?( zGcw)dC5zAO>JER8>VX58bjlBVwGPNKH)}A8T4h`pM`TJX5ucp%S_RjF^I!9Ca0m>M zD`j{l`pV*!V@`0|=B?-}6k4(4)WMQ zkI2;cB_t!uN4R#(;X|V*5@3KD?x3F5a1e>XCIoqZ(y#tgGBWw|$gr*DxKqeeuT1;> zSOyL4QvYBpl#J>W;)7RiBQ_bH0{e-=F;PE3jC5$&7VoKdN`7-TJhb$g39``89Al@4 zlWWE}ayGzAO8W$y0>a=_nYMfENLmpNL4wMkug;1($sIs(RBB^sM$^S6S)E@Pt!gad6gNYMhV26>S#X=zM0f%e)F@1L_SH)uyoN2?_T6iW%Ya!#^v`po+(6mS>o)x`@dpPx=d$j$&%s`{|Mi$M z-aSmLv$*Ny$NBbFI0uwHTsZdtqO+h)WqTbc*gRvzZ-^93dvfrBqhK7}S^xx#zpN*H zXth_!G$N!9*jI7c_B z|J4i0rm^B^ydMj$+~$fkQ=HF*YV}F;+(55s9)Z7F@>15qlv%H|>FRmD%g8)q2pkm9 z+4apywyQy_p3H_u+C9Ygha|$(?8b!k^T$&nj_^KnVuwITulmoKXug$UJj{@OJ_deN zd!X(=+m(t}Rn>x)gSrnCL&P3^(lI%BR`xPBP+rFNWBe7TKG-GMuJ(qiBVJy$IYCZa z{o(;fQd?|Ykdg&GW!b;H*qxRVp*s!dX5eIA)v95f2oWFgQcSZ}WF?`Hq%cLXlc3nO zIKaKz#FRU&K{gQwmM+U<<|@YHvaPk{7aVCvYOh=@fZ5jqjQMmpDtSvMPrd3-bjT-p z!Qzt&Rac=bf@`Emv}^rew6h`6P8PcO4ZHPJERXHVGfvyd(rIWdcm9xEo2RH;~HHuwn ze5uf5qUaV8MaEL+6;1p0MZ7jh~~8>?yQjG!39+`4X|a1Pr3|v~Zq`RILRN00fIctsq6~NPZN5 zN>Rc!IL4fFW+1uoY9d=P18j33R!2-LCmUqo(~Bp;kAfq#yM!Xg6GKgNY7C-c)OT&p zSmu<896>3A2U5DAFb0cI2D^#J`lVU+SekQNd^u_kRrun`cUTEfP!!kh7NHt?l} zSkOkj;DeF(4RF~y4kx};HaqPY%x6{n^{5ewN&g|sE6)g*)6xOvn$K1AR1|IOXzr7? zX3&z6&4B@7xrKS-7ZVNM}+b}l0{Rv zrmADUJ;c?W@~cV?ze(zA$#p6G+E((bux|KoReO}?(zXOSgqwrk=Z$Qq+%*q%f=12Y zFD;$i?@SJm&?!VhCk_dnSeJi6r?L>VwSNHxXILY3RNrP}_37qNy(~b#etZtEoo@#c zdRX5b9s4hdxX!|JKIi_y+no~EH59zK@cxi_#>XR>3K8J&DF*y}i6DUCXC8^XKLSj_HcDok|}0GawFU5G^_F7kg_zJk6G_M;jj79ldE;?C%I@ zy7mGPdpVvJK>lK4Kcuv21ANj) zxP?ebTep5vo-?63i5(E!1bwq#1_BTl(Ls{DPk)dy!Py z9xg4;jGmf3zjRMp71Y!7k@uFgbJP+~se_H{$z1W3!^Vkv-ja0bkO)~$L>wZbb$Kz5 z?z5)6bQ+6t3GhA!wtu)F1Cy8k6smS?D&3mUg`W2xb07f_X815KNs zGN*em6NeA_Clh~5V_o+VSLWGMaW^V};uH}H?iwOAw9gJ02V?e3%NM+=I*W%KlUvA8 zCqL>q%b|P(^!Lf;RoEJWYJV4IcsS@=^>=5bAk74}(Ler5`Ylb0btcXHn>qu2?OO7) zEvrHQJ)mFUZ1 zT}G3{8|Ss7-{yL8D#v!LrQwa9OPk{-Rn4Yn;{r=b8zpH+bNLUlFYJ6Qx!uj&KNoCA zT#pq3%i9W*_aA3@rTM&l9Qq{l?JU1VOWIv~L9odjSM(!+h=Ok?QDHKS3HT^{*Fm9u zIpXx>%{DsA4s)FFBG)-Sr8X-qXOaFTmLc!>Jed$u!G7T&SlYHVj(J3bj zN7_9VF#>?emC+5iZCO$=x3JY(oUdHz1|{dniuK_raOmV*g`#RMqRB&_vGcU)M@J&# zs4rZf_5Hl~ZLe|-hFuqC5{4akQ2e&;&#A(-#PwQf*6tIx2N3cT(oSqL@OrIn1M5g$ zDzgl77mN7DNn-FMH!YomcXsA7whsC22OzaQU^=awttxUL^K-JVZt2C9gPu2hwRy4h z8wp$2+aqY#Hgfk^C^hcUPBTTH_2ay_o_t)DK@31H^yRNy z^?&$Ru3kfldY})ii}U*S(At%7pVVd&$*RPwK^bNefW|=MeVPOdkl~+ zgz9Oj@SX*O@D1QtT$f%iKHP7Rov<~)?kGM>R^frmX`X$uCd-+psl|ZNj9~MN>$cF< zGU3@Or=`z$#b6?kBU{f}p{3J+x1BgvEJopibN*4|r_CVTy$VT>c(XLwm#mm8{~#s& z;ufMAU{uVHyTX-e_IVq;J;1MXL_{7xbj`xMic1?4#{tlTT_GFr1nsoM$`Eil>|s^! z8M#J(4rM`x}Lt%VqV|@Wr z%?6s4x!4hY%2hYU`YR$+e1@st$Qf@?F2M zD*z3x2E5Vxk9`q*+;7>4-)kMe-zrtJpFfh<`YTmXPlVtk1IuJKHwV2t>kbL_;87=` zG3t1Np5E}SHGN-*Lny@%tXb<8+sir zwV>l988t+<{V^!h*q+@tQ#lc?q)BAconeDrneMSHpqCyO7XAaa0(LcI$njdS1=@zyqU>dvK_Dke-*aThwUWK=AvHkZfGP8WRnRQVC@?-u z_S!cn-pYtGFnH2Ig<+V4N2YeVsM)?QkTQ``D(ryD+Q^CRUB=*$QR!=<^65L*7?%v4 zv0$?rc*K^n?=QrrP+?D~cfRAns2@0hNM!8&0JT?a+0oB1nPg&MkM5QLG#P=MAVTCh zKluGKWoV7MHy&JMra1~jd*(lpnk#&MPC}0!_Eyb3^E>!WG%C#C#ztzO3PTPBV?k}5 zAifZim5B;nMg^Hv4;ODO#h@nvmpI(j?3uXZP;U)>OJ1OA} zbhxnGX&-vXK@ugJ3kL;f*&irU41+&x@aq`p4C$r5GC6DX}*Q3Feh9CUbPP z)%Z{)*?)%~x#J;)?QFBSN)jfczWjeZE((6bo6~5t1FPBvy?26-^ z|B*a)!gXtQY~EWj9J-wjnB3nfI9yPd-%h`{+@(BpkB@chk-YcL%HvU<{$ZNjw#fFq zW27*tJm=oKwX7;2*(})oO8^5+!x2lcJN<}(Uox~uThKJmYmv!rJ|Ax)DBhPyJD#wW za%T=>TahQK_lJyW*)wgKmelnOlqKoaV+;S$=#QYb-}E#}ezB*(Uc<50I97*V1aq{} zNqOg1!o%eX0Zyca+F0zoadcFJoPLFOrz2bEYX%l&s5sdQ}v8jF26ZD4&fCCxOF9u zV~f^9&%%#3g9q<$v2=nYjqhl6ktv$0K#9<1Zg#^rV#%2_^%J$GF>Hkw6=(gq6BmP?j6bvs)evBC z-WEW8gHRSuP)H>qAYER7azyw`3QPrxgr2J^QIiN5ce0M3n;KtaP>(P%MLdy;&RLW# zgRe+h^4Z~#8l4wQj4gNnV5H+ZWe9?$a6SC0)T3NNu4Hf!CVHIe3<0??JQMkN0dDUS zO^k*Zy%o5EH%87nOZq7Z!1`FF^#S(~GetKRN~QJtmRGwxd?WWaHSa8^Xdo>gEO+M> zeO=K76MIPi5PAKe3cRBaMx|zTtymnXf-I{-#d*)cb$D$7L)6~*`lHrzv`jprO|o`+ z$J$I`G2m8JO39>!;J7ED1@Ih&3)zRYO3%zE6liPyuu!?3aR(s4o?A&o_`uu222k@0 z4NNUc39L(eI!>&;5n@YMzx7Mj0CmYJd2Pgj3Oy=zK$Sm&XK2CCDd`sSHwwYhSpw>=j@<^rz=$X8=FCox-B>JL$}I>fl>O#VLJz4@a{K z9MQG7y6hQw=E4jFNN$2(9_1YEq+Y{g(JgSf(jam%=pHKIPkp;3_g107-7#SIW6LXL zJ$VuZaKXYL(U9K$o%xrd?0H~ZqX4y2K9YHfa?aI#O7=}mrkVx8=e82<7F3fNn4EEyG6(4uhReR8mgj^ zp3ydUoXW46Zg?I%bSBUGB~Mb-SMIkq@s6DP`%O7|`)B z76P~B(HaN*>a1V35+EY_55yr6F2hn5J>owoHPfx6nsX@`#!J=0x^g8os)*y)8}9cX zfHyQKhFFq&y)sstX>iegvhDK^z8R?Kmav`{CP2QXD1C=~O|hpzK|Q{!L2&|&B7}Ed z+M~CUz>)j0M!MCEY_Lu4+YPpAOM`6zc0du^loS&W;7*pc;OIDvsXE~6ef^l^VXZPK z3)XN^c`YMK0zKB%LcU<~V8!3jd5;Fqi!j2%l}*dH&gDUFD!6XQm3lZl3S^oplbb@% zu+=>^pMHSpg_CG`0js8{(7^OkQr>#E&=u9-HO6UToc@Y68oBxp{^jBq79dGcBfSw) z?+KCQ+FdDPjoxqx?V>(?c^|#oyCfR&ieD35&QIgT9<|z|$oyz9tO`GXVOZT7l28+I zSi>@CkUgy1dq=D0A!*eE*65~ER$QwaRy`RSJrelbn3X&-i7Jtd6eG4sm9+>)`qdQ6 zmZ_J|G)b(J?fme{uK*b?D-Py}MfK<$f9*aA10qd|zfl|4m5vsb&$F3l$vg!|+9B{a zAc2d1xCk~E;Z9AE`^Xv@7@iZkch;RCubHeFjFyk2Yd421m3HHigyF$w0V5`!p7pYn^aF6ys-(mbfO9P0B>tshJRO!Tfe z=(OQ;kiV@)n|=O3v;YOmCus-tV5F|&*Bq|-FovXP;J|^M5Cw}$^nO+_2_}*3NxZ5` z5B%Yzltmd#)`AeA!?JEobug1&UYI<3K~YJOt?1S+V15|kD%O&$(Al2pD#J9n%GvbY}T(0Wcrp?)G3q4Z7Labv}M@MSD!ge_hyLS82nhJ%8 ziV1P$c@Iyv+3qiJrGZ>^V6Tj%?I%kSA2$-If^sE)fn0z&x*ohkTE-k5K!0`yin4H( zBU#xA?2Wzs5ZD`X`|Ez?tow#kP7c3$mA@TQrDmdTK#{`wMA6BbV?~q{ib`fy&O=ps z?hrE`*gI45#f^i%dSii(nkGlQh`s(<8oO0H`CuB9jR0sFhV6eLJ)KRHY{P&2Ka^VN_wDeaEZUYmT$3p&?S=uijqgRZ%)N$+9p@ ztnl8zttpRMLEq5f=jZ9``Jo>T9}%VBHxz#8$q&$6Dd|5M;(e|^pisfE z>PI2xHBjH4{Gz;)eH^3uG%eAoxAe58fA3L@Sha|Lb9S`5zYPD)Fc5Ga&QDW(e*MM_ z!SVwU8od1LY3`77{~J-AJt8{p+1zx|u|?2YY|j4DVAI9qV%U63+76WHitt7BkBfo} zzo_d3sU7OwiUgK>pMMv3T%X+xhY6c0Y`Oapnu$WL5kz{8Z*H-4rODB1>Y3t*B=~sf zdG3D5#|7NKE_g)wcm)3;6X56jr%XsNDTH4P5c*#eynF(@!vE8q7x^)az&~U{yh8t! z3H|TzynMp{A6;P~{=alp{^L9HTN~d0^%E%*;{C@UqP+ifT2xRpDNF#0!y|&r%q*`C G#r=N-hAegf diff --git a/figures/plots/scatter.pdf b/figures/plots/scatter.pdf index 157b58d34a9c9a597a3aa46639fb5abdc66de411..7abca1c80310277a2c30f21a18132d2ae614dd05 100644 GIT binary patch delta 8082 zcmZX22UJsAu(nhM6%_<7(n}d^`K$PAIK|nxyl@cH{DbhQl z2!twK0g-^x%M0lJ|F_=FS~=@v_SxC{?3ww#Z>GVv*EUuRBqk^+bVu}#umprsK*QNm z$=uzF4Xh+3EGT?O@c)g)A)q88{K9eJWyC@f3DQIYfc@8n#^R^Vd1BNguP>EQSax4vNwq5AK<>S>wJz!V$%Z57hL@(Y`6@^klCu zCXnP%9LQ)2I{xl_(xDsE1vKs-vzA;ndr$cH%F~ly)={e``vcbuj+gOkq@K|>lPm}z zdRu7l7#FQ7f8fTzdI|5h%%5pu#%jC_3@g zIYz!2p2{vCHIscmduXMJ(NA`_?Yq;ZE(^yWUd9S$LzN7TvSrJ^#?z0554`!=oB4*; zJWEFLjuhU{H|=NN`)BX&7;7|T486PPCv$CWxv$Tzl zgc)*KB92O?i_S!VyRPnpaFyNMEo#=$obOe8&#R?xid{%nb{}{Wk>Bu;N$#e=wc$~f z)LY$8j%}PGhh8j8-aY6lADZtrkER+CMh28f)+?6h0`2rIeL?OE^lp|FIJXj^NsU88 zx{q=b6rzoE_SD6r>42M z$*dTElmn*U0N!-s9RnmV^Efk>tZSYZQnU_*D?>}+=uJQ5L6Cdpj3kYi;irx1rAnPD z#sI&Er&|8~9QD=}ooRW^e0_6{bbhTu5#$PEDNkgKNGBh9-r1D+m|pJDtFVP$xKeqK z&m4Yt%3YQnqs%tshjLaB2%u zE<8Kk+WVX@rRWRm?)>nSEo?*D>WKSV&T@+)U3A8n-i#OlZ${SWSUgNJ>)gxl!5EA( z`!i{c++38^ZiMoOvxyF7Rd2W5vKgM(XO99ZH|FNlqbyiP>t=6PqcWDOq$wsHAKFko z{GLLhE8!C*QvWTd?%q&#=|@pj=hBbc1-w;@Unymu|IDqsCtf?Iyx{rR4H4d1A&|5E z-tJ95#{v4)EgRQn*QZRc+b7@5@Ox0MxL=O>Eta3C>{f66fom{1OVY1%x9@P>_FMGs z1TdVe2ooNl`tCMF0eQ8+(k4n9eeu`mNQ{<5FTu~5W{>%zOA}iIyQd#GY$knr0CkFd zqh$SnyTURv;nGFXwDvFMmR-YrY1XXnqKEEoBA36%D?63m;;VOOp7z*&Ra)$ccLJBC z(w8wI=k2K^@!Wp=lK9~`an)^TvPPa&AcSrKH10Z=58 zaUj9gZ7s$imy0g(6UPKT6FyEXjx?|vmit7%G$A8-T)dhPYOt?Gk-fsir9xEO7s2+5Qag3`gYG$vXf(y%+X+7Wk!ptUpM7eKouJFT z5YvSkRWz(#v_TyChVp>vGxhAvfJGgav^`{&n-VKPzH=3%c z#T-jLXDhpEJ}mo;rl7z%YG{A)X9YbQTno%tm;7SR6%T zyy4wGI0e(d$rgMmVU)3aTHJZ4qB8&5k6uc{mo;6)N?94#J%D0Cn*f*rn#%RJ95t)A zGiup<+6J(%x4&nbPGsHuig%g1UxKl& zRXT0Hx{MJgydq;PK-ZA5!7_TYoz-vhdArkhOm7=I?eG0tGThCoTCXu#Sq_WVYU{sq z);IY=pwb21@9IC$^LNj<;84u(T}hoU1H5% z@=>VtWc;%PO=Hx%+UlLut~kZ865~U4oJ`Bc=!VSO0Q40o@WLCWG56WvsIA#r2mJ{3 z?h6xSK#>x{I2;3&1}|DqFVfIk@j5xfsvF0K2TG%5Pf`@uvGgf#tA7R9_eRTxFK<6q z{w8x=AksA;mx8v~Qf$nF7J?D=>^kJpJz8SXskE$)426tYPFrCS7~0(tH~Ca-72Ogo zqW9*6;$b)d&wq!KmUAp*wOdO7T8wn*9#!P31-OoMPicGMZGR%-}x+uj^l0nh)J$l@bVPJhT%h|;v07#xE4Ptg&mZ_p}MHu9+*M@gkt+*nVm%d%L?3Zyw)9@Gt zO7!N}y=Rc?@T~$W{o8m{<9zql2=qp3*o0###C2JiSPjN!2p@3*c)0l{opz^1Xe;I?iAzg51 zP+=LUXElT^_!Bra4t7ir8vT-l$!5__jg^BxAGR;7J}21J{UaoEe(LBfHJ~6z{JGF( z7<}i5UvK(|AH?sUq(5aFA)V!_=DPiTS4WqM?4|Fb>UzdXsj=1#k~LRjBVWK5p{75Q zSC*?jmn_lS8P019m)KxVJMpm%XOv&oXUe4ExJOXs*j*+xAt6P^+D9Atcq=*?DmD zt*U#4n5aTy)YQdRro{}en8&FjS~q!gvp(yJz4`RNbG5^Z3Ke_4X*|mQhW?MHLG0!m zqph~aq>=&$SRQFQ>S_tx0}EOsH>+g>&LFRfS(SnF6sAEDA{V`Tw=pFW8&(rGD5v1} zJ=QrR-}UQ7Octe7UVB%7rE2Xd591>OEX@Cua%V7L<8fM*)D0wW4N@n08cpHq@MR#p zYz#!jZRf-C$d;q|rC7YP)c8%X@&lz0-4QRDwP&uspO2>KBfA0Deph>ew2l&Zu2Sge zZgUxRU62D)bvZNoT0foV8?Qb(G2{A7aF|0NH8C++VBFFmUAd7fk*z_2)Lie7j75Ht zFDsZ#?&OAa_|~ZxN1-q4n^4w>*<;dtOJ$811;{2IbotcGN;5{`S&)EnY~$RhZ@8x8 zkET&n-?h>MNM>kcmm`Fkmnav=ls{Np_;`*61V)$+@{h0P*-5uLBQU z6M*>(5$*4~=d}<0U&avt8WlAaH_yrzp#(U}9qfCdt`NhspCALnkmB8}biG$AUkdJd z4_gX@vX-2>|52IMOtWQARM|{9>Web&jlV@nY^}y5s)h@#vcVOjb^fMfUa~8VuRSbJ zNx}b4gB8Mk&l%gH`Beb6<>eTO$-?j@Q+w)obeJfi(;R#1d$`Y7WXKU<*?lpUG>uHP z2~?q~^yqll9tsPT;v`jDK9%8?=P$gkKT;TLp0_?u{sMaKZKZ?O#4^7ASxoF1SMt%; z;PnSVS3_RbCR5X`Zj$g=d>Mvb*HTDwz8lxjNhi8=U#ZXz!53JQaBm}W;VP829f4ct zth<{~lR~|AbZgu@unc&4dFFWYC-}l*T?NJLK+5bZ4sj}n`&3G;2KQRtGF`xoP-Hv{ zzkCP9YiI>s5Jb;bO4bSG-t;9t>^No_tCVHRw{TK`^v~aZc=XRQ3PfwgDXtrQPlA2& zMyoka0rHjR+4EN+rOyb0qQmfDZeC>2XB1WYX5DPeb@1s2J-#=?6Nte*Q z4F)c|pnlXU_h1_29#{V5NL6wN*BAA9?PN_?QaUY%-mZS+Mv|NK*xW^9a`gQ;Y}AZz zC$2N(OX-Of+$yI%U$L2ToS>pWo0cDM+L+Kj8o*k9O!OtUgur+v;~|0nDwdkPv)4If zDX#|pQSU>&zCf7*9PH=v1r@!0+M@4EY6*fB4P{HIhcy&HB2pK5e3Dh48UB$uhjit# z3o5Q*V=4TW{ew+8>I4&uk?Ft#uG7F!AITdVZ}v|1{D&$oubUav)CuZ!*2OH0Sy6C zmRBl))yCDI#H199=Q-DT;5s<8`rI1DzbHx2g#?$*wqYdN5N(*{+H?_ULucKOD{1^G z7a5E|tmXzi@!=SILyr$3I zNOY2*O5>pAMe8`aFMpK^-8?w6Kk|0%dt#U7hs*~>MM&6sx5k$tYup;nJf=3ul5nj z#F}nxFj>5}zo^#P99Lgac+=q_k6?s=0GYMTYn?#&<_+(}k^+XjN0egRpol0~qO7=> z9wAc5aEu#>fYG#^9|c54*DANQ13y&OBTVm-=n+@-#9Sz`BnWTNT}Asiz+4!*XVEC3pWXksSUP#urvzfxkH9H1*j}>UjN9iIqA3rp6OiK4ehfU z!hAi@kRcW9Xuf{Ww_TJA_7tlOF&TW^O@Kh^cHNEv}X;!drJQ7p9JV zej$ux_g50-qI_bcb^9ZA@LLB<6aoq^qNNW28+>hw9&qa76Z1Kk+hj=16;gKa7%tIq z@W6>2P}W5bAYF&X9oz$-RGu_xc3RU9aCcY-6^}b0HVnvXAYNA)q_&O3j?$AwJQZrt{-$kHJd4_*je@8_6HzyfH zD-p5(SG2E>~zxH!v62*%u~bi$wM{F-&kmyUiSK@~abVTi4%0ldvW zhl$U#6%x~O7R}0wP@rdEGz#Klqz5ob{ryB2}nzzpyQ|95U1%`kNEkeoT-5@8y92^cR_sum;c7i zx{?C$AvR%TLE@F=f=~jq0TQ43R}7q#kC$YT`!eR`wF6D598}}ab(M%8MR~*y*^*l+ z(;i3NG`Qu{o1ST+w<9Sku|~&Q%SAwdZ01V7WX!7UM40J~FdUsSv|aU%BH;)ssPg{5 zSj7jz2;Se#VwU+wD~oxdXgJ#M-Eg!CJ|fa%-9E}AvRFl0zeC_Pb?PC*kI<;oNWf$9 zLqMcQ?0Jo5#+)}0M3hVz^m;2QcCy`K&vK=kq+&7S!1@|o$4f8(GX6Cmj6GR#gM?(& zX;;(+K0wlL6}wecGN=*VOtT_WPkWX&@DXmf<@J!j71koL_av@^KKhO*58UmTC%@gJ z*S*VGwPc&D<02DTVZ_30Q37C@_$*5FkMi;j4&@5Hs?;G8EGlrq2yV%iVeK|SA^NmM zVQoW#%N!C>ey;gi4bq3Jyu2CT;nl&_(p+*L(Oj1)t;W)Ff4Dq{VGy;u1bH|Or_Q(! zA)MHAv^4@nEHg1pH0rAfp3ItN$Im@{KsQ3bT_#ToMhqkC9Y6_J$ybE%<;ZpD@SB7e zr5Hj!=8aFs<$XD{)TOa_4&LXUKSqJ{sK(PKoSS!0xp`JNWg0@1ou6|^$J6Z;Cz%pt zgQ?8Cq@oFadVV9g^n$aOIn3Z?&&qG}?(kK@J%l%1-(iAvs_)BE6){}=9BwuY#v%y* z1gJf?!hZ($bn&*}(C<0a zue`8)Fn|SHeS81Q>$ZNPETLJ_@KEUDIh)x~KQ3Mqp*V5>ov2CcfPa>me+9eG+?G$( z-X(G=9wgdqfucUmaX3n_T8Zrc+eF~)eF@m;?vdMQ*lE4?L|R3ux^0V&WEPKkW0ZYm zRyg1~3_tr&WDQcWCj9HGZwcZrAC6tTOmFF2`*WzVpU{11t00!C!z@d0P8K*0C-BWL zNYeMj_hK!3U-*&!3vz$znnO1M&tCk_6ig;5b*QoZceU?^l_}~o;3lxRT}R&=pF+R1 z>?Ww1G1O?B>Gpf7FHs6U)EJvF^zkL3n%5rU+Pu-_4I0Z(o7pmg3Vhz2fPy6l(^`=m&V4POj@`ID3=P z{>F!w@6`fwP`Hk^7?RVfVY65yc6yVM&U;53yr?ogPm*NGyK(rLc>AV#Ui*W)CO^9% z3MQu+`Immn1nhZA-Noc=5k}Ts%Go1yU+nj~rgJ>{Tf=3TA%0B-gjKux-2M>m@sr?? z*JRQMzKN3HacP-nY$qRK_Yb8AO%km8SHON>MrX81I)<;09E7+f&qn;G1M#Ic@(7*! z!`aRFh;~$x|LtmRXvXnOZ*3&*bA$e@M4m&NbIu}c;sCAtS2bjqBYxHPv0U5vg4(+p zKGP}+KkHe%x)*(Z}D(a+g+-;yg-LD zam{)$tY>x!p<<@+rJH|0a!N;kB&8u(4ZeAk(^7WE(;D~osinbree4&rK8UrX)O-~A zLN4R_5Ed-OE5}h24@x9dpURqzCLw>M%AY!_N!epx5c$Bisjp+a$gg|DEA*;T*eVdt zsDS7R4Ap_%-+kXYe{Vx(fdbm}gi$GMC7jWg%bTX7)cb0B-`MNMi|Aie zQ@q*kKtvPpQpotR<(22y0HL#|0m9?O3s`nHX(z@Eu-BrfoY~fW2mhJHIQW*=8#+P* zZk@XI_uKK5GPXUMDj3TWNbg?l{CNsEvri3n+1abT8U1X}&*?v+QdSqh!#=gRAWFyq zqnwt+7vNw^pQGsPUUI^mo{RiF`cpO6oe1nk%F?(@T;e@og;fU1>Oy12jJXIx{eOu& zJ~@{M*!dcB!164=>xeQpl8=CMrOI4U2^|qEiv<0vV!hpZ!Q0nHNi)E z{)b;|na~lJ$)JV%j>^_c#>_qR$wzy`St9!!3G&`Ydu4kZ2|q7dA8}pT=V%`L zn&ARURfulra?+Id03DJgxe2tgR>x0=?q=Ty7S-fjAx57FzmE-?Mon`7TTnBw*37zk zJhpu}dZ<5l^`fo(XJ+opQxtAY0>JTnpS=9Z&Onlx@9|DU+37*jbMu3SRe{74%ws!^ zvYoP2lD%VrlbwC>3>kmBqJ8>VzvG>?X1|l9pDfO&dt+s@{`)_;0Z-qZhHKy1-gOyV znmx{5Iw5hfaYH;wm%OOjbg3< zrqEE_mbdR&W}h`PO-w)esBE#=?b-QIOLygR9;<4LS9bsucr9H86ky6lpQa>9U8nf& z=!Q}>DE3GV*cR`{%g4Vb+`5xR%U(i7HW(6 z7hO8XqQ*{tm^B}8Y@9avZx5CdF>ZYsJ3S7=?RE(?9nEJ)cKV$9Y|FnaQ?9OD<>I1} zlMrSTJo~UYu~DAi2ug~IiT=kXE-d+9o0xEdx`>pusQCXG5)u{{y7M2KgpkmGZDNA| fu}KR4-=&0v334%ET}0`!h@hZ2B?zRfsY3Zb`~c&; delta 8041 zcmZvA2RzjA|9=RDD6+E0WuEPjoe;8jMv1dWC^A1r2-)X^jI7LX>X5ALJ3Hab?7jE? zAN3#K-}iry$F2AKUiW^#U+?Gp`FcLzDBBv_SO`c&NK{ZvNKjaai-=#{#p;p8D=0JP zqX$9)LLvhHr!UTg7f+B+I4+_TS3o#H5{DnK`!iEpWM7{r6n~9?oH;IlgoNxZ{4RPk zksd9#eZ*2d_`>zsw~_JCKY;|_!Bmm@Mds{cE3SiWE=(&kF)x{F%13o+Hg|ceC3vZBX<$DX_OfzE{Yj1Y5 zWs|l7Z%p!BbU8MVlkd4WHCR?AS$4PXg!6inv);M)aCTb%4 zQFhNl9i5-_LDQ1+wS*72Ob=F^ZhCvJI!ylX&py~URjQzXRkt(uZeW#AI zekpbFG2G);r^4XtwDqQh83Jh>mI@X^7aYJ{S9e6{tKHN+^2y&hOJ9d=zMW7EUL-5L z>abAzT%$xUbDN*6Z}4F%OGQ3ut{DrxDxEMmyVHGv z4t_P$hFB=`axWI_QQt5k?Um{NX=i5d)I*=+Kza}|27ULY6&?Jws0mrae?$hx(({^avM`Ks(8hcqclqA*_dgia2yW{@w(jthEll zfvsyjsa{zGZ&0MY_8P*jpprgvCFNx#mB;5_GGqhBVW~oDRk#F?J#SQ0cx3dW` zQ!5_XlBVI05ay>c56=of&*~vHr6PrWG)$s1W!1+?tHulv)-t6@Bhbo33mYDvXyjRqeCT39V`rO($xx zj33px8{s$h(MF#Wi`Xvt5JXjz7$s4sNm|`#OZT5$rzyL^r6T&2IuPN!LpO=~q;TDX z3ClYQ52cQ9MkN^n2^BbgK`f&E1!63$_OWk4SVt2uYcYHK`kL`tmVUpk*E-snZzG|-F zu?c>wQE?6g113Fwt#HTf8iz6%eh4Bil}ov~gVRAlV!xU(_@VQFdR^?p{m_%O)E?e? zmkY3oQ>0En+lpK z+S;VJFZG0?++<%z*Mm$%bcKM-awLAME>2L~dNPEHVlr&uiLtge#7#Ii3I7{S(JqR=A3x z`mHYGs*RjhGVwj>Q9mFzGPhTyVH4fj-aPw<4%5VaXu$q zbu#%0KbNF09W$Dg+pWlXGlUo+Hulj{^5+0HbL4MW>werQqA$b1tp{Jf?}$OkcJVP?N5HjzYa!E4IOi-|Y=7f&!;V$LHCZGv} zfheI4qROf}FXd25y{u?bqH1HYqh>`+#OSDw)5hyp6-AaY_4E7R*Fx^wf*HnB5SEjf zM^Z+2x-}M-I$bMbKFxhN)hfj<7X^2I2Mf&RfkojxZxJ+6S`N)m^_|XgY?0m@+}chc;@^yO_tk$&f{68vOyTWzGH#Tb22V9lMX~A82T3t>nByrB-(8) zs++a<_tl9#jr|A8j&TmlgVIz@_b~%qSMWAfnsZf9rWODCyFu!SR?NnWq znq;Xd|C4Hi;wl24!LsR19OAc^OU2DVj`Oj+2>+|^$VC{8PAs=x;xeTIRHv#kg{W{o|x^kc#`O<^Tf&qj58 zvl;h4IZO@!PM^RG4!-g&;lhgqd56~D{G#8BdfYkgA2aJ4g5+1moiq3nx7$#L8P_wY(`&RYM2<5D#US0ERxij#gT^{^Rr!w|4TQ`ykOGg%z{ ze~%d)9A%kRxuKKyXSDyoe(-2K+{pBdFUm5g!$mM`*G#q0H%TaDG8Z=I$%Rgz%q>iv z%x!Mj=#QXta#a@5Aa#~GGq*>53W$w4eOIzRW=hQc;<8v2 zW+tW%;-?!hv1cg-Os)>d*$BLj=l~9LRkL zJ{9iScS*Px@er!|${Q7`n*t>a87uyK#7e|_0?84Vla@+%1OF^s8^BPdGpPAj9mG@q zCo%jXrnYja!oZO;Z!Bg3pnz+tw~dQF-!P%Cbpz6e;mGG_*k0QaA)S{q(35Y29n!?E z-Z*#$4UQNqHaRsTWSM#TutiL=_y|Otq|PaxRa}0Tvh1(uczI?AK<^2kfc84UQPaa}3;1nt)>EddlzEm-_BPzg`w2qIVTd(e#F^@e#m^kn70=?}PgI|Wj zXjrI4h4noKw`P2(;O%)t=^x2+R$7!;|XAPHnrLr~M-B z8J-FMX5!1n8m4mw`~WF{6f#v~ymO*5+p2P_Bk=K|>J1tgo*ewaRE36XV5l2@(R(O8 za=V6~Lo5##r1`vFb#eyvHEiaUFPJmLc&s^`sp%skC(Cix8l`oE!%Fg@IF@ekeY4iT zZ1#TBz(my0sL9^;Z7t|ZA*T*KTY1o}3$Bt0!3)alH$ryAWRC+}YRgY-)UjpSiM_@p zt{oxIBDPF*RZP%Eka~K(%hF#y@ieI!Rq;L-Jn_qKJ4U>a*;uCnOJZ0?HbH7Lb9HQr#El({Lvq)d_kOdAnEH5vLBzXqscUUEet7{us| zJ%#{6$HBjf;=DZotRXi9tI|u7|Mg*&R*&KMQFkWK5lOoae4tQs_DMK#dc5>Zdmf5F z9}jso;<_QusH0lm^>mVyX_BI#LJXYbdwfp-FH{9;A_UreOl73M?-==SIg@zEeX?M? zWO52Qd0pB1uTU*Yy<5XTL`*=jmCBx%UKp3FtO(S##std0lTW&G|0)NQxJa9^ZSdC9 zUm@=?m?W6p>AS{1os)bkodHEvcOnz%-*J9>WZx6INjz)tc3;c>Hs}ctU$54e=d~@i z#J=~h5oG0DZ>I>Skaaqe1duEt?;f@$U_z_ihu82zQH;NCOAF5lNETB8Hjt!7$6z3qDVfXFPXsz-n~5?P9r!0{P8i{7moh24m$MA+qdfvB z+1Tp)zBL1>agWAx`Su0;hkRQ1F`^3dbTAsEA94?c4vA-{efp)i!>mv=%eKmWsX=a` z`3`$CSIM2BjBKQZq0qnF!9)e`(vo+^z_flMdIXw2*`GL_a5deXwYd>2*FCzgoA>JN zAI{OcJd2fUz~_Jbea77^I4Y!YGg{e#;JH4{3#6h2(YEm0JBS=>6Bw>kJibAq^k~-i zM(RtsI{A|oJry+tg8}+mxtYqK(Yr|z2z2L*afuzJWMj$$ zFt3RiBr&K<1>bsWF1KNoM3-BDSmt4AJf&z8Oto=7$Ew`50I|@Tw~R#!vrm&Baqstp zuMDDcD;6H+yG&9=HFJJKS|GX9(;E${26`|v+^<_HjqwP}wgenbI=!2{5m%ZEP%?3I zzh5>6Nld9>sxJTK$GXs}sa2^C+OD>HgwZwOaje#n(s9pqB$MK*kK29w)OVWm`G4@z9ZIcC?gsQ5jG}zS0a{f{UooSyfD|9d+IzF;sF=W zDT>N~GCIv)MPvcR`EFmaKq&x0XU_fL40_{wpk;d;(zHmXV6QrpvV@U*b4iuL>W)R!7ya)$9z+%UE>zDq0*@26_zYgRT@ zxZ1-0Ig-%-a1eDZ(uMen0B$EU8`;;`?-J5rPB!XkWj%!~<~)2TF4*D*N=d1@B`E8u zE6zJN{PmULRgbY36`z?%Gs=zc={8$acX78^CrM|Q#-fLLsSn` z4vpZ-Kg!fl8XgoArcpz|s@A*5>4AIY&%Bms^)M=`=lDgc#yVjyRFV;&^;-mE$_lY4 zz^DMKS9YpyjXe=*g~CsusgYKlwC-a_I%5YSl@;q|4u4*jJ+jCy zZuU zK`b#8O}tB6q>TIEq}}|`_mMC=HU%rsiT*FcGF7o--3RALV=xYAMDWYKwc#P`ld`_u zwykc^wynzktEK*F5LbD1Mql*sVUMjcBS|}U;ATbO63iAOQZ%uborF&d4PpH+t#l@Qx z__mRdqKt+G%~-&&J#e#Vc)9cT6}{RnEYJ&C^6RFu{sw$sw|9gJ!erY#p9 zBQmu-#*IfvJuBI3LVQXQ^1stSTj4Nl#{6vTkcwOl+2ah)yQl*Z;1LFLrgfN=d>3|u z1BVJ&YPRC!VV7Li!l8o#48|(catp%~ah)Q@8SI|cvOz=zJ;LGv)(K5xow)JPa& zy1FUvtg+VAOoJ!(pU3lmT`P{=&p?jSL+1~^KX8GMP8~cyXLm%%{}s&~US2pVol5X{Uc7`ENyEf6zkzZ7d%$`WR`Y0{y|SOSQTz=TTi&6pWC!%dp18|8>Z+2SV8F z-D)G+oJFO$RJdd#6*ZCXViruQ&w@Fpn!PyyOYLc-_meqr3+evt&TqJEC4KPFUq*of3nXRfga&Gw zl|)mZ?hV?Z{Ku0fhl*7C z@O|ARdUWqag^syM>ESV(kzd^`rQSdr&(GJ%A6PoK^Vs-g zT?)wy(7_z-PV))>GL+$r(6)&)vcAi8t5R>FDxh$N)qmfIXGX|cbK|Ms%m>7+dXxAK z(-jkqfA9nXo%&|@e53{qb*DY?vB@Pawce~w;@L1=!oEJ69nM?*#Ie+U;?D1f9+L4M zp7QTTgt6BAZ!t8(3e?CX$Cy}owuBQY8Tzp}()!`+v(?NByQvv#JFMU1hBi{XkKKDS zlMD?mj5FH`kG!Nrgd`ROQrdz<$G#hqqgji^hwsKz;g?>6N_x+O;pbKLtvtNELOP+s z%n`SNJ{|E(c1a(9X3gz{xjj zHg4$9w}Y+x$SXnFx?{ohX*;%tsNEO}e&J9T^hdLs&kL5McXR5mV5tjvDVdR7!sYr@ zv$XpC;V0VKdECOcFB4*huxG$N4%d(i!>X(6gxwV4vSraE z3{~wstBkpO9B+Xm@c-zQf7-@hqM>u&ALzoD)y~E{CghLYsXoe7y&-lae#&pGqs$G- zUG>%${kearP7l9zIkcu}jk#M|RV^6kiD_Q`Hb!=uJK%XCV}C|DNUf8-cG#HQQ4`ErEeA998th30=cMIsUtlh4(njy2*WzVDQA`}YSi%e?cK+F0aA5O=( zoF7Hz&z>R9*&uI{EO}YRkTG!5RL{`e{jJw08#cot9!oQp%SxS%Kbl*QW?7ldR1BUKbIUrr?QtJ#*jxQ6g*;LK22TGcmU!V`3`JH`k%edIyGn}kDTdAA`4(_2R z^*swV4e^y7yc2-&AKR0D+Y9^`+eH5^080I`nW9hx@?^Xq>>hLsN5|6g$Ct$zIfa-yl~Fr!T^N57hilasr<=YLt* zNbVmK)x@-%r3uDFY!moix=HM9#EDsSuaD9e_%;RdJ zqoyR_XxWcwcA@n2eGGN!N!UjH(TXh6*}dTlPQn&VSLMY)LGorB$1`BB5i*^2QF{b< zhw&~{9I>;L%0R@J1uh@V&dfwt9|R;sMFsw=O-xYazuQE_5=?|2Sc?e#k1jzWh=jy{ qwTVeU{;N$~MD)MgB>v}6f