From dcffb2fc638e51ff8763b1a52ea032805c21d53d Mon Sep 17 00:00:00 2001
From: wendtalexander <wendtalexander@protonmail.com>
Date: Wed, 22 Feb 2023 13:40:57 +0100
Subject: [PATCH] now its the right amount of individuals

---
 code/plot_chirp_size.py             |  33 +++++++++++++++++-----------
 code/plot_event_timeline.py         |   2 ++
 poster/figs/chirps_winner_loser.pdf | Bin 0 -> 18064 bytes
 3 files changed, 22 insertions(+), 13 deletions(-)
 create mode 100644 poster/figs/chirps_winner_loser.pdf

diff --git a/code/plot_chirp_size.py b/code/plot_chirp_size.py
index 92522a8..20ec7f7 100644
--- a/code/plot_chirp_size.py
+++ b/code/plot_chirp_size.py
@@ -67,6 +67,7 @@ def get_chirp_size(folder_name, Behavior, order_meta_df, id_meta_df):
 
     size_losers = [size_fish2_row[col].values[0] for col in ['l1', 'l2', 'l3']]
     size_fish2 = np.nanmean(size_losers)
+
     if winner == fish1:
         if size_fish1 > size_fish2:
             size_diff_bigger = size_fish1 - size_fish2
@@ -76,11 +77,8 @@ def get_chirp_size(folder_name, Behavior, order_meta_df, id_meta_df):
             size_diff_bigger = size_fish1 - size_fish2
             size_diff_smaller = size_fish2 - size_fish1
         else:
-            size_diff_bigger = np.nan
-            size_diff_smaller = np.nan
-            winner_fish_id = np.nan
-            loser_fish_id = np.nan
-            return size_diff_bigger, size_diff_smaller, winner_fish_id, loser_fish_id
+            size_diff_bigger = 0
+            size_diff_smaller = 0
 
         winner_fish_id = folder_row['rec_id1'].values[0]
         loser_fish_id = folder_row['rec_id2'].values[0]
@@ -94,11 +92,8 @@ def get_chirp_size(folder_name, Behavior, order_meta_df, id_meta_df):
             size_diff_bigger = size_fish2 - size_fish1
             size_diff_smaller = size_fish1 - size_fish2
         else:
-            size_diff_bigger = np.nan
-            size_diff_smaller = np.nan
-            winner_fish_id = np.nan
-            loser_fish_id = np.nan
-            return size_diff_bigger, size_diff_smaller, winner_fish_id, loser_fish_id
+            size_diff_bigger = 0
+            size_diff_smaller = 0
 
         winner_fish_id = folder_row['rec_id2'].values[0]
         loser_fish_id = folder_row['rec_id1'].values[0]
@@ -107,6 +102,7 @@ def get_chirp_size(folder_name, Behavior, order_meta_df, id_meta_df):
         size_diff_smaller = np.nan
         winner_fish_id = np.nan
         loser_fish_id = np.nan
+
         return size_diff_bigger, size_diff_smaller, winner_fish_id, loser_fish_id
 
     chirp_winner = len(
@@ -126,11 +122,12 @@ def get_chirp_freq(folder_name, Behavior, order_meta_df):
 
     fish1_freq = folder_row['rec_id1'].values[0].astype(int)
     fish2_freq = folder_row['rec_id2'].values[0].astype(int)
-    winner = folder_row['winner'].values[0].astype(int)
+
     chirp_freq_fish1 = np.nanmedian(
         Behavior.freq[Behavior.ident == fish1_freq])
     chirp_freq_fish2 = np.nanmedian(
         Behavior.freq[Behavior.ident == fish2_freq])
+    winner = folder_row['winner'].values[0].astype(int)
 
     if winner == fish1:
         # if chirp_freq_fish1 > chirp_freq_fish2:
@@ -174,6 +171,7 @@ def get_chirp_freq(folder_name, Behavior, order_meta_df):
         loser_fish_freq = np.nan
         winner_fish_id = np.nan
         loser_fish_id = np.nan
+        return winner_fish_freq, winner_fish_id, loser_fish_freq, loser_fish_id
 
     chirp_winner = len(
         Behavior.chirps[Behavior.chirps_ids == winner_fish_id])
@@ -197,6 +195,7 @@ def main(datapath: str):
     id_meta_df = read_csv(path_id_meta)
 
     chirps_winner = []
+    chirps_loser = []
 
     size_diffs_winner = []
     size_diffs_loser = []
@@ -208,7 +207,6 @@ def main(datapath: str):
     freq_chirps_winner = []
     freq_chirps_loser = []
 
-    chirps_loser = []
     freq_diffs = []
     freq_chirps_diffs = []
 
@@ -259,6 +257,16 @@ def main(datapath: str):
     scatterloser = 1.85
     chirps_winner = np.asarray(chirps_winner)[~np.isnan(chirps_winner)]
     chirps_loser = np.asarray(chirps_loser)[~np.isnan(chirps_loser)]
+    embed()
+    exit()
+    freq_diffs_higher = np.asarray(
+        freq_diffs_higher)[~np.isnan(freq_diffs_higher)]
+    freq_diffs_lower = np.asarray(freq_diffs_lower)[
+        ~np.isnan(freq_diffs_lower)]
+    freq_chirps_winner = np.asarray(
+        freq_chirps_winner)[~np.isnan(freq_chirps_winner)]
+    freq_chirps_loser = np.asarray(
+        freq_chirps_loser)[~np.isnan(freq_chirps_loser)]
 
     stat = wilcoxon(chirps_winner, chirps_loser)
     print(stat)
@@ -299,7 +307,6 @@ def main(datapath: str):
 
     ax2.set_xlabel('Size difference [cm]')
     # ax2.set_xticks(np.arange(-10, 10.1, 2))
-
     ax3.scatter(freq_diffs_higher, freq_chirps_winner, color=winner_color)
     ax3.scatter(freq_diffs_lower, freq_chirps_loser, color=loser_color)
 
diff --git a/code/plot_event_timeline.py b/code/plot_event_timeline.py
index 32c0c18..cb75cd9 100644
--- a/code/plot_event_timeline.py
+++ b/code/plot_event_timeline.py
@@ -48,6 +48,8 @@ def main(datapath: str):
         # Associate chirps to inidividual fish
         fish1 = (bh.chirps[bh.chirps_ids == fish1_id] / 60) / 60
         fish2 = (bh.chirps[bh.chirps_ids == fish2_id] / 60) / 60
+        embed()
+        exit()
         fish1_color = ps.gblue2
         fish2_color = ps.gblue1
         fig, ax = plt.subplots(5, 1, figsize=(
diff --git a/poster/figs/chirps_winner_loser.pdf b/poster/figs/chirps_winner_loser.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..728ff7ea6bc12664c3dc0698e6cbfaeab5fcc83e
GIT binary patch
literal 18064
zcmd^nc{o+y7q^f(8A38$^X$%-Ysx&2*F2BckojIDrBFhKj2S|anG6vblQBbv6d|Q3
zDn%nodC#GK-v+<;+w(r}e?Ol6oOAZs!`f@Fz1L^0v)2_gQqz!xOQC24D<48CnrYw=
z7{t%{2#vfv1ZsNB9}j`5It4m;`?*7)Mo#W{A_M^xz(EuhY4E<TU_qImDrorm211bA
z1)!FO&PVVrfe`7PNA-|EO_M;UKs*GslVRi(7>FnMLeStN4b;TL$<@=>9fJA(mf+`N
zf)9k)0Ntu-0HokU0wGXsAAkgv?^o6DS8a&hPvpS=8~{i+NId}$P`5)q6o)7J1rc0;
z@qqN>AW(h0tEZESUkK164ERT5QBu+f2plT|v7>>i0^<TI6Cv20!pgqBe!vqPEcr7p
zz}G*Rr-%1-5A=Y*zjvtN=?x4Q0@d&ac%X)N@pHw4R1*UUcqbp4kevHAx_+F}Ozmsc
zOD+Tk{cG#;%e;xq&`3K8PL7MJNvfA`C(`NaY(5~_`1r+aL{B=-4w=1P`-a6??Wt*M
zkF!Gi`TOtZ#_pf0f>yVm+VtQ2{9@H{<JG%6bG;esYtHrax|?5I0_R;ekHalDgOJ^a
z^ggYBy#K_)g7#a@m#r;Xx{+u~hbIxz459LLqS<rhJY|Rc$lO*RWzEF}y@h=$=ORuW
z8bNY#=TH`xo#JVsw$tmK#@W`_HjgDEF(q$5-e2zukyBhBs1rK5Nq;;f_S7l5i3jyt
z{VrRY+KSO0wdpH7=X33B;r$e^JmXz$uVWY+=5cUE1~1cG?ME4NmzP%qLlC|yR&134
z!+5sB2JA&^k+b@u9H(Qks&-`+#0j+(kscpUFMS%A_abFhe9(gKboRAcQ<|Wn%Kg0F
zdKP}RrPvIqj&rdSzQZ-c_i{?z%Dm5L)UK`;r`_o_(W_dEbUe4O=$?=0D}8FUiE0TO
z5yr9NH2p_ON2vs|nl4N9By5mmEmt<?mKPn0Q(0~-@_xFR=GZ$-d9I&@4{tKS*LUq}
zVKLR!H19{Y9vcUmF6-b^dEF}>vz`$hzOi6|6+dKC<}Go9<eVy9mie1QK96#k52&(4
zHH=i#Lxz0z^>Od9S-;cMx!1>JfFjMRg~Ep3s>0X6ii@4^<%wasPE3(R4s)7N;3xw=
zGG&mXt9Dpcb+7Ov&z8BFwfh%@BKJ|aQ%}hj%ld~TU0LR*znjf)Eun~HIicEt^OHc{
z!Sjl6N5e&J=QmH*ZY7++l1MJ*GIQH_5auEY$B>HaC0snScb*1I^v_*6bI($=f{9Hq
zTZ`8!o8E$jSFaOuHn^E&n6~m-!SqytPDjb*Moq~8%@Y5t(g%myq!=SURXqN>WI`@r
zbJ(XVMa)?Lw#;e@``eejuzoBL^+g4YdLb28*%HM~e;H;ALrK!fLY48XCAp71Y5mic
zXMGJE8=Aj7?BjQ*STk5S%6eb@T0d{i%Bp&tkg}G|!qIO*LrC1{zCCLcm6aj{`=hNS
zDakIg`j55bH<COO+9?(Z=VB&m;vn&JsanyOq^uW(d#@NSU-?MaFJ--dgGwq^E6=Z5
ze(2gvjg8#?5UHvDFjK~NmQu8jmxS}_L`;_tBpR^Pmgf4h3EQ)*a(ci0ASRPJH}Jv4
zj4`^K>F#jb)8wxdr}sYnkZPKIF8uu`Ber|m*U)Ma^9kDVRT9C=<D#jRMBAGj0?+(k
zt|ny%NCmMy&}_KKU?jR|H)f%pR4sb0^mRtl(L^Dk<1SV29_w^@mor~?6x=f!Nj}MO
zgzb8?66aA~H~l9N4xab7=PIAYjKqA}sA}el6EZw{zpE=^?3PrbRLGPqBkPGW6>Q5`
z>XcPT$hvpPs^rt5^CNkNlZzuK*-vubTz@h+<}xiPcJoGiTU+SYXY1qZ%PmhXx*G-G
z_4Y3fV_^|3(zxQbd_zOZj6^@0Kb-f?<Pwj=yC$<sK2u(SYsD5z5)x+~2AL9;jK`;l
z4RX5dZ%*I4ko=JEoD>!0*@HKdx7)(Pj=rC{Kl7?_@@V7yY)fOtQE3ZjmxnJHPWo%T
zPBtrVu}N8&JL;@gC#x9MlRteh$;Nwt=qc+K(xUjG>s=FMO^@B*>LmO4h3La1CQM?!
z6qMHw>~|t~fAlb{T=B-pkQYzJk>0sn9MM<8xVR`WdLxR$-}KbLTLy;`QxYvDYPZVq
z(i*xv(>>|8d5wzsk`2rJr)ky&d#I%`ViKo?&B9x+-g@)yGd0WR0L>-o2O-RCCb@V^
zCsj`7Mu*#%;_Do$9kwv`kJ?(oz7)?6ZAH%H40<`nJUL;HG*t^2;>MhdPY|0YucD^T
zi{f0!eq!+&Hq2!jk>6k8Mjuw6`_v-1c<uV6rc;$}Q>vW(*`sVZJl)oZuNBAHMHLUW
zy>+)}Aa!pMwUb!ux}9|SdXOBeUD^TZH+6YD{duRN$a8hWKRJa7?IAs*@?ic7(>iIr
zY~vb%0pf5354{xaaJ?&Pao?i=a_Lk~c3<I*&E&QYjoL-G_0!O>!%^4be4?)LJP8Xw
zUn<#=_@4WIKrOGAO$W{T@F|<(+;ksrk@~Nw2Fa#namlXNc9$t?1PAGk-<=f2dWxo!
zAFlB=3txOPa`g!F`V?7M2}7Hn*x+iVSKVf(@<MCHWQlpTx>MChaWm<c;${pAuYxKy
ze7BazHkJ#&G=&5;JxDS7c;ibuo3@0EYxv9A{I0$JsH#HdxsRXx9yMML+9I|2awuo^
zCWUkd*M=Q6MQj~ogmB5lDI6<5#-;97w)-<LpXGgn(R;3Y1v`#ZyOCr(-G6e%fH?PF
z`SsU?@O-kNmHSi>k60CN=8nue*8ClfYZ#a*Y3)@w@AI_0aE5lJlR8ZT?Re7}GL>e?
z=h{Fg{WKdCmYDaZH<>mzasth8mpsLe%BqyyK%#;^lHvJaZ~W#<RS|2XGL7m!!`S0{
ztzL+0<%T&C_xgohGxkDgy(OR7TWP<63k!Q9bY9({%~$bunbUcrw>%E~S2-AkU9T6P
zIUt;~SAT6FYo*Sdi>3_oDP&csBuHG=3?`kTby+D!yeNpb+t}c;z|iyjr0L$JL`{(^
zA5hH+k%JuqNt!`6Ytn^e4cB}sRJHnFU!CE$^LE08pj-AfRDV#~BwYOLwDD5D{^Tif
z*QLyk<JVQF>{Hf@`9s!kOwjNiEBEQVJI$xks&A=tTFqcqurjDlEhaAh`kRZl%kaIp
zY!RI8bBUt-yRnT{x}H^|vG)p4)@N=ES-;Ao&DDw3A=9`$Xj}Nf*+6@3ulntQ`s==9
zqp_y>AH0@AY`xkzu*WlQOeAp1=Bch+Y#p>F@)gWC`PhXnshOJ}I*jO4o5p1;1-dcc
zp6p(^z#5Tsv36dP&0#(4vz8XBay~gP&c%qa_(7iE6SlbxxzpB*CHj&#-%_*(v7N87
z|F*eu)bLb$7%_kt_T|mKV3cAOv@-VV#Pgg(Ju|iQ=a`4HX$F|C(LbGk<d#gZF`zD@
zVNWzXW#7s!R~u3(bY&o7hMBO~7S|X2cH(npiBM7x%;#x*`{?+q!_~bGRhJGgYOvi%
zXnZ?bz5ev&r*}s--@bcZ^3j}w{4plx1nKI_=#qWX`|BlhpGQ<xj1$f@P>#!wH9Ade
zPUmTV(VE0~@^~KSe9don(WRg)=AgC_TDEYE)4vC;+iYybgq%idUpSE4|1#N1TH73c
z&a@BK!zz$?e#o5cR!e3S8grgXUuiLuRZdE^t}R=5(l+QD?8?csEw8FW2tIzNI9{3d
zFLuQH5Gs;WFZR7nSXR;IuTc5SM@GeBDn)6!Xx1#YkSr!iBPn!VtaT|$m2|GNHof$%
z8E-;e0J{-arW9*bb#SQgSt7~xv=j}?C=^rFGn=|I<u(U;r}bagYO<-;D0v)mVR>ph
zSaD|$16vw?AIBMgbehV*_2xNL=d>4>(IwJjLX_g2V(C&RzU;-3<nU<WE5D9~nDRl|
z%16z3X&xjUd3I9aa$Ki*5t>BZNK{OOlqpul{B(DD_*_h9O53p6Gj3JoG)g7^N2w}i
zs?mqV=fl)OahXS!ZV8jPE{OB;?6pWIPQ;=~h@Tl)s>FoVm3t#M)mhJ6RWa*|+)$LT
z^Ep>eDsNEzsn9X&gqSH>$2_7&c+uzW-uFjr&2FD&IeVK*fu-cyeGx9YX0{yBPAkC*
z9Gm|kN1TR~5{$-<<!<<ZR1^7K%J!mKGu~X`d-x57;1$z`=`vgicc&W76%7j)w%|RA
zq(-sb#dq&ND6(GgUrt1x?X@V9V$Y{|raaj6q5$CNc}GGiHrC;qP<qY43-Oz*2Xdod
zN8)gWs|Q5-`LDk{uIDu{jb%Hf43(==7CA+Oe;r24!m4>-<Kw}EYcv6oEKQf+7R75=
z%9$l$o0M}?CB<CQ6JIr@FEM9ae~6nfTb*;gao#Aut9e~ZSXrOM<=H4f44+P<RFh(=
zBhwZ;gimnCnG1N5pF81K%Oz*fPWt%nj8V?Tmwjyu;+}pw`DYGw@IEH_=(g1TR7oR-
zp<7FhqV?3<N=fdi4dX}7s`9`>8Pp6m&wOM|ex&)-249AsbPAJ&e@>`6PBA+4<ojyP
zkv^WSZzU>Fj?vk>gB?cIJb4ZG<wx%hK95MWZ>VW6%DBU<(B{qK`2zn%OA|VyX&pMH
zNMLQdw;tE4*YU=@{<+2#G=qF9vvz-$%}b(hiuKtT^27%kY5G%})h2`Q`=?7Fo~IXE
zM6;znUY0J`%#j^qPdbYzpK?lyJ}6z1%Q+B!mE=(=H3@QV;^3WL(yIlma~m(>hFRHj
zqm6}EIra~W?<qc`b3*e`;)G|}?E5>A%g3(;dAna_Jr+%1IJA0^L}dD8wv#bi7u}so
zj+rpCn~x{2ZGB}pe8*#NrZ}1qFS9;aL+I>S31pe6P7IBvKEaTk(ky!SlVN$;1uKpU
zsY5Z$miK(ldreKAsH%8>^I@AFx~;LTI{Y<W=XTDiqx&H&lO`*(c^~<{H7#wKT>ciZ
z;&z(t(SFvUXZu($QXVUH9?8Eh)L75=E`#j?We{pOQtrm{v=INB11ECsJb3c3XY`Ts
z^oE<mkqf>LtHuZcW5LqL9qd2JAD6#{yzKEr@ry){L&YW8hS%OFj`bg9AUFwluUlVL
z(Zpd?Un!(IXXD#@&%VFg`LOQ7>wbxrEW;YP#c_N3;Etx}oIE$TD!2^KZjtgI7P$-B
z#@{Ww?^g1*83>w&NHh}t!*t%ZAHhr;z{vbw5CNLczuAxoBpNUr|HXhr8*BTLz^Pgm
zS-vS1)d;U-mA}mYx^?}mF_et#!2^MnTN(n2-AQjmEKUr(L=CRWSuOO8=sa0yi=inL
z>!k6d93hw^>GC+TNHdUct=tr|Smjua^t%Is7f`2_V?1p&2VcO(<(H_N6b}Ztlvq_c
zYrUGqwY5*ayCgQCIWa5IXfSn8XSTq7S=%<)-q6R>;b71K>AlJHJaVm+8+^|tq-2wQ
zKV1}!PA}w@kqqaF%{LvPuVxfF=_U{|kY})R1=adVdHiGBab-<kAy=_RBMq^<N>p!r
z{^tY_p|X{&h`f(&?cYdL^-Y`pWG(`9Dg4GWI0p6aJX1)z#i+_lrM6;TD;OTGq|on5
zSiiaDrIc#x%H_CLCKr)p9F6R9f7*X~ycXr}#}e;FRn7l_N}n{fqx)8N2p#R75~pMr
zsmEp)%u3R$vJ>(q+gz0IduE{b`AkQ=yL$d`C;HKq)kDkNDwMS%fxdWa+AK?XuN=9d
zD?^JX4Q7ri9K2e5uTq*YHK`ywc+be*`pcPNkNU5|lh78pJx|mDEY&W91sx{84ORyJ
z?_<@!)3UDvLPtKDCn`@mx20&D#y52{VIO0D?WI8BT*qh{-vCRxdO6Ok>-S7D!<M+2
zQjCbRYL-aB=`)26p9f)GnU=R3s--Qe2<@y7IpcOA40PT6Mwm4A-w8uXqj4~!ntBo@
zw_+ih=NcDiM9<nPDaND^;ypa+-JXS83H4GaT_`IPWuoPy9jA4=wcuD^8;iY1PEIR8
z*?^&Rv!S=3usB3Ak;X~43t6B8>Nm15f59&#8f~#dR&lY$u4J+HDp6>$=QxoR-pVFA
zh7xWdZrmjqWMLTEg@~WMWZykQ2>9QuY9rk$GC0-D;(>OhA|)oO#jNtsi;qryftAKF
zPR#MVeU&M5HbY{aY2CaPU-<xLa-t6>Ss3j2N-5}xZ4YvRCec1df4P3nwxLP8unMBe
zr^7!(P|Yx^7nKz}<(F;KEZQinpY1z7&~a}tsg6J9pbmwu(fD<b&Vtuv>$&gZZ_388
zo$R9*4}AD!+=TI?LYKQ@jBZ=+r_!k`wPxt;r15NBYJJqy*vjn#!}8LFNy89joNc5R
zKI!X=j4NbJFnFVcQm{oqj-Y4X-J8geyNf(dwx$WY@Mjm#84~puvxLUU-My$x_3Ea{
z5^0{*qrqc2{e|nqEexp~i$skI>0MIFtdLmVa^tFXK|`3m@aoYv`GiY-ElJj^riDB=
zFh28NE$f76E@((ZP4%_&2rh&i+2kNqVDs#Z*VCXEAPcxTu{aW%DjNMb9OYxcBdEQ<
zU~h{}S0a=ap%sejO3B@<soHDKmQPJ|(HGz^F)x0;@c?Pt%j_ahn!Q)x1@l^2B%^1q
z%VMSAegvoC$EQ8jjc=`=!#`dtSL(&w;`A0qJ~^OH{**0EKsta;pQ~qV$dcGHYJ^*8
zZ;4J#wO`9fE$tBJ&FQi3gsU1_J$EfprZ&BrLxQZFY=;?)Yfv9O(%dHg>Ert$f%qBD
zU~`?~TUh!x_*879Yqgw^cjr9A`$r0q$(hnO<0K=+Gq97cpKE#4IbNK8`Z=dUw2nn^
zQ#D<__UWOmOOL;7yx;sr8832@|IgV$f_~25mJk;4?<-p&;nuzm7!}!;-YqHy((W&&
zW;`#xJr-hoNy$Sob<okN4n22ksw3m>v~*jcdg?Tl`E~VKqr8>_MV5A&MS7ivS#qav
zft;`X>1Q5KJd@F%UvgL*{kmtr4ma&CH2&=O{l3Dak^fF3vH@iTEOVI2C>n(^g<V6L
z!=#Z$TJQx@nh47#TUae8Y;cIF8-bEmb$*e}u7Xdtgv=9eI8rQj3=KXRnA0I6;f`fZ
z--RI1@%tOoFn_Zh?hu5s*dYi5GyYCcB%F9lqQ;w^O^@pgp2GPcdp4wYD7-g9@CjWK
zJJx@EXaD(J#;U#xF}tFQznDT0RdmDXTE4Jsc92ZGX%+6d9Q|o)5;t&wN=dNchM=<x
zKPeQxj*WS^o`Sw2bkdz*u_<8~rP+CPPWwPa#uNMC$U%-?{9$?|PiuzxaA9vM-EDP@
zf;=PQfgpn{tc4=zf!lO!JMWV^e%u|63*(baAMT2g512GOei6y!Q}LPjrIOa?%QUk=
z@bQ=0&zNHtFI0Bm=UGyFEGZO>n{?q}NiV(gh4WiQE}iyPJQAapA>uw4d+V0SDE)HN
zYxk#SydHM5B|a{t8ajRBH}~o5r%-LEQM)&w-mrc<H3GnwFZ&I6O&ND#*e;L)f&7~d
z$q<`D&Pz9VEJ_h(c5p`NW=Tl{=E26P=9Dahq*dPjab~WpkSu=w8n?a~YgdTdoQYn~
z%}@3FnakJaxcSoJXr0a``%hHP`{pIGZ3U-~SW>QF*3VhlGd$?jo3LPYmzvq@;+fNp
z?e~qhW3EjIx>{mlsR_%cI}_N@l?EZ27JleYL_A-TNfyqFTlq925sW=a|B4oqz{|3@
zE*2VzIbE=PY=!g?L3zt=cAZ<jWl0diLCrd0?EEI@*%^nn;6U!{B1}Y;Y-K5C-a9YE
ziW#zvD(kLX#u@vUDV8i%UOM&6!B%n?e(i!v5U9TpcLgK4L^WQzSJ)`%s=zybe{PKm
zNw<~9Rb<(sMtUb3kOHHO5t?-HwCZ`mtCy=HZpNv+^?my4bSI=j6$cp(-5(RUT#QL(
zmgUMp9(Iy;N$gOEQ{=n|ELNZt!hY$bIYZ>YFPm!#$CmRYzC<R9=7}A@6v8U;`ZmT!
zI7MLr@&t#A?Jv5NHgMQyAdN=WZH!!kscFuYSzl(bMH?5-o`!hIT^!68Rw;k`itLiB
z{%|FK*3`J+Ox8N5f1|^+t3i5SAotLEPW|v{EKBp+d@P1iKVf=ld5zgq42o;vvd%pj
z#FRdJaPc5hhWjb@$F5r&l;5JQ`F7#pF7OG7{R@bofYbHuQ34@??d7|x7e|ia%w}7^
z_<5<~%EVe@;k~f1fuyEqx^_$fQT>hfRW|a)C!}nwvyR2I#3cwSHayID!`omNC{-)N
zrTLB)kt?QFdZ&CWe_iE56lIdZYCdge81~ydH2l_N&AuB?6cEZa0?$QNdVP)Avze3`
zofeZ>N(+rHS#>@(e)DP#vG<5yd=ae+=LhQgy(&Ih-J~y-NYs|oOs^&-x(hJhWSMn$
z=5$;fySAr_-Bj(u@#Vglvf<MBR<&+LL;cSp7xW~}dwosV9X8YA@0|$TD^`Yjm-cSd
ztZtOnGWHnLC+w!y`r?E)?^#*jiFdKjc#u5^reQ~2tYd!0yyo1DRS;tpSJRE*zzCl7
z34gu%S0`_OH2fAgV|^#6l|WUX^h{rQwS~S)qqe&np_rQ>8d@Mxb11u219{7dfpG%x
z>yjP%N@!(3m3Z;hTk>7J?yQz@&YP&vMn0jp(u1AoV(#=_$(_i*hD3_-Mou@MJ@QH^
zxUT)Zf5WB2$fW|GUFLKb>?!>hc1azq0Z>XX6wCyMlE%TLF*q0;Z44?U>WEB0Dcu5;
z(!o**N{Ki}WmjZ}H@}HVDK`UcK_5dE)2+zcZvI7kNQ~RrS$Cmm7a$5pV*XuE>YrDp
z2Y@Kq%CV?RN)V-x5p%QsGg}V24}!;Ox%<PtjVU4oqeG%Oo{Y`xf2@vSai7m|z==BZ
zW*Yimf|+n8h#r0Cay~pD?TAE)M~J=duznYt<WXkHqN^FH?0z@iWUOhX`y}sEACaV0
zne<NCKtxcUIGCv!H1Ao|H~Y!E`qL)Kus2G3mjUg99RF(L(KpsUM+RuNMV8M>y!NTB
zC8}Afi(8We^if!H>gtr#3z?R*`{8GJBox+@d5O|rIZhT#wOECw^iOSNF7lZeuj5Lk
zoK@^%I7adrR~7QTn5=5s-eM)LWD(mbd7^HRBlGvTglQw1N}s3mG8@nZanf5CcHZ!E
z?s{rHCZF)o_lt%bUg*re^R6ln4~EA;UXflqD3ON2Fe9vn$zJd4-)HSOS;k=;7b%`L
zK^<6N-d#R4SND7(SXx!`-CSwd+s_{68IDn7x*sh=9HJC6%d&#sq^N{G(U;>J;-)D?
z$VrM{Ki>V)H{$c*LEP%s{<U?k>iu^%ZFBzQAZizM_gB`ldQxA&UO1iVN;#<@ES%<J
zW#9Y#ufA2E{V=W^dQgSV*|4+5)*`?uu<LxIBy!&=P4+jvpAjyt_bnm9OJsD(u@z&5
zpM!hY%RemKcdoipU<;MyV;V}CuTbJYA$zSawV?s|HPD`CnRkVGOH6Fh*J+Pb`5UDT
z9M9W(=ckNcX6K(56ZhgiGHr3jj!|buQK|%CFxy6E{iH3uLScI4Fn6%*yIlsm8%+EQ
z9E3zPf-n)ahv_UUV5GWHY#XWVWwA{vQfyju-2;2-=)&s-b0TBh!!sd13_TE5wc4y$
zb~Ub?J;{S0N)()7?car-T~HzdE%Wc|Q9W6gKn|f&Tg!{&58qsnPu0Q0o3=3fhFvki
z8Kqy(=<pXaNh}^X1`F$dcgMBi#Zj!e9LL@R>LQu_W=plC6>=usSzfV6Ub)~GZj8%O
z^l%K@Var?tuEYnK28sFWjii4xp4|wxN(<AZ$J0EhByxJs4fu2&6xO(S^$q=6EW}XQ
zM6>5LWNMXCI3X;(Lfk*s`_q^%jGp|f3UlT{oYE`0h~7IU9~PLF=jVeLBOLvXF=|aQ
z+)m~_$Q2o!YrFMyVxw&<;H=!pb@nRGm8VBDCxZJNoKskcd*U97mM(fpSx>UdK1xOQ
zEfjw_JKW@$(Q;aEW#jO~$mefl+@fxucVXqv7<cFB%8xTBP?I3%!0p2<8h!*H2vlX~
zpob(JhJZsPQ5YlyE{#P2XEl&W;4lgFyZL~`()_E}Z5{z`(ipVa?PmyS3~3#qLxP{k
z51|#<-kAY^0Ly!t6T8s53%rB_Wc+J?>Pd1;WW02$#|`rN*ZIj5N8Avim^D9_2i~}j
zaZAYQ7Usx%(~+@Dk<kr;vwK&rlDoTXzQFjRnHI={!p#fkJ<@PtO+B9Jfp0x9L~+`y
z%v^5tp=npxFEhxpJbGQQKq!7QyrJs<(c$Aw=UME<*ogiCE~+~(O&)Gw<O-Ub`@a?T
zJ*=EODZXsx<Y}clNZpv*;-*CX&;&bOp>UXoM)uP#W8DoJ!v8|V6^xn^$l-K@c2S2c
zJ3?QVxwFRU7H@sf+W5pb#;T4JLuaLkLW9x-txO)tp19CWf9&itUssHvm^w7IMsRJE
z=51K)^8;(M_b=3`qB1<NCuxyl<s4k&x|nQ#G3MsW54(_`_xS{pnMZTi+CWAN9VZ(g
zgueX;>275)$MMnW=u>P}5zb$-kK<76daHUMIfec?ddu?bJL0}-lk>vKzJBdb?Q_cL
z$rw3PvRnD>ENi8h7(LO-`Xe3IlHo0l_TM<_`J3sKKE-wiWaAR^Z+iI`jX<9Xn>^Ts
zMY~`>nZJ>Sc)+Izqf$Fw>_~o6=o8@6>wI)`AA_Lp4dGb@TQ2G>QD0=2Dn%Zh#96z5
zhz}QFw6>Z?iMh*HobI>3Bv1HQjC;SbKP_*$^6A-Edq`Cdx$Q#QE)WKR{0m^OZ`AJt
z%-)w48Z@k9iL-^9GEV;5*Hfwk)F~Mrx=-`prP*8<z_0OR{ra1qsYV;5d%p7>y;EKL
z?Ft=f@RpG&!&Sjrogv%1zO^WA=v0?XLWbA%#Kp4XXRNM-tj<3e7mawKfVkyW*M@~v
zC=Aml8QkctemO3eE<hIBob+C$h~^ddJ;8$W;xn1JTm#8x2Qw}$9Y2vsqR{gw;o(@{
zAX~h~b;;7PBk3~xVaz6D=>ZkV5oyd6s_vQ=DB7ai1EDJD&ID|JR$G_9BiUfxnZly0
zxYqsnJxio^)3Wm@oybIz&c`D0A3sx+Z`w=m!nfUk76$e=D;&!N04;K2Tcog-ZRA?z
z;PqRdx5_9dGz(skqBxsb^#?)^cqqQ8)yQui^;(k9Nm0ASn!;sN;u0e#qNS(hlA_hX
zeVRYpvd%6u%T%lm_oUvlr}O5N_5Isw3H}M2u`ei0Hw|CfrmgCqu1<k}K4{u1cRF5s
z-Cq5nHT?BAHCQeCCu%=+p}l8fYs;qEM;y-%t<g-%^u6Dw7b*Q}mx1mAz>vtl0J)Ml
zT|y6pZiYx>cNw`(xgh@{usNh{t05&-U4!271N_<9Se+AAkj5($-Hi!DN>NX$D~6&*
zv}SJ0#tbIK-AIv)I<6+M0PmLXj$EmD%gwz}In-ZA*$bUIx?lX$y)$^;`5yLZCzd`Z
zmIce2Al~xhAtr@gWdaZ8<F*ugh?3-zCp1RG9TqNCC(do8uzId#`ixhF*67KV7uIlk
z8xrk>9IV}&WVuHu`#)7ZV8peM9*Al}4ile<jGJrJvj$_dj@WwpC4RNk8IAI_)i}qW
zbg7I_{&HdekPyw&^EuoEF*2u@*;RoA4zqxgydm_(xJIP5+pBi1+N{oFb~1I}=uF(M
z9}^xS6x=JLTB&u{9UaA1<+WX~@(8@h+*){Mq}RaEpxpxbiWkd5UFK|Uam({SmNWGX
zlY@z|Np>jvw=eJByp=ooZI1-;1o6*lLGOa<r2j(vkVuRP=no@hDlR_zi`#oI^WdJk
zA&H>~34Yq73!GA_v@S3oXYWEsH)jf0D+nCY4Enz_vzf_up=KA@g+%-XIcUJR-3$1?
zM?`NZ6|EW$zb`M3YX^*5<~9~RnB%-C)atHo`#r&#qOJuSo;AT*l}G1aXAj|e)>!Kf
zj7S0>R2ds3X~>P46)dx$EPVGyJ}k+W`s-YZu%K!@z|neIo8<J7bxQ6cB*bV@skg-$
zfg_6)ctG?Fpo%33#hxl`&3i0-wJ`sKf>qYw<5b<EQ7(Ohr9C|ZJ=&M;v^OqiW$LJ<
zWV{NiRrzqXBSl_O6W$-bjEVIr*ekTiT_iZveL%LRkF?b-`@Y$dV<&gCXuJh&T`9Ns
zI~(ToLawXx>{GBugD?~4{3kVg?>%?MN4|MBkFM^Fd2O6(*X94_!kwd3V;<i)SO?}@
zSW?0@OKO?+dCGLebmZB~^w_^=#?p7hxa4L_(7AXBm4qh7c8uoL6*=|xaa(_Uzu)x5
ziRfJ>Vi&ymH@X;!sW*a!cO6izLXmJiC&TU1^3U+=$Kh81PyIWjIj>7k4t#`L%PC0g
z=;K)oK)i?|ZoLlA^<{QZh%$|6ip*D#@*>M+;o)U(c-vTZAyM5*(KHlizQ2Diy0UkW
z_u*xU)ktlzv%E%enTL`k@xG1WGnjqXuBDq4CDz1Uft~+U@6E#PP+-FslydlVw;1%{
zd^6YWqKTGA1##`C3irUNcNxtt01biu3;&;_u{MEBnW|+({KFZylE*F?;?WCRlT-tg
z6#L?o0<E83j=uzvuxNdjT;b967zvMjyf^x!_@s~Jf|4mj+r<c-AN)4-1NVb~OzNT3
zi5fb$%Ke-1q7RpPp1|4N8z=Gi7?$m-6|$83D{Z?4nc4L2JMp{IAI?HZ)e{w@A&VN8
z8BR&Q=7*3>iKr+o-l+16g~Q~^Oz~A32^hkGFGbMvfoIJUO5w{yNgsz1eHpLOWy{KU
z3u)YT0aC(j&`Qn^A2DQX#-5e$l(ag&`tJWoyq5YnM0atRb$tH#o&I;7rj9)qCw-ML
z_n$s)jhX)vU%j}vu}S9OXcMyw4|f56@V{|?=qG8@s{$Z*EvJ*M^bE3NQT}Fmv&YKK
ztE!Yi5t|DTe|h=Q00F5$*QNMG@r#`E9d!_wr(-f({Wlc5;mWfGl9n04eN#hB#7+d9
zv_YePO1_|Df1fQClWPV~dk2Ki?)--Nd#ltBBaL#i=bP)@925#VVwb6UYW7%!e{JrI
z(Cyr?A79826LevYei!kQ4~LP<XSi;<OejBBH+xaspXI~mnChF?d$04$Ki*~7yMQ_A
zzwodq#LLl#!?s7Q>2XqkQY9wh(WxTQ`nazj1ZH~19IcZ{8$0?J6wAg5BLZ?K7%uWI
zSa4WhuNt_UO&6CPe{<jCS^Y`cVi~rMU1-`3#r@UyM`J8KU~eq10jCNC;-Mji1tNF8
z-MkY&X3l-f_H2C`Z`M0pYANqC6UERhshsAAlS4;Pix=u8iQzXITA^1{I%3^oPT%s(
zys;-+evj(pSgElp8~5nFA`_9ChvbUq8h!WAx?GRQXy9Fpg*c0_ODNwdQ)g!N84)f-
z-vp5Nqm0+Hw@r$AFC2`Z*kdlkIk_~snyz3dG#1E$xL1K5JV?!~&L_W-SbU%Pl{M-v
zg6Uo_W91RDR=KY*0}>MEfI(8%OFC(F3*COTAL=M9K7ZSDNz?SgpZo#Lz<)m0?AVUm
zu?$cZCnEm)oxZWLnz^_d{)m%#kcpEo5!~4b#6VRKCjw9zs^|3m9)U&!b|VlmpoS;9
z5Ip?@{Rj|1)on)^{KpOj`vCSv11BFKG(qL}rg!T9i<Js?q3;w(@Z64Bkb=R$F#!Mn
zxdQ_i0OE{*PzoSy1Kjw5I16=OAoKtXwD>2$<3|<^RL#@P4G)A#fN>Ev5G3FiA^JPH
z03jQIwHAse_(7n4z5u;}NBtmZfO<F2VBm=e5CA3OgYmu)DBjcEBM<_9UdIOnIe9}c
z(m?zO&@Ke34gta`pjr@cCFnuGa{-1B(8Mu;$N*(5AXuP|Gr`Hl3m@o>cMF7Ifp_1Z
z3EPcgfhsOEP*<R7pl2ZV0)Utv5U3|`2fP7{dZ2rNKz|4v4!j`%Bmsm%ph19D4|E&~
zoUj4HEua>ju7Msv<OT>@e~#SPo}eFJ>;4gdpQdnoO8?7z0ixjT<W2-Q_E$5fOmqQ5
zQLqSraIgpPK@tIaVO5>{wRU(2qk&rfxCH`TfDjx{7iC{}Z#+Oa)Fcq^V-DWvJB4hI
z2m(i<5I;x3`~c-P5DX0V|0B`-zqW$}-~|ehO%!1K$RK4P7$k-U3!qbJ85jg1jfNnx
zu<e8f=mg(`<&l5{ZC{bnVEO-BqWSS236upB*v5~92EN$AVITu!P(T720^)h30T_${
zWE%~1E)6C$0)oJzX@CU5pfM0M8fXOY#Yh8EFAWGM1_s=NSMWV>5(@)_7)b;5!8Wln
zSm3&o2VS=agOsKLuOJcNpumKIBe(0Kfr0O|0h9&$-M$BmC4{uhc0vOS1H8uqY7ztF
zW8es2fII!7fPT@yMuh}Ey@T)2z*Pod7I@uG2(TPxn+f0*?0<XGKn4Od01O3Ap8;io
zI@?V7K^%D9mLISu6!-uGSp<>|4h_5lX?s4v?sgIl8W`esO_0C;+yDgbGyxoB1FNAy
zssPHt>Yy0>n64jJV9Edo!A~uqq=7=faXFBDP!4_}70lcgCZHv-)_+<7iGsmsw%>hk
z01JF=|B*mJ*pUHH4t`DByawa~l!>3xwmg6`@gr@^%6FLn<>BW9(zX2zv;oS-j|7e!
z149FI2+#`3z%B_Z1I!617l1zgIsK4{A6Jkmz`WCJU%~sI()L{Mya&ttCxLPS&MJ8Q
z=ZN9YHh?N*`@|s**oh$9=jebtbqFwN+xkiif&kGBkf#U8<<6ZU1o1s)&II5Mm<MdD
zknJ-fK)IjyaDeICb${8%0q$&Px&VT@bLR>HDBf-c4=jqEyq!%DklgLEJIVxD+}lrX
z5P;h|cc20SoKf2k9)M(kcQk-b0FDlA=lMW@ZE^d~7lPXE1XL2J@120k4E4PWe+cS(
z^8`@cZMQ>&0AH}TYwxJHZJQbxBdGBJTVUtSwt@tf)%H^$s26wcKppvgSV6%2?&O9-
zfDLB*{kA#<b{g>EpJ2+LZ72;mWE+<J{QhUFL>d0G4*yShB~(=zXcX$=1X$#t##2W8
z_hM3j%K4`t^1l`Ybo2LOKhX*f>AS}KX8_>n{>TJJ`9rtrgL@6|46sojSlr*|L>~yX
z{LUhMz#0FYNBV$E{`;ipqanZ&2dV)I7B2Y@2$S^j^bG>$X{X75=-BS#N1r?T{NFwX
z$Ur0k%>n%X^92BogaiHk%fVZpBxU+D5IXfW1$F4Ewyrw<)m^U}S$a>uHj)J>>CK`>
z+vQIhe@j_4DC$@jpw>00TL>Lby4$Q5fAIPJYa)G@jU^|BgJeoaM*Up^<7y8S-<wzp
zxpMIZ=27b5&ox%O8izUQtIr%~RTZ)og;PX7riIRQZ87HA-|OHhNI2<Yt@`y(IJ?c(
z|J~95B$?l1nE#jN{--7QKREiI_nQB`bq{Wl|Fi4<XIQf{^*g^<0M+$$1?^b?i-CFv
z*xptBf&i0iTR(mey8u49|HC;|ASNHQ+C>2K+TYtR(A(1)5{#6>NFgL39)W@WL|G`*
z=RdEc{0Q!%G*H0Ib`5d?Li~TO<?reSaR!ZAU)S#~042AZ_w@5s1Ah%cL`@a}Lm*)Y
z1OmX|Fa$~*h7g9qghhYSv>^n%8w~)Fu`~eI{{9D;DnKL&#0~PFJWy+be-PjAd2kFG
zfT_Rbp`<~u^=mnlH0F=GC@gH3JP-r^Ru>DRrC;)Z2Mn+U|CWaY;oh%#NGuSf@mn4Y
z0sEspBo@@~ztn}JVStJDdma*m3%`{^1LnZ*c`_ir{k0qh2AE#I=lw}D9EkvIqu<IQ
zccBIF0{u}Ag#@woFa5%hSil_lEf0kN;QH@*Xb@}vS{H=@Oo!j|Wd7tc3j4?W!T`JM
zkABe@%pbf*%YcCW*M8Ah>EHPb11z0C>S9oTju#{ShYSJMCWy0t?L!*<$6NyE;Qp9%
z7+`h%L7xm9ILz@|don1{8u_&k8PEmrd)}XN55oez{9X==1VPrX?O_2^?~gob!0+*E
zIl!X*gJw7k`A1#A3<pfQ-|Nc!F>i1<a60|>a&XWW@>^Xv`p<O)mjQfqzm-FPF0|kC
z5STw`!C-#wH;~}u>5V6B|8k6pXDA-nuOLvIpC7o<ZsUDzUpK%e0yw#2@(=@^2!Y#M
S8yq781Ac%;P~A|2=Klah!-enw

literal 0
HcmV?d00001