From b632d512f7c72059bb6981f9fd51bb089423d5fb Mon Sep 17 00:00:00 2001 From: hadley Date: Thu, 18 Aug 2016 08:37:48 -0500 Subject: [PATCH] Function proofing --- functions.Rmd | 82 ++++++++++++++++++++---------------- screenshots/rstudio-nav.png | Bin 0 -> 15430 bytes 2 files changed, 46 insertions(+), 36 deletions(-) create mode 100644 screenshots/rstudio-nav.png diff --git a/functions.Rmd b/functions.Rmd index 8c12c51..5df09ef 100644 --- a/functions.Rmd +++ b/functions.Rmd @@ -2,7 +2,7 @@ ## Introduction -One of the best ways to improve your reach as a data scientist is to write functions. Functions allow you to automate common tasks in a more powerful and general way than copy-and-pasting.. Writing a function has three big advantages over using copy-and-paste: +One of the best ways to improve your reach as a data scientist is to write functions. Functions allow you to automate common tasks in a more powerful and general way than copy-and-pasting. Writing a function has three big advantages over using copy-and-paste: 1. You can give a function an evocative name that makes your code easier to understand. @@ -13,9 +13,9 @@ One of the best ways to improve your reach as a data scientist is to write funct 1. You eliminate the chance of making incidental mistakes when you copy and paste (i.e. updating a variable name in one place, but not in another). -Writing good functions is a lifetime journey. Even after using R for many years I still learn new techniques and better ways of approaching old problems. The goal of this chapter is not to master every esoteric detail of functions but to get you started with some pragmatic advice that you can apply immediately. +Writing good functions is a lifetime journey. Even after using R for many years I still learn new techniques and better ways of approaching old problems. The goal of this chapter is not to teach you every esoteric detail of functions but to get you started with some pragmatic advice that you can apply immediately. -As well as practical advice for writing functions, this chapter also gives you some suggestions for how to style your code. Good coding style is like using correct punctuation. Youcanmanagewithoutit, but it sure makes things easier to read. As with styles of punctuation, there are many possible variations. Here we present the style we use in our code, but the most important thing is to be consistent. +As well as practical advice for writing functions, this chapter also gives you some suggestions for how to style your code. Good code style is like correct punctuation. Youcanmanagewithoutit, but it sure makes things easier to read! As with styles of punctuation, there are many possible variations. Here we present the style we use in our code, but the most important thing is to be consistent. ### Prerequisites @@ -26,7 +26,7 @@ The focus of this chapter is on writing functions in base R, so you won't need a You should consider writing a function whenever you've copied and pasted a block of code more than twice (i.e. you now have three copies of the same code). For example, take a look at this code. What does it do? ```{r} -df <- data.frame( +df <- tibble::tibble( a = rnorm(10), b = rnorm(10), c = rnorm(10), @@ -36,7 +36,7 @@ df <- data.frame( df$a <- (df$a - min(df$a, na.rm = TRUE)) / (max(df$a, na.rm = TRUE) - min(df$a, na.rm = TRUE)) df$b <- (df$b - min(df$b, na.rm = TRUE)) / - (max(df$a, na.rm = TRUE) - min(df$b, na.rm = TRUE)) + (max(df$b, na.rm = TRUE) - min(df$a, na.rm = TRUE)) df$c <- (df$c - min(df$c, na.rm = TRUE)) / (max(df$c, na.rm = TRUE) - min(df$c, na.rm = TRUE)) df$d <- (df$d - min(df$d, na.rm = TRUE)) / @@ -55,7 +55,7 @@ To write a function you need to first analyse the code. How many inputs does it This code only has one input: `df$a`. (If you're surprised that `TRUE` is not an input, you can explore why in the exercise below.) To make the inputs more clear, it's a good idea to rewrite the code using temporary variables with general names. Here this code only requires a single numeric vector, so I'll call it `x`: ```{r} -x <- 1:10 +x <- df$a (x - min(x, na.rm = TRUE)) / (max(x, na.rm = TRUE) - min(x, na.rm = TRUE)) ``` @@ -85,8 +85,8 @@ There are three key steps to creating a new function: Here we have just one argument. If we had more the call would look like `function(x, y, z)`. -1. You place the __body__ of the function inside a `{` block immediately - following `function`. +1. You place the code you have developed in __body__ of the function, a + `{` block that immediately follows `function(...)`. Note the overall process: I only made the function after I'd figured out how to make it work with a simple input. It's easier to start with working code and turn it into a function; it's harder to create a function and then try to make it work. @@ -108,7 +108,7 @@ df$c <- rescale01(df$c) df$d <- rescale01(df$d) ``` -Compared to the original, this code is easier to understand and we've eliminated one class of copy-and-paste errors. There is still quite a bit of duplication since we're doing the same thing to multiple columns. We'll learn how to eliminate that duplication in [iteration], once you've learned more about R's data structures in [data-structures]. +Compared to the original, this code is easier to understand and we've eliminated one class of copy-and-paste errors. There is still quite a bit of duplication since we're doing the same thing to multiple columns. We'll learn how to eliminate that duplication in [iteration], once you've learned more about R's data structures in [vectors]. Another advantage of functions is that if our requirements change, we only need to make the change in one place. For example, we might discover that some of our variables include infinite values, and `rescale01()` fails: @@ -151,7 +151,8 @@ This is an important part of the "do not repeat yourself" (or DRY) principle. Th ``` 1. Follow to - write your own functions to compute the variance and skew of a vector. + write your own functions to compute the variance and skew of a numeric + vector. 1. Implement a `fizzbuzz` function. It takes a single number as input. If the number is divisible by three, it returns "fizz". If it's divisible by @@ -179,7 +180,7 @@ This is an important part of the "do not repeat yourself" (or DRY) principle. Th It's important to remember that functions are not just for the computer, but are also for humans. R doesn't care what your function is called, or what comments it contains, but these are important for human readers. This section discusses some things that you should bear in mind when writing functions that humans can understand. -The name of a function is important. Ideally, the name of your function will be short, but clearly evoke what the function does. That said, it's hard to come up with concise names, and autocomplete makes it easy to type long names, so it's better to err on the side of clarity, not brevity. +The name of a function is important. Ideally, the name of your function will be short, but clearly evoke what the function does. That's hard! But it's better to be clear than short, as RStudio's autocomplete makes it easy to type long names. Generally, function names should be verbs, and arguments should be nouns. There are some exceptions: nouns are ok if the function computes a very well known noun (i.e. `mean()` is better than `compute_mean()`), or accessing some property of an object (i.e. `coef()` is better than `get_coefficients()`). A good sign that a noun might be a better choice is if you're using a very broad verb like "get", "compute", "calculate", or "determine". Use your best judgement and don't be afraid to rename a function if you figure out a better name later. @@ -195,12 +196,12 @@ impute_missing() collapse_years() ``` -If your function name is composed of multiple words, I recommend using "snake\_case", where each word is lower case and separated by an underscore. camelCase is a popular alternative. It doesn't really matter which one you pick, the important thing is to be consistent: pick one or the other and stick with it. R itself is not very consistent, but there's nothing you can do about that. Make sure you don't fall into the same trap by making your code as consistent as possible. +If your function name is composed of multiple words, I recommend using "snake\_case", where each lowercase word is separated by an underscore. camelCase is a popular alternative. It doesn't really matter which one you pick, the important thing is to be consistent: pick one or the other and stick with it. R itself is not very consistent, but there's nothing you can do about that. Make sure you don't fall into the same trap by making your code as consistent as possible. ```{r, eval = FALSE} # Never do this! -col_mins() -rowMaxes() +col_mins <- function(x, y) {} +rowMaxes <- function(y, x) {} ``` If you have a family of functions that do similar things, make sure they have consistent names and arguments. Use a common prefix to indicate that they are connected. That's better than a common suffix because autocomplete allows you to type the prefix and see all the members of the family. @@ -228,9 +229,9 @@ c <- 10 mean <- function(x) sum(x) ``` -Use comments, lines starting with `#`, to explain the "why" of your code. You generally should avoid comments that explain the "what" or the "how". If you can't understand what the code does from reading it, you should think about how to rewrite it to be more clear. Do you need to add some intermediate variables with useful names? Do you need to break out a subcomponent of a large function so you can name it? However, your code can never capture the reasoning behind your decisions: why did you choose this approach instead of an alternative? What else did I try that didn't work? It's a great idea to capture that sort of thinking in a comment. +Use comments, lines starting with `#`, to explain the "why" of your code. You generally should avoid comments that explain the "what" or the "how". If you can't understand what the code does from reading it, you should think about how to rewrite it to be more clear. Do you need to add some intermediate variables with useful names? Do you need to break out a subcomponent of a large function so you can name it? However, your code can never capture the reasoning behind your decisions: why did you choose this approach instead of an alternative? What else did you try that didn't work? It's a great idea to capture that sort of thinking in a comment. -Another important use of comments is to break up your file into easily readable chunks. Use long lines of `-` and `=` to make it easy to spot the breaks. RStudio even provides a keyboard shortcut for this: Cmd/Ctrl + Shift + R. +Another important use of comments is to break up your file into easily readable chunks. Use long lines of `-` and `=` to make it easy to spot the breaks. ```{r, eval = FALSE} # Load data -------------------------------------- @@ -238,6 +239,12 @@ Another important use of comments is to break up your file into easily readable # Plot data -------------------------------------- ``` +RStudio provides a keyboard shortcut to create these headers (Cmd/Ctrl + Shift + R), and will display them in the code navigation drop-down at the bottom-left of the editor: + +```{r, echo = FALSE, out.width = NULL} +knitr::include_graphics("screenshots/rstudio-nav.png") +``` + ### Exercises 1. Read the source code for each of the following three functions, puzzle out @@ -306,7 +313,7 @@ if (NA) {} You can use `||` (or) and `&&` (and) to combine multiple logical expressions. These operators are "short-circuiting": as soon as `||` sees the first `TRUE` it returns `TRUE` without computing anything else. As soon as `&&` sees the first `FALSE` it returns `FALSE`. You should never use `|` or `&` in an `if` statement: these are vectorised operations that apply to multiple values (that's why you use them in `filter()`). If you do have a logical vector, you can use `any()` or `all()` to collapse it to a single value. -Be careful when testing for equality. `==` is vectorised, which means that it's easy to get more than one output. Either check the length is already 1, collapsed with `all()` or `any()`, or use the non-vectorised `identical()`. `identical()` is very strict: it always returns either a single `TRUE` or a single `FALSE`, and doesn't coerce types. This means that you need to be careful when comparing integers and doubles: +Be careful when testing for equality. `==` is vectorised, which means that it's easy to get more than one output. Either check the length is already 1, collapse with `all()` or `any()`, or use the non-vectorised `identical()`. `identical()` is very strict: it always returns either a single `TRUE` or a single `FALSE`, and doesn't coerce types. This means that you need to be careful when comparing integers and doubles: ```{r} identical(0L, 0) @@ -321,6 +328,8 @@ x == 2 x - 2 ``` +Instead use `dplyr::near()` for comparisons, as described in [comparisons]. + And remember, `x == NA` doesn't do anything useful! ### Multiple conditions @@ -339,7 +348,7 @@ if (this) { But if you end up with a very long series of chained `if` statements, you should consider rewriting. One useful technique is the `switch()` function. It allows you to evaluate selected code based on position or name. -```{r} +```{r, echo = FALSE} function(x, y, op) { switch(op, plus = x + y, @@ -383,7 +392,7 @@ else { } ``` -It's ok to drop the curly braces if you have a very short `if `statement that can fit on one line: +It's ok to drop the curly braces if you have a very short `if` statement that can fit on one line: ```{r} y <- 10 @@ -427,7 +436,8 @@ if (y < 20) { ``` How would you change the call to `cut()` if I'd used `<` instead of `<=`? - What are the advantages of `cut()` for this type of problem? + What is the other chief advantage of `cut()` for this problem? (Hint: + what happens if you have many values in `temp`?) 1. What happens if you use `switch()` with numeric values? @@ -456,7 +466,7 @@ The arguments to a function typically fall into two broad sets: one set supplies * In `t.test()`, the data are `x` and `y`, and the details of the test are `alternative`, `mu`, `paired`, `var.equal`, and `conf.level`. -* In `paste()` you can supply any number of strings to `...`, and the details +* In `str_c()` you can supply any number of strings to `...`, and the details of the concatenation are controlled by `sep` and `collapse`. Generally, data arguments should come first. Detail arguments should go on the end, and usually should have default values. You specify a default value in the same way you call a function with a named argument: @@ -474,9 +484,9 @@ mean_ci(x) mean_ci(x, conf = 0.99) ``` -The default value should almost always be the most common value. The few exceptions are to do with safety. For example, it makes sense for `na.rm` to default to `FALSE` because missing values are important. Even though `na.rm = TRUE` is what you usually put in your code, it's a bad idea to silently ignore missing values by default. +The default value should almost always be the most common value. The few exceptions to this rule are to do with safety. For example, it makes sense for `na.rm` to default to `FALSE` because missing values are important. Even though `na.rm = TRUE` is what you usually put in your code, it's a bad idea to silently ignore missing values by default. -When you call a function, you typically omit the names of the data arguments (because they are used so commonly). If you override the default value of a detail argument, you should use the full name: +When you call a function, you typically omit the names of the data arguments, because they are used so commonly. If you override the default value of a detail argument, you should use the full name: ```{r, eval = FALSE} # Good @@ -590,7 +600,7 @@ wt_mean(1:6, 6:1, na.rm = "foo") Note that when using `stopifnot()` you assert what should be true rather than checking for what might be wrong. -### Dot dot dot +### Dot-dot-dot (...) Many functions in R take an arbitrary number of inputs: @@ -626,7 +636,7 @@ If you just want to capture the values of the `...`, use `list(...)`. ### Lazy evaluation -Arguments in R are lazily evaluated: they're not computed until they're needed. That means if they're never used, they're never called. This is an important property of R as a programming language, but is generally not important when you're writing your own functions for data analysis. You can read more about lazy evaluation at +Arguments in R are lazily evaluated: they're not computed until they're needed. That means if they're never used, they're never called. This is an important property of R as a programming language, but is generally not important when you're writing your own functions for data analysis. You can read more about lazy evaluation at . ### Exercises @@ -644,7 +654,11 @@ Arguments in R are lazily evaluated: they're not computed until they're needed. ## Return values -Figuring out what your function should return is usually straightforward: it's why you created the function in the first place! There are two things you should consider when returning a value: Does returning early make your function easier to read? And can you make your function pipeable? +Figuring out what your function should return is usually straightforward: it's why you created the function in the first place! There are two things you should consider when returning a value: + +1. Does returning early make your function easier to read? + +2. Can you make your function pipeable? ### Explicit return statements @@ -704,11 +718,11 @@ This tends to make the code easier to understand, because you don't need quite s ### Writing pipeable functions -If you want to write your own pipeable functions, thinking about the return value is important. There are two main types of pipeable functions. +If you want to write your own pipeable functions, thinking about the return value is important. There are two main types of pipeable functions: transformation and side-effecty. -In __transformation__ functions, there's a clear "primary" object that is passed in as the first argument, and a modified version is returned by the function. For example, the key objects for dplyr and tidyr are data frames. If you can identify what the object type is for your domain, you'll find that your functions just work in a pipe. +In __transformation__ functions, there's a clear "primary" object that is passed in as the first argument, and a modified version is returned by the function. For example, the key objects for dplyr and tidyr are data frames. If you can identify what the object type is for your domain, you'll find that your functions just work with the pipe. -__Side-effect__ functions, however, are primarily called to perform an action, like drawing a plot or saving a file, not transforming an object. These functions should "invisibly" return the first argument, so they're not printed by default, but can still be used in a pipeline. For example, this simple function that prints out the number of missing values in a data frame: +__Side-effect__ functions are primarily called to perform an action, like drawing a plot or saving a file, not transforming an object. These functions should "invisibly" return the first argument, so they're not printed by default, but can still be used in a pipeline. For example, this simple function that prints out the number of missing values in a data frame: ```{r} show_missings <- function(df) { @@ -733,16 +747,12 @@ class(x) dim(x) ``` -And we can still use it in a pipeline: - -```{r, include = FALSE} -library(dplyr) -``` +And we can still use it in a pipe: ```{r} mtcars %>% show_missings() %>% - mutate(mpg = ifelse(mpg < 20, NA, mpg)) %>% + dplyr::mutate(mpg = ifelse(mpg < 20, NA, mpg)) %>% show_missings() ``` @@ -756,7 +766,7 @@ f <- function(x) { } ``` -In many programming languages, this would be an error, because `y` is not defined inside the function. In R, this is valid code because R uses rules called _lexical scoping_ to find the value associated with a name. Since `y` is not defined inside the function, R will look in the _environment_ where the function was defined: +In many programming languages, this would be an error, because `y` is not defined inside the function. In R, this is valid code because R uses rules called __lexical scoping__ to find the value associated with a name. Since `y` is not defined inside the function, R will look in the __environment__ where the function was defined: ```{r} y <- 100 diff --git a/screenshots/rstudio-nav.png b/screenshots/rstudio-nav.png new file mode 100644 index 0000000000000000000000000000000000000000..927fceac5b51bde9b7be1954de117315a27f71d4 GIT binary patch literal 15430 zcmch6WmH_vwjk~ffdqGVcL^OdI01rNkl@gGaEIXTZh_$L!5tcR4bZqXZqwg=Z|;4w zX3gI@Yn`s@s_nIF*WPs^RF!2h(TLGtU|=xi2B65HcB0NyKnxLzf}-OPQ0PJ>RcmRY z&?*t%;Z}lwh1;?Dvocdh+Tot~m`Q#5N`%g@`6>?wZu1j*6SLg-aX#7U-sQbHF2V!D zC{z>DYpeOepH!-gQFu5tIW1|NS{?@03g%NDyzO$9XugW3rl3Ff%$vc60djalw<%rp z&D7-^Gqb)SIWi0_+K#3@IRo-eAIz36X+9$b%txX*lF4x*g-gBgc|aT`44{snbgZ_C zpmel0dYqnr1>4jQW>}M7F9{3AfH?O{k9?RD_+z*of_@jw4jYcQGxw5yI2}&x*ho5T zE#5oR`nmI+qBk~+BdgE+B0dCTF7c?#xw#DE7E`YhW0&1@40@v6VX(Gy7;26Uyp}r6 z248u4i{3lS>LQfm#^Y3c{E;%1ysUV6is9tAR^vm#FonHzVoedLWCU#0 zJZ_ZLlotJW&f@yUSHi~FJj?lN&aX=%Ccc`LStKz*j$6pju;}3?iikXt^jq$?e7AoD z9^@#E;c?}VMUD9?KU~ViMWbx-AQ1#gQ0)YCs7V+tGh!1+zND)Cezr!OF)u{qj%B?- z*df(Nzw2MWbi9XzPGGtJF_3-MQJ|^EgAaf0H#VgN{YXMauA2;1Lb!jEAwu#dA6Fjz z#F^07>=DcmY~L&0%NH2zh~6JGRL9t-HT&?XfFiBF2S6DBf(xrx4B_&YwROZFUctQ^ z5To=PCbQhp?Ll+*Pru(XyeHF}2y19G6I6fo<4w(>^7a9#Z%pDCe z`fkpLqd(P9u+;eL6`ah^-19s!QhN)@%5wNBFakqPDiY z`LOZfId*xu4+Btlj(O!-e^pG<>k#VL6rP(CKXn!0XHL!n=&2`9W>LND^IM4-37jJY73Qb=CNQ*m`&JjVrLdpMr7L^4OXvmO!ze^ne*Qz^VOq>3C^XUgMrV|E<`1<4OOs0)GLYd+Y*xAz!2 z;_wC%I^t!)5=j!ZAf5Fzg5VAPeJ`*us8K$Cv_gRg0U$VL{dl#~`Zy}IFsQNx@wBvv z^0Elw@6)4LWk%!BiinpZ=cu}=W#bcvg^Ea8LqsLL;unV*2&Bu)h<^oT4QRhV-Q{6|gTIrtzZ8v88CW`*-O651m!T5^(yxm@V~5I2m`-fh1?>LBY5l<0Z7tVSS~ z$~2aDqzT5sLp=;m561gQXM((qvdz0DA{OipppKX4r4#P!Ob3k9PEw`d3*g=1>qV^h z%l5If$4M7kCIxn$ML7b0e{Ww{Ph7z7*@KzLaT{nI1VGDSv0VA-utimf^|woK3KPXxM8I zUSV4CSwR*KvRA^+cR54{^Ffdvzdb6_A<~)AQ50J%geW8{ge9>v-7v8*W->l8)h8dY zEHYAj>B<+)>s77&5}`Sw#-zRc+gW)pomqY!z$ckppLT18Sw_1y09Z;(j-T~_hrmC3jIUpl9izCLy@#y+*LS})2kQt%o`3kE?3z~BYQHpMHnKsxK6aKSL$wvP z?cY`Eu^WtxH_hKh|I7p&)QoJG(DXSK?BWga4{<0@DpQyGD0fcEOfpT9WI42ywE|j` z+Ou5gS{~cj9qpW6!25qw&fTU@Iow(9Snt-)rs&g@x0IH|x8uVCkXOwr!uk{U(sx=1 zG%FLOg{46~N8Q;J3;M)Bvkn|HqUzSqyXFI;bwZ|ZN8F!TN` z{zv{#0Z0Ld0XG3JQZZ64fx1Dce@y;V1UdwRgB?E}eSDP8mOcn+3RCMd?fun99;z1{ z)tkvsfG^EZmE$4MEBFxy^d6O3np#_SBZ@nEH2TpF&q;l*C_`SQz^?E>#=W1@v}EWc z#*^rnq@DU!IqPedU=^Doldr-LgXbrNx8O(w`JMZI&j_M~hJ<^0rIDh)_kYPo;1l8L ztFza52AmW|l8!Z$N+&bC>@yY`POOJYSd=VPB$7VL$;&NgU!{|cRDHuP3N5miJ3Y8} zT5vLS>}<7Z=@i8_L_D|7kk$RPQV)!U%pZ1sF#M4G!RLpBFW=i23&;Ka@w>0_Uq$Le z>+tLD4&q$Yonz265x)lBerV&dm;e*%qA3O^=r^x5w!Rc_f?`Bz9ljfCbP}l*-mTt|xpKoJ~^nJ@$^BW6o zhf=cxvl9$J?JLJw$4729=Z0#r3sJ}3cb}JM8yEc?X2mOf^zL}>QpsSD`;l|rGOmtq zf=rNCv06!2+?VbzLTP@|6l2QuorJy){*;l|P}jHvO&8itoK493*1y_WM)@jEaNU@( z#Aw72QBVl0d6?#XSz2-De)(J2bZF+8s_(G0+m(BermQ?~Ws#@=Bx-Bj7D~&z(2D zCvi;K0cM8~4hL(BG1HhQG!4!B@KR)fB=7dN&`#%S=H&EmZP^8(^gSmpLH|WH1yo15uP!h! zgmnMBu=47(7cele+18pmKpmye!e$P3>?Y<8rr+2->>Q!cFfgJX!qBeWH=qfnhn=mx zi?D|n^}is5q5Xe`IjAZB1p>4Yqt;PUrIdDX{zl2i&dbh8EsjP>Nh#`VZXv7=koj+P z=sz)PDVFvh`}r?9zj;{yhm*a_f2##mkmH{j4lZ_1j(^jJ zqKf`ADy(Yl@y%8TU~Tu!-UZ4-oSTnR^k4A**UW!7{twKr|H0%E{68@N$IO3YigNs; z!2eL_zpVAIQK-4Z(L_1^ZF+GuM5S0i=xUH!15`AjE!;m#1ATO%FXMl=&^}hrES9Aw zwEZa$kka&kJ;_FOB=|x+^hUK99B11D5F&8H#m9`x%li||U3DoLu*Vu~0MJ@fMgPc} zs;bBMwV1b2an68{NqYg_bRLa>wL}#i6}Au{PXPOg7nQPR7`IB3Tbh0ZaMzjBk&|OV z?n3)vc6L)!%H41M{P_H(`}Ibo%@Xex|H)sx2!;t~1}zvDOqFf)izu!*snj0ILoN&Q8Gty0}UnF6E8F=Yp zV`FR6PM^(wr#{imEWODKVAr)cRc~1Nap)+RbT4?HbCMItK*D6qi-?eBI-+M=JDCc1 z5+05s&~CLnpsqJagy~B-Ep#%$#1&92?P=p6ci9lcfy`473Kv^>Ve0U+(_@+h_zO~>Hl&Fw(9Ie zRD6;TP}TBpUT~@je>RZuPHMt}143 z4m${#5JLWkK|Pm^x@$DM&DJow0ch-Q@&0dDotmphER8(+{j2)7`d-uTSmyz7)~ZL1;KF9>bTC)2(f zOnHbeO$}6|YhP?_E4itL+ArWd?4sCn+M&Q!d8?;x#h(uHu15qn3a;C6qd}d;;kTKe zU}1$RHex<^%fpJ%(+@uAXR`@$6ns1Iw;4Lx;6)<_D=8_VfjiKrq&QgKvi>7CcK5e9 zy{Jh&H>|d_EPRJa*LK2eNr8Dj#biv$fH;_Si+UzDVdR`<@ZP(pAKyXU74^uJAKp-l(Mulmt^WPvlT2i}v%8y|YTmDt?Y`$ZHp#HtVCGsrdrKyD+=dU{P?vg| z$b$)BAK`8UfcM-a>acJI3H)#d)EMcQIjGZs;lGG|b&GdD=Cf)uyul%Jcfg?ZUe-UnS^f}B-w?Sf zNC|BiIEM9`-EpsW`H&T!8I~b6Hb8jE<&P1j+lQ-(8Gj=bK~;;NV3>9#a-F=+-AR^M zW4>JK93O=3BoA}b-Zo$V`h5PR3u%3dO%ts%y2G1=0YMqi{-az+n#c>%&Qh2VJh$4i z?*@h&kAY2GW;s!LU|XhS>hStro}S9CKe1Q^78DRT5M|W>(L<)4Nmkm(8Cc-$(c1}v^ERmdzUI?N}*3b?ps`kRZ zUTy@^{&h$o`~xX)CK&IM))@SHmY7N|eHQuqD)7|%8Gb!d@|l0JD;n3#xf*Xwzw&LL zLQVy^0fx9xo*~xPz7F#j)winny|g`X4G}{9%|$!Fx9pzl^Mv2fX(dy6pQ;ZyFO#j% z7}g79 z-)T)*;MsR77>(%}bXmEgm+n6y%TJvwQEcv~nz#K|nRPU+Y=MTZ8et&n%($I;njwIq zUl-{2!dyIlLjg6DySBffhA0{SwSk9{pQ}VpWTQj0urn~zO+vlRJ-IhOO~YRo4SES;yAu>O#DJh(hCsOJj76%`O~7-d|GTra#vo6mrqof?w_Cl32R}2mJ3@2uf-E3lgR9%_ph7h% zGuriCP5}{rHoLHrzAM2HJnVz(JAnaWXtZ)!4DzM>QA1(AKDRn=&y3ue5D?Q@5q}wX z?D~cLD%f+3(>WY5x8sI)bcJD0A7|;((AVKFQPk++mdU=wZQT?~_vX*z@J$yF?@>vD z?^IJ65h2(X)Q@X!9WCk@*1%dXL@~+kkHm?rP){omp#X06beEb9^6Ug6p)IpMm`|_@ zEt#*=o|>)TDzNDx)RjaL^ov7V0 z>MDXBCsr@_yR+I1BKg8kvFeZe!Da4R$>@P8Cle&@iAS5>=4ggjd0H1Z8R@v-2DsLr4kw_ zMUmlz-IICj=-UNZXIyog787sU zV*zVJgnvLPknt*4yFAqCiJ6W1GL+6P|5ma7>DxVb>Y^&6n5?-Ldwlkb196}Q4T}07 zD_6TdImb--R;yArU_x!!9T3Rwgbe!KG;jaseCgWxCRsTaF*LuXfgo|UTdBnFWWP3? zv!pA0)Brt6TmG(&$il+)4>}Lj=ST-Okue5{aZ6K+Z*A)9h=56EGaN`Ki;5^zjnJc@ zXqIR?(`LCvavY7lgHs5p%|a+9Sx4#2^)J)C+^2|@I{c`c9T3hG!ggJ^J8~H@w@wJbe)TlQCmmzc4Jya(92=rY zE$Au##Yi%h;AVNv;`h+ujzS$Owyk>kmE^22R0DY9`bxX!{N$+$9J#Np7ME$btD+m zbfxq>Ow)(Z6CF<;_5WXj^6&g-JG?>7M&p! znMrUV4W*7<3szO{DT5L+h;JRFS!!s+V_Oq=a=33fe4q38Rpi9LX$9+{14721{39ZT zk=%tKR3$5TyR+SMoa|RkD;hLuyt0#s$zH{-vB|~Hd*ak{xE;cK+~9->O4ympjy8J5 zzofaIJMBJg4X69iIe&WAV%ne9UKq2^S5k%W@&+wW-`;w&%b#IjD6IdfOHU2Ze`pTO zOJv@W#HB{Up)pefKm}S7$3$v2Hhj#C1-PoVl10SU^GGi+ z<>3gmjF`rO^pgJmNW(!WnYM0;>wj;vTJKIayzkP9Xf5{dR3y@_$-P<_-Nn!xnpA0v%Y#hcdRtku34bLy?LME8oj(SqUefkNq60~YE zw)nISj(@2sy<$y^_r%WnfID3f$FrIXmgu=v=?;F8jqQE)@y^?f<}to{XBlc6OJHk+ zKnn`mvZj}`FoYeymv<1nU0F_SE;6ahgNB`-BA~F3I0@VfIJdwKgZnPc4D)w zPAb(P%h_B%k)3EGqEF=cp{9}PE8@i=O=?qExwm|Km-73bCQ z(m2gvSV7%{1<%K$3{WkGVPyY|(Z%{ePWtf^Mi!qffqwRt#``2I;QW`TNaKkpYfiiP zcJAQ~RT%9Xgm4jWL5xEZ;oS-a{yy2`H_T)m>(8Xx9=l`EBR&ve0X* zYw7qF2P2%VMJFczd26#rxlaBI?Z*lHSClSIKG3dfIf)?&`J z8`bW5JYp#>DgrREv3-dl6{2pDChfgpM?AQUoIMoE34o|?#O8u;Khha2qQVq0A>M8Q zIEIz4sXsb9*JARB14bwcJ+7^W4w7;-KQ~tP4n?j28Fr%C-HmxK(#mm0-PhedKM$qEC_wmO*IhcrdXn*)l z2ZRn~I!m^@meG&x_}`bUG6^jGQr25rzJ-%Mh|eHfDGM4!o;$9Rv~NIwm4{GTrOQ_^ zG+Bx$jj%w**C)P-7CA9{et4wvima1z-6v*x9hnwEZtMl3?j(2ltA+lrEs>HAGg5x2 zAQS3UU^zafe$_W~x8!iDXmB)=w+0J?Iv8+ER0qR68-Ga+C>G}ycs0{?M zWQM=1A#$VQc~b^qGFS%IFckdpqF$$GS-()OEgi7&C~jP+5CdZ(xt%MFrk={O|J zpSDfHyaD)DbYsM+n*RH%m$QF@hV7Byc$Di0&okVhc<_>Zg9$k28?$~F`z&YJ!Q_Zx zhK-!c6lkJR@g{XW5NDlN)}hhxt!69><{v+T39EyPUoZkmepAr?K(v|=Kim2QjCsUdu(?I1CXgQS@TEPUil-bMsk7%oLKticr;4 zvk;e6J87(0aYNiNoUMurjWU5r55eYe2+V{npUV^Sjcof)=R8{#q)hVb*vDG9vEr)0 zmS_Wp<`b$~dwWct55)oFYS-g5h$hmJ2^qhPNFf;Kt=aocFj%VxmX<8()*E}=E(=WG zw4zV?cS-qHd7E#(tW|u@j~xLj)h1+{WL5J7sPoh%e}s5-eZ*W} zGwo87ZRS3&_{lI)_1BW56id@`+*~Hb<~gramvV;=V=2c3$NT82MDw@z{-W1z~O& zMM;3XKIl@OJBEKl#0$`f^{0I-mq409&cI252DoCI@$2r1U(Y2%dm=7#;I;qe6x2oONG@$d)n?oL8{{8_%P zPWlL6yvg@pKcxb;WThaFk6xF2U+5nPClF@}-A88Ch}YFnz!Bol&-KBwDxW_S7=4b8 zio!S2;PVDP_I_vFE)s)MARpb1<6i#&D+c0d+p>iR9lDGFj+9k^kL<$YzSx+Un5HpZ zCbRs9p!VypY^!G>F!*HMY?+W8dG#ZfV_4Z&Ehscz zdI+N8r`N@pPv80rqM@Nd1TmFJTIUBRHXSKTyIp|+z_3DCw{tf5%n7V_^L$&*e61tO z+>gc^iw7fi5t7=n{=*YgqM{gHH*&#COIGi>Jzx?fo{Ip>u2?d5L?n|cJYsTy` zHR6yaylmKnM}#c4*TnlqY|t`3&$e|gs*1UVsfPH0ny`C=e$LECR0g}4$>49M_$G%# zF78kQPO_3*i1f5W3E+HEt=ONm7OUaQ8~0pX5oH|{2>WSQE%3&%%AOE2Z(%l2dKqaZ z>J-8OmL!L4DO;(A+3H6(B)Ow&AQIbh!hPfEbK!L=+>qCfI`P}E%QgRimynO=C^}sh zZ|r3M{Lp8e$mJjxWXVAd;j!ay($;AbA6_CxUGEby=?so+GKWs=p=^apZlQ>ut>}fb zgOVl%ohmu-M?Jk}!tLpq&+%)Sro~Isn}v)}nl58uciL9UCVk?Hx)!LG*JF?i`%f2! zx|-7OA4M0WxJr*_f7sN!+C;Whh>gS3c= zr`NyF#V$fFYlq=yvHq_!6PRdCQ!w&Af8wsPIWDiG6%{o;$U^7aJahM6ie0*dng~vq z?N`4?el1)}MS(^PGPsn2^OXVi_zX3-BmS9^u*?fzwP9a9)Mdv^KmpPqur}Vf!GrEj zICyY6v2w{$rKl-!!7C^_l}%pIVRAt$N?|tygZ7XauLyqt%ygj*U=iFk1b_$&fEvBO zvt{J9d)N5a-=h4u+wR-?Rbho+go|3*x7^I`KiofisU@Vf*KJ6(-Ir4p=WBYWrPFUv^8x}AYv06SYL$}W1kfb zs8+|udMk`3MBm6#iQW3d5I$!i`#ZeP6b~^ylXOT5L@pH_evV!mn~ofh>T?-pPDB+P z+7|OiTdyk={~qj>Yogs-278(Piol_J#VK0bL`=D@UGu#DlqV~K(`i+!e2@?^9dwZh zSX8FL+;Z#r%6MYx#R{%nr`LdY%iISyQMy2UTFq(+x~P!kV+5{HG47}E$YDVdmXhXc z){ZN~tMN+{s*Psp!)GGj11az0IS2#myJU*CqpJh|dhAkQw7Tm!e_%qoX_AxqS-XRA0Bx8kC_@iZJ^!LKRwl?8=@U9;ko4_W- zMNUqYOon9Mt$jxhh+LP5qL8C;DFs1_uiEDk7k&_^-Un0*vQp!jdS%1MUv}nC4eVa{ zZdjB(PvKE5o-^P}{Xy34dOM+WT#4~a?le9kN*la!AAi904-7?kiz%0g0m zRl}wB9eo+7vDWLDd8X(60pCqljr(#dRO2}F$@i0pq70WOqxv!*QLpwxfG7h^ z4(~_d-F=~pMWfswq9ch8<8P!#N4ED(Tf*C+J8^$eQ7;|W`QG3!XX=~5noUcjU{HDA z{Kt_l_|PwRUO&t#P`ttlQ>VHjquR_Abp!$ZH=!3`4>{3ehW8AT09y`P;V$l9P5Q&1 z7_j^;wBes?4$K^vRkmXHkB(?6XgO_(J%Jd$M|?#hZl2-3$ap@~UO@KJ@0TYeGj^)C zN>nFeK)9?VauEC^MLLG`(<~{)U3oUCjh8LoJ4f&qRaseWApSK?Z8nkN4^&7gZZWne zvyoZ{KH8Wg8~pX|V-|d&tP#wY3d=up?&)CR^RXMJ%K76aX~k+rscd9bhO?m?ne5cU z$uDKs4jGfps-+|3t6U8)A`P9yc0XleHdA4S`cc(XAk?7aCIbxU}uiMAIGiT)R@fDZ_i)KD7jNQ z_^$IUQiFDX_IWIfaiwng#I#$tA#TBW*(90NY2N zf?bZ}+PkU=U5|28Srz6@I=%5r?2e>w#&V_!qP!DrT~ zQsitxhIO#}FXDO_cTS5qp=+FoIX?IF;NOa0^pUq454(e6nO!5=CXjP|NY04=NbCPL z)JClIx|R3vE>&WP$b6X$0@V;T#j|G>c6WpADeYA`S4RZfXVnfyaa{5IK1r_D!!3y; z+EG+58#Z_dKi|GEvs;HJ+WRUPyjj5#8s6Fv!xwY-or}TdO=yM9oTaoY4i4>pWT#ho zFR`PB6~r39hX4Oi@CMQQ_V}(F-`yOko0p2W1Xc z((r);@fyyRsc8P@v9zv)aB|KvooEX8FjNYa{DnEiV))(7kegL9%%kb-8tGBL_9~VSc3ya+QT+Go| zJQr)?GvwZd=|1_QGfZc2Ny(byicgCoqOaQajY;%j2iepL<0UQ$48rGy2FuTfL+)un zivndQebob^GM6r81~+ovOq@RvO$Yg>i-o6Re5RQ>e?0zNwrRv%a%3>|%Cj0BK>BBk0V>iqGcpG=GGEBG6Rc#mN(S&( z5Y0SogArK2>d?HvrTquEdNrVF3bd^JC}$W3iz|kc^AmoV3BpA9LhcT8KgfQdEqbJAFf<8DcyYE45*IP`cQFzrEUdG;8y$8dZqH{D;YP% z$fKAftY-Lq?Y)fj-3HBI1Y(P%^c}svoy6^Fm?_v^tmZ=o@$q_FsLJLqyQiu}pWqR4 zQ7+33p`l$613svcrtF>NgPBnWynAN-7%LppVf zUYfBx7)a?es&K)R&AHS4YLJFzw&p$?D-Eyb-zXQpY+i*tVWv&Z^UFx?CBGWdU#IMj zDfy=ru1)9#%@|X=$zhs47OW}Va5R{JyQOl0%-E=9<#cw{-tOB!6FZ!^FC_RNl0cjn zl3c)f?Bc}v$AUZF@sSY&K)A@?3liHKEvw#q?}@Y6IB>Et{o^M#+D9|t0!w5{ukD?5 zP8U8ev5wK_eiFX`669%WDkL&0v{Z27N!j91SzxtK`DTpH-OwKsK!aMkHt;yHA{Ufx=nJZO(1K5(y5K5vbdy zM+y-}32zEcb3=^VW55BI7_jL56JU_Ego%yK{0%9$k)M8lS5;?3#8R^L_jMBbtZ!|R z8E2lrdHAyxjqf@?ZGgWh)E(hVQVMA(o7_5mLl_XHO7gvR%y|v3J}`^6`6^=x7M%yV z4>KvTGY*OD=r?`0a$blU`@Wz56du++2Wg8fYN}qasR;|skON*RNe_v`JK%xV|1cV5 zJ-r{tp+!o&9OZ8DVVu0i_yd>w=E*$%)hMK#8N?wCm-MoWn07VDuT;me#?AGv>H{CC z@f7h1r#Mkv^~FduAoqy=5vX%RfVXZC_Fcv4>>PwOUwD8Ex-m0uX0@}rr2$0?S1$s+ zE~pABgwKOTZ^32xpPDq1r7l>XhOWhHT2N|rPNNZ5&rN9O{r`CA&++n!xeuUC+57+Y z9@-3PeZ49BBz9l+w64)*5HUFfSi^N+Y4u?1WPp%fbyGih9Lk(z#e+S<^y|kC?_vZO zKxzO9!-n^RlLekkJxiWKoIXkXz=KGqh?8vO^RQyb59MAC^Z9k4pJ0@OdNi3#rFPCO2LJsBR(02gqX|I=V6BI7YAr=suGyfz!vCL(3l3U;!d0# z$xxn(I7tMvHD=pl0DwdG{iZ3{tgPV!-Wxg%f#HhURsEAju$kFe{-r+V8Fb@tm$e?z zaJ$}TZHelPAA;1BODf7Zu?Y32IVppyQ@A~VrMK3!>=rRmgdU>82O+WpQnC3=J28T7 zYN07N*&dx)aO+fyYgTvOvZTlFg@VL{OA*2^`^Vo4=-q(vbAmqIVO_JFFh0Tf3hh>Z z7sm=%l$gj-V^P0EUE}|j-$w1bZGJzb*!&U03kjPoD*7?}fFa5pr=qUK&WZu(U*99>s8;S9ic{g|bh$U^wETMVB@q=l^}JVL5qv0iH`xgh zv~NC@EO?Wmv*@+*I8&M3PB3I&x*pgxans6KVKPp4s^1JeOsz;g8eNzSLEb1CK?pJT z*q;wGrcGgb1R{=q{VeUn7W^4V8rUpj$zqF48sFi#YF1;ICP-{J;X|6>ZIu|i0N!m1 zeLwqzc1V5~`&NvXEh>U|^*0bKKndF2{Zs7oOao|K*BuYcrCShf6GmWO>s;Id8L;0d z^J6QB!c6}q(mxxSY3l#hG^|7VtMF$F3l`7)tyWg)>~T3k$2|A6&94MY&UKy6T0`Y` z()Y_)8)2eI<_>mRzn)RKZQahaZucCpYEE@YC3RVkzn@PW z+M?5Uce)GJ3c_9kDBmeN@PML-HY4?ILS-mo%iCh@!byMv8U!2W(%F3>iF2L}rB?p% zApk#3k~Z(Zoz0V(7s#Xwy9Kz2t5iByqhoRP$y$pyCX{QrLiqHCVseUeM~$aJi@85Y zH(0KR=`qRS>)+1B9jmf2(<#1tJ(f~CF_cHyG^u;|oauX1>!_r4Ili5gZk5g}Qu$E1 zSwJty(IHVzm76NwcHU|SJmSE~8KC_vkQ=0iy?pr>U))wjdi4No+l^_LTjqvS)mQU3 zoq@_SZ0MH9>3y#Eo1BtSna<}0;VjD2D-djOSeSK=0nH<(91!wXNeO~#CviUv#I&Lj zC(H`6VfsF$qf<8o1|gR57Hbc*2&!CVFZAtPE#=?T_wyY75Y)q;%?4U--ji<7Y+@HA zOcF3Z&#xLDVlw>7Z(BeF;j)*x4)lDmyul?O&V9v~?G6*pdKb)a+iwknb&UI=r+qwz zObiW68=?R15L3*|Dy?Ei-YbTJEsU5IKSG=KF^gQMyW4Bxa442MltR@~6To|Pg88F6 zv&#@p)vse1i-i0PrJ=LCCYv@R@6ig8x22Z~*6y6jwiGz;`@I_eJXZC>EixyK0^h`z z=pzFpdqY~kaY$V`e=0^Nt}BP@XoADy>8oE2#Nz5?MFiK>)C_iU3ozZR*tl@_@Zqyx zGiX+;06#8zfkw5>g1JlmsLYE%_d%sPj@Wr9Q|^jKHRvFi+I4;3--<2y^ae+ObW;_O zkE*>z=cG3hMX#M0S#~>`pFJ=r?T3mW0^H1mp{lgXJs~5Zc7yQd}t%V>fQF3vB}cl$>yavYXUmUdBJos{napG{(zU1E~pM zgICT`Bsl~)8e^uV*GAg^G5z2gPwSN`@v{yVFQ5lSCsOfNp#XVMZXZfo6fR}mTLCKE zvq=(W914O=5A-Xilsu&8{fK^qWU6dW-CGZpJcQlv{Wo}8xB%RH4p6xAM`iJ4ap>1~ O82L}ifU1wif&UAFceNY< literal 0 HcmV?d00001