From dd6565296c9abbf8009a6984e6832af147b319a4 Mon Sep 17 00:00:00 2001
From: riinuots <riinuots@gmail.com>
Date: Fri, 15 Sep 2023 16:01:20 +0100
Subject: [PATCH] {} removed, multilevel variables lumped, more numerics
 grouped

---
 R/app_ui.R                                    |   2 +-
 data-raw/appdata.R                            | 146 ++++++++++++++----
 .../rots/gecko_shinyviz.dcf                   |   6 +-
 temp_histogram_ggplot.png                     | Bin 0 -> 24151 bytes
 temp_matrix_ggplot.png                        | Bin 0 -> 43708 bytes
 5 files changed, 118 insertions(+), 36 deletions(-)
 create mode 100644 temp_histogram_ggplot.png
 create mode 100644 temp_matrix_ggplot.png

diff --git a/R/app_ui.R b/R/app_ui.R
index aa5c561..16967a4 100644
--- a/R/app_ui.R
+++ b/R/app_ui.R
@@ -18,7 +18,7 @@ app_ui <- function(request) {
              
       ),
       column(8,
-             p(paste("Preliminary data updated on", format(Sys.time(), format = "%d-%B %Y"))),
+             p(paste("Preliminary data updated on", updated_date, ":", nrow(appdata), "patient records")),
              mod_render_barplot_ui("render_barplot_ui_1")
       )
     )
diff --git a/data-raw/appdata.R b/data-raw/appdata.R
index 992f6da..175eec4 100644
--- a/data-raw/appdata.R
+++ b/data-raw/appdata.R
@@ -19,39 +19,103 @@ patient_data_orig %>%
   ggplot(aes(age_years)) +
   geom_histogram()
 
+labels_keep = extract_variable_label(patient_data_orig)
+
 appdata = patient_data_orig %>% 
   select(-redcap_data_access_group) %>% 
-  mutate(ALL = factor("ALL")) %>% 
-  mutate(age.groups = case_when(age_years < 18 ~ "Missing",
-                                between(age_years, 18, 30) ~ "18-30",
-                                between(age_years, 31, 50) ~ "31-50",
-                                between(age_years, 51, 70) ~ "51-70",
-                                between(age_years, 71, 120) ~ "71+",
-                                TRUE ~ "Missing") %>% 
-           ff_label("Age"),
-         admission_prior = parse_number(admission_prior),
-         admission_prior.groups = case_when(admission_prior == 0 ~ "None",
-                                            between(admission_prior, 1, 1) ~ "1",
-                                            between(admission_prior, 2, 3) ~ "2-3",
-                                            between(admission_prior, 5, 10) ~ "5-10",
-                                            between(admission_prior, 11, 10000) ~ "10+",
-                                            TRUE ~ "Missing") %>% 
-           fct_relevel("None", "1", "2-3", "5-10", "10+") %>% 
-           ff_label("Prior admissions"),
-         post30_los.groups = case_when(post30_los == 0 ~ "0",
-                                       between(post30_los, 1, 1) ~ "1",
-                                       between(post30_los, 2, 3) ~ "2-3",
-                                       between(post30_los, 5, 10) ~ "5-10",
-                                       between(post30_los, 11, 10000) ~ "10+",
-                                       TRUE ~ "Missing") %>%
-           fct_relevel("0", "1", "2-3", "5-10", "10+") %>% 
-           ff_label("Length of Stay (days)")
-         )
-
-# appdata %>% 
-#   #filter(admission_prior < 10) %>% 
-#   ggplot(aes(post30_los.groups)) +
-#   geom_bar()
+  # remove REDCAP internal varnames that are between {op_code}, etc (matrix variables only):
+  mutate(across(where(is_character), ~str_remove_all(., " \\s*\\{[^\\)]+\\}"))) %>% 
+  # remove HTML tags:
+  mutate(across(where(is_character), ~str_remove_all(., "<b>|</b>"))) %>% 
+  mutate(ALL = factor("ALL"), .after = 1) %>% 
+  mutate(
+    # age ----
+    age.groups = case_when(age_years < 18 ~ "Missing",
+                           between(age_years, 18, 30) ~ "18-30",
+                           between(age_years, 31, 50) ~ "31-50",
+                           between(age_years, 51, 70) ~ "51-70",
+                           between(age_years, 71, 120) ~ "71+",
+                           TRUE ~ "Missing") %>% 
+      ff_label("Age"),
+    # n admissions ----
+    admission_prior = parse_number(admission_prior),
+    admission_prior.groups = case_when(admission_prior == 0 ~ "None",
+                                       between(admission_prior, 1, 1) ~ "1",
+                                       between(admission_prior, 2, 3) ~ "2-3",
+                                       between(admission_prior, 5, 10) ~ "5-10",
+                                       between(admission_prior, 11, 10000) ~ "10+",
+                                       TRUE ~ "Missing") %>% 
+      fct_relevel("None", "1", "2-3", "5-10", "10+") %>% 
+      ff_label("Prior admissions"),
+    # LOS ----
+    post30_los.groups = case_when(post30_los == 0 ~ "0",
+                                  between(post30_los, 1, 1) ~ "1",
+                                  between(post30_los, 2, 3) ~ "2-3",
+                                  between(post30_los, 5, 10) ~ "5-10",
+                                  between(post30_los, 11, 10000) ~ "10+",
+                                  TRUE ~ "Missing") %>%
+      fct_relevel("0", "1", "2-3", "5-10", "10+") %>% 
+      ff_label("Length of Stay (days)"),
+    # symptom onset diagnosis ----
+    pre_symp_adm_day = parse_number(pre_symp_adm_day),
+    pre_symp_adm_day.groups = case_when(pre_symp_adm_day == 0 ~ "0 days",
+                                        between(pre_symp_adm_day, 1, 2) ~ "1-2 days",
+                                        between(pre_symp_adm_day, 3, 7) ~ "3-7 days",
+                                        between(pre_symp_adm_day, 8, 14) ~ "1-2 weeks",
+                                        between(pre_symp_adm_day, 15, 28) ~ "2-4 weeks",
+                                        between(pre_symp_adm_day, 29, 182) ~ "1-6 months",
+                                        between(pre_symp_adm_day, 182, 10000000) ~ "6+ months",
+                                        TRUE ~ "Missing") %>%
+      fct_relevel("0 days", "1-2 days", "3-7 days", "1-2 weeks", "2-4 weeks", "1-6 months", "6+ months") %>% 
+      ff_label("Symptom onset vs diagnosis (days)"),
+    # diagnosis vs decision to operate ----
+    pre_diag_dec_day.groups = case_when(pre_diag_dec_day == 0 ~ "0 days",
+                                        between(pre_diag_dec_day, 1, 2) ~ "1-2 days",
+                                        between(pre_diag_dec_day, 3, 7) ~ "3-7 days",
+                                        between(pre_diag_dec_day, 8, 14) ~ "1-2 weeks",
+                                        between(pre_diag_dec_day, 15, 28) ~ "2-4 weeks",
+                                        between(pre_diag_dec_day, 29, 182) ~ "1-6 months",
+                                        between(pre_diag_dec_day, 182, 10000000) ~ "6+ months",
+                                        TRUE ~ "Missing") %>%
+      fct_relevel("0 days", "1-2 days", "3-7 days", "1-2 weeks", "2-4 weeks", "1-6 months", "6+ months") %>% 
+      ff_label("Diagnosis vs decision to operate (days)"),
+    # decision to operate vs operation ----
+    pre_dec_op_day = parse_number(pre_dec_op_day),
+    pre_dec_op_day.groups = case_when(pre_dec_op_day == 0 ~ "0 days",
+                                      between(pre_dec_op_day, 1, 2) ~ "1-2 days",
+                                      between(pre_dec_op_day, 3, 7) ~ "3-7 days",
+                                      between(pre_dec_op_day, 8, 14) ~ "1-2 weeks",
+                                      between(pre_dec_op_day, 15, 28) ~ "2-4 weeks",
+                                      between(pre_dec_op_day, 29, 182) ~ "1-6 months",
+                                      between(pre_dec_op_day, 182, 10000000) ~ "6+ months",
+                                      TRUE ~ "Missing") %>%
+      fct_relevel("0 days", "1-2 days", "3-7 days", "1-2 weeks", "2-4 weeks", "1-6 months", "6+ months") %>% 
+      ff_label("Decision to operate vs operation (days)"),
+    # anaesthetic type ----
+    op_anaes.grouped = op_anaes %>% 
+      str_replace("Total Intravenous Volatile Anaesthetic \\(TIVA\\)", "TIVA") %>% 
+      str_remove(" \\(e.g., midazolam\\)") %>% 
+      fct_explicit_na("Missing") %>% 
+      fct_lump(8, other_level = "Other")
+    # end ----
+  ) %>% 
+  ff_relabel(labels_keep)
+
+  
+
+appdata %>% 
+  count(op_contam, sort = TRUE)
+
+
+appdata %>% 
+  filter(pre_dec_op_day < 1000) %>% 
+  ggplot(aes(pre_dec_op_day)) +
+  geom_histogram(binwidth = 1)
+
+appdata %>%
+  #filter(admission_prior < 10) %>%
+  ggplot(aes(pre_dec_op_day.groups)) +
+  geom_bar()
 
 # %>% #select_if(is.factor) %>% 
 #mutate_if(is.factor, fct_explicit_na, na_level = "Missing") %>% 
@@ -67,13 +131,31 @@ appdata = patient_data_orig %>%
 # levels(appdata$diagnosis_timing) = str_trunc(levels(appdata$diagnosis_timing), width = 25)
 
 
+appdata = appdata %>%
+  select(-record_id, -age_years, -admission_prior, -post30_los,
+         -pre_symp_adm_day,
+         -pre_diag_dec_day,
+         -pre_dec_op_day,
+         -pt_comorbid,
+         -pre_img_finding,
+         -pre_img_finding_cbd,
+         -op_anaes)
+
+
+# pointblank::scan_data(appdata, sections = "OV")
+
+
 allvars = appdata %>%
-  select(-record_id, -age_years, -admission_prior, -post30_los) %>% 
   select(-contains("_day")) %>% 
   finalfit::extract_labels() %>%
   select(vname, vfill) %$%
   setNames(as.list(vname), vfill)
 
 
+
+
+updated_date = format(Sys.time(), format = "%d-%B %Y")
+
 usethis::use_data(appdata, overwrite = TRUE)
 usethis::use_data(allvars, overwrite = TRUE)
+usethis::use_data(updated_date, overwrite = TRUE)
diff --git a/rsconnect/argoshare.is.ed.ac.uk/rots/gecko_shinyviz.dcf b/rsconnect/argoshare.is.ed.ac.uk/rots/gecko_shinyviz.dcf
index 04ede41..b76661e 100644
--- a/rsconnect/argoshare.is.ed.ac.uk/rots/gecko_shinyviz.dcf
+++ b/rsconnect/argoshare.is.ed.ac.uk/rots/gecko_shinyviz.dcf
@@ -5,9 +5,9 @@ account: rots
 server: argoshare.is.ed.ac.uk
 hostUrl: https://argoshare.is.ed.ac.uk/__api__
 appId: 639
-bundleId: 3265
+bundleId: 3273
 url: https://argoshare.is.ed.ac.uk/gecko_shinyviz/
-when: 1694698060.84986
-lastSyncTime: 1694698060.84987
+when: 1694790034.70899
+lastSyncTime: 1694790034.709
 asMultiple: FALSE
 asStatic: FALSE
diff --git a/temp_histogram_ggplot.png b/temp_histogram_ggplot.png
new file mode 100644
index 0000000000000000000000000000000000000000..e315db3ac12bdea4b9208b6fc77281b492ec4da9
GIT binary patch
literal 24151
zcmeIad0bOzn=Oi3YL%!|oKO*{3S}W8L}iQ&(K1LG1QcXYCItbRhX7$rL`#W`r64i{
zP%$7uKx7V(piH6=kSWYTPyz%JV+d0qId9Oa?(d%aNB2GVobLWk^ZOwfc6RoDhv!*q
zJ?q(b&zl%-+p>R)goMPlv%e$FBqTNtNl0w?<NHnEm1=EX1NgDY(Z~=9euA&1)73ce
z;|Jg0FQFwQR29VE>w<}hAPI?|B+epFT?{774C4ZZT_Y41+j$QTh)@LiWI{t;{rx|?
zw*Qck1D|OA(`%>8(ZJ`*J1#@sZatEyS30a}mr`!yvS;A^$=tzV^D8Eg1Jaa?-p%cL
z*7>7R{P3m$ZKM9*7LO&Ftv`n{Uw<rlS7OuIBQ|7)s|f#bMEjTI0_L=`L#;h?C|jT_
z#3I_8UlwwDY#j8o)~{Xif80OXYQ!8F3CrZXyu6rv68XwiTRMMcbaa%SKJ!Cp|EpV@
zVH^&}fzDvyMu*Vtdffqtj+mUB99kY}x|3FQI%@kcu5q{30WXr`r1NGextbSLD()v@
zz~W6w)kT-C)swKd@LX_|9u@7~mZrhzFLmpxH|Bp=JEOx`qBUuF^sY!j?rHgqtUSmO
zTy@S*^Omr_WQOu7U@gCstiUfKyYx=M@=4PPEVA&+%*;%i@xuDr!_06mCfdf1?w~JP
z;lVt}ut^IgBF>!7e>gZhAC}JE7}|fg-ODgRE+PWFzi67p3j46)G%IPm@cONy?`H96
z-&p@0du`Ou->8#a`(?|Il^evroLsjav-ZNFlcU?$e%bi_&^qxik0c^>)?V24(?8cq
zIv4>PBwp-=$q(0-A$ke)0C}R`QT`%XBf>&=Za6gt6VHyocC9AOzsYG9vSB!3k6{{H
zJXp!M>SS-I{YF|I>5pGd&Ji4%^{R*n`M`!5?xqf1gAMz#omlpJAz5BzDR8v)_VT2)
zf&5QB_`e%#o%BgfOIx1BR~Ap{Ti-W^6hC$vipS^O35^`b3fhJ9L)69Px#b+eX!lB#
zmEF>3o~mwGIXINxd!}JYvDMhH!FEIK$@X{ftVFYc6tsSt*|Im-3sTjSe{etGgKO@d
zc|<2B(cZ!x{xo|FXWlYiI1Qn1+rAxJAah$#eIdzWwkCKvQQx(k^4aI@mJX|>ecaS%
zsN$+AxzT-F)rTeKNJ*Wqo69}~{&FN2J>hfmDx7&cK}olQ{|AXgPL7crQaVhW%f#>X
zAWNAz(3X<%eLcrCGR&5za+V`EF35+^;NtKPHHJ6mPIzqIPUU}*Y_Fn%%O9R53qK@v
z4xV6w9T<=I(dgbWS4u7vS)itJ@Qa?PDSA44bLa(>Qcqxn4bzskv_K-IT^H|sC1O5!
zuapw8A0KSowR*nc*6`sYzPZmLBO}8U#zw$M%^91Tir&07^=8WM+q<{2HV=E{RXk><
zKk!SowVUv>C~TBX7sbf0cy>j+ARC*Q$Suir$;*QyefrZQEw1?cf4R-WNCnIOFnkYf
z($&>vnrfsJ(c-k7?Bwn~T@8W2nW4VE9pkaCpN!4_xi{-SS?u{=xt$xr9GOXVb@e%p
zIm(Ac&<UNJT*3?o(#MBwG1snre6C;^Ka^xPup0u4!M8#DjuHc`?C3jL^dg3ire?3b
zaoO5R(biC&*s}29!v~J3Qrww_w4L6=_859=Yij^%N&&(p<12rG_8K&x{b@eul;YZ3
zTXnnL%gT$A2i9@{fxH`1nGhR0_lcmG95xJ2qxFfffNX#(!UK$uLpti}-P8PUhh~C9
z<@G}69zw^Wmp*KM=+tJJ&mSZ3!gO-BBO8>Vj^hmtNDVD5j5jUh7pR<__D@80iU+;$
zU1*u++K1&9arK5Uvjc&bXiM~oC{6)sS|QL*#8QQ&L&vgZ$9-~=3ir#)BWjKo`!Fbm
z3S;Fm-jyF<%gf7@J&<%LltfxSkF*`U_4UTBZ=Nagm6w$*Ad&WG$(zXQhc7-Hm<?jj
zlTV*M4U5McNjE7GEnTk<yKzQp8F{Sb9$9mkSJ~zE*|M)Ubi_D2Qouxrmb45T9`Ne3
z`lvWjtN-(Kvrwbv&qdA6&D1&1W^f^^;}Gue%}W(?m%Y|*BRQwIuC~YBZa;*RDUS_y
zbR5fczI?)<A+C9|QZ472Vr+bTFN7<+?MGXhkjvjDeox&FjUAN<@$r@rzvh+Ccau6z
zVe|oTnKlYz&*Sq<VG}hZ5}Uo(zSyZbKCaeGM*N-|M2}y5MMF&e1#^fWf8}#od#O~2
z*3+c&z<5_zS9$B!kU%?JExP7v{he3#d0#)KQzdLvUstzAgE%)Rjw+_-&Yg?Rnjzfx
zr%7Yt+w}!A&C=+Rp`js#jPaoA)zZ=Ba|Pet$%pYmm6+j!E2}}(!&oom$4{P*NF*bF
zeagp~zCP=)V@$8X^4T|$8#c!$plQ;4)a-qLq=Xy-kLF!?3Il?wA1V(${iLO(h3ZX`
zxBg5qU70<YMG<|-864eJyV|6TSb3b4s@oYL|J)nSM4JSA#}IgIaCr8%vhVa%hyhoT
zNZL|pk-}pgaN@BBW;uS5p{S_n<Mu>;G7!#;<0Vv;fGiQ4UhdiNkJH)$^BC0!yXOu*
zxMk&Y?!xS#RvCe1i0lrmrHr4Qdf%8?5y<Kda`m;d=pozZ+4Yqtbz0cj-}+S6L{(CL
zH6UynVWuCn?y|BI>QVOn;mqIaaLa0O;14!RDPt$Z-thI8Ds5r_w?%40@=q<lhn}{u
z@Tj(|wz+K@&}FRUbyjKhAs(O9$t<dxNRv{f4`i63GH}c;3Yk*K%8g*ty1y;5VgK7j
zZr<MBN1)44=)}YXHasAp)IO)obPn0PY8<2ICK9fiF4<aG@N2Pp!oQ|thswL;;{}#1
zMTIfLi6~JJSi(5!0TvbpXKs+C>WBBfdGkBNasD&|AkyV{y>10!zzZvP=^#^-Th2R=
z4m(&Y96@zj2}Zd+FZRANf{{necj7CbOzCTDB*>)>Cz_>slcg+cjr?ga4&isDh$B$p
z<Hy)Xbh|_}_7qY##ca7GCq2^s0UBHXgL#&*(yH;!+LJ{_N_H3N`~+gs0EFw1&W<9!
zgCjZ%2JEDUFf)C{$DmjrY_NzN$iik1dSo_K=xhcccKAm~lxdI|YO01dN-kuz$dT+?
z0Yb;c<3&Sl8k2=nQ&XZ{kc>#dV~}of!xs@Rt~Be+hyU~)+m(QsU7gkuezdxS7!<FX
z#Kc4$a}*45b=hmfh7AiEoj=Z1g7xEIs#N<ERG3b`*C7jMs_W^k-r1Rh&^>zebLOwh
zk@i{|8i5hW6>D0EH^{MQnRXH@Qqb?3RqGZM6x6ozIYdoO&G2$TZ%IZc9d2D6G>a&`
zTyT5Vm7SiFvY2Z$m0MKAu?bvFsO~QvoMj>#9-0Q(+tE#b2<a`bSt!*C*E9E6L*6d}
zjB6q%(h`SEX7Y;0+mlx|+{`bcEDu8bJg?eX^u+85nGTd)pdQoIoVPCO7Lz}TD^qFT
z56NhpIB_{QsYx8Jv{c4QcR>(?$AsDTq;DHDFv2Kx(1*HuQg`2qDk{H+tO?bxI%=Zi
zB@UHS_MdW}410@QEG!TZ$UU_1u@&XZgLX!3*n(B&)Z*jEj~^N|OvGpEcX)h@nrXW+
zx5H>v!|=@1)V!>zzOJ$e!Mv$KL*#8gkX?8_>}7Qw4Eq4c6@lj>T*AXwg!bQJSm*P~
zy48TP%LNJLmkXSD71W4wCqqL+7P-)eokI4U9Vl~fE4z#!qnWfa9$yoT;$xvuD1eTw
z+};3~4J)kpYZ@@>Qi}}Jt<s}CVc7C%I-Hejs3~ck9Mw4dD)>?>uv9i|7>68#$lAR8
zNHCmd4-5GA-=fqMi&(ds6qk@dRfc*;;}uQ(Co$$II=jDwS#DqSb-<>kz%f;l7&DYW
zI;^vyp+OdnMoT|UvX_R)$K$8)*P9c=h`6K76;geal+QkSc~N(sW<XBw+Kh~LX9Tbw
zd@6Ts&!ooU0}eB7K}Kdl$EI)lp{8bY{e_eWLHwBp`Bp~siEs1n>mE9`Wg$*R%f8AN
z!CbjFT#Z5WT|<cHSK_(89;DLBB5?KfDi{hJ{Tl4+%O>Du!v@n<T+M2{4m<c!%bB&;
ztuttgbrDQjhKItK0R*B%tk4|#O5*anqr<wo%{4?#k8|gK&0jnE{J=ik?XQ&X`FH=<
zvA+Km<n^DB3ICbqePiS}93FxOJGn1F{z%I^G;L6)DIq6mgl?^Mlp~7?ccxp(ht6eb
zim6YVPUgZZX+59OlVZy8i2zH$r<<lwG>GcMNJIx~noQ9r&fUi^a+HY0If6+0ixw6j
ze2K>Ym@<8u#P)jn^JGMO3taRLe?2!zSd^1~b!uRC-p4-p!|l!07i>mNE|{B}gY80B
zAKu&%=3Oh0@oTrrP4m8GK1Y+M9#<Xo59;y)tFX?`!osdBgf~H*Ix{Qu-sQ|kDXmrz
zcyQx|SNC-h(9FSO`a2-Fk@GO$AtyC7VlQMcjeX0UxA?va#2f{16iI6@1g3kHorw>A
z9%=uG`H3J_U_y`wV)5-7!%}lE>=a|E;1qPELI=wXm5!#tOq6=y%n0L!`c4aro-)ox
z+Yj4y8|UWcaBtz%DyQjJ;LAwH<NJx(;O0yzWP=L$l)Str*biANZ$G=Hvv2iy`TP63
zqTGC#wy3G4wN%oxbMEE0RH6ZG+PrzQL%&!n@~63lhLZKKE;GSPjIQ5LcZ>sUQkb9!
zJEMC-Tf4Gw3Txe8<92v$KmN1nD*t3V?El}}S$>Ev&d|YZpXc-KQq9OGpjd6<oT)!;
z5b9`-Y(+LiK)6)^-tt*=oWG~%OrQUp#_{9+z4sIhy-`ZD(RguM#%rgTk@E{#cfnB{
zSs59RqMUD3E6ut3HJ}*V_zMfBIZ5nzJjSqr3Ivy(5F4nY9g-jof!)Pt0S&>?#l<D`
z`t?aZ<Q}=7{`vFg!qn8%siaOr=vW}cuO%D-HaYbeQC(k6jdwo3XMl3qeGRPmPFsFS
zKF)lagvEh(W3r7N3&5Mx{b^y)F+$=2@<AwDcK?2NeyDhpOGucszud@HB39x1WFd&1
z5ZI*>m$qJzK8vUmX&%eJ{`xB^21Z$YUR71)#lvAM$7^ivL`8|5r?3E5#h_`N*_0DM
zLW=PCxScz9_5k`(1!vyF2bX3)<l0pS+4bgsqp}h2%A*8}i;G#vhAe_MpkqQjzLMW-
z<Ek)b1L5jujRcPy0wfUu@w*`A3bonU**I;0o}}Fs4H^t}ni%Jl>@@WCvAUCt3Li!`
z#BXq$*3f6^h=>S>7{!Iu1QJPq7lb~4<?Br)Uor4jLO07|bj|2|nG-p*>V=W8iU%YF
zAQfZy%7Ci`Gxxwdy%}Nzl0hzX7(T}D1dMm+U+~bIx`>h5j9WklwlyR9qNY@2>#1@8
zP;{2aQUR(;Fl<neWimSd1$xO;(gTOt2bctT<%|D{F}C_R^MvvC-t6LUq~U$^wjuiS
zOT3#jE@1ME<Y$10muf(1(<#eyjx?Ysv*SPP)(HZ{O!*tEBk$f<=^(}`{xtc?lu(f!
zW|$v*cR-1lJF^j?qa&yTZ(}X9GRC*bXl|WHeH(#7>eQw^fUjEO_*S?;A&{u~78bRN
z%#{c+`he*#Pi9lkyJrx1jY`Dv@sf7^unVVzV*0qcDq5U!uNb@q;huzbtGWS&>rPVV
z0w!O7_#v9I{6`)d#7UvCLfgqG98O?WSR)!;2+z81`9jMc)9)c-b$6G6s61^+=Od=V
zWi=2B2xDVoS<hz&Cj$qW2zz$mF9FSSa{<e8U5FDJ8f7!Jb?$$nvGyn>GgL(8U}*h7
znFj(;pPZ=b3)O3EeAN<J<A?2ASZc>nqVc~g8PMR+CqOjpiBTm_BH|Mg+LJF=06seb
z7-DODTwGb!?Bd43axWeK(Yna-aRSev)~xzzQYXNA8pF1?oH6#un(C=vNAL6Kx9e70
zc@gPalusWGPhl*2x^Q`F#JS#0mm_;7-hMu@Be_eN2~<o;X`hb-K|R{WG$kXWRvAj(
zQK?l<Jcr7Q#%u5@b_%o!Xt_&(Ug-GL@_QYxh9Qx5be7JY1l9SP0-Fl`gQJBcce}rd
zf#pXzHpI^aH8l@5w9hsvDF=s_pTj$-KOx2=6)({FUHTYT_rGSb_;Vlt>bTUDl;3$3
z(_=9(eGgEsurWs=_Cau3xIP0o6CuBUf9QBg&j|z5fFVX{U*bH6Axdhuu&^);efABP
zmzs^JYqa~wqJM(w9|vDSbCd5yU@<I<iEe?`kZ-+R*aCZ1a1fTni;UN5JK=Hd+Ebrn
zg2z4srEU?;a?|(BLI=`4GxExe_V3*r?1ii;PX}C73bFDz49HLHlu+L*S4QD+8sG8Q
zOVh4cEvUB}(0JQ)5c8j;jcdFnOT=qJ6U?_7F~F4xII)$z@2iaq8sh96C68ScwZf_5
z^t@@4w<lSC@&ZWC%MIs7r8jNbG}86ER>!y;cJTuXO~dVh7=nX(v*y0DEPJq{645AD
zFo33(3<mR|VM8MX0%8L7M~KKK@(Y7~;DBCez-c65&Ti2R#SGj;541C#KR@eehRUrS
z-lkj+phefjO7+VfGXPNJ*o>;T_P_;*KXQQ3t3T^;)YT!qLj`ctW`-?od@DJdCs?7m
zCBP@+K~-uzK_x-W#HEqV)Wh-TeFD$bl90m`3dJWd9t5%{g+riJ6iwp-&YTlU{5B00
zjvE$27&Dc)vv+S=IYjU)6vp~qSCUB5Dqp%;C)hYAkWY&LR(uNP>Q*;{R|WeLPHY4L
zF%abz4kA}*7u3WE!Yz><M@$p{k|ID<F+6jI7M**iQGn|4t{9}QQI^w1;h&VVM4ru3
z+qTsnfoj5H%~A4`f2m|`(OZXb)45ppHV}&B4H`lu(Mm*yk_ApMKbsYe-yk_oE>r}e
zsF)ijTMfA1tE5ir=5E)f=4KA(uX7lS{Er+5x>Q*q3qwr$6*?x#>MhXm6)rImeyoBD
zn(BcPF>3@^^1j)CAJ~zYGYxgGvpx};6@UFzf5FoOkjaEI4aU6s&EN%9T+SNt1i51v
zFl&(xBCsC>1Zt;ynLR>FYhen+URj{@$w0@PA>1&CpYV0nY><H8-rJmd`J1+ODY-5p
zD5!jz!@(INYkq}-MEdLY1KK|V_6dCw?@H|y!oUROOBbKRA)ilGik`D;cr?@PQc@yE
z+LG2p?l-n78KGxP3%Lg583#_!x8zNE@PgQ?#{c{~h)BcaK+(9EGX^<*;BvW`#zd<|
zgx~$s)O0g)P<y6c#lOn1-?71n3M_Fr(MjYAkHCWwle}Z_u*Ss4p|X%>=$H_ng2z`F
zHb_dz>+5PCJI1Yca&nqAP~{EX)=?s^#>{#FDKZ}9U`aDnQ4C(nupt0hI^I&-iAX_X
zdp^`hRm3JF;M|Kf+ULGqDcQN-3`~ZTyr8ch!e#O*o)lKI$pAR;UjsBVZw|)^w=5%4
zsJjFIrpM3!P;pN79B+!VQXOiC<E|~Rhov<8#q}`>iF{Mj@7Q1tslHSyRm9Nfm;Cy&
zId~bu-N}>pp^nGv$X$O4S$o|wS54oSQH5p>;DbRvM@<!Mt5^dkb%)D5ebNS0cs>Wx
zzE<iaZuHE?QeS^9by}|fCoN?D+YqX5z{M|rc>%=Z_+PWodlJ{Rt;1%!oNn+{Bw<+0
zcdHa*1ybPuAqp8Mqq0_4VVUS)m+~o1Z39@%%nC=adTX;3sQTrPG%A_^6<k&~__OuS
zHnAv~p@qX_i`75CS(S)$NmWUmE+GD<6^KF1N55VfAQ*qSg3427NIgxIUmpV$K|UyC
z8h1rV)G~XksgQ;6K$7@_I4fM`@s}#t?O@^7y}e-3W09`qBTO&sW}msRwq&inY(JP>
zh}tfDVxh0t$pdi7E;^mQkPqQDiRJ#s<iaBF7YFfpe3q~lMq$1wuA;B8^|wZ}kVYF!
zgh(2^#e$Qyz~`@CWg;CM7-h}zvI35=#&@;PNu-M0B&;fPB{$HHN}0kwk|rjZSy(XO
zSC_2-Ip?U0^-droJ)9u;m@3Jpr>7701|W_M-<GxyK15Wf#|_MbqNSv2<*0lbJCc~1
znTdrm!@~iJbqIKGkhrb<I|!tBTsF)F{M8)8s1SXq7D_#K1DPv<plwCx{|xe&Ecl-(
z68SVKK9q0@^O4X6a^Kyqgpk%#m~+O)dOLu`NRx9V#Hu*2jmJN&r{LeV760pQzki2q
z<c9Fm<!3$aUcW2p6&13iFIw>H;Y_As0}tx>len@7^1ZY_E$bP{-BFwjRIPjz6ck)#
z;r<_Uk~$0V!IzMpnH**z013TC4pfY%=kwaMZzR+``O}sKtD3Rh(#VE!IFkyrqf@jc
zroveBSWDe%P4IDFoC$3SRPPJOh17ocGT~xQy7$1)><D0&7yw7k<t71j3FJDX6U!k^
z0Qk7LOgy&CFN-Mqs;`fB7kpH4wg)xIg$01T$<J*hw@OJ-m6$b(P#V1Du?=)=1W4qT
zMf9)Y@0>cMk|7U4WC0x=fHa=CwnMpl_wHe+=9vZ~GgQ`$c94^WMb`d(`z|?<@BSOg
zpDFSZguK7m?aQbz91cf=Q^7!jv8rBRW#EDU0YxH#Ox>^(!kxk|z6WRZfN*`7QdQT%
z?YKh6^fua{Al&!yTrG2y70|xr?Nj>xCwhxFhCnPVA`U?dXJ37Td#_&Fe!wy_Fjx!#
z8VHPIpK$9J4Z1B4)YaUZhip*uYX5ekt2ZQ8`tiY5&D}&$9<a3`VNrQ%bs`o3H88zF
z?qG$FC>W{R=_gdFwfqZW-HQYyhqff7`Be|AmaQ5`<7E=0%t9G0CMlwIbpoBx>Z<0i
zE8zuFS<|<OG!2)n?*|Bkb$3T#)SWx!XBsZanrGQxLtrpli$K6_Zu%G4+gGQ@9I9Gv
zv4h+*>NZyLVjXCw(ez)W@TSSww?dwbrsv+72j51*F3;<7^nN)x;VS|!;Yci+M!i4^
zqAiUjb#l3&612MLNCVL{>?u_wkFc8!3Ofpwo4h@6P;ai|v<#^7ROHz3rzk<y5o3ab
z7sbK9;-1T+dKOjyYftl%Ug}qVH`&=HLz~BaZwyW&u%~{&#1|P)#^0tV!kHR!t{t)o
zx?g@j;ipP0r2)@J_OQ24ETE{jn=koj0}GhEky2Upg|NnrY_M@yt}ZU~LFb0cKm{TU
z)dR-lJ|SVVW(k0LG0P3!|IZR^FAxbLK_xfMyYSV^>Zilr2>&njkE(_JyEptT*diP-
zOIa0zb|Kp~#Po)|e|?j{a}D?U4FK_=P}v{e0fOfN&@mk#t<X!fkOF4;i0@>_<9z6t
z5BWlvftQjRP*}ib3s;A#re=kXfq6jeyIr1m7ypbE8qxguTcI^K+|Y8Z(7H_=73(?+
zp4RC_e2{clEWlq^IR~P~OT=myLhoY0y`&Q&TCVIzdf;kl`|#g@+(zkxhkMy&4h=49
z<=#ySni;7TsM%Sw?#2t!_k?u{R|Q#Lsz{{i-=%uK*3Rx)zLsO>^|;sSKu0={EkuKr
zn1FFACID(0q1z6&=#)(>b~0Vl>oTJ_D)ed6_IcFoYwtEV?g=Q68UUH&SLO=rM9y{~
zIOi4LxiZ2qMR4!nD}fVATSMv`LYUt!JG3ZP*QaZLBF;_F?gB(^*bwQeEoYr3Z>lsK
zi-%1R81w+Oj+oBGb+2$BO=rWIh1Zl`3=osBKS8m_pc5xxdna0Z075ul;jSfHnG1pq
zBNneEoAtXJ6X;xacDB61j9b};3Ep^oy&^#I0rsF?P3K>e>QQy^_5IAoS5}KDG_VZ8
z@OWC?Dp$p(Jp9o&`_Yo?&%TK9@Xwh(v5b{$s|c})Yxq}1ML?d+*8_KxFc{En%ALzM
zGc`eB!2mfEi}wLD{>_7=xQU&`O}WKRNlPQ^XKo!D1bO{jqRrc<zyz`YiUPQKww`Wh
zK4s?$C#f4#9XGqDDJ9+K$E|y}s&AL^xvuKg<^3T12if0t?7gy~<;VMN*AlOOxwd{w
zf168C)y(=9_O)c^^}pDrX_TXew{3sbcJtDHPuJfzovZs1T|D?qX0Xx?lTQq77yjIm
z^e*DUS)Zu~B>gOAh&&yBf?0#GsFKU0hO%m;qGk;))Pw`ZYd)+yRaaLhb$PJOh1bGM
zwc6F#STI}WL6@6+@J$1&L9JhS|Ni~uT6-LhWq36DN$^skUA6x3EZQNc-x)}##Z^HC
zmet{z<`t#7LK-bY2iK|mt#})3Bqug#^m<p%mQz99qb$IoYS&f4%u-BJXP<2~e~TdT
zcL+BpdA>Q!?W^pv$LH(@@y87mRE#WZLAW012qj`c+?nk3uyO>dq(8iU2*^t4Ldx>A
zk%=ZYIC$|O#BbhKhZyj>HqCPSTWF%)-^d0+-rma9`9Vao+483<3W4{osi`S>W%~d?
zpOk~p`S!AyC)Ygzm8|DBfe7Grh*(-e^kDS$^<~tSe-#a4N`8qH#LMbz5>NjG!Qp-9
zDu``dDvWi2AG+!UzWZjlvx#xCva%{*T%3*Q8$EiTFNK1kdSC7}*xOoIL9s8AlG?s~
z1OjoTZS&MjFPdunHm-zMMD3RM*@juN5PEqaj$~id9u}9>c`&!~HxMc<%XIbg02PnG
z0)b@#h`lMZp=QV#*#ie+clJC$j}!yV$Y=hq`yIN^P%mpFnZ!RiWY)hiHZHDK%|-a*
z=MJSiBUxr-M<8SEtfhRtCrNYXQTuWpQ|@WfJPl)({l4-l(0iUIb>=IKDF8ndH0@Z4
z)$(S48eR58z+^FqkQD)!J5Y9C(+^FhCHZ*Vnelo2q17#XxJ|j|B^k)WwlZ^+h?2Mz
zHb0XZi<fGbBVNw?x&k^<Gb>+5)Zf|Go)Sht3n$~WMiy5Ov%b3VM03bd`}`hz-GyMV
zE2Z+EZxjc)8^Z3LAjQ=<&Qs|O-}-GBi%Z{&Uzb>QYuwBR6J@g#^3V`~=@BzO_}U-l
z@(xzHxofqJXqS3uUHf{HC2OwLSMP3}|6j0k$6ub^zw5e@-+aEP|2_na<fH#x<8stG
zaNUdhUNXb0e@vY8^Xdsqn!4ps$?9hkiswJfUo6q~Vr9jPf^HvtpvoM}$1U{UYqzjl
zvt-vAjBbzlb90)sfe9y%9Dh7wI)9-Yld%t?$JfiT4f8>Mozg=O{-RF*o?AFYB1?q6
z1Kz8V{Sa6IIiBc0Xt$s@VHZXTbn6N1NL9`xPqA^Uosv|$K-qu^oS?fRbNK<TisNTN
zKjkXC{?S$V;XaS1eKAmW4pe6^DG^=E0~z}vQWh7a<|26tV|U?9d14a9GVoYOJ(g0i
z3f#lh(Z4;C()|^xYOJ&y&V*gG+Zsdn)N~M3HtVq+$i?$lwDu5-^z`IbPg7kB9A;HY
zcN942`Pw51R^Sp`R?{L$MpxGL7VfWnKE3}Z;<&xjmla~+hGk>a)XZcyPNwd-?}wRn
z0cUAdd$8sFKP>Vs3pl+4oxh<Tm@QvNt~m7s9(hpBLHM6S);yj*G1)r2Aiv~@%DW&w
zxhNq)==l8`OF|efd=93n%bb(HXt#N_;81ifG0A?gca@#giQLUUS@aH0)rH#cf?%E|
z>01Ztb!3~-mf8oFgNRACA(D1~nvowJiUyXt9EGuj{KBY45QwemLX|5<`>g+XC~|W$
zy4TeP1~AIMW+A6yo5-Zcxg)i}M;?bm)=<DTol}XH4xcViaL6NXUmgSq9yrX{y%Z1`
z0qP)UlhPhNG(LB3B?J_PAjdS-)dK>1_-ic%K*`wSHB8k4_?wN|a#Ici6i;cx23G|~
z|LfPUd-d>N)Kq;dDkuGA<3IS9!T*Q9tp9(>VE=uHIJTbr=NXs1wr>-?`d>PKeqKV`
z@2;87ls4eb;e)-AHAkVD5Qtk)pQ21$lk&rtg|a$9&-ML3z9~~8&O*5N@!C32ET}2x
z9g1!OC50f)P*q94>0kNnVD;+=;5`glp2{@^)wW^tU5zRbF3zAK9Xj?BZ!HVScS$(7
zs!;dl&6{}uN#)KdMW{pVWg$$12B6UT)0VaYFA4a=h<OoRYvJ8j025GHRFSlA-@Xtu
zEkqY;f@th$KGX0pDk_TgcN#D0juZ?-xLo|=t)xz{MH+#|1oo+ZWkVpc)?U;P9K<%{
zdEmhXoeJ--cjx5hw!eXM@c7NlO;3}W%6%BvQ$>YeS;cFW{&^Kp<S;<$eh7Auk(0|b
zQIh2t4c`YMw#qDrn3M_E)C{^_=mZ$A)R_jwLh>{Twsq^)6*ST_77tVm`B=PisP8ub
z_ZS8)9k4=R`|kv}?4r-M0-IStJ4T)y!sZcK%xV7e(%%59>tSRVY6@eMCQ=6uzMjN5
zMokEhVbHKah3CNV_=Y02SNB%t;wvXUoI!DPdj9C2^%eRTW}vao3EvECP`;hmpu4h&
zw;t}~6PH@^d)&1u)L|oFy6)It`U+hdS=DuV#JPsv`R?v+yEUk*Witqj-}2LcL)}*-
zNo#jtqculKRha6`me%u9wSe=nZgkW%tI(wc4q_hv^?>sU-_bg13*6y-0smn-49BEA
zUKk2+1##^0JJ<;pNy_xxTDYFb2&3?aaEt6!i-)J8gjWPzqWaA^^lM4O9-15HQ1los
zD<wsn-NuYG2y+3PK7<FNoyZ(yi<q^d$-Y<5=_cpS&9zaRuPS{aM7c1~`755RZMSaS
z5>0?rylijP%Wit|?oxrJ4veWwCk#}fUD5iy%*6&pxqeXSub(CG#zUQF?*pzCVE3@S
z*wfRtf;7Fu`i5F>0U21|uZq)L9BY|vXlbhLH9H!Uz)Lvw4ET4XNrme|?NC!g_{sFO
ztm1~6>LkV7kh5at4&^`mL|_aQlnK0zQlA{^EmM~#aB7+nLT<wMl^$gS1OGVU-OGB+
zo!fP9KRG1FA1dhni&bD(!iA?*60~TT&XL+MkK8Qma{cV(f}p2Sy{cV6lW~~&0#A-p
zXijrQHmH#bZ_y5<9-apTiYQK7_ng{dpku}2H2iCBOi@k}*6eEzPK9vi#JK{2L{{Kn
zHz`SItUr?Q(2NX&QZ9g0mG4gz8QTa7!3UXI@J}Y_C@@!IlguCx2NHSj9v@&!@(}Hr
z1Pg^(Lv@Dr_lQYEN;`H^M$6ll{!&bMv3<%O58Z&h!wbSK_s(>Jhb(1jmLs30=Ywku
z7s_nb=QbDBtYE~K>=_|g6ex`G*WkiTQ_<>j?#PKR4@j&>RFJ|vjg)4+S;8LUbj{%l
zXHRvE1ILF@mpGS9;7Q-PbEmM7EH?>i%mJTAUp*;<f#w}>U~RQH|H*3mjy6KLV2^iq
zbaVi$B=yb);SKZ0csKzM)@t~##3Z?{0o7>7CsAJwm{;yCtNMb@nZ;8c)%~9OA^N=N
zCzX|DwL!B38){#EBRh0Vq`>*1T@VI&Al-xRUiOX&4>Wo=&Y2mi;$0BDFZ_GzD3W5`
zYK_9_01h8=4wOXyfx{JX)-T}kE+tEuRn0)vw2k2RJS2Fo056OYbZn+72sO3(uOon9
zqAlI$E`Aj~`{dAmn+oq?-<o28^JRdg$8%~x3&bARqHjBWZ0_8seGwf2#KnX*Ifj0J
zAK;2l4n<p)2iU!EX-g}jv-r}$Qt#C^6Ut)omp=njv@FDLVyF&WO`5&!w@|PT<N{Fr
z7rCRRHZ(7Fl9HS0JuJS{>;``tRZU-Ch|9dT^eQ>fWGMppZ&^H`AW^f>*k(_9BGi#H
zI(ns>#UH1&z&U-u5v&y!_$~5zpDpq=RqAp4b(arfK-Q=AWIun6ia`tG@KGkm=7dR|
zyT`}JvnKi9HbF|JxKc>V3j5P+K!<S_@Txn&xLpDYjq&&dWWy``^`=()Dfu^W2Mo}T
zKodjUaHrcT<!_&NqT;E|TlM{MQU(oCZ|yPfl2w4g^`&@p8JOcZ5UUd*+z<WP7ymxK
zYo}GEPv46V`_Q1*K)#mSBtvy@*857KvNFFP0(&=xDF@i9@aluIwU;}xLZQA1c&k?7
z!m_w)r0}rhSL%i+0G3)X%pTCOFJd2q&0*EAEtL9YPqmC^@AS(=<&aow=XgXdSXBPn
zEO`a~-*5|z$@%}P81u@D>GhK;z-@+b$nT*Q0mjt17!`({P!KT|NUFEC+s##WWeO`c
zO9}Hqj8nrca&hFv>djJCcAI2FhhtQSgdafjuo&0d(3U#K102vGwn-leIo1L6hIT^4
zV#6&g^Jd}V`vh!Kd{=FM*mm;76r732uj7=A96fr}cRbZKbLpg3@19C?(C$bDs;0k~
z(MZZyKNe>-$vc`gN0=)w#UF0kU#kBdTe;m_vj8_Z1Iu~w{3AgGNYf)J+hnw`-9ma@
z#77G_)5^l`ySd6!h{4d)7k0JBejll${JIvxEg;2MAMP>$?sn5><n1TW&r{{`_zmJX
z9&aoS+7+6_A7S1$CV=>rf`sG)X95`4Trm&WUOHOe(!#i7T+^T-70#?zB8HJvH-+{!
zdd+TlPyGD-vK`9Jj<)oJ5|MzWluU<O*SM7Qm*FH2a1?&8)2&@Gm9dE?co%{XO&YJ&
zs`XmLfRX)&OyYi_X4KqAb_oTNHv{C1av(X2KUr1WNG)s8^Ifecuq?EL_P4wEVxW!Y
z=NEzwO&T23c3)g)32VVi7%MYS_1v&|dyO_!Uj~v0nx;GeRr}QTSEJ;=herPgiYQwx
z?bc$<i!WaLYD0XF%3=X=Pnt!{uF^_i)Xy_-;B@<W#c?Z3JJWu5B|YCRn+i@oIo|%~
zuf5uPyvE6!>9QFI7B7Ogd2wVnJA!W@$orBcBxD)l=K?@2d(7sB7Np@g(1=qweLVE`
z{}-j<KX8ctGihGp>5{Af-FQ8Q@dv6jiKTi_`yCKh!v-mHl-6#8_tqfyvw{CtQ(1MV
zcwQgDfOvQt6dn4pz%+MDi3r?EK&Gt0TgyO9LBa25qPo~k7X~QBguC*Z@HZ*7eqr}@
zpqaM8nOwZ~>aNrBG7twjIk`gEsY`2XCBk@=@){)YtTs5cg6M%m+U|qM?2kL+1!T}z
z0xupP?CX6J<>m=&L()eA3)bGV+91bPHj2h4l$BB+nuSgPk@#^{yk6M6E@))0oF?T{
zKy>=t4guX-Pb~8zjzSBqJKfhZ<uYD*$m&lZ{%6JE6~#I!zxPFL52#XdRHDZoUM6D;
z?TfM)v-`gu{GcvkBj|ziHsWOgv6>S4k}I}KRO5pmISX0U7HjC`Bz0$H(E%lBVMyh@
zPJxNWRlq*+cq+kzNpR@T=v=$rA1jM~tgo+6NWk}~*g2%N*8seWHsXzes+;$rvZHIq
zKz|#hI6h6v>>J1gwF(b0BwTdy2wPp|pFWtpHkNhI0Ux@8^`|j&lH|Ozf(G28m}&qc
zFGs$5l|^Z7b3q@1vJY2N;9#i2ERebdHY}JHBPDoJ$zgHcRGbhpFioU)q5uZvy*OTu
zJ5-BahRG91<gF*>?&s=HY{R|uVa$U+QHLWIhNXS_KEIJ#lGP7iteBr2#JqTN|4}Z$
z$+pFz{ORkO>868BjOt`Gu%Qu4XeU^Z2*7x)*n<_E*l5&EJYM^FkL{{fh49(*@yU)U
zk4>3<KJXS$+rWhjd0AAdj7P4y)K&Q2HB!A)Thf+yCcAj~=}wngZk&&88E{2Ap#fcb
zv$z<Xt!!|n0d$B~afFM_2&`E}?;5tTjs*|9SV(J2*JAX=m3b+Z{I<)WxRAIs8qU-M
z+{dT9TugpCqt_naF?toOLnyI(`SPdRh330%(+`LT{gLD!Mas?pLVM)7`vK5aAX<DU
z{%c+Rkb2;4z|wuH+@`023b7SH?~;(0a0o=qKOXW-VR%!e)fYww&mbH2xf1-^0KcV8
zO~5A$#U-%4pff&)$f~Pd7!jns1p?TDSX6{T$8LfhSIeOY*n<{;0>x;04^@VMY4-7%
z)d`sF@HqemJZvp(bmbTq$s9PaLIhx$2(rV<S+TT+QX-;ha94~61JsIz)8IUUYZe-(
zV1t!@<0bSW2N4r;ERa*X5P#85N-Snl`wGf7vXz^M3*tcY2Or`$c`?xpmEg8*DCSAZ
zA!=z>C=H#-lNz?oV(SEd9(`vc%DRish6TTp)VHwE>tF%y>w5_52(Hk9I{74~hW<zr
zWI3~Uw<Uqo2D`f44veu=@nFs~HIE%DcAnSM@&lGfLH{PZKhY-Oc_ZOS1{fF0hlR6~
z@!%%T(e#g8`-Kp#;%tcr&5-h?<>_q`iA=SkL(wrdgP@aI@4U;`EYd38Hh=*24}<Q)
zLc9DP)huSWib?`iY@^Idx+E>0BndTuDjFsyuirD<pYD;HG@_f}rUV`7?i2=3-Co<p
zKSIFN6Z<_Wu5935f{V?j2Nh)yK%W3oQdnVzStw5W-I^tD(gPfMjW^L&fou>~1zp$`
zC*GjiBYl9uYy(-nJYZHvMy9%CP-AUP%_&C;wj2j_<6Q>ax5<T);xaWqv{v{899{ql
z$8^^JOnHcfYghd!>-!<A6mdNXG%yKN>`bHawLnHNWbs4e@%51R_~e3X{xluX5xd%U
zP<6NXJ|7i<w>&B?E(V3qG6?ry*OH9R?-3ty(g8%G^Z&1`MUL|iti{ajvG@%V5sm;s
zT33R=dG1U@l<QZ6<Kd56zCB#!KX;lGXRjAEyA1q(A<os=>M&Qv53?fm!BUk2QjT@c
z?UO~opPix_Jr>zZA)qgam$DfOw2VaJwQU1cfQ_14qC!-^WbtzcXea0@xpa47-`b%3
zAO0smT=sahpC1|#ftq^q<WQzXeotmPn(_-&H5&L`iZ$27PCe#IDrhsQTqB8IjS??{
zK|b*@Qd+&>!T9*`W2)A0K-UM7`}!NV#py5e#L=TKKr%l7fyLsxfT@<iqXRo~+2z>@
z(8dw&44N+u53B5Y)c*a@V=23kPY+xR>}!_|WlgHFpz;@&wenF>Cu6h3yVpP^?eeEF
zu5+>o1YP6N?R6->HsG%5^VCG_hqPxke+~m0ZPjq$4>z(%7GOFYZN>WB@a$lP!}-9U
zB|7rNQS&1~#)7_1BIcN`Zn#%Te~sFNboqw!>mqf0U6CvHU`vALG*^92NW08+_5Ee|
zz8@hlpvqF5XYPPB`3}^3do9clOB`th_eqVz$GE1;KYO`H!Ws{<cqnRWIk^zndN+e}
zFpI#3GR}VgtN{G<1xi7}StUKNE+&8#$l%lZITf3?C*=o#PdOFzoC39Z1T?YYc)RM&
zfj(P;U&N?HQ0dpzM}`}Btb3#mygt-N06T=0FO6I-m<OVUcdCF4z@{}@Y#&DF=Hy_1
zMAZ6zc?Uk(Y6KKUr7nVqlrXUcGkIEU!@_|oZv_)lZixKI)Z>KUc0mflMugrEQ5$ZU
z#eo{Ud)YvrCvWue=_nx>l-37(HcH5*rKN>Btf<#l2|6C1hV3ZdaKv;yW;E#xP<Qab
zu154kP{l6KA%#j@p@>WPB7JD)vBB_0AlM~XUXqdsa0fu^3BJ->9Mm$hTqWK;tdmW7
zNZ@4@%mTI4+CoL*s642WpM`J>`zH+d$2BoY&o@i#nY2(zkOFE(M9xc;g!aKmoue?Y
zR9aNp^3uJ@Q;*+9ue)Ks7d(T)pkZ4{Ta+1!dI883FV4f6y;)aryCLOgHdr^=KHk!5
zL|fvkw57j>g-9H6su&4n$(W-awrS#=fwSQ}apStW^|7(BgiN%45s{UH$7@lzarx{j
zW+uRp(&i}8slbAc4Jr|B0P^WG#ePYI&43vqN2je@&QKnTGzY-5N*#*DinV6L+n*%r
z&c?;Y>Vg1@Ln1#C5M9^Rn@lt|OC;x^FU!E0IzXO>jR4hE$_xb>>Zu3Bz;+oQYzdt0
z0Nw(9bBa;DJe?0Kdg^#XA3So(#2vU)K`5xUA$x<a^NP)?67MYD!f}rCekJ{$Le%2L
zJrc?nFYGeC?gc<XZcdWW1F1W|%e);RM<Xz*7P)&jpKSS#UH5}{dgJBvmW%Nzg73#q
zgVTmte6K2DodX<w4l)oJ!AsbbsTa2WqV;x(BPi5VX)o|;HHnK8QP44~m?zr#phD1F
zXr}`#?Y|*w3cT)8Ye5TI=1cXS^*>&f!i+X_&Q>sl(vS;TEFky8mgaI15e~OC*46!B
z7dBg#kdc`wG!@%|2QvjjOzflg*Du=bkvL+F%1h=c6k<;3LFd8pj|IDA7(oZtbeaot
zNEYHF4}g6JTKIN|p4K)Yk_8QLq;Eq%eiq!H9h}JsuT&rwh&$9Mhhx^?Nc<gGc61bg
zL(&nTP-(FU8SSzDM&pUjk$YwX@AsEih|{G8Q+J?~KNea`$&+bZ_xbvI%p*GzoyCd9
zGrXp%#eIVu(&6tVOq_!f2C8JjMb+YMS@AYh;?Y)v3+GA$^^Ss4M5B3|O>HROp)Anu
zvb}`_A+}ET`U5bWjGZc1!%%KEz&DZiI!$8qX}R|>`Ddt*=gFf1+YBZ9wZFt8b<obg
zYX3`COvsQqXb&r*YC$g}ruebPS?fm|W<G3hJqr{?Tt7MP+=X2*26$U;{*f&bCZ9^X
zOa`?<qmFxd$a~-nNdPP*(Rkl=ujiLaCYnLOR!yqRet%uJS|)J)Ne~Zmu<F{2_{j-&
zhzWR-4Gu}r8yv9)ScZ@SR<C&f08Mj2a&X-wB;UlZRJ2{U<BNv$Q;8RzGyNs5$D!UC
z0KU*NnBBV|uH~%)zZT*;-}UP;lW0&1DjNj*y+~m!P-WO^%OHR?$+r{Q^xnXI?xr4H
zeJU;hp`RS&DMgP_gNVl^<QjMDf(n(Pp@(=AGJr@PQV1iK_G~j-XMMJiDqLPYmhkzi
zNv%B<3e`u@!BZ9v$jIn)#|=naDK<mNzlDp6O4`(bT>A-t2_waHH9)2qW#yndSY$-s
z^4%x|^dd)ycguzf23F8~xqcqFdH$G{m?IO=C1nwS0*)b_);)3uW)%d8<c=U54OD^1
zwUgZfTCV)vkShr)<A{*@ZCbh=pN_9fIp`auUkX5RB*3Q`Km#fbpay{G$x&7Y&l9{c
zE{=-k$a4l@lV=^()(#eTfR=$7%WdC&R;){d#RTMFVRLd`8)!yiSbw>>Zn>`p5h^(s
z575{UzOtYojFhNhhI$C(-{DJO`MgjU??mwYk1^XY(9~Bd4s*T_!2Zb-J7<9r?siI5
z%R<dZf`hKUKBXeeVqFJqH@l8tk#WG_!s0E^p;5O({7eOZq^=5-qLJ{y!2|H2ogiN8
zt#xhGZRxnKE+m4r0iO6`fUHr2f@k}H=9cToa5!<U5F~j~4M+-|;+0;~>9__{|2M-9
zvI@x05)yjK1`UbOF=^=7ZHS*6(smbwr7-4;4A+5D_CZ9zhSerqc!fHwLI;uPfy!*I
zK?C*TMaj8n9q8EW<oSMx-?Jx|#aDWz<taY+Z0jFWRrTY_5MU|Xd-)6L9q(T7$Ay4X
zWuNB$01iMWm1dQm<kK4|yA#44GVHy8JE++zpA6m|`j(9E+X=zRL!@C>gut0m1a*W$
zH!+V1ySQ(DnjhLoNc`mU6EIyJ;4R^9eR>tPOJY)4{P2gjFAuwv->G_8fOA0F0x=Me
zFFVH5xvt_;c)s=*sQpe!s$+=7;JPEQk!%2m{hWbQ$$1InPD`V(!!EBE(oK&F6_1s`
zc(TVzI7@JG+Z0#u;q{N=JCa!L4&hF&_EjA5Er^3~X(^1xpR>Gla%T^RcKI|aI4~o?
zy=;CPnlMl=D`73(0um?p0UL%e964h*dKnaAMK4ku!Ffod+29_NuG-iLZaugN=Jx#0
z;zi*cd=(VB4o6SkkHjxl;TNSKTq0h3FQnNeA^+j%iV{?_JyX=1vh-cu2Y^>vXVR(8
z<xP#3cSGoB8XOv1S(sm-*u83DpG8-ZwzjX~mT9l>i#`ftxM+NPLo3hk7bti<RWP#V
zAT+c0HZQbIr{rS$8r2A}i#FsigfVuy3_iE3u6nhmOT|HFgv(O2=Eax$p#`Tw5b-Tr
z)kQa#HD6+#_VUO;BygsXOzJOCO2~a)yX-=1yxlX0XuN#;xGtcEC9cesP@A^lmGOEr
z&IkP;AwIFGGJ?!$c+0V!Gs1IG%*jvm8X3;}owFY(;%lg@sRTVuHyYnUj=k@#D2rjr
z#HK2vtz?-V8^jml_k@AlJlI4`={FQk;QIpZUyiRB;r5<*Fb{Y(Sau%AsQIz&8M0xG
zx=rQ4BF@GyTDJ~8pTB6oAHukcS8Q%w7@1nA+CkV20c~TrSo}87RrWmPhYno+E#BN0
zkCC{AN$!zG-bRTI$l5$C+S&-{K(#u&MXANdutCQ~9@NX};52G`)u^X|X&IOdvQ;xg
zfuNe--;BWH9b)j!Wdu=dN)>r=q&x&C4ap7k%@o0**op$5BT&^Cd|(OlvJ9jh54?8R
z2q$ECNO=)`t3d<8uz{!&f!5-JC-?L=w4z`B>Qa6`#B^|F0oW1t1bY{SxF%08iMiyB
zgBT0lP{}Cyrq&VWZ~%G`p7o+GaG)YTMJK57++yF|gwg`PS5TVhG2S}jen_sYE8B*7
zi8(P&;8>X)V>&B7Q!{MP;OTOHJJOhZJRdqH<HfAfH;H&4Ks&i`c)cOgZR6pqU(zrF
zy@597TF2Vq;=@pie6v&o1*HCXd~P$hBE>qky*z_8w6PFfbCe_RG~wQF@?du*?!-ie
z@jIdyZ+mkdq9x67Lg10Hpu^L)0?$K>H}`;j$M*5iY0h;NW_>U0FRv1DbSLKCd5p9z
zT-VB4dlbyG?t_r`u{#~M>0rV&pu&gvy?0ad?W9W<Fk1CBMy4mm8{Tdi(pj=+Ur!fw
zo5$i#4K9W~DUtXxKi*w^f5Sp^X#wtt%e+6Qf$!y^ySi&!@JWL+5muTmUR*a9^Yz7>
zJhb(XZcLWrPziKQ3Od%5Qgv83f)5TmqKjwUDQaw$-$e3{K+bTktT4*Ua1*=69fo^_
z+4w^zv>Eb2M}rP7_LR@-D2!3I(}Tlo!VIH$wG$ar^I*}u_-QTf%$4D_5iaF_R^c+y
zaxxIy0X3^3?4ji5ZQBs%xl1B)Jy+D%OZohZlR0uFznn5Go5|tH#yhYFd;-XS(nVPV
z14SHP3cW^cIGiFOvF99kZ0X7V9OxJf>UhYdym+xhn{UBEW4(fAO33B|HN$4yCF=VR
zi(!KUih|O*OAx<MBf+bbs#ji$a?z=X!nvFLV<oa?`J5wHGF)`Oz(f3zTLX?VJVsKz
z^&KQ6B!2~h`Xie(2=_rsKY8qjV>=i1bU=87fR&!Dm{M^J3VJRTRZqn!cD7jV9`=1v
z&rS52EmPb#D5pNs`=M$gJS38roQ@tCa7e&UZ7d8lw5d*v7!{Angm_G@h{q)Pt*N7n
z!<<xO>--)W$B8YA?FU?ENGbi<A;+{X+QrIS1&qKmHc_f>70=?{ELuh}(|EHRX*S-O
zoxOnEZm&3CI1t~GJZ?S6jV>?K3^*V4kE3G|ol>=42ohBqSh-6OZq31h$=nlFsn@gu
z)DGy@77S0donq>~==&}C?zr7;3gwgaU=Vobnq`^nY&;zUwdc{k?^t)k`79XeI)7g$
z-Q}Ek2&Bkn`Zt#!cqRo_nb~4<$*Is!fH6$&7YBxSjZx$8tvV=Aqkjt>G94I5xHNh~
zGNWRs6d^AkpIQ+6R(lUQ(5Sy=LPEj_xpvkI$cE;@W$dq{!RA`4)kg&p-1YP--{RXx
zc4ErPVEIq*s=Txf5=Yj;&bt<^_%20X7+qLC|Jy<fPUt9k=~oJ)7keh4d|??i3NFJM
zx;Ctd8eX&67go2O3u{Y&j`<o1jzNxff7lPWQD^62QtvIUr~>2jbC8|0M6<}a^YVl5
zH<eM4>i42=Zh5e0+f53dZ@Cd>>V7KeH^EOg!GEO=<AYBj8=|6vPb0(6f5(yju|U6X
zb3yM<hYqcS87OyJs&lWlb)@Pg65wjg9conhBd-V0z4i<O$5VUhC2w|__p$H6muj#3
zXo(xUFMxG5>hJGzT=XiIF2WpNSPd%Y^|&fF+v_snIU&oqh2oZ#0*||nA-$<+`XLG5
zXNxa4NYp`VYRMBfnA-NWj{Y+?YP~DY$8mG(oV1VRG}0F|7>sQr)H-X`%@~sjpJ$on
z@h6k7+rcq$>m{P1Z-R-5@wLg8o;Mo6weu4={xfQ7y+b;wSwqZ`s5VnlSb1yK#aibC
z51%~m%fU!iVBzwck9}Gs!-AiaoiAC@?6M{59zPLJ*Av~x7oM(M3UmsC)i`urpW-&P
zaHT94+N{Iq*Mwv}c1WMyb=P8kVvmQ*^RQF=9wM?7%OTOcQ?p3y!on<`uA&x=Kew^5
z>=5g1mB%hIWS9sjztAywLTFuE5^5j#BuetLmBs){H~W3`gwU>6SbBmS6_LPsIUnDG
z72Fjq9`MMeXv67wyJtktTbA8Szr-y0>M*D<4SzJh$s=APnwB;2C5>vCmGV5Shj6TY
zJif)g3(9_?OYy4^%;_!t#$O4-78Vv{TI*GW8uHUsJY4jy%)R)y;$L-UL4Xy*$hx<)
zmc9(mqCX!GpT{m7r17(+88%Jn3wuRZl`7@kTp@wIi&lZ;_h+uNxHc29A#z<e1-HwY
zfg{{DqZxk6bX4_>7N^oylFfRv^vO5wFfeT^KPGM!&%){K#rdi;?>!L|2ff@-FMh`X
z)`(5S92B85Q0aj|FAA%Nm!{{Z`EgICk6<)}>po!g!nNjdHcQJ4oU@YmwQXJPk(9||
zB917|m0dRcIx~Kr!{D-g4HjMbWpREU{G>D7l{;r(tXxY|@6yrW7t7}3>*5zI8Tr_e
z8SJbyni6WmkDMm<w8hQu&DENBDsR^ky1r7Dk)qUAkqa+`Q*xJs_z$KZ##P_)VAjLh
z3x4VBsq9!hf~&pYiM=ph^_}Z}>ViW4=KWa$blm(|Y9Y3a%j}=#C?|yK=eK^F*3li5
zuK7K(2{u03&u<Xcq~RldWnQov{XZ>#vN{3ZEYY{Y;ZKmYLux~hQ%l}!{-ONvufv>V
zqy~pVD|m4pad_uM@Njc**gkvv+UF<Ao9S63EMxvPE8*deeea-5@!uqjKIfcS$CvL3
z+mBJndKT}!vht>B<s5&iT-r`c5bRLViH9Eo9|*pzZ?=uGuuTy?{Cug@4dU&1@Jq)<
z0LAU)NhX$5JIqhdW5HF-ys}`4AE#!Wt`G;EF3C?PxmxDRsGK%YXS)k+_%YM$>FS;Y
rQdBSGh~&>-uR(nF|Jbj(<111V5q68aeysHrU&&bm6J+6US8o0fF|fzK

literal 0
HcmV?d00001

diff --git a/temp_matrix_ggplot.png b/temp_matrix_ggplot.png
new file mode 100644
index 0000000000000000000000000000000000000000..85be6a423221397da07409a3ac34f285358cffe1
GIT binary patch
literal 43708
zcmeI)eM}o=90%~@#qpvPh>NlXN|9-n5r^7EK^di30?m~8GQwbxB5_57x<Z8&3bd;(
zL*2$Qr>1GQE}3Mek+-mRA{53<2E+>4A5cQ2OzCK#&<<HjN=MI@WZ9)l{9~D0*L*Lz
zyFPcldw$R7-k!TWdA@olKQ}6JYa~GsQQSPvdj#>6**k`B@+`?&5@dR2o67g(ay)a-
z;C!=J;+aJV@`{86u`}5_28mi(wFJQ+xE%Igk#f=4^-XKxmdk&w>f<_JdYf<i@vSYP
zii?9Ugnshn>5osp7|MNqGMt_FbYAR}?9Pi_Q6D9lQ}Yr}U43r6Jh}H!SyF4EWisu<
zF4thf_eU)AruUVr4Lnj2w2?g~5P$##{78U)rT$v<^&2<MBi|Yf#RCd?Pha2danEqc
zVyeoZvnIs@jVxYC$?e&-wYBp?yIfwkc56PdVs2bBkw`oG+Ygq0(&2DS_xASYj;Ru8
zM|tfUca5!dXm;5kwqN4Uttb?VvC&4dyJF0gb|ycOC8@2ebIyBmm>JVqMH%kT?dxuR
zgHR}38R@Z;GbdO}I$n>;P%KH0r%A-p-3M(tiMXn2S!K0ac}kaCE;kyD(RQ1yNyO^V
z^k>G?1|-d;6w(HZICU+kZ&y?2MEx$N@bxoIOlC!NHGgiXvQ0G;Ew*>?=h!l(X2Qi3
z?yFp{;g8bIdc9(mw3=pzT390KZYI<Iv&Li;lS1!0W&Q<&o80d(NzTf<!_q6;693Sj
zmpjY(^lig(xvEXmFJUd}2IStAV!KDN#bF}P%EKb69VS*$rEcZVYZhkmFy#>-Dtp!w
zAthwu<}0zD=hodJi$VQAA<a4&a(5*`9A<=3lfi5nH5n?4rzXQg)2Ydr@M>z35rd4W
zhXT|%xo|n*2130@;+Bbr1fIvJfl%)WsB@u~iF!5R{=Gk8ffMt^yK(DA*=&;IyVuyu
zwC$UtrNKwCWI;JG+x521ArZqVTeGWASes4qV{Xbt#PE5`k@~#W_?%%rFD}?v3NRx0
zzfTguKZexHq#%PnGZg1it6#o4^!4$HshR>^qspR8O{wkI?o=-<3eqy#e~}oqE|be`
z>7M@VP)S!oL1CNbg8J8@Drekm&8XA5e@;1;qE#=nYqT@kPe&xqYE$tW>zTJq2Mdqs
z(|PrFSGq;}w_B~Am^x4}O-?rK(>y^r*YOhq5P$##AOHafK)@dbdd#cRJ8c2{F+Riu
z0SG)=0oCH*;f)ARl9kCnsWCtR0uX=z1Rwwb2tWV=5P$##AQ0#Ronc|Y-j4xU?vCT<
z0==*(9|Rx(0SG_<0uX=z1OivUK6Ps8PFsM$EhS0@0SG*Bfh+X*@Qp?S&Vg5D4}9;~
zAp{@*0SG_<0uX=z1Rwwb2tWV={w5G97BDux1(?aHaDNL9@jw6q5P$##AOHafKmY>1
z0@;qLJ!lKytAaEHAaLITJJ6B9ecNF}2tWV=5P$##AOHafKmY;|fWTuZfQ|$n%fp7C
zApijgKmY;|fB*ynUZAL@`3Bko1pHw}0U_WA0_aG<4-NuCfB*y_009U<00Izz00bZa
z0SG)q0dyqr5V4UJ1Rwwb2tWV=5P$##?h)7)%pE3(NCua~-h0nVEPwz6An?Bn>{uF1
x3i9+YC_P7ZlztGh(bpg*5P$##An@n~u0Jzw?{F~=#O8Qo;pXIX+Ox}K{{n)cefj_Z

literal 0
HcmV?d00001