From 46c7cd13a960fb105b7c9dc5e05fc2293a370cb7 Mon Sep 17 00:00:00 2001 From: fawda123 Date: Tue, 8 Sep 2015 15:05:23 -0500 Subject: [PATCH] new options for lekprofile for different grouping schemes of unevaluated variables in each facet --- DESCRIPTION | 4 +- R/NeuralNetTools_lek.R | 88 ++++++++++++------ R/NeuralNetTools_utils.R | 55 ++++++----- .../figure-html/unnamed-chunk-6-1.png | Bin 15616 -> 15090 bytes .../figure-html/unnamed-chunk-7-1.png | Bin 6325 -> 5922 bytes .../figure-html/unnamed-chunk-8-1.png | Bin 5800 -> 5743 bytes .../figure-html/unnamed-chunk-9-1.png | Bin 11140 -> 10900 bytes man/lekprofile.Rd | 26 +++++- man/pred_sens.Rd | 12 ++- 9 files changed, 126 insertions(+), 59 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index 5c292f6..a7cd9b5 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,8 +1,8 @@ Package: NeuralNetTools Type: Package Title: Visualization and Analysis Tools for Neural Networks -Version: 1.3.12.9000 -Date: 2015-09-07 +Version: 1.3.13.9000 +Date: 2015-09-08 Author: Marcus W. Beck [aut, cre] Maintainer: Marcus W. Beck Description: Visualization and analysis tools to aid in the interpretation of diff --git a/R/NeuralNetTools_lek.R b/R/NeuralNetTools_lek.R index 40f17df..5921d9f 100644 --- a/R/NeuralNetTools_lek.R +++ b/R/NeuralNetTools_lek.R @@ -12,17 +12,23 @@ #' @param ... arguments passed to other methods #' #' @details -#' The Lek profile method is described briefly in Lek et al. 1996 and in more detail in Gevrey et al. 2003. The profile method is fairly generic and can be extended to any statistical model in R with a predict method. However, it is one of few methods used to evaluate sensitivity in neural networks. The default method of this function attempts to find variables names from a generic model object. +#' The Lek profile method is described briefly in Lek et al. 1996 and in more detail in Gevrey et al. 2003. The profile method is fairly generic and can be extended to any statistical model in R with a predict method. However, it is one of few methods used to evaluate sensitivity in neural networks. #' -#' The profile method begins by obtaining model predictions of the response variable across the range of values for the given explanatory variable. All other explanatory variables are held constant at set values (e.g., minimum, 20th percentile, maximum). The final result is a set of response curves for one response variable across the range of values for one explanatory variable, while holding all other explanatory variables constant. This is implemented in in the function by creating a matrix of values for explanatory variables where the number of rows is the number of observations and the number of columns is the number of explanatory variables. All explanatory variables are held at their mean (or other constant value) while the variable of interest is sequenced from its minimum to maximum value across the range of observations. This matrix (or data frame) is then used to predict values of the response variable from a fitted model object. This is repeated for each explanatory variable to obtain all response curves. +#' The profile method can be used to evaluate the effect of explanatory variables by returning a plot of the predicted response across the range of values for each separate variable. The original profile method evaluated the effects of each variable while holding the remaining expalanatory variables at different quantiles (e.g., minimum, 20th percentile, maximum). This is implemented in in the function by creating a matrix of values for explanatory variables where the number of rows is the number of observations and the number of columns is the number of explanatory variables. All explanatory variables are held at their mean (or other constant value) while the variable of interest is sequenced from its minimum to maximum value across the range of observations. This matrix (or data frame) is then used to predict values of the response variable from a fitted model object. This is repeated for each explanatory variable to obtain all response curves. Values passed to \code{split_vals} must range from zero to one to define the quantiles for holding unevalauted explanatory variables. +#' +#' An alternative implementation of the profile method is to group the unevaluated explanatory variables using groupings defined by the statistical properties of the data. Covariance among predictors may present unlikely scenarios if holding all unevaluated variables at the same level. To address this issue, the function provides an option to hold unevalutaed variable at mean values defined by natural clusters in the data. \code{\link[stats]{kmeans}} clustering is used on the input \code{data.frame} of explanatory variables if the argument passed to \code{split_vals} is an integer value greater than one. The centers of the clusters are then used as constant values for the unevaluated variables. An arbitrary grouping scheme can also be passed to \code{split_vals} as a \code{data.frame} where the user can specify exact values for holding each value constant (see the examples). +#' +#' For all plots, the legend with the 'splits' label indicates the colors that correspond to each group. The groups describe the values at which unevaluated explanatory variables were held constant, either as specific quantiles, set mean values based on clustering, or in the arbitrary grouping defined by the user. #' #' Note that there is no predict method for neuralnet objects from the nn package. The lekprofile method for nn objects uses the nnet package to recreate the input model, which is then used for the sensitivity predictions. This approach only works for networks with one hidden layer. #' +#' Finally, an alternative plot of grouping means... +#' #' @export #' #' @import ggplot2 nnet #' -#' @return A \code{\link[ggplot2]{ggplot}} object for plotting if \code{val_out = FALSE}, otherwise a \code{data.frame} in long form showing the predicted responses at different values of the explanatory varibales. +#' @return A \code{\link[ggplot2]{ggplot}} object for plotting if \code{val_out = FALSE}, otherwise a two-element \code{list} is returned with a \code{data.frame} in long form showing the predicted responses at different values of the explanatory variables and the grouping scheme that was used to hold unevaluated variables constant. #' #' @references #' Lek, S., Delacoste, M., Baran, P., Dimopoulos, I., Lauga, J., Aulagnier, S. 1996. Application of neural networks to modelling nonlinear relationships in Ecology. Ecological Modelling. 90:39-52. @@ -79,6 +85,20 @@ #' lekprofile(mod) #' #' } +#' +#' ## group by clusters instead of sequencing by quantiles +#' +#' mod <- nnet(Y1 ~ X1 + X2 + X3, data = neuraldat, size = 5) +#' +#' lekprofile(mod, split_vals = 6) # six clusters +#' +#' ## enter an arbitrary grouping scheme for the split values +#' ## i.e. hold all values at 0.5 +#' split_vals <- rbind(rep(0.5, length = ncol(x))) +#' split_vals <- data.frame(split_vals) +#' names(split_vals) <- names(split_vals) +#' +#' lekprofile(mod, split_vals = split_vals, xsel = 'X3') lekprofile <- function(mod_in, ...) UseMethod('lekprofile') #' @rdname lekprofile @@ -90,50 +110,62 @@ lekprofile <- function(mod_in, ...) UseMethod('lekprofile') #' @method lekprofile default lekprofile.default <- function(mod_in, xvars, ynms, xsel = NULL, steps = 100, split_vals = seq(0, 1, by = 0.2), val_out = FALSE, ...){ - # subset xall if xsel is not empy if(is.null(xsel)) xsel <- names(xvars) # stop if only one input variable if(ncol(xvars) == 1) stop('Lek profile requires greater than one input variable') + # standard lekprofile method using quantile splits or clusters + if(inherits(split_vals, c('numeric', 'integer'))){ + + # quantile approach + if(all(split_vals <= 1)){ + + grps <- apply(xvars, 2, quantile, split_vals) + grps <- as.data.frame(rbind(grps)) + + # kmeans approach + } else { + + # sanity checks for integer, one value + if(length(split_vals) > 1) stop('split_vals must have length equal to one if an integer') + if(split_vals%%1 != 0) stop('split_vals must be an integer greater than one') + + # get means of cluster centers + grps <- kmeans(xvars, centers = split_vals)$centers + + } + + # use matrix or data.frame input for constant values + } else { + + if(ncol(split_vals) != ncol(xvars)) stop('split_vals as matrix must have ncol same as xvars') + grps <- split_vals + names(grps) <- names(xvars) + + } + #use 'pred_fun' to get pred vals of response across range of vals for an exp vars #loops over all explanatory variables of interest and all split values lek_vals <- sapply( xsel, - function(vars){ - sapply( - split_vals, - function(splits){ - pred_sens( - xvars, - mod_in, - vars, - steps, - function(val) quantile(val, probs = splits), - ynms - ) - }, - simplify = FALSE - ) - }, + function(vars) pred_sens(xvars, mod_in, vars, steps, grps, ynms), simplify = FALSE ) - + #melt lek_val list for use with ggplot lek_vals <- melt(lek_vals, id.vars = 'x_vars') - lek_vals$L2 <- factor(lek_vals$L2, labels = split_vals) + lek_vals$L2 <- factor(lek_vals$L2)#, labels = 1:nrow(grps)) names(lek_vals) <- c('Explanatory', 'resp_name', 'Response', 'Splits', 'exp_name') - + #return only values if val_out = TRUE - if(val_out) return(lek_vals) + if(val_out) return(list(lek_vals, grps)) #ggplot object p <- ggplot2::ggplot(lek_vals, aes_string(x = 'Explanatory', y = 'Response', group = 'Splits')) + - geom_line(aes_string(colour = 'Splits', linetype = 'Splits', size = 'Splits')) + - facet_grid(resp_name ~ exp_name, scales = 'free_x') + - scale_linetype_manual(values = rep('solid', length(split_vals))) + - scale_size_manual(values = rep(1, length(split_vals))) + geom_line(aes_string(colour = 'Splits')) + + facet_grid(resp_name ~ exp_name, scales = 'free_x') return(p) diff --git a/R/NeuralNetTools_utils.R b/R/NeuralNetTools_utils.R index f66e107..e22e9fa 100644 --- a/R/NeuralNetTools_utils.R +++ b/R/NeuralNetTools_utils.R @@ -270,15 +270,15 @@ neuralweights.nn <- function(mod_in, rel_rsc = NULL, ...){ #' @param mod_in any model object with a predict method #' @param var_sel chr string of explanatory variable to select #' @param step_val number of values to sequence range of selected explanatory variable -#' @param fun_in function defining the method of holding explanatory variables constant +#' @param grps matrix of values for holding explanatory values constant, one column per variable and one row per split #' @param ynms chr string of response variable names for correct labelling #' #'@details -#' Gets predicted output for a model's response variable based on matrix of explanatory variables that are restricted following Lek's profile method. The selected explanatory variable is sequenced across a range of values. All other explanatory variables are held constant at the value specified by \code{fun_in}. +#' Gets predicted output for a model's response variable based on matrix of explanatory variables that are restricted following Lek's profile method. The selected explanatory variable is sequenced across a range of values. All other explanatory variables are held constant at the values in \code{grps}. #' #' @seealso lekprofile #' -#' @return A \code{\link{data.frame}} of predictions and the sequence values of the selected explanatory variable +#' @return A \code{\link[base]{list}} of predictions where each element is a \code{\link[base]{data.frame}} with the predicted value of the response and the values of the explanatory variable defined by \code{var_sel}. Each element of the list corresponds to a group defined by the rows in \code{grps} at which the other explanatory variables were held constant. #' #' @export #' @@ -294,25 +294,38 @@ neuralweights.nn <- function(mod_in, rel_rsc = NULL, ...){ #' mod <- nnet(Y1 ~ X1 + X2 + X3, data = neuraldat, size = 5) #' #' mat_in <- neuraldat[, c('X1', 'X2', 'X3')] -#' pred_sens(mat_in, mod, 'X1', 100, function(x) quantile(x, 0.5), 'Y1') -pred_sens <- function(mat_in, mod_in, var_sel, step_val, fun_in, ynms){ - - mat_out <- matrix(nrow = step_val, ncol = ncol(mat_in), dimnames = list(c(1:step_val))) - mat_out <- data.frame(mat_out) - names(mat_out) <- names(mat_in) - - mat_cons <- mat_in[, !names(mat_in) %in% var_sel, drop = F] - mat_cons <- apply(mat_cons, 2, fun_in) - mat_cons <- sapply(1:step_val, function(x) mat_cons) - if(!'numeric' %in% class(mat_cons)) mat_cons <- t(mat_cons) - mat_out[, !names(mat_in) %in% var_sel] <- mat_cons - - mat_out[, var_sel] <- seq(min(mat_in[, var_sel]), max(mat_in[, var_sel]), length = step_val) +#' grps <- apply(mat_in, 2, quantile, seq(0, 1, by = 0.2)) +#' +#' pred_sens(mat_in, mod, 'X1', 100, grps, 'Y1') +pred_sens <- function(mat_in, mod_in, var_sel, step_val, grps, ynms){ + + # exp variable to evaluate across its range + chngs <- range(mat_in[, var_sel, drop = FALSE], na.rm = TRUE) + chngs <- data.frame(seq(chngs[1], chngs[2], length = step_val)) + names(chngs) <- var_sel + + # constant values exp variables not to evaluate + const <- grps[, !names(mat_in) %in% var_sel] + rownames(const) <- 1:nrow(const) + + # iterate across rows of const, combine with chngs, get preds + out <- apply(const, 1, function(x) { + + topred <- as.data.frame(rbind(x))[rep(1, step_val), ] + topred <- cbind(chngs, topred) + + preds <- data.frame(predict(mod_in, newdata = topred)) + names(preds) <- ynms + + x_vars <- topred[, var_sel] + preds <- data.frame(preds, x_vars) + rownames(preds) <- 1:step_val + + return(preds) - out <- data.frame(predict(mod_in, newdata = as.data.frame(mat_out))) - names(out) <- ynms - x_vars <- mat_out[, var_sel] - data.frame(out, x_vars) + }) + + return(out) } diff --git a/README_files/figure-html/unnamed-chunk-6-1.png b/README_files/figure-html/unnamed-chunk-6-1.png index 4bfe7ebc723ef92cd2edc30887e1107a7c56a1f1..dee05e0608e6daafdd12ff7c6512cccff7ffd72d 100644 GIT binary patch literal 15090 zcmch8S6oxi6K@g-kkGq;h)7520R$n0NDDPIQ4moP2q>Te(gH>d(nRT1NP?7rf`uX? zpd!5`6qO>+GP08;{7 zR?>ngMRl-B7y^6^(f9iZUF39DRyir%@B|_ zd6X3e%Nkp%F1TF%{z3IvY4un~_1IE_FQDP5q2cJLwY7CJ&2aLlVe;tcsKNJWGL3aR zI{NQTV2S5H(IL03E`%fL+1l8vwvb0LskNCF0pe5%-6MhCAi6%1?~w)^sy> zf8#eHSl`icbYXhx(a~}PYR(Dl0t8eN=$F7xzIR=IxRaX?0 zTy;~0SXEPot7^&tRjs=JKPP1kaRt6w zQ=nVmeB;C8x4tx875c@d>>5#_0e^W~k*0C-(Cpl7617@ zg`ZY}YUU)CwVey2n|taeL23Cm{In`mGnly8b}obt_S8!X2h?4NP*R7#T2pic8i%Za zJ@x-rPJ`BjY8ECE`pMIWu>1Q5T_2ydZpS6c|J#~pjD%M>%vA#p52uTdJ`UhxnY#m` zKwi|rV%PHh;HAl0B{~NVetRW`dA1-s+#5iwP=hO-N+OsXXqHF{M}{63Uwwl6r@6(P z+Rp|aevY}n!&@+-NRP#}zIW#bMs3}?nz4IOBUX36W4O?WOD9_vG-}j2Z`l8p%L=B_Uy zKuP8bOS-2e>z6B^UP>36Gqk8w+9(SUan6{k`lW6MILEpL14}TH|`dhkdk=YHezLOpVSAKt8ayqNM@`q7<%;8z78H z*%>b{{N)6?+m^T@KDeD}%*vcmc+XQQ?vpg6P0)~qIk4yy*$jb)cI2053@mP9TW>66 z8r^{-gSnK|_%XsHp63Z{z1BeEYDGlgIg-mm7vPjVlsv+&tW8oYOOgl^Vd2AwG&S&G z?Yiv06)smqk!=ui(pp&48Z^;`#Dce9YU3e+YW8m9Cd!C&ml9v$ie358ACjpC$D!09 zu{lC28a&DkXi7$R#wl{s-zA6Dico8NU)i)NhN2!=Ab2mA^8xC+T2zog;}<-E$QNIV z^EsT~5jlMidRW|AjEEAX)WL%hN(Lb6;^1kb(pMu%hdgQn-_RQw1Ut>~(tiwV!DY7O zX7n;7+LJ9x8!1XXdAUpD)b@%G3J_L1Gm3|n;yB_Bi#17aox-qDNCU)|GXvh*UO?kY zo|^ed2ro8S>C9=8s|wLH2o|eF!a4M)hRq3TWPqaVQ~Q>9_%UH+8PcDmwc$^D?t`~W zttClzczca|r^?BL;5 z#upU(hDCt_CM;;}ezAcqv3@7sey}$qgCupMQQi1tuA*rJi^UsL*bK^;CzPfF2Kv}x zt~44EvRzk|@wqG;RU{g=T9fbp3Y3+zxpkL++>HgBMay|9F5S&fy;v8dfeH~};fpk- z2emKv^_c>V%dm1TwlrUCYtu;)s&g_&meMHWJxX1bR4Rf#$w_`OZ=(piXCBzGoRY#) ze_7kJ#VSaB(CY9#Y-__AR*l|hR1&1@ZcpImdpowQnBG`JZW2^E%T(G<0pa;;Ifui= zPFVpgq6$Oa#1Bp<1TwMy{>nr@5pGO|!{+2`AFzXwuUiow4HoHXMD^u|nK2ib3tcn0 zmRBcdpiT&i^25W30&4d%A9Lp|$dUDaf;x6OnJJPL67!ZTIqqOM%#f4C(CHTbuHvW+ zmtV@_WkR)WJ4vlKr3_b$zGkiLhZH2-Ig^xwfg1}ZEn)sM8PX2dW(b#sVIOYay}43q z1Al9ECOPVZy(t-#EGk9Hf+#TSY#nv_3>DHZ3Nh}f5S3i=gnIL+m-5*Nxw@$p0z==k zl3C68q@rnxg1qPyK43M-1;Xouc+)H;v$i|anxLrs+3&)b5kXq83Rmp#)27!Jb3jDz zat|jOtLpVr7yd!k1!zoYG>jFi=j2K|{W}%limV0|BmqRisSm~cYnUA1&~qxpBy9vQ zCa%0$NK)gDNSmEVd7$9Oo9v-S!;T1tinEd?!q69?ZnL?JTSUO&V%PNHc7)Ef9;gbh zO3a^|i;%*ef(2HB$kHbah^!&}_rY(DhRp0vUq{IWCxnBNOV^hxWxK6~*5!h@LC&89 z-%-@FU+2WrKv9NZgi@~YHth3`mKu25xsIr-!!uW_DUJMNPX-l0VMtyY6pX!acfbrW z6ia`)Govt`uo903Z{LFA2DuE7j^1@yfpU>`;FzcIVQkqG8SBWsKw`Jzffsu1&ig#&)11R4=hQ})Y~?Yg`)sRhQhLRNQkMKmgn^390a zj;=>VRMgh76N8r$I63OXS>3-@g%d^?HS@Xrqgsb^o56X*WGa?^*X)-KOHN45dYv4R zQ~~kHLHz=$edVFP8!sa4o6g75J^!I%|J>XodQLNr$|Mss=VYhs6;?HzePx0SO)qhxI zN0wRM{u*vi!^$s@^zNjX6tnbr1{&7zV&2gMgd>tD zmqhhqA~=PraDO#!NhRF*c^QH}3MV%2`v{_uDk{X<-;IlnqjKNq9Ia;KOy=e!VtZQS zNu0otak4FN<8l&1_(#OrEqM;uW#fgNnxytSER@BZUV}?MYLppRue3zkwJqPfM8+>U zA=J_#=V?_tO7mIe5nB47CPhU{R{!~X8Uz85ArDcEj!!0cr})KeuPb~v);h5H+8V_!?(ipZl!f?3AcG&2?WF7tUd z;vq*Sd1o~E*;CH{npIR?A-)S?551zb>=$DZUVIFvaibLT)=JVqA@o}eP>n9pQtDr$ zy$Uha7&(tly8Walocxf(OYu&KaE5 zS<5>lE?kAXcdCh+ur29rz-@h&KHK&9LtLpK13u`JA)f>PTAggKiYC*C_Y!zI*| z72V#6gG!cs*CmC6vQFRNDcX4eAc_H?mW8y6bpWl8r2u}BTo>3U$TOm*pM<@TguLRY zW8Qk^jkybN2XL*)EM8~7GTCkBZfTOf0Q}Ur0rd+P66fuIT&$IZvX0NrrwdhXv>hK3-9b{x48 z5fIjBt5rpZxQTR~{~|CIaX$b|3?>R1hO*qaDE!#<)Y{D7Ug0C_3h%>Y{!Kth6sc20 z-<-*fvq4HXWxl)?8Jc7MJ?Z?%f+h{{8OgJM+9hr;y?p|mNqPS~t~m9$N5=<#2V}xE z!IFKJjan~bg%TEwn(pF^Nu9$uMZNF+hX*{wPpR+D{cpk$h4Xe*_yBxtY3xr}?$xn{-e-u9_SNN4k(JfzEM=8d z{qFHxnUTKS!po6)D+p>`mzmCR@IFJ{i;1`>0mN?s(p0%PL8F${Dy zD=o;T8ns0{n18TuvBSdo6M}00jZ8I8hq%*k9*ReYe7mz4V5h`*q zmm%@W-hMF^?O=90)7nA!+ppQ#6agPfJ$qK?jgA(NwG5nC@_@j53LVn4%DFCXfUw)R z$Y?RTnEO(9%=#Hr_cUVC4D`di+OF0Dd7gWYZDaU3}FWErmLaXxbfV87Mwje1jPBB42DB=?ln_ouG-FnSoMA3|7}g--t0 z1QV(Bi3GHYvIh$*LF(EHSJxBX$J2@u>})P}eR%@Ny!2E0?pv_qPgg^IQZ6uSNe7!j z_oue;rT}$`&sBEm^~YY0JYWq8lgZWz{tG&`X_g7ciKI21P>mJlSQOI2PAp}z6MC8-X*`c#)BkPWSuSli6my@KSZ$F7e>3xD zSTl^Iwg3>e2+L(H`FDDNp*k! zHT>z*MONE+-S|x)Erm6OyHy3idmyJ2?+2e(c@4nmRmbkN=c9`$VV?F#vsdS7NkE0Y zQQ<{~QL60E(@1u)kC~o-=nx*F0^Dm@r59#4nd|}de4a%)NXC|A43yu|eju@}1+KHF zJ;Al6p~0c?TUrd4cC)g|KC;^DQUb~s*@vxBaCt;C=Q2J`tmnIDGsvD!Px@%L`GMOM zK#Pb@rF_Nr8k(b6zwE~4XJp){FR<(nACVxA*FH{Fgh2{nqH59TtKBR)D5(p^5IN_g;I^7#BxcFFASv4ia-PS})^ z%6e`cgs8#tc|D2OAgs3d9STjK3=K>*tF{Xqb?4yqb*#N=G2h?DVb}g(?4R1-aOKu)N(h^ipT&pwFGRKo8;wVtTyf0&J;M@ zl+b|SB=Y`9R~ry1Q1~xj?MbaB84{FyQE?wP7$uccO#a;ijH?dY&Oo~I6>DBxdv__1 zZ-KS7@?j>H9x^5+>#kGM9x`Z3iv*}Kb+i~#%mrh$V~qrL?RQ!#CMVwAU5Q1PKaV|a zaXL|5aeEVIy?;64K8G*%G@-uYb5-+zsQ9X{NwIVv06P97vPz1 z?oD%nLy%t9wCff=_@7Yf$D;$z{k8VRk-S)Q9NY;por(g)eA-I@k5uBUCV-!!S6mu3 zHTf}yi<=!d_3BjvxhM`~_}z@y{eMT;Bfp*ZJhTtYO=IvV!1V-7lb|-E!z@fb{@vkA zZ?pUy{gq0)nqc;S^EZKSm~haQ437DS0xDc-{zM*B8j&mTr6gv5$Rn!hD%}yG%?%epx=8-&dDEBf#}-pD z2W>wmkG-p|<0gKLF>$9H?NuB-YgGe;o*LDxHuX9)Bl4%f8^w@0u2H@CImoDW?<`BU z7paayLp6{oHtPL@Fyz)R6AGKP4WOtq{W#RNJ>jCB*}RF=LTMaVS6}F(p9&tpZ2;0^ zfOZ%8ij6j6S`eMrTJZiwD=#(DX`->P-0m78IdD|=wCq9Aybj)_C)%XX;2ypKxOzYI zLuHa9)G=cJ53stDhfIJO2(G3c|4d86#s!XwxM@pnVg@3j!nwcGAHIHKsr%(gmY(Y$ z9#Z|kMsF@-0D2a`n2!#!Gw=07Z`z9Pz8sNIIvD?JSZp_WA^?%4Spq_eqA`d-dDtmy zM^vQaCgs_R5g|j&Te%L5%ZZKsU^R$&f+?cz15_`=f(}Xu$1LlLwV5KDy#>$B-jo$x ze2nOODNVf#e)F{5)!5tY&~huMKHa1AJU~?aT=9<5+eSMBHIIj0Y1%0eE;b)Ov(0Of6kFOgoV(U@i4sv{ zFgrPPC_-N38>5)74XX$_<+MD_aAniIDN^^9ZbW!>(&ggUOng%z`Q@bQJvBopI-gU^ zLv6-#GR@A+A?X9#)DuxR;Rlo5Z~KE5fn{8*4P+r?{NwV(zm_Z>y<(|SyL)^_0a^cc zKOfQ1qnZMoUpmOC*1Ej6&db%YA*x?i17`*}_J6R>Xh(=!eAR(MBPos@tixM8Q#2 z;%(u04VkXOR%0G&fGa_;AUt_Yj?0*nIH*wY@yojtB`Dd~IR)zSD0L%p1d2RY{qFgj zPnm|rfY6D6O$4`LWumfn$#Iq#3YL$8@7pJAsdoQ4+@*|L)zY*&7^S&p5r_OBYN%5K zQ7T3}`&fLUHp{^%_5FM6<;L>Ug_jF-!^~lo)hj}t4_f5WUD4=-F?RF9p4k0vfKPPY zJ=a-MzF^V6Ih((!mus8>lAq05eWQ{_5VyM3-aPqgCGYuX?5WH%o))izilaY9bZK5% zpJ}$NHESu9q0KiSYjO4sad!Tosy(#8-WY2te_L+U*auf&0|%{k?;lRLVvmRC?=L`f z<%zG@vtyiCyK|~(`4hx9=V5^loA{l7`tlu!Esd1p_}+uA;lV{BQY9=`l3jI2;Hp{l z$!&QYAMlA*Z@xJH%^xox#tO5xlgw}DQ~k#WZpaW!K5aCbr)@$5L>yOhjR&ZC6wq*T@M6cv9Rqe3ZE9Nk7l`q3bWI~}=}Odh2(*y63Hv<}0uP+HMkro9 z;r=FXv*4M%jXTw0-tw1x_RjN2WAd4jfARrm9c8(!@A&`t3%MZwL#O+S-ei?Q6|*_5 zQi$@z9eeep#A3@u$^IXyoG*!-0_ZcV zh{{s1`0lpto^b1|e_PLE3QRe#UQzP9WZtgW*;a_P3m=B@C}fr%{jlJty)e|d>EE(j z>PW-qBn3UShZMnI$C)Z}vGl@^nx~JQ@_+Eb-^vN<%QI9+{RCpFL*OwtciG>JrH?LH z@$rlu1%UG$MgYD{O7!i!5|m>&fFB51tG8AlNTZj@r*E0A{a?Pm?!H&2yLWS7gv| zAJoJn>YtSxzYEQImNI{DL7e>yK0khHV-tOKZ8vPcLu06%sGT@Cz5Hd)>2hw>Mx4-# zEbXiYSx=Sq$V5fh+SR+nw|2H^YpFwz7T366u+!GOPUN(QD3Eu<^8{7ZJF+Xs&&aWw>g^Im%DLhjk4xU-abl7r1uz zA9!F;!e3@c(~s`GN|HhPrWE}~&;nWL8amb-G48m^Hg4w5RdwzWBs*tr&nj5t$>p|E z-fMS``*XPx9ZIe6D?4%DZq* zGR-ZBIt?=|yGF5nCSjKXeX79R#iOAPG)FZe$gILQf{hqt@;5U0)c)khLw#Xs#SD+~ zvD2I4M^8Wg8@UfRJ!Crz8P*Ic(X39{RktbFFYEO2&&LK2J&TEq5=bxc59<*u>12bq z<=bISX_1SElt`r4@cWwA#QNtMGDmUNcfLLT@+v5}|2>IpM%FEIT4;v zuzswSTBj)ZO1-3w8|Rx!eIc+a57^FpA$o8b2LMR?{Vx_^?`LPYbKb9?{s;o#V4AQ} zlAETrxTreaXanShgnk$eP zF}ChDzcykI0PlC)?`W|dy;3x0xn)OX-Yrk&&&-|S*$)du-#WK4-KhqQXU@g8o)<;` z;FoMkga23KsNi=#+g|auwsE;Zu1dCm3#ZcY9IuTIVb^mV$}0+EOIOq3^o%Ntq-&MA0O@Y?y;0EAW%QfBB0 z#OGMLH|O~hcIrE7Ukb3krhpR{o7oZUQX)2XT$pl*AMSX zmjNOy%&k3Brt%j6_L%2=O=Q%sx8^8G5=RKQ1k;>`J7-!SSKMcB#2U_R_S@Gw7wl4> zy6JKfk2bd^a>Z833Bl#GW7GlxVsOY?Zxf7FzEWua6CMPWBdi8{>-pO!=f3{BQWhc> zadT<=>(u%nT@m(5(lA4l{%A(0oQF4z1M#kKFA6{$uB=I;-hz`&k4%M^mQH6n*khGt zi{L+)iZ?q_(eS{8Dut`4*WaVfEO@v7cg{9WP)t7vORkdCFfGmE?uu! zG&6V=eh2xgdHrH523zxmtt5m^vw48;{rTGgnLfh4-)8+=?82|Tx7-2KdPU-@q&afv z%!0C`E`X{R>r5#$4NIj)YNoK&wJ$H>gV+>#PV?t4FnTrkYlqF}JNF7= zR?|2Ud&zp2b9rsrDE=ztiz*ztQK)}P=V3@M5H06w9H8>I1bpM4%@IGlOC!ceAfQCF z9_a4&>)|FU;-m&e#c63l)1u|sUTQJxV&1>WZ?+AFf*WP`3V&-b)7C=(}r39b7 zDFTeC0dp7^Z;dNo}pANgD{o+%% zB@HbNNF7e1rEoS;$(vrdeNdB9-VAV5r|ap%>aF+B2LYIu?Yy^h+E8_UmfEjMFEng9 z@e`MD<>1DzvGflQA%r%#K@m6odJ%~FUP!QX1;A%is`aHd8}M7 zy#uS9^N=9E0t&-l{KlFHE{mpD@-k8;{;{X)#6vq~6*>9_@A>do}gMj%$@G3VbuzuX_q4cyzB{QKM| zVJi}fU@K_l4fve35Fx5r>_7cJTRs#*qLY5#c@K^$6u&60pBa+`O-)v38nF>SKiH{n z1piT`{YFXc^u(coq_+M~A;n3>Xb=@**X(t$<{Ge7Zui`N*9by6zsRYit23)5PwGHAB-ikpwGBki74q)Pui`BXvdt>%(Q zvE$#2Gho`mlRUN;yQhjG?;~fNfz_ho3px$wy*ZWpJu=W2ziMcfupe*(imtS1g+Dqm zBbU7*|9d@wqFz3k{Wuv2H`MVX0F>awOxr)^%6|U4{b6(u* zE1MhYWEKe2yr|BVvfn~?akNWLJHS@G^P zJlbpP0kdvKXYeABRI^kS0qKbVga`7pN>jf}Q1={CU%!f>M;;?;k6cGgHv{vQg3|R% z%#q9SW<7G4l=t0n$y_m_gt!8nHP#P%f}LFVR|iLl(CiAyI$vjtfrZvrHU9`mdC+z8 zcDi+D<5T(2j7wB0LJXTE^YNfyziTpe;o@HmN@~QU$oPCDN*fVArcsg_;Y^34b?!iA(FtcX6Kxv_$z?|dTg|{ofQ6b;lqDcM*H&&0fY^9?%m{UQm>YS zY0T4ygMwvh<>^|zPr2>+8T^Tqn)2Fz{jM1o_(SDXu(Sk94N&(JRIyUjRIg=qpW;%U=K8Zr4p@G*kwpmt>ehg*MoBE*Dn+gS-Tr1{SdNjREktQWEry6AW@G7;BdQqmI%4*sqGw;HGc)qBvZXSj z?y_>EN$YFDJy1Qh=ib|}zGb)=Y*XAM@b12-Ak->!UtXd)d!38?piPjE^I8gw&B_d; z9sDY*^ktdULt`OoPSw=bQmq#}ylU=CEPX_McBF#HPQ6nT`F31uPR9&J@z`V>NXx38 ztaJK0?<5+H`cRcR*u-2JdW$8UKw7q#dv=W;x-8l zn7Cy*(U`t7{2(XKwUiIYj9{_u+Z${)zP-J3v7r-McPMkXLZ=?zqdx z2qg2MBIJRJ*3O4FdYqwt!cF)=p;1g@wF7(y_@COJC*5l-1EdDSZsY8Y^W`(gO@r7+ z74du<_`9wx1|@<`j|ZK2Lpk8#K#_IvqnB$(%giErM)^h~J&q2bjb&a6NbGb5A&pBm z@Wb8lXIr3DbJq`QbCG=To&WhKbo<{ubuhn{J7uPjYNZS)_QbiVpsG2O{Kn~BHoF2?h^U8;h+U7j^u8>i{~5SV^`pZ#u9%bA z#h8qW8B4_;KcNm-1&w3bl%Lg8=}@aW6G5{O*{MJ_`@q)w27*z zZ*(JRf+c9;HO8p1n|AfM0->p^BJ}nf<7Cr zZX%vcPNlA!+cQVj1-$)Vpe-#G9f#FQ5Y3*w%A7b4WP2}*N}Te$0T zy=P}!@+E94!hV|$`CRBP?#DOH3!=r-A=DKGlh(ND)fkeHk;A|;6VMT|@c!mj+zWwT z1DmpgSwPv}Rok2k;>f`C?g2wxf)YFIQpe}ZnR2S9vS=ri%}OCuhu@1LSR$`bzQdd) zRuN9&QhpK=^IC{6Jcr9iK8$yv4|egj%@BFzG#qR0Tnov7^f0}%0!GXT$dQ` z*ed{Z?1{UxpzqfR+P#fpwCZ!pS(3!eZ%vza5mz2lP6MokdOcchr%Ard6-eGcnhm_3 zI1%ax*8pNzXMl29h%D8vIKgUNq&`}z*lnI6iOUGl8;_DWF1M8VTEUj%oe5DbLOau& z(!XpD@Oex_8wH-R7)Vme;tC2z6!tHRGMjmb^-0BH$Ug5skdmS863###j8yTDxj1J= zq%=YTsZ%=aA@u5*Wr-C{{mQT}D~8!Ku@-yafBv|JX0kwPpH;ateKSWncQkR?oY^cy zyzE|u*E?lB)oe_xj<}LVu>tm~ec}l=Roo;xyncLLE+CvORE`lD3x5U1bR@pTIXT;( zwU8}&w3VUPJ#@iHqks7;zy~6omuW;gUXn*;r`p+)w;-Xqh;ep1^RE`lUCM~jTV7rI z^O64WSDea(5(T{WV-OOI)3X!v?_^A1uLw{uEPXQ{D6258W)!*(FzyhYodI!|o|Wv~PH}h2^-)pBrt&36BR1h-ox>@SE z3-dZtw)bq{8rm#e-(m6Fx|vzVl5Wtl6sO)^+gSmMQgr27{HE9k0ap*@^J)>f>CK{q z@0@i=#w zphJ0$5+#w!1mqN$YDmhtysIJ$17}UVV2tSBkHr2H&@&2pjAg|wt}8H``H4cvN}sU* zpfeaP-F0YK=uI)1BCsOlYYa&XZ=2@Mtq&=H5z8%4=AU$Jym>Vf?h60KuDsUoEtl`* zqLAf}_Z8cYc@M6phF^jQgOv;2%`pB0!`Bgar0QZ!3MRKAbrAUrieCx_@Y<UfXfwL z7;1tj!!O;OJ5&l#yL->~1+>+Ax}MoAN3;!hw#Ix;T=HwjHEfD|#s8Y)*v`ReR}T6< zWSP@oE%`~X~T6dpMWv5bcfV$;wR0kjLM}*}H#IloD&2;lTbpM&&@#xlbsIrQp<^idsB= zNBSXpHWf3_d(!f>YJS%8fS+~(`%@CF{pZoY%d%OFj2vs`c9T9)4KMVnN0w`fL()`P zf0nU{(-W|6>4r2juGt`bScA$tG&kw|J)`5}Lg+jeXJkld;0l^cSB5&HAaKweb=YoB4j zWHB2j=F2S~iHL9=A{68DZOTPaq{TV(Brq`FLF3EjJr-|tD&9SVx;5;*In=8p2NtQy zr+y~L1AKDAhJi^X;pgF^=rdnp7PnM{3INvgJtwQMZ&@+p6AeBZAhRiPB2V&P6z{uEt`0%KOzaVJ(4*X zSu4*^KGx*HsKjj{r|gizh6XlIk5gF(YHFCD7S5|&qB`NhtQjyv67uJCIs`FNEg;id zfresGf%vG*KDL0X(i+24Ocz97z$q9fY+d}owYz0>vBa8|8kE@X!yEIcC)oDKi? zm&1j1!W&No5$CTIO|!C+^n<#pS;}@b^=n09AXG=f(O84@JC;VQ%LKO={kP*E~X`qQxNm5V2p2Ctd)65AmiC_NV^FFuLTLfkt&my%uA*ABD zszHXR{w3u_VH;pzwqKfj7JD7*IQkSTlnIRd4ywYIJ0Zj|CM*E3Fe?N^>VsW)F^b^O z^ovM0)^x!Nm(ER+3*{yDD#?OFNwze|6coLf$%-Gk_}Oye5$iDv)O1;%WWzF_zY6iR znP8IuiwS?bh(GEWWpF6B#8sp;5msza z((88#N>c=}j_Rwh7fNwK&$H0Q4rZ(vsZ3B)`eQo&2`*=s9j{7YrPLJmG)}beQooc$ zvw}n?dj$7bYh%A3>lcC5*1##W4yki=ce}7hl@O*v^awcR0wr&$u>PV_(FDc8dSv}@ zbw37KD~koe=NbK>vkw#ITQUL>NhoBy6d$bl2j= zg)6jU#gFd538cAB(-bOXtWZm$0}t0weplINrdX@((VilOanDN&7!2T z8j}AXa`g|lJyt{~@ literal 15616 zcmdse=UY=v(04))NN>`6m7;VhA%rG9bm@XpMHB(01c(?!y3&h~gc_NOMU&H_aF@S;?z)|duf>>x5 z@eWmhhw8g2B8@IUK{sdB^$x(*7U0^j>OrJ@_s-+XHQW{e|5F72v-(aI@NVilEmR@x z&uUeUV^zbesy{_lf4ZvvtbXwTe7XMe<@&m{wRP&9>eTg@sq5?OFCN!Z?}(S{>;GB| zt2h1=9rWy;T>t<;*YV#E*z*+o2ms&*80ctO1%KGdqsd+Q;$N~Lo>UpU@<_0%CI8RH zQG%(`h+Ylh`g&T)z_A=5>3HuQS?~D~c8!~fjOKe7H9HA%NYK2;!j_AlsaO~h*O&p2 z+_WdbVjqLh1=S+|>m3>l`0<$PFAEQ5#A9(JN6~@xiZSbc25=_yN0=?wi1Ufj;%h)6 z^&Z_yOi+g~m=adL|3(&EsU}L0Kn%GnN1Op?7 z9LZ)(JZ`JnfCZv2qk(P`B`_jNXr3-F1~filLIJ$0_bFt+6&IrB0G@`U54*WKz`*69 zfX1gxD6(}=#VzGX$F}gZOP4pY8e)R6nwO`9jf;iUQ*lV-1&S|n&IZ~RalV>_XcWFW zTIxY5AMdOizyWAk$50o|Uw$i25`gtjU7k9^gn^fuKe6XJ&QF%~G#K%G9^kgfHhHvW z*52M1EqPxi6yps+%>1gq{Ai!^Q?o694w2MT zI=yT%d=hKgEEnl|3hQe*LlFY)!P(GjCe@PBAnKWcb2F-S$6vDTr9c6Zl+%)kU?2#6 zDd~%@B3}cJ-a;qg#}Z^tRBZg2+Tbz&N&%C#z8e_*JUV>Cen~@W5vx?EXhQ3|zg^bP zg;1T!a$J%C3$aABgB7nXsS)C&QdRj(+1n|X_*^E$8*vHOfD0gV=?hsWV+LB~uG zHTFFRU)es!oc>G_{ZK}9pD;E0a)8~jH#SN}LEkH5PytoSz$s=?XJMzq$y>-5@|4vv(WfNqpPsh}Q+knOfQ7~z=VbXK%k2olIu#gr<28WK+11Me(091U*m z)+xM*(g46?H2@LcgIdz33xW?TM#77BEMTI*K$rlef-3QBL}9&EA(SxhltEmo-eQPSx}N8wNrwWkwbso%z0xUt^L)B$c6) z%$Uk^4kVCroL#>I9sy~c#I>=`q)PqTz_cOp7x{$@u-Rx%DAz%Vi~`B}JA0sFkox`y z9ZheoIw$hBP6=c;XXnq_&PM!x_P>KtY1DsQ=_JUA>Pf$=Quc)-SP(-bdmaRzhrSjG zT*cQdbg6iQFqrO72MVr`_cX^y$;hI+ur)@aFj#V!Rgb=m6%afyJDrv%MH0mMKg5Z6 ze01ZE=Jr>EsaCjP^lgI+%NTT)MY^Ygj7vy5j2B8M9fL1>uv&5!~Th&DUt zVK#r=vp)Tz)Lur2e%e!{o{ZIn_u-#hxq7+UO8Jkql4k-`Jnh}`?to$^@eZ+AmPAjE z0x_Z-bXbE!grvGQi-%Lf2rtA;kZ3LDt&E%_gst2?g+9*&7F;QVeKF8j-PVz3Pjvh} z{LQ`M)9o0dejz{zhUT~frGNrjGyyfD|Lf946zJ-={z(2vcCYx11$t|wEn1A`C`BDD zlgR`pDS8fnp^XKo=G#E-V5o0IzY#Kqz36SFMT^Y~5heMSWczGK|9MJ{Kpi6*ie-a5 zj7Bsrh*-9lX12>_a#AmdQV)E2J7-H88c3i*FZGf5P&r-y0K!%^T=fp?@1Ti+Ay>0? zcI0aig4gW7X)^S%Uu0N}uL605ecqR5);I(~A#L!Jj-D9hPkFxCNc`-MnhQ3w$Z=-Z zGw~5%lm)%@>V#g|o^pKwNsaF41D5Zo=_%YqQzmDr5D*c6c@=)CB<$QrlHm!Ey9=m!LB#wYC${G&9@%ehM&Hi(lL9-6#VOqGuiU$4 zM1zUmoN6>e0_&eqn7zCA?FNbk{d6;iokC$aEn2lMg`LRKA=DR8Ux;33w74<;XiiiQ zaK56ptm&eZL`TI+9CqQN=M};jn5e#ad;n*<=xVF_H5tS`vX?20^cb!8XZN=F#^|G; z!%_gMRrLPud4^o!CDKt2Z)=3};54;JTuKqRq07Z4exiLC_5ZeW<_i*pkg{(&s}_)g z<4@v0uK(oIBkDXY?G0iQ^kN|P+r-8d=&bVo3n`@NXzdeQubppGde}UV($|x-y+ajl z{}GL^jQ%qxYI>7M|8wX0u?4fRcQy_K$(3Dok}^Z5k=pfA>7OYU(A=5*Pd6RP96EkmI;t35_nLbBK%93O^p9TB0;L0c=>_X-n|@n4c4+35!PT2Z)egtHd05A{|KXeMVO zWc65HmW59C>~~7rr!sFS)*5A#Cn&h_Ec-nri`343ZR!NM;paHf_|A++j8BY+gxT1% z*LuZ!A|}t!eHt3nXP6Ml>?-4m1bog-Tta2udEF*i@!1jL+yLYELe7M~-EsjQss~-K zlu1<8Qb1ARi>ZSGpo;6Q@`yJv{h4rlKl9{&2Cin2WLE3cB$?B*}qCdQQjQqYjiU+44OHO@fvfGp~ zJUM}tai@-J@*$ARVUz31Sc(yDU#j(3qE(DnAaOw#+VRCCGeY>ghK+ay)C`f2vy5s3 z8H)28Y%aWqIZNC(R62?-1sDSGDXw8v0)VD0Zvi zqsMo!D^x8Rw^krA%T8VaklMU6z4T^4XmM82n_*Q_$U^0c1X-HE2<+x|YSBNNno$He z6Tf>y)i?w^zG1kV`vRFWQ0fi+1*C)d-0BPw+ojFr8Gd8oL-f6&+0>_=205=IXiV;)9( zT98YCM&pOg-T*}FjyF;;z`LI46k4E;)nVW#5I8Hp!t;6+lE7Y4AGh00JQq8UjvsST z^)M3^!mpHwu}nCtGgB+=OVku`T8Bdhr~Gkqzhb&x`9TGSv)eE3$})(ap><3#L% z*oo}{$$FpRA&7 z+e6My@FqL8$J-nRBxT<9^+wO$6F$&%IoPC!xRm*2!l(sd_mjTgSB_i3l+NWzD2Q_0 zAK9xVoc;$*9r&gRgxip)+x*jJ=hIHtcIy{Urb)13$$$}XUJu;H^gbJud9y#OdvTcU zW54jJBQSaghk!My^mS4^&FqxtEdHo3IFUVPKxC0J1b_aDReaWcgRqQx0fh{XoKJ)= zPL1wo*VU_o_`f(MfY58=ePZnQwZOk#7_guR)84Q9b+I-ZkIl03*2 zG-KUIgYKy)zx39w_;g=P+3G2Byj2CiZye?`L8kfoGa z1qj75?bfE4eU&kYp-%^j7jg}#>NmUJDJ`BlUW=aI_z{MdMpq}j9M+zn;cUViLM@Q_&n~n_69uF@Tc;?l zUVp^1pjy=wqY$fpR;AVc{`{yW8n@U+bqs{1}9I;&Md%f(9`r}oqtAr)(o z0(7uWgjIW@^3OV_b=HD3sSD0Y6yDt$*1v0_lXg_ec9;<0ub0`ejblWYtfGn-<#;>N# zAAda@+huWPP)~sdb!^oM0+y1$ctV!p+|r5Yv4KC&RISX4Zov#G9|XrM9Lq z&_|p^_`1ukI3RTTd*8rjb?_BP)98TI4>1YUK z#$80_=q7}QCSmvIIPtySwfY^N0ykZOSMy1s7wSR}GGsXOaWlva7CUC*T7^>^i@kYE z<+d~-U5Tti-{HgogIHk?S#zxkkjaFqGYaVSnYp;S zQsNmfFy_;Ol`BIBMi)*oxeQRdRev%@i>lHyUJ_$x+>D3TT&~0f_#WRbQvq#-@=E|*=Nu|JW+bW>6fL9gS9NWJ1qt z8|n?2i@0b5qH|pgnRBz^z37S7T5R6wVYoCMcu!xMK$)9U_i<_O7Nhb`zJL@T5gm5gh2zn>Hb@x|4JR>pU=__SLC;y zYcsSvX;1=pRH%UKw>F1^b8ZbGNQOiYu`Z0fN58hmK|EgCdh(xr#jF5T?`jnlk2qL$ z_tTi1&=bEI&lTUb2ve0L_Hf4mFR}9-t|en{OUx`MEZFwCkT=M%YpJSRj~~O$<8ZMr zQCg2lj#AvpMTK6fQyl$rME;F`x}|q$s&K<4&=Oc@HD(d&eD7$w*D-B1wM(A+r^W0; zi$oZ1?Ky&%vY28sS1u>XO$d`rSr-Vy6FBnPlcV-^FV2om@rn4}@J)5Jmu61M+Y-|x zgZ1>6kpd#G8|i57+~SBxI#WO`4*cmZePLoc)qSL|KoMB5Ke8M_n4f#rR`kYCfumgs z%15()%Suox#`AVSEjfxzf&$+gt*RarXpr{J*l_Ajj(|cx29BgBudAH5c}3tw-~L*0 zmH?e9SEn`wY;GhFDraPQIX|X1mAq&r|{a#|IYzGv zg3!1;tB&#~H0Uep^>Ltg3$n}{R%2V6R%<$rQ{L?bzAj9#Q$YP|`Fz~O=#H~uK39QQ z!YPoAuKCtdj(1S|@~8Qnj*9^g{WBwsLApaSFw!Xi+OgLaSHTQdxxIuPIEB+h> zs-Xh+9hlHXh{`=2(i*MqXjS9cqBx1Q_?&G1xA8qy%lw!(SO9|$FOl@fn$yab78tc$^FO8$`D1m6c ze)3M7?DOi6U7YfbTOkX&fer8Q%`sIo`V+33y`w3Ed@|2`*Ng_+VkCfRs!JUl1Zgb( zZo67QU3p%ecpk51*Q(0|jyw>X zE+cK5*~4cauclB}?@5m?l_Nh~pB-GZbs{9?9pvRyq{`lz-5iT?b&H3RXsNc>3I|`L z0%kV!vW=k7t%oYwJ=14x;H-ZaE|JQS|8&c4pj|Y>nSrXBYQk@nVizUmhs|v|#CQ+L zw2yn^hN)(H#9cJMpQb8*Y;nKE6R^~Fx$*LREqZ|7^PGm{h?V5L+6u64)R9Ly?9TLJ ze&MdN&Rr2M8XeIGs){-IqVGj3^1#theQPt9+vE+lVV4}L>)phMzW%p>TMgjEzdXJS zqgq!W;QBuwz(d7)zj_p^6kdUk(Dp9m>%7W2A6G5SXX)_CbVF^E$sF|O1T!kE;mhe3 z$T#fIxU+-;;Wve%JM$2B@?>CK!faVpLX9a)w`0gLNOiq6^19+_$QVhvOKRkLp46cf z3LGT|o@k1CWdk3bkO^$@*4Ci`wQ!-YkIk<4%-ZBpk$uOD+VtgX)V%MyUgazf|C&sQ z+Pr=9ZPK-gaaQbO;LNt(SLRT8NThkB(nFI_8Z-@t5SoO@N0|A(f;u5cz{4_An0F*c z#!zKBL-L(mytd!P7N5u7nalM)-@sQ|gwNr#x8=Ul`W#?$KJ^~YYORHf-S80tzVh;| zsmiVRk7~O5(|`9Oo!m@$XdrV5UQ!#|JT@ z|LVU?*yJeA^IG&hIBjp7ul9!Q6!q1ooS4eI?d17PPO6d{cP4f7)@>2Z2_qB5xkA0lJqHOWnka4}Z_G%sn+Fr0lEBlQ;C zBX%nKbQ>~5InQ1DR~!efJ`aHu&gmxMhN+}^LIfudtVFKwK>~8>@_tv(h3Xm{o7NY&e8808Q@OjueYa+a?CazcX4 ztzC#&2t)6h>mDcHk$YqTyCo7u@~9zSZCYjMSqUjagFt3Gbt^6?8c z%NO5N2RutM$@~m~=YMx@k)+iedbu-nC58n6!>@#Cf$Qfac04y^-;>D|=W~?Hd}k&# ze3_cz1?K~ioS#L8e>l&lIC6x^`>}l83+i3oKqir3cQkVyov2wnGPp!VW;|)Y%t9D+ z@0vpxPFgp7WAv_ z9T0tEudN%61JafcZq3C%CAc@ zQAa8;OFtH+A{o&@D#)6}0Y~hBuSj>ghzw5Rw#HiY^ zcu)h+!OiF1o!^2-nI)DeE(F{KOGp+vGilL|#-lVRO)Ww-b5gb3XS`@JXBN}Gvv&sB zK*=|rW4PZl=X3+~`WE85qvBa#3dV}rrZUe2$56+`Viz?d9S^lZsM`x-uljFO=UQDE zXL-68f1CNTdM9I(HI8moW7Ym<#E-AisQR@c(2ZEnm-F0%`)=O@{k?98FO$Y3GJP#% zF#{?_i;G$@YIc!=T5jNnq7D_$?I~vMax7k%W$YI-y3Dt$OU7|$$&heWSV&e*4gjmB zp3&t*iem%M^O#4AXlAYt2XaQ08r{+x*~Y_P+VL3|9;UIXYvjObT*C`oC`y{9^|s!q z7ZwvT>MRY9FQXpmKsan2VH?sJuEJ_4-EA*NW0p1-_welQ&Y z6ViLaWh8NtXy6hC{RmJe1CcBbQLW>S+wdCcy^`h{TI3pS^#8MnI6+uX>;}}8z z^+*1`X3tpk%cJ*_aTk4(kK@1or}-sZA(|&Z&RV$EJz>#>l2SZ3lf5Rkj(mjwp;+zi zYIFK%3yG@4HCDWr_!9q8#Fd&c9=?0~MLbcwG2E&v6EK-8siphY9vVZzvdeyb(iTy) z>9R0UQ{<{{e)Ov8!@pBXV+Pavab|hX3KuSy-1)LqTRmA{b9%hwL3?Qa4SFhg`5&0_ zq%YH!8{PYl!Bs~j9C5DpjH$?VR>#$joe-X-&{kkL~Byz}_{9(V$n?f@iJrpymtCFH_9`05{oxwEzI@ob1)OJ8lVu zF~5mNMf^h?nP3!?Ih(a+X?qJ@2tMC9;6pjvR$)3c`y*hly6)`;^rB!TS>UPJ!#HsVB4_ImLX~Ze(u=n9LF+ zrkWBgY!|^!%?}oNU#bA*z|gRXm{_p&s3>FxINdKkV}DTRq<5 z8<~8%O988MQ6JZTj%Y1(9%^)f551PueA>YmPwIDo8s5}y0Br3Xrvma=14+2iC1H=o zDo7xNhqb8mDl;`8alUIBUiX5(o1NiB<8mpF*!$n3CFiclaP!(XcNABKHxP$XWlJ-8 zubh-Ettfmpy3Fh#1QqOF^bt^;n6XMJGHzL6!f z)e&PI<|@09flT?`epq&BtsYk`4{47e%zkx^z#85F?a1hdqKcGWMH#8+Jt&KY#E8BKk3d|4SVnq5+rJ?Ln*}0pSkAyFFS? z{h~WDv*u4{cwBR_#(Kj6+0TP>p6#aAYOyLnm&!5Btgk??H$|95zm9d4fHXD(5npo8 zLv{05Y&MJ~stXBYkmf7+kJ`EE;UeP_N`2}LvX_;Qkq6#Q#;gyb_6K@`raPTr_v4yq z{n#uUQb5`l{)JiveK=3bx^F)n(g`JA;PAV7-x0<_ED_+kc}Z&;(PpnRb#gdw%@L{w zV`t1NX_9RhPmB3cWzlw9qkLD)zlLtSHer5KZBv+a)92WR_ z3O>2ltsP~pTO#nNP4cFd_KiXUXX=Uk@rDxZ>D7Pa<;O!-y?FVz=N^_B zrswRXdKD}9MmE0U%l)j5vHJShx)CZP*6!QYv zzZ)S|rAzOd6(4{KM%s}*Qe=VQ`ah$s;a`lDj{DY*MVI4@Fy}LmQZj`qAiIS?K`Y9Mw%Zfkb z_L`^#Y{oP@h9eC|8~fwr#uv|yLY|1dLDJEEsE<>hM}Dm|2>uxuPfLaYYIqL@8%8IW zm$aTW_a+QYnq^an_CMCQ^lC!=d89is$e3nqo$p5@C19+={wWR5Y<14BXdL{E zE(4r+09zJUT#KzQD!ffgG(8942quP(YG|7kghr;$?c-&n+{{1kaxEugTaS2!g51euRE zROL9y_8M`i(Edpv0-fCHt0*dR;hgXkmRZ9w-(}X}zA$5Ieg)~+FY>|VX@(ZhUONM- zDkl}!z)dTM{9)4(gaMg_0kKj5KBPM2)5K>UiW-idSD@OQ=l%@D&lJZm20C}hfoS$o z!f0TMfsW>{*;Dmi{$>|-V#OoS=?mhZqR)m=NcWMpf*)8 zdMF-`k05BuUj*O9aggMY4ohfZ^8kkA(YDd2{(bnQ?)cRqr^=d$hhH;3fNX@R(9 zYV>~USxP^%q1jhTTXW?FM&0*LgH$s*VwvZ=XX^dppIn;9zHyn^5VkxTT2t$i#r9yH^{DcoXX;({H%hGxc7{PcHu|U9JItmqh`zjfBt;G%pJ| zjbwjlLB$Te7-#0up>gB=_PSvz{6nJ13xAT|FdLjKvQXhXLYipwfQvZ#=#%^7&EvXF~*?(E;rs?KEZ#Ca- zwNl6GBAJGFUjNCGFX9RsHpP(yJ|3+;TRpl7=XGVDmc8moj#%D8Eo^*$9r1>ydmNl1 z|0L%?mpQ~L3D)u+mSb)`Xmp_2m^f)FhFa7*ZGO=>W6>396VXOs&rvA*iPLldgYgG8 zG4lz@vd_P+KRTcr3@V{vcS?8|A|T@T{V{#sKd(-4_+NEbj{4)FVr2Pq7NzoQOJoUS zmuAlah0)J`QOodyTsdGP>SaVwK7mMx;wR7a>V>4?Z-zKjV^8tTjJQQpuWud_euWu=o2AHdD&j_`1k?ea z{UMkqaf2FN88^Kq_)yr|wnIn}saK)HjqO6_C}%dtFFa8OPNIFr?bIok4lPJEavS1s zOomEgm#^N!LKQ~1F{$}}g%svqUmwtu=B%cgK4w-;t-p8+cYf;HYOWOTAi$mG^RGEZ zjZvn8=AkqoNRG|qFVjmLH)`5B?_Y~y`iEs6&6<|xqeIT`oz_x5Sl_$Q*yp(lsL-yj z=Kf?`9cA88>2O8#lkyi&nO;8E0mppuPk;DMOJQ@na6knbACv0_TY#*eTX-0Gv7&#i z|7=FM8H0OnKesNxks3Cn!05~~IOxm6Nj_aUg&9B2?1NFA7|7X1bA{9w%AZ7}X`OL* z_FLz}%B4vB3WyAu^NNe^RFChqnypshMj{BcKQ0%itanr`o013C#67A>S_F6pPYpyZ z&t8|Jmw(l4|CwdPUS^hGbpfcq#78Oy<%{eFz`!w|1FYTd^^x%#M#3{w=Z9PBLgQoJj>8y?wFb#KCxiQ(!Oe8D`Lthta`Kit<(kM>AmW2~rEDSj} zmmHxEI`LE@dgpy`H=qW*?X=)!=}_|2-?@KQi5?$qAB4Pt^@z8W$me))&PwTeO;c$Q z@eh;w(;sWW6k*ZGl+=_#@}IS#f?jE|3Yg$82H7-EzI&RZ!}jl@yPPpQ8%}IvMt?l# z88MRoOdWg12|dJ&<|cTreeJDim~QtPXocgKNh!LzoHFnvf>6ABC<>?ZuC+xPyBV1U zv-XKE_!Vxg8;L`bCuCcD!)rdinSHtUSfOYTV)t{MO{`D+j>pZPDz4T1p$4@%<9qJG zvvD}BDvCG|tznEnTKQ&e$D11%D_L*zhbMM7w@u+mK1vLOLG9!^^l+Ki1paM&!P}I8I#nC`+}K0_SiT@|uNgk%>%c?nh7dfbpn_e0g(+Ix&TD67m$h z_YMj1jZ#`qSTQ(nAby#R3bc=U(;Zp*a-ic0d$(+S>H3>9v`& zW_>{tg0~?vr!!@gg>_39pSUnf^yF<+D=#kWFeKHA%&>X)u^|JeEOXcU?zb|p(n5rg zG#+M8?;n`QtdX8J8_7XMK6v?D!c$;TgFzp>@3N8hOKllrhr&zBPo`M5orV}C-&cMAGS-0)p zDgBOL|3@sW5v}9I38R!H$1aM5EmH8Go^4U-CdgFeGGLCMew~W?kCUZ^U@?_RrO9vE z4CVHJEbt|r9fCbXxdG?7A5u}LW2pWjK~HDls~l9i(8H#Ti^2MACrQXbGQIQA**&zP zv(u1Sm7FqbAm#R-a%YJrnjRU#`|D==z1$lGk(NX@kM>xT#_>n;h7SC(p?-1Cwgv+5)pWT4yth(rAb2wT{zf zmnb6H%9v<{$|JcP1@zOlv!$&W=A7iTrg7q^Yn&U=efiAyJ{w${4P(!*eHvKwLd{t- z;)?AXzG_!~hJTj9bFR@Pj_PxVi9q(kGJ~L(Aln7LmJDwa>`{K=kg-aesv}{vj%;Jj zBzVnB;Eo~)TC8n5K?(HBhZ)`lX+k4{jWPIV`8<=dS56Xb9E)FGN^8Eyg^9*lWWa{Y z%pAKH8D$;|s;{p-y*fBQm^lgKUKcuR_AVD^WiXVkxFw)Kzuq|hdG2HFSn#D4crQoF zFa+A?9@I&bo_&)EQ60u}!INB(hzpu~>hSKR3A;$Od@^5fH}@k@UKHX3W%%5-7T>^H zII~_k62Hv065|Q25`{rG zFaqc&C-ISTsF$XkRdP@v!6Bp#(0&WT_iL?cSg+zO%&q+nGYgNTfYoW`P!+R HU zHwg++in`GV!o=>yO#J?AS)neyvaV5Duo!2a1ZCaxw$q-&kdx&*{mh)EUwbnBkp;nA zBYPq0PhAi%&G-t5v*_4CB_Jnz?6Pn$|EmOQkh`Zcc~7ir#Mt0WAtlS)rDZQR3;!TT zgEfUba6E$WW37@>_&sC#a}vCl098mDX`s1AZX-7A=ECXwFBdAdjGXt>-PH21PJNz32H42t&AYP$>*)}#m33caw_A-JkXpYay z7e~yh6v*S4%C*HNC!?uRuXxj}oXtg_r8jq*vgG@bdNJlFOV(40C}KFFz=C&pQ6j7N z`DtTF3*^CJssKhBzzTk-A^90(%@p_S!vi!Aha{XS)W(OoT*zvi!Cu8q=d|eP)3Rt9 zf0|>jcTaj62tr~Hc4(-)xN0T|vkzr9+-5f`SYjCU0$KVvF+~O}qz*1IeXtP0a58d| z+J9ayWWZ6fjUQTkuu#OLwj==3?8Xdwgy3?a4^k|tYlJAV8r46N!e<&!Ex@Sz_VmWj zEsl|e4xQ1+xsi1#xx~4V4Y=C;WjCTU8{8Hy2$ORScex zfu_A=nIf@L4e2!bVxHpU{1jrb&_AuUM!OVXT@+5AeJFJ}pRvIOPKTd}ZFPgRMdW-% zhvYhKlvAstob-}S;ikD%qs&0BZ>(;=A_zE{0?$%#sUh^21TU=3wRztT8zZJ!stOf? z%tR18zCNV!>T1eBIYZ3MB2CFKc96an)lAcrLv-%$^*s$R>GN^{1>kzK?V9rS)zkVZ ze5I7P@VYamV;RU;*X?)P!^*1QaI9~NDrs=0&1lN z9-UfYWuNH3@wi4jEnzQd!6s!JYa553K^VSzU=ZO8i^8`c3`59LZcGS5yc{6KDo}$Y zUnZk+KiSO970+L@w0A>QT!mxmNfQ5lr@j->>>r(uOm8Opgg-GhE5mDJTT zcuqWK2301N6wxxQaZi@Q&^lsS$fWx+odl_7D6nayLVa1`aiBz?be5E0*LQ>G9;Mq6 zgd?>q8qZrs5}#5kC#uq6)cB#ZBMYi#;irFTm~)(!dz@t5gtLbs)Yb>>YP$zWLk?kL z6Ds)vqma*1U_mBfu&cl*T*m~HcvU1EHc7?M_;2Ox230!AJ)?ut_h~iA-;ow@L8!>J ze98?U(Rx2F7jo>QNYfPUon?^8;RYyu8fRF8#oWMHCYc#wA4NE)=u6i#*;8Ye!(R!f z=J*iN9|NrKoXJ3@KK=X|#6NE}6G7>JEpCS(5tGs4`R1n*uWpfV1w&{Nc2M8_d4ByQlmij`5 zEvMH+LFgW%Khpw0tC4y#nE5#Iw2PAwfHg?^|7yq9HECFTVxDMeob!JtnGJMJb*iBc GqW&K$$XOl$ diff --git a/README_files/figure-html/unnamed-chunk-7-1.png b/README_files/figure-html/unnamed-chunk-7-1.png index e7114c8dc43ed731ef1f7565b7f9f6e2f74d7303..521bacae308cdc8705d2026cade19d54171c6c7a 100644 GIT binary patch literal 5922 zcmeHLdpMM9*B^$NWRMPwoRTyZF(NVyyPOj$nV}pW)KEr3l4)KEhx)7>oqN;b6E#7^!HD zs_{Av3n}bgDNQTs9s8uUY|z@aGCL2*=p0PZ=BZ@vwbA1EMoVwyy<}yplbiPWs#yD}*!XX@J-zur5a!?+ zj9tix5*spfQ1%{gkjYxH*Din?2@hI^cbM*{X}BKFu- zP11Eu&l@|ul6IU((el2jO}?e=L*MDk&^h(Hu3v_({~x;ovUUaD(F?k#cP1M=5_ldQ zj07DE`X%B*n7GhfTq0LLEZ-oUWpK6-f9^h>^1zTOXL2;NV@}!F$7YIiKv9{~XF4JCy&~zF?eC_?0vd zItfN1;SxzjdQ0_a^nA6(@7Mcjc1B67;YF#rY=uY~>rwZ8X`!eDE&ElqYG z3QimDKr-ZB+t+kCCYm5T{1NcY{c>VO0&Z(ep6m6e?e@MybgAqSwS7$kcA-4D}yW#jJ+QW@U?$jFJ!QT}8KVu@1z}l)S9JT*Wp! zo)h`HDT`pU38AysPQLM+`UC%snydq<|6X-0VkTK+mD!Rx$OxQYL{g2^UjHFmDgqX){Z&*h63d|cPg zF*ZaTkx<}V6Hz9qW}7nD&kZMIsd-4a)#lg~4V7G5K;Jf4%FlYTk&dv#6^5uY2u*D{xA75C{0L>8&gcHUhE1t{K$4C>8mY+^kBG(<4Mr#^FeLjq9b z#tRQ*O#!PByVoVM=Y%i?C#B5;$*4x%#q2SboLxic2p)~dE-^ndn=j!!wRaED|5d%h zY;b5a$ZrhzI;J7hF0@D*`77@g3kPfRcwMzlQy-5HsESb6MJ0Ehkm)Wx)XHgSLPI&a z7z%OIes#L_8%o zJYPZ8SN}F&!Lem#Op;34+|Rz%)G85D`YF z#-*}?yb}1^gLbl7C>^gva?2AO-kJ7T#@l~a4RU%drxPH;xo^}?3vx?~D$YFKRqjY1 z^s;G*DvQDtjBY~0WwqRf_X|sY~V!thx+o3%hMZF zXQ=ZBQ+Vo{!0?wdCtXo$GK+bb9yc4-o|rH7UOzR&KKta7c>Q&iA}#x(_55UfpZG5* z$?wb5Flc^9Ut6AFzV61{ZCwpT7e9XqlYLQV@{7L+ z@c6&;)~HVf(WI6ul$VtwvYWhYmNX#{$5$*6N8@qnp$A>k5UZAt>V1Bu1@kpy}7QL-j34BkE-wf5~{f$y`)}4FU8V|ol9`mlf8r< z??5g*=8t71Am2xe|Ki%1YVWH8C(@U0zm};PcMaV7U+|DTA~rPHSN*NDMY>bo?`ysm zYpgYTta8^c2nZIfcL+IwH(0BGwkbk&Mk?wUBg{SPmr#VDgb5V!szmP=Mkmuo@#N^` zzMlXSKe9HkcT?@(>|pz?mzPuspB_eaj_j?LFG|rL6gwDvWjt$oMNRVFrBxT?hn7R& zhb-@w5&My=Z?AW5do9yGtCc!s5}>*aLZi59otN1>`DZU~rH=nhF*}ERcE-TS+j8)Z z^1brOPGJfKj%%Zrc0~|nBdVB1rqVgH0OT)#?7EXkY;df4`l_4vv#MIwZnxIpuI8V( zw*n-vpaz}3oO=@P#-)dUOW-_V%6NSISrd|_l#QD7v1W2rkV?J4UOT=t#%An`)|FraS}#{I?U7F{04b!;hFh?V{FYd$?7X6&aK15ND394Hoha zmmo>D?><`|N!6=mukVbr*Ufjq&AcE^BZER`+QvLsS!{82dN%;%YqDz2xV6SHlb@d~ zzqYFViWNpLDZv%hxD(MfX_NcNXFTU1v3!R%%j)ictxJbX=|!qu4XBYyE6xKo))NuP zpox3o86o{{g?6{~dya=9bVwpp;m-0uX6sJ?y;=7kz`bkPavngOZzTcFRHr{SJ*u1?n9<0w{;Tcm1cmPI;rBRExX#YWSADf+Kxnih+ z1@3yUsyB(w<;)B6+W#)vYY?+#!SyVtv|loin7C90y6&Y85n6cuDu`(|n_CUFCBtU5 z_mUu$l-s4D@*OR$a11~-(m-mRKPHkZvEYX!ct;u7^SbLkO_xcr@a-L3nOp8I-g!y} zJ0Gb%xuu0Ag66@Pu`bcCF_j-h{ZAuntPnyASiz3G#8zB-*;Vbvn|()8auv@}o$qR2 z`W8;qR9nx@caa{Q=mRTe!Q~Zj)ozWPr%(_Z41F!9)~$7T`1VcANPMWgnK#gkPVqqH zj+UTTs2Pw~GTcv3-&T@V$ZKl(tl2Jx=(E`aFe391MMTAtOU}43-kZ61)X=mvS^igdNJa+ z_O$w^&UG7lHexKW^Au@;I#jsAV&v@Ri{Bz&?k?qvE)<8{@~%Ha3N>8o{1zO?csL63 zz-$F`%l8S$r$l(z-02|C;<@U&$c{d#1y@#s1#+f8(wSS(Nhk`?oxPWaFRf&w+wZ9e zO% zhOYCaq%Bj0buR(!0mXY900a?zDEMcvmT671gPU)ufK0xMY(v5W=0CjGfV4-aXrLNp z5PhQXfJ@Y=C+svgp#0y7Z|(5SXXq47OGMDyZuUbxZD4|IJS&^+lb!r+Ii|Y8%i{u5)$Ie!|4{Veu2OlqlOkg zKk99*P>q_1A_aFK{oekz7q*@%qY}m7Nm3i2UNwL7{Zv1`S2TB{0_RUzPH|!_G@nuR zuv)b^h#}S_ZCISRLUhDvwxAko!NP5y#6ROy`dGVY~fjDjq!aPDe#hg|Kr)taSLFghv==hvz2oPN7SEp zhny4e0OTJyH~LSa&+_Qzr#}s{uC}nnUCt2>`OZC>mEwuqL+Xbd(+P9`_H0W|KMwl% zOXBypbZN8tZ5eAZ1$faBfyY|e?7m1kKJ=iR^>NJeeq2ak*{;;KD4)<(IsyK5b^h#|hC ze6eEcMSa;`qq|ZL`1Y4N-QoGyiFZyn5&Lvlk*_Pq7}rH}6@~zJ3H!m4Np$1I`x3KF z$ojz4!@xLzTr(s6&7A{@BcrN!hg`177S>TDpQtS$anLEqXZDXUQtjcaYItFGO4GVe z3b!BWWY;DMGQi2#`cF;Co&tB)fi@{LhQv<#{pc8B{gh{Sg2lquwhrmfD%+l1xB!3TZ{$KIwsU`pb literal 6325 zcmcIp2{hDg+n+Gh2q}@h(r!uCjIBcUwaiR0W>DxgmaIiggObRQcSN>QQH>Z6O^IQK zLP>hUP?IfNjOEErnD;;AVRXLt{l4#_@u3o%to zF*WNQ>NesUcH)|+B#xez(1J-q;ZoXXq;$?o>pDv7A!YPYGRM%e$6aJkxX2ou+i7@y zr;*#PlODT_FYGq)+--`HGxL(O@Y!SOyVuHZueJX^o4|dxLHq449XNIQfPLse*p-8@ zutRX{p)*(I&ql~QL@GF5Q$Su>96p4zGMZ@rc`P759Itco0-C z+)?#RQT0q!!`xN#N>}&JQ1{K!xOh*~FI&_9;n9FcM+0-Uf}Q{_0scTiZfv*Mp9H9z z8#MSS^l~0FfuSApOeeHJ=Srb&_#e91V!f+mJsep-;^nc(SI44Ck6)u4kFGdzgL)#S z(%@#5LF^mDxEjOwS|j|slL_xnCe|A#H5mW%qY2@YX>yb4ohGvsx>;(gdD>_5bcRJn zn`LIZ<-Ja;`(4)A-PR9!Z65X6Rk$0Z~AL z#iX%VfDX1eiv?`5aL-L2PRQmJzN0X^ydR?%AiS8KKU(~B@_qRDZX**(%8g=n$jox`R?%e(ps&qiXjt*0H3 z3r*$7`tgO7qGEs1pRWdkDUGUaMe5rVG32s}4T(oip?8ZPR8EW&NijJIvDhOA-gi5Z zYuj5BT+F$LR6`ZJ|GYi9X{fMt!xjF^ae>OL5|3U(J9MNPoi^W}P1t6N6ZnMS8_$Y` zSv}Ol=>}?E_1T3l6ztO%dzNkjjre8w_wNbAm9|5npFVJ(3vyc{~t zM1KEmF3COMSoTZfok!WCLFiJ=O4N9VBODFEe2jZuXJ^Z= zN=gdVuIr{3Vd=s{g4oYLmU`z;g{!1C^N(zYS1WC!R5`W`@QCyn7+5Wd%P## zGZ(+@TKq6=^PxN;%huRQYD|OYxZ^CR^}Y{JVUsxp_CIjj3FvNlN);4AaIeDX4zQtUOd2j?dA{i= zK?0Edg_gP`+oHaHbf`m1+J` z{{$4VZH@i+@52NNBEcd1jn}PyH-0_0om3?@>GYa+HE(P!ou=`?G87@Xs;0*JafVgn z{(qBjL9d>QnvjQuQC!-T68k?*;#qDzHzsZl1^pJ{Fex-VIMNiIv7bmhrnT^d{^THR+ z-Gv>1%3mQ2V276|t3xT-TGseZX3CU5E;UnM6IIK!*~9eHsLvZQqwjXfk^VX{e*9&TPZ;WTztFxGSvp=P z&BCe87T})RxaS7K^IAcAq;718A!O)Fz#?D9XNvG5B)A>Y>6iv?B(9eXxc{^_Tm)i^d3y0g^jSA4l|I1Su)& z@`dJ6*-whlWL-+QfyJC8skBuplfKz0MAE@|`6xEEAeF2eqh*Eg?IGDz4vFbP9h_xD zR=0kM1Lb2`s5`VR%L0}eagel9ErLdWg(C8)eZf=eza<@+XAS3nb^h3z(*%I(C_3!b z*{g6hX1?9qgwYAmZA;%J}{%$Vjn6oCiz~Vl-+^xMF2ig5w^O6X> zTJprep|u16sEuzH%K!SU?_U1y9$S=r>JG%hkN2T{^^6(q7xLib&Fqr|AYMSJ1dxO) zqw`%tZBMO0GKLid89wq%lsVCklOX|>H|_BmrE~K0LE2-1-*TuE;`O4dH4gY#fNpl} z&qwKL52lyPgsS@+prH|GR-HvQx-IsqH^Vi!wXggb(rW!tci1TXp2sGO2AaYug#ds@ zRcde(!}Iv9-^8N^{Odp&Mrv!0>f%JKumP^#N4c37##9@NQVz&k(R&+_}m~bly$Btw!XtYgH1)t`kLi>`G zf0K-TqtOJ(%qg&z^>(Ns<&LD6gE_*{fTnaq-j|v^@068z^+|XwFE1eb*K|(+=)9Le9vH+4av}7%)CVAVz4E3w zd(CMZwBl;^rN7j6;G2cK+b43SjcZ{o4ugA)RXJRoc5Ze}Yl)m(+%nOamUww=J)1h^ zg-n?*!wYMh$O%S6-#V|b%o4=0nq+ZASYui_ZeUOJ znl=dV1Gxh-G;}9QR(SG<5ocX?#MH;SWcA#rULS33(j4tpTKj~)WdU_F2X1|82Kh97 zLR5c3;Zwtfc?a)TQ`b6n3-6;H#xqr+D5X|v?vXf;xme`5D4_~miI8W)R zpDoUxxwJ2J(c2#x*T9*TRwJ;(YNPF`Bd}bJ08v|zpE>5iAb)9Gqpq<}nk(K}H^=cana7#((JKuy0Av&nV8V***}>5&|SRHwWJakNOr`;7m*8Qy_c3e_AY%;0D+DnuBYF+SlvgTe48!;jt3kUWoH-6Ozf z<|HGtafh!8DDe@o*5X(73-si`oxgFw8z4M$l!4s3(2fZRM#4$3k-uJG1q{F?25=<` zh>TB25wgvs`^$X0CWWg}IPVsRrTm4h*U!>|{eb(1Hb=pVp#e_%uHeE+0{yPQN%xV7 z<~)zGq9;vUIh%eDxe9(}V)gB?Zd)*!#*x+%0ROK#*E5n(f6rIS#yPob@jXuqs-PH2e(Ae5JCU$FpL3H|14Hir3$o*(G<63Q5}1 z>dro_YeQ~V0Y-sgDN$5CX~58}5KMSx`LsX#z9)ik-pch;4e0O2^X=qOru>+npB_-^ z+_7gz)7k;>F6+s4DU#3_In+F@dbdh3tV$qYDx`y4|F=>w(X;g>v1xE`N(>H3Tm+2^ z*`mRx&`iC{mXm}*W?C=-I!zeCU6aJ0w0}y4&h1p96 z?|0cr8qBw&aPwPtxrv#{D=%BUfw4L76;Wgt$pn%IbbSNeg@c9lrDl_QZ|_rK2T7HC zc&07_2&qOs3dzcWqlwF2A(p*5Vx1+y?cnJ3B-pQMX(M#MD8nX~`pSlFq`7S5QbNU- zpYg~umN2)>5(*kfTsY#o{I&%L!qeQd z3jUC2{ifCW&}Fd;K+U}pb{X)qyC6!xW*N8nOp-=wK_rd(4`{@83u9##b+O+IoI^FM z62uIPYBraCucG+mdG2-KO2Z!O8e-`>PN6T-eYV)JMjk!JD;Ie2iaITP(a@jS zwUmpC8r=M18nP9)L+Y-cLEG*(6#=4&NUx_pwYUu&NBjzdIx+0{#}QXco;=K@S=xl?QEF0;wB%4`mz5KL7e%2 z|2`?wdMXxR6%QpggG8`LC_|>e&?SIJ_4`Ixim)_bFo1wJHfW&Wc7S1`|IObi-?%~- z(S8eE)rYykFlBJiks&JkF#f<&S&FAv)#jm|_74-^r>gub@(z%I#y{4XSM(<2{Idl= zyCW=Lz=w~HZ0Z5Cr$=82V)BZlyO(+Dyg?l&WZ2`c{5Zwqz_h9wn!<7vcZgV-gsG5u z&g;HSWOXD@YbMmO=V|h;<~H=kiCLD{Wh{5e(%?-cu&GvX2$;j7oH%Fq4 zt28FhU;Wh5T(24mS6!Oc{img(Ke)GvA6++3JfuciXvwr-#_Gcl7s}vffm@r!zG+H@ zc=oJg!@eYV@DOY6l88j>INx-xv5)Sima9&(_UIy-)ku4mo@Ap{dVL_!ElhlV343qj zfo^m!(EY6Y6E(B&_Gfh2hE>VO5u{HjVZOu z@nbYvzs0m7-uFTDWn?JoPy3Uvt8j-d?wr>_Wq^0U0X2(7b5|b!NH=m+Xk9#~LU#2p nRpH)tFzWv{fc+oFw-eY86qfjNyqFjJ-w&onR)+Z}&=LOvD(dKu diff --git a/README_files/figure-html/unnamed-chunk-8-1.png b/README_files/figure-html/unnamed-chunk-8-1.png index 0062a134c00cf82c49be7521202ffdb0ba192ba4..66a8312243cf7b1dd6e51c65c492a5c63eaffff8 100644 GIT binary patch literal 5743 zcmeHLdo>UgK>lp}eMPC>oDSh+>jFN6e%+ zMP}wG5|V>4r5bWo9_1A#G{5_s>6{FmKkgrQt$Wv9*IHw(?|yv0dw=%+?*04jL_3@P z@YRUb5C{ZrX|eY(1OmlCAj=e2K_L*x0I9bFTuj+nIhcb>2t*$O!9w&eL9kWpcbG}4 z@0ZlDkkYi2(mJqVr_}~+YiS)@X{6mo-9w07_K4j_Wb_CzNbYDcOZ7 zABs@6KfTQ%a@&!p?T#_qonlqcXH}dDJB|`}xWuWto>z5CP;*aI!;sWHk~BP%HM~+a zy)S9{T+za&?ex8>?VqU~a7`yLTPOHBG9(9y%hf%8Q};yPu9Nw@LT~RrRj3zUq!&>F zI;IcN*N0&Bu~_}n6#Yv*24~6)qRI`UD-2_3sMxzG{5_+y_l*eE#^-+B6IZh*zSiV? zooNEYG?8g`vEGdIcyCgJdGb^9)TVuxn)h99*`M}&|J4^3>FpMo9hTQR4`i_pWU~)m ze|hjmk5$ea>)f|CH~Vbz-r3&zhu!V>b_E{}6%HOM8aiA&;z*mp3Q)#Y0eJpjIv!8K zr!(;s)+K$&r5*rWr{mMp@s#v*W;%tE&K#x)Bq(?ag~_C_fGRAc3ILetcqRqU90phv zJPTm40HEiReos|Z4*(1d3=A{Vhgqy)0GODV02M$XkPQqo0RU9x%K`w%1js%e&%v;Z z9=C7~g+RpX`M+h2xMFVz*yW4(&6Fq!Kl*DF{3V{sA{9OQPIx& zPR$bgnH%MwOLNPKCVPLArQ*&$z2T49OwiH1ya#d3LJ@AY^%6|!7Ypcm(c6T?6~8m~ zF8>$U;hI#9%5j`Y+?m2#jT19(_Cbfc4DH*!zRrbI5($qO4fK`k-e+BuZ{|eEpw3MN zxQI*w89Zm@jwP@kvmYLzWJ&zeh$ zftvd9@$tgmbAHK4k+LNloy#!{4UOLF;(EjV-};WotTlyF(nO_^x2>PSrss1` z9BTC4^4EM-cjDinE~@5l4gHf)62{Z0(}{r?-|T4bH?t;*;9(!vmDwm^T9R}-7?FZ# zBxSW@GYEs0DH#GYwX|Qge9^y_TZpp|9C`PanibUGK~;48nJ_K?lVdKU+SBk=dHgtv zIk!vAI+nQokHvQ_lFQKqbYSo@cR^luV&JNwWEyVIDmVz8VgMDE@t}>pOEHe0z-y6c zX!GkI5lTRm5qv)XKxi+V=hR}HYyut5{4UIOU$k>CzM+sZ%6kro)i>M98=g9G$Obwr zxEP=jWNO#_newExhrFlUwo}3!1zP0}35mBcv#Y{Kjg7R8uG~?7B|frJm#3D%uzZqx z!Ck4hPztDrX3$TDKgkS9B|VJNz1tP*X_s9X+%z>5T;}9o*8n%D`3+TX{qocsEqiUqF0MjKLIhn z;aVViy(Trv$>@~A$OagT(8!%gEeTiuCr^3e}RQvC$^ z`;0zCjFKquY^&D@4Ju}m8ADF@u^p8OA|e&rV)WE*)`BHY?&?L91igip2LtpRQKey# za%ECc^0C__!GBp^W^mr__vD6~gf6eP`h}eo>gPM8b;W|BBm18~AHAsM3omGBmqj`Z zzhFGiHL8s$7;`aQ{wy$`&}CA>o*ZkIl?fOldrnb~R`!KaFT08KfTtTTm>}v%FGb3@ ziO38Z=v9GsFC2+a>Z%BU&;D(EM3P#QIuy08+G7zrj?a$1(*HR>Qcjh$Xwyvszq0EV z`?c+_{9=NBjWOI}Coq<_Cxhy9FDM+~r`db99nm=HOn8`hoj7*R;QF$G`|XUOn+_vq zxIQUxB-aexbyJ;KK}?U8Xv{i0@d_+-b}G|`^E(^QBj!c;%$^C;(2}K`If0sn&h^-! z6j)Wo50%=*`ErAj>RUhcPoaReAj9 zcX-&i6t=C~ z6@BsaL`!QQw}flZmL0;ZiTuS?$hHl?Qq+$q<+MHUD0h7$uLKeKXa{plo?K| zu-4F_@8)@P;Km!qJu#$$*YtNgrsb1BM8CsYp28UQl*Atns9)3C;j(!XYriJwPN{rf zEtkA84%29w9IGsaq>kfsa8U|d;bo%X>xopaynTSKZYz@;tMZbdhIeF4HR{G*^ozMd z1oK!*gpApxwaP#$;0w;f*D0H`PSuz?F7bNW$fBWW@a5!K;EAz6b}9YzC0}WtVxjP>=SF_W26{?9g+aEsV)dS}6!WJ+=#xF9~pV}#@L!Ztgw->SS^^Xb0LTR;@Iq$ zz60vNfHzT1+=vtt$mQTs zMT1U9Wuwo%d8{~sE@m2^2g^@#Y^j5=;i-Z;H?Yo1L&sKF_8chN`hLcbwWjrh5FN+l zm**aipZfg4)}H{Ao(u`JtRPDI{YYRKuvSpQB>fji^o^6YUd55A@0;D8V!*PWskwJL z_KmH;!~}l*G9sbptz{2BNZrT}2={3=245^8$V#QMkVBatMN*}NPk=a97RAM$3Kw>3 z#_I_2p;!z&e8#mf0{?U=O8d$J>@t}IR4I96?yOK4Wi;E$pEGUO1>=-1FHM#?YN?59 zp|YNEi0K<#0JY}oRZ(A?27JRPq#^#NTeA>e9O9jCDN#-hmUrH&w0ueNh`f@&N9Z<9 zv7g{}8=jsS%JH^UXbKuVQ&O#=^5R|ZdO^5fb-o*oWlT7!m(TzH>BJ#Qbzs$_z{-L6hUq>Jid|9GxT{OCaG_0c?#i+Uiy|akk6mqNEW%_Z7~U6*i$XU+ z5b?z%_X}T(wEQsU!dS)qM5=i5)6M9{fHwnyzAO|qik9F5(LxAVnV(5VuHg|xnB#4s z%LxCOH1#tn>}S&bpGioNWcHPjLXHfP{>(ue#-7k>TAX&zzICz1-ybogrUmxj=6AO& zb)_X&sK!0>f54} z4|3{$>1iSde?_QnWM{o@6a(AT(P(dLs&S_XSULt+jHe4=zwRW+sYy8@HJ1z8)hBNx zEMvbWs;T`KL38lL-pd*c(?nTw)(Y09j`!mI!RJy&g=WQ4qxrnLBc(rblK+B*m;HbF zy=DFf(YP)rHWJUzf|H_V-lN^$9;fg}1M*;x`cYB^vjWE-{pr7UX-y_*G@ z^G7tL4;o(ii9^>9W-OT!eVSsVEHR$VUMNOo1t_2FF0$vO?ElPVqR+usf9q0U*;Y;p{>ucw}OoNkTIDxN7TAOQD0d3|oaY?)zbKF#2t-#E9Sc*UDz+7)RrR?3au}^eiF> zcO$xwMl}XrwJXKQSk^VCr@#AZKnT>K}p{jxH&!AN3+mpmRio3G4^(8_MEOt3FCm_)G8l;uD1bR!=Rs%`_Z%Uy!q hzc69@i}~CaeBSb{64M70_wWbDmgY8l%S|!o{s|O+`k?>- literal 5800 zcmeHLc{r5&+n+EcLP{t!X(fBv6|&2&><%ND$TCS0Vj8l=SWZWTs3z30#}HynPPC}8 z6cZ*nP9w&Ub!71To=4lz`_J$Cz3+9s?{&R%UGvBDJooqh-uLJ8x$p0iU~Os4$FqY6 z27~dLni$x^V4O%8Y~3zyP8baKG4Vqe493A{ZGO@ae1yR?VK5X-GZlua-YTQNO~F)H z@wl*(x$q%#k;9gv%2uK(HrrKgcc`7%pshj^1l=ZlC4(eO5vHtuOAk2|i$Z>A;Cl zsgq$+c9*5?BM&-WIp~CyagLI4xhm@#BkL9?=N>PQye99Fq;LwS=$WGEm8yjLUFppA zL*9QJ^0{&NY^Jhbma>1gO291@bgt^T+p2+g)PnG87w)NFEKm0YhZjd^@5u2C<( zNiX4v{S0HhAA(MQePUSy)yotY;wKbB)!x0#v9Yjx5sbx9M7Vd<@B20 z>Nn4&TI3B_-X65NJ7jhDAM1NG>-=Gx!V#PMqqYymY>UQE6n{QhGIg?a#;$DEzJiHb z0UZTHp{fB4^h?KJiZJOUOi_2LCM*>Ux(@(frDM|5F-71%J)M+ZRFqDl7p>S8VTy`K zq@r%18U?Ec08%=J1X>hfNOYjP2-6L8cLP9Qs%9SmeEj&4PD-bDchdo2YHA7`09t_G z06GZ(z@d;K0DvEXo8DFPpd2=yGjR%n!T6h@|GHN6{nId*&@oekV<#`=FsKRPt|1p& z0tE~6;k})3o~kOZ5&l02yuNqf;wGld7=pF4H@3R{hQtO7KN(B$ty%`E=@sI3X(dmY zXanZ!-O`c9!^xQ8oPnY~=JeJ1&?7D6oW<7;hj;jt-q<0yvGOE~KYr6qj;$N-H5qV7 zE2io4nD0&fcOOq~p&Xu7>W|yCLvBX?v<ory61_I=%0NNG5W4 z=6C2#=n2DHj^ktA$H?*|r5o`!@Hz>J%HV{@{h+OKtBC&tiD!4U;6P z9FQGPalFz!7hT8|cAHY{1zu=2>GItzsR8RCe{lQ1L_6hnK~pwe6VdBydc~QPA_*?VRy$g>skSdiT;iG;q^YV|)n82b$)zftfY2 zsZ{6Payavab7jQj2T9qFk5XFCA)=(936)gVHE`5VM-%2UR7+JK7Hq$$H<9_qR*Mfo z0Aws@dHk3!w&IPu`QfTJrY2`DZDLBs?<9dS#lmM@8#wwRFPFTF7>(t4c@GM7<)wSU zfPxv;Yb9O%Q+mQ1%@xj{{{9pO52ladP^``Zw*vCY7^B$>8IdVSG5okpTYbiC3{80@ zY-`JquA1Zj8E7Dt5IiO|S0J)DWKR8eF!6ufY5v~wGjvEef8y*xl2;Tr!`*T`1yW~+ zHbzWF>Yfeq@d|ELWS+4Ah;azIggGZH6*s8OW%JFFS8U)CYo@dFLH4WnFFGR}wox{( zD9SsJn1Dm~*AZD>v&(Y`xGD@HbZ_r48aX89fH(^@On;7{&BY0aa-aA{P!-v|wXSb& zXG8mPXUkiof+kCCs07Zxx5O~w$3yA{214hX&B*n)o7NjZ1!KTkv@5P>%yEij-I{Q- zLpI-m)OR;x8Q~Yqige=fQ~wq}fRo?{<+myBYX)D4<>S>QtWwn)#N!7UBXRoQoNe>F zJhWiT+(x@~ctLMk_j?*i+k`r$6y3gs5B$Z$<0B-h07bTiF^xgOW0Zq9`Y?AxM``on z`a+^}%xIW#(XCj4p~ZG>K#PW;bWh6{=hf8)wH$22_=^-YU2Ez%TR7OPt^+3gq@T0T zCG`F~fnui`H~C5Tz8J6poU7{fnOA%*z?c?MVm1Z8mJPxBOP7BE2XQ8^(Dmg`#4-nM zCuK&^?;?QV@k<8nui)QWHpc}~hL8kB|3hm*kUT9?FA%UD9$Zh68$PFWM} zPcTq*XNK}+(c!qcFVl{?$fIA%7*9ape!dNIJSK_bsCOWAeO~u7GhWrkiiS_Grp9{J zVz2}Sw4dAv@2l4PtNi6vj91U)&z3Bav#Ny`>=!h!ICt%S4bP?DmVN6L;r3iV(~<3_ zR=@r2p?_@!k@%dy7)M}f&0bua^MSE*qSt)~{(yH!N)+c${G)CJM*a;=VyM>A%{_Gk zH>ur`z=!!>H7~5nP}vo#cFBZcPFe5{?GN5y0vuNT4{peqY!&8?4uE z$`dIr{2k$CVNgn?G4B6>yrx_oGUBdQK}+-8pyHxaS6K2JSLd9k6A!W+C^u*Ma)^Wee6{yieIo+m9ZQx>RxZ{8$ensw4Hly@ zoA-&|#cI$Wl^ z2yz*V!yK7j$Vse`E*Jj}FuJv$GkfW(S;KGn*V)oRc#LDeEG5yviZW3TxcWV40WgXg z?%{QF_nO(nS*}|vI&LAr`7v(YF>o?jyj2*J-t>c#BJHs_l@cm z1mIC6;aSLA8O0N9bD&A6ni>&as*6*A0=FI59jNhh(YG}W*YJiK)9{g)^J)=ez)<$pp6R~-o%9TfAY=FFjLtSfky5N9VCev7 z!?QpWHaBG}{cq$;?4e&KRdF76tn#3aZ`P%^D1)qltZU*ccB1VCqCc#gv&`&&Vf3RgH15$Clh*}kx-Y{s3 zmH5FDs1;yJkUG&8knPGSkS%{dQGF6g`qpgoFdMU}+O@s7our(7Kbk3Pk*B!c|8ht6 z_*O+Yja>{Gp2o{ zyD~=swxS_&|2k3+RJc7;{&`$MPvhh!KUPJ@r+!q>aPOxvKKPC3udSzg5l%_d>54R7 zvY$ZbPe@Qorl|@&2y3vAl<@unQ1w2(CbHAra=eY zn?DwcS~z%`F8EP?(AhTghtjg3^DpZEhjQWncfxez_oiF|xMzsdKy6~_Oe7d8LW~`a zy17L=cQ+S;2JSHaC+XJXOoZ+zy#LcHm{oT;W~N9kBf^BPI8K!ty1o(Pf-ZCTqM)mt z^~a!#p6J`qwa+p2L2Sag6Ug0Ip-huq@cmm%unAn653v3@eSTT=;cP~$5x&rOV7q%) zM;40pH?B9NrRPGG3%+BIFJZzw`T7pZEnh3Oz6{k61s!Vt4vLe7k#83D<(Un_l+#s= z)}AaE{cL~G%rS`GPn#|5!-N#7N)2^n;X0VHbKI&D{sZIb9ubZ?!+)*sb_F6?HjO?y z?EfjYoy?ah4_F)$>5pbH3t;+>Y(`zrSx1Pk&0{0iVG vF?pux}E&G-)9qQbJ20fkaTGi&OzYLN7uHU7CXQ7DZY>r56>cp-3-+6d@=< zsvuH=ND&1=5T!ST7ryVk`|g{0bLYO7nN0GZ^V>ap&e=V?|J_ZvVWi7?hW`uz0ARhQ zr-cCkXs`ePP?(Vh001l|&VQsn*xfKR)22QG05Sl85&%F2xdLST0ZK{$C0Bq_G60hT zaBTv#UIuh_QxP%%85saZM#f(TgOTyCp%PGWN&weW_~|EE$<C)3GF)2DMbO`l9ppPVc%F0OWEtWHdvW8 zsGvemR=ZA4s9aC!PEM$h6EQ|{Bz1S0ZtGe50|0DqPrtx-zGcn;fH2^imWEmIv!A&M zS?Au}>EW<3f`N<#Y6O%EoX$C~o3}T=fjZEOFi?`#tsAtEgt7epH=sCR?SP_f@ zwG{=)4(rdMq_slWyQ2Mp8B^Nqo7{`R+)*uI1VBS#p}^w>*&t)qGqsX)e>V@d*yo1q zBerxqZ`Du&i@s?@ZoTgN_Pr|OlX<#>u|(%y;)K39Y4Gspk$Tke&|ZV@(sW+9TO&S0 zBl0LMPyM8I$g?I5)_L&OLFJ@%=U4r8D&Ox8QF{aI4Gue!Y3hQ$%og-Hf1@IKno5ic zoey`1r;qtOOy7dm4&v0ESSeBAiVHGO_Dx04vUki@=UrxPy$nv*L{AS{#e9o;#D@}e zX2`M+BX3YA5uAy(9>eW8=e`x94)LP|zolLBHBw4UUGMU`RTKC;>gZ$CiN~>b5|!#( znnl#7HUf31@^nc3pleIgD?M_IDtXg%WB|{iKWOm9?yWdS^YHVc^p!02tsI^Pxmta# zadY8$^omKhs{A64)K*`94?6xzm&9GtHESH$`!qeG0RJ9-|F!c9oJzZPq-9ej%&Bs4z{(~Zvk@Cyl{|P0J z@8+V=)H#M0e$`IFTeMDc&?y{O4@ zEA9g^xHA0pfcwjA0v{b#+@|Q(-|G5ds{skFo z87f0`5!z(4_@4l5Y~0QFL7?0=|>Ur55GFqWKX%s0cj^xUBY>;34I(FzY_Kx_QvpI&d zAvS+^4u9vCu0_q>J&O%Yv`v@$gWY<-(;%pT9j~(7=P6p?p@l@RY`!0|T$uAGy+0Xz z%tlo!lfq}v*eMoOh|)E5KJNE5!dnJI$S1}Q+})0=fVOsjgF|KC*YEdyG=$!r>e;Py zQvvk06$B^36Ais|n^c_gmL#pSD(L!z9%umw2>H%_>~~@6JCYJ59Pe;KT+pL&UWp z$`HzR1mcDQ`V~(2P3#%Sw%Jr04$M4**_8*ytymA3g92@7*OXSF80O3okS^0*uEc2I zgFDTStb)ma0W$p=$YqbWHMmXdb1#U8u^^yRghmyoZ@|fDz z!b;=)T`y03?hLkDec8!?EJ^#7YJF*gh`{fa>`aPU*#Vtsc{Fz}+L7mcT_WE7-FiL! z6X<$yhiPB+f^F#NedLZ#l91u5#2FtVYhop&h36cOIrh5Jc6l}(mZS4YGI&jFm3}su zwBB``@Yh*scft=5}Z28deeN%c6&k|RM)R8(ywdGeX}eoY*QWgvrEQX*jCjX_?IrA6zQoOFS!c z?Y=4{oxveS`^U>+vqV|!2dakwd#3-9HfOgGDY3omm>}iqy}Nz3c41xm2f-vyS}rl% zROO>l5>>WC28g7gdiF7^zxb4Q+r+2kmGju}+oDT^4@=&a>F-WPem#}O_KJd{IH(G( znNDVxyj{Cs5R1@OKwFa`&C~1Wu;UcReIA;@M4LXRI4pZEm1!$IOZG++B2u|48Xh!t zf2m~iue7Twg!H5CU__e!Nnjbf4}eIu5Q9m)cb9J#2fdJ2EI98X4o*d}kh7jt0_)OtD5SHJ)xw8ZhXNs#_^Mjsd)vDTKnpd!Q{{ZHLyL-7DDX|+)Y*zD zS#Z*O>(_1v2PmI>F%43p`g*zqPJT(9)`US2S7sJ_0GFtgj$|b-UcC2PzL)&4*q`F6 z}$$u5Fb6Hp0~*%op}hl1F=Ba zZ#5z*ecz1B=?px(iW-e$xsW5jq-7mM5!E=*1;LicEDgM41tz?3OimG|?t0cCKQISBL>2s5eil)$>($*#hiH?JK!%DQMIF5z~%(tU8=U1>O}2jXSf{B<%wCk2NT zht9F1q+X;ay%0>gi*laii9zjHDoDU3 zzQNC~veZ5#e?3c7vKT;&6v`46@Xo(~Dq?-}R9(G4l``yYDr|<;)9k(Jsi6BkdzTO% zd`V8gJi00E}Q3UA;k9~B8qOND%cl=@mIq~u1W=n_mO0_RLVa2z&k(Bzi zo-&5$z>YzeUZsJyODj@6Xkzwi(2D&Z2_;?*dD&Zu)%xXvA!WlVoXuGKO`5gs$OE4l z4>8gi4ch^W2)-H6oWi#=@0AimCRSwtTzLqF$iMc99l>x?W41kw1goFf+jF$p)pU1j z7eWDuV4Lsrl?|^LhT#+3*v;;TwKa38QrclY4u;BdSycDc?RfYHj1VgC+)sKQ?E$K1 zM=*B_CyG@0V<)Ul(;+xkV2O$~+reR-$JS?%dSmnTXX_$ZOZF{>p~?-Z>8g)3$ICf@ zG+NmJdn1lhvKtfD{^?EWI5^PM&1XOA@9*b;`*hC&V-eM>3Rnw8v&lyQdN+KUDT2_? zk8b>fuWZPKA>m}ZqDEe!yfsf?Kr`7+YrZk0u^$95kf2)i27zdK zofHa*SHEFf*0WUkB9?}-nAAquva;k zg_Ur(AMZE^bWH@}8&0KOtADNe3p#jRpH z{HPej5GObB?n^W9!`glv1z{?J&>&xdlM4aU@@a{Ttw*W)DNpta02Z^+xrVzlY(VA8r+0_Ptk4CzN}gIWXtr7l7r}7S*r%0toGin5NG>yr?Eb7*U$q% z6I7)AQHl8hkT{4OHl>LS0Tre{5+A=snS4BN{`mNcLg~B!@3xy5V{J^ynOd&L-OFdO zx(W`tzh*uO&hV>Y1(4A$BFX0YFQ3=h{xHMWN$WM=?v`ind``kF^*oEX6Y(8DN*Aii zuMx7}c8QF0KAcas+O~-KZOL zujC?E^W_@4!E}uQdnHS`Pxf)T6y1e_f#7`1#E%3r&C+v|LR(NF-n+bZnImYE9er~AErSBcCbrwRApb}0NY zC#zqUshT+KduvRA%%7uC-C_)1U10CxWThO0sr}@=)QfJ}m_HmzsZndhx)iz8BevNb z>t^al+(n3D37pnzNiy7W=?nF>BOBdUZEN?<;J<&rNKRz5HC?gSzvv+D!W<~?_$kBj z@S)1~BbEz{&P?t|lg2U3s1}oKpiUZY!2oT!d^6G(FCtp6#zUhjEbj#|*Sns9LuhB; z$qCJ#)zze95eTtFUbLAG<&KUvw=2H>{@V@Y{Fa$lLAD~i#LZ%`7+=t^q$;IjoN>YY z_o%V8aRv^iZL>F_K&ocna=kq|tikp;&W7yyuo8Zp&#TgK*W~^+^cNp!11X+|Xw@X? zuYx{-u_j`kN$cckUNYw;x|gyL^)(Evof2mcfs~ZF%S;eP{x+!FK*>5Ae~l=Yjp<(- zA6`^$?Prh(^>TRAkHOD+Z{e+d?&kJ0_;COesJ2HkgImFW^f0Q@42Hn{;i$8SKyx;5 zRo{KpURdm`drPjB9NuH5d$cgm!2`%&*UJ69IHJ7p^JT$-8y$AQ80~Atr8eOo-$O15 zt}tN*8#aDA5bSU2t+mpvu=Wc&O3+Kw=~QI2aW%shr7P>_6^}@)UIu;f7X^soJSnE+ zE9DpxQsDtkAR2)r%~}KT;#{xhM+5L3$zf~wc_IDV_SH)%u?#yVES~5dE>9NJ0ki^n z#_>Kmxxu)D*Hd2qI6w27NM$8Yt%GGwokZs7a1klu@dH|n)%2aazfp1g(* zEEZ(0CR}994$ff$w>*W}uqYnYYWY8sToNE^mT-rWm*1A{$C`}imKei(>A$GbARK(A z)?v^wXGN8R zZUk9-7_gOLx<*l-kI?XTw+t>Lf51g?_5KLXd~U1@F(kluxpqhpQkT9xfdd6NQee#;XOpIVzk(^8yW^)yX!slM6MAvq=>3h!?XzY*DeCT)J_k9^l zXEp0xZx>pRm(+S`c)7*Y7h^=0AKscf@z?Vc$Eg z8`4?f0`An`ZX`)v%1nJzsFPsxM&g0UjF#7~Sn7HiqQ^}caxiFDv|Z-B?2(tmdjs91 zE0d}nTE9k5U6pzY5+iD~?}oZCP(2Z4gi8&Fad)kSdq6lu3djSRP-EFb#_q{Ni^qsdHXba$a z(>N14j%mjb%h)p`es#^VdUsLzINwxNnzawd$Xm&(Sf6Sg!ETxeci4%E9HswuNUA%(_I%KUejn^5tg{0U0d7@eaO^*VWsEsRB9Q^pn^w=}|(%GE>!jHe0ht*8Gf0 z)g*JRWWn2Q$@c?>C{i$u)#rZrWf$2yBvd1z&uE^I3ZyBSOCW>jzpt78>efm$%(hGV zjyI7a<(tEo;YC3Yds?+hhK^ohg*wjBXP2@dt7zHdVlSucedbF>EluVjk~C{x85BD-Sopp_gY4Bmj*-GIJsQXiV98cM9h5ves#ZIdaL514e2xZ=unGEg3QfiiX2JS zAZ|}r@IQ(nLSMXd*;V`>PkK>87hdjaFczE^|$ zXCh2T)xM5x-^IIpuRM{67Yw$S4aMC+>wgM-C&vsz=SAOVrFL(`_ zybJat@Ta_`(qu->xUT@OUbu6b?dQVJ#Sb|Jx9p{iy~V#=JF-}wF#?hBxuQsAqJKIZ^RW*cU(3@<;hXw$KPy$jr21|Tv-@GO+QyB9QBSECht46i_z)@c zwd%}cP}7@CTr4&58DBL*%(nGX2ig-nqPQkuF>)40)wNfm&tcM6{5Zb2#oADF3B+N% z$OE!jQmmEF#sI-O1E)uv=XZ(d34OuswRR(ERGOCqRnr_sdLOkvmdobVwRrmOV(My$BG zzI;V<>7!Y-d&l`O_8!~XBCW8GZ9UX=s$=#k&DW8@;ss?xf~B%$#(bO zej_%Pb3Zb969!@gFDu^rF0)*@qpvDmVurvj-g4d(TiO;F^kGWEEK8KUfQY#X{QwGW zXJSk>sJEnxd_~Wt+IP2O7cXGpA-`%%@n~{0mU#T`=jC#0Hl2kI1*~)CX#gN6Dzk6H zH$yq%V7ym0OzVJTu=%0ZO3HQP4_5$Rk2-D&jlysTQY!rqJIQuuo^gN7U|E#JX z`71_E994Wj;L)h;R92w|Ak(<8v%MkeYP7&kMHx)t80-)KeF5!cF}M#~yg9N8jk$_n zsG9D9gB`shEAfk5#ntRl@^@}lQmYSaP1kOsmvL@mPJ>hL9QXAB z*wQg6D}n=?accK`D6G;y_Hf`|zRjhbI2dD`o@wVWxR9|X{k9@kIu@{1&%%#-nU`IN zTk~*pVWuhyn3%rGEU-0Ks67AHq5pa>|MMDHqt*8+ zqW$rAZeqB4+&a_zo&NCKPI$bcZ|^F#f>!kPCs{ItYgfXvX;v-l*>_DzQaD?QEygW( zh$&hn7-L9;sOZO;JG{_Nv5CAGNDY*hp7@cfPr3s5vP7{gGW2}ke}4Je>~P*7Pmg-- zxp3o&!%9yBxHQ;TXYYYs5}3e?W;qZ#fdOiqVyb#+&WS`V18lD6y^_sf0jzCDPg*`! zQ+<`yc_(?`H7YfKNa8?5u3Pwe-SDvMgr07ruV|InFCvVO)H8=g!NdB0Ad=>hvK{s*;)*-J31&LpUyvtfZT|Qg>O5 zucWhmb@sslFfgP zq>U#gaXKXM*`VjS-|!H!Mf9D$G!e5@f03pncPkiQm&E#A-RMA|Qn17=;;o^tmBNNx zDBIEyK z>5tJrKF9%=?xKpstucfBG{{tWF*3rT&*T$^1IVR-D1D|z6-)dY?+0CZ^L`%>CH>^u zG6JyS{rwrOI9QjN)<`?-D_kMvx-99mPFY05($1Aw0UUX96pY#td4z$_PpQ@|Tpc2v3 zQVZ@r-GCCa4%1MgO+iIRt*L#@xkB-Oru#~3_?pMS1U(Wssju*XTfF_lbidpb@9KFr z_F89Zaz3-EFrCo=gKpAe=j{X*Gq>E{MZc(DcGK_5OO7NdYRAtR`6_9D(b z9cCYy(e87By5tyY38aTPJKn0E;26IGMoYUK^b0P2&R%7%{oTb1Lg;u~PG$oujuM!p z=?AAu#mS92)`4s@i8b%om}(&}Z{GY$kAT7@egOSjh{NtOG%)b) z7SHr0lKHJRn#qp`Q-x!B(zc1dvQw6Laeh?Z_@Lpx&(bOL4}pm-PVIfXjoguK_bN3FXLKtsL7+Pshjza5f#mZ~dl98;dIP zIiW$T8SAES;(uW(s5IJ!8V04o446j<+Bz)&Q-X#`XbH3*L8(s{X~IIVY17kkj~{X{ zQBejYKzCu$15-k|39cb@W5L#Vmj`=ofZ(^vhP;+c-!ST*wFptIs+`8fB!*~g!UT3@=$NyKrSHQ-*v ziAXx%dB}O$VoKpND~S^A%KuI8){2<^P<@?*^TwkKng}4Z+B@k3Sc3nvvHBzoH1k5g zwmf!~&6G&{&k4E#r8!nF>#tail0<;?@l<_LO5m)^-qdbn%f29Ej0c`7pi!X~avw3D z3H&G+#(~A3#)N4wRxo$nB_PdnsUUka8`zSRdRkyqN*HyY^19OdiH7A_w<#r)>!1sYiq%-JC9-8-27^ z>7S})K!z5q)`gd!S zsl+_B6}s=H6OHZ9$i@WuQB8n0yU0zn&-Kp3)YDNE_2ua)&QQ|1Q*`(We z6Ez!h+WAmMx3w#lC$(r9KFK3T0ayOBrqk968LFKBwrp6Q0E#r&K!~ROYUz?m{0EdO z5Nu0@fAl3)3obp*a0s`5jX0UV)z; z_X;TqFs9lPmql3mDm?z8f22RhOoUh;=Z&>~%77lyWE zKL`^3KM3mZQ@koq5x$#wHSZ2^a+A1z_9N~FWJhEex|E?)=DZ*3rxxGYZ$4T%7xWU5 zEjVC<@04Rb%4KI~SNiQcbK?WYVVu#_q%U$QsOG@(Kl5&+q*jEH>JM8_Qw6F{%lAwF z)H%C(gZ%)Ntabqo3)Q{y!D!lBd{Owq**q- z&hC*>ErrqC|NIbUO7GDqpI;4BUhLFVQ!T*b)+dcTRYL!hghfPPfO2oj$K{qO_;g@S7B}B_bOpH{>SBDg@h>Yyb1N}Aq`~c!P-#%53kSAMocc44 z&u!n5|GKn$5Oz|pkn>rF52xB2v3vMy{hMsR!u!2Vqs6J>{ckf-&)c^qn<;3Y0LYm> z4O**#7r-{P_SdQ{2U1RM>Y{Gy3~Ml&rEfjaYXQ4aYubNxbZj@uZKqBkUVM3zabwjx zc2t^$XaQ@@ySwkBW+IX=8ciAK*z}CJ?C!0x#(wzuC+2rgq>=Md>xS@&qv!OQ14H&4 zeqq?_e*HV%H*_+v{(?ZA^p9W?tDPB*pzRLp7r`W8V8p)Y50PE_e-`hQ?F^t6Xo5+6 zQlDP*d6&~k&%HuEf-WhzNpk&9&_YqXf*@+Fn%y_|*WaJ^2W!%ui~Yw1UPavs7xI1! zR05KE&N2vtDm{6EEHvH=nX=dIK}oF^m$8y#nroemwbeXlC4dG({gwi ZmA_!C?x3Up`Lw(AnzoTv^;PWs{{hCo8{z-} literal 11140 zcmbt)2{hF2_y3HsZ<%Q9gbXqm+4o{7J7Xt|v1Cn1kuX^@)(Dv_B_V6ZHnIW7DgUM006+E ze+_{I06-`J0C=7r1ONc$6Tgm78Wv_omb#QC03Zzjr~&{KA`gIc5I|KGpy~loO$8v+ z03HNDlQN*4M4^xdNJ|5d($YcFNThU76@`I9rwZ^;J$+JDJv>xB2&$=56e&o6G!lSB z1_dG0$S%1j9;b3uJx)neJv;y&lqWUSgFx^goCHY&f{=iqDstrMlNlLV9+}w@Sw5Dk z3ZO_#^#D)^Qwgc5Q-sv1snn^HG$bI6oR&r=Q%Ex-GczO0DZk9jj?D7%%#M}vzx?Hq z<>g~1RYB5KX=zntGQk5tNL3|Fr4mjEgp(=4$w|X=Ku2a|2Ze3rgq(&XR|S#D6p<8K z@>qG~*vZLMs_N8fpM8Irq}=M~A<9T2Bf zmpe$Et6V&L;WpocUWV>`(eH(5N2tUl?Ob|}3r{aRzU3kfoSnUYV)Jyp@t|b+wWh}% z55lq5C+e({($e)%!iOb$4j80dNop({1RKlk!c^O8BmiBj?{R?Mp?Ab-le(78nYj`_ zzT|>K%CE|zVAhY&MB561|2M*mAk+sBDb@g1gUDXjE=(v%XyD3{X=nu({E;$AMf^zO z(K*lpI8K?aNNc_tO;n#{`;ou>O5K_nca?9JFXK3)+glr>v45QP;aT{0^wul?lSLaH z75?2NAMHQ8afw}FIx3E)*e@_h`(EOU+)F5qqun1H8i%`A-FNQqjcoLD;8jT*5eGB1 z{tv_cjDcqxB7Y~|_`+c6He$dgRxcuE_rX z{hKT8s@ukRFw>o##cOpQJA!L6ZWlpiJ@V0D#a95d8?N=g{0?JRlk@7l*a#p=Ol8%uKGi*B zx4DD+qi`n!>hbkN4I2h3Yizm~;N8e}-~8?_VW)2bPoCr3PD$$JgMdw#)4$uVRvB>A zd8ipXgLeNsumz_}l6M4?qYPQ3DZ`lrg4kBhGu=wL*A1S66YK8yyu*ZY*m~P?n;DNy z>mu0Lv#7v{VY@Ex73tySrtdL+^r&D-4{d8lM$c~a^Kkm=WFRwFAdRXDtP7L;QZWyL z!UE9!MKBWbSvY%jGJyH%2l{|O+wvm}@MTUC5H19g-9hr8^*^*dn4kqbG~pc1tGgAz zL|C6R3OkS|Z5v8pi-5%1BxH(bh(C>>bpj7}qY&o(h1+q&FJ8xAbV-;3*!j8sVg=4C9 zrahp2NjL$r8dB|W%xOSJ;lpcxjsiPKfjb%99-)1)=Upsl$j@+R?x2;N9gNTE_6rNZ zAPtNvtZMA3rI@bU?A#PfunAwoz2Ufqlc!}!WF3O~W@E0bI6W}I4B6Od1dMDK%$EV1 zUsBz1yaR5e>075}2A?-_V)Q0?yJGOa*(&vXA7dW$^TZ#zF4b*=!! z{TX{|0cAupxLlehiNrI94+=!pLp-zs5-=Ujid4Cym02Dytc;XVQTg6wZPDcihm|!O zZ^9xDO~m=R0S@)Zh{B7X-B4oJur)bn?0O}_Hjy8Ds9a_O^hr@LGqIj}t~p0}XMK~h zF~clt$*HM4hDI7S>E;JY>xr@xTJc&_!T9K|r^Fj>T~ zilvBO$EIT(LacY9yne9w`F0z^jI&ytF}PaAm=b_M^l3&yT z9As#HXECKQOIb*Ixm?K@clzq5t2j24)k@|p3HOu}VjtNmbr?e2iU#Mi#tckG&EAkVU2V7d%edc|N?v8bADAFd~yd+JW+%Q=|8tYCC%$)(zN>RI50Ab(Q z2)$IcwoWgCc)wwU*l`S$T@EtCwz6C=?(NA56RUl!e%O~yIU_}99UE^D2vG=d{7j6V)J_qJpK%IAS_*#TW z78!V%QhQbER;L~arzs~mM;$6zRDW4%2#ykhuB;2AV;+omII{<}uJQ%6OLTt{8yD^n z{>e2Y`UhDA+nb~k_onqtz`LK(PB zxUgI^*vcsO8s!a6s4p`H$g8BCrwmn!p}GE^c|sltqMMSe085&6T2JKQs`r?{id4k> zaB~-v?dXM~rePqlKb_y428y!l|a|B2;Cq zlQcbNp5-(!Hgrn5GTqrfKAQbw+Q#@Qnlg|1LX*6Ao}70Xrn*qEds7K5vB#VCQ@jxp zz;t^`{geKo@tB=0Z&+pGs4fp^qEwUthcUbifzgs(AHFy*o<^|F$Bd%t zbLN3k1(4qk5(qH^?gq%hJEs-G%%Ky0{gX@=71`Tm(&6%mguOHumirL4UdVNVe zOt2fz*%WRFmYPM$3GC3=2419tr$^UUkH3+ER~rX%R$nDVy8Bi66)ocbBIbQNaxS%F!Y zieddv3il6p(sp^T6hX`fdmtKAW8h*(gQI6wQk#Xi$(Ly_%(RgSJ3yS6Rl=8e1!D&k zZ`e@EWfS5-PVqP(>~Oom?W|F{;j1O!F)lV%yHp?yNhS8~zmd^1E=L+RqDyz>yC+4rnM4r`R)O;txla zVL-WAQZ7$)4?`}`OxxW6F~#C)Wn@r=*AdcUSNd({NVTtmNSM-*?Ld0yjb$0lAHi_l zJG)4UeMSfSHAmiecrj8^Gs|S+{05z$Nw<9H?9=xBC@AMQb9q zK*Kx}%4H_nzsNKl_^@qPH-Him2X21WYTMu1l1_tPEmg>-M)mSC{xp3)jXmcnLg=JV zH`%(7S80s-ykrN|QQRFar=AG(XnTS_m$~t0teWahYve0&pkzwCgU%*0QQWg#|BQ)c zol)W-+Y61`@2?r*lv#a%Jk_4WC%Z}ECd8=p@TVqC+CD~nJ#F5hrsoAEDVpCEMcv>9 z;C4OuJ;4IPcqYD^q{8^QH?m1XbIQrz_BfQnuF^w~a!dCHZyQ0Pf@=#W= zLk7&um8L=(Ehiau7gt1gL`xA1C#`XRwu8^4 zpk0h6B%!vGc}2&{MWC@{QXs&Z)8l?Cotp^bhSCpeqtKW86@03W8x}cuV*v98*c%#G zrg}F`h+s1Wlkc=-|7}-RHz)Xk4lpK9$@J<)c^o!&vmeeO2=C$*8nQX@I);N!`3lax zWS}yt!0^jJM&$?1F1yh1yPZXUw7w0g zf2!S4&*T=ju!-SMf2C|jWHopLZAb?>C=hHE_&u$9<>S^ujG5(2pGywL2TeJ`QGFR~ zdGm$=>$uM)4GnDkJ@X0(iObz;e9`uma6u{oC^JoY8= zL$ywoga*GlgUmC5KQ1Vc-etYsa*2v`Y!iE^42}^pxVTSbThnIUxtSGDDnD+#XEX3v zQnT6}O1_UV{tc#;W><{8jLO;Z9CBaLWF7L;fj${aa zl=s1tP=_9&PWE9#(}mkMCPowbq;LhTA)T0JhHpNc3tM5o4;8%Sdwgtx;`jRa zpN4&h4{ZR@Be5;&XAH`&Uu7h-$sfPY)Uy$=Y@a>4|Te}+_9oN=KdKHnx>$TSk zwN)(l%Z&1q%@2+)20YGP3|6-%z7P9ZM8#yEbr<1RZ@;c9^=13~&$fF2vOv5Q6yTD{ zFo>-;F(TcgGcxHvl!;YUispoVeJB1g9&(}RyR)l|p9Qc_W;3=AFE4e4#*a}p;5L1; zG$(guMm^>`+c~B=AsPuNzvH zE~+FB_X#?T&1xsc~8TlkIU7rgusC`WEWg5*PC?w;W7)nzM(MN6I%)$92v*L zyPgQ3ZNvgF!2T>B@3ufHn67djG57;U23;h8$D_83Ej;YEaNLB^25q_PNf3BI+0_p>0z=~L-Rk&EgxC9OuxG#xCC^nx z3Ps8AB{ooQh_ws5&ZxH-9~>YfWy6SvRCn&`eiMCk$MGHG#3yZ7%uaAwMhfK5Zy}M& zyu1?K${`>v|6A!%dGr%~9j@=B)<@zJh`A{l4mf>P5sz?Pwm{4DV|jV5bEw(t2%%-4 zq2o&><+R@Qmh9c+W?rd}@y+?APz{iEtapyDm-9whzSV^NKpvh!$UZOeGLfyN?z+T# zuZi1YtSxNBw!=?V=+h`KIjQt^7T^QPmTMbIx=Sya&Q=##K@(zpK+8M@xUTUaHU2&i zSnNSqfB#{W4(VanfT4uy?RjBsff!!ZN0txOts_AnFf|2R8IL*Z0yIwhI`?GVct!JD{pINK005VSOU1SKQ0qQMJXbcy!l<{jen*E*8&1-MLvqjs(Xee__bnU z#m@fh&j-bvv(*YH5*WLw_JBT-iX7ZNU?{bdY@QjrW1SFs%X)L^tNk979zS7+putB^ zlk7LW!|Ht&p+eM720&Answf|rk_I&$|ts9L^ZJ_Ef+`2f7ug%5r6sMBI#JsmSukxpTR%1;Pn;H?Hlw znKHKFS*`aUfTU(yQ?v`>2g3WCVwijQb@NGFv`>7QYNl=Qtd-XYqvCGU4tgDT`F@fL zW#oNEw!Nfm(^U!#h_3*5e)gCL>~t|XGkMh68QPS`%A8H9qhT_=y6WFik!I_MendH; zl!I4tamf=HDV)Z$>z|@!M)Wt2?9UjRG3{Om%yv>CNO9vmY6|T=6>nXPBmw1mg!3~4 zhEaI^ltKwb5+Y2|fju2R*10OiqwyMl*bCy4cZ@demkgiGhcyaCGvE_{51p9>Bxb&( zDTibyOC)bq{k}l=*;B=#Qi6(B5gm^aQaMaY!=ZtiRhin|@^*mepV}ih&i-WDWbO`4 z)Ay>l!My^4fjf0Y#LGrI&voehmp@y&0YCR8?`OB!QhBOCFFj$WVb1O`p=BQU6nBGT z1k!a|>_Kgy3Pn<<(#`x8MrpClRCKC-N-n0}U*1MW!C51tuNR$g>}Gmw8;J&aJf27N z!J;uytaXd7{A6JmKm6)Io7A9osJF4$bgdda>dr5#>7BAx)1aAJBH^yeI%jACl+L!K z2}9KDYxGvPBT|e>OXO2I{F;QGPWMq`uzjPevpwT4irTjNvQ^e-IHg8tU}|;VcDNS} z01Xyo0qZB;QQ?NL?m*_APva>h7?Y=6DI|OaEKqO~ku|YcUXR%fYt9C&72l=H1l?)xDdWa}eFeEh4_{#6n` zdK~KxVh%Gmy1Du=fs zI)S?nn%Rp!$Amg}JgU_792>E<*Yn4&7V-y>>F|RbHD;H=?w1uPl|hE?nI@x6J!S&- zJd7KTnjY|NtDv76;vCW+e)fEbO7R6p&5@c=ObE~2Th0F86ffLxexQh!tjNJJ;#z-^ z1j4z&7t?v=#g=gqdKel;M)3mX66hQVc4>#+#s9hA12*gYY0|kdtSGSfmyA??4~vHO zj^;J10^_H$s2nCP!rPg(&nweV-`G5{^DtCr40|e`n#h{GYh{Z`=k<;`FlH;e5IXK{ zP0VneLeSalm{zoHZU_|Gn0Om_fbE|G0b}D(3lj)P@as`e7y0UE$OW%vg#zYcJ@zdwQwoq?n58qf zZ@9>#^J4XJv(&saG%?r~t9U7?h9ip=`+-36Hrm1zlQ;dx^dODzHDvYt-ly&LOZnMJ zTyQ+tklHA$_#9Inyo>7oJD<`=_l&TEj5-l5u|*B2J_1+{)^xML52PKDhNH zoo2`V94&L3rJ=1RG2hVuIbjGR@`@ynX6;T3Q6*AmxBFfAl;z>l17BKL_;Ciu)y6+~ zmq`vaQRLkjNcHeL3GszJsAbuR#XNw1Ehffk^`yjk(<0h4Tw@j3{k%e3G4q^MT&>aXG2gNt$IR0(bH|?>xLMH6 zz*A7YE7}shgseLh*d0)>E}wiQ19f@B;X!{Ff3-4v=*2iL0&uYb2Nnblyhniuzr z&ZzItW%FW+s~a)F|LBqn*^E0cGGXDA!sG(O(2-LwLC??CS4!%zjCx{1?JiL`|PU>j}v6Q?Yz0Sr@Hy~44G(V|){*#S6)tb87zgm0;=*Yc#73=wYr7niP z+;}yjKn}6tlet@zpJzUBk%7G1g^@)q7JWS}h-Xq#t_u|kvJ62n*uq1kZT*H4EN7;M zICJ=)F>wH+87PJDxy#jX9@#Dg9f$>Nyf?q>N9!??`2Me`*M{!;E57!1;w2DOm$LYG z51HU^jl(;v!dVU@Tx*tuPAJE&j36-Oh^@4I38k6%osBn<&nr)g>HJYAas%oO9UOhG zM1oA5YpYgL&k^_|Ge=Vq_6;e(m^mNYXW;*w!2MK)Qs)f2e)O|Lk`}g!ynERuC?KnobcF-0Rr9Jp|!6$PwSmsU9V~w z!YmXW{GV$|Khr`H(bAdE#<=TliJ_~FR#{?F*u0k7wAUNf`BNb9!o(6qhxuDcZQL_V zr`6rD&E3tPQ&o@ZzWj~o5}@LZmG=;!;XD!O6*Rh#*-w!Ld$*c0crO$ z*zsgnI6pq!xs6g1R=N1PEoA6e?{82lbb!BzsT+7hGigk1y#0w$8F4tsGKYNdY1_%i zu+PG$;~?kzAV96r<|Mz)<0$tA1s~9NiUsgu=Gs9f&yJ|7!LCmuE-!~2qp+7Z%=QFu z#xg1tWs8^Qhqm^cu)P=JeA~J3i-z@OSJ9Lz8|y03DR?)AeqG(NWyR1r=1{?=4N8X# zltwSaF5Jv^+Ga6q1y7Scg!03FBNHyTJOxJc%8ENa|3+dQQcv&SZgYi6g>5KZB~&SZ zAd*@rn=|v;zkbnm1*yFH`9}GZPu4STJjT!~!C+;z8UtOmu@!mv5GxQm5pt%#m~L&w zow(S08cv@A7VeF!%)}ex;6R2T;MgB}s5b+-n9e`v#=jJe0PAE@XU)Lb+{<2h9TCDp z2-hV;C&o#x`6g$#FG3(ao6|)~rp3=Qiokl}hbW8AK@mN}Mp<#^2vy?&?Gpqi+N9Ti z-3g59qm(xCF|ggU#^AHtG$w#aQH_}I`f4~2_N(2dgzDxaa=sY2T0NBu#tPo#WURS zrXD_^98_&fxs0eP?UP&b9VJLm;%CmJ#`^3BM}Q^GW%<#b8wpKTob@1hNYe zqUfi2*GQMyB}oS@#$}WK2LP7q8|k#Lku2tqZpTtM_ydP(d9izv6eqj69NXF+!P4OL zCU!l~yQWdf_H^jJv?w^b)Wni>9RrrOlHAGG6VC_Gx>kT^n;8R<+os8$wHk$ zTXR)+DffmoPcci5T`eY9p*F>^hY=euOaHzt#i)Q+1mWZNwkX0FZ6p5a=@$?RGjZuX zpsD!8G~dqgQRUys|Kv-Yypt4SM9Fsyx^Ka(&VGwIRX`Pmgt{F8Jd&E;Aw|I^T<|6a zkhNhfObOMyV)kiWY3RS98`Vn07?TcG*L-#n1lwz)y&Tkk134xnQQEIwNt1=eDh%rzYY46@Cg+)`VTAq|Zl;(%d?b^Y(JpOoDHXNzR! z>YxJje@_PTA|i;-BE^p_A}uw+r|W{97WSUOfCr?AZNWHds9u!T41e*5!j^fuvohNz zu+#ebVi+>0k<*W+#ln!SUqli!OX8`lB~G0n#sUA?%cj$0}>lIUL9IXv%37^x^={H{!#Vu%fEX*9Mzb&%{lA298SiB{@ee}f8a0A zi?G>`}Rv zW{569JW_$!!yAh409BMAKtBZg4XhW#1^xz&djl<^T8I=Jao^I5NYr?;(K(nq z+RiguSn>?`JbjjexT$~eztAtrt{9hM_Dwl%J*Cg209<^!t5zVC4pddZA28keKbvig zH5Sui;0r*=rV6{K@Jd%I-EB=CJz+qy-9?#}pmvFZ0Xdfa~x!h4J~VGXZDOL|RKa9I=pk@y~AZA&ZWcL;ahV^{5@X zyQ6jp6=gN2{}9O~hRn5HUb}ZPw!hZi>{u1vBV@CGCYZM