{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": { "tags": [ "hide-cell" ] }, "outputs": [], "source": [ "from IPython.display import display, HTML, Markdown\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import seaborn as sns\n", "from myst_nb import glue\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Why do we have to learn statistics?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```{note}\n", "This chapter is adapted from Danielle Navarro's excellent [Learning Statistics with R](https://learningstatisticswithr.com) book {cite}`Navarro2011`. The main text has mainly be left intact with a few modifications, also the code adapted to use python and jupyter.\n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Introduction" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "> \"Thou shalt not answer questionnaires \n", "Or quizzes upon World Affairs, \n", "Nor with compliance \n", "Take any test. \n", "Thou shalt not sit \n", "With statisticians nor commit\" \n", "-W.H. Auden [^quote]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "[^quote]: The quote comes from Auden's 1946 poem Under Which Lyre: A Reactionary Tract for the Times, delivered as part of a commencement address at Harvard University. The history of the poem is kind of interesting: http://harvardmagazine.com/2007/11/a-poets-warning.html" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "## On the psychology of statistics\n", "\n", "To the surprise of many students, statistics is a fairly significant part of a psychological education. To the surprise of no-one, statistics is very rarely the *favourite* part of one's psychological education. After all, if you really loved the idea of doing statistics, you'd probably be enrolled in a statistics class right now, not a psychology class. So, not surprisingly, there's a pretty large proportion of the student base that isn't happy about the fact that psychology has so much statistics in it. In view of this, I thought that the right place to start might be to answer some of the more common questions that people have about stats...\n", "\n", "A big part of this issue at hand relates to the very idea of statistics. What is it? What's it there for? And why are scientists so bloody obsessed with it? These are all good questions, when you think about it. So let's start with the last one. As a group, scientists seem to be bizarrely fixated on running statistical tests on everything. In fact, we use statistics so often that we sometimes forget to explain to people why we do. It's a kind of article of faith among scientists -- and especially social scientists -- that your findings can't be trusted until you've done some stats. Undergraduate students might be forgiven for thinking that we're all completely mad, because no-one takes the time to answer one very simple question:\n", "\n", ">*Why do you do statistics? Why don't scientists just use **common sense?***\n", "\n", "It's a naive question in some ways, but most good questions are. There's a lot of good answers to it (including the suggestion that common sense is in short supply among scientists), but for my money, the best answer is a really simple one: we don't trust ourselves enough. We worry that we're human, and susceptible to all of the biases, temptations and frailties that humans suffer from. Much of statistics is basically a safeguard. Using \"common sense\" to evaluate evidence means trusting gut instincts, relying on verbal arguments and on using the raw power of human reason to come up with the right answer. Most scientists don't think this approach is likely to work.\n", "\n", "In fact, come to think of it, this sounds a lot like a psychological question to me, and since I do work in a psychology department, it seems like a good idea to dig a little deeper here. Is it really plausible to think that this \"common sense\" approach is very trustworthy? Verbal arguments have to be constructed in language, and all languages have biases -- some things are harder to say than others, and not necessarily because they're false (e.g., quantum electrodynamics is a good theory, but hard to explain in words). The instincts of our \"gut\" aren't designed to solve scientific problems, they're designed to handle day to day inferences -- and given that biological evolution is slower than cultural change, we should say that they're designed to solve the day to day problems for a *different world* than the one we live in. Most fundamentally, reasoning sensibly requires people to engage in \"induction\", making wise guesses and going beyond the immediate evidence of the senses to make generalisations about the world. If you think that you can do that without being influenced by various distractors, well, I have a bridge in Brooklyn I'd like to sell you. Heck, as the next section shows, we can't even solve \"deductive\" problems (ones where no guessing is required) without being influenced by our pre-existing biases." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### The curse of belief bias\n", "\n", "People are mostly pretty smart. We're certainly smarter than the other species that we share the planet with (though many people might disagree). Our minds are quite amazing things, and we seem to be capable of the most incredible feats of thought and reason. That doesn't make us perfect though. And among the many things that psychologists have shown over the years is that we really do find it hard to be neutral, to evaluate evidence impartially and without being swayed by pre-existing biases. A good example of this is the ***belief bias effect*** in logical reasoning: if you ask people to decide whether a particular argument is logically valid (i.e., conclusion would be true if the premises were true), we tend to be influenced by the believability of the conclusion, even when we shouldn't. For instance, here's a valid argument where the conclusion is believable:\n", "\n", ">No cigarettes are inexpensive (Premise 1) \n", ">Some addictive things are inexpensive (Premise 2) \n", ">Therefore, some addictive things are not cigarettes (Conclusion). \n", "\n", "And here's a valid argument where the conclusion is not believable:\n", "\n", ">No addictive things are inexpensive (Premise 1) \n", ">Some cigarettes are inexpensive (Premise 2) \n", ">Therefore, some cigarettes are not addictive (Conclusion)\n", "\n", "The logical *structure* of argument #2 is identical to the structure of argument #1, and they're both valid. However, in the second argument, there are good reasons to think that premise 1 is incorrect, and as a result it's probably the case that the conclusion is also incorrect. But that's entirely irrelevant to the topic at hand: an argument is deductively valid if the conclusion is a logical consequence of the premises. That is, a valid argument doesn't have to involve true statements.\n", "\n", "On the other hand, here's an invalid argument that has a believable conclusion:\n", "\n", ">No addictive things are inexpensive (Premise 1) \n", ">Some cigarettes are inexpensive (Premise 2) \n", ">Therefore, some addictive things are not cigarettes (Conclusion)\n", "\n", "And finally, an invalid argument with an unbelievable conclusion:\n", "\n", ">No cigarettes are inexpensive (Premise 1) \n", ">Some addictive things are inexpensive (Premise 2) \n", ">Therefore, some cigarettes are not addictive (Conclusion)\n", "\n", "Now, suppose that people really are perfectly able to set aside their pre-existing biases about what is true and what isn't, and purely evaluate an argument on its logical merits. We'd expect 100% of people to say that the valid arguments are valid, and 0% of people to say that the invalid arguments are valid. So if you ran an experiment looking at this, you'd expect to see data like this:\n", "\n", "\n", "| | conclusion feels true| conclusion feels false |\n", "|------------------ |:--------------------:|:----------------------:|\n", "|argument is valid |100% say \"valid\" |100% say \"valid\" |\n", "|argument is invalid|0% say \"valid\" |0% say \"valid\" |\n", "\n", "If the psychological data looked like this (or even a good approximation to this), we might feel safe in just trusting our gut instincts. That is, it'd be perfectly okay just to let scientists evaluate data based on their common sense, and not bother with all this murky statistics stuff. However, you guys have taken psych classes, and by now you probably know where this is going...\n", "\n", "In a classic study, Evans, Barston, & Pollard (1983) {cite}`Evans1983` ran an experiment looking at exactly this. What they found is that when pre-existing biases (i.e., beliefs) were in agreement with the structure of the data, everything went the way you'd hope: \n", "\n", "\n", "| | conclusion feels true| conclusion feels false |\n", "|------------------ |:--------------------:|:----------------------:|\n", "|argument is valid |92% say \"valid\" | |\n", "|argument is invalid| |8% say \"valid\" |\n", "\n", "Not perfect, but that's pretty good. But look what happens when our intuitive feelings about the truth of the conclusion run against the logical structure of the argument:\n", "\n", "| | conclusion feels true| conclusion feels false |\n", "|------------------ |:--------------------:|:----------------------:|\n", "|argument is valid |92% say \"valid\" |**46% say \"valid\"** |\n", "|argument is invalid|**92% say \"valid\"** |8% say \"valid\" |\n", "\n", "Oh dear, that's not as good. Apparently, when people are presented with a strong argument that contradicts our pre-existing beliefs, we find it pretty hard to even perceive it to be a strong argument (people only did so 46% of the time). Even worse, when people are presented with a weak argument that agrees with our pre-existing biases, almost no-one can see that the argument is weak (people got that one wrong 92% of the time!)\n", "\n", "If you think about it, it's not as if these data are horribly damning. Overall, people did do better than chance at compensating for their prior biases, since about 60% of people's judgements were correct (you'd expect 50% by chance). Even so, if you were a professional \"evaluator of evidence\", and someone came along and offered you a magic tool that improves your chances of making the right decision from 60% to (say) 95%, you'd probably jump at it, right? Of course you would. Thankfully, we actually do have a tool that can do this. But it's not magic, it's statistics. So that's reason #1 why scientists love statistics. It's just *too easy* for us to \"believe what we want to believe\"; so if we want to \"believe in the data\" instead, we're going to need a bit of help to keep our personal biases under control. That's what statistics does: it helps keep us honest." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## The cautionary tale of Simpson's paradox\n", "\n", "The following is a true story (I think...). In 1973, the University of California, Berkeley had some worries about the admissions of students into their postgraduate courses. Specifically, the thing that caused the problem was that the gender breakdown of their admissions, which looked like this...\n", "\n", "| | Number of applicants | Percent admitted |\n", "|-------|:--------------------:|:----------------:|\n", "|Males |8442 |46% |\n", "|Females|4321 |35% |\n", "\n", "...and the were worried about being sued. Given that there were nearly 13,000 applicants, a difference of 9% in admission rates between males and females is just way too big to be a coincidence. Pretty compelling data, right? And if I were to say to you that these data *actually* reflect a weak bias in favour of women (sort of!), you'd probably think that I was either crazy or sexist. \n", "\n", "Oddly, it's actually sort of true ...when people started looking more carefully at the admissions data Bickel, Hammel, & O’Connell (1975) {cite}`Bickel1975` they told a rather different story. Specifically, when they looked at it on a department by department basis, it turned out that most of the departments actually had a slightly *higher* success rate for female applicants than for male applicants. {numref}`tbl:simpsontable` shows the admission figures for the six largest departments (with the names of the departments removed for privacy reasons):" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "tags": [ "hide-input" ] }, "outputs": [ { "data": { "application/papermill.record/text/html": "
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
Male ApplicantsMale Percent AdmittedFemale ApplicantsFemale Percent admitted
Department
A82562%10882%
B56063%2568%
C32537%59334%
D41733%37535%
E19128%39324%
F2726%3417%
\n
", "application/papermill.record/text/plain": " Male Applicants Male Percent Admitted Female Applicants \\\nDepartment \nA 825 62% 108 \nB 560 63% 25 \nC 325 37% 593 \nD 417 33% 375 \nE 191 28% 393 \nF 272 6% 341 \n\n Female Percent admitted \nDepartment \nA 82% \nB 68% \nC 34% \nD 35% \nE 24% \nF 7% " }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "simpsontable" } }, "output_type": "display_data" } ], "source": [ "simpson = {\"Department\": [\"A\", \"B\", \"C\", \"D\", \"E\", \"F\"],\n", " \"Male Applicants\": [825,560,325,417,191,272],\n", " \"Male Percent Admitted\": [\"62%\",\"63%\",\"37%\",\"33%\",\"28%\",\"6%\"],\n", " \"Female Applicants\": [108,25,593,375,393,341],\n", " \"Female Percent admitted\": [\"82%\",\"68%\",\"34%\",\"35%\",\"24%\",\"7%\"]\n", " }\n", "simpson_df = pd.DataFrame(simpson)\n", "\n", "glue(\"simpsontable\",simpson_df.set_index(\"Department\"), display=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "```{glue:figure} simpsontable\n", ":figwidth: 600px\n", ":name: \"tbl:simpsontable\"\n", "\n", "Admissions figures from the six largest departments in Bickel, Hammel, & O'Connell (1975) {cite}`Bickel1975`.\n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Remarkably, most departments had a *higher* rate of admissions for females than for males! Yet the overall rate of admission across the university for females was *lower* than for males. How can this be? How can both of these statements be true at the same time?\n", "\n", "Here's what's going on. Firstly, notice that the departments are *not* equal to one another in terms of their admission percentages: some departments (e.g., engineering, chemistry) tended to admit a high percentage of the qualified applicants, whereas others (e.g., English) tended to reject most of the candidates, even if they were high quality. So, among the six departments shown above, notice that department A is the most generous, followed by B, C, D, E and F in that order. Next, notice that males and females tended to apply to different departments. If we rank the departments in terms of the total number of male applicants, we get **A**>**B**>D>C>F>E (the \"easy\" departments are in bold). On the whole, males tended to apply to the departments that had high admission rates. Now compare this to how the female applicants distributed themselves. Ranking the departments in terms of the total number of female applicants produces a quite different ordering C>E>D>F>**A**>**B**. In other words, what these data seem to be suggesting is that the female applicants tended to apply to \"harder\" departments. And in fact, if we look at all {numref}`fig:berkeley` we see that this trend is systematic, and quite striking. This effect is known as Simpson's paradox. It's not common, but it does happen in real life, and most people are very surprised by it when they first encounter it, and many people refuse to even believe that it's real. It is very real. And while there are lots of very subtle statistical lessons buried in there, I want to use it to make a much more important point ...doing research is hard, and there are *lots* of subtle, counterintuitive traps lying in wait for the unwary. That's reason #2 why scientists love statistics, and why we teach research methods. Because science is hard, and the truth is sometimes cunningly hidden in the nooks and crannies of complicated data." ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "tags": [ "hide_cell" ] }, "outputs": [ { "data": { "application/papermill.record/image/png": "iVBORw0KGgoAAAANSUhEUgAAAcMAAAFgCAYAAADHDhKeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xd4XNW18OHfmqreZbn3jm1sPGDTi4EAoSb0EgiETkJCCEluGilfQtolhIQAF0gIxcS0QCBU08EYy71btnGRm3qXpu7vjzMysq0yI82MNNJ6n0ePNGdO2XI5a/Y+e68lxhiUUkqpgczW2w1QSimlepsGQ6WUUgOeBkOllFIDngZDpZRSA54GQ6WUUgOeBkOllFIDngZDpZRSA54GQ6WUUgOeBkOllFIDnqO3GxCJM844w7z++uu93QyllOqM9HYDVPclRc+woqKit5uglFKqH0uKYKiUUkrFkwZDpZRSA54GQ6WUUgOeBkOllFIDXlLMJo1EMBiissnHlrJGPiwpp8EboDDDzemHDaYgw0V+hru3m6iUUqqP6hfBsKrRx4vLS3nw/a2U13sPeO+Pb21iypBM/uesKcwckUNmirOXWqmUUqqvSvph0qpGL997biW/fGX9IYGw1fo99Vz16Gf8a8lO6pr9CW6hUkqpvi6pg2Fds58/vLGJhevLItr/V6+u59OtlRhj4twypZRSySSpg2G9N8D8JTuiOuaXr66josEXpxYppZRKRkkbDP3BEM8W7yTaTt7OqmZ21zTHp1FKKaWSUtJOoKlp8vP8stJuHTv/sx1MH5aNzdZ/UglWNnip9waoqPeSk+okO81JQYYbkf7zOyqlVLwkbTA0xlDT2L3JMJWNPryBEKkue4xblXj+YIjPKxq5Y8EK1uyq2799ZF4av/nKdGaOyCHdnbR/zUoplRBJO0wK4HZ2r/luhw2HvX/0mHbXNHPBXz8+IBAC7Khq4spHF7NmV20vtUwppZJH0gbDFKedw4Zld+vYuWPzcdqT9lffr8Hr5w9vbqTRF2z3fWPgR/9eQ0VD+0tOlFJKWZI2ImSlOvnmyeOjPs5lt3Ha1KI4tCjxGr1BXlu9t9N9Npc10OQNJKhFSimVnJI2GAKMyk9neG5qVMecO3MomSn94xlaIGQIhLqeTtvkb7/nqJRSyhLXYCgi3xGRtSKyRkTmi0iKiIwRkcUisllE/iUiru6evyDDxWNXH0mqM7KJMGML0rnrS5NIc/WPYOi0CTlpnaeXE4HsVE1Bp5RSnYlbMBSRYcC3AI8xZhpgBy4Ffgvca4wZD1QD1/XgGowpTOPFW4+hsItE3LNG5jD/hrkMykrp7uX6nJw0F187elSn+5wwoYD0fhL8lVIqXuI9TOoAUkXEAaQBe4BTgOfC7z8OnN+TCzjtdiYOyuS/tx/HP75+JLNG5NC6tM5pF86aPphXv3Ucj3zNQ1E/CoQALoeNq+aOZvLgzHbfL8hw8avzp5OlPUOllOqUxDNPp4jcDvw/oBl4E7gd+DTcK0RERgCvhXuOBx97A3ADwMiRI2dv3749omtWNfoIBEOEAJtAmtNORj+vVFFe7+X5ZaU8+uHnlDd4SXfZueyokXzj+LEUZenCe6USRP+jJbG4BUMRyQWeBy4BaoBnsXqEd0cSDNvyeDymuLg4Lu3sLwLBENVNfoLGYBPITnHijvBZqlIqJjQYJrF4Pkw6FfjcGFMOICIvAMcCOSLiMMYEgOHArji2YcBw2G0UZmoBY6WU6o54PjPcAcwVkTSxxunmAeuAd4ELw/tcDbwUxzYopZRSXYpbz9AYs1hEngOWAQFgOfAw8CrwjIj8Krzt0Thd33p+GDL4AiHAmnBitwm5aS7s/ShJt1JKqZ6J6wSaWInmmWEgGKKq0cfCDWU8+el2Nu6t378w3W4TJgzK4GLPCM6aPoS8dBcuR1LnHVBK9R36CTuJ9atgWNPkY/HWKu58diX1XaQgS3Xa+dX50zhl8iBy07u97l8ppVppMExi/aZbVN3k4963N3Hjk0u7DIQAzf4g3312JXf/Zy1Vjb4EtFAppVRf1S+CYaM3wJOLtvP4J5GtRWzrpRW7uW9hCXXN3auNqJRSKvn1i2BYWt3EvW9v6vbxj3+yjc1lDTFskVJKqWSS9MGwpsnHT19aSwTFGzr1o3+vplLr/iml1ICU9MGwwRtg8edVPT7P+j31VDfpUKlSSg1ESR8MV5fWxuxcn26tjNm5lFJKJY+kD4aLYhjAPt1aiVcL4Sql1ICT1MEwEAxRE8OhzYaWAP6ePnxUSimVdJI6GNptQlZq7DLKpbrsOO26blYppQaapA6GIsKcMfkxO99RY/JwO7TskVJKDTRJHQwBDh+RE7NzHTMudoFVKaVU8kj6YJjhdnDEyJ4HxHGFGeSnaz1ApZQaiJI+GOalu/j5edN6fJ5fnncYBVocVymlBqSkD4YAY/LTuf74Md0+/qtHDOOwoVkxbJFSSqlk0i+CYUaKg5tPGs/5M4dGfezJkwr5n7OmkJ2mZZyUUmqg6hfBEKzh0p+dcxi//ep03BEU7HXahZ+cPYU/Xnw4+Rk6PKqUUgNZ7Bbp9QG56S4umDWcEyYW8sKyXcz/bAel1c0H7DMkO4ULZw/n0iNHkpfhJNXZr/4IlFJKdUO/qnTfViAUorrRRyBk8PpDGCDFacNhE7JTXbgi6D0qpVQUNGNHEuu33SKHzUZhZkpvN0MppVQS0O6RUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAS9uwVBEJonIijZfdSLybRHJE5G3RKQk/D03Xm1QSimlIhG3YGiM2WiMmWmMmQnMBpqAF4EfAAuNMROAheHXSSkQChEMhXq7GUoppXooURlo5gFbjDHbReQ84KTw9seB94DvJ6gdMVFe72XNrlpeXrkbgPNmDuWwoVma8UYppZJUooLhpcD88M9Fxpg94Z/3AkXtHSAiNwA3AIwcOTLuDYzUtopGLn34U/bWtezf9uLyXQzNTuGZG+YyMj+9F1unlFKqO+I+gUZEXMC5wLMHv2esLOHtZgo3xjxsjPEYYzyFhYVxbmVkyupbuOz/DgyErXbXtnD5I4spr/f2QsuUUkr1RCJmk54JLDPG7Au/3iciQwDC38sS0IaYWL+7jj21hwbCVqXVzZTsq09gi5RSSsVCIoLhZXwxRArwMnB1+OergZcS0IYeC4UM/1m1p8v9Xlm1OwGtUUopFUtxDYYikg6cBrzQZvM9wGkiUgKcGn6dFGwRVCsT0ZJmSimVbOI6gcYY0wjkH7StEmt2aVKx2YTzZw5jQXFpp/udN3NYglqkVOzVtfhpaAmwr64FERiclUqG205GirO3m6ZUXPXb4r7xMKEog5F5aeyoajpgu8MmnDx5ECdPGsTw3FR2VDXh9QcJhgxOu41Ulx2HXUh3OUh36x+56nsCwRCl1c388pV1vLuxjFB4WpvDJnzpsMH88KzJDM1OxRbJ8IhSSUisCZ19m8fjMcXFxb3dDABKq5q48tHFbKtsoiDDxRVzRnHypELe21TOB5vKWbenjhb/oQvxCzPdzBiezTXHjGbS4Ezy0l04bJoNT/U+YwwlZQ185YFPaPAG2t0nN83Jv289llG6dKgz+kkhiWkw7Iby+haqGn0EQob739nMW+v2EQxF/udYkOHihhPG8pUjhlOQ4Y5jS5XqWnm9l4se/IRtlU2d7jdjeDZ/v+ZI8vXfbEc0GCYx7Zp0g80mLNlWxVce+ITX1+yNKhACVDT4+PV/N3DJQ4vYsLeOxg4+jSuVCHtqm7sMhACrSmupafInoEVKJZ4GwyiV13u57enl/Pjfa/EGepaXdEt5I2fd9yGvrt5DXbPeZFTveGVl10uGWr23qTyOLVGq9+hsjihUNHj5xj+XsHJnbczOGTJw13Or8PpDnD9rKJk6a08lWJM/8pGJJh3FUP3UgAmGvkCQ2uYAIWMIGYMANhFSnHayUrsOQNWNPu58dmVMA2FbP3lpDSPyUjl+fAF2u3bYVeIcNjQ74n0nDc6MY0uU6j39OhgGgiGqmnxs3FPPc8tKWVVay7bKRlrnDOWkOZk2NJuTJhVy1vQhZKY42u2Z+YNBPiwp572N8R0i+tYzy3nz2ycwODs1rtdRqq2TJw3CaRf8wc6ffWe4HcwckZOgVimVWP12Nmllg5dXVu3hzwtLqGz07d+e6rSTlerAGwhR2+zfHxjtNuH0qUX8+MtTKMx043LY9x+zu6aZU//3fZp8wZj8Pp0547DB3PPV6eSkueJ+LaUAGr0B/vJOCX97f2un+/3oy1P42txRuJ32TvcbwHQ2aRLrdz3DQDBEaU0ztz61jLW76wArjdopkwdx+ZxRuOw2Khu9pDjtDMp088bavSwoLqWq0cdra/by/qZyfvOV6Zw4sZCcNBdef5BHP/o8IYEQ4PW1e7nrjEkaDFXCpLsdXH/COBq8QZ74dHu7+3zzlPFcOHu4BkLVb/WrnmEgGGLjvnoueejT/YuH89Nd3H/5LIq3VTP/sx0HVJ1w2W2cMW0w1xwzmr++t5mF678ooHHjiWO5+cRx+AIhTv7DezQmKBgCXH7UCH5y9mGkuvTGoxKnpslHdZOff3z8OUt3VCMIx4zL58q5o8hOdUb0bH2A055hEutXPcMdVU1c/OCi/YEr0+3gwatm8+tX17N8Z80h+/uCIV5euZu31u3jgSuOIBg0+6eOP/T+Vsbmp5OZ6kxoIAR4ftkubp83UYOhSqicNBc5aS5+eNYUmnwBQMhwO3A5dEKX6v/6zb/yygYvNz6x9IDAdcfpE3nw/S3tBsK2mv1Bbnt6Gd8+dQLuNv/x1+yu4821e+PW5o54A9bEH6V6Q4rTTl66m7x0lwZCNWD0i3/p3kCQBcWllJQ17N+W5rIzfVj2AUOfnWn0BfnPqj2cPWPI/m3Th2ezsjQ+Sym6snxH5wFcKaVU7PSLYFjb5Oev724+YNuXZwzh5ZXRFdp9dulOLpg1fP/rUXlpbKtsjEkbo7VoawW+QOfDsw0tfsrqWiira6E5wUO5SinVn/SLZ4ZLtlUfkm1/XEEGb0Q5xFnXHMBp/+IZuAF6a35RbZMfXyB0wBKPVg0tfnZWN3P/OyUs3lqFzSacNqWIG08cS1FWCik6408ppaKS9MGwvsXPv5bsOGS7wy4EokygfbBETw3LTXMybVg2g7NSGD8oo93acfUtfl5dtYcfvLD6gO1Pf7aDZ5fu5B9fPwrPqFydAq+UUlFI+mHSZn+Q1bsOfa5XXu9lWG50mVzsNsHeNgAlIBo6bMKFs4fzzA1z+d2FM5g10srwsbm8gQVLdrJmVy1l9S0EglZS8KpGHz98cXW75/IHDd94vJgaTfqtlFJRSfqeoT9oqG6nrMwrq/Zw97mH8fqayIdKT50yiA9LKva/rmr0UZDhoqIhPjM7pwzJ5O5zD+Od9WXc8EQxdc0HDvU+G/5emOHmlpPHcd7Mofxz0fZOh26b/VbquAtnj4hLm5VSqj9K+p5hsIN8irtqmjHGMCo/LeJzXTFnFPM/+2LIdXVpLdOGRZ7EOBonTSzkp2cfxncXrOShD7YeEgjbKm/w8vP/rOONtftYtr26y3N/vLmSUA+HiJVSaiBJ+mDY2VDmvW9v4vcXHk6mu+sO8C0njWNVaS1l9d7921bvqmVWHBITzxiezfUnjOW6x5dQWt0c8XFl9V5SIliIn5niaPd5o1JKqfYlfTB0O2x0dN9fv6eee9/exGNfP7LDbPt56S5+ds5UBmen8Me3Nh7w3uKtVZwwsRCJYVxxO2z89Oyp3LFgRdT5Tj/YVM5Z04Z0ud/FHh0iVUqpaCT9M0O7TRhTkMGW8oZ231+0pZLv167iG8eP5YdnTubNdfsor7cSdc8dm8egzBTmf7aDV1cfWu3bFwyxaEslJ04ojFmF78vnjOTfy3exr87b9c4HWbq9mru+NInhuakd9ihnj8plSHZKT5uplFIDStL3DLNTnMwZm9fh+0OyU8hJc/Lc0lJ+9vIa9tQ2k5niIBAK8dD7W7ny0cXtBsJWT3+2g+uOHxOz9p49YyjPLSvt9vH/8+Ia7r9sFocNzTrkvRMmFvDglUeQn+HuSROVUmrASfqeodNh45qjR/P04i8mvhw2NIurjh7FhEGZ7K5pprS6CWNgaE4qw3NT2VndzBOLtrNxX32X5y+tbqbZF+RrR4/in4vaL28TqXGFGWwua6DFH+r2ObaUN3D7Myu447SJTCjKYOn2apw2G0ePzyc71Umuln5SSqmoJX0wBMjPcDF1SBa7a5v52TmHYYzh7x9va3f9IcCkokyuPmYUN5wwlp+8tIby+s6HLMcPymDWyFzeWrfvgBJQ0Zo+LJtVpT3PObqjqolv/2sFv/nKNL529Ogen08ppQa6pB8mBcjPcPPny2by92uO5NnindyxYGWHgRBg4756/ufFNTz8wRYevmo20ztZPnHFUSMpyHBTmOnmyevmkJXS/c8PQ3NS2FnV1O3jD6bJvJVSKjbiGgxFJEdEnhORDSKyXkSOFpE8EXlLRErC33Njca1Up51vzl/OJ1sqIz5m2Y4abvjnUn52zlTGFWYc8n5RlpvbT5uwv6jp6Pw0nrvpGHLSulfkVESI5eq/kIFkKM6slFJ9Xbx7hvcBrxtjJgOHA+uBHwALjTETgIXh1z1S2eDl+n8ujWrNXqvyBi93LFjJry+YdkAqtky3g6e+MYfCNpNR7HYb4wdl8Oo3j2fOmI4n7XSkrK6FIdnRpYjrzMSiTCSW6z6UUmqAilswFJFs4ATgUQBjjM8YUwOcBzwe3u1x4PyeXMcXCPHvFbtZt6eu2+fYUdXEG2v3csWckYCV/uyFW45hTH76IcHGZhOG5abytyuP4I8XHR7VsOma3bWdDslG68jRMelUK6XUgBfPCTRjgHLg7yJyOLAUuB0oMsa0rmXYCxT15CI1TT7+vLCkRw0FeGrxDp6+fi6N3gDfP2MyBRnuTrO45KW7OefwIRw/oYDlO6v588LNrN3deUDeVtGIZ3QuDlvPK2pkpToYlhO7XqZSSg1k8QyGDuAI4JvGmMUich8HDYkaY4yItBsVROQG4AaAkSNHdniRkrIGamNQpcEbCLGzqomfnzuNjAh7ey6HnUFZdr502BA8o/LwBUNUN/pYsbOWXTVN+AIh0t0Opg7JYkJRJmkuOy67jXNnDuWFZbt61N6vHzNGl1EopVSMxDMYlgKlxpjF4dfPYQXDfSIyxBizR0SGAGXtHWyMeRh4GMDj8bQbMEMhw3sb2z28Wz4oKees6V2nO2tP60L3IdmpTB3a+VDoXV+axML1Zd0O4sNzU7ly7iicjn4xGVgppXpd3O6mxpi9wE4RmRTeNA9YB7wMXB3edjXwUnev0eALdDk0GY2Ne+upb4l/LcDCDDePXO3B0Y1k2mkuO49dfSSFmZplRimlYiXeXYtvAk+JyCpgJvBr4B7gNBEpAU4Nv+6WUMjgC3Q/m8vBvIFQTJc+dMRutzFtaBbzb5hLdmrkyzSKsty8cPMxjClMj2PrlFJq4IlrBhpjzArA085b82JxfpfD1u01f+3JSXV2q7fWHakuB0eMyOGNbx/P717fyMsrd3c4qcbtsHH5USO5+eRxFGa4dTmFUkrFWFKnY0tzOZgzNp+318fmueHsUblR9dR6ym63MTg7lV+eP427zpjEsh01fFRSwY6qJmwCYwszOG5CAdOHZZOV4iQ1glqGSimlopfUwRDgxImF/L9X18fkXF+aNrhXel3pbgfpbgdnTU/l9KlFNPuDCFaw1yK9SikVf0k/HTEv3cXkwZk9Ps+wnFRG5Pb+uj2H3UZmipOMFKcGQqWUSpCkD4YFGW5++9UZPT7PPV+dTn66ztBUSqmBKOmDIcDYwnRuPGFst4+/aPZwZgzP1p6YUkoNUP0iGGamOLnpxHFcflTHmWo6ctb0wfzwrMlkp2o2F6WUGqj6RTAEyE13cdcZk3joytkRJc9Oddr540WH86vzp5Onw6NKKTWgSTLUw/N4PKa4uDiifQPBEFWNPt7fVM5Ti3ewfk8d3vDCfKddmDQ4k0uPHMFpUweTm+bCpSnNlFKxoc9ZkljSL604mMNuY1BWChfOHs68KUX4gyH8wdZgaMNhE3LTXPp8UCml1H79Lhi2EhHy0vU5oFJKqa7pGKFSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg2GSqmr0sbummR2VjeyoamJvbTO1Tb7ebpZSqp8SkfdExNPb7TiYiJwkIq/09Dz9trhvf+QPBqlu9LN6Vy33LSxhVWnt/vdE4LjxBXzrlAmMLkgnP92FzSa92FqllLKIiMMYE+jtdnRGg2Enapv9+IMhbCJkpThw2HuvI93oDbB8Rw23zV9GTZP/kPeNgQ9LKviwpIKh2Sk8es2RjCtMx+Ww90JrlVK9RURGA68BHwHHALuA88Lb7jTGFItIAVBsjBktItcA5wPpwATgD4ALuArwAmcZY6rCp79KRB7Bih3XGmM+E5F04H5gGuAE7jbGvBQ+71eADMAOnHhQO/8GHAmkAs8ZY34W3r4NWACcCTQDlxtjNovIP4AWwANkAXcYY15pcz4bsBE4xhhTHn69CTjaGFPe1Z9bl8FQRH5rjPl+V9s6OHYbUA8EgYAxxiMiecC/gNHANuBiY0x1V+dKlNbe14qdNfzjk21UNHhJddqZN6WIC2cPJzPFQWaKM6FtavIFWLihjNufWY4xXe9vswkflpST7rbjstux2yAvzYW9F4O5UiqhJgCXGWOuF5EFwFe72H8aMAtIATYD3zfGzBKRe4GvAX8K75dmjJkpIicAj4WP+xHwjjHmWhHJAT4TkbfD+x8BzGgTTNv6kTGmSkTswEIRmWGMWRV+r9YYM11EWq99dnj7aOAoYBzwroiMbz2ZMSYkIk8CV4SPORVYGUkghMh6hqcBBwe+M9vZ1pGTjTEVbV7/AFhojLlHRH4Qfh3pueKqyRtg+c4abn5qKXXNB/boV5bW8qe3N3HNMaO57ZQJ5KW7EtauDXvr+XYEgdBlt/GL8w8jP93Fk5/u4J7XNhAyMCjTzXXHjeH8WcMozHDr8KlS/d/nxpgV4Z+XYgWRzrxrjKkH6kWkFvhPePtqYEab/eYDGGM+EJGscPA7HThXRO4M75MCjAz//FYHgRDgYhG5ASsODQGmAq3BcH6b7/e2OWaBMSYElIjIVmDyQed8DHgJKxheC/y9i997vw6DoYjcDNwCjBWRVW3eygQ+jvQC7TgPOCn88+PAe/SBYOgPBlm+s4arHl1MqIOgEzLw2MfbqG328+MvTyU3AQGxssHLT/69psM2tbIJ3HfZTN5et4/nl+064L2yei+/eW0DD76/hQU3Hs24wgwNiEr1b942PwexhiIDfDFpMqWT/UNtXoc4ME4cfCcygABfNcZsbPuGiMwBGttrnIiMAe4EjjTGVIeHQNu2yUTw8yGvjTE7RWSfiJyC1YO8or3rt6ezcbOngXOAl8PfW79mG2OujPD8BnhTRJaGPwEAFBlj9oR/3gsUtXegiNwgIsUiUlxeHlEvt0eqG/3c/NTSLoMOwPPLdrG5rCHubQKobvSxdnddl/udP2sYG/fWHxIIDzhXk59LHv6U8gZvh/sopfqtbcDs8M8XdvMclwCIyHFYQ5m1wBvAN0VEwu/NiuA8WViBslZEirBGGw+5Tvj7ojbbLxIRm4iMA8ZiPSM82CPAk8CzxphgZL9W58HQGGO2AbdiPfdr/SL83C8SxxljjsD6RW8NjzMfcAEOjfSt7z1sjPEYYzyFhYURXq77VuysOWRotDP3LdxEbTsTWWKp2Rfg4Q+3RrTvRbNH8PePt3W5X1Wjj3c2lGEiefiolOpP/gDcLCLLgYJunqMlfPyDwHXhbb/EmjizSkTWhl+3S0RWABhjVgLLgQ1YHa+DRxtzwyOStwPfabN9B/AZ1mSgm4wxLe1c5mWsSTsRD5ECSEc3RRF5xRhztoh8zhdd4VbGGDM2qguJ3A00ANcDJxlj9ojIEOA9Y8ykzo71eDymuLg4mstFpbbZz01PLmXRlsqojvvkB6cwNCc1Tq2C8novFz34Cdsqmzrdb3huKnd9aRLfemZFp/u1GleYwb9umEtBpjsWzVRKWfTZQwyEJ156DpprQngo9RVjzHNdHO8B7jXGHB/NdTvsGRpjzg5/H2OMGRv+3vrVZSAUkXQRyWz9Gesh6xqsqH11eLersR529ip/MERlN4YO/cFQHFrzBWMM9S1d91YLMtzsrm3vA1L79tY2E9KeoVKqnwlPynwe+GG0x0a0zlBEZmDNRtq/vzHmhS4OKwJeDA8jO4CnjTGvi8gSYIGIXAdsBy6OttGxZhMh1Rn9ejybxPmDoEBKBO1q9gVJc0Xe/jSXLi9VSvVNxpjRHWy/JoJj7wHu6c51I1ln+BjW1Nq1WDOLwBo27TQYGmO2Aoe3s70SmBd1S+MoK8XBvClFrGyT0aUrI/JSSXHGd92e225jRF4qu2qaO91va0UD04dlYxMimgD05emDyUpN7FpJpZTqyyK5m88NT2S52hjz9fDXtXFvWQI57DYunD2caFYb3HziOAoy4vvMLTvNxU0njutyP3/Q8GFJBfOmtDsx9xDXHDsmoh6nUkoNFJEEw0UiMjXuLellmSkOrjlmdET7Ds9N5bSpRUi8h0mBqUOzyE3ruhf3z0XbuPWkcQzP7XxCz8/PPYyCjMQlDFCJU9Pko7zeS2WDN+7Ps5XqbyJ5ePRPrIC4F2shpmDNJp3R+WHJJTPFyW2nTKC22d/pWr3huanMv34uhZkHr1mNj9xw7/A3r23odL+KBh93PreK+y+bxT8+2cZrq/fia3NDHFuQzo++PIUjR+eRkeB0cip+fIEgNU1+irdX849PtlHZ4CXFaefUKUVc7BlBZqqDLP377rM8Hs844A7gSqyEJvVYa+T+t7i4eEtvtm2g6XBpxf4dRDZj/WWt5otnhhgSn393AAAgAElEQVRjtse3aV+I99KKtqobfWwua+C+hZv4eEvl/hRoI/JSueWk8Zw6pYjCBC9JqGr08YPnV/Hmun1d7pud6uTSI0dwxrTBVDX6sIkwLDeVvHRX3Id1VWI1+QKs2FHDTU8upa6dWccicO2xY7j15PEJTR84gEU1VOTxeM4EnsNao9f2E4s//HVhcXHxa91qiMgZwH1YCbIfCU8safu+G6ujMxuoBC4JrysfsCIJhouMMUcnqD3tSmQwbFXb5KfJH8AfsKpWuJ02CjLcCRkabU9lg5dfvbqeF5d33Gtt6zunTuTqY0eRk6o3wf4oEAxRvL2ay//v0y4nTV125AjuOnMyuWn6byHOIr45hHuEq4C0TnZrAmZE20MMJ77ehJVXuhRYgpW0e12bfW7BSqB9k4hcClxgjLmk3RMOEJEMky4XkaexErfuX4wXwdKKpJad5iSbvjO8lJ/h5qdnT+WquaP48zslvL+p/JDE3U67cPb0odx00jiGZKfojNF+rKrJx21PL4to9vD8JTu5Yu4oDYZ9yx3Q5Q3GiZV95bYoz30UsDk8ox8ReQYrJ/S6NvucB9wd/vk54C8iImYAp6aKJBimYgXB09ts63JphYq93HQXueku7rt0Fg0t1nOi8novDrtQlJXCrBG5ZKY4SHfrOsL+btPeBioafBHv/5d3NvP7i2YkvPyY6tCVRBYMryL6YDgM2NnmdSkwp6N9jDGBcKWKfKCCAarLu6Yx5uuJaIiKXHaqk+xUJ8NyOxthUf1Vsy/AE59G98j+rfX7aPIFNRj2HZkR7pcR11ao/Tor4XQ/HSTRBjDGfCsuLVJKdcoXNNQ0Rd4rBAiGDMFIxlRVotRjVW7oSnfK4+wCRrR5PTy8rb19SkXEAWRjTaQZsDpbZ1iMVRQyBatacUn4ayagDx+U6iUOm3RrKFxrWPYpT2LNGO2MH3iiG+deAkwQkTEi4gIuxcoJ3VbbHNEXYlWqH9Cfljr8H2WMeRz2F/k9zhgTCL9+EPgwMc1TSh0s3e3gglnDeGdDWcTHzB2bR4ojvukDVVT+F7iGzp8b+jmwyntEws8Ab8OqM2gHHjPGrBWRXwDFxpiXgUeBJ8JL56qwAuaAFsnHy1ys7nxV+HVGeJtSqpfMGZNHustOoy+y2qW3z5tIjs4m7TOKi4u3eDyeC+l6nWG3Ft4bY/4L/PegbT9t83MLcFF3zt1fRfJR8R6s5RX/EJHHgWXAr+PbLKVUZ7LTnPzivGkR7Xvc+AImFuk8jL4mvKB+BvAUVtX3UPj7k1jrC7u14F51T5eL7gFEZDBfTM1dbIzZG9dWHaQ3Ft0r1dfVNvt5ZdVufvzvNYesOW11woQC/veSmZp9KDGiWXRvB74M3IWVBaZtQVU31nyN3wGvFhcXR9b9Vz0SaTDMBSZgTaYBwBjzQRzbdQANhkq1r8EboKbJxzNLdvKvz3ZS0eglxWHn+AkF3HbKeIbnppKXroEwQSIKhh6P50Ss4VE3nS+xqAdagIuKi4vf73nzVGciScf2DeB2rOm5K4C5wCJjzCnxb55Fg6FSnfMFQtZyCwEMpLrsuqYw8boMhuHnhP/ESmYSqWbgquLi4ue72zDVtUieGd4OHAlsN8acDMwCauLaKqVUVFwOG4OyUhiUmcKgrBQNhH1QuEf4ONEFQsL7PxE+XsVJJMGwJTzzCBFxG2M2AJPi2yyllOo/ws8In6PzxNydSQWeDZ9HxUEkwbBURHKAfwNvichLQMLKNymlVD/wZaxnhD2RApzV1U4ikiIin4nIShFZKyI/D28fIyKLRWSziPwrvCAfEXGHX28Ovz+6g/OeISIbw/v9oIe/S5/TZTA0xlxgjKkxxtwN/ARrseb58W6YUkr1I3cReT7SjmSGz9MVL3CKMeZwrIxhZ4jIXOC3wL3GmPFANXBdeP/rgOrw9nvD+x0gXBbqr8CZwFTgMhGZ2sPfp0+JKiWFMeZ9Y8zLxpjoEiMqpdQA5fF4hmEtn4jR6TxDO9vBWFpzmrYu6DfAKVhDtWA9u2zt1JwXfk34/XlyaOHW/WWhwvf/1rJQ/Ua38jOJyCuxbohSSvVTIzhwHWFPeDkwCXe7RMQuIiuAMuAtYAtQ05pWE6us07DwzweUcwJayzm11V5ZqGH0I91NVnh9TFuhlFL9V0+HR6M+nzEmaIyZibUk7ihgcozb0O90KxgaY/bEuiFKKdVP1ffW+YwxNcC7wNFATrhcExxY1ml/yadOyjlFUhYqqXVWz3A17dczFKxh6Rlxa5VSSvUfO+n5TNJWbg4crjyEiBQCfmNMjYikAqdhTYp5F6tc0zNY5ZteCh/SWs5pER2Xc9pfFgorCF4KXB6T36iP6KxqxdkJa4VSSvVTxcXFuzwez1Lg2Nicrnh3F/sMAR4PzwC1AQuMMa+IyDrgGRH5FbAca2UAdFDOSUSGAo8YY87qqCxUDH6fPqOzeoa6llAppWLjd8B8ur/oHqApfJ5OGWNWYWUKO3j7Vqznhwdvb7eckzFmN23WNbZXFqo/6fKZoYjMFZElItIgIj4RCYpIXaQXCM9qWt46A7WjhZ9KKdUfeTwewUpp2ZNASPj4o8LnUzEWyQSavwCXASVYKYG+gbX4MlK3A+vbvO5o4adSSvVHPwe+E6Nz3QH8LEbnUm1ENJvUGLMZsIen6/4dOCOS40RkOFYaokfCr4WOF34qpVS/4vF4zgS+C6TH6JRpwPc8Hk9E92AVuUiCYVN4KHOFiPxORL4T4XEAf8JKHxQKv86n44WfBxCRG0SkWESKy8vLI7ycUkr1DR6PJweran1Ph0cPlgY86fF4smN83gEtkqB2VXi/24BGrLUmX+nqIBE5GygzxiztTsOMMQ8bYzzGGE9hYWF3TqGUUr3pR8Q+ELZKB34cp3MPSJEEw/ONMS3GmDpjzM+NMXcQ2bKLY4FzRWQb1rqWU4D76Hjhp1KqnwiFDNVNPsrrWyivb6GiwUuzL9jbzUoYj8eTAtyIVWkiHlKAGz0eT6zWLw54na0zbHU1VhBr65p2th3AGPND4IcAInIScKcx5goReZb2F34qpZJYXbOfRm+ADXvr+XhzBct21FDV6CVoDCkOO2MK0jl2fAFHjs6lINNNfpoLu727GSH7vEQlsT4f+NfBG8Nl9x4BpmElT7kW2BjedzSwDbjYGFMdnstxH9YyiibgGmPMsnbOORv4B9ZEyv8Ct7ezOD9pdZaB5jKsDANjROTlNm9lYS3M7K7v0/7CT6VUEqpq9LGjspH739nMuxvLCHVweywpa+DNdfsAKMx0843jxnDezGHkpbtwOfpdUJxH7HOSHiwTOJl2giFWcHvdGHNheM5HGvA/wEJjzD3heoQ/wLofnwlMCH/NAf4W/n6wv2HlpV6MFQzPAF6L6W/Ui6SjwC4io4AxwG+w/tBa1QOr2kyCiTuPx2OKi4sTdTmlohYKGZr9QRw2we0cGMXIfYEg5fU+vvfcSj7ZcnAqy8ikOG388rxpzJtSRF567y45rmjw4vWHCBmDy2EjNy3qIL1//Z/H41mLVfcv3tYWFxdPO6ARItnACmBs256biGwETjLG7BGRIcB7xphJIvJQ+Of5B+/X5tghwLvGmMnh15eF97kx3r9gonSVgWY7cLSIFGEtGgVYn8hAmEyafAHqWwJ8sKmcdXvqaPQGKMhwc/LkQYzOTyMv3Y3dputl+5P6Fj+1zX6eLS5lVWkt6W47V8wZxYRBGRRk9t/HOXXNfj7aXMGdz66kqQfPAlv8Ib733CqOGp3LfZfNYkh2agxbGZmqRh8flZRz/zubKSmzygBmpTi4fM5Ivn7MGAZluTm0vF+XRsa8oe0b1c62MUA58HcRORxYirXeu6hNgNsLFIV/7qg8U9uCDMPC2w/ep9/o8pmhiFwE/AF4D+uTz/0i8j1jzHOdHjiAeANB9tS0cP+7m/nPit34gqED3n/gvS0MyU7hxhPGcvaMof36JjmQ1Db5eHLxDv7w5kbaDrC8smoPx08o4N6LZ4JAMGQIhscO7TbBJkJ2qjNphwbrmv08v6yUn/9nXczO+dm2ai56cBH/umEuw3LjNQHzUFWNPn704mpeW7P3gO11LQEefH8rzy/bxfM3HcPI/KjbFMl8jFho7zoO4Ajgm8aYxSJyHweO7mGMMSLSb573xUIkf2E/Bo40xpTB/ozob/PFwvkBrckXYMm2aq5/vPiQINjWntoW7v7POuZ/tpPHrz2KwdnxmmSmEiEUMry3qZzfv7Fx/za7TTh9ahGXHTWS2mY/97y+gVWlNWwpb9wfDEVgTH46M0bkcM6MIUwflk1O9MNxvabJF+D1tXtjGghblVY3c8nDn/L8zcdQlBX//x/BkOGVlbsPCYRtldd7uf6fxTz1jTnRfohN1OhZe9cpBUqNMYvDr5/DCob7RGRIm2HSsvD7kZRn2hXe3tk+SS2S/4G21kAYVhnhcf2ePxhkza5arv3Hkk4DYVsb99VzxSOLKatviXPrVDxVNvoOCITjB2Xw1DfmMLEok+8uWMk35y/nuaWlbNrXsD8QAhgDWysa+ffyXVz3eDGn3vs+979TQll9C8kwMW93TTM/fGF13M5fWt3Mj15cTXWTL27XaFXZ6OVv72/pcr+N++opq4+6UP2ObjUqeocUVDDG7AV2isik8KZ5wDq+KNUEh5Zw+ppY5gK1B9esDb+uC+eqFuBr9LOVAJEEtddE5A0RuUZErgFepR9nLo9GVaOfa/9RfMDNLhJbyhv409slA2rdVX/T6A1QWt0MwDXHjObucw7jB8+v4r6FJZQ3RH7jrGsOcP87mznn/o8o3l5NXbM/Xk3uscoGL7c8tSzqf+/Rent9GUu2VRGI8ANmdzV5g+ypjexD6Usrou4EfRx1g7rnow62fxN4SkRWATOBXwP3AKeJSAlwavg1WPfzrcBm4P+AW1pPIiIr2pzzFqzlGpuBLfSjmaQQWTA0wEPAjPDXw3FtURL57PMqGrzdGw15YVkpdS1998anOtc6EnDn6ZMYlZ/G1x5bzLbKpm6fb1+dl4seXMQTn26nJgG9omiFQoa31u1j076GhFzvzmdXUtUY3z+HaIJtffT/VxcC8f6L9GIV7D2EMWZFOIPXDGPM+caYamNMpTFmnjFmgjHmVGNMVXhfY4y51Rgzzhgz3RhT3OY8M9v8XGyMmRbe77b+tMYQIguGpxljXjDG3BH+ehFrXcqAVtHg5f53Srp9fIs/xPubypNiaEwdKtPt4MYTxpLqsvHz/6zrcG1dtH7/xkYe/2QbtX0sIFY2+vjLu5sTdr265gBLtlfH9RrpbgeRThKdOjQr2tO7wl/x5Aaccb7GgNFhMBSRm0VkNTBJRFa1+focWJW4JvZNzb5gjz8lL1iyk9o+PCymOpaZ4uC4CQX86tX1Xe8cpXvfLmFxAoYJo7Grumn/sHCi/OWdEiqjGHKOVprLwXHjC7rcz2ETTp1S1OV+rTwez2CiK3PXEw94PJ7IG6c61FnP8GngHKyHq+e0+ZptjLkyAW3r0/wxuFFVN/nwB7VnmIya/UG+/9wq4tWxv3PByqiePcZTIBhiQXFp1zvG2Po99TT7Y/9cvbLBS2l1E5vL6vnuaRNx2jvvHl5/wlgyUqJaKfF3rF5bIrixUqSpHuowGBpjao0x24wxlxljtrf56kkqtn4j+jW4h7KJxOQ8fUWjN0B5fQv76looq2+hsZvPU/u6Fn+QRz76nN0RTr7ojrqWAL9+dX2feK5c2+xn+c74Dll2JJa90coGL2+t28vVf/+M4377Ll99cBGPL9rO/ZcdQYb70GAnAtceO5objh9Lhjuy0UiPxzMDOJH4D5G2cgEneDye6Qm6Xr+VqIWh/Y7L0fOUW4WZbpxJnqi42RekrsXP+5vKeX3NXtaHM++kuRxMHpLJ6VOLOGVyEVkpDtLaueEko7pmP08sOmRGe8y9unoPPzxrClkpvftYKBAylCRo4szBPt1aydyx+T06hy9gzRq97enlrN5Ve8B7Ly7fRU2Tn8euOZKNe+v5oKQcfyDE+KIMzpkxlLx0F+nR/bv9Dol/jucGvg1cl+Dr9iv94+7UC9KcdmaPymVpDx7yf/3Y0WSnJu/z74oGL898tpP73ynBGzhw2LiuJcDeuhbe21iO076Wm04cx9XHjKYgI/mz73yypbJHKcgiFTLWrOObThyHoxc/NPkCIQJxXk7Rkc1lDfgCoW4nJfD6g6zZXcsVjyymxd/+o413N5bx7sYyjhiZyxGjcnDZbeysauLihxYB8Pi1RzFzRA4pXeSc9Xg8GcAlJP6+agcu83g8txcXF/fOp5Z+ILm7Jb0oN93Ft0+d0O3js1IdHD48J4YtSqy9tc1c9ehi/vDmxkMC4cH8QcP972zmkocWsbsmsZMwYq2u2c+C4p1d7xgjLy7fTU1T7w6VhnpxxrMvEOrRusbPKxq5/P86DoRtLdtRzSMffs4D723hP6v24A2E8AZCXPXo4kh7xkcBvfWX5eeL/NGIyO0iskZE1orIt8Pb8kTkLREpCX/PDW8XEfmziGwOT5I8or0LiMhsEVkd3u/P0o2ErX2ZBsMemDQ4k6Ks7vV0rj12DDlpvZulv7vK671c8chi1u+pj+q4LeWNXPzQIsrqkjf7Tku4p5Eon1c0xH2Re1dsvXjPc9ptdLdTXF7v5fonirv8sNYVf9BwwxPFkWSN8hC/Yr5dSQlfHxGZhlVq6SjgcOBsERmPlZJtoTFmAtY6yNZ8pW1LON2AVaqpPa0lnFr3PSMuv0kv0WDYA4UZbp68bg5pruieHx49Lp+vHT0qafJRtlXX7OePb25kS3ljt44vrW7mV6+u75MLyyPhC4aoa07cxKCQgXpv7/YMnQ4bvVVsZWReWreeq/sCIZ78dDs7q2IzErGntoVHPvicls5ntyZy4szBXOHrA0wBFhtjmsIVht4HvoJVcPjx8D6PYxUGJrz9n+HF958COeHcpfuFX2cZYz4NL7b/Z5vj+4Xkuxv3ISLC6II0Xrj5mIhrsc2bMoj7L5tFXnpyPjsrq2vhmSU9GyZ8eeVudiV4zVqs+HrYy+iO2gQG3/bYRRhTkNEr1z5mfH53yidR0+TjsY8+j2lbnvh0e1frgnu7pFFrIu01wPEiki8iaVgV7EcQfQmntrSEk+qc025nYlEmr91+PP9dvYeHP9jabr7Do8fmc/upE5hYlNGjQBgMGaoavTT5gmwua2BXTTMCjClIZ1R+OqlOe9xKRDV5A/ztg60xOddf3t3M7y6cQWYvz5RMBqFeHibNSXMyc0Q2W8oTPzdjdH56t45bt6eO+hgv7Wn2B1m2vZozpw/paJfefu7hAjDGrBeR3wJvAo1YhX4P6NJqCadDaTCMAZtNKMpK4WtHj+bsGUPYXdPC9spGmnxBctJcTB2aRVaKo0fPCI0xlNd7eWnFbh75aCv76tpfkD22IJ1bTx7PiZMKuz1zs7LBS9AYMJDitJMVnvHa4AvwZgflbpx24fSpgzlz2mCyUp3UtwR4Y+1eXl+zt92KHu9sKKPZF4woGNY2+2nyBfAHDU67kO5y7G9TovVGcebUKIfhY81pt/HV2cN5flliK/aMyk/rcgZne4KhEAvXl3W9Yzcs3FDGqVOLOhq67e2x//3XN8Y8CjwKICK/xurJaQmnTmgwjCG7TSjMTKEwM4XDR8Rupmh9i591e+q47anlXWYl2VrRyHefXcm4wnQeuGJ2VDeUigYvxduquP+dzWzYW4/DJhwzLp9vzZvAmIJ0fIFQu5+2jx2fz/dOn8Sb6/bx6/+up6LBR266k/NnDmPBjXP58zubeWfDgTcnbyDUZXaRumY/m8sa+ONbG/l4cyVgLYQ+cWIh3z1tImMK0slIcM8yxWnHYZOELjXI7QMTrcYPyqAww53QrDi3nDSegozof/fa5gArS2vi0CJYVVpDbbO/ow+au7AmrPSW/cOYIjLIGFMmIiOxnhfOBcZglW66h0NLON0mIs8Ac+ighJOI1IVLPC3GKuF0f7x/oUTSZ4Z9XH2Ln/+u3sslD30a1Y1oS3kj59z/Ect2VHf10B+wZt5d9/gSbnpyGWt31xEMGbyBEO9uLOeCBz7hD29sbDd4HTMun5tPHM9Vj37GA+9tYXdtC75giH11Xh76YCuXP7KYq+aO4pTJgw45trMMNa1LGL7yt0/2B0Kw6gG+t7Gcc/7yMS+v3E1DgjO0OGzChKLEPT/LS3d1mS4sIe1Ic3Hd8WMSdr0Up42TJhV263lhMBSKWz3EmiZ/Z8PW79N7vUNf+PqtnheRdcB/gFuNMTVoCadOaTDsw4LBEMu2V/P957uXF90XDHHNY0vYXtn5zM/aZj8/enE1K3d2vGTgycU7Dpk8YhP43pcmcctTSzt8PtPkC3Lr08v41injD7mp220d//PbWt7QZRLs/3lxDTsTPBEnO9XJCRMKE3a9uWPzo82AEhcOu42vHjGcYTmpCbnez845jJy07vf67XFaDmITgY5PXQz01rqhlvD1ATDGHG+MmWqMOdwYszC8TUs4dUKDYR9W3uDjW8+s6HrHTviCIW5+chkVnfQq65r9vLV+X5fn8h60cPnkyYP4sKSCupbOJyo0+YK8uW4fZxw2+IDt7eWDBKht9vG/b2/qsj0Af15YktDeocNu4/I5IxN2vZtPHNsngiFY6QMfuKLd9dgxdcSoHL502GDc3Ux5mOK0Myw3PkF7eG4qro6XenxG75VUcgJLeuna/YIGwz6qxR/kwfe3xKTE09aKRt7ZUNZh7cQ31u6NqPrCjqomRuR9cZM5bWoRL6/cHVEbXlqxm9PbBMPCDDdOR/sfsZt9QT4sqYjovG+u20djAlKjtZWd6uSoMXlxv87o/DSGZCemJxapsYXpPcq81JW8dBd/vnRWxEuV2pOZ4mTOmJ7lM+2IZ3RuhykUw6nQ/gUkei1MAJivqdh6RoNhHxXrtF9/fXczlQ3tP84o62Bm6sFeW7OHK+aM2v860+2MePF8VaPvgJ7gFXNHdjgxxB80EZdGCoZMwvNm5qS5+H/nT4v7zNLfX3Q4+d2YQBJPmSlOrjlmNF8/dnTMz52d6mTBjXMZGoMPAKdNjU+JvzOnDenqOea9JP6+agP+lOBr9jsaDPuoHVVNMU0Gvb2y4/MNj3BI6e31ZZw9Y8j+oFbb7It4+UZhpnt/LzfFaeMSz4gOM4s47ZFnPHHZbTh6YbnDsJxUvj0vfj2kS44cwaSizG5NIIm3nDQXt8+bwPfPmByzzDRjCtL5963HMqYgA1sMTlqY4WZcYWwnOo3MS2NoTpfZ1txAop+lGRJXP7Hf0mDYRy3+PPZlI7dWtD+KcurUoohuakOyU0hz2rnnq1bptP+u3ssFsyJLQvGVI4bx39XWbO1fnDeN/E6CaJrbzrwIK4t/uU1wTqQ0t4PL54xkztjYD5dOGJTBnadP6rW1lJHISXNx9dGjeO32ExhX2L2F8WBNwvrmKeN59sajGVOQHrPedkGmm5+dMzUm52r14y9PieTD3x0xvWjkeuu6/YYGwz4oEAyxcW90SbAj0VE+0Qy3g/MiCGp3nj6JnDQXx40v4OvHjubjLRV4Rud1uRYsO9XJSRMLWbihjEs8wzltSlGneVmzUpx857SJXRY+ttuEW08e32sTTPIz3PzlsiPwjMqN2TnHD8rgievmUBinLEKxlOZ2MGlwJv+68Wj+cvksJhVlRnys22Hj0iNHsPC7J3HTiePikjVpxvBsvtxxtpionDJ5EEeOzuu0p+7xeHKw8nUmOkuCHbggfH3VTX1jmpo6QLzGWDpaH5WV6uTHZ02hot7b4cSV2+dN4ORJhdhssn+YbEh2Cr96dR1/veIIbn1qGRXtPJPMSXPywBVHcM9rG/jWKeP52tGjyY1gcsTI3DT+eNHhfPfZle0+P7TbhL9cPiuSYau4Ksx089BVs7lvYQn/7GHB3y9PH8Ld506lMLN3f6doFWS4OXvGUI4em09lo4+PSir4ZEsFa3bVUd3kIxAypDhsjClMZ9aIXI6fUMDhw3PISnWQ6orfLSgnzcXd505lc1kDG/d1/8Pl2IJ07vnK9Ej+3R4NeOmdyhVerIX1rwOIyGPA2UCZMWZaeNvvgXOw1iRuAb4eXn+IiPwQqzhwEPiWMeaN8PYzgPuwAu4jxph7OIiIuLESd88GKoFLjDHb4vabxonEa6mIiKQAH2CNZTuA54wxPxORMcAzQD6wFLjKGNPpLAyPx2OKi4s726Xf+d3rG3jgvS0xPecDlx/BWTM6/qRc2eBla3kj979bwrrdddhtwgkTCrnhhLEMyko5ZBZdsy/I7ppm5n+2nTOnD2X5jmpeXL6LigYveekuzjt8GEeNyeO/q3dzkWckw3NTo+rFNXoD7K5p5q/vbubV1XvwBw0uu41zZg7llpPGWcO2cbyZRqO+xcqW8835yymNcu1jQYaLP158OIcPz0nasl4Ha/IFaPAGrE92AhjrA0xOmrPT9aXxUFbXwk1PLWXZ9uiz0kwfls3/Xe1hcFbX8c3j8fwU+Cm908nwA78oLi7+FYCInAA0YFWjaA2GpwPvGGMC4dylGGO+LyJTgflYJZ+GAm8DE8Pn3QSchpXdZglwmTFmXdsLi8gtwAxjzE0icilwgTHmkvj+urEXz2AoQLoxpkFEnMBHwO1YY9svGGOeEZEHgZXGmI7qZwEDMxh+WFLOVY9+FtNzvvPdExkbwaSC2mY/Xn8QBDJcDtK6CGB1zX4afQEq6r00+gKICMYY0l0O8jPcpLvtZKd2/ybf4A3Q5A0QCBkcNiHD3XWbekt5vZdtFdYHik82V3Y409UmcMTIXG6fN4FJQzLJT3f3St7TgaKq0ctrq/fyi1fWRVTf0O2w8YMzJ3Pu4UM7fb7dlsfjWQjM62FTe2JhcXHxqa0vRGQ08EprMGxLRC4ALjTGXBHuFWKM+U34vTeAu8O73m2M+Rr9LY4AACAASURBVFJ4+wH7tTnXG+H9FomIA6siRmGyLcqP2x0l/AfROmPDGf4ywCnA5eHtj2P9oXcaDAei8YMysNskZoVd89JdEU80yU51QhSTN7JSnWSlOuO2Ji7D7eiVSTLdUZjppjDTzf2DZ9HiD9HkC1Cyr4HKRh8YyE5zMrEog8wUJy6HrU/kHR0I8tLdXDh7OKdOLeLdDWU8uXg7G/fW4w9+8f/LYRMmFmVy+ZwRnDZ1MNmpDlKcUf276+2SRkOj2PdarDWRYLX70zbvtS3PdHBppzntnGt/Cahwr7MWa+QvssXCfURc7zAiYscaCh0P/BVrnLomXHASOqmJJSI3YFVdZuTIxGX86Csy3Q6+dFgR/13dfpWIaH392NE9WsisopOd6qL1s0Fv1QJUB3I77RQ57Vxy5AhOm1pEIGRo8QcJhgx2m+xPwp6b5uru8o7e/g8WURdWRH6EtVD/qfg2J7nENRgaY4LATBHJAV4EJkdx7MPAw2ANk8anhX1XRoqTu740mbfW7Tvg02t35KY5udgzAkc3KoYr1d+ISMRDn1Hq7RJOXWbPEJFrsCbWzGszjNlZCaeuSju1Pb40PEyajTWRJqkk5O4YnrH0LtZsq5zwHxj0w5pYsVSUncIPz5zS4/P86dJZ3a5tqJSKWG/fyzrNjRieGXoXcK4xpqnNWy8Dl4qIOzzBcQJWntUlwAQRGSMiLuDS8L4HexmrJBTAhViTdJKuAxO3YCgiheEeISKSijUjaT1WULwwvFvbmlrqIKlOOxccMYzLjxrR9c4d+MnZU5k1IkcnZygVf+9jzersDX7gvdYXIjIfWARMEpFSEbkO+AuQCbwlIivCExgxxqwFFgDrsJZm3GqMCYYfZ90GvIF1714Q3hcR+YWInBu+3KNAvohsxpog+YO4/7ZxEM/ZpDOwJsjYsYLuAmPML0RkLNbSijxgOXClMabT7v1AnE3aVnWjj+eWlfLb1zZEnIcz3WXnDxcdzrHjC/p0JhOl+guPx3MW1hKF7F64fC1waXFx8eu9cO1+IZ6zSVcBs9rZvhVrPYuKUG66iyvnjGLe5EHc89oG3l6/j45iotth47yZQ/n2qRMpzHThtCc6GYZSkWnxB6lt9lPX7GfN7lo2lzUQCBpy013MGpnD8JxUUpz2eD3fi4dF9F6OUDcHzghVUYpbzzCWBnrPsK36Fj8N3gDr99Tzyf9v787j46rLxY9/nlkyk0ySZm2a7gst3ShtCaVAKWUTlLJ4QWSTRQREQdArXLhXf/Lzyu+nl6ugKAgiiCwVFEXRK1spSwu0TRdoaelemi5JszV7Zsv3/nFOadpmmWT25Hm/Xnl15syZc55mJvPMdznfZ2sNFfVtOAQmFPs4ZUIRE4Zmk+uN78oeSkWjxR+iutnPE0t38KdVu3tckH76iFy+deZEZo/Oj8uSbTEmZWVli4Avkdgl2cLAC+Xl5Vf2uqfqln5ippkcr5scr3VN34JJxQTCHQjWtHGlUl1Ns58nl+3gkbe2ddu70dn6PY3c9PQqxhf5ePjq2YzKz0qZYsfd+ClwIZCVwHP67fOqKGjLMM21BcI0tgXZ19jGml0HqGn2k+N1M3NUHmMKs8j2uMjx6pihSr69B9r4ym9XsK26fzVoHQL3ffE4zj+uNFXHwQWgrKxsCXAqial6HwSWlZeXn5GAcw1oKf0VS3UvFO6gusnPQ29u4cXVe7pcYsohcPaUEu75whSGDfGSqa1HlSSVDW1c9uj7fV63tbMOA/f8eR3hDsPFM4eTnbpf8r4CfEJikqEfuDoB5xnw9CrsNNTcHmLZtlrOffAdnltR0e1aix0GXttQxeceeJsXVlZEXJVeqViqawlwaz8WMO/O915az/aarsuRpYLy8vLdwG1AvINsAb5VXl6e7OsbBwRNhmmmLRBi6dYarntyBY3tod6fAATDhh/87WOeWLqDhjZNiCpxAqEO3thYRfnO+pge9xvPrqa6qdcFV5Lpd8BviF9CbAF+U15e/mScjj/oaDJMM5WNfm5btLrLGn+9+cWbW/l4T2O3dQ2VirUDrQH+8+UNve/YR7vr23hl/b6UfS+Xl5cbrAvQ45EQW+zjanX7GNJkmEZqW/x8+/m1Ua1V+q0/rKG2JaW/UasBZOXOepr8kfVg9NWv395uVQNJUZ0S4m1YFXyiXZ0maB/ntvLy8m/bx1cxoskwjdQ2B1hb0fcCpZ3VNAdYvSu2XVbJYoyhpsnP9upmVn9az7rdB9jX0EZTe7JWxFKdHWgN8OSyHXE7/p4DbdSn+Dh4eXm5sbsypwDLgFas6wL7Imw/bxkwWbtG40Nnk/ZDuMNQ1+InYE9cyXA5KPR5+lv2JSLGGF5aE5tx8qff38UpE4rS+pKLmmY/Sz7Zz6+WbGVn7aE1h50O4ewpQ7n97EmMLsgk25O+/8d0Fwh1sH5vQ1zP8WHFASaV5MT1HLFgT6o5o6ysrAyrtfhFrJmgWXQ96zSIlQA9WBV/flau15fFlSbDPqpu8vPyh3t57J3tVDa2A1A6xMvN88dz/oxSinO8cTlvQ1uQlTvrYnKsdXsaaAuG0zYZVje1c9PTq1iz6+hWcrjD8OrHVby2oYr7Lp7OBccP7/b/2R60rtEMG6ueXX5WBm4tcxUz/lAH7cHeq8pH471tNVw8a0TavG52QruyrKwsD5gLlAELsArzerAS5F6sRbfLgQ/Ky8uj6w5SEdFk2Af7Gtq44rEPDmuJWNvbufflDTz9wac8+7W5DBsS+4QYDHewr6E9JsdqaAvSEd/PqLipbfZz66I1XSbCzoyBf//LeobmeDlz8tCjWu2VDe08+s42nl9ZQWsgTH6Wm2tPGctVJ42hOPWX/UoLwXD832T1LUH8oY60SYYH2QnuFfvnR0kOR6FjhhGrbfbzjWdWH5UIO9tW3cKtz62mLk6D+g6JXTdsDA+VMKGODnbVtbJ8e+Qt5Htf/vioCUN7D7Rx0a+W8uSynZ+ti1nfGuTBN7Zw1eMfpPqUfdWZ2Mu+KBUlTYYRqm0JsCaCySvln9ZTH4dk6HU7GVsYm+UOS3I9aZMMa5v9bKps4kd/38Br66v4xeItfXr+7vo2KuoOXezd4g9x/yubqGrsOuFtrmrmhfIKQglo1Qx0Ga74f7yUDsnE49aPMRU9fRdFwBjDouW7It7/+fKKmMeQ43VzyjFFMTnWzFF5ZKf2YscYY9hV18o1T6zg3Aff4fGlOyjMzuDtzdV9PtaLq3Z/drupPcTf1/VYEJwnlu6gLsVnKaaDDKeD3Mz4vs/mHVOIy6EfYyp6+i6KQDDcQVVT5ON1lQ3tcRkvWTijNCYtuq+dNp6sFC/xtK+hnUsefo+P9zYetr0/11gfaAt+NvM3EAr3ep1mbUugX+cJhDrY39jOngNtbNjbwNqKejZXNlHZ0EZ1kz9lLxCPl8wMJ7NH58f1HNNHJKOOrhqIUvsTMUW4HA4KfBkR71/gi8+sxByviy9ML+Uf6/b1+xgTirMZV+SLYVSx1+IP8dPXN1PdHJuxu2yPE7fT+haR4XLgdAjhHhJTbqaLvlwl0xoIUd8SYNGKXfxhZQU1zUe3KscV+bhh3ljOmTqMQl8GrjSb8NEfOV43Xz99Am9t6ntrPhLThuemfA+HSh8D/y8yBhwO4Yo5oyPe//I5o+ISx5DMDL6/cApD+lm+xukQHrl6NkUpXjm8qT3Ey2uP7spcv7eBk8cX9vl4F80agdhNap/HxdlThva4/1fmjiE/M7IvP7XNfp5ctpP597/FL5ds6zIRAuyoaeF7L33MOT97mzc/2T9oFk2fUOxjWG58Ljf61lkTKUzx97JKH5oMI1SS62VSSXav+00tzaU4jn+gRdkenv3aSWRl9K0ck0Pg4StnMzIvM06Rxc6+hjYCXXQzP/vBLq6aG/mXEoDiHA/HDD30uuV43fzH+VO7/UIxMj+Ta04eizuCyR/7m9q54amV3P/qph5bmp01toe46elV/Ncrm+Iy0SrVFPo8PPDl42N+3JmjhnDi2Ph2warBRZNhhIqyPfz22hMZ2sM1aMNyvTx2zQlx/bbqcjo4dlgOL982j4lDe0/OYCWE528+mXkTi8hKg26lUDdjettrWsh0O5lamhvxse75/GQKsw5v5Y3Iy+Tvt83j/ONKcdn9oV63gyvnjOJPXz+FkghaMrXNfm78fTlrK/q3wspzK3bxwBubaWgb2EvHORzC1OFDuGT2iJgd0+Ny8PPLZ1Hg01ahih2tdN9HVY3tPLF0B88t3/XZAsS5XhdXnTSa604dF9EHaazUNPt5d0s1v3xzK9uqj14Yf2iOhxtPG89Fs4ZTnO35rKswEcIdhvrWwKFJIwLZGa6IkvGntS2cfv9bXT5Wkuvh11efwHf/+FGvFdO/ffZErj1lLHlZXXd5NrcHafaHCYY7yHA5yPW6yezU4m5uD9Ie7EAE8jLdOO1xvlZ/iN+8u50H3ujbZR5dee5rJzF3fGFcl/JLBbXNfm55ZhUroizl5HYKz9xwEjNH5eFJvWLVA/tFHOA0GfZDINRBfWuAtmAYwboGMD8rIyHXVXWlpsmPP9RBdXM7DW1BstwuSvO8eFxOCrIOfYgnQos/RFN7kH+s28dzy3exv9FPsKODXK+b+ZOKuWn+eAp9GT22nutaAlz/5Ao+3N11q2tkfiYPfnkmb36yn+dXVhxVuWD26HzuPPdYpg3PJbcf46t1LQF21bbw8Fvb2Fbdgtft4OKZI7jg+OEU+NxUNvpZcP+Sfs04PVJxtof/uX1e3JbxSyU1zX6+/9J6/rm+sl/PH5Lp5rfXlTGtNJfM1JwNrckwjWkyVDFT1djO//ufjfzjo32EesgUU0tzue+L05k8LPewllhnK3fWcdmj73dbt9HtFM6dNoy7zptMXUuAqsZ2PC4H44p85Hrd5Pdh9m9nNc1+7vzjhyzpYgakx+Xg5VtPZdHKCp5ctrNfx+/KE9eVcebkkpgdL5UdaA1Q/mk933lhLY1tkZd2+vz0YfzggmkUZ2ck9MtdH2kyTGMp+fVKpZ/d9a1c/tgH7K5v63XfDfsaufTX7/PQFbM4fVIxvi66TqeU5vDAZTP5zgtru2yBGWNdd1noy2B0QWxW5mlqD/L///lJl4kQrIWnKxv9/Kl8d5eP99fj7+7ghDEF/Z4lnE7ysjI4Y1Ixb3zndJZvr+ORt7axsbKxyy89OR4XF80czlfnjaPAl9Ftd7dSsaDJUEWtqrGdK3+zPKJEeFC4w3Drc6t55oaTOHFcPm7n4S3EbI+bc6aW8NadZ/DUezt5cfVuGtqC5GdlcFnZSK4+aQyFORlkumP3Fm5qD/GX1T0nOq/bGfNitet2N+APhem6kk9i1bUEONAaYHOVNR47pTQnqpZ2V5xOB0NzvFxw/HBOmVBIMGxobAtS2dhOuMOQ7XUxLNeLx+UgN9ONN/XGBtUApMlQRaU9GObhJVvZVdf9Aubd6TBwy7Oref3b8xmae/QHns/jwudxcde5x3Lz/PFWJ5SxWhfxGJ99e3N1j+OARdkZcVnEu8kfIhhK7nCFPxhmR00Ld734ER8dMVZbNiafn1w6g9EFWTFfTOLg2PGwIV4mDUv9uoRq4ErZzneVHhragrwQRbdhQ1uQdXt6vjzB43YyNNfL0BwvQ3O9cZuoVNNLovO4nLSH+lqkPDIdSRy7N8awZX8zF/1q2VGJEKzF5y98aCmf1h49Y1mpgSJuyVBERonIEhHZICIfi8jt9vYCEXldRLbY/+qVs2ls2dYa2oLRJYifL95CbYyWXovGiPyeFyRoC4bxxWkWozOJl1bUNge49bnV+EPdr6fbEgjznRc+jFt5MqWSLZ4twxDwr8aYqVgVnb8pIlOBu4HFxpiJwGL7vkpDB1qt9Tij9dHuhqgTaiycMqGIjB66AetaAhRmx34SR36WG5czeclwX0Nbj3U6D/podwP1g2QZOTX4xC0ZGmP2GWNW27ebgI3ACOAi4Cl7t6eAi+MVg4qvYNjErKXQWyWJRMjxurh+3tge92luD/Vp0fZIzBqdT1YSJ4ks21Yb8b4fRlDTU6l0lJAxQxEZC8wClgMlxpiDZRcqgS4vsBKRm0SkXETKq6vjs+q9ipaJeE3OSI6VbD6Pi5vnT+DKbhZlL/BlML7Yx9V9XB+1NzfPH0+2N3kzSfvyGvZ0/ahS6Szus0lFJBt4EbjDGNPYeUkwY4wRkS7/uowxjwGPgXXRfbzjVH3ncjisa78i6GKL5FipoMCXwb+dN5mbTx/PM+9/yobKRjLdTi4rG8XM0XkU+jxcfuJoHn17e49jbJEaVZDJ+OLI1piNl77UHJw2PPJ1YZVKJ3FNhiLixkqEzxpj/mxvrhKRUmPMPhEpBfbHMwYVP7mZbhbOKGVtlF1n44p8KXUt2ZAsN0Oy3Nz1+cm0B8M4HXJYMeRCXwbfXziV7720PqrziMAjV51AcQ+LvyfChKE+irM9vdaPHFfko2QQLBunBqd4ziYV4LfARmPMzzo99DfgWvv2tcBf4xWDii+nQ1g4Y/hnlR/66xsLJlAUh4kp0XI7HeR43YclQrAu9fjCcaWcObk4quPfee6xjC2Mzeo50SjIyuAnl87ocR+nQ/ivS2fEZQKRUqkgnn1TpwJfAc4UkbX2zxeAHwPniMgW4Gz7vkpT2R4n500f1u/ne90O5k8qTmhFjVgo8GVw/6XHc+60/v3f7zz3WK6cMzqpY4UHuZwO5ozN5/FrysjPOjqe4mwPz9xwEtNKc9PudVIqUrpQt4paRV0r5z/0bp8WXj7ox5ccx0Uzh8d0WbVEqm8J8O6Wau758zpaAr1fHjJ8iJdHrj6B8cU+clIgEXYWDIepbw2ycW8j72+vRRBOm1TEMcXZFPpSeoHsVKHfFNKYJkMVtXCHYVNlI19+9IM+rdt5+1kT+eq8sQzJTO+ut2A4TH1LkLc2V/PUezvZVNl02KzLTLeTmaPyuGXBBKaU5iZ9jFDFjSbDNKbJUMVEKNxBRX0btzyzik8qm3rcNzfTxb0XTOOsKSUDqlKDMdZ1l6EOgz8YpsNYY20elwOP25H2SV/1SpNhGtNkqGKqttlPTXOAX7+9jX98tI9A+NDlB9NH5HL7WRM5fmRe3BbbViqJNBmmMU2GKi5a/CFa/CGCHQZjDC6HkOF0UNBDhXul0pwmwzSWnrMWVMo7WH5JKaXSgfZTKaWUGvT0q7tSXQiFOz6bDOMQq8J9XpZOgFFqoNJkqNQR9je28/zKCp58b+dnVTlOHl/Idz43iUklOQNqBqxSyqITaJTqZHd9K5c+8j6Vje1dPn7H2RO5/tRxmhBVV3QCTRrTMUOlbLXNfm5+elW3iRDgwTe2sKmyMYFRKaUSQZOhUrbqZj8f7+090f3s9c0c0IrvSg0oOmaolO2f6yoj2u+D7XW0B6OvZRiNA60BWgNh2oPWeqhet5OsDJ3ko1R/aTJUytbcHvm6qh1JGmuvbfazpaqZBxdv5oPtdYc9NmdcAXecNZFJw3Io0sUNlOoTTYZK2aaPjKyKe36WO+oajv1R2dDOjb9fybo9XXflrthRx5WPL2dqaS5PXFfGsCGZCY5QqfSlY4ZK2U6ZUIQngvVSrzl5LPm+xHZH7m9q56rHP+g2EXa2YV8jV/xmOft7mAiklDqcJkOlbLleN3edd2yP+wzN8XDlSaNxJ7C2X6s/xE9f28y26paIn7OjpoWfvLKJlj6U1FJqMNNkqJQtM8PJJbNH8r3zp3TZQpw8LIcXbzmFklxvQuNq8od4ac2ePj/v5Q/30tSHcVClBjMdM1Sqk7ysDK46aQwXHD+cpVtqWLenAZ/HycIZwynO8SR8YooxhsUbq/CH+j57NRDu4NWPK/nK3DE4kjDGqVQ60WSo1BEyM5xWK/GEkVxywsikxtIaCPP6hv39fv7ijVVcMnsE2V5dMUepnmg3qVIpLBQ2tAT639XZ7A8T7Ej9JReVSjZNhkqlMJdT8GX0vwMn2+PErV2kSvVKk6FSKSzT7eT0ScX9fv5pE4vJjCKZKjVYaDJUKoU5HMJ500twO/veunM5hIXHl+LUlqFSvdJkqFSKy/a6WThjeJ+fd970YeR4tFWoVCQ0GSqV4rI9Lu7+/GRGFUS+vNrI/Ey+t3CqziJVKkKaDJVKAyW5Xp6/6WQmDs3udd8Jxdm8cPPJDEvw4gBKpTPtQ1EqTQzPy2TRjXNZW3GAny/ewro9DYc9Pm14LrefNZFZo/MpztGqFUr1RdySoYg8ASwE9htjptvbCoDngbHATuAyY0x9vGJQaqApyvFw9tQSZo/JozUQptkfAgM+j4usDCeFWrpJqX6JZzfp74Dzjth2N7DYGDMRWGzfV0r1UYHPw8j8LCYPy2VyaS6jCrI0ESoVhbglQ2PMO0DdEZsvAp6ybz8FXByv8yullFKRSvSYYYkxZp99uxIo6W5HEbkJuAlg9OjRCQhNqfgLhMI0tYfoMAYQPG4HuTrjU6mkS9oEGmOMEZFuF000xjwGPAZQVlamiyuqtFXX4qehLcTKnXUs21rD5qom2gJhnA4Hw4Z4mDOukFPGFzKmKIv8zAzcERQYVkrFVqKTYZWIlBpj9olIKdD/5fiVSnF1zQHWVNTz4BtHz/w8aFt1M8u21vIAkJvp4pq5Y7hq7hiKsjNwO52JDVipQSzRX0H/Blxr374W+GuCz69U3PlDYXbVtnLD71dyw1Pl3SbCIzW2hfjlkm2c/dO3eWV9FfUtgThHqpQ6KG7JUEQWAe8Dx4rIbhG5AfgxcI6IbAHOtu8rNWC0+EO8v62Wcx54mzW7DvTvGIEwty1aw7+9+BG1zf4YR6iU6krcukmNMVd089BZ8TqnUsnUGrAS4Y1Pl2NiMMr92oYqmv0hHrpill42oVSc6Ui9UjFSUdfGzc+sikkiPOi9bbX892ubaGwLxu6gSqmjaDJUKgZqmv3c8swqwnGoKr9oRQWfVDZhYplllVKH0WSoVJRC4Q5eXLWb7TUtcTvHrc+tpkbHD5WKG02GSkWprjXAo+9sj+s59jf5+XhvY1zPodRgpslQqSht3NtIXQIug3jwjS3UtWjrUKl40GSoVBT8wTAvrt6dkHOtrTiAP9SRkHMpNdhoMlQqCo3tIT7cHdlF9bGQiBaoUoORJkOlohDuMHxa25qw832UwMSr1GCiyVCpKIQ7EtttqS1DpeJDk6FSSqlBT5OhUlFwOhyIJO58xTm6LJtS8aDJUKkouBwwvsiXsPPNGDEkYedSajDRZKhUFHIzM5gxMi9h58vLcifsXEoNJpoMlYpChsvBZWWjEnKuOeMK8Li14K9S8aDJUKkoTSj2JWQs746zJpKflRH38yg1GGkyVCpKhb4Mbj3jmLieY2R+JhNLsuN6DqUGM02GSkXJ6XSwcEYpU0pz4nJ8EXj4qtkU53jjcnyllCZDpWKiMNvDr66cjccV+z+pr80bx7gEzlhVajDSZKhUjIzIz+R315+IyxG7Cw8/N62EWxYcQ45XZ5EqFU+aDJWKEY/LyczReSy6cS7ZHlfUx/vSCSP58b8cR4FPJ80oFW+aDJWKoUy3i5mj83jt2/NZMKm4X8fIz3Lz+6/O4T/On0KBT1ecUSoRov/6qpQ6jNvpYHheJg98eSZbq5v5xeItLN1agzE9P29Yrpeb5o/n/BmlFGV7cMawu1Up1TMxvf2FpoCysjJTXl6e7DCU6peG1iDNgRCb9jWydGsNG/Y10hYI43Q4KB3i4dRjipk1Oo/C7AwKfZoE05i+cGksLZKhiFQDn3bxUBFQk+BwepNqMaVaPJB6MaVaPKAxRSLV4qkxxpyX7CBU/6RFMuyOiJQbY8qSHUdnqRZTqsUDqRdTqsUDGlMkUi0eld50Ao1SSqlBT5OhUkqpQS/dk+FjyQ6gC6kWU6rFA6kXU6rFAxpTJFItHpXG0nrMUCmllIqFdG8ZKqWUUlHTZKiUUmrQS9tkKCLnicgmEdkqIncn4fxPiMh+EVnfaVuBiLwuIlvsf/MTHNMoEVkiIhtE5GMRuT2ZcYmIV0RWiMiHdjz/194+TkSW26/d8yKS8MU3RcQpImtE5O+pEJOI7BSRdSKyVkTK7W1Jez+JSJ6I/ElEPhGRjSJycpLjOdb+3Rz8aRSRO5L9N6cGjrRMhiLiBH4FfB6YClwhIlMTHMbvgCMvsL0bWGyMmQgstu8nUgj4V2PMVGAu8E3795KsuPzAmcaY44GZwHkiMhf4CfCAMeYYoB64IUHxdHY7sLHT/VSI6QxjzMxO184l8/30c+AVY8xk4His31XS4jHGbLJ/NzOBE4BW4C/JjEkNMMaYtPsBTgZe7XT/HuCeJMQxFljf6f4moNS+XQpsSvLv6a/AOakQF5AFrAZOwlo1xNXVa5mgWEZifXCeCfwdaxmtZMe0Eyg6YltSXjdgCLADe4JdsuPpIr7PActSKSb9Sf+ftGwZAiOAik73d9vbkq3EGLPPvl0JlCQrEBEZC8wClpPEuOzuyLXAfuB1YBtwwBgTsndJxmv3IHAX0GHfL0yBmAzwmoisEpGb7G3Jet3GAdXAk3ZX8uMi4ktiPEe6HFhk306VmFSaS9dkmPKMMQbrAy7hRCQbeBG4wxjTmMy4jDFhY3VtjQTmAJMTde6uiMhCYL8xZlUy4+jCPGPMbKyu/2+KyPzODyb4dXMBs4FHjDGzgBaO6H5M1vvbHsu9EPjjkY8l829Opb90TYZ7gFGd7o+0tyVblYiUAtj/7k90ACLixkqEzxpj/pwqcRljDgBLsLog80TkYPmwRL92pwIXishO4A9YXaU/T3JMGGP22P/uxxoLm0PyXrfdwG5jzHL7/p+wkmPS30dYXxZWG2Oq7PupEJMaANI1Ga4EJtozADOwuk3+luSYwIrhWvv2tVhjdgkjIgL8FthojPlZsuMSkWIRybNvZ2KNX27ESoqXJjoeAGPMPcaYkcaYsVjvmzeNMVclMyYR8YlI9OUHFwAABAlJREFUzsHbWGNi60nS62aMqQQqRORYe9NZwIZkxXOEKzjURQqpEZMaANJ2BRoR+QLW2I8TeMIYc1+Cz78IWIBVRqYK+AHwEvACMBqr5NRlxpi6BMY0D3gXWMeh8bB/xxo3THhcIjIDeArrNXIALxhjfigi47FaZQXAGuBqY4w/3vF0Ed8C4LvGmIXJjMk+91/suy7gOWPMfSJSSJLeTyIyE3gcyAC2A9djv4bJiMeOyQfsAsYbYxrsbUn7HamBJW2ToVJKKRUr6dpNqpRSSsWMJkOllFKDniZDpZRSg54mQ6WUUoOeJkOllFKDniZDpeLErkRRlOw4lFK902SolFJq0NNkqBJORO4UkW/Ztx8QkTft22eKyLMicoVd22+9iPyk0/OaReR+uzbiGyIyR0TeEpHtInKhvY/T3meliHwkIjfb2xfY+x6s0fesvWLPkbFdYNc1XGOfo8Tefq+IPC0i79u1827sdNx3ROQfYtXX/LWIOI445g9F5I5O9+8Tu9akUio1aDJUyfAucJp9uwzIttdUPQ3YjFVb8EysGognisjF9r4+rOXTpgFNwI+wlnj7IvBDe58bgAZjzInAicCNIjLOfmwWcAdWDczxWOuUHmkpMNdeoPoPWNUtDpphx3Uy8H9EZLi9fQ5wm33cCcC/HHHMJ4BrAOxEeTnwTM+/IqVUIrl630WpmFsFnCAiuVgFgFdjJcXTgJeBt4wx1QAi8iwwH2upuwDwin2MdYDfGBMUkXVYtSXBWtdzhogcXGd0CDDRfu4KY8xu+7hr7ecsPSK2kcDz9qLPGVh1/Q76qzGmDWgTkSVYSfCAfdzt9nEXAfOwFrcGwBizU0RqRWQWVomhNcaY2j7/1pRScaMtQ5VwxpggVpK5DngPq6V4BnAMVpHb7gTNofUDO7ASKcaYDg59sRPgNmNXRTfGjDPGvGY/1nmt0TBdfxl8CPilMeY44GbA2zn0I/8rvWzv7HGs/+/1WC1FpVQK0WSokuVd4LvAO/btr2MtkL0COF1EikTEiVWl4O0+HPdV4Ba72xURmWQv8BypIRwq33TtEY9dJCJee3HoBVjVUwDm2BVUHMCXObq1CdZC3Odhdd2+2od4lFIJoMlQJcu7QCnwvl2brh14165afjdWSaUPgVXGmL6U5Xkcq9zQahFZDzxKL8MB9gSXC+279wJ/FJFVQM0Ru35kx/UB8J/GmL329pXAL7HKU+3gUAWKzxhjAvZzXzDGhPvw/1FKJYBWrVAqAiJyL9BsjPnvI7YvwC4D1cvzHVhjo18yxmyJV5xKqf7RlqFScSYiU4GtwGJNhEqlJm0ZKqWUGvS0ZaiUUmrQ02SolFJq0NNkqJRSatDTZKiUUmrQ02SolFJq0Ptfen1Ffv4/5QcAAAAASUVORK5CYII=\n", "application/papermill.record/text/plain": "
" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "berkeley_fig" } }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcMAAAFgCAYAAADHDhKeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xd4XNW18OHfmqreZbn3jm1sPGDTi4EAoSb0EgiETkJCCEluGilfQtolhIQAF0gIxcS0QCBU08EYy71btnGRm3qXpu7vjzMysq0yI82MNNJ6n0ePNGdO2XI5a/Y+e68lxhiUUkqpgczW2w1QSimlepsGQ6WUUgOeBkOllFIDngZDpZRSA54GQ6WUUgOeBkOllFIDngZDpZRSA54GQ6WUUgOeBkOllFIDnqO3GxCJM844w7z++uu93QyllOqM9HYDVPclRc+woqKit5uglFKqH0uKYKiUUkrFkwZDpZRSA54GQ6WUUgOeBkOllFIDXlLMJo1EMBiissnHlrJGPiwpp8EboDDDzemHDaYgw0V+hru3m6iUUqqP6hfBsKrRx4vLS3nw/a2U13sPeO+Pb21iypBM/uesKcwckUNmirOXWqmUUqqvSvph0qpGL997biW/fGX9IYGw1fo99Vz16Gf8a8lO6pr9CW6hUkqpvi6pg2Fds58/vLGJhevLItr/V6+u59OtlRhj4twypZRSySSpg2G9N8D8JTuiOuaXr66josEXpxYppZRKRkkbDP3BEM8W7yTaTt7OqmZ21zTHp1FKKaWSUtJOoKlp8vP8stJuHTv/sx1MH5aNzdZ/UglWNnip9waoqPeSk+okO81JQYYbkf7zOyqlVLwkbTA0xlDT2L3JMJWNPryBEKkue4xblXj+YIjPKxq5Y8EK1uyq2799ZF4av/nKdGaOyCHdnbR/zUoplRBJO0wK4HZ2r/luhw2HvX/0mHbXNHPBXz8+IBAC7Khq4spHF7NmV20vtUwppZJH0gbDFKedw4Zld+vYuWPzcdqT9lffr8Hr5w9vbqTRF2z3fWPgR/9eQ0VD+0tOlFJKWZI2ImSlOvnmyeOjPs5lt3Ha1KI4tCjxGr1BXlu9t9N9Npc10OQNJKhFSimVnJI2GAKMyk9neG5qVMecO3MomSn94xlaIGQIhLqeTtvkb7/nqJRSyhLXYCgi3xGRtSKyRkTmi0iKiIwRkcUisllE/iUiru6evyDDxWNXH0mqM7KJMGML0rnrS5NIc/WPYOi0CTlpnaeXE4HsVE1Bp5RSnYlbMBSRYcC3AI8xZhpgBy4Ffgvca4wZD1QD1/XgGowpTOPFW4+hsItE3LNG5jD/hrkMykrp7uX6nJw0F187elSn+5wwoYD0fhL8lVIqXuI9TOoAUkXEAaQBe4BTgOfC7z8OnN+TCzjtdiYOyuS/tx/HP75+JLNG5NC6tM5pF86aPphXv3Ucj3zNQ1E/CoQALoeNq+aOZvLgzHbfL8hw8avzp5OlPUOllOqUxDNPp4jcDvw/oBl4E7gd+DTcK0RERgCvhXuOBx97A3ADwMiRI2dv3749omtWNfoIBEOEAJtAmtNORj+vVFFe7+X5ZaU8+uHnlDd4SXfZueyokXzj+LEUZenCe6USRP+jJbG4BUMRyQWeBy4BaoBnsXqEd0cSDNvyeDymuLg4Lu3sLwLBENVNfoLGYBPITnHijvBZqlIqJjQYJrF4Pkw6FfjcGFMOICIvAMcCOSLiMMYEgOHArji2YcBw2G0UZmoBY6WU6o54PjPcAcwVkTSxxunmAeuAd4ELw/tcDbwUxzYopZRSXYpbz9AYs1hEngOWAQFgOfAw8CrwjIj8Krzt0Thd33p+GDL4AiHAmnBitwm5aS7s/ShJt1JKqZ6J6wSaWInmmWEgGKKq0cfCDWU8+el2Nu6t378w3W4TJgzK4GLPCM6aPoS8dBcuR1LnHVBK9R36CTuJ9atgWNPkY/HWKu58diX1XaQgS3Xa+dX50zhl8iBy07u97l8ppVppMExi/aZbVN3k4963N3Hjk0u7DIQAzf4g3312JXf/Zy1Vjb4EtFAppVRf1S+CYaM3wJOLtvP4J5GtRWzrpRW7uW9hCXXN3auNqJRSKvn1i2BYWt3EvW9v6vbxj3+yjc1lDTFskVJKqWSS9MGwpsnHT19aSwTFGzr1o3+vplLr/iml1ICU9MGwwRtg8edVPT7P+j31VDfpUKlSSg1ESR8MV5fWxuxcn26tjNm5lFJKJY+kD4aLYhjAPt1aiVcL4Sql1ICT1MEwEAxRE8OhzYaWAP6ePnxUSimVdJI6GNptQlZq7DLKpbrsOO26blYppQaapA6GIsKcMfkxO99RY/JwO7TskVJKDTRJHQwBDh+RE7NzHTMudoFVKaVU8kj6YJjhdnDEyJ4HxHGFGeSnaz1ApZQaiJI+GOalu/j5edN6fJ5fnncYBVocVymlBqSkD4YAY/LTuf74Md0+/qtHDOOwoVkxbJFSSqlk0i+CYUaKg5tPGs/5M4dGfezJkwr5n7OmkJ2mZZyUUmqg6hfBEKzh0p+dcxi//ep03BEU7HXahZ+cPYU/Xnw4+Rk6PKqUUgNZ7Bbp9QG56S4umDWcEyYW8sKyXcz/bAel1c0H7DMkO4ULZw/n0iNHkpfhJNXZr/4IlFJKdUO/qnTfViAUorrRRyBk8PpDGCDFacNhE7JTXbgi6D0qpVQUNGNHEuu33SKHzUZhZkpvN0MppVQS0O6RUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAS9uwVBEJonIijZfdSLybRHJE5G3RKQk/D03Xm1QSimlIhG3YGiM2WiMmWmMmQnMBpqAF4EfAAuNMROAheHXSSkQChEMhXq7GUoppXooURlo5gFbjDHbReQ84KTw9seB94DvJ6gdMVFe72XNrlpeXrkbgPNmDuWwoVma8UYppZJUooLhpcD88M9Fxpg94Z/3AkXtHSAiNwA3AIwcOTLuDYzUtopGLn34U/bWtezf9uLyXQzNTuGZG+YyMj+9F1unlFKqO+I+gUZEXMC5wLMHv2esLOHtZgo3xjxsjPEYYzyFhYVxbmVkyupbuOz/DgyErXbXtnD5I4spr/f2QsuUUkr1RCJmk54JLDPG7Au/3iciQwDC38sS0IaYWL+7jj21hwbCVqXVzZTsq09gi5RSSsVCIoLhZXwxRArwMnB1+OergZcS0IYeC4UM/1m1p8v9Xlm1OwGtUUopFUtxDYYikg6cBrzQZvM9wGkiUgKcGn6dFGwRVCsT0ZJmSimVbOI6gcYY0wjkH7StEmt2aVKx2YTzZw5jQXFpp/udN3NYglqkVOzVtfhpaAmwr64FERiclUqG205GirO3m6ZUXPXb4r7xMKEog5F5aeyoajpgu8MmnDx5ECdPGsTw3FR2VDXh9QcJhgxOu41Ulx2HXUh3OUh36x+56nsCwRCl1c388pV1vLuxjFB4WpvDJnzpsMH88KzJDM1OxRbJ8IhSSUisCZ19m8fjMcXFxb3dDABKq5q48tHFbKtsoiDDxRVzRnHypELe21TOB5vKWbenjhb/oQvxCzPdzBiezTXHjGbS4Ezy0l04bJoNT/U+YwwlZQ185YFPaPAG2t0nN83Jv289llG6dKgz+kkhiWkw7Iby+haqGn0EQob739nMW+v2EQxF/udYkOHihhPG8pUjhlOQ4Y5jS5XqWnm9l4se/IRtlU2d7jdjeDZ/v+ZI8vXfbEc0GCYx7Zp0g80mLNlWxVce+ITX1+yNKhACVDT4+PV/N3DJQ4vYsLeOxg4+jSuVCHtqm7sMhACrSmupafInoEVKJZ4GwyiV13u57enl/Pjfa/EGepaXdEt5I2fd9yGvrt5DXbPeZFTveGVl10uGWr23qTyOLVGq9+hsjihUNHj5xj+XsHJnbczOGTJw13Or8PpDnD9rKJk6a08lWJM/8pGJJh3FUP3UgAmGvkCQ2uYAIWMIGYMANhFSnHayUrsOQNWNPu58dmVMA2FbP3lpDSPyUjl+fAF2u3bYVeIcNjQ74n0nDc6MY0uU6j39OhgGgiGqmnxs3FPPc8tKWVVay7bKRlrnDOWkOZk2NJuTJhVy1vQhZKY42u2Z+YNBPiwp572N8R0i+tYzy3nz2ycwODs1rtdRqq2TJw3CaRf8wc6ffWe4HcwckZOgVimVWP12Nmllg5dXVu3hzwtLqGz07d+e6rSTlerAGwhR2+zfHxjtNuH0qUX8+MtTKMx043LY9x+zu6aZU//3fZp8wZj8Pp0547DB3PPV6eSkueJ+LaUAGr0B/vJOCX97f2un+/3oy1P42txRuJ32TvcbwHQ2aRLrdz3DQDBEaU0ztz61jLW76wArjdopkwdx+ZxRuOw2Khu9pDjtDMp088bavSwoLqWq0cdra/by/qZyfvOV6Zw4sZCcNBdef5BHP/o8IYEQ4PW1e7nrjEkaDFXCpLsdXH/COBq8QZ74dHu7+3zzlPFcOHu4BkLVb/WrnmEgGGLjvnoueejT/YuH89Nd3H/5LIq3VTP/sx0HVJ1w2W2cMW0w1xwzmr++t5mF678ooHHjiWO5+cRx+AIhTv7DezQmKBgCXH7UCH5y9mGkuvTGoxKnpslHdZOff3z8OUt3VCMIx4zL58q5o8hOdUb0bH2A055hEutXPcMdVU1c/OCi/YEr0+3gwatm8+tX17N8Z80h+/uCIV5euZu31u3jgSuOIBg0+6eOP/T+Vsbmp5OZ6kxoIAR4ftkubp83UYOhSqicNBc5aS5+eNYUmnwBQMhwO3A5dEKX6v/6zb/yygYvNz6x9IDAdcfpE3nw/S3tBsK2mv1Bbnt6Gd8+dQLuNv/x1+yu4821e+PW5o54A9bEH6V6Q4rTTl66m7x0lwZCNWD0i3/p3kCQBcWllJQ17N+W5rIzfVj2AUOfnWn0BfnPqj2cPWPI/m3Th2ezsjQ+Sym6snxH5wFcKaVU7PSLYFjb5Oev724+YNuXZwzh5ZXRFdp9dulOLpg1fP/rUXlpbKtsjEkbo7VoawW+QOfDsw0tfsrqWiira6E5wUO5SinVn/SLZ4ZLtlUfkm1/XEEGb0Q5xFnXHMBp/+IZuAF6a35RbZMfXyB0wBKPVg0tfnZWN3P/OyUs3lqFzSacNqWIG08cS1FWCik6408ppaKS9MGwvsXPv5bsOGS7wy4EokygfbBETw3LTXMybVg2g7NSGD8oo93acfUtfl5dtYcfvLD6gO1Pf7aDZ5fu5B9fPwrPqFydAq+UUlFI+mHSZn+Q1bsOfa5XXu9lWG50mVzsNsHeNgAlIBo6bMKFs4fzzA1z+d2FM5g10srwsbm8gQVLdrJmVy1l9S0EglZS8KpGHz98cXW75/IHDd94vJgaTfqtlFJRSfqeoT9oqG6nrMwrq/Zw97mH8fqayIdKT50yiA9LKva/rmr0UZDhoqIhPjM7pwzJ5O5zD+Od9WXc8EQxdc0HDvU+G/5emOHmlpPHcd7Mofxz0fZOh26b/VbquAtnj4hLm5VSqj9K+p5hsIN8irtqmjHGMCo/LeJzXTFnFPM/+2LIdXVpLdOGRZ7EOBonTSzkp2cfxncXrOShD7YeEgjbKm/w8vP/rOONtftYtr26y3N/vLmSUA+HiJVSaiBJ+mDY2VDmvW9v4vcXHk6mu+sO8C0njWNVaS1l9d7921bvqmVWHBITzxiezfUnjOW6x5dQWt0c8XFl9V5SIliIn5niaPd5o1JKqfYlfTB0O2x0dN9fv6eee9/exGNfP7LDbPt56S5+ds5UBmen8Me3Nh7w3uKtVZwwsRCJYVxxO2z89Oyp3LFgRdT5Tj/YVM5Z04Z0ud/FHh0iVUqpaCT9M0O7TRhTkMGW8oZ231+0pZLv167iG8eP5YdnTubNdfsor7cSdc8dm8egzBTmf7aDV1cfWu3bFwyxaEslJ04ojFmF78vnjOTfy3exr87b9c4HWbq9mru+NInhuakd9ihnj8plSHZKT5uplFIDStL3DLNTnMwZm9fh+0OyU8hJc/Lc0lJ+9vIa9tQ2k5niIBAK8dD7W7ny0cXtBsJWT3+2g+uOHxOz9p49YyjPLSvt9vH/8+Ia7r9sFocNzTrkvRMmFvDglUeQn+HuSROVUmrASfqeodNh45qjR/P04i8mvhw2NIurjh7FhEGZ7K5pprS6CWNgaE4qw3NT2VndzBOLtrNxX32X5y+tbqbZF+RrR4/in4vaL28TqXGFGWwua6DFH+r2ObaUN3D7Myu447SJTCjKYOn2apw2G0ePzyc71Umuln5SSqmoJX0wBMjPcDF1SBa7a5v52TmHYYzh7x9va3f9IcCkokyuPmYUN5wwlp+8tIby+s6HLMcPymDWyFzeWrfvgBJQ0Zo+LJtVpT3PObqjqolv/2sFv/nKNL529Ogen08ppQa6pB8mBcjPcPPny2by92uO5NnindyxYGWHgRBg4756/ufFNTz8wRYevmo20ztZPnHFUSMpyHBTmOnmyevmkJXS/c8PQ3NS2FnV1O3jD6bJvJVSKjbiGgxFJEdEnhORDSKyXkSOFpE8EXlLRErC33Njca1Up51vzl/OJ1sqIz5m2Y4abvjnUn52zlTGFWYc8n5RlpvbT5uwv6jp6Pw0nrvpGHLSulfkVESI5eq/kIFkKM6slFJ9Xbx7hvcBrxtjJgOHA+uBHwALjTETgIXh1z1S2eDl+n8ujWrNXqvyBi93LFjJry+YdkAqtky3g6e+MYfCNpNR7HYb4wdl8Oo3j2fOmI4n7XSkrK6FIdnRpYjrzMSiTCSW6z6UUmqAilswFJFs4ATgUQBjjM8YUwOcBzwe3u1x4PyeXMcXCPHvFbtZt6eu2+fYUdXEG2v3csWckYCV/uyFW45hTH76IcHGZhOG5abytyuP4I8XHR7VsOma3bWdDslG68jRMelUK6XUgBfPCTRjgHLg7yJyOLAUuB0oMsa0rmXYCxT15CI1TT7+vLCkRw0FeGrxDp6+fi6N3gDfP2MyBRnuTrO45KW7OefwIRw/oYDlO6v588LNrN3deUDeVtGIZ3QuDlvPK2pkpToYlhO7XqZSSg1k8QyGDuAI4JvGmMUich8HDYkaY4yItBsVROQG4AaAkSNHdniRkrIGamNQpcEbCLGzqomfnzuNjAh7ey6HnUFZdr502BA8o/LwBUNUN/pYsbOWXTVN+AIh0t0Opg7JYkJRJmkuOy67jXNnDuWFZbt61N6vHzNGl1EopVSMxDMYlgKlxpjF4dfPYQXDfSIyxBizR0SGAGXtHWyMeRh4GMDj8bQbMEMhw3sb2z28Wz4oKees6V2nO2tP60L3IdmpTB3a+VDoXV+axML1Zd0O4sNzU7ly7iicjn4xGVgppXpd3O6mxpi9wE4RmRTeNA9YB7wMXB3edjXwUnev0eALdDk0GY2Ne+upb4l/LcDCDDePXO3B0Y1k2mkuO49dfSSFmZplRimlYiXeXYtvAk+JyCpgJvBr4B7gNBEpAU4Nv+6WUMjgC3Q/m8vBvIFQTJc+dMRutzFtaBbzb5hLdmrkyzSKsty8cPMxjClMj2PrlFJq4IlrBhpjzArA085b82JxfpfD1u01f+3JSXV2q7fWHakuB0eMyOGNbx/P717fyMsrd3c4qcbtsHH5USO5+eRxFGa4dTmFUkrFWFKnY0tzOZgzNp+318fmueHsUblR9dR6ym63MTg7lV+eP427zpjEsh01fFRSwY6qJmwCYwszOG5CAdOHZZOV4iQ1glqGSimlopfUwRDgxImF/L9X18fkXF+aNrhXel3pbgfpbgdnTU/l9KlFNPuDCFaw1yK9SikVf0k/HTEv3cXkwZk9Ps+wnFRG5Pb+uj2H3UZmipOMFKcGQqWUSpCkD4YFGW5++9UZPT7PPV+dTn66ztBUSqmBKOmDIcDYwnRuPGFst4+/aPZwZgzP1p6YUkoNUP0iGGamOLnpxHFcflTHmWo6ctb0wfzwrMlkp2o2F6WUGqj6RTAEyE13cdcZk3joytkRJc9Oddr540WH86vzp5Onw6NKKTWgSTLUw/N4PKa4uDiifQPBEFWNPt7fVM5Ti3ewfk8d3vDCfKddmDQ4k0uPHMFpUweTm+bCpSnNlFKxoc9ZkljSL604mMNuY1BWChfOHs68KUX4gyH8wdZgaMNhE3LTXPp8UCml1H79Lhi2EhHy0vU5oFJKqa7pGKFSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg6FSSqkBT4OhUkqpAU+DoVJKqQFPg2GSqmr0sbummR2VjeyoamJvbTO1Tb7ebpZSqp8SkfdExNPb7TiYiJwkIq/09Dz9trhvf+QPBqlu9LN6Vy33LSxhVWnt/vdE4LjxBXzrlAmMLkgnP92FzSa92FqllLKIiMMYE+jtdnRGg2Enapv9+IMhbCJkpThw2HuvI93oDbB8Rw23zV9GTZP/kPeNgQ9LKviwpIKh2Sk8es2RjCtMx+Ww90JrlVK9RURGA68BHwHHALuA88Lb7jTGFItIAVBsjBktItcA5wPpwATgD4ALuArwAmcZY6rCp79KRB7Bih3XGmM+E5F04H5gGuAE7jbGvBQ+71eADMAOnHhQO/8GHAmkAs8ZY34W3r4NWACcCTQDlxtjNovIP4AWwANkAXcYY15pcz4bsBE4xhhTHn69CTjaGFPe1Z9bl8FQRH5rjPl+V9s6OHYbUA8EgYAxxiMiecC/gNHANuBiY0x1V+dKlNbe14qdNfzjk21UNHhJddqZN6WIC2cPJzPFQWaKM6FtavIFWLihjNufWY4xXe9vswkflpST7rbjstux2yAvzYW9F4O5UiqhJgCXGWOuF5EFwFe72H8aMAtIATYD3zfGzBKRe4GvAX8K75dmjJkpIicAj4WP+xHwjjHmWhHJAT4TkbfD+x8BzGgTTNv6kTGmSkTswEIRmWGMWRV+r9YYM11EWq99dnj7aOAoYBzwroiMbz2ZMSYkIk8CV4SPORVYGUkghMh6hqcBBwe+M9vZ1pGTjTEVbV7/AFhojLlHRH4Qfh3pueKqyRtg+c4abn5qKXXNB/boV5bW8qe3N3HNMaO57ZQJ5KW7EtauDXvr+XYEgdBlt/GL8w8jP93Fk5/u4J7XNhAyMCjTzXXHjeH8WcMozHDr8KlS/d/nxpgV4Z+XYgWRzrxrjKkH6kWkFvhPePtqYEab/eYDGGM+EJGscPA7HThXRO4M75MCjAz//FYHgRDgYhG5ASsODQGmAq3BcH6b7/e2OWaBMSYElIjIVmDyQed8DHgJKxheC/y9i997vw6DoYjcDNwCjBWRVW3eygQ+jvQC7TgPOCn88+PAe/SBYOgPBlm+s4arHl1MqIOgEzLw2MfbqG328+MvTyU3AQGxssHLT/69psM2tbIJ3HfZTN5et4/nl+064L2yei+/eW0DD76/hQU3Hs24wgwNiEr1b942PwexhiIDfDFpMqWT/UNtXoc4ME4cfCcygABfNcZsbPuGiMwBGttrnIiMAe4EjjTGVIeHQNu2yUTw8yGvjTE7RWSfiJyC1YO8or3rt6ezcbOngXOAl8PfW79mG2OujPD8BnhTRJaGPwEAFBlj9oR/3gsUtXegiNwgIsUiUlxeHlEvt0eqG/3c/NTSLoMOwPPLdrG5rCHubQKobvSxdnddl/udP2sYG/fWHxIIDzhXk59LHv6U8gZvh/sopfqtbcDs8M8XdvMclwCIyHFYQ5m1wBvAN0VEwu/NiuA8WViBslZEirBGGw+5Tvj7ojbbLxIRm4iMA8ZiPSM82CPAk8CzxphgZL9W58HQGGO2AbdiPfdr/SL83C8SxxljjsD6RW8NjzMfcAEOjfSt7z1sjPEYYzyFhYURXq77VuysOWRotDP3LdxEbTsTWWKp2Rfg4Q+3RrTvRbNH8PePt3W5X1Wjj3c2lGEiefiolOpP/gDcLCLLgYJunqMlfPyDwHXhbb/EmjizSkTWhl+3S0RWABhjVgLLgQ1YHa+DRxtzwyOStwPfabN9B/AZ1mSgm4wxLe1c5mWsSTsRD5ECSEc3RRF5xRhztoh8zhdd4VbGGDM2qguJ3A00ANcDJxlj9ojIEOA9Y8ykzo71eDymuLg4mstFpbbZz01PLmXRlsqojvvkB6cwNCc1Tq2C8novFz34Cdsqmzrdb3huKnd9aRLfemZFp/u1GleYwb9umEtBpjsWzVRKWfTZQwyEJ156DpprQngo9RVjzHNdHO8B7jXGHB/NdTvsGRpjzg5/H2OMGRv+3vrVZSAUkXQRyWz9Gesh6xqsqH11eLersR529ip/MERlN4YO/cFQHFrzBWMM9S1d91YLMtzsrm3vA1L79tY2E9KeoVKqnwlPynwe+GG0x0a0zlBEZmDNRtq/vzHmhS4OKwJeDA8jO4CnjTGvi8gSYIGIXAdsBy6OttGxZhMh1Rn9ejybxPmDoEBKBO1q9gVJc0Xe/jSXLi9VSvVNxpjRHWy/JoJj7wHu6c51I1ln+BjW1Nq1WDOLwBo27TQYGmO2Aoe3s70SmBd1S+MoK8XBvClFrGyT0aUrI/JSSXHGd92e225jRF4qu2qaO91va0UD04dlYxMimgD05emDyUpN7FpJpZTqyyK5m88NT2S52hjz9fDXtXFvWQI57DYunD2caFYb3HziOAoy4vvMLTvNxU0njutyP3/Q8GFJBfOmtDsx9xDXHDsmoh6nUkoNFJEEw0UiMjXuLellmSkOrjlmdET7Ds9N5bSpRUi8h0mBqUOzyE3ruhf3z0XbuPWkcQzP7XxCz8/PPYyCjMQlDFCJU9Pko7zeS2WDN+7Ps5XqbyJ5ePRPrIC4F2shpmDNJp3R+WHJJTPFyW2nTKC22d/pWr3huanMv34uhZkHr1mNj9xw7/A3r23odL+KBh93PreK+y+bxT8+2cZrq/fia3NDHFuQzo++PIUjR+eRkeB0cip+fIEgNU1+irdX849PtlHZ4CXFaefUKUVc7BlBZqqDLP377rM8Hs844A7gSqyEJvVYa+T+t7i4eEtvtm2g6XBpxf4dRDZj/WWt5otnhhgSn393AAAgAElEQVRjtse3aV+I99KKtqobfWwua+C+hZv4eEvl/hRoI/JSueWk8Zw6pYjCBC9JqGr08YPnV/Hmun1d7pud6uTSI0dwxrTBVDX6sIkwLDeVvHRX3Id1VWI1+QKs2FHDTU8upa6dWccicO2xY7j15PEJTR84gEU1VOTxeM4EnsNao9f2E4s//HVhcXHxa91qiMgZwH1YCbIfCU8safu+G6ujMxuoBC4JrysfsCIJhouMMUcnqD3tSmQwbFXb5KfJH8AfsKpWuJ02CjLcCRkabU9lg5dfvbqeF5d33Gtt6zunTuTqY0eRk6o3wf4oEAxRvL2ay//v0y4nTV125AjuOnMyuWn6byHOIr45hHuEq4C0TnZrAmZE20MMJ77ehJVXuhRYgpW0e12bfW7BSqB9k4hcClxgjLmk3RMOEJEMky4XkaexErfuX4wXwdKKpJad5iSbvjO8lJ/h5qdnT+WquaP48zslvL+p/JDE3U67cPb0odx00jiGZKfojNF+rKrJx21PL4to9vD8JTu5Yu4oDYZ9yx3Q5Q3GiZV95bYoz30UsDk8ox8ReQYrJ/S6NvucB9wd/vk54C8iImYAp6aKJBimYgXB09ts63JphYq93HQXueku7rt0Fg0t1nOi8novDrtQlJXCrBG5ZKY4SHfrOsL+btPeBioafBHv/5d3NvP7i2YkvPyY6tCVRBYMryL6YDgM2NnmdSkwp6N9jDGBcKWKfKCCAarLu6Yx5uuJaIiKXHaqk+xUJ8NyOxthUf1Vsy/AE59G98j+rfX7aPIFNRj2HZkR7pcR11ao/Tor4XQ/HSTRBjDGfCsuLVJKdcoXNNQ0Rd4rBAiGDMFIxlRVotRjVW7oSnfK4+wCRrR5PTy8rb19SkXEAWRjTaQZsDpbZ1iMVRQyBatacUn4ayagDx+U6iUOm3RrKFxrWPYpT2LNGO2MH3iiG+deAkwQkTEi4gIuxcoJ3VbbHNEXYlWqH9Cfljr8H2WMeRz2F/k9zhgTCL9+EPgwMc1TSh0s3e3gglnDeGdDWcTHzB2bR4ojvukDVVT+F7iGzp8b+jmwyntEws8Ab8OqM2gHHjPGrBWRXwDFxpiXgUeBJ8JL56qwAuaAFsnHy1ys7nxV+HVGeJtSqpfMGZNHustOoy+y2qW3z5tIjs4m7TOKi4u3eDyeC+l6nWG3Ft4bY/4L/PegbT9t83MLcFF3zt1fRfJR8R6s5RX/EJHHgWXAr+PbLKVUZ7LTnPzivGkR7Xvc+AImFuk8jL4mvKB+BvAUVtX3UPj7k1jrC7u14F51T5eL7gFEZDBfTM1dbIzZG9dWHaQ3Ft0r1dfVNvt5ZdVufvzvNYesOW11woQC/veSmZp9KDGiWXRvB74M3IWVBaZtQVU31nyN3wGvFhcXR9b9Vz0SaTDMBSZgTaYBwBjzQRzbdQANhkq1r8EboKbJxzNLdvKvz3ZS0eglxWHn+AkF3HbKeIbnppKXroEwQSIKhh6P50Ss4VE3nS+xqAdagIuKi4vf73nzVGciScf2DeB2rOm5K4C5wCJjzCnxb55Fg6FSnfMFQtZyCwEMpLrsuqYw8boMhuHnhP/ESmYSqWbgquLi4ue72zDVtUieGd4OHAlsN8acDMwCauLaKqVUVFwOG4OyUhiUmcKgrBQNhH1QuEf4ONEFQsL7PxE+XsVJJMGwJTzzCBFxG2M2AJPi2yyllOo/ws8In6PzxNydSQWeDZ9HxUEkwbBURHKAfwNvichLQMLKNymlVD/wZaxnhD2RApzV1U4ikiIin4nIShFZKyI/D28fIyKLRWSziPwrvCAfEXGHX28Ovz+6g/OeISIbw/v9oIe/S5/TZTA0xlxgjKkxxtwN/ARrseb58W6YUkr1I3cReT7SjmSGz9MVL3CKMeZwrIxhZ4jIXOC3wL3GmPFANXBdeP/rgOrw9nvD+x0gXBbqr8CZwFTgMhGZ2sPfp0+JKiWFMeZ9Y8zLxpjoEiMqpdQA5fF4hmEtn4jR6TxDO9vBWFpzmrYu6DfAKVhDtWA9u2zt1JwXfk34/XlyaOHW/WWhwvf/1rJQ/Ua38jOJyCuxbohSSvVTIzhwHWFPeDkwCXe7RMQuIiuAMuAtYAtQ05pWE6us07DwzweUcwJayzm11V5ZqGH0I91NVnh9TFuhlFL9V0+HR6M+nzEmaIyZibUk7ihgcozb0O90KxgaY/bEuiFKKdVP1ffW+YwxNcC7wNFATrhcExxY1ml/yadOyjlFUhYqqXVWz3A17dczFKxh6Rlxa5VSSvUfO+n5TNJWbg4crjyEiBQCfmNMjYikAqdhTYp5F6tc0zNY5ZteCh/SWs5pER2Xc9pfFgorCF4KXB6T36iP6KxqxdkJa4VSSvVTxcXFuzwez1Lg2Nicrnh3F/sMAR4PzwC1AQuMMa+IyDrgGRH5FbAca2UAdFDOSUSGAo8YY87qqCxUDH6fPqOzeoa6llAppWLjd8B8ur/oHqApfJ5OGWNWYWUKO3j7Vqznhwdvb7eckzFmN23WNbZXFqo/6fKZoYjMFZElItIgIj4RCYpIXaQXCM9qWt46A7WjhZ9KKdUfeTwewUpp2ZNASPj4o8LnUzEWyQSavwCXASVYKYG+gbX4MlK3A+vbvO5o4adSSvVHPwe+E6Nz3QH8LEbnUm1ENJvUGLMZsIen6/4dOCOS40RkOFYaokfCr4WOF34qpVS/4vF4zgS+C6TH6JRpwPc8Hk9E92AVuUiCYVN4KHOFiPxORL4T4XEAf8JKHxQKv86n44WfBxCRG0SkWESKy8vLI7ycUkr1DR6PJweran1Ph0cPlgY86fF4smN83gEtkqB2VXi/24BGrLUmX+nqIBE5GygzxiztTsOMMQ8bYzzGGE9hYWF3TqGUUr3pR8Q+ELZKB34cp3MPSJEEw/ONMS3GmDpjzM+NMXcQ2bKLY4FzRWQb1rqWU4D76Hjhp1KqnwiFDNVNPsrrWyivb6GiwUuzL9jbzUoYj8eTAtyIVWkiHlKAGz0eT6zWLw54na0zbHU1VhBr65p2th3AGPND4IcAInIScKcx5goReZb2F34qpZJYXbOfRm+ADXvr+XhzBct21FDV6CVoDCkOO2MK0jl2fAFHjs6lINNNfpoLu727GSH7vEQlsT4f+NfBG8Nl9x4BpmElT7kW2BjedzSwDbjYGFMdnstxH9YyiibgGmPMsnbOORv4B9ZEyv8Ct7ezOD9pdZaB5jKsDANjROTlNm9lYS3M7K7v0/7CT6VUEqpq9LGjspH739nMuxvLCHVweywpa+DNdfsAKMx0843jxnDezGHkpbtwOfpdUJxH7HOSHiwTOJl2giFWcHvdGHNheM5HGvA/wEJjzD3heoQ/wLofnwlMCH/NAf4W/n6wv2HlpV6MFQzPAF6L6W/Ui6SjwC4io4AxwG+w/tBa1QOr2kyCiTuPx2OKi4sTdTmlohYKGZr9QRw2we0cGMXIfYEg5fU+vvfcSj7ZcnAqy8ikOG388rxpzJtSRF567y45rmjw4vWHCBmDy2EjNy3qIL1//Z/H41mLVfcv3tYWFxdPO6ARItnACmBs256biGwETjLG7BGRIcB7xphJIvJQ+Of5B+/X5tghwLvGmMnh15eF97kx3r9gonSVgWY7cLSIFGEtGgVYn8hAmEyafAHqWwJ8sKmcdXvqaPQGKMhwc/LkQYzOTyMv3Y3dputl+5P6Fj+1zX6eLS5lVWkt6W47V8wZxYRBGRRk9t/HOXXNfj7aXMGdz66kqQfPAlv8Ib733CqOGp3LfZfNYkh2agxbGZmqRh8flZRz/zubKSmzygBmpTi4fM5Ivn7MGAZluTm0vF+XRsa8oe0b1c62MUA58HcRORxYirXeu6hNgNsLFIV/7qg8U9uCDMPC2w/ep9/o8pmhiFwE/AF4D+uTz/0i8j1jzHOdHjiAeANB9tS0cP+7m/nPit34gqED3n/gvS0MyU7hxhPGcvaMof36JjmQ1Db5eHLxDv7w5kbaDrC8smoPx08o4N6LZ4JAMGQIhscO7TbBJkJ2qjNphwbrmv08v6yUn/9nXczO+dm2ai56cBH/umEuw3LjNQHzUFWNPn704mpeW7P3gO11LQEefH8rzy/bxfM3HcPI/KjbFMl8jFho7zoO4Ajgm8aYxSJyHweO7mGMMSLSb573xUIkf2E/Bo40xpTB/ozob/PFwvkBrckXYMm2aq5/vPiQINjWntoW7v7POuZ/tpPHrz2KwdnxmmSmEiEUMry3qZzfv7Fx/za7TTh9ahGXHTWS2mY/97y+gVWlNWwpb9wfDEVgTH46M0bkcM6MIUwflk1O9MNxvabJF+D1tXtjGghblVY3c8nDn/L8zcdQlBX//x/BkOGVlbsPCYRtldd7uf6fxTz1jTnRfohN1OhZe9cpBUqNMYvDr5/DCob7RGRIm2HSsvD7kZRn2hXe3tk+SS2S/4G21kAYVhnhcf2ePxhkza5arv3Hkk4DYVsb99VzxSOLKatviXPrVDxVNvoOCITjB2Xw1DfmMLEok+8uWMk35y/nuaWlbNrXsD8QAhgDWysa+ffyXVz3eDGn3vs+979TQll9C8kwMW93TTM/fGF13M5fWt3Mj15cTXWTL27XaFXZ6OVv72/pcr+N++opq4+6UP2ObjUqeocUVDDG7AV2isik8KZ5wDq+KNUEh5Zw+ppY5gK1B9esDb+uC+eqFuBr9LOVAJEEtddE5A0RuUZErgFepR9nLo9GVaOfa/9RfMDNLhJbyhv409slA2rdVX/T6A1QWt0MwDXHjObucw7jB8+v4r6FJZQ3RH7jrGsOcP87mznn/o8o3l5NXbM/Xk3uscoGL7c8tSzqf+/Rent9GUu2VRGI8ANmdzV5g+ypjexD6Usrou4EfRx1g7rnow62fxN4SkRWATOBXwP3AKeJSAlwavg1WPfzrcBm4P+AW1pPIiIr2pzzFqzlGpuBLfSjmaQQWTA0wEPAjPDXw3FtURL57PMqGrzdGw15YVkpdS1998anOtc6EnDn6ZMYlZ/G1x5bzLbKpm6fb1+dl4seXMQTn26nJgG9omiFQoa31u1j076GhFzvzmdXUtUY3z+HaIJtffT/VxcC8f6L9GIV7D2EMWZFOIPXDGPM+caYamNMpTFmnjFmgjHmVGNMVXhfY4y51Rgzzhgz3RhT3OY8M9v8XGyMmRbe77b+tMYQIguGpxljXjDG3BH+ehFrXcqAVtHg5f53Srp9fIs/xPubypNiaEwdKtPt4MYTxpLqsvHz/6zrcG1dtH7/xkYe/2QbtX0sIFY2+vjLu5sTdr265gBLtlfH9RrpbgeRThKdOjQr2tO7wl/x5Aaccb7GgNFhMBSRm0VkNTBJRFa1+focWJW4JvZNzb5gjz8lL1iyk9o+PCymOpaZ4uC4CQX86tX1Xe8cpXvfLmFxAoYJo7Grumn/sHCi/OWdEiqjGHKOVprLwXHjC7rcz2ETTp1S1OV+rTwez2CiK3PXEw94PJ7IG6c61FnP8GngHKyHq+e0+ZptjLkyAW3r0/wxuFFVN/nwB7VnmIya/UG+/9wq4tWxv3PByqiePcZTIBhiQXFp1zvG2Po99TT7Y/9cvbLBS2l1E5vL6vnuaRNx2jvvHl5/wlgyUqJaKfF3rF5bIrixUqSpHuowGBpjao0x24wxlxljtrf56kkqtn4j+jW4h7KJxOQ8fUWjN0B5fQv76looq2+hsZvPU/u6Fn+QRz76nN0RTr7ojrqWAL9+dX2feK5c2+xn+c74Dll2JJa90coGL2+t28vVf/+M4377Ll99cBGPL9rO/ZcdQYb70GAnAtceO5objh9Lhjuy0UiPxzMDOJH4D5G2cgEneDye6Qm6Xr+VqIWh/Y7L0fOUW4WZbpxJnqi42RekrsXP+5vKeX3NXtaHM++kuRxMHpLJ6VOLOGVyEVkpDtLaueEko7pmP08sOmRGe8y9unoPPzxrClkpvftYKBAylCRo4szBPt1aydyx+T06hy9gzRq97enlrN5Ve8B7Ly7fRU2Tn8euOZKNe+v5oKQcfyDE+KIMzpkxlLx0F+nR/bv9Dol/jucGvg1cl+Dr9iv94+7UC9KcdmaPymVpDx7yf/3Y0WSnJu/z74oGL898tpP73ynBGzhw2LiuJcDeuhbe21iO076Wm04cx9XHjKYgI/mz73yypbJHKcgiFTLWrOObThyHoxc/NPkCIQJxXk7Rkc1lDfgCoW4nJfD6g6zZXcsVjyymxd/+o413N5bx7sYyjhiZyxGjcnDZbeysauLihxYB8Pi1RzFzRA4pXeSc9Xg8GcAlJP6+agcu83g8txcXF/fOp5Z+ILm7Jb0oN93Ft0+d0O3js1IdHD48J4YtSqy9tc1c9ehi/vDmxkMC4cH8QcP972zmkocWsbsmsZMwYq2u2c+C4p1d7xgjLy7fTU1T7w6VhnpxxrMvEOrRusbPKxq5/P86DoRtLdtRzSMffs4D723hP6v24A2E8AZCXPXo4kh7xkcBvfWX5eeL/NGIyO0iskZE1orIt8Pb8kTkLREpCX/PDW8XEfmziGwOT5I8or0LiMhsEVkd3u/P0o2ErX2ZBsMemDQ4k6Ks7vV0rj12DDlpvZulv7vK671c8chi1u+pj+q4LeWNXPzQIsrqkjf7Tku4p5Eon1c0xH2Re1dsvXjPc9ptdLdTXF7v5fonirv8sNYVf9BwwxPFkWSN8hC/Yr5dSQlfHxGZhlVq6SjgcOBsERmPlZJtoTFmAtY6yNZ8pW1LON2AVaqpPa0lnFr3PSMuv0kv0WDYA4UZbp68bg5pruieHx49Lp+vHT0qafJRtlXX7OePb25kS3ljt44vrW7mV6+u75MLyyPhC4aoa07cxKCQgXpv7/YMnQ4bvVVsZWReWreeq/sCIZ78dDs7q2IzErGntoVHPvicls5ntyZy4szBXOHrA0wBFhtjmsIVht4HvoJVcPjx8D6PYxUGJrz9n+HF958COeHcpfuFX2cZYz4NL7b/Z5vj+4Xkuxv3ISLC6II0Xrj5mIhrsc2bMoj7L5tFXnpyPjsrq2vhmSU9GyZ8eeVudiV4zVqs+HrYy+iO2gQG3/bYRRhTkNEr1z5mfH53yidR0+TjsY8+j2lbnvh0e1frgnu7pFFrIu01wPEiki8iaVgV7EcQfQmntrSEk+qc025nYlEmr91+PP9dvYeHP9jabr7Do8fmc/upE5hYlNGjQBgMGaoavTT5gmwua2BXTTMCjClIZ1R+OqlOe9xKRDV5A/ztg60xOddf3t3M7y6cQWYvz5RMBqFeHibNSXMyc0Q2W8oTPzdjdH56t45bt6eO+hgv7Wn2B1m2vZozpw/paJfefu7hAjDGrBeR3wJvAo1YhX4P6NJqCadDaTCMAZtNKMpK4WtHj+bsGUPYXdPC9spGmnxBctJcTB2aRVaKo0fPCI0xlNd7eWnFbh75aCv76tpfkD22IJ1bTx7PiZMKuz1zs7LBS9AYMJDitJMVnvHa4AvwZgflbpx24fSpgzlz2mCyUp3UtwR4Y+1eXl+zt92KHu9sKKPZF4woGNY2+2nyBfAHDU67kO5y7G9TovVGcebUKIfhY81pt/HV2cN5flliK/aMyk/rcgZne4KhEAvXl3W9Yzcs3FDGqVOLOhq67e2x//3XN8Y8CjwKICK/xurJaQmnTmgwjCG7TSjMTKEwM4XDR8Rupmh9i591e+q47anlXWYl2VrRyHefXcm4wnQeuGJ2VDeUigYvxduquP+dzWzYW4/DJhwzLp9vzZvAmIJ0fIFQu5+2jx2fz/dOn8Sb6/bx6/+up6LBR266k/NnDmPBjXP58zubeWfDgTcnbyDUZXaRumY/m8sa+ONbG/l4cyVgLYQ+cWIh3z1tImMK0slIcM8yxWnHYZOELjXI7QMTrcYPyqAww53QrDi3nDSegozof/fa5gArS2vi0CJYVVpDbbO/ow+au7AmrPSW/cOYIjLIGFMmIiOxnhfOBcZglW66h0NLON0mIs8Ac+ighJOI1IVLPC3GKuF0f7x/oUTSZ4Z9XH2Ln/+u3sslD30a1Y1oS3kj59z/Ect2VHf10B+wZt5d9/gSbnpyGWt31xEMGbyBEO9uLOeCBz7hD29sbDd4HTMun5tPHM9Vj37GA+9tYXdtC75giH11Xh76YCuXP7KYq+aO4pTJgw45trMMNa1LGL7yt0/2B0Kw6gG+t7Gcc/7yMS+v3E1DgjO0OGzChKLEPT/LS3d1mS4sIe1Ic3Hd8WMSdr0Up42TJhV263lhMBSKWz3EmiZ/Z8PW79N7vUNf+PqtnheRdcB/gFuNMTVoCadOaTDsw4LBEMu2V/P957uXF90XDHHNY0vYXtn5zM/aZj8/enE1K3d2vGTgycU7Dpk8YhP43pcmcctTSzt8PtPkC3Lr08v41injD7mp220d//PbWt7QZRLs/3lxDTsTPBEnO9XJCRMKE3a9uWPzo82AEhcOu42vHjGcYTmpCbnez845jJy07vf67XFaDmITgY5PXQz01rqhlvD1ATDGHG+MmWqMOdwYszC8TUs4dUKDYR9W3uDjW8+s6HrHTviCIW5+chkVnfQq65r9vLV+X5fn8h60cPnkyYP4sKSCupbOJyo0+YK8uW4fZxw2+IDt7eWDBKht9vG/b2/qsj0Af15YktDeocNu4/I5IxN2vZtPHNsngiFY6QMfuKLd9dgxdcSoHL502GDc3Ux5mOK0Myw3PkF7eG4qro6XenxG75VUcgJLeuna/YIGwz6qxR/kwfe3xKTE09aKRt7ZUNZh7cQ31u6NqPrCjqomRuR9cZM5bWoRL6/cHVEbXlqxm9PbBMPCDDdOR/sfsZt9QT4sqYjovG+u20djAlKjtZWd6uSoMXlxv87o/DSGZCemJxapsYXpPcq81JW8dBd/vnRWxEuV2pOZ4mTOmJ7lM+2IZ3RuhykUw6nQ/gUkei1MAJivqdh6RoNhHxXrtF9/fXczlQ3tP84o62Bm6sFeW7OHK+aM2v860+2MePF8VaPvgJ7gFXNHdjgxxB80EZdGCoZMwvNm5qS5+H/nT4v7zNLfX3Q4+d2YQBJPmSlOrjlmNF8/dnTMz52d6mTBjXMZGoMPAKdNjU+JvzOnDenqOea9JP6+agP+lOBr9jsaDPuoHVVNMU0Gvb2y4/MNj3BI6e31ZZw9Y8j+oFbb7It4+UZhpnt/LzfFaeMSz4gOM4s47ZFnPHHZbTh6YbnDsJxUvj0vfj2kS44cwaSizG5NIIm3nDQXt8+bwPfPmByzzDRjCtL5963HMqYgA1sMTlqY4WZcYWwnOo3MS2NoTpfZ1txAop+lGRJXP7Hf0mDYRy3+PPZlI7dWtD+KcurUoohuakOyU0hz2rnnq1bptP+u3ssFsyJLQvGVI4bx39XWbO1fnDeN/E6CaJrbzrwIK4t/uU1wTqQ0t4PL54xkztjYD5dOGJTBnadP6rW1lJHISXNx9dGjeO32ExhX2L2F8WBNwvrmKeN59sajGVOQHrPedkGmm5+dMzUm52r14y9PieTD3x0xvWjkeuu6/YYGwz4oEAyxcW90SbAj0VE+0Qy3g/MiCGp3nj6JnDQXx40v4OvHjubjLRV4Rud1uRYsO9XJSRMLWbihjEs8wzltSlGneVmzUpx857SJXRY+ttuEW08e32sTTPIz3PzlsiPwjMqN2TnHD8rgievmUBinLEKxlOZ2MGlwJv+68Wj+cvksJhVlRnys22Hj0iNHsPC7J3HTiePikjVpxvBsvtxxtpionDJ5EEeOzuu0p+7xeHKw8nUmOkuCHbggfH3VTX1jmpo6QLzGWDpaH5WV6uTHZ02hot7b4cSV2+dN4ORJhdhssn+YbEh2Cr96dR1/veIIbn1qGRXtPJPMSXPywBVHcM9rG/jWKeP52tGjyY1gcsTI3DT+eNHhfPfZle0+P7TbhL9cPiuSYau4Ksx089BVs7lvYQn/7GHB3y9PH8Ld506lMLN3f6doFWS4OXvGUI4em09lo4+PSir4ZEsFa3bVUd3kIxAypDhsjClMZ9aIXI6fUMDhw3PISnWQ6orfLSgnzcXd505lc1kDG/d1/8Pl2IJ07vnK9Ej+3R4NeOmdyhVerIX1rwOIyGPA2UCZMWZaeNvvgXOw1iRuAb4eXn+IiPwQqzhwEPiWMeaN8PYzgPuwAu4jxph7OIiIuLESd88GKoFLjDHb4vabxonEa6mIiKQAH2CNZTuA54wxPxORMcAzQD6wFLjKGNPpLAyPx2OKi4s726Xf+d3rG3jgvS0xPecDlx/BWTM6/qRc2eBla3kj979bwrrdddhtwgkTCrnhhLEMyko5ZBZdsy/I7ppm5n+2nTOnD2X5jmpeXL6LigYveekuzjt8GEeNyeO/q3dzkWckw3NTo+rFNXoD7K5p5q/vbubV1XvwBw0uu41zZg7llpPGWcO2cbyZRqO+xcqW8835yymNcu1jQYaLP158OIcPz0nasl4Ha/IFaPAGrE92AhjrA0xOmrPT9aXxUFbXwk1PLWXZ9uiz0kwfls3/Xe1hcFbX8c3j8fwU+Cm908nwA78oLi7+FYCInAA0YFWjaA2GpwPvGGMC4dylGGO+LyJTgflYJZ+GAm8DE8Pn3QSchpXdZglwmTFmXdsLi8gtwAxjzE0icilwgTHmkvj+urEXz2AoQLoxpkFEnMBHwO1YY9svGGOeEZEHgZXGmI7qZwEDMxh+WFLOVY9+FtNzvvPdExkbwaSC2mY/Xn8QBDJcDtK6CGB1zX4afQEq6r00+gKICMYY0l0O8jPcpLvtZKd2/ybf4A3Q5A0QCBkcNiHD3XWbekt5vZdtFdYHik82V3Y409UmcMTIXG6fN4FJQzLJT3f3St7TgaKq0ctrq/fyi1fWRVTf0O2w8YMzJ3Pu4UM7fb7dlsfjWQjM62FTe2JhcXHxqa0vRGQ08EprMGxLRC4ALjTGXBHuFWKM+U34vTeAu8O73m2M+Rr9LY4AACAASURBVFJ4+wH7tTnXG+H9FomIA6siRmGyLcqP2x0l/AfROmPDGf4ywCnA5eHtj2P9oXcaDAei8YMysNskZoVd89JdEU80yU51QhSTN7JSnWSlOuO2Ji7D7eiVSTLdUZjppjDTzf2DZ9HiD9HkC1Cyr4HKRh8YyE5zMrEog8wUJy6HrU/kHR0I8tLdXDh7OKdOLeLdDWU8uXg7G/fW4w9+8f/LYRMmFmVy+ZwRnDZ1MNmpDlKcUf276+2SRkOj2PdarDWRYLX70zbvtS3PdHBppzntnGt/Cahwr7MWa+QvssXCfURc7zAiYscaCh0P/BVrnLomXHASOqmJJSI3YFVdZuTIxGX86Csy3Q6+dFgR/13dfpWIaH392NE9WsisopOd6qL1s0Fv1QJUB3I77RQ57Vxy5AhOm1pEIGRo8QcJhgx2m+xPwp6b5uru8o7e/g8WURdWRH6EtVD/qfg2J7nENRgaY4LATBHJAV4EJkdx7MPAw2ANk8anhX1XRoqTu740mbfW7Tvg02t35KY5udgzAkc3KoYr1d+ISMRDn1Hq7RJOXWbPEJFrsCbWzGszjNlZCaeuSju1Pb40PEyajTWRJqkk5O4YnrH0LtZsq5zwHxj0w5pYsVSUncIPz5zS4/P86dJZ3a5tqJSKWG/fyzrNjRieGXoXcK4xpqnNWy8Dl4qIOzzBcQJWntUlwAQRGSMiLuDS8L4HexmrJBTAhViTdJKuAxO3YCgiheEeISKSijUjaT1WULwwvFvbmlrqIKlOOxccMYzLjxrR9c4d+MnZU5k1IkcnZygVf+9jzersDX7gvdYXIjIfWARMEpFSEbkO+AuQCbwlIivCExgxxqwFFgDrsJZm3GqMCYYfZ90GvIF1714Q3hcR+YWInBu+3KNAvohsxpog+YO4/7ZxEM/ZpDOwJsjYsYLuAmPML0RkLNbSijxgOXClMabT7v1AnE3aVnWjj+eWlfLb1zZEnIcz3WXnDxcdzrHjC/p0JhOl+guPx3MW1hKF7F64fC1waXFx8eu9cO1+IZ6zSVcBs9rZvhVrPYuKUG66iyvnjGLe5EHc89oG3l6/j45iotth47yZQ/n2qRMpzHThtCc6GYZSkWnxB6lt9lPX7GfN7lo2lzUQCBpy013MGpnD8JxUUpz2eD3fi4dF9F6OUDcHzghVUYpbzzCWBnrPsK36Fj8N3gDr99Tzyf9v787j46rLxY9/nlkyk0ySZm2a7gst3ShtCaVAKWUTlLJ4QWSTRQREQdArXLhXf/Lzyu+nl6ugKAgiiCwVFEXRK1spSwu0TRdoaelemi5JszV7Zsv3/nFOadpmmWT25Hm/Xnl15syZc55mJvPMdznfZ2sNFfVtOAQmFPs4ZUIRE4Zmk+uN78oeSkWjxR+iutnPE0t38KdVu3tckH76iFy+deZEZo/Oj8uSbTEmZWVli4Avkdgl2cLAC+Xl5Vf2uqfqln5ippkcr5scr3VN34JJxQTCHQjWtHGlUl1Ns58nl+3gkbe2ddu70dn6PY3c9PQqxhf5ePjq2YzKz0qZYsfd+ClwIZCVwHP67fOqKGjLMM21BcI0tgXZ19jGml0HqGn2k+N1M3NUHmMKs8j2uMjx6pihSr69B9r4ym9XsK26fzVoHQL3ffE4zj+uNFXHwQWgrKxsCXAqial6HwSWlZeXn5GAcw1oKf0VS3UvFO6gusnPQ29u4cXVe7pcYsohcPaUEu75whSGDfGSqa1HlSSVDW1c9uj7fV63tbMOA/f8eR3hDsPFM4eTnbpf8r4CfEJikqEfuDoB5xnw9CrsNNTcHmLZtlrOffAdnltR0e1aix0GXttQxeceeJsXVlZEXJVeqViqawlwaz8WMO/O915az/aarsuRpYLy8vLdwG1AvINsAb5VXl6e7OsbBwRNhmmmLRBi6dYarntyBY3tod6fAATDhh/87WOeWLqDhjZNiCpxAqEO3thYRfnO+pge9xvPrqa6qdcFV5Lpd8BviF9CbAF+U15e/mScjj/oaDJMM5WNfm5btLrLGn+9+cWbW/l4T2O3dQ2VirUDrQH+8+UNve/YR7vr23hl/b6UfS+Xl5cbrAvQ45EQW+zjanX7GNJkmEZqW/x8+/m1Ua1V+q0/rKG2JaW/UasBZOXOepr8kfVg9NWv395uVQNJUZ0S4m1YFXyiXZ0maB/ntvLy8m/bx1cxoskwjdQ2B1hb0fcCpZ3VNAdYvSu2XVbJYoyhpsnP9upmVn9az7rdB9jX0EZTe7JWxFKdHWgN8OSyHXE7/p4DbdSn+Dh4eXm5sbsypwDLgFas6wL7Imw/bxkwWbtG40Nnk/ZDuMNQ1+InYE9cyXA5KPR5+lv2JSLGGF5aE5tx8qff38UpE4rS+pKLmmY/Sz7Zz6+WbGVn7aE1h50O4ewpQ7n97EmMLsgk25O+/8d0Fwh1sH5vQ1zP8WHFASaV5MT1HLFgT6o5o6ysrAyrtfhFrJmgWXQ96zSIlQA9WBV/flau15fFlSbDPqpu8vPyh3t57J3tVDa2A1A6xMvN88dz/oxSinO8cTlvQ1uQlTvrYnKsdXsaaAuG0zYZVje1c9PTq1iz6+hWcrjD8OrHVby2oYr7Lp7OBccP7/b/2R60rtEMG6ueXX5WBm4tcxUz/lAH7cHeq8pH471tNVw8a0TavG52QruyrKwsD5gLlAELsArzerAS5F6sRbfLgQ/Ky8uj6w5SEdFk2Af7Gtq44rEPDmuJWNvbufflDTz9wac8+7W5DBsS+4QYDHewr6E9JsdqaAvSEd/PqLipbfZz66I1XSbCzoyBf//LeobmeDlz8tCjWu2VDe08+s42nl9ZQWsgTH6Wm2tPGctVJ42hOPWX/UoLwXD832T1LUH8oY60SYYH2QnuFfvnR0kOR6FjhhGrbfbzjWdWH5UIO9tW3cKtz62mLk6D+g6JXTdsDA+VMKGODnbVtbJ8e+Qt5Htf/vioCUN7D7Rx0a+W8uSynZ+ti1nfGuTBN7Zw1eMfpPqUfdWZ2Mu+KBUlTYYRqm0JsCaCySvln9ZTH4dk6HU7GVsYm+UOS3I9aZMMa5v9bKps4kd/38Br66v4xeItfXr+7vo2KuoOXezd4g9x/yubqGrsOuFtrmrmhfIKQglo1Qx0Ga74f7yUDsnE49aPMRU9fRdFwBjDouW7It7/+fKKmMeQ43VzyjFFMTnWzFF5ZKf2YscYY9hV18o1T6zg3Aff4fGlOyjMzuDtzdV9PtaLq3Z/drupPcTf1/VYEJwnlu6gLsVnKaaDDKeD3Mz4vs/mHVOIy6EfYyp6+i6KQDDcQVVT5ON1lQ3tcRkvWTijNCYtuq+dNp6sFC/xtK+hnUsefo+P9zYetr0/11gfaAt+NvM3EAr3ep1mbUugX+cJhDrY39jOngNtbNjbwNqKejZXNlHZ0EZ1kz9lLxCPl8wMJ7NH58f1HNNHJKOOrhqIUvsTMUW4HA4KfBkR71/gi8+sxByviy9ML+Uf6/b1+xgTirMZV+SLYVSx1+IP8dPXN1PdHJuxu2yPE7fT+haR4XLgdAjhHhJTbqaLvlwl0xoIUd8SYNGKXfxhZQU1zUe3KscV+bhh3ljOmTqMQl8GrjSb8NEfOV43Xz99Am9t6ntrPhLThuemfA+HSh8D/y8yBhwO4Yo5oyPe//I5o+ISx5DMDL6/cApD+lm+xukQHrl6NkUpXjm8qT3Ey2uP7spcv7eBk8cX9vl4F80agdhNap/HxdlThva4/1fmjiE/M7IvP7XNfp5ctpP597/FL5ds6zIRAuyoaeF7L33MOT97mzc/2T9oFk2fUOxjWG58Ljf61lkTKUzx97JKH5oMI1SS62VSSXav+00tzaU4jn+gRdkenv3aSWRl9K0ck0Pg4StnMzIvM06Rxc6+hjYCXXQzP/vBLq6aG/mXEoDiHA/HDD30uuV43fzH+VO7/UIxMj+Ta04eizuCyR/7m9q54amV3P/qph5bmp01toe46elV/Ncrm+Iy0SrVFPo8PPDl42N+3JmjhnDi2Ph2warBRZNhhIqyPfz22hMZ2sM1aMNyvTx2zQlx/bbqcjo4dlgOL982j4lDe0/OYCWE528+mXkTi8hKg26lUDdjettrWsh0O5lamhvxse75/GQKsw5v5Y3Iy+Tvt83j/ONKcdn9oV63gyvnjOJPXz+FkghaMrXNfm78fTlrK/q3wspzK3bxwBubaWgb2EvHORzC1OFDuGT2iJgd0+Ny8PPLZ1Hg01ahih2tdN9HVY3tPLF0B88t3/XZAsS5XhdXnTSa604dF9EHaazUNPt5d0s1v3xzK9uqj14Yf2iOhxtPG89Fs4ZTnO35rKswEcIdhvrWwKFJIwLZGa6IkvGntS2cfv9bXT5Wkuvh11efwHf/+FGvFdO/ffZErj1lLHlZXXd5NrcHafaHCYY7yHA5yPW6yezU4m5uD9Ie7EAE8jLdOO1xvlZ/iN+8u50H3ujbZR5dee5rJzF3fGFcl/JLBbXNfm55ZhUroizl5HYKz9xwEjNH5eFJvWLVA/tFHOA0GfZDINRBfWuAtmAYwboGMD8rIyHXVXWlpsmPP9RBdXM7DW1BstwuSvO8eFxOCrIOfYgnQos/RFN7kH+s28dzy3exv9FPsKODXK+b+ZOKuWn+eAp9GT22nutaAlz/5Ao+3N11q2tkfiYPfnkmb36yn+dXVhxVuWD26HzuPPdYpg3PJbcf46t1LQF21bbw8Fvb2Fbdgtft4OKZI7jg+OEU+NxUNvpZcP+Sfs04PVJxtof/uX1e3JbxSyU1zX6+/9J6/rm+sl/PH5Lp5rfXlTGtNJfM1JwNrckwjWkyVDFT1djO//ufjfzjo32EesgUU0tzue+L05k8LPewllhnK3fWcdmj73dbt9HtFM6dNoy7zptMXUuAqsZ2PC4H44p85Hrd5Pdh9m9nNc1+7vzjhyzpYgakx+Xg5VtPZdHKCp5ctrNfx+/KE9eVcebkkpgdL5UdaA1Q/mk933lhLY1tkZd2+vz0YfzggmkUZ2ck9MtdH2kyTGMp+fVKpZ/d9a1c/tgH7K5v63XfDfsaufTX7/PQFbM4fVIxvi66TqeU5vDAZTP5zgtru2yBGWNdd1noy2B0QWxW5mlqD/L///lJl4kQrIWnKxv9/Kl8d5eP99fj7+7ghDEF/Z4lnE7ysjI4Y1Ixb3zndJZvr+ORt7axsbKxyy89OR4XF80czlfnjaPAl9Ftd7dSsaDJUEWtqrGdK3+zPKJEeFC4w3Drc6t55oaTOHFcPm7n4S3EbI+bc6aW8NadZ/DUezt5cfVuGtqC5GdlcFnZSK4+aQyFORlkumP3Fm5qD/GX1T0nOq/bGfNitet2N+APhem6kk9i1bUEONAaYHOVNR47pTQnqpZ2V5xOB0NzvFxw/HBOmVBIMGxobAtS2dhOuMOQ7XUxLNeLx+UgN9ONN/XGBtUApMlQRaU9GObhJVvZVdf9Aubd6TBwy7Oref3b8xmae/QHns/jwudxcde5x3Lz/PFWJ5SxWhfxGJ99e3N1j+OARdkZcVnEu8kfIhhK7nCFPxhmR00Ld734ER8dMVZbNiafn1w6g9EFWTFfTOLg2PGwIV4mDUv9uoRq4ErZzneVHhragrwQRbdhQ1uQdXt6vjzB43YyNNfL0BwvQ3O9cZuoVNNLovO4nLSH+lqkPDIdSRy7N8awZX8zF/1q2VGJEKzF5y98aCmf1h49Y1mpgSJuyVBERonIEhHZICIfi8jt9vYCEXldRLbY/+qVs2ls2dYa2oLRJYifL95CbYyWXovGiPyeFyRoC4bxxWkWozOJl1bUNge49bnV+EPdr6fbEgjznRc+jFt5MqWSLZ4twxDwr8aYqVgVnb8pIlOBu4HFxpiJwGL7vkpDB1qt9Tij9dHuhqgTaiycMqGIjB66AetaAhRmx34SR36WG5czeclwX0Nbj3U6D/podwP1g2QZOTX4xC0ZGmP2GWNW27ebgI3ACOAi4Cl7t6eAi+MVg4qvYNjErKXQWyWJRMjxurh+3tge92luD/Vp0fZIzBqdT1YSJ4ks21Yb8b4fRlDTU6l0lJAxQxEZC8wClgMlxpiDZRcqgS4vsBKRm0SkXETKq6vjs+q9ipaJeE3OSI6VbD6Pi5vnT+DKbhZlL/BlML7Yx9V9XB+1NzfPH0+2N3kzSfvyGvZ0/ahS6Szus0lFJBt4EbjDGNPYeUkwY4wRkS7/uowxjwGPgXXRfbzjVH3ncjisa78i6GKL5FipoMCXwb+dN5mbTx/PM+9/yobKRjLdTi4rG8XM0XkU+jxcfuJoHn17e49jbJEaVZDJ+OLI1piNl77UHJw2PPJ1YZVKJ3FNhiLixkqEzxpj/mxvrhKRUmPMPhEpBfbHMwYVP7mZbhbOKGVtlF1n44p8KXUt2ZAsN0Oy3Nz1+cm0B8M4HXJYMeRCXwbfXziV7720PqrziMAjV51AcQ+LvyfChKE+irM9vdaPHFfko2QQLBunBqd4ziYV4LfARmPMzzo99DfgWvv2tcBf4xWDii+nQ1g4Y/hnlR/66xsLJlAUh4kp0XI7HeR43YclQrAu9fjCcaWcObk4quPfee6xjC2Mzeo50SjIyuAnl87ocR+nQ/ivS2fEZQKRUqkgnn1TpwJfAc4UkbX2zxeAHwPniMgW4Gz7vkpT2R4n500f1u/ne90O5k8qTmhFjVgo8GVw/6XHc+60/v3f7zz3WK6cMzqpY4UHuZwO5ozN5/FrysjPOjqe4mwPz9xwEtNKc9PudVIqUrpQt4paRV0r5z/0bp8WXj7ox5ccx0Uzh8d0WbVEqm8J8O6Wau758zpaAr1fHjJ8iJdHrj6B8cU+clIgEXYWDIepbw2ycW8j72+vRRBOm1TEMcXZFPpSeoHsVKHfFNKYJkMVtXCHYVNlI19+9IM+rdt5+1kT+eq8sQzJTO+ut2A4TH1LkLc2V/PUezvZVNl02KzLTLeTmaPyuGXBBKaU5iZ9jFDFjSbDNKbJUMVEKNxBRX0btzyzik8qm3rcNzfTxb0XTOOsKSUDqlKDMdZ1l6EOgz8YpsNYY20elwOP25H2SV/1SpNhGtNkqGKqttlPTXOAX7+9jX98tI9A+NDlB9NH5HL7WRM5fmRe3BbbViqJNBmmMU2GKi5a/CFa/CGCHQZjDC6HkOF0UNBDhXul0pwmwzSWnrMWVMo7WH5JKaXSgfZTKaWUGvT0q7tSXQiFOz6bDOMQq8J9XpZOgFFqoNJkqNQR9je28/zKCp58b+dnVTlOHl/Idz43iUklOQNqBqxSyqITaJTqZHd9K5c+8j6Vje1dPn7H2RO5/tRxmhBVV3QCTRrTMUOlbLXNfm5+elW3iRDgwTe2sKmyMYFRKaUSQZOhUrbqZj8f7+090f3s9c0c0IrvSg0oOmaolO2f6yoj2u+D7XW0B6OvZRiNA60BWgNh2oPWeqhet5OsDJ3ko1R/aTJUytbcHvm6qh1JGmuvbfazpaqZBxdv5oPtdYc9NmdcAXecNZFJw3Io0sUNlOoTTYZK2aaPjKyKe36WO+oajv1R2dDOjb9fybo9XXflrthRx5WPL2dqaS5PXFfGsCGZCY5QqfSlY4ZK2U6ZUIQngvVSrzl5LPm+xHZH7m9q56rHP+g2EXa2YV8jV/xmOft7mAiklDqcJkOlbLleN3edd2yP+wzN8XDlSaNxJ7C2X6s/xE9f28y26paIn7OjpoWfvLKJlj6U1FJqMNNkqJQtM8PJJbNH8r3zp3TZQpw8LIcXbzmFklxvQuNq8od4ac2ePj/v5Q/30tSHcVClBjMdM1Sqk7ysDK46aQwXHD+cpVtqWLenAZ/HycIZwynO8SR8YooxhsUbq/CH+j57NRDu4NWPK/nK3DE4kjDGqVQ60WSo1BEyM5xWK/GEkVxywsikxtIaCPP6hv39fv7ijVVcMnsE2V5dMUepnmg3qVIpLBQ2tAT639XZ7A8T7Ej9JReVSjZNhkqlMJdT8GX0vwMn2+PErV2kSvVKk6FSKSzT7eT0ScX9fv5pE4vJjCKZKjVYaDJUKoU5HMJ500twO/veunM5hIXHl+LUlqFSvdJkqFSKy/a6WThjeJ+fd970YeR4tFWoVCQ0GSqV4rI9Lu7+/GRGFUS+vNrI/Ey+t3CqziJVKkKaDJVKAyW5Xp6/6WQmDs3udd8Jxdm8cPPJDEvw4gBKpTPtQ1EqTQzPy2TRjXNZW3GAny/ewro9DYc9Pm14LrefNZFZo/MpztGqFUr1RdySoYg8ASwE9htjptvbCoDngbHATuAyY0x9vGJQaqApyvFw9tQSZo/JozUQptkfAgM+j4usDCeFWrpJqX6JZzfp74Dzjth2N7DYGDMRWGzfV0r1UYHPw8j8LCYPy2VyaS6jCrI0ESoVhbglQ2PMO0DdEZsvAp6ybz8FXByv8yullFKRSvSYYYkxZp99uxIo6W5HEbkJuAlg9OjRCQhNqfgLhMI0tYfoMAYQPG4HuTrjU6mkS9oEGmOMEZFuF000xjwGPAZQVlamiyuqtFXX4qehLcTKnXUs21rD5qom2gJhnA4Hw4Z4mDOukFPGFzKmKIv8zAzcERQYVkrFVqKTYZWIlBpj9olIKdD/5fiVSnF1zQHWVNTz4BtHz/w8aFt1M8u21vIAkJvp4pq5Y7hq7hiKsjNwO52JDVipQSzRX0H/Blxr374W+GuCz69U3PlDYXbVtnLD71dyw1Pl3SbCIzW2hfjlkm2c/dO3eWV9FfUtgThHqpQ6KG7JUEQWAe8Dx4rIbhG5AfgxcI6IbAHOtu8rNWC0+EO8v62Wcx54mzW7DvTvGIEwty1aw7+9+BG1zf4YR6iU6krcukmNMVd089BZ8TqnUsnUGrAS4Y1Pl2NiMMr92oYqmv0hHrpill42oVSc6Ui9UjFSUdfGzc+sikkiPOi9bbX892ubaGwLxu6gSqmjaDJUKgZqmv3c8swqwnGoKr9oRQWfVDZhYplllVKH0WSoVJRC4Q5eXLWb7TUtcTvHrc+tpkbHD5WKG02GSkWprjXAo+9sj+s59jf5+XhvY1zPodRgpslQqSht3NtIXQIug3jwjS3UtWjrUKl40GSoVBT8wTAvrt6dkHOtrTiAP9SRkHMpNdhoMlQqCo3tIT7cHdlF9bGQiBaoUoORJkOlohDuMHxa25qw832UwMSr1GCiyVCpKIQ7EtttqS1DpeJDk6FSSqlBT5OhUlFwOhyIJO58xTm6LJtS8aDJUKkouBwwvsiXsPPNGDEkYedSajDRZKhUFHIzM5gxMi9h58vLcifsXEoNJpoMlYpChsvBZWWjEnKuOeMK8Li14K9S8aDJUKkoTSj2JWQs746zJpKflRH38yg1GGkyVCpKhb4Mbj3jmLieY2R+JhNLsuN6DqUGM02GSkXJ6XSwcEYpU0pz4nJ8EXj4qtkU53jjcnyllCZDpWKiMNvDr66cjccV+z+pr80bx7gEzlhVajDSZKhUjIzIz+R315+IyxG7Cw8/N62EWxYcQ45XZ5EqFU+aDJWKEY/LyczReSy6cS7ZHlfUx/vSCSP58b8cR4FPJ80oFW+aDJWKoUy3i5mj83jt2/NZMKm4X8fIz3Lz+6/O4T/On0KBT1ecUSoRov/6qpQ6jNvpYHheJg98eSZbq5v5xeItLN1agzE9P29Yrpeb5o/n/BmlFGV7cMawu1Up1TMxvf2FpoCysjJTXl6e7DCU6peG1iDNgRCb9jWydGsNG/Y10hYI43Q4KB3i4dRjipk1Oo/C7AwKfZoE05i+cGksLZKhiFQDn3bxUBFQk+BwepNqMaVaPJB6MaVaPKAxRSLV4qkxxpyX7CBU/6RFMuyOiJQbY8qSHUdnqRZTqsUDqRdTqsUDGlMkUi0eld50Ao1SSqlBT5OhUkqpQS/dk+FjyQ6gC6kWU6rFA6kXU6rFAxpTJFItHpXG0nrMUCmllIqFdG8ZKqWUUlHTZKiUUmrQS9tkKCLnicgmEdkqIncn4fxPiMh+EVnfaVuBiLwuIlvsf/MTHNMoEVkiIhtE5GMRuT2ZcYmIV0RWiMiHdjz/194+TkSW26/d8yKS8MU3RcQpImtE5O+pEJOI7BSRdSKyVkTK7W1Jez+JSJ6I/ElEPhGRjSJycpLjOdb+3Rz8aRSRO5L9N6cGjrRMhiLiBH4FfB6YClwhIlMTHMbvgCMvsL0bWGyMmQgstu8nUgj4V2PMVGAu8E3795KsuPzAmcaY44GZwHkiMhf4CfCAMeYYoB64IUHxdHY7sLHT/VSI6QxjzMxO184l8/30c+AVY8xk4His31XS4jHGbLJ/NzOBE4BW4C/JjEkNMMaYtPsBTgZe7XT/HuCeJMQxFljf6f4moNS+XQpsSvLv6a/AOakQF5AFrAZOwlo1xNXVa5mgWEZifXCeCfwdaxmtZMe0Eyg6YltSXjdgCLADe4JdsuPpIr7PActSKSb9Sf+ftGwZAiOAik73d9vbkq3EGLPPvl0JlCQrEBEZC8wClpPEuOzuyLXAfuB1YBtwwBgTsndJxmv3IHAX0GHfL0yBmAzwmoisEpGb7G3Jet3GAdXAk3ZX8uMi4ktiPEe6HFhk306VmFSaS9dkmPKMMQbrAy7hRCQbeBG4wxjTmMy4jDFhY3VtjQTmAJMTde6uiMhCYL8xZlUy4+jCPGPMbKyu/2+KyPzODyb4dXMBs4FHjDGzgBaO6H5M1vvbHsu9EPjjkY8l829Opb90TYZ7gFGd7o+0tyVblYiUAtj/7k90ACLixkqEzxpj/pwqcRljDgBLsLog80TkYPmwRL92pwIXishO4A9YXaU/T3JMGGP22P/uxxoLm0PyXrfdwG5jzHL7/p+wkmPS30dYXxZWG2Oq7PupEJMaANI1Ga4EJtozADOwuk3+luSYwIrhWvv2tVhjdgkjIgL8FthojPlZsuMSkWIRybNvZ2KNX27ESoqXJjoeAGPMPcaYkcaYsVjvmzeNMVclMyYR8YlI9OUHFwAABAlJREFUzsHbWGNi60nS62aMqQQqRORYe9NZwIZkxXOEKzjURQqpEZMaANJ2BRoR+QLW2I8TeMIYc1+Cz78IWIBVRqYK+AHwEvACMBqr5NRlxpi6BMY0D3gXWMeh8bB/xxo3THhcIjIDeArrNXIALxhjfigi47FaZQXAGuBqY4w/3vF0Ed8C4LvGmIXJjMk+91/suy7gOWPMfSJSSJLeTyIyE3gcyAC2A9djv4bJiMeOyQfsAsYbYxrsbUn7HamBJW2ToVJKKRUr6dpNqpRSSsWMJkOllFKDniZDpZRSg54mQ6WUUoOeJkOllFKDniZDpeLErkRRlOw4lFK902SolFJq0NNkqBJORO4UkW/Ztx8QkTft22eKyLMicoVd22+9iPyk0/OaReR+uzbiGyIyR0TeEpHtInKhvY/T3meliHwkIjfb2xfY+x6s0fesvWLPkbFdYNc1XGOfo8Tefq+IPC0i79u1827sdNx3ROQfYtXX/LWIOI445g9F5I5O9+8Tu9akUio1aDJUyfAucJp9uwzIttdUPQ3YjFVb8EysGognisjF9r4+rOXTpgFNwI+wlnj7IvBDe58bgAZjzInAicCNIjLOfmwWcAdWDczxWOuUHmkpMNdeoPoPWNUtDpphx3Uy8H9EZLi9fQ5wm33cCcC/HHHMJ4BrAOxEeTnwTM+/IqVUIrl630WpmFsFnCAiuVgFgFdjJcXTgJeBt4wx1QAi8iwwH2upuwDwin2MdYDfGBMUkXVYtSXBWtdzhogcXGd0CDDRfu4KY8xu+7hr7ecsPSK2kcDz9qLPGVh1/Q76qzGmDWgTkSVYSfCAfdzt9nEXAfOwFrcGwBizU0RqRWQWVomhNcaY2j7/1pRScaMtQ5VwxpggVpK5DngPq6V4BnAMVpHb7gTNofUDO7ASKcaYDg59sRPgNmNXRTfGjDPGvGY/1nmt0TBdfxl8CPilMeY44GbA2zn0I/8rvWzv7HGs/+/1WC1FpVQK0WSokuVd4LvAO/btr2MtkL0COF1EikTEiVWl4O0+HPdV4Ba72xURmWQv8BypIRwq33TtEY9dJCJee3HoBVjVUwDm2BVUHMCXObq1CdZC3Odhdd2+2od4lFIJoMlQJcu7QCnwvl2brh14165afjdWSaUPgVXGmL6U5Xkcq9zQahFZDzxKL8MB9gSXC+279wJ/FJFVQM0Ru35kx/UB8J/GmL329pXAL7HKU+3gUAWKzxhjAvZzXzDGhPvw/1FKJYBWrVAqAiJyL9BsjPnvI7YvwC4D1cvzHVhjo18yxmyJV5xKqf7RlqFScSYiU4GtwGJNhEqlJm0ZKqWUGvS0ZaiUUmrQ02SolFJq0NNkqJRSatDTZKiUUmrQ02SolFJq0Ptfen1Ffv4/5QcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import pyreadr\n", "import numpy as np\n", "result = pyreadr.read_r('berkeley.RData') # also works for Rds\n", "\n", "simpsons_berkeley=result['berkeley']\n", "simpsons_berkeley_small=result['berkeley.small']\n", "\n", "import seaborn as sns\n", "\n", "simpsons_berkeley['number.apply.sqrt']=np.sqrt(simpsons_berkeley['number.apply'])\n", "\n", "g=sns.relplot(x='women.apply', y='total.admit', size='number.apply', sizes=(100,1000), data=result['berkeley'])\n", "figure = g.fig\n", "glue(\"berkeley_fig\",figure, display=False)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```{glue:figure} berkeley_fig\n", ":figwidth: 600px\n", ":name: \"fig:berkeley\"\n", "\n", "The Berkeley 1973 college admissions data. This figure plots the admission rate for the 85 departments that had at least one female applicant, as a function of the percentage of applicants that were female. The plot is a redrawing of Figure 1 from Bickel, Hammel, & O’Connell (1975) {cite}`Bickel1975`. Circles plot departments with more than 40 applicants; the area of the circle is proportional to the total number of applicants. \n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "Before leaving this topic entirely, I want to point out something else really critical that is often overlooked in a research methods class. Statistics only solves *part* of the problem. Remember that we started all this with the concern that Berkeley's admissions processes might be unfairly biased against female applicants. When we looked at the \"aggregated\" data, it did seem like the university was discriminating against women, but when we \"disaggregate\" and looked at the individual behaviour of all the departments, it turned out that the actual departments were, if anything, slightly biased in favour of women. The gender bias in total admissions was caused by the fact that women tended to self-select for harder departments. From a legal perspective, that would probably put the university in the clear. Postgraduate admissions are determined at the level of the individual department (and there are good reasons to do that), and at the level of individual departments, the decisions are more or less unbiased (the weak bias in favour of females at that level is small, and not consistent across departments). Since the university can't dictate which departments people choose to apply to, and the decision making takes place at the level of the department it can hardly be held accountable for any biases that those choices produce. \n", "\n", "That was the basis for my somewhat glib remarks earlier, but that's not exactly the whole story, is it? After all, if we're interested in this from a more sociological and psychological perspective, we might want to ask *why* there are such strong gender differences in applications. Why do males tend to apply to engineering more often than females, and why is this reversed for the English department? And why is it it the case that the departments that tend to have a female-application bias tend to have lower overall admission rates than those departments that have a male-application bias? Might this not still reflect a gender bias, even though every single department is itself unbiased? It might. Suppose, hypothetically, that males preferred to apply to \"hard sciences\" and females prefer \"humanities\". And suppose further that the reason for why the humanities departments have low admission rates is because the government doesn't want to fund the humanities (Ph.D. places, for instance, are often tied to government funded research projects). Does that constitute a gender bias? Or just an unenlightened view of the value of the humanities? What if someone at a high level in the government cut the humanities funds because they felt that the humanities are \"useless chick stuff\". That seems pretty *blatantly* gender biased. None of this falls within the purview of statistics, but it matters to the research project. If you're interested in the overall structural effects of subtle gender biases, then you probably want to look at *both* the aggregated and disaggregated data. If you're interested in the decision making process at Berkeley itself then you're probably only interested in the disaggregated data. \n", "\n", "In short there are a lot of critical questions that you can't answer with statistics, but the answers to those questions will have a huge impact on how you analyse and interpret data. And this is the reason why you should always think of statistics as a *tool* to help you learn about your data, no more and no less. It's a powerful tool to that end, but there's no substitute for careful thought." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Statistics in psychology\n", "\n", "I hope that the discussion above helped explain why science in general is so focused on statistics. But I'm guessing that you have a lot more questions about what role statistics plays in psychology, and specifically why psychology classes always devote so many lectures to stats. So here's my attempt to answer a few of them...\n", "\n", "\n", "**Why does psychology have so much statistics?**\n", "\n", " To be perfectly honest, there's a few different reasons, some of which are better than others. The most important reason is that psychology is a statistical science. What I mean by that is that the \"things\" that we study are *people*. Real, complicated, gloriously messy, infuriatingly perverse people. The \"things\" of physics include object like electrons, and while there are all sorts of complexities that arise in physics, electrons don't have minds of their own. They don't have opinions, they don't differ from each other in weird and arbitrary ways, they don't get bored in the middle of an experiment, and they don't get angry at the experimenter and then deliberately try to sabotage the data set (not that I've ever done that...). At a fundamental level psychology is harder than physics.\n", "\n", " Basically, we teach statistics to you as psychologists because you need to be better at stats than physicists. There's actually a saying used sometimes in physics, to the effect that \"if your experiment needs statistics, you should have done a better experiment\". They have the luxury of being able to say that because their objects of study are pathetically simple in comparison to the vast mess that confronts social scientists. It's not just psychology, really: most social sciences are desperately reliant on statistics. Not because we're bad experimenters, but because we've picked a harder problem to solve. We teach you stats because you really, really need it.\n", "\n", "**Can't someone else do the statistics?**\n", "\n", "To some extent, but not completely. It's true that you don't need to become a fully trained statistician just to do psychology, but you do need to reach a certain level of statistical competence. In my view, there's three reasons that every psychological researcher ought to be able to do basic statistics:\n", "\n", "\n", " - Firstly, there's the fundamental reason: statistics is deeply intertwined with research design. If you want to be good at designing psychological studies, you need to at least understand the basics of stats. \n", " - Secondly, if you want to be good at the psychological side of the research, then you need to be able to understand the psychological literature, right? But almost every paper in the psychological literature reports the results of statistical analyses. So if you really want to understand the psychology, you need to be able to understand what other people did with their data. And that means understanding a certain amount of statistics.\n", " - Thirdly, there's a big practical problem with being dependent on other people to do all your statistics: statistical analysis is *expensive*. If you ever get bored and want to look up how much the Australian government charges for university fees (*The original author of this chapter is Australian*), you'll notice something interesting: statistics is designated as a \"national priority\" category, and so the fees are much, much lower than for any other area of study. This is because there's a massive shortage of statisticians out there. So, from your perspective as a psychological researcher, the laws of supply and demand aren't exactly on your side here! As a result, in almost any real life situation where you want to do psychological research, the cruel facts will be that you don't have enough money to afford a statistician. So the economics of the situation mean that you have to be pretty self-sufficient. \n", "\n", "Note that a lot of these reasons generalise beyond researchers. If you want to be a practicing psychologist and stay on top of the field, it helps to be able to read the scientific literature, which relies pretty heavily on statistics. \n", "\n", "**I don't care about jobs, research, or clinical work. Do I need statistics?**\n", "\n", "Okay, now you're just messing with me. Still, I think it should matter to you too. Statistics should matter to you in the same way that statistics should matter to *everyone*: we live in the 21st century, and data are *everywhere*. Frankly, given the world in which we live these days, a basic knowledge of statistics is pretty damn close to a survival tool! Which is the topic of the next section..." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Statistics in everyday life\n", "\n", ">*\"We are drowning in information, but we are starved for knowledge\"*\n", ">\n", ">-Various authors, original probably John Naisbitt\n", "\n", "\n", "When I started writing up my lecture notes I took the 20 most recent news articles posted to the ABC news website. Of those 20 articles, it turned out that 8 of them involved a discussion of something that I would call a statistical topic; 6 of those made a mistake. The most common error, if you're curious, was failing to report baseline data (e.g., the article mentions that 5% of people in situation X have some characteristic Y, but doesn't say how common the characteristic is for everyone else!) The point I'm trying to make here isn't that journalists are bad at statistics (though they almost always are), it's that a basic knowledge of statistics is very helpful for trying to figure out when someone else is either making a mistake or even lying to you. In fact, one of the biggest things that a knowledge of statistics does to you is cause you to get angry at the newspaper or the internet on a far more frequent basis!\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "## There's more to research methods than statistics\n", "\n", "So far, most of what I've talked about is statistics, and so you'd be forgiven for thinking that statistics is all I care about in life. To be fair, you wouldn't be far wrong, but research methodology is a broader concept than statistics. So most research methods courses will cover a lot of topics that relate much more to the pragmatics of research design, and in particular the issues that you encounter when trying to do research with humans. However, about 99% of student *fears* relate to the statistics part of the course, so I've focused on the stats in this discussion, and hopefully I've convinced you that statistics matters, and more importantly, that it's not to be feared. That being said, it's pretty typical for introductory research methods classes to be very stats-heavy. This is not (usually) because the lecturers are evil people. Quite the contrary, in fact. Introductory classes focus a lot on the statistics because you almost always find yourself needing statistics before you need the other research methods training. Why? Because almost all of your assignments in other classes will rely on statistical training, to a much greater extent than they rely on other methodological tools. It's not common for undergraduate assignments to require you to design your own study from the ground up (in which case you would need to know a lot about research design), but it *is* common for assignments to ask you to analyse and interpret data that were collected in a study that someone else designed (in which case you need statistics). In that sense, from the perspective of allowing you to do well in all your other classes, the statistics is more urgent. \n", "\n", "But note that \"urgent\" is different from \"important\" -- they both matter. I really do want to stress that research design is just as important as data analysis, and this class does spend a fair amount of time on it. However, while statistics has a kind of universality, and provides a set of core tools that are useful for most types of psychological research, the research methods side isn't quite so universal. There are some general principles that everyone should think about, but a lot of research design is very idiosyncratic, and is specific to the area of research that you want to engage in.\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## References\n", "\n", "```{bibliography} ../../references.bib\n", ":filter: docname in docnames\n", "```" ] } ], "metadata": { "celltoolbar": "Tags", "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.9" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 2 }