06-data-wrangling.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "## Assignment 6: Data wrangling I\n",
    "\n",
    "In the last two assignments you've explored data with both basic statistical methods and different types of data visualization.  While this is a key part of the data analysis process, datasets need to be in the right format before analysts can start drawing meaningful conclusions.  The process of preparing data for analysis is called __data wrangling__, and often takes the bulk of an analysts' time during a data project.  Possible issues might include: \n",
    "\n",
    "* Missing data or problematic/incorrect values in a dataset; \n",
    "* Data are formatted incorrectly, preventing the analyst from working with the data in the right way; \n",
    "* Data are spread across multiple files or data tables; \n",
    "* Data are in the wrong \"shape\" for analysis and visualization\n",
    "* All of the above, in varying capacities!\n",
    "\n",
    "A major reason why you are learning to work with `pandas` in this course is because it can flexibly handle all of these tasks.  In this notebook, we'll be going over some basic examples of how this works, which you'll then put into practice with the Exercises at the end.  \n",
    "\n",
    "In this notebook, we'll be working primarily with the US Department of Education's College Scorecard dataset. Read in the full dataset with the code below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(7804, 122)"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "\n",
    "full_url = 'http://personal.tcu.edu/kylewalker/data/colleges.csv'\n",
    "full = pd.read_csv(full_url, encoding = 'latin_1')\n",
    "full.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "We see that our data have over 7800 rows, and a whopping 122 columns!  At this point, while you might not know yet what insights are contained in the data, you will have some sense of the research questions that you are interested in.  As your research questions don't likely require all 122 columns, you can safely restrict your dataset to only those columns that you need for your analysis.  Our core research question is as follows: \n",
    "\n",
    "* How does the proportion of non-traditional students - defined as undergraduate students aged 25 and above - vary among comparable bachelor's-granting colleges and universities?  \n",
    "\n",
    "To accomplish this, we'll need to identify some columns that we need and subset our data accordingly.  The columns we'll be keeping are as follows: \n",
    "\n",
    "* `INSTNM`: The name of the institution; \n",
    "* `STABBR`: The state the institution is located in; \n",
    "* `PREDDEG`: The primary type of degree granted by the institution; codes include 1 for certificates, 2 for associate's degrees, 3 for bachelor's degrees, and 4 for graduate degrees; \n",
    "* `CONTROL`: The ownership of the institution, coded as 1 for public non-profit, 2 for private non-profit, and 3 for private for-profit; \n",
    "* `UGDS`: The number of undergraduates enrolled at the institution; \n",
    "* `UG25abv`: The percentage of undergraduates at the institution aged 25 and above.  \n",
    "\n",
    "## Subsetting data\n",
    "\n",
    "As we discussed in class, `pandas` includes many different methods for subsetting data; I encourage you to review the corresponding lecture notes for the full set of methods that we discussed.  In this notebook, we'll be focusing on the methods that allow us to accomplish the task at hand, and you'll be learning a few new methods as well.  \n",
    "\n",
    "To subset our data column-wise, we can specify a list of column names that we want to keep then use the `.filter()` method to restrict our data to only those specified columns.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>INSTNM</th>\n",
       "      <th>STABBR</th>\n",
       "      <th>PREDDEG</th>\n",
       "      <th>CONTROL</th>\n",
       "      <th>UGDS</th>\n",
       "      <th>UG25abv</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Alabama A &amp; M University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4051.0</td>\n",
       "      <td>0.1049</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>University of Alabama at Birmingham</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>11200.0</td>\n",
       "      <td>0.2422</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Amridge University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>322.0</td>\n",
       "      <td>0.8540</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>University of Alabama in Huntsville</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5525.0</td>\n",
       "      <td>0.2640</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Alabama State University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5354.0</td>\n",
       "      <td>0.1270</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                INSTNM STABBR  PREDDEG  CONTROL     UGDS  \\\n",
       "0             Alabama A & M University     AL        3        1   4051.0   \n",
       "1  University of Alabama at Birmingham     AL        3        1  11200.0   \n",
       "2                   Amridge University     AL        3        2    322.0   \n",
       "3  University of Alabama in Huntsville     AL        3        1   5525.0   \n",
       "4             Alabama State University     AL        3        1   5354.0   \n",
       "\n",
       "   UG25abv  \n",
       "0   0.1049  \n",
       "1   0.2422  \n",
       "2   0.8540  \n",
       "3   0.2640  \n",
       "4   0.1270  "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "keep = ['INSTNM', 'STABBR', 'PREDDEG', 'CONTROL', 'UGDS', 'UG25abv']\n",
    "\n",
    "df = full.filter(keep)\n",
    "\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "Our data are much simpler to work with now!  We can already see some noticeable variations in the data; over 85 percent of students at Amridge University are over age 25; however, the university itself is quite small with only 322 undergraduates.  These are the types of things we'll want to account for in our analysis.  \n",
    "\n",
    "You may have noticed that I decided to assign the result of the subsetting operation to a new data frame.  Certainly, I could have assigned the result back to the original frame: \n",
    "\n",
    "```python\n",
    "full = full.filter(keep)\n",
    "```\n",
    "\n",
    "Or: \n",
    "\n",
    "```python\n",
    "full.filter(keep, inplace = True)\n",
    "```\n",
    "\n",
    "My personal preference when doing data programming is to assign the results of major operations to new data frames, creating a data frame object that represents each step of the analysis, and do minor operations in place.  For example, had I assigned the results of the subsetting operation back to `full`, and I later decided that I needed an additional column from the original data frame, I'd need to read the whole thing in again rather than just add a column name to the list `keep`.  Python will hold all of our data _in memory_, so it will be accessible to us throughout our Python session; given that our data are relatively small, this won't cause us any problems.  Bigger data workflows might require different methods, however.  \n",
    "\n",
    "There are small things, however, that you can do to your data frame in place to make your lives easier.  For example, I don't really want to hit Caps Lock or the Shift key every time I am typing out column names, but my column names are capitalized.  Our column names in our data frame are simply a list of strings, which we've learned how to work with already with string methods. As such, we can use __list comprehension__ to convert all of the column names to lower case: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>instnm</th>\n",
       "      <th>stabbr</th>\n",
       "      <th>preddeg</th>\n",
       "      <th>control</th>\n",
       "      <th>ugds</th>\n",
       "      <th>ug25abv</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Alabama A &amp; M University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4051.0</td>\n",
       "      <td>0.1049</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>University of Alabama at Birmingham</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>11200.0</td>\n",
       "      <td>0.2422</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Amridge University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>322.0</td>\n",
       "      <td>0.8540</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>University of Alabama in Huntsville</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5525.0</td>\n",
       "      <td>0.2640</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Alabama State University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5354.0</td>\n",
       "      <td>0.1270</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                instnm stabbr  preddeg  control     ugds  \\\n",
       "0             Alabama A & M University     AL        3        1   4051.0   \n",
       "1  University of Alabama at Birmingham     AL        3        1  11200.0   \n",
       "2                   Amridge University     AL        3        2    322.0   \n",
       "3  University of Alabama in Huntsville     AL        3        1   5525.0   \n",
       "4             Alabama State University     AL        3        1   5354.0   \n",
       "\n",
       "   ug25abv  \n",
       "0   0.1049  \n",
       "1   0.2422  \n",
       "2   0.8540  \n",
       "3   0.2640  \n",
       "4   0.1270  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.columns = [x.lower() for x in df.columns]\n",
    "\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "You've also already learned to modify the column names of your data frame by passing in a list of new names. Individual names can be modified as well.  For example, I think `'state'` makes a lot more sense as a column name than `'stabbr'`, so I'm going to change it.  The following code gets this done: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>instnm</th>\n",
       "      <th>state</th>\n",
       "      <th>preddeg</th>\n",
       "      <th>control</th>\n",
       "      <th>ugds</th>\n",
       "      <th>ug25abv</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Alabama A &amp; M University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4051.0</td>\n",
       "      <td>0.1049</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>University of Alabama at Birmingham</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>11200.0</td>\n",
       "      <td>0.2422</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Amridge University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>322.0</td>\n",
       "      <td>0.8540</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>University of Alabama in Huntsville</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5525.0</td>\n",
       "      <td>0.2640</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Alabama State University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5354.0</td>\n",
       "      <td>0.1270</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                instnm state  preddeg  control     ugds  \\\n",
       "0             Alabama A & M University    AL        3        1   4051.0   \n",
       "1  University of Alabama at Birmingham    AL        3        1  11200.0   \n",
       "2                   Amridge University    AL        3        2    322.0   \n",
       "3  University of Alabama in Huntsville    AL        3        1   5525.0   \n",
       "4             Alabama State University    AL        3        1   5354.0   \n",
       "\n",
       "   ug25abv  \n",
       "0   0.1049  \n",
       "1   0.2422  \n",
       "2   0.8540  \n",
       "3   0.2640  \n",
       "4   0.1270  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.rename(columns = {'stabbr': 'state'}, inplace = True)\n",
    "\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "The `rename` data frame method takes a __dictionary__ of values, which is a Python data structure that we haven't discussed yet.  Dictionaries, or `dict` objects as they are often called, are enclosed by curly braces (`{}`) and are made up of __key/value pairs__.  In turn, they can really come in handy when working with paired values; we may return to them later in the semester.  As this is a minor change to the data frame, the argument `inplace = True` makes sense to modify `df` directly.  \n",
    "\n",
    "We now want to subset our data even further.  As mentioned above, the `preddeg` column designates the primary degree granted by the colleges and universities in the dataset.  We're interested in comparing primarily bachelor's-granting universities, which have the code of `3`; in turn, we want to tell `pandas` to keep only those rows where the `preddeg` column is equal to 3.  First, let's check our `dtypes` to see if the columns is formatted as a string or number: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "instnm      object\n",
       "state       object\n",
       "preddeg      int64\n",
       "control      int64\n",
       "ugds       float64\n",
       "ug25abv    float64\n",
       "dtype: object"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.dtypes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "It appears as though `preddeg` is an integer, so we will work with the values as numbers.  To subset rows in __pandas__, we can use the `.query()` data frame method.  `.query()` requires an __expression__ to be evaluated by Python using boolean/logical operators.  The method will then return those rows for which the result of the expression is `True`, and drop those rows that return `False`.  In this case, we want all of those rows for which the value of the `preddeg` column is equal to 3.  Let's create a new dataframe, called `ug` for undergraduate, from this expression.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(2133, 6)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ug = df.query('preddeg == 3')\n",
    "\n",
    "ug.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "We've gone from over 7800 colleges & universities down to 2133.  \n",
    "\n",
    "### Missing data\n",
    "\n",
    "As mentioned in class, missing data in __pandas__ are designated with the value `NaN`, which refers to \"not a number.\"  Data analysts need to take missing data seriously, as they could be representative of a systematic flaw in the dataset.  You can check for missing values in __pandas__ with the `.isnull()` method.  Let's query our undergradate data frame to check to see how many rows in our dataset have missing values for our column of interest, `ug25abv`.  Note that we need the argument `engine = 'python'` in this example as we are using a __pandas__ method within the expression.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(36, 6)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ugnull = ug.query('ug25abv.isnull()', engine = 'python')\n",
    "\n",
    "ugnull.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "It looks like we have 36 rows with null values for the `ug25abv` column.  Let's see what universities they are: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>instnm</th>\n",
       "      <th>state</th>\n",
       "      <th>preddeg</th>\n",
       "      <th>control</th>\n",
       "      <th>ugds</th>\n",
       "      <th>ug25abv</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>104</th>\n",
       "      <td>Frank Lloyd Wright School of Architecture</td>\n",
       "      <td>AZ</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>2.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>649</th>\n",
       "      <td>Yeshiva Ohr Elchonon Chabad West Coast Talmudi...</td>\n",
       "      <td>CA</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>131.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1171</th>\n",
       "      <td>Rosalind Franklin University of Medicine and S...</td>\n",
       "      <td>IL</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1893</th>\n",
       "      <td>New England College of Optometry</td>\n",
       "      <td>MA</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2598</th>\n",
       "      <td>Beth Hatalmud Rabbinical College</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>47.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2599</th>\n",
       "      <td>Beth Hamedrash Shaarei Yosher Institute</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>49.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2668</th>\n",
       "      <td>Yeshiva of Far Rockaway Derech Ayson Rabbinica...</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>57.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2707</th>\n",
       "      <td>Kehilath Yakov Rabbinical Seminary</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>120.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2733</th>\n",
       "      <td>Mesivtha Tifereth Jerusalem of America</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>50.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2767</th>\n",
       "      <td>Ohr Hameir Theological Seminary</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>94.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2781</th>\n",
       "      <td>Rabbinical College Bobover Yeshiva Bnei Zion</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>263.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2783</th>\n",
       "      <td>Rabbinical College Beth Shraga</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>47.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2785</th>\n",
       "      <td>Rabbinical College of Long Island</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>96.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2854</th>\n",
       "      <td>Talmudical Seminary Oholei Torah</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>321.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2859</th>\n",
       "      <td>Torah Temimah Talmudical Seminary</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>180.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2875</th>\n",
       "      <td>Webb Institute</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>82.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2882</th>\n",
       "      <td>Yeshiva Karlin Stolin</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>119.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2886</th>\n",
       "      <td>Yeshiva Shaar Hatorah</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>75.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2889</th>\n",
       "      <td>Yeshivath Zichron Moshe</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>164.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2920</th>\n",
       "      <td>Davidson College</td>\n",
       "      <td>NC</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>1782.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3657</th>\n",
       "      <td>Talmudical Yeshiva of Philadelphia</td>\n",
       "      <td>PA</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>124.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4268</th>\n",
       "      <td>Washington and Lee University</td>\n",
       "      <td>VA</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>1846.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4622</th>\n",
       "      <td>Yeshiva Gedolah of Greater Detroit</td>\n",
       "      <td>MI</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>58.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4915</th>\n",
       "      <td>Yeshiva Gedolah Imrei Yosef D'spinka</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>125.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5066</th>\n",
       "      <td>Rabbi Jacob Joseph School</td>\n",
       "      <td>NJ</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>86.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5152</th>\n",
       "      <td>Yeshivas Novominsk</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>118.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5433</th>\n",
       "      <td>Yeshiva D'monsey Rabbinical College</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>48.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5585</th>\n",
       "      <td>Yeshiva of the Telshe Alumni</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>126.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5884</th>\n",
       "      <td>Yeshiva Shaarei Torah of Rockland</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>116.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5909</th>\n",
       "      <td>Franklin W Olin College of Engineering</td>\n",
       "      <td>MA</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>343.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6061</th>\n",
       "      <td>Beis Medrash Heichal Dovid</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>96.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6566</th>\n",
       "      <td>Yeshiva Toras Chaim</td>\n",
       "      <td>NJ</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>182.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7520</th>\n",
       "      <td>Bais HaMedrash and Mesivta of Baltimore</td>\n",
       "      <td>MD</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>49.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7528</th>\n",
       "      <td>Be'er Yaakov Talmudic Seminary</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>298.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7689</th>\n",
       "      <td>Yeshiva Gedolah Kesser Torah</td>\n",
       "      <td>NY</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>70.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7691</th>\n",
       "      <td>Yeshiva Yesodei Hatorah</td>\n",
       "      <td>NJ</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>64.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                 instnm state  preddeg  \\\n",
       "104           Frank Lloyd Wright School of Architecture    AZ        3   \n",
       "649   Yeshiva Ohr Elchonon Chabad West Coast Talmudi...    CA        3   \n",
       "1171  Rosalind Franklin University of Medicine and S...    IL        3   \n",
       "1893                   New England College of Optometry    MA        3   \n",
       "2598                   Beth Hatalmud Rabbinical College    NY        3   \n",
       "2599            Beth Hamedrash Shaarei Yosher Institute    NY        3   \n",
       "2668  Yeshiva of Far Rockaway Derech Ayson Rabbinica...    NY        3   \n",
       "2707                 Kehilath Yakov Rabbinical Seminary    NY        3   \n",
       "2733             Mesivtha Tifereth Jerusalem of America    NY        3   \n",
       "2767                    Ohr Hameir Theological Seminary    NY        3   \n",
       "2781       Rabbinical College Bobover Yeshiva Bnei Zion    NY        3   \n",
       "2783                     Rabbinical College Beth Shraga    NY        3   \n",
       "2785                  Rabbinical College of Long Island    NY        3   \n",
       "2854                   Talmudical Seminary Oholei Torah    NY        3   \n",
       "2859                  Torah Temimah Talmudical Seminary    NY        3   \n",
       "2875                                     Webb Institute    NY        3   \n",
       "2882                              Yeshiva Karlin Stolin    NY        3   \n",
       "2886                              Yeshiva Shaar Hatorah    NY        3   \n",
       "2889                            Yeshivath Zichron Moshe    NY        3   \n",
       "2920                                   Davidson College    NC        3   \n",
       "3657                 Talmudical Yeshiva of Philadelphia    PA        3   \n",
       "4268                      Washington and Lee University    VA        3   \n",
       "4622                 Yeshiva Gedolah of Greater Detroit    MI        3   \n",
       "4915               Yeshiva Gedolah Imrei Yosef D'spinka    NY        3   \n",
       "5066                          Rabbi Jacob Joseph School    NJ        3   \n",
       "5152                                 Yeshivas Novominsk    NY        3   \n",
       "5433                Yeshiva D'monsey Rabbinical College    NY        3   \n",
       "5585                       Yeshiva of the Telshe Alumni    NY        3   \n",
       "5884                  Yeshiva Shaarei Torah of Rockland    NY        3   \n",
       "5909             Franklin W Olin College of Engineering    MA        3   \n",
       "6061                         Beis Medrash Heichal Dovid    NY        3   \n",
       "6566                                Yeshiva Toras Chaim    NJ        3   \n",
       "7520            Bais HaMedrash and Mesivta of Baltimore    MD        3   \n",
       "7528                     Be'er Yaakov Talmudic Seminary    NY        3   \n",
       "7689                       Yeshiva Gedolah Kesser Torah    NY        3   \n",
       "7691                            Yeshiva Yesodei Hatorah    NJ        3   \n",
       "\n",
       "      control    ugds  ug25abv  \n",
       "104         2     2.0      NaN  \n",
       "649         2   131.0      NaN  \n",
       "1171        2     NaN      NaN  \n",
       "1893        2     NaN      NaN  \n",
       "2598        2    47.0      NaN  \n",
       "2599        2    49.0      NaN  \n",
       "2668        2    57.0      NaN  \n",
       "2707        2   120.0      NaN  \n",
       "2733        2    50.0      NaN  \n",
       "2767        2    94.0      NaN  \n",
       "2781        2   263.0      NaN  \n",
       "2783        2    47.0      NaN  \n",
       "2785        2    96.0      NaN  \n",
       "2854        2   321.0      NaN  \n",
       "2859        2   180.0      NaN  \n",
       "2875        2    82.0      NaN  \n",
       "2882        2   119.0      NaN  \n",
       "2886        2    75.0      NaN  \n",
       "2889        2   164.0      NaN  \n",
       "2920        2  1782.0      NaN  \n",
       "3657        2   124.0      NaN  \n",
       "4268        2  1846.0      NaN  \n",
       "4622        2    58.0      NaN  \n",
       "4915        2   125.0      NaN  \n",
       "5066        2    86.0      NaN  \n",
       "5152        2   118.0      NaN  \n",
       "5433        2    48.0      NaN  \n",
       "5585        2   126.0      NaN  \n",
       "5884        2   116.0      NaN  \n",
       "5909        2   343.0      NaN  \n",
       "6061        2    96.0      NaN  \n",
       "6566        2   182.0      NaN  \n",
       "7520        2    49.0      NaN  \n",
       "7528        2   298.0      NaN  \n",
       "7689        2    70.0      NaN  \n",
       "7691        2    64.0      NaN  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ugnull"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "Interesting!  It looks like primarily we've returned Jewish _yeshivas_ as well as a few graduate colleges - e.g. New England College of Optometry which is designated as a bachelor's granting institution but has no undergraduates - that appear to have been mis-coded.  This reveals why careful inspection of your data is important.  Missing values often aren't at random, but instead may be systematic or clustered amongst a particular type of record in your dataset.  Additionally, it looks like we've uncovered some errors in the original data as well.   \n",
    "\n",
    "In this case, let's go ahead and drop these records from our dataset.  Recall from class that there are other options available for missing data as well; you can __fill__ the `NaN` values with some other value if appropriate, such as `0` or the mean of the column.  We'll drop these values using `.query()` again, but instead with the `.notna()` __pandas__ method.  This method will identify the rows that are _not_ `NaN` in the `ug25abv` column, and retain only those rows.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(2097, 6)"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ug_nonull = ug.query('ug25abv.notna()', engine = 'python')\n",
    "\n",
    "ug_nonull.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "I can also use the `.dropna()` method, which we learned in class.  One thing about `.dropna()` is that it can operate over the entire data frame if we don't pass it any columns.  Let's check it out: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(2097, 6)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ug1 = ug.dropna()\n",
    "\n",
    "ug1.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "We get the same result, which means that we've also removed any of the rows that had null values in the `ugds` column.  \n",
    "\n",
    "## Group-wise data analysis\n",
    "\n",
    "In class, we discussed the \"split-apply-combine\" approach to data analysis.  This approach involves: \n",
    "\n",
    "1. __Splitting__ the data into groups based on some common characteristic; \n",
    "2. __Applying__ some function to each group; \n",
    "3. __Combining__ the results back into a single dataset, allowing for group-wise comparisons.  \n",
    "\n",
    "While conceptually simple, this approach to data analysis is extraordinarily powerful - and extremely common!  For example, let's say you are working in analytics for a business, and your supervisor wants you to compare sales results of your stores by region.  This is a very common request when working with data professionally, and fortunately __pandas__ can help you out with this with minimal code.  \n",
    "\n",
    "Let's start with a couple guiding questions.  \n",
    "\n",
    "1. How does the proportion of undergraduate students above age 25 vary among public, private, and for-profit universities?  \n",
    "2. How does the proportion of undergraduate students above age 25 vary by institution size?  \n",
    "\n",
    "Group-wise data analysis in __pandas__ is conducted by the creation of a `groupby` object, in which you tell `pandas` exactly how the data should be grouped.  To address the first question, we'll want to group our data by institution type.  The code below gets this done.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.groupby.generic.DataFrameGroupBy"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "groups1 = ug1.groupby('control')\n",
    "\n",
    "type(groups1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "Notice the object `type` - Python knows that our new object, `groups1`, represents grouped data.  In turn, anything we calculate over this object will apply to each of our groups, defined by unique values of the `control` column.  Let's try:  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "control\n",
       "1    0.215475\n",
       "2    0.235799\n",
       "3    0.670381\n",
       "Name: ug25abv, dtype: float64"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "groups1.ug25abv.mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "Interesting stuff - whereas public and private non-profit colleges and universities are pretty similar in their age composition, with 21-23 percent of their student bodies above age 25, private for-profit universities are quite different, with over two-thirds of their undergraduates aged over 25.  \n",
    "\n",
    "Let's visualize this with `seaborn`, which also sets up very nicely for group-wise data analysis.  Box plots and violin plots, for example, can be extended by specifying a second column for the plot.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='control', ylabel='ug25abv'>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEJCAYAAACUk1DVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAeVElEQVR4nO3dfVRUdcIH8O+dARUekgGcgSxDt3WxFKyndBEQOrmBi44rppVarI+G6drhWSzDfG1rfcG1tbJdT3pMjwklRxTlqQhfEkIw0z1p+RqaiiIw8lYoyDB3nj8YJyd1GHDu3DvM9/NP/riXyxeG+M59+13BbDabQUREHk8ldwAiIlIGFgIREQFgIRARkQULgYiIALAQiIjIgoVAREQAWAhERGThJXeAu1FXdxWiyNsoiIgcoVIJCAj4rzsud+tCEEUzC4GIyEl4yIiIiACwEIiIyIKFQEREACQuhMbGRowePRoXL168ZdmJEycwbtw4JCQkYP78+WhtbZUyCpFLHDt2FNOmTcbx49/LHYWowyQrhCNHjmDixIk4d+7cbZfPmTMHixYtwhdffAGz2Yzs7GypohC5zJo1q2E2m/Hvf78rdxSiDpOsELKzs7F48WLodLpbll26dAnNzc145JFHAADjxo1Dfn6+VFGIXOLYsaO4du0qAODatavcSyC3I9llp0uWLLnjsurqami1WutYq9WiqqpKqihELrFmzWqb8b///S7ef3+dTGm6tv37i1BcXOj07TY01AMA/P01Tt92TEwcoqNjnb5dZ5LlPgRRFCEIgnVsNpttxo4KCvJzZiyiu3Jj7+DmsVZ7j0xpuraePX3g7a12+nZ/+qkBANCrV5DTt92zp4/ifx9kKYSQkBAYDAbr+MqVK7c9tNSemppG3phGiuHj44OmpiabscHws4yJuq7w8CEIDx/i9O1mZLwFAJg9e57Ttw1A9t8HlUqw+0ZalstO77vvPnTv3h2HDx8GAOzYsQOxscrelSJqT79+D9qMf/Ob38qUhKhzXFoIKSkp+O677wAAK1euxLJlyzBy5Ehcu3YNycnJroxC5HQnTx63GZ84cUymJESdI/kho71791r/vW7dLyfYBgwYgK1bt0r95YlcRhRFu2MipeOdykREBICFQEREFiwEIiICwEIgIiILFgKRk6hUKrtjIqXjbyyRk2i1wXbHRErHQiBykrq6WrtjIqVjIRA5SUBAgN0xkdKxEIic5Ncz9nIGX3I3LAQFqq+vw/Llb1qn4iV38euJFjnxIrkXFoIC5eVtxw8/nMLOndvkjkJEHoSFoDD19XUoKvoSZrMZRUX7uJdARC7DQlCYvLztMJlMAACTqZV7CUTkMiwEhdm//yu7Y1Kunj17/mrsL1MSos5hISiMl5eX3TEp18SJts/0mDz5zzIlIeocFoLC3O65vOQecnNtn++xbVu2TEmIOoeFoDDBwffajENC7r3DmqQ0VVWVdsdESsdCUJg+fR6wGd9//wN3WJOIyLlYCArz/fdH7Y6JiKTCQlCYYcOiIQhtL4sgqDBsWLTMiYjIU7AQFEavT4KXlxpA2xVGY8aMkzkREXkKFoLCaDQBiImJgyAIiImJhb+/Ru5IROQheJG7AsXFPYkDB/bjiSdGyB2FqNOysjahvPy83DE65MKFtrwZGW/JnKRj+vQJxaRJye2v2A4WggIVFu5Fc3Mz9u3bgxdemCp3HKJOKS8/j9NnT0Ht303uKA4T1W3Txpyp+VHmJI4zNbQ4bVssBIWpr69DcXEhzGYziouLMGbMOB42Irel9u8G/9jecsfo0hqKKpy2LZ5DUJi8vO0QRREAIIomTm5HRC7DQlCY0tL9N812akJp6X6ZE5GjunfvbndMpHQsBIX57/9+3Gb82GNDZEpCHXX9+vVfjZ13bJfIFVgICmc28zGM7ouvHbkXFoLC/Oc/h+yOiYikwkJQmGHDoqFStd2prFKpOXUFEbmMpIWQl5eHxMRExMfHIzMz85blx44dw9NPP40xY8bgpZdewk8//SRlHLeg1ydBrW57WdRqNaeuICKXkawQqqqqsGrVKmRlZSE3NxdbtmxBWVmZzTpLlixBamoqdu7ciX79+mH9+vVSxXEbGk0AhgyJBAAMGfJ73oNARC4jWSGUlJQgMjISGo0Gvr6+SEhIQH5+vs06oiji6tW2J4I1NTWhR48eUsVxS4IgyB2BiDyIZHcqV1dXQ6vVWsc6nQ5Hj9rO7T937lxMnToVS5cuhY+PD7KzO/bIwaAgP6dkVZLa2lp8880BAMDBg6V46aUXERAQIHMq6iyt9h65I8jG21stdwSP4e2tdsrvmmSFIIqizTtcs9lsM25ubsb8+fOxceNGREREYMOGDUhPT8fatWsd/ho1NY0Qxa51ad9HH32E1tZWAEBrays2bNjE+YzcmMHws9wRZGM0muSO4DGMRpNDv2sqlWD3jbRkh4xCQkJgMBisY4PBAJ1OZx2fPn0a3bt3R0REBADg2WefxcGDB6WK4zZKSoqt9x6YzWaUlBTLnIiIPIVkhRAVFYXS0lLU1taiqakJBQUFiI2NtS4PDQ1FZWUlzp49CwDYs2cPwsPDpYrjNoKCguyOiYikItkho+DgYKSlpSE5ORlGoxHjx49HREQEUlJSkJqaivDwcCxbtgx//etfYTabERQUhKVLl0oVx23U1NTYHRMRSUXS6a/1ej30er3Nx9atW2f9d1xcHOLi4qSM4HYee2wISkq+so4ff3yojGmIyJPweQgKx7mMpLF/fxGKiwuduk0fHx80NTXZjJ395K2YmDhER8e2vyJRJ7AQFIZzGbmv3r3vw5kzZTeN75cxjfwaGurRWn/dqQ9woVu11l9Hg1e9U7bFQlCYYcOiUVj4JUTRxLmMJBQdHSvJO+1Zs6ahqakJAwcOwiuvzHP69omkxEJQGL0+CV99tQ+iCKhUKs5l5GZ6974PFRWX8OKLf5E7iuz8/TW40lrHR2hKrKGowmlT3HC2U4XRaAKg0wUDaLu7m3MZuRcvL2888EBfvm7kllgIClNfX4fq6moAgMFQjYaGenkDEZHHYCEoTF7eduuVRaJoxs6d22RORESegoWgMKWl+2Eytc1lZDK1orR0v8yJiMhTsBAUZtiwaOskgIIg8CojInIZFoLCxMU9aTO53RNPjJA5ERF5ChaCwhQW7gVwY5pwAfv27ZEzDhF5EBaCwrSdM7gxXYWZ5xCIyGVYCAozbFg01Oq2+wXVai+eQyAil2EhKIxenwSVqu2QEe9UJiJXYiEojEYTgJiYOAiCgJiYWN7xSkQuw7mMFEivT8KlSxe5d0BELsVCUCCNJgBz5y6SOwYReRgeMiIiIgAsBEW6cOEcZs2ahvLy83JHISIPwkJQoDVrVqOpqQlr1rwndxQi8iAsBIW5cOEcqqouAwAqKy9zL4GIXIaFoDBr1qz+1Zh7CUTkGiwEhbmxd3BDZeXlO6xJRORcLAQiIgLAQlCcwMAgm3FQUNAd1iQici4WgsLU19fZjOvq6u6wJhGRc7EQFObG09LuNCYikgoLQWEeffRxu2MiIqmwEBSmW7duNuPu3bvLlISIPA0LQWH+859DNuPDh7+RKQkReRpJCyEvLw+JiYmIj49HZmbmLcvPnj2LF154AWPGjMG0adPQ0NAgZRy30PbENDUAQK1W84lpROQykhVCVVUVVq1ahaysLOTm5mLLli0oKyuzLjebzZg5cyZSUlKwc+dOPPTQQ1i7dq1UcdyGXp8Es7ntmcpms5nPRCAil5GsEEpKShAZGQmNRgNfX18kJCQgPz/fuvzYsWPw9fVFbGwsAGDGjBmYPHmyVHHcyi+FIHMQIvIoDhXC22+/jfLy8g5tuLq6Glqt1jrW6XSoqqqyji9cuIBevXph3rx5SEpKwuLFi+Hr69uhr9EV5eVthyC0vSyCIGDnzm0yJyIiT+HQE9PMZjMmT56Mfv36YcKECYiPj7/laphfE0XR5hp6s9lsM25tbcXBgwexefNmhIeH45133sHy5cuxfPlyh8MHBfk5vK67OHBgP0TRBAAQRRMOHNiP2bP/V+ZU5Chv77bzP1rtPTInkd+NnwVJz9tb7ZTfOYcK4dVXX8Xs2bPx1VdfYdu2bVixYgVGjhyJefPm3fFzQkJCcOjQL1fMGAwG6HQ661ir1SI0NBTh4eEAgNGjRyM1NbVD4WtqGiGKXeu4SmRkNL78crfN2GD4WcZE1BFGY1uZ8zX75WdB0jMaTQ79zqlUgt030g6fQ1CpVBg0aBAeffRRBAQE2Pyxv52oqCiUlpaitrYWTU1NKCgosJ4vAIBHH30UtbW1OHnyJABg7969GDhwoKNxuqy4uCdtxk88MUKmJETkaRzaQ9i1axdycnLw7bffYuTIkVi6dGm7f7yDg4ORlpaG5ORkGI1GjB8/HhEREUhJSUFqairCw8Pxr3/9CwsWLEBTUxNCQkKwYsUKp3xT7uz//m/HLeOZMzu250RE1BkOFcL69evxzDPP4J133kGPHj0c3rher4der7f52Lp166z/Hjx4MLZu3erw9jzBoUNf24y/+eYAC4GIXMKhQvjkk09QV1eH0tJSeHl5YfDgwejZs6fU2YiIyIUcKoTCwkK89tpr+O1vfwtRFFFeXo5Vq1ZhyJAhUucjIiIXcagQ3nnnHWzevBn9+/cH0HZT2cKFC7FtG6+RJyLqKhy6ykgQBGsZAMDAgQOtd9MSEVHXYLcQ6uvrUV9fj0GDBmH9+vW4evUqmpqakJmZicjISFdlJCIiF7B7yCgyMhKCIFj3Bv7xj39YlwmCgPT0dGnTKdz+/UUoLi6U/OtkZLzltG3FxMQhOjq2/RWJyOPYLYQbN42R69xcwG1jPrKCiFzDoZPKLS0tKCwsxNWrVwEAJpMJFy5cQFpamqThlC46Otbp77aPHTuKt9/+ZT6nV16Zi4cfHuTUr0FEdDsOFUJaWhrKy8thMBjw8MMP48iRIxg6dKjU2TzSwIER1r0EHx9flgERuYxDxyNOnDiBbdu2YcSIEZg3bx4+/vhjPt1MQr173wcAmDXrr/IGISKP4lAh6HQ6eHl5oW/fvjh9+jT69++Pn3/mbI5S8fO7B2FhD3HvgIhcyqFC8PX1RV5eHgYMGIDPP/8cp06dwrVr16TORkRELuRQISxatAgnTpxAdHQ0VCoVnn/+eUybNk3qbERE5EIOnVTu27cvXnvtNVy6dAmvv/46goODpc5FREQuZncPYfbs2QCAy5cvIykpCWPGjEFiYiImT56MyspKlwQkIiLXsFsIP/74IwBg2bJlGD16NA4fPoxDhw4hMTER8+fPd0lAIiJyDYfOIfz444/WcwaCIHAPgYioC2p3crtvv/0WoaGhOHPmjPXjly5dgkrFKRWIiLoSuyeVExISsHz5chw/fhxGoxEffPABduzYgaVLl2LBggWuykhEbsrU0IKGogq5YzhMbDYBAFQ91DIncZypoQUIcs627BbC3LlzAbTNZVRdXQ0AGDBggM3DcoiIbqdPn1C5I3TYhQvnAQAPBLlR9iDn/azbvez066+/RkVFBYYPHw4ACAsLAwBs374dSUlJTglBRF3PpEnJckfosBtTzaenL5Q5iTzsngjYsGEDFi1ahM8//xyJiYk4cOCAddmmTZskD0dERK5jdw8hJycHOTk58PPzQ1FREdLS0rBx40aEhYXxEZpERF2M3T0ELy8v+Pn5AQBiY2ORnp6OmTNnoq6uDoIguCQgERG5ht1CCAwMRE5ODq5fvw4AGDt2LOLj4zF9+nTOdkpE1MXYLYQ33ngDOTk5+Oyzz6wfmzt3Lh5//HFUVLjPpWRERNQ+u+cQHnjgAWRlZd3y8fT0dEyZMkWqTEREJAOHZjt9/fXXbcaCIMDHxwf9+/fHhAkToFa7z00cRER0ew7PP3H8+HGEhYXhoYceQllZGSoqKlBcXIylS5dKmY+IiFzEoT2EM2fOIDMz03rF0YQJEzB16lRkZWVh9OjRkgYkIiLXcGgPoaGhwVoGANCjRw80NjZCEAR4e3tLFo6IiFzHoUIYPHgwXn31VRw4cAClpaWYM2cOIiIiUFhYCB8fnzt+Xl5eHhITExEfH4/MzMw7rrdv3z48+eSTHU9PRERO41AhvPnmm+jduzeWL1+OlStXIjQ0FIsXL8bVq1fx5ptv3vZzqqqqsGrVKmRlZSE3NxdbtmxBWVnZLetduXIFGRkZd/ddEBHRXXOoECZMmIAvv/wSJpMJLS0tKCgowNNPP43ExMQ7znpaUlKCyMhIaDQa+Pr6IiEhAfn5+best2DBArz88st3910QEdFdc+ik8sKFv8z8ZzQa8emnn6JPnz52P6e6uhpardY61ul0OHr0qM06mzZtwsMPP4zBgwd3JLNVUJBf+yu5IW/vtst4tdp7ZE5CHcXXzr15+uvnUCEMHTrUZhwVFYXnnnsOM2fOvOPniKJoM9+R2Wy2GZ8+fRoFBQXYuHFjpx/HWVPTCFHsepPsGY1tD+kwGDg9iLvha+feuvrrp1IJdt9Id+o5mHV1ddYH5txJSEgIDAaDdWwwGKDT6azj/Px8GAwGPP3005g+fTqqq6sxadKkzsQhIiIncGgPQa/X24wrKirw7LPP2v2cqKgorF69GrW1tfDx8UFBQQHeeust6/LU1FSkpqYCAC5evIjk5OTbTpNBRESu0eFzCIIgIDAwEA8++KDdzwkODkZaWhqSk5NhNBoxfvx4REREICUlBampqQgPD7+75ERE5FSdOofgKL1ef8vexbp1625Z7/7778fevXs79TWIiMg5OnUOgYiIuh4WAhERAWAhEBGRBQuBiIgAsBCIiMjCoauMiOSSlbUJ5eXn5Y7hsAsX2rJmZLzVzprK0qdPKCZNSpY7BsmMhUCKVl5+HufKTiLEzz1+VX0hAgCaK2+d2VepKhtb5Y5ACuEe/5eRRwvx88L/RATKHaPL2nC0Vu4IpBA8h0BERABYCEREZMFCICIiACwEIiKyYCEQEREAFgIREVmwEIiICAALgYiILFgIREQEgIVAREQWXX7qCnebHA3gBGlEJI8uXwjl5edx6ocyqHto5I7iMNGkBgCUlV+ROYnjTM31ckcgorvU5QsBANQ9NPANHSF3jC7t2vk9ckcgorvEcwhERASAhUBERBYsBCIiAsBCICIiCxYCEREBYCEQEZEFC4GIiACwEIiIyIKFQEREACQuhLy8PCQmJiI+Ph6ZmZm3LN+9ezf+9Kc/YcyYMfjLX/6ChoYGKeMQEZEdkhVCVVUVVq1ahaysLOTm5mLLli0oKyuzLm9sbMQbb7yBtWvXYufOnQgLC8Pq1aulikNERO2QrBBKSkoQGRkJjUYDX19fJCQkID8/37rcaDRi8eLFCA4OBgCEhYXh8uXLUsUhIqJ2SFYI1dXV0Gq11rFOp0NVVZV1HBAQgKeeegoA0NzcjLVr1+IPf/iDVHGIiKgdks12KooiBEGwjs1ms834hp9//hmzZs3CgAEDkJSU1KGvERTk1+463t7qDm2TOs/bWw2t9h6nb7PZqVuk25HitXNHN/5eeOrPQrJCCAkJwaFDh6xjg8EAnU5ns051dTWmTZuGyMhIzJs3r8Nfo6amEaJotrvOlSs1MDXXc3pmiZma63HlihcMhp+dul2j0eTU7dHtGY0mp7927ujG71tX/VmoVILdN9KSHTKKiopCaWkpamtr0dTUhIKCAsTGxlqXm0wmzJgxA3/84x8xf/782+49EBGR60i2hxAcHIy0tDQkJyfDaDRi/PjxiIiIQEpKClJTU1FZWYnjx4/DZDLhiy++AAAMGjQIS5YscWoOf38NDD+18gE5Ert2fg/8/TVyxyCiuyDpE9P0ej30er3Nx9atWwcACA8Px8mTJ6X88kRE1AG8U5mIiACwEIiIyIKFQEREACQ+h0B0txoa6lHX2IoNR2vljtJlVTa2IqChXu4YpADcQyAiIgDcQyCF8/fXoHvTFfxPRKDcUbqsDUdr0YOXDBO4h0BERBYsBCIiAsBCICIiCxYCEREBYCEQEZEFC4GIiACwEIiIyIKFQEREADzkxjR3e2Ka2Nr20EiVVw+ZkzjO1FwPoJfcMYjoLnT5QujTJ1TuCB124cJ5AMADfdzpD2wvt/xZE9EvunwhTJqULHeEDsvIeAsAkJ6+UOYkRORJeA6BiIgAsBCIiMiChUBERABYCEREZMFCICIiACwEIiKy6PKXnZL7q3SjZyo3togAAL9u7vNeq7KxFX3lDkGKwEIgRXO3m92qLTcV9gpxn9x94X4/Z5IGC4EUzd1uLORNheTOWAhE5Hb27y9CcXGh07d7Y9qYG8XuTDExcYiOjnX6dp2JhUBEZOHv7y93BFmxEIjI7URHxyr+3bY7cp9LIYiISFIsBCIiAiBxIeTl5SExMRHx8fHIzMy8ZfmJEycwbtw4JCQkYP78+WhtbZUyDhER2SFZIVRVVWHVqlXIyspCbm4utmzZgrKyMpt15syZg0WLFuGLL76A2WxGdna2VHGIiKgdkp1ULikpQWRkJDQaDQAgISEB+fn5ePnllwEAly5dQnNzMx555BEAwLhx4/Dee+9h0qRJUkVyOl765r7c8bUD+PqRtCQrhOrqami1WutYp9Ph6NGjd1yu1WpRVVUlVRy34umXvrkzvnbkziQrBFEUIQiCdWw2m23G7S13RFCQ390HvQtjx47C2LGjZM1AncPXjuhWkhVCSEgIDh06ZB0bDAbodDqb5QaDwTq+cuWKzXJH1NQ0QhTNdx+WiMgDqFSC3TfSkp1UjoqKQmlpKWpra9HU1ISCggLExv5y7PO+++5D9+7dcfjwYQDAjh07bJYTEZFrSVYIwcHBSEtLQ3JyMsaOHYvRo0cjIiICKSkp+O677wAAK1euxLJlyzBy5Ehcu3YNycnuNZEZEVFXIpjNZrc95sJDRkREjpPtkBEREbkXFgIREQFgIRARkYVbT3+tUnXsvgUiIk/W3t9Mtz6pTEREzsNDRkREBICFQEREFiwEIiICwEIgIiILFgIREQFgIRARkQULgYiIALAQiIjIgoVAREQAWAiK1NjYiNGjR+PixYtyR6EOev/99zFq1CiMGjUKK1askDsOdcC7776LxMREjBo1Chs2bJA7jixYCApz5MgRTJw4EefOnZM7CnVQSUkJiouLsX37duTm5uLYsWPYtWuX3LHIAQcPHsSBAwewc+dO5OTk4KOPPsLZs2fljuVyLASFyc7OxuLFizv8fGmSn1arxdy5c9GtWzd4e3vjwQcfREVFhdyxyAFDhw7Fpk2b4OXlhZqaGphMJvj6+sody+XcerbTrmjJkiVyR6BO6t+/v/Xf586dw+eff46PP/5YxkTUEd7e3njvvffw4YcfYuTIkQgODpY7kstxD4HIyX744QdMnToVr732Gvr27St3HOqA1NRUlJaW4vLly8jOzpY7jsuxEIic6PDhw5gyZQpeeeUVJCUlyR2HHHTmzBmcOHECAODj44P4+HicOnVK5lSux0IgcpLLly9j1qxZWLlyJUaNGiV3HOqAixcvYsGCBWhpaUFLSwv27NmDxx57TO5YLsdzCEROsn79ely/fh3Lly+3fuy5557DxIkTZUxFjoiLi8PRo0cxduxYqNVqxMfHe2Sp84lpREQEgIeMiIjIgoVAREQAWAhERGTBQiAiIgAsBCIismAhEEnk/fffx+7duzv8eWFhYaitrZUgEZF9LAQiiXz99ddobW2VOwaRw3hjGtGvbN26FRs2bIBKpUJAQAAyMjJQVFSEjz76CCqVCr169cLChQvRr18/zJ07F35+fjh16hQqKysRFhaGjIwM5Obm4vvvv8eKFSugVquxZ88e1NfXo7y8HE888QRmzJiBv/3tbzh58iQEQcDw4cMxe/ZseHnxf0mSD3/7iG5y8uRJrFy5Etu3b8e9996LjRs3YsqUKRBFEVu2bEFgYCC2bduGWbNm4dNPPwUAfP/999i0aRMEQcAzzzyD/Px8TJ482frfp556Cnv27EFzc7P1c9LT06HRaJCXlwej0YiZM2fiww8/xPTp0+X89snD8ZAR0U1KS0sRExODe++9FwAwZcoUjBgxAomJiQgMDAQAjBs3DlVVVdYn2g0fPtz6DITf/e53aGhouO22b54bp6ioCM8//zwEQUC3bt3w3HPPoaioSOLvjsg+FgLRTdRqNQRBsI6bm5tRXl5+y3pms9l6fqBHjx7WjwuCgDvNBnPzA1dEUbT5OqIo8nwDyY6FQHST3//+9ygtLUV1dTUA4JNPPkFhYSE+++wz65U/OTk50Gg0CA0NtbsttVp9xz/yMTEx2Lx5M8xmM1paWpCdnY2oqCjnfjNEHcRzCEQ3CQsLw5w5c/Diiy8CaHss5q5du7B79278+c9/hiiKCAwMxAcffACVyv77qSeffBL//Oc/YTQab1m2YMEC/P3vf4der4fRaMTw4cMxY8YMSb4nIkdxtlMiIgLAQ0ZERGTBQiAiIgAsBCIismAhEBERABYCERFZsBCIiAgAC4GIiCxYCEREBAD4f6HtGQmihL7bAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import seaborn as sns\n",
    "sns.set(style = \"darkgrid\")\n",
    "\n",
    "sns.boxplot(data = ug1, x = 'control', y = 'ug25abv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "We can get a clear sense here of the variations of the distributions between the groups.  While the mean percentage of undergraduates above 25 at public non-profits was below that of private non-profits, it appears as though the median for private non-profits is lower; there is simply a longer tail of private non-profits with large proportions of their student bodies above 25.  For private for-profits, the distribution is quite evident - with the median value above 70 percent.  \n",
    "\n",
    "`seaborn` includes even more functionality for group-wise visualization - incorporating, in some instances, multiple groups!  We'll learn more about this in a couple weeks when we focus on data visualization; however we can take a look right now.  The `catplot` function in `seaborn` allows you to split categorical plots such as point (the default), box, violin, bar, or strip plots into separate smaller charts, to facilitate comparisons across groups.  \n",
    "\n",
    "Let's create a new subsetted data frame, `sub1`, by indexing our `ug1` data frame for only those colleges and universities that are located in New York, Texas, and California.  We then call the `catplot` function, specifying how to divide our plots into __small multiples__ with the `col = 'state'` argument."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<seaborn.axisgrid.FacetGrid at 0x7f3415dc95e0>"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABDoAAAFqCAYAAAD7mzrpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA2MElEQVR4nO3da3hU5aH28XtmEg7Z1AyHBGxNo1UEi0R76RZCYui2VRSYSBAtggTEolK9sqWtDSqgggdCaT1gtxa3RVFTRQ6a1F2KgiQ7IeDhupRdKyeVEEBIIJmRADHJzHo/8DIaUZyEWVmz1vx/X8zKDM/cziRPJnee9SyXYRiGAAAAAAAAHMBtdQAAAAAAAIBooegAAAAAAACOQdEBAAAAAAAcg6IDAAAAAAA4BkUHAAAAAABwDIoOAAAAAADgGBQdiGtTp05VfX191O7XUV988YUeffRRjRkzRldffbV8Pp8WL16sr1/9+fnnn9eAAQP0/vvvm5YFAKItFubaTZs2aeDAgaqsrGzz+blz52rRokWqqanRRRddpLfeeqvN7WVlZRo2bJj2799vSi4AiDar59zPP/9cV199ta6++mpdfvnlysjICB8XFRXprbfe0pAhQ7Rnz57wv6murtall16qDz/8MOp5EJ8SrA4AWOnrb3hP9X4dYRiGfvWrX+mss87Syy+/rK5du6qhoUG33HKLjhw5ojvuuCN835deekk+n0/PPfecLrzwQtMyAUA0xcJcK0mJiYkqLCxUSUmJevXq1ea2tLQ0zZ49W7NmzVJpaal69eqlhoYGzZo1S0VFRerbt6+p2QAgWqyec0877TS99tprko6VzPPmzQsfH3fdddfpP//zP1VcXKzm5mb96le/0m9/+1sNGjTIlEyIPxQdcLzDhw/rrrvuUnV1tdxutwYNGqS5c+fqnnvukSRNnjxZixcv1pYtW/TnP/9Zzc3Nqq+v15gxY3THHXforrvuanM/t9utuXPn6rPPPlNLS4tGjRqlW2+99YTHLSgoUHV1dZvPnXHGGfrTn/7U5nPvvPOOPvnkEy1evFgej0eS1LNnTy1YsKBN071p0yYFAgHdeeeduvzyy/XZZ5/p9NNPj+pzBQAdFetzrSSlp6crIyNDd999t5566qkTbh8zZowqKip033336fHHH9e9996rvLw8XXrppdF4igAgauww557MHXfcoffff18LFy7Uvn37lJ2drauvvrqDzwbwDQzA4VatWmVMnTrVMAzDaG1tNe655x5j586dhmEYxrnnnmscPHjQCIVCxg033GB8+umnhmEYxr59+4zzzjvPOHjwYJv7GYZhTJo0yVi7dq1hGIbR1NRkTJo0yXj99dc7nO+ZZ54xCgoKvvN+BQUFxvz58w3DMIxp06YZCxYs6PBjAkC0xfpcu3HjRmPUqFHG4cOHjSuuuMJ4/vnnDcMwjPvvv994/PHHw/c7dOiQcdlllxkzZ840rr/+eqOlpaXDjwkAZon1Ofe443PvN6mtrTUuueQSY9y4ccy1iDpWdMDxLrroIj3yyCOaNGmShg0bpsmTJys9Pb3NfVwul5566imtX79ef/vb3/Txxx/LMAwdPXq0zf2OHDmid955R4FAQI899lj4c1u2bNHIkSPb3DfSxtvtdp+wF8fX1dXVae3atVqxYoWkY391vO+++3TbbbcpKSkp8icDAEwS63PtcUlJSfrjH/+o/Px8XXLJJSfc3qNHDz3wwAOaNm2a3nzzTSUk8FYJQOyxy5x7Mh9++KESExP16aefqrq6WmeffXa7xwC+DT+94XhpaWl64403tGnTJm3cuFE33nij5s6dq8suuyx8nyNHjigvL08///nPdfHFF+uaa67Rm2++eUIBEQqFZBiGXnrpJXXv3l2SVF9fr65du57wuI8//nhE+S644AI999xzCgaD4VNXJGnz5s16/vnn9fvf/17Lli2TJE2fPj2co7GxUatWrdLEiRPb94QAgAlifa79qkGDBmn69On6zW9+o4yMjG/8f0lMTFS/fv3aPTYAdAY7zbnfZOfOnZo5c6aeeOIJbdiwQQUFBXrllVf4Ax6ihquuwPGKi4t11113KTs7W3feeaeys7P1r3/9S5Lk8XjU2tqq6upqNTY26o477tBll12mTZs2qbm5WaFQqM39evTooQsvvFBLliyRdGxX6euvv15r167tcL6f/OQn+tGPfqSHH35YX3zxhSTpwIEDeuCBB3TGGWcoGAzqlVde0f33369169Zp3bp1Wr9+vW655RYtXbr0O1eDAEBniPW59utuuukm9enTRyUlJVEbEwA6i93m3K86fPiwbr/9dt188826+OKLdfvtt6tPnz6aM2eOKY+H+ETRAccbM2aMgsGgRo4cqbFjx+rQoUOaNGmSJOnKK6/UpEmT5HK59NOf/lRXXXWVrrrqKr311ls655xzwkvzjt9v27ZtWrhwoT744AP5fD5de+21Gj16tHJzc08p4/F2fOzYscrNzdWUKVN0xRVXqKCgQG+99ZZCoZB8Pl+bfzNlyhQdOHBAZWVlp/TYABANdphrv8rlcqmoqEj/9m//FrUxAaCz2G3O/aq77rpLZ555pqZOnSrp2Gncf/jDH7Rx40YVFxeb8piIPy6DPwcDAAAAAACHYEUHAAAAAABwDIoOAAAAAADgGBQdAAAAAADAMSg6AAAAAACAY1B0AAAAAAAAx0iwOsCpOHiwUaEQF40BgEilpHyv3f+GuRYA2qcjc63EfAsA7XGyuZYVHQAAAAAAwDEoOgAAAAAAgGNQdAAAAAAAAMeg6AAAAAAAAI5B0QEAAAAAAByDogMAAAAAADgGRQcAAAAAAHAMig4AAAAAAOAYFB0AAAAAAMAxTC06GhsbNXr0aO3evfuE2z766CONHTtWI0aM0D333KPW1lYzowAAAAAAgDhgWtHxwQcf6Prrr9fOnTu/8fY777xTc+bM0T/+8Q8ZhqFly5aZFQUAAAAAAMSJBLMGXrZsme6991797ne/O+G2PXv2qKmpSRdeeKEkaezYsXr88cc1YcIEs+IAcIjKynJVVJSZMnYg4JckJSd7oz52dvZwZWXlRH1cAAAAAG2ZVnQ8+OCD33pbbW2tUlJSwscpKSnav39/ux+jd+8eHcoGwL5OO627EhM9poz9+ecBSVKfPr2jPvZpp3VXSsr3oj5uZ2CuBYDOwXwLANFhWtFxMqFQSC6XK3xsGEab40gdPNioUMiIZjQAMW7w4H/X4MH/bsrYRUXzJEm//vXdpoxfV3fIlHHboyNlC3MtALRPR4tt5lsAiNzJ5lpLrrrSr18/1dXVhY8PHDig1NRUK6IAAAAAAAAHsaTo+MEPfqCuXbvqvffekyS99tprysnh3HUAAAAAAHBqOvXUlWnTpqmgoECDBw/WwoULNWvWLDU2NmrQoEHKz8/vzCiIc2ZtaGnmZpYSG1oCAAAAwHcxvehYt25d+OOnn346/PHAgQO1fPlysx8e6FSBwLHNLM0qOgAAAAAAJ2fJZqSA1bKyckxZGXF8M8vCwtlRHxsAAAAA8N0s2aMDAAAAAADADBQdAAAAAADAMSg6AAAAAACAY1B0AAAAAAAAx6DoAAAAAAAAjkHRAQAAAAAAHIOiAwAAAAAAOAZFBwAAAAAAcAyKDgAAYDt+f4Pmz5+rQMBvdRQAABBjKDoAAIDtlJau0vbtW1VSstLqKAAAIMZQdAAAAFvx+xtUUVEmwzBUUVHOqg4AANBGgtUBAAAA2qO0dJVCIUOSFAqFVFKyUpMmTbU4FQAA7VdZWa6KijJTxj7+h4DkZG/Ux87OHq6srJyojxstrOgAAAC2UlVVqWCwVZIUDLaqqqrS4kQAAMSeQCCgQCBgdQxLsKIDAADYSmZmlsrL1ysYbJXHk6DMzCyrIwEA0CFZWTmmrYwoKponSSosnG3K+LGMFR0AAMBWfL48ud0uSZLb7VZu7liLEwEAgFhC0QEAAGzF6+2p7Ozhcrlcys7OMeXcYwAAYF+cugIAAGzH58vTnj27Wc0BAABOQNEBAABsx+vtqZkz51gdAwAAxCBOXQEAxCW/v0Hz588NX3oNAAAAzkDRAQCIS6Wlq7R9+1aVlKy0OgoAAACiiKIDABB3/P4GVVSUyTAMVVSUs6oDAADAQSg6AABxp7R0lUIhQ5IUCoVY1QEAAOAgFB0AgLhTVVWpYLBVkhQMtqqqqtLiRAAAAIgWig4AQNzJzMySx3PswmMeT4IyM7MsTgQAAIBooegAAMQdny9PbrdLkuR2u5WbO9biRAAAAIgWig4AQNzxensqO3u4XC6XsrNzlJzstToSAAAAoiTB6gAAAFjB58vTnj27Wc0BAADgMBQdAIC45PX21MyZc6yOAQAAbKC4eKlqaqqtjtEuu3Ydy1tUNM/iJJFLS0vXhAn5pzwORQcAAAAA4AR+f4OeemqRpk8viPvTPGtqqrXtk63yJHexOkrEQp6gJOnjg59anCQywUBz1Mai6AAAAAAAnGD58pe0bdsWvfLKX/XLX063Oo7lPMldlJzzfatjOFagfG/UxmIzUgAAAABAG35/g6qqKiRJVVWVCgT81gYC2oGiAwAAAADQxvLlL8kwDEmSYYT0yit/tTgREDmKDgAAAABAG5s2bTjpMRDLKDoAAHHJ72/Q/PlzWYoLAADgMBQdAIC4VFq6Stu3b1VJyUqrowAAEHOGDBnW5njo0CyLkgDtR9EBAIg7fn+DKirKZBiGKirKWdUBAMDXjBs3Xm73sV8X3W63xo0bb3EiIHIUHQCAuFNaukqh0LEN1kKhEKs6AAD4Gq+3Z3gVR2ZmtpKTvdYGAtqBogMAEHeqqioVDLZKkoLBVlVVVVqcCACA2DNu3Hide+5AVnPAdig6AABxJzMzSx5PgiTJ40lQZibnHQMA8HVeb0/NnDmH1RywHYoOAEDc8fny5Ha7JB077zg3d6zFiQAAABAtFB0AgLjj9fZUdvZwuVwuZWfn8JcqAAAAB0mwOgAAAFbw+fK0Z89uVnMAAAA4DEUHACAuHT/vGAAAAM7CqSsAAAAAAMAxKDoAAAAAAIBjmFp0lJaWauTIkbriiiv04osvnnD7hx9+qGuuuUa5ubm65ZZb9Pnnn5sZBwAAAAAAOJxpRcf+/fv1yCOPqLi4WK+++qpefvll7dixo819HnzwQRUUFKikpERnnXWWnnnmGbPiAAAAAACAOGBa0bFhwwYNHTpUXq9XSUlJGjFihFavXt3mPqFQSIcPH5YkHT16VN26dTMrDgAAAAAAiAOmXXWltrZWKSkp4ePU1FRt3ry5zX1mzpypqVOn6qGHHlL37t21bNmydj1G7949opIViJbERI8kKSXlexYnQUfw+n0z5loA6BzMt0DsOv4+EeZKTPRE5b24aUVHKBSSy+UKHxuG0ea4qalJ99xzj5599lllZGRoyZIlKiws1OLFiyN+jIMHGxUKGVHNDZyKlpagJKmu7pDFSdAR8fD6deQHB3MtYpHf36Cnnlqk6dMLlJzstToO0EZH36Qz3wKx6/j7RJirpSUY8Xvxk821pp260q9fP9XV1YWP6+rqlJqaGj7etm2bunbtqoyMDEnSL37xC7399ttmxQEAAA5SWrpK27dvVUnJSqujAACAGGNa0TFs2DBVVVWpvr5eR48e1Zo1a5STkxO+PT09Xfv27dMnn3wiSVq7dq0GDx5sVhwAAOAQfn+DKirKZBiGKirKFQj4rY4EAABiiGmnrvTt21czZsxQfn6+WlpaNG7cOGVkZGjatGkqKCjQ4MGD9fDDD+uOO+6QYRjq3bu3HnroIbPiAAAAhygtXRVe3h8KhVRSslKTJk21OBUAAIgVphUdkuTz+eTz+dp87umnnw5/PHz4cA0fPtzMCAAAwGGqqioVDLZKkoLBVlVVVVJ0AACAMNNOXQEAADBDZmaWPJ5jf6vxeBKUmZllcSIAcCa/v0Hz58/lFEHYDkUHAACwFZ8vT273sSu5ud1u5eaOtTgRADgTGz/Drig6AACArXi9PZWdPVwul0vZ2TlcXhYATMDGz7Azig4AAGA7Pl+e+vcfwGoOADDJN238DNgFRQcAALAdr7enZs6cw2oOADDJN238DNgFRQcAALAdNsgDAHOx8TPszNTLywIAAJhh+fKXtG3bFi1f/pJuuulWq+MAgOP4fHmqqChTMMjGz5IUCPjV6v9CgfK9VkdxrFb/Fwok+KMyFis6AACArfj9Ddq48dgS6qqqClZ1AIAJ2PgZdsaKDgBAXPL7G/TUU4s0fXoBb95sZvnylxQKhSQd2yCPVR0AYA6fL0979uyO+9UckpSc7NWB1gYl53zf6iiOFSjfG7X3ZKzoAADEpdLSVdq+fSu7yNvQpk0b2hwfX90BAIguNn6GXVF0AADijt/foIqKMhmGoYqKck59AAAAcBCKDgBA3CktXaVQyJB07NQHVnXYy5Ahw056DAAA4htFBwAg7lRVVSoYbJUkBYOtqqri1Ac7GTduvFwulyTJ5XLr2muvtzgRAACIJRQdAIC4k5mZJY/n2H7cHk+CMjOzLE6E9vB6eyozM1vSsdeSc8cBAMBXUXQAAOKOz5cnt/vYigC3281u8jY0btx4nXvuQFZzAACAE1B0AADijtfbU9nZw+VyuZSdncOKABviSgAAAODbJFgdAAAAK/h8edqzZzerOQAAAByGogMAEJeOrwgAAACAs3DqCgAAAAAAcAyKDgAAAAAA4BicugLAFMXFS1VTU211jHbZtetY3qKieRYniVxaWromTMi3OgYAAAAQMyg6AJiipqZa2z7ZKk9yF6ujRCzkCUqSPj74qcVJIhMMNFsdAQAAAIg5FB0ATONJ7qLknO9bHcOxAuV7rY4AAAAAxBz26AAAAAAAAI5B0QEAAAAAAByDogMAAAAAADgGe3QAAADgBJWV5aqoKIv6uIGAX5KUnOyN+tjZ2cOVlZUT9XEBAPbCig4AAAB0mkAgoEAgYHUMAICDsaIDAAAAJ8jKyjFldURR0TxJUmHh7KiPDQCAxIoOAAAAAADgIBQdAAAAAADAMTh1BQAAAIhzZm0+K7EBLYDOx4oOAAAAAKZhA1oAnY0VHQAAAECcM2vzWYkNaAF0PooOAAAAALAxs0494rQj2BWnrgAAAAAATsBpR7ArVnQAAAAAgI2ZdeoRpx3BrljRAQAAAAAAHIOiAwAAAAAAOAZFBwAAAAAAcAyKDgAAAAAA4BgUHQAAAAAAwDEoOgAAAAAAgGNQdAAAANvx+xs0f/5cBQJ+q6MAAIAYQ9EBAABsp7R0lbZv36qSkpVWRwEAADHG1KKjtLRUI0eO1BVXXKEXX3zxhNs/+eQTTZo0Sbm5ubrpppsUCATMjAMAABzA729QRUWZDMNQRUU5qzoAAEAbphUd+/fv1yOPPKLi4mK9+uqrevnll7Vjx47w7YZhaPr06Zo2bZpKSkp03nnnafHixWbFAQAADlFaukqhkCFJCoVCrOoAAABtmFZ0bNiwQUOHDpXX61VSUpJGjBih1atXh2//8MMPlZSUpJycHEnSrbfeqokTJ5oVBwAAOERVVaWCwVZJUjDYqqqqSosTAQCAWGJa0VFbW6uUlJTwcWpqqvbv3x8+3rVrl/r06aO7775beXl5uvfee5WUlGRWHAAA4BCZmVnyeBIkSR5PgjIzsyxOBAAAYkmCWQOHQiG5XK7wsWEYbY5bW1v19ttv64UXXtDgwYP16KOPav78+Zo/f37Ej9G7d4+oZgZOVWKiR5KUkvI9i5NY7/hzAXMlJnpM/3pjrkWsmTJlkioryxUMSh6PWzfemK+ePZl37YKfld/OqfMtr7l98dp9ife2nSNa721NKzr69eund999N3xcV1en1NTU8HFKSorS09M1ePBgSdLo0aNVUFDQrsc4eLAxfI4uEAtaWoKSpLq6QxYnsd7x5wLmamkJtuvrrSM/OJhrEXsSlZWVo/Xr1yorK0etrQnMuzYSDz8rO/om3anzbTy85k7Fa/cl3tt2jva8tz3ZXGta0TFs2DAtWrRI9fX16t69u9asWaN58+aFb//JT36i+vp6bdmyRQMHDtS6des0aNAgs+IAAAAH8fnytGfPbuXmjrU6CgAgTgQDzQqU77U6RsRCTcfKGXc3e6xGCQaapd7RGcu0oqNv376aMWOG8vPz1dLSonHjxikjI0PTpk1TQUGBBg8erD/96U+aNWuWjh49qn79+mnBggVmxQEAAA7i9fbUzJlzrI4BAIgTaWnpVkdot127qiVJP+xtk+y9o/c8m1Z0SJLP55PP52vzuaeffjr88QUXXKDly5ebGQEAAAAAgFMyYUK+1RHarajo2BkVhYWzLU7S+Uy76goAAAAAAEBni6jo+MMf/qCamhqzswAAAAAAAJySiIoOwzA0ceJETZ48WX/729/U3Nxsdi4AAAAAAIB2i6jo+O1vf6v169dr6tSpeuONN/Tzn/9cDz30kNnZAAAAAAAA2iXizUjdbrfOP/98ffrpp9q5c6feffddM3MBAAAAAAC0W0RFxxtvvKEVK1bo/fff15VXXqmHHnpIgwYNMjsbAAAAAABAu0RUdDzzzDO67rrr9Oijj6pbt25mZwIAAAAAAOiQiPboeOmll/Qf//Efqqqq0v/+7//q888/NzsXAACm8vsbNH/+XAUCfqujAAAAIIoiKjrKysp05ZVX6r//+7/1X//1Xxo5cqTeeecds7MBAGCa0tJV2r59q0pKVlodBQAAAFEU0akrjz76qF544QX1799fkvThhx9q9uzZWrmSN4cAAPvx+xtUUVEmwzBUUVGu3NyxSk72Wh0LAAAAURDRig6XyxUuOSRp0KBBMgzDtFAAAJiptHSVQqFjP8dCoRCrOgAAABzkpEWH3++X3+/X+eefr2eeeUaHDx/W0aNH9eKLL2ro0KGdlREAgKiqqqpUMNgqSQoGW1VVVWlxIgAAAETLSU9dGTp0qFwuV3j1xu9///vwbS6XS4WFheamAwDABJmZWSovX69gsFUeT4IyM7OsjgQAAIAoOWnRsWXLls7KAQBAp/H58lRRUaZgUHK73crNHWt1JAAAAERJRJuRNjc3q6ysTIcPH5YkBYNB7dq1SzNmzDA1HAAAZvB6eyo7e7jWr1+r7OwcNiIFAABwkIiKjhkzZqimpkZ1dXX68Y9/rA8++ECXXHKJ2dkAADCNz5enPXt2s5oDAADAYSIqOj766COtWbNG9913n2688UaFQiHdd999JkcDYGeBgF+t/i8UKN9rdRTHavV/oUCC3+oYtuX19tTMmXOsjgEAAIAoi+jysqmpqUpISNCZZ56pbdu2qX///jp06JDZ2QAAAAAAANolohUdSUlJKi0t1cCBA7Vs2TL96Ec/0pEjR8zOBsDGkpO9OtDaoOSc71sdxbEC5XvZWwIAAAD4mohWdMyZM0cfffSRsrKy5Ha7dcMNN+imm24yOxsAAAAAAEC7RLSi48wzz9Tvfvc77dmzR3fddZf69u1rdi4AAAAAAIB2O+mKjl//+teSpM8++0x5eXnKzc3VyJEjNXHiRO3bt69TAgIAAAAAAETqpEXHp59+Kkl6+OGHNXr0aL333nt69913NXLkSN1zzz2dEhAAAAAAACBSEe3R8emnn4b35HC5XKzoAAAAAAAAMemkRYff79f777+v9PR0ffzxx+HP79mzR253RB0JAABA1Pn9DZo/f64CAb/VUQAAQIw5aVsxYsQIzZ8/X+Xl5VqwYIEk6bXXXtPYsWN18803d0pAAACArystXaXt27eqpGSl1VEAAECMOelVV2bOnClJam5uVm1trSRp4MCBeuGFF9S/f3/z0wEAAHyN39+giooyGYahiopy5eaOVXKy1+pYAAAgRnzn+SebNm3S66+/rm7dukmSBgwYoP79+2vVqlWmhwMAAPi60tJVCoUMSVIoFGJVBwAAaOOkRceSJUs0Z84c/f3vf9fIkSO1cePG8G1Lly41PRwAAMDXVVVVKhhslSQFg62qqqq0OBEAAIglJy06VqxYoRUrVmjx4sVauHChZsyYoa1bt0qSDMPolIAAAABflZmZJY/n2Nm3Hk+CMjOzLE4EAABiyUmLjoSEBPXo0UOSlJOTo8LCQk2fPl0NDQ1yuVydEhAAAOCrfL688PsQt9ul3NyxFicCAACx5KSbkfbq1UsrVqzQ6NGj1bVrV40ZM0ZbtmzRzTffrEOHDnVWRsSx4uKlqqmptjpGxHbtOpa1qGiexUnaJy0tXRMm5FsdAwAi4vX2VGpqqvbu3aOUlFQ2IgUAAG2ctOi47777NHPmTLndbuXl5Uk6diWWoqIiPffcc50SEPGtpqZaO3dsUb8eJ/1SjRlJCkmSmvbtsDhJ5PY1tlodAQDaxe9vUG3tfklSbW2tAgE/ZQcAAAg76W+PP/zhD1VcXHzC5wsLCzVlyhSzMgFt9OuRoBszelkdw7GWbK63OgIAtEtp6Sod3yrMMAyVlKzUpElTrQ0FAABiRkR/Jr/rrrvaHLtcLnXv3l39+/fXtddeK4/HY0o4AAAqK8tVUVEW9XEDAb8kmbISIDt7uLKycqI+Lo75pquuUHQAAIDjTroZ6Vf961//0oABA3Teeedpx44d2rt3ryoqKvTQQw+ZmQ8AAFMEAgEFAgGrY6ADuOoKAAA4mYhWdHz88cd68cUXw1dgufbaazV16lQVFxdr9OjRpgYEAMS3rKwcU1ZHHN80uLBwdtTHhrl8vjxVVJQpGJTcbjdXXQEAAG1EtKIjEAiESw5J6tatmxobG+VyuZSYmGhaOAAAgK/zensqO3u4XC6XsrNz2IgUAAC0EdGKjgsuuEC//e1vNW7cOBmGoRUrVigjI0NlZWXq3r272RkBAADaGD78Mm3cWKmf/vRnVkcBAAAxJqIVHXPnztX3v/99zZ8/XwsXLlR6erruvfdeHT58WHPnzjU7IwAAQBtlZevU1NSk9evXWh0FAADEmIhWdFx77bXhj4PBoNasWaM1a9aotLTUtGAAAADfxO9vUEVFmQzDUEVFuXJzx3L6CgAACIuo6Jg9+8uN2lpaWvT6668rLS3NtFAAAADfprR0lUIhQ5IUCoVUUrKSy8sCAICwiIqOSy65pM3xsGHDNH78eE2fPt2UUAAAAN+mqqpSwWCrJCkYbFVVVSVFBwAACIuo6Pi6hoYG1dbWRjsLAADAd8rMzFJ5+XoFg63yeBKUmZlldSRLFRcvVU1NtdUxIrZr17Gsxy/xbBdpaemaMCHf6hgAgAhEVHT4fL42x3v37tUvfvELUwIBAACcjM+Xp4qKMgWDktvtVm7uWKsjWaqmplrbPtkqT3IXq6NEJOQJSpI+PvipxUkiFww0Wx0BANAO7d6jw+VyqVevXjr77LNNCwUAAPBtvN6eys4ervXr1yo7O4eNSCV5krsoOef7VsdwrED5XqsjhNltBY9kz1U8rOAB7K1De3REqrS0VE8++aRaW1s1efJkTZw48Rvvt379es2dO1fr1q3r0OMAAID44vPlac+e3XG/mgPxp6amWlu375Cnm9fqKBELBT2SpB01ByxOEplgk9/qCABOUYf26IjE/v379cgjj2jlypXq0qWLxo8fryFDhuicc85pc78DBw6oqKjIrBgAAMCBvN6emjlzjtUxAEt4unmVlP4zq2M41pHqtVZHAHCK3GYNvGHDBg0dOlRer1dJSUkaMWKEVq9efcL9Zs2apdtvv92sGAAAAAAAII6YVnTU1tYqJSUlfJyamqr9+/e3uc/SpUv14x//WBdccIFZMQAAAAAAQBwx7dSVUCgkl8sVPjYMo83xtm3btGbNGj377LPat29fhx6jd+8ep5wTsS0x0aMmq0PEgcREj1JSvhf1MWE+M167r3PqXHv8a9Ts5w8wG/Nt54iV+ZbXu3N0xuttB/ystLd4fv1MKzr69eund999N3xcV1en1NTU8PHq1atVV1ena665Ri0tLaqtrdWECRNUXFwc8WMcPNioUMiIam7ElpaWoNUR4kJLS1B1dYeiPmYw0BxTO9V/l1DTsa83dzd7vIkMBprV0rt9r11HftA5da49Pr9E+2sf6Gz8rOwc7flZ2dFfKiKZb3m9O4cZ743siJ+V9ub01+9kc61pRcewYcO0aNEi1dfXq3v37lqzZo3mzfvyklIFBQUqKCiQJO3evVv5+fntKjkAxLa0tHSrI7Tb8cvf/bC3TbL3tufzDAAAAJjJtKKjb9++mjFjhvLz89XS0qJx48YpIyND06ZNU0FBgQYPHmzWQwOIAXa89nxR0bEytrBwtsVJAAAAAHSUaUWHJPl8Pvl8vjafe/rpp0+43xlnnKF169aZGQUAAAAAAMQBU4sOp6usLFdFRVnUxw0E/JKk5GRv1MfOzh6urKycqI8LAAAAAEAsMO3ysui4QCCgQCBgdQwAAAAAAGyHFR2nICsrx5TVEewTAAAAAABAx7CiAwAAAAAAOAZFBwAAAAAAcAyKDgAAAAAA4BgUHQAAAAAAwDEoOgAAAAAAgGNQdAAAAAAAAMfg8rIAAMAUlZXlqqgoM2XsQMAvSUpO9kZ97Ozs4aZcPh4AAHQOVnQAAADbCQQCCgQCVscAAAAxiBUdAADAFFlZOaatjCgqmidJKiycbcr4AADAvljRAQAAAAAAHIMVHQAAAABgsuLipaqpqbY6Rrvs2nUs7/FVdHaRlpauCRPyrY4RETP3szLz9Yv1/awoOgAAAADAZDU11dq5Y4v69bDPr2BJCkmSmvbtsDhJ5PY1tlodIWYkJydbHcEy9vkuAwAAAAAb69cjQTdm9LI6hqMt2VxvdYR2MXM/q3jGHh0AAAAAAMAxKDoAAAAAAIBjUHQAAAAAAADHoOgAAAAAAACOQdEBAAAAAAAcg6IDAAAAAAA4BkUHAAAAAABwDIoOAAAAAADgGBQdAAAAAADAMSg6AAAAAACAYyRYHQAAAABAZAIBv4JNfh2pXmt1FMcKNvkVCPBrEmBnrOgAAAAAAACOQVUJAAAA2ERysld1n7cqKf1nVkdxrCPVa5Wc7LU6BoBTwIoOAAAAAADgGBQdAAAAAADAMSg6AAAAAACAY7BHBwAgKoqLl6qmptrqGBHbtetY1qKieRYnaZ+0tHRNmJBvdQwAAICYRdEBAIiKmppqbd2+Q55uXqujRCQU9EiSdtQcsDhJ5IJNfqsjAAAAxDyKDgBA1Hi6ebkSgImOVK+1OgIAAEDMY48OAAAAAADgGBQdAAAAAADAMSg6AAAAAACAY1B0AAAAAAAAx6DoAAAAAAAAjkHRAQAAAAAAHIOiAwAAAAAAOAZFBwAAAAAAcAyKDgAAAAAA4BgJVgcAAAAAAKcLBPxqaGzVks31VkdxtH2NreoZ8FsdAxZjRQcAAAAAAHAMU1d0lJaW6sknn1Rra6smT56siRMntrn9zTff1KJFi2QYhs444ww9/PDDSk5ONjMSbIbm23y03gAAAOZLTvaq69EDujGjl9VRHG3J5np1S/ZaHQMWM63o2L9/vx555BGtXLlSXbp00fjx4zVkyBCdc845kqTGxkbdd999WrFihfr27avHHntMixYt0qxZs6Kao7h4qWpqqqM6ptl27TqWt6honsVJIpeWlq4JE/KtjoE4UFlZroqKMlPGNvN7Lzt7uLKycqI+LgAAAIC2TCs6NmzYoKFDh8rr9UqSRowYodWrV+v222+XJLW0tOjee+9V3759JUkDBgxQaWlp1HPU1FRr6/Yd8nTzRn1ss4SCHknSjpoDFieJTLDJb9rYNN/mo/X+EivKAAAAAPszreiora1VSkpK+Dg1NVWbN28OH/fs2VOXX365JKmpqUmLFy/WpEmT2vUYvXv3+M77JCZ65OnmVVL6z9o1NiJ3pHqtEhM9Skn5XtTHTkz0qCnqo+LrzHr9zDBmzCiNGTPK6hhxJZK5Vjr2dQTz2en71UzHv954Lvje6yyd8b0X6XtbmM+M15v3tZ2Hn5UwregIhUJyuVzhY8Mw2hwfd+jQId12220aOHCg8vLy2vUYBw82KhQyTnqflpZgu8ZEx7S0BFVXd8iUcWE+s14/xJ6O/NCPZK6V+H7tLHy/HnP8643ngu+9ztKe772O/oLFe9vYYcZcy2vXefhZGR9ONteadtWVfv36qa6uLnxcV1en1NTUNvepra3VhAkTNGDAAD344INmRQEAAAAAAHHCtBUdw4YN06JFi1RfX6/u3btrzZo1mjfvyw3+gsGgbr31Vl111VX61a9+ZVYMAAAAwFGCTX4dqV5rdYyIhVqPnbDhTuhmcZLIHNt/ro/VMQCcAtOKjr59+2rGjBnKz89XS0uLxo0bp4yMDE2bNk0FBQXat2+f/vWvfykYDOof//iHJOn8889nZQcAAADwLdLS0q2O0G7Hr2r2wzS7lAd9bPk8A/iSaUWHJPl8Pvl8vjafe/rppyVJgwcP1pYtW8x8eAAAAMBRJkzItzpCux2/bHth4WyLkwCIF6YWHbEgEPDbbnmf3QSb/AoEHP+lBACOVVy8VDU11VbHaJfjfyE+/guUHaSlpdvyl1QAAOyG304BAIhzNTXV2rlji/r1sM/bgiSFJElN+3ZYnCQy+xpbrY4AAEDcsM87mg5KTvaq7vNWJaX/zOoojnWkeq2Sk71WxwBgMVbQmc/MFXT9eiToxoxepowNacnmeqsjAAAQN0y7vCwAAAAAAEBnc/yKDgBA52AFnflYQQcAAPDdWNEBAAAAAAAcg6IDAAAAAAA4BkUHAAAAAABwDIoOAAAAAADgGBQdAAAAAADAMSg6AAAAAACAY1B0AAAAAAAAx6DoAAAAAAAAjkHRAQAAAAAAHIOiAwAAAAAAOAZFBwAAAAAAcAyKDgAAAAAA4BgJVgcAAABAxwUCfrX6v1CgfK/VURyr1f+FAgl+q2MAACLEig4AAAAAAOAYrOgAAACwseRkrw60Nig55/tWR3GsQPleJSd7rY4BAIgQKzoAAAAAAIBjUHQAAAAAAADH4NQVAAAAmwsGmm2zGWmoKShJcnfzWJwkcsFAs9Tb6hQAgEhRdAAAANhYWlq61RHaZdeuaknSD3vbKHdv+z3PABDPKDoAAABsbMKEfKsjtEtR0TxJUmHhbIuTAACciqIDAAAAADrBvsZWLdlcb3WMiDU2hyRJPbrYZ2vHfY2tOtPqELAcRQcAAAAAmMyOpz/V/v9Tzfr0s0/2M2XP5xrRRdGBmGen5pvWGwAAAN/EbqeZSZxqBvui6EBMs1sbS+uNeBds8utI9VqrY0Qk1NokSXIndLM4SeSCTX5JfayOAQAAENMoOhDT7NZ803ojntmtMAtf+SHNTsVBH1Oe50DArwYbrZ6zo32NreoZ8FsdAwCAuEDRAQCICopJAAAAxIK4KDrstJRast9yapZSA4C9JSd71fXoAd2Y0cvqKI61ZHO9uiV7rY4BAEBccHzRYbel1JIdl1Obs5QaAAAAAID2cnzRYbel1BLLqQEAAAAA6Cj7XAMTAAAAAADgO1B0AAAAAAAAx3D8qSsAAABov8rKclVUlEV93ON7kR0/VTeasrOHKysrJ+rjAgDshaIDAAAAnSY5OdnqCAAAh6PoAAAAwAmysnJYHQEAsCX26AAAAAAAAI7Big4AAKB9ja1asrne6hgRa2wOSZJ6dLHH32z2NbbqTKtDACdh1p4sEvuyAOh8FB0AAMS5tLR0qyO0W+3//8WpTz97ZD9T9nyegWhgXxYAnY2i4xSwGzkAwAkmTMi3OkK7Hf8ZWVg42+IkgDOwJwsAJ6HoiEG03gAAAAAAdAxFxymg+QYAAAAAILbYYwcvAAAAAACACJi6oqO0tFRPPvmkWltbNXnyZE2cOLHN7R999JHuueceHT58WBdffLHuv/9+JSSwyAQA8CX2QwIAAEB7mNYq7N+/X4888ohWrlypLl26aPz48RoyZIjOOeec8H3uvPNOPfDAA7rwwgt19913a9myZZowYYJZkYAwO/7iJPHLExBN7IcEAHAKO7635X0tzGRa0bFhwwYNHTpUXq9XkjRixAitXr1at99+uyRpz549ampq0oUXXihJGjt2rB5//HGKDtgavzgB0cd+SPZl1htviTffANAZeG8LuzKt6KitrVVKSkr4ODU1VZs3b/7W21NSUrR///52PUbv3j1OPSji0pgxozRmzCirYwC2wFyLjjrttO5KTPSYMnbv3r0kyZTxTzutu1JSvhf1cYHvwnyLjuK9LdCWaUVHKBSSy+UKHxuG0eb4u26PxMGDjQqFjFMPCwBxoiO/vDHXoqMGD/53DR7871bH6JC6ukNWR4CNdbQoY74FgMidbK417aor/fr1U11dXfi4rq5Oqamp33r7gQMH2twOAAAAAADQXqYVHcOGDVNVVZXq6+t19OhRrVmzRjk5X57v+oMf/EBdu3bVe++9J0l67bXX2twOAAAAAADQXqYVHX379tWMGTOUn5+vMWPGaPTo0crIyNC0adP0f//3f5KkhQsX6uGHH9aVV16pI0eOKD8/36w4AAAAAAAgDrgMw7DtiYCcxwgA7cMeHQBgPvboAADzWbJHBwAAAAAAQGej6AAAAAAAAI5B0QEAAAAAAByDogMAAAAAADgGRQcAAAAAAHAMig4AAAAAAOAYCVYHOBVut8vqCADgeMy1ANA5mG8BIDpchmFwsW4AAAAAAOAInLoCAAAAAAAcg6IDAAAAAAA4BkUHAAAAAABwDIoOAAAAAADgGBQdAAAAAADAMSg6AAAAAACAY1B0AAAAAAAAx6DoAAAAAAAAjkHRAQAAAAAAHIOiIwY1NjZq9OjR2r17t9VR0E5PPPGERo0apVGjRmnBggVWx0E7PPbYYxo5cqRGjRqlJUuWWB0HnYC51r6Ya+2N+Tb+MN/aE3OtvcX7XEvREWM++OADXX/99dq5c6fVUdBOGzZsUEVFhVatWqVXX31VH374od544w2rYyECb7/9tjZu3KiSkhKtWLFCzz//vD755BOrY8FEzLX2xVxrb8y38Yf51p6Ya+2NuZaiI+YsW7ZM9957r1JTU62OgnZKSUnRzJkz1aVLFyUmJurss8/W3r17rY6FCFxyySVaunSpEhISdPDgQQWDQSUlJVkdCyZirrUv5lp7Y76NP8y39sRca2/MtVKC1QHQ1oMPPmh1BHRQ//79wx/v3LlTf//73/XXv/7VwkRoj8TERD3++OP6y1/+oiuvvFJ9+/a1OhJMxFxrX8y19sd8G1+Yb+2Judb+4n2uZUUHEGXbt2/X1KlT9bvf/U5nnnmm1XHQDgUFBaqqqtJnn32mZcuWWR0HwEkw19ob8y1gD8y19hbPcy1FBxBF7733nqZMmaLf/OY3ysvLszoOIvTxxx/ro48+kiR1795dV1xxhbZu3WpxKgDfhrnWvphvAftgrrUv5lqKDiBqPvvsM912221auHChRo0aZXUctMPu3bs1a9YsNTc3q7m5WWvXrtVFF11kdSwA34C51t6YbwF7YK61N+Za9ugAouaZZ57RF198ofnz54c/N378eF1//fUWpkIkhg8frs2bN2vMmDHyeDy64oor+KEOxCjmWntjvgXsgbnW3phrJZdhGIbVIQAAAAAAAKKBU1cAAAAAAIBjUHQAAAAAAADHoOgAAAAAAACOQdEBAAAAAAAcg6IDAAAAAAA4BkUHYJInnnhCb775Zrv/3YABA1RfX29CIgBwHuZaAOgczLewE4oOwCSbNm1Sa2ur1TEAwNGYawGgczDfwk4SrA4AxJrly5dryZIlcrvd6tmzp4qKilReXq7nn39ebrdbffr00ezZs3XWWWdp5syZ6tGjh7Zu3ap9+/ZpwIABKioq0quvvqp//vOfWrBggTwej9auXSu/36+amhr99Kc/1a233qr7779fW7Zskcvl0qWXXqpf//rXSkjgWxJAfGCuBYDOwXyLeMRXHvAVW7Zs0cKFC7Vq1SqdfvrpevbZZzVlyhSFQiG9/PLL6tWrl1auXKnbbrtNr7/+uiTpn//8p5YuXSqXy6XrrrtOq1ev1sSJE8P/vfzyy7V27Vo1NTWF/01hYaG8Xq9KS0vV0tKi6dOn6y9/+YtuvvlmK//3AaBTMNcCQOdgvkW84tQV4CuqqqqUnZ2t008/XZI0ZcoU/exnP9PIkSPVq1cvSdLYsWO1f/9+7d69W5J06aWXqkuXLkpMTNS5556rQCDwjWNfdNFF4Y/Ly8t1ww03yOVyqUuXLho/frzKy8tN/r8DgNjAXAsAnYP5FvGKogP4Co/HI5fLFT5uampSTU3NCfczDCN8jmK3bt3Cn3e5XDIM4xvHTkpKCn8cCoXaPE4oFOKcRwBxg7kWADoH8y3iFUUH8BVDhgxRVVWVamtrJUkvvfSSysrK9D//8z/h3aJXrFghr9er9PT0k47l8Xi+dYLPzs7WCy+8IMMw1NzcrGXLlmnYsGHR/Z8BgBjFXAsAnYP5FvGKPTqArxgwYIDuvPNO/fKXv5QkpaSk6I033tCbb76pyZMnKxQKqVevXvrzn/8st/vkPeFll12mP/7xj2ppaTnhtlmzZumBBx6Qz+dTS0uLLr30Ut16662m/D8BQKxhrgWAzsF8i3jlMr5tLRIAAAAAAIDNcOoKAAAAAABwDIoOAAAAAADgGBQdAAAAAADAMSg6AAAAAACAY1B0AAAAAAAAx6DoAAAAAAAAjkHRAQAAAAAAHIOiAwAAAAAAOMb/A50CQhMX15g8AAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 1080x360 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sub1 = ug1.query(\"state in ['NY', 'TX', 'CA']\", engine = 'python')\n",
    "\n",
    "sns.catplot(data = sub1, x = 'control', y = 'ug25abv', \n",
    "               col = 'state', kind = 'box', order = [1, 2, 3])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "Now, we want to look at how the percentage of undergraduates over age 25 varies by institution size.  Certainly, we could create a scatter plot as both the `ug25abv` and `ugds` columns are numeric:  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<seaborn.axisgrid.FacetGrid at 0x7f3415cd5d60>"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWQAAAFqCAYAAAAtESBQAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABaHklEQVR4nO3deXxU5dn4/89ZZslksm/si0RFwSCKhapFrSIg4EK19SvVtlprbX9PrX2KBW216tf6lLq2VZ+2+tTnW8GlKlAsota2tFWr4MIiCBJ2AtnXmcx2zv374yRjdhLIMpNc79erxZk5M7kygSv33Oc616UppRRCCCEGnD7QAQghhHBIQhZCiAQhCVkIIRKEJGQhhEgQkpCFECJBSEIWQogEYQ50AP2tsrIB2078Sr+sLB/V1cGBDqPbkineZIoVkiveZIoVBibevLy0Th+TFXKCMk1joEPokWSKN5liheSKN5lihcSLVxKyEEIkCEnIQgiRICQhCyFEgpCELIQQCUISshBCJAhJyEIIkSAkIQshRIKQhCyEEAliwBJyQ0MD8+fP5+DBg+0e2759OwsXLmT27NnccccdxGIxAEpKSli0aBFz5szh5ptvJhAI9HfYQgjRZwYkIW/atIn/83/+D3v37u3w8cWLF3PnnXfy2muvoZTihRdeAODuu+/mmmuuYd26dUyePJnHH3+8x1/7tife5s4n/3084QshRJ8YkIT8wgsvcNddd5Gfn9/usUOHDhEKhTj99NMBWLhwIevWrSMajbJhwwZmz57d6v5jcbAiKElZCJFwBqS50H333dfpY2VlZeTl5cVv5+XlUVpaSnV1NX6/H9M0W91/rA5WJE8DFCHE0JBw3d5s20bTtPhtpRSapsX/bKnt7Z7qqutSIkj0+NpKpniTKVZIrniTKVZIrHgTLiEPGzaM8vLy+O2Kigry8/PJzs6mvr4ey7IwDIPy8vIOtzx6ory8/njD7TN5eWkJHV9byRRvMsUKyRVvMsUKAxNvUrXfHDlyJB6Ph/fffx+A1atXM3PmTFwuF9OmTWPt2rUArFq1ipkzZx7X19pcXNHlY8tWfMBtT7zNshUfdHmsEEL0hoRJyDfeeCNbtmwB4IEHHuD+++9nzpw5BINBrrvuOgDuuusuXnjhBS655BI2btzI97///eP6mo++uJk7n/x3u2S7ubiC5W/spCYQwec1qQlEWP7GTknKQog+pSmlEn98Ri+64f++Tll1Y/y229TJ8LtZNOskiibkArBsxQfUBCJ4XJ81rw5HLTJT3dx2zRn9Eqd89Os7yRQrJFe8yRQryJZFwrFshWHorHt3f/y+itoQbrP1W+M2dSpqQ/0dnhBiCEm4k3r9zTT1dsk2N8PbboUcidnkZngHIsResbm4gnXv7qeiNkRuhpc508fEPxEIIRLDkF8hp/tc7ZLtnOljsCybcNRCKUU4amFZNnOmjxnASI+d7IkLkRyG/Ao5ZimgdbJtXjn2ZEWZyCvQde/uxzD0+Irf4zIIN92fKDEKISQhEwxFufnyye0SU9GE3G4nq+YVqGHorVagza8zEFr+gqhtiJDhd7XagpE9cSESz5DfsghHbZ55fedx1Rq3XIFqmobHZbQ7Udif2m5RaBpU10cIhqLxY5J9T1yIwWjIr5ABNI3jWtVW1IbweVu/lQO5Am27RZGZ5qGytpHahggpHpNIzO61PfHmlXhVQ4RsvzuhtmqESDZDfoUMUB+MYlk2tQ0RHnt5S49Xy7kZXiIxu9V9A7kCbVu2l+IxyU73opQiGIqRmdq67vpYtVyJp6XIyUIhjpeskIFQxCIUseK3S6sbW62Wm1eBhyoCWJbCNDRG5KbGV4Nzpo9h+Rs7CeOsjHtzBXosOirbMwydCSMzevXClo62auRkoRDHThJyGxpQH4yQ5nPH94CXv7GTaMwmGI6BgnBUtUva0HFVRsuTa163AUoRitp9WonRX78gEm2rRohkJwm5DcPQAI1gKEpJhc1v/7SNcDSGUs5es6Hr2EqjMRwjM80TXw12VJXRsvpC0+BwZQDQyEpz92klxrGU7R2LwXgBjRADSRJyC7quoTf1WI7EbGzlXFqNAlsBCizbRgNspR11NdjyI31pfRhN0wFFfWOMYdnuPv1435OyvWPVciVuGlrSX0AjxECThNyCbSuU9lnyTfWaNDRGaXO+DgUoW1EfjFKQldLp67X8SB+zbHRNQymINb1gsn+8b7kSr26IkCVVFkIcF0nIbVgKDF3DZWqkp7oJhGI4Kbi9hmCUr805udPXavmR3jT0pqsCnf4ZMDg+3jevxJOty5cQiUjK3lowdI2TRmXw3SsmM25YurNtYXfendTrMbpcDbbsiZHmc6GUjVKQlmLKx3shRDuyQm7B0LVWH7mXv7ETQ9ecrYwWx2mAbmiMzE3t8vXanlwbnpMar7LITJWP90IMJS0vonrqjlkdHiMJuQXLtttVPrz4t13tJlQrnP3miWMyj/qax3pybeP2Up5//ZOEbFYkhOiZlhVXaSmdp11JyE1MQ8c0tHgPipalbHc++W8OVwWxmk7uuUydVK/JJ/truLQPYtlcXMFzb+4CjT5tVpTIHeqEGEzaXkTVGUnITWyliFmKI5VBjlQGufPJfzNtYj4bPymLr5Bdpk6m343P60Ip1WcVEuve3Y9pahi6s8XfF1fAddSh7n/+vJ30VDehiCUJWohe1NFFVB2RhNxMqXh5mwYcqghy8F970Vv8MovGbCrrnCRsGHqfVUhU1IbI8LvjVRnQ+yVybRsQWZZNIBQjFLEYluNLiBaiQgwWHV1E1RGpsmjSspiiqQw5fr9pfJaVbRtqGyJ9WiGRm+ElHLVa3dfbJXJtGxDVBaNomvNJIRFaiAoxmLSdQtQZScjdoGsapqHRvPWjlOqVbmmdmTN9DLGY6tMRUm071MViTkmeaXz2VyLZL1wRIlEUTchl0ayTyEx109AY6/Q42bLoJl3TQNfQNThhRHqffowvmpBLRoavT6ss2jYgMnQNy7ZJT3XHjxkMF64IkShaXkTVGUnIR6E3fYx3PmUovF53v1zMMe2UAsbm+vrs9dvWSOdneqlrjKHrGkqpAW8hKsRQJAm5CznpHlLcBqXVITRNUZCdypXnTxg0J7na1khLGZwQA0sSchvOzoSzSrx29snxhNScrJ55fSe5Gfs7TVZHS2qJnPT6o0OcEKJzkpBb0ABN0zB0jfwsX6tk3J2p0kc7LhGnUwshEseAJOQ1a9bwxBNPEIvF+NrXvsaiRYvij23fvp0lS5bEb1dVVZGRkcErr7zCypUrefDBB8nJyQHg/PPP59Zbb+21uHIzvRiGjmXZXHn+hPj9bWt2O7tQ42jHdfd1hBBDU78n5NLSUh5++GFefvll3G43V199NdOnT6ewsBCAU045hdWrVwPQ2NjIVVddxU9/+lMAtm7dypIlS5g/f36vx6VroBQdNv3p7qiiox0nI4+EEF3p9zrkt99+mxkzZpCZmYnP52P27NmsW7euw2N/85vfcNZZZzFt2jQAtmzZwsqVK1mwYAE//OEPqa2t7bW4NE3jqxefxJzpY1j37n5ue+Lt+PTp7k6VPtpxiTadWgiRWPp9hVxWVkZeXl78dn5+Pps3b253XH19PS+88AJr1qyJ35eXl8f111/PGWecwUMPPcQ999zDgw8+2CtxaShW/WsvoXAM09TI8LtpCEV57s1dfPGs0fx1wwEs23a2GaIWKPjKxRNb1RR+5eKJ/OblzZ0ed7TH2+qqXjERJVO8yRQrJFe8yRQrJFa8/Z6Qbdtu1e1INV2q29af/vQnLrroovh+McBjjz0W/+9vfvObzJrVcU/RntA052SeQuNQWT05mSkYujPdw9B1YprFB9tLufrCws+mR7ucznCP/fGjVpUSY3N9rY5rfmxsro/y8vqjPt5Ssk3gSKZ4kylWSK54kylWGJh4E+rCkGHDhrFx48b47fLycvLz89sd95e//IWbbropfru+vp6XXnqJr3/964CTyA2j60YdPWHozrw7t6kTDEWpC0aJxWxMQyfQGI2XhB2tUuJopWNSWiaE6Ey/7yGfffbZvPPOO1RVVdHY2Mjrr7/OzJkzWx2jlOLjjz9m6tSp8ft8Ph9PPvkkmzZtAuCZZ57plRWyapoobehQkOWlLhChqj6MZdnomjOcNBS22FxcAbTvaypNeIQQvaXfE3JBQQG33nor1113HZdffjnz58+nqKiIG2+8kS1btgBOqZvL5cLj8cSfZxgGjzzyCD/96U+ZO3cuH3/8MYsXLz7ueDTAMDRyMlK48oJCZ6iparo4BKehkN/niifctl3SQColhBC9Y0DqkBcsWMCCBQta3fe73/0u/t85OTm89dZb7Z43bdo0Vq5c2aux5GWl4HUb1NSHWffufmKWDcppv+k2ddJTPXjdBhW1ITYXVxAMxaiqD+MyddJ9Lnxel1RKCCF6xZC/Uq+iNoTPYxCO2tQEIrhNg5hlo2mQnuomxeNMiPa6DZa/sROP2yAcjRGL2VTVh4nGbFymLk14hBDHbcgnZNtWNDTG0DUor25E1zVspdDRqG0Io+salmVjanrTiTwDl6lTF4gQjVlEojZfnztRTtQJIY6bNKhvYje113TK8Jy95WjMJjPVzaJZJxGK2vG94xSPSUG2j5F5fnxeU5KxEKJXDPkVcku2DS7TqYnWNY3CERncds0ZAORm7G83EytR944TuaOcEKJzkpBbUHzWjD5m20wck8myFR84F4O4DYKNUcA52ZeoDdylo5wQyUu2LNqIWQo0yErz8NbWI9QEIvi8JjFbgaZhahAMxeJbGYmW5KROWojkJQm5A8pWpHjMVonNthWhcIzS6sTeBpA6aSGSl2xZdKIuEHGmPdtOrw3LVjR33EjkbYDcDG/S7HULIVqTFXIHbAWNEYuYpdA1jVjMxrYVoDBNPaG3AeZMH4Nl2YSjFkop55dKAu51CyHakxVyJ1K9JvWNUeckX9N9lg3ZPhdw9G2Agap0aDtNuje/tlRvCNG3JCG3oWtOs/r0VLdzAUgwSsyyAKc22ed1EnJX2wADXenQFx3lBvp7EmIokC2LNvwpLrLSPBypDFJdH0YDMv1uDB1cpt6tbYDBWOkwGL8nIRKNrJBbMHQNl6ERCEWbTuY5V+tFYxYel0l2mptgKHbUj+uDcXbeYPyehEg0kpCb6BqMzEvlSGUQy7bxp7gIRW1iMRtD18lOc3PPN2d067UGY6XDYPyehEg0smXRxDD0eHkbaISjNhlN+8iWbVNa3RhvUn80g7HSYTB+T0IkGlkhN4nGbIKhKKapY1k2kZhNVV0I0NA053/dPYnVl5UOA2Uwfk9CJBpJyC2U14ScFbHlFLopXUPTnLK3DL87fhKrO0loMM7OG4zfkxCJRBJyG7GY3TyGGnC2Mpongyil5CSWEKLPSEJuQwFel4GuaaT73XISSwjRb+SkXgfCEQvT0OQklhCiX0lCbqN5t2JEbiqLZp1EZqo7odttCiEGD9myaMMpeiNeQSAJWAjRXyQhd2DGqfnHlYilCY8Q4lhIQm7BbepcMmMMl557wjG/hjThEUIcK0nILfz3D88/7tdo2YQHwOMyCDfdLwlZCNEVOanXRNPo9qXRXZERSkKIYzUgCXnNmjVccsklXHzxxSxfvrzd47/+9a+54IILuOyyy7jsssvix5SUlLBo0SLmzJnDzTffTCAQ6NW4lr+x87iTcm6Gl0jMbnWf1C8LIbqj3xNyaWkpDz/8MCtWrGDVqlU8//zz7Nq1q9UxW7du5aGHHmL16tWsXr2aRYsWAXD33XdzzTXXsG7dOiZPnszjjz/ea3G5e6m/rzThEUIcq35PyG+//TYzZswgMzMTn8/H7NmzWbduXatjtm7dym9+8xsWLFjAPffcQzgcJhqNsmHDBmbPng3AwoUL2z3veIQjFrGYRUlFgGUrPuC2J95m2YoPerxiLpqQ265++ZzJw1j37v5jfk0hxNDQ7yf1ysrKyMvLi9/Oz89n8+bN8duBQIBTTjmFxYsXM3bsWJYsWcLjjz/OokWL8Pv9mKYTcl5eHqWlpb0aW1VdGE2DhlCUDL+bhlCU597cRUaGj2mnFHT4nI3bS3n577sorQpSkO1j4fmFXDhjPBfOGB9//Dcvb8Y0tW6/ZrO8vLRe/f76WjLFm0yxQnLFm0yxQmLF2+8J2bZtNE2L31ZKtbqdmprK7373u/jt66+/nttvv51rrrmm1XFAu9vHHZsCU9cwdJ2YpZw/NYvnX/+Esbm+dse3LHHzunTKq4M8/uJHra7oe/71T0Cj26/ZLC8vjfLy+l79/vpSMsWbTLFCcsWbTLHCwMTb1S+Aft+yGDZsGOXl5fHb5eXl5Ofnx2+XlJTw4osvxm8rpTBNk+zsbOrr67GaBo62fV5v0DRQqvV9XVVIdGfOnFRdCCG6q98T8tlnn80777xDVVUVjY2NvP7668ycOTP+uNfr5Re/+AUHDhxAKcXy5cuZNWsWLpeLadOmsXbtWgBWrVrV6nnHKyPVhcvQMfTWq+6uKiS6k2yl6kII0V39vmVRUFDArbfeynXXXUc0GuXKK6+kqKiIG2+8ke9973ucdtpp3HPPPdx8881Eo1HOOOMMvvGNbwBw1113sWTJEp544gmGDx/OQw891Gtx1QWjoJxeFvuO1GMaGmk+N5ZlY+oatz3xdrvLoLszZ27O9DEsf2MnYZxkHYnZWJbNxDGZLFvxgVxeLYSI05Rq+yF9cLvh/75OWXVjl8c0d3zTALdLRymwbGfYqddjcv0lEymakNtqD7llsm3bFa5tb4uJYzJ5a+uRLp8ne3F9J5liheSKN5lihcTbQ5ZLpzug684Jw5ilCEdtTEPH0HVsBYHGCC/+vTieOD2mTml1CFAUZPu4+ovtW3S27Rq3bMUHcnm1EKIdScgdsGxosQtB87ayroGlNEqrgq1Wx8NzfURizsUg3VFRG8Lnbf3Wy4k+IYT0suhE1Op4J8eptNO6VWHRGTnRJ4ToiCTkbrCVQinV9CcUZHmPq5xNLq8WQnREEvJR6BromoatnD9TvSZXXlB4zKvc5hN84YhFXUOE6vqwjIcSQgCyh3xUtoKYZWPoGumpbrwunWde34nXbRBsjAKty9m6WuW23HfOTPO0eo4kYyGErJC7wWUapKa4qKoPEQhb+LwmMVuBpmFqdHsI6vHsOwshBj9ZIXdDJGqBBhoajeEY6anueMma3+fmnmvO6NbrSHWFEKIrskLuBgVEIha65mxfNOtpMpXqCiFEVyQhd5PCKYWzbUV1XYgjVUEOlgcIhmLd7m8s1RVCiK5IQu4hpZy+F5GIhYbC4za6PfqpaEIu50weRl1DhINlDdQ1RDhn8jA5oSeEACQh94hpaDS3YNY0yMlIIT3V3e0Tc5uLK3hr6xHS/W5G5ftJ97t5a+sRmSAihADkpF6PxCyFBhg6gEaKx3n7uruX3LLKAtr3sGiuUS6pCGAppwZ6ZG6qlMUJMURIQu4hhdPrQtMUB8saME2dFLeB32setZ1mV1UWzTXKMUsRCEWbr9HmSFWQ5W/sBJCkLMQgJ1sWx0g1rWBjMZuahghV9RFqAhF8XpOaQKTDfeWuqiyaV8+N4RiapmPqOhrQGLGkVlmIIUIS8nGIWgrLdpoQhWPWUS/46KrKork3Rsyy493lNJyEL7XKQgwNsmXRSyxL0RiOtdpXLqkItNvGWDTrpFbN6pu3NnIz9lMTiGAaTcNQNWd7xGy6LFtqlYUY/CQhd8HQnf3iZs2TRDpTXR8mxWPSGI5RWRfCshShkjoy/O74NsaiWSdxW5sr+zYXV9DQGKWsKoiua1iWwjl9CCluQ2qVhRgiZMuiC1br7d4ukzFANGYTDEWprG3EshS65jQnqq4PY9uqw22M+Mk8W5GT4UXXNNDA1DV8HpNh2T7pBCfEEDHkVsi5mSnUBSKEIt2b7tHS0VbIGlBTH3bGPdlOhzhNA1tp1AUi5GeltNsLblUK5zLweV2EoxZ5WT5uvbKoxzEKIZLXkFshV9eH8ftc5GZ68biNoz+hhaOtkEfm+sjwexiW48PtNuLHN/fA6GgvuLNG92VVwR7FJoRIfkMuIVuWTUVNiOq6MF63QW6mN34i7ni1bFyf7nOhcKaNWLZC17RWe8GbiytYtuIDahsiHKkM0hiOxV8nErPJz/b1SkxCiOQx5BJyM8tW1DZEqKgJoZSzf5vmc+Eyj+0tyc3wUjQhN17aZhg6WX43uqahlCI/KyW+F9y8b1wTiJDhd2HZisraRoKhKHWBCJU1IfYfqWPZig/ksmohhpAht4fc3IuipVDEItTUXjPF6yI91U00ZhOOWERiFupoexVAXSDC5uIKiibksvdwHa9vOEgoEsNl6mSleQhFrPgJvZb7xh6XQSRqU98YpbwmhKY5lRW5Gd54ZQbIVXpCDAVDLiHnZqYAGpGoRSgcIxy1sZsyrq0g0Bgl0BhF1zVS3AYZqSkozWlSH43ZRKJ2q57IzaKWzfI3drL3cF28gVCqZVJVF6K6PoyuQWVtiB37a1A4fSzSU10ABMIxDF0jZjlbG+GYTTAUa9fr4ng198ro6vJuIcTAGXIJuby6kcraEG6XgcdtkpbqbFEEGqMEw7H4ati2FYFQjEAohqY5Y5y8HgN/igvDcK6oizbNxItZiphlU1MfYvVbe9Gbjneu4tOwlcJus8oORy3Ka5xKD0MHNB0NhaE7x9c0hMnP8vXaVXot5/m1vLwbZPUtRKIYkIS8Zs0annjiCWKxGF/72tdYtGhRq8f/8pe/8Ktf/QqlFKNGjeL+++8nIyODlStX8uCDD5KTkwPA+eefz6233trjr2/ZzlV1zSfSzKYkVZDtc6ZBByNY1mcZVClnhRyJWtQ13WcaGi7TwGXqeN06puHCMDQ0NJSysWyIxCxiloonbcu2O9z+cBbcdjw2XYNI1Ka0Kkg0ZuFxmfHtkGN1tE5zQoiB1+8JubS0lIcffpiXX34Zt9vN1VdfzfTp0yksLASgoaGBn/70p7z00ksUFBTw6KOP8qtf/Yof//jHbN26lSVLljB//vxejSlm2dQFItQ1NQfKzUghGrOpbQjHe1W0f44iZsVoDLe+X9PA0HU8Lg00HdPQ8LhcmIaGYegopbCaVtSW3fSnpYhaNrbtVGTYNJXYxWzQNNwu/bhXszLPT4jE1+9VFm+//TYzZswgMzMTn8/H7NmzWbduXfzxaDTKXXfdRUFBAQAnn3wyhw8fBmDLli2sXLmSBQsW8MMf/pDa2tpejy8YilFaFSQciZGXlUKm34Oud3AmsBNKOQk+ELIINkapbYhQVReirLqRwxUByqobqW0IE4o4DYbcLoO0VBf5WSkMz0klLzOFzDQP6T4Xaalu8ptuH2/HN5nnJ0Ti6/eEXFZWRl5eXvx2fn4+paWl8dtZWVnMmjULgFAoxG9/+1suuugiAPLy8vjOd77Dn/70J4YPH84999zTZ3EGmhKz3VSy5k9x9fg1FM4l0C0rO2xbEYnZNIZj1Aej1NSHqagJcaQySGlVIJ6sR+SnMjzbR06ml0y/h+HZvnhNs8JZiW/dU8mDz33IbU+8fdQSOZnnJ0Ti6/ctC9u20VpkKKVUq9vN6uvr+e53v8vEiRO54oorAHjsscfij3/zm9+MJ+6+opRTzhZojJLh95DvNakNRAj34LLr8SMz+NypBax4fcdRy+ds5axaidkUH6wl0++JX7QSiVpkZaRg6zq2Bjv3VfPKW3txuw3GD08jHLNZ995+srJ8nHFyQbvXvjAvjYwMHy//fRdlVUHys30sPL+Qaae0P/ZY5eWl9dpr9bVkihWSK95kihUSK95+T8jDhg1j48aN8dvl5eXk5+e3OqasrIwbbriBGTNmcPvttwNOgn7ppZf4+te/DjiJ3DB6dunzsbJsRVVdCI/bINPvIRK1qA1EsDvZX27p0wM1fHqghrQUF/WN0aMerwE+b1PHuNoQGakuDNPp+Pa5k/OorAoAsPJvn1LXGMVtOu+BrmvYts3zr31CiqGh6zqmrmGYOoauoWsaI7O8fP/K05wGRji/cMrL64/9jWkhLy+t116rryVTrJBc8SZTrDAw8Xb1C6DfE/LZZ5/Nr371K6qqqkhJSeH111/n3nvvjT9uWRbf/va3mTt3Lt/5znfi9/t8Pp588kmmTp3KlClTeOaZZ/p8hdxWOGJRWhUkPdVNflYKtQ2RVpc8t6U1/Z9SdCsZA2SmudE0jUjMIhpTVDdEKcgymHfOeE4ekwXAjv3V7C+tx1ZOv2R/iguv20QpOFwZJGYpsCwiAOHWsWiahq6BqesYpnPS0dB1DA10Q2uVrLtrc3EFb764mcPlDVLfLMRx6PeEXFBQwK233sp1111HNBrlyiuvpKioiBtvvJHvfe97HDlyhG3btmFZFq+99hoAkydP5r777uORRx7hpz/9KaFQiHHjxrFs2bL+Dh9wtjGCoRhZ6R58XpPahrCTBNtQTf+nAy6X02j+aIlO0zRqG8Jomoapg9W8jdFkx/5q/vTWHjRNQ8Op2KhtCIPfWSVnpXk6fN3mWJRyqjhilgXRz7ZetA6StaFrmE1JWm9aZWta62TdXN/scRtS3yzEcdKU6slaKPnd8H9fp6y6sddeLzXFRZrPRaAxSn2w81Wwz2OQ4fdQH4zQ0Nj5qto0mhvUfybN5yIvw8s3F0ziyTUfU9cYxbJsagMRpyWoc/0JOhoet0FBVgpfmDIivqI+Xk6LZicZG0brrZCnX/2E2kAYr8sgHLWbaryjpHldLG7TiD9RyMfqvpNMsYJsWQw6gcYojeEYmX4P+VkmtQ0RwtH2J/2CYYvMNI30VA9KOVUcHelopR0IReP71dX1YUDRELawW1axKXC7NCIxmz1H6tlftpPzTx/BF88cfdzfo1KgUKDAsltvhew+VENqiivem9Rt6njdHmIxRSRqoxs4K+wOVtdCiNYkIfcCu+VJv7Sjn/TzuI1OE3LHrw/RqN1Ut6zHmxCZupMHLdvZI45YNhpgaE5Mf/+whJF5/l5bKXck3eemrjHaVH/tfL+RmEV6iouq+tBnWyE4q2uj6QKZ+FaIpsUb+UuyFkOdJORe1PakX30g0irxRmIWLkOnpiHcxat0LGLZPPjcR011xE6S+yyJOVscOsRPyhmaUx3yz00lfZqQvzBlBH96aw/hqIWuaUQtp7/HF6aMAJpW18371rYFLXZ1Wm6F6IaGqekYppOgDU3H0I/9RKMQyUgSch9wTvpFyUzz4PO6qK4POTP1NOcxu32zuE4ZutNP2VZQVf9ZIrdsRduNkZbV3LYCw9Catjj6zsljsrgUeGdbGeVVAbLSPN3ev269FaKIYuPshzjanWiMr7A/W107JxslWYvBQRJyH4lZioqaEKlek9zMFBqCURpCToOi7tB1Z7Wr4bS0yEl3MyLXz/Z91R3uM9N0nMugqbOcIsXj6rTqojedPCaLz58+iqqmGuneoo5WFdK8utY1TF3HNJ3tEL0piRuyuhZJRhJyHwuEYoQiltOPQtcIR46+d6w1/S9mKVyGRqbfhabrHCxv4Lwpw/n7hyV0kpOJWgpDdy4uMXQtvnUw2PR4dd2ijC8UjmErJScaRcKRhNwPnBFNIdJSnNVyfdBpgt8ZBaT73Ri6Hr8SD5x8s+dwPV6PSSgc6zQpWzY0BGOMKfATtRTRmH3Mo6mSVVera0yTmqYTo7qh4dI/O9nonGQkXnfd/FpC9AdJyP1EAxojFo2RRjLTPKR4DKrr2rf31DVn37g+ECUvM6XVYy5Dp7o+TEFWCgfKGqBplQdO0tc08LpNwlGnD/O+0gb2vbETj8tg0vhsphTmEIvZvLXlMNX14R7t9w42tmqzum7S0RWNelOybr4E3dBA052TjyAJW/QeScj9RPFZWVhFTYjUFBd5WSnxq/6a2QoyUkzqG2NELbvVCjlq2fEk+v/W7XC6yGnOczQN0lPdoOA/rz6dbXur2bSrguJDtYSjFh/sLOeDneVomtOcPtVrUhuM8Ke39nApDMmk3JGOr2hsfUzzdohGiwtlDOcEo6HpUnstjpkk5H5iGlqrk3GBxijhSIysdC9et0lNQzhet2yYBgXZboKhKBHAilk0NMaIWQqjaUWcn5VCVV0IWykMQ8ef4kbXNdKb+lqccVIeZ5yUR0NjlK27K9m0q5J9pfUo9dlQV0N3mt+/ufGAJOQeaN4OgRYXyrTQMmG3L+drqgzp5FJ0MbQNrY3FAaLrGiPz/O3uj1mK8upGIlGL/KyU+EQPy7K54vxCLj1nPIYG1Q1Rok3JvKouxIt/L+a0E7Lx+1xkpXvJSfei61qr+t9m/hQXMyYN46bLJjlN730uTMP5sTuXOVscrAjyyB838fcPD1FVJxNEjpdSxKe/RKM2jZEYDcGmYQX1YSrrQlTUhCivaaSiLuS0eA1FaYzE4kN0baWaEvtAfzeiP3Vrhfzggw/y5S9/mdGjj/8y3KFIA0qrgp0+3tAYJRSJkZXmJd3nYs70MUyekEtVVYB17+5D15rLvJy9z2AoxvqPSnC7DKyQhWFo3epfkZvhpa4xSprPTbSpSX4wHMO2FWXVjby+4QCvbzjA6Hw/UwpzOO2EHNJ87t5/Q4a4DitEWmh5daNu6OguV3wyuZx0HNy6lZCVUixatIjx48dz1VVXcfHFF+N2yz/U7rJshXWUpvYxS1FZ28iXzjuBs04piK+MKmvDTjVA/B+fc1Ve1FLkZrpbXRl3tG2H5qvqIjgnCL0eA9OAsycPp7YhwubdldQHoxwoa+BAWQN/fmcfE0ZkUDQhh0njs+PN8kXfalkhgm3RGIlRH3A2Rjpro2rocsHMYNDtbm+2bfPPf/6Tl19+mQ8//JA5c+bEm8cnk97u9tZbmvcSXabOgs+PZdG8Sbz57l5+s/rjeCWGBq26wGWluUnxuOK9I765YNJRv86O/dX8c1NJh1UWtq3Yc6SOTbsq2bq7klCLXyKGrnHymEymFOYycUxWuzK67OzUXr8wpK8kU6zQs3ibL5ihqVqn+YKZ/jrhKN3euvc1O9PtJY+u60yePJk9e/awd+/eVlM/hEPXmq+S6zmlwGVo2LZizdv78PrcvPXRIYblpFBd71RitH3pmoYI9YEItoKauhB/ff8Au0vquixpO3lMVqcraV3XmDAigwkjMrj0nHF8eqCGj3ZV8sm+aqKWzba91WzbW43HZXDquCyKJuRQOCoDQ5dTEYmi5XaI3cUFM80nHNvWYOtNVz5KSd/A6NYK+Y033uCll17io48+Ys6cOVx11VVMmnT01Vgi6ssVsselE4523qii7Qq3LbfZfLLNxutxTsChFFbTP66a+vZ1ywCG3vwP0emdnOp1xbcyLm0xaaQ7OlpBjxuezva91WwqruDTA7VODW8Tn9fktBNyOHfqSLJ8rvjWSiIbzCvk49XyCseOxoAZOhi6Ht9Sa5s9ZIXcva/ZmW6tkJ966im+/OUv88gjj+D1ytj4znSVjKHrZAwQs50LPUxDJxiKEotZzgUfgMdtkpfV8VV+iuaZeopQxMKf4sZtGkSgR93emqeRGIaO12NS1xh16pTPGc/pJ+Zy+om5BENRtuyuYnNxBXsO1xMMxXh3WynvbislI9VN0YQcphTmMjzH1+HwWpHYWu1fd3MMWMv962jM6UbY/KFJVtg9062E/Nxzz1FdXc0777yDaZpMmTKF9PT0vo5tyLFthY3Ctp2/yOGoHT+h11yJkZnmwZ9iUlkbitc1N/e9AIi2GPfUfGVfd/1zUwmG8dnl2h0ldZ/XxfRTC5h+agG1DWE2N9U4l1QEqA1E+Ofmw/xz82HyMr0UTchlSmEOuRkpXXxVkSy6MwYMw6C6thFdx6m/NpqudJQKkW7pVkJev349t912G4WFhdi2zYEDB3j44Yc566yz+jq+IanlroRSYDX9rY1Ziuq6EP4UF7mZKR2OjXIu/IjhdZvxK/u6q7o+jLdNJUVXST3D7+ELRSP4QtEIIgr+8f4BNhdXUlEborwmxJvvH+TN9w8yKi+Vogm5FE3Ica4mFINS87aZbTuLihhtEnbT/zWvsA29aXVttr8sfaheNNOthPzII4/wzDPPcOKJJwLw8ccf85Of/ISXX365T4MT7Vk21Aai1Dc6Y6Pyskyq68LQopa1tiGMnt7xhSJdyUrzUNcY7fBy7aMZlpPKRdNGc+GZoyipCLCpuJLNxZXUBSIcLA9wsDzAq//ex/gR6UyZkMOk8TnxC2HE0NDhChviWyIgl6V361+EpmnxZAwwadIkhths1ITTPDYqxWOSm+lttVq2bDA0jXmdnNDrrPStbZ1y2+kf3aFpzlWJI/P8zJk+hr2Hm8ro9lTRGI6xu6SO3SV1/OmtvZw4KpMphTmcMjYLt8s4+ouLQa9Hl6U3lfW1vSxdS+I67C4Tck1NDQCTJ0/mqaee4uqrr0bXdV5++WVmzJjRH/GJo2gMxwhHLTL9HkbkuKisC2PbNqles9Nk3OGJOz6b/tFZnXJP6ZrGCSMyOGFEBgvOGceug7VsKq5g+95qIjGbT/ZX88n+atymzsSxWZxemEvhqIz4pd1CtNU6YXdS1hevw27aFjHaXjgDuqaTiNWaXSbkGTNmoGla/A34xS9+EX9M0zR+9KMf9W10Q4ChOyvanmpZQte8Wk5PdZOb6UHX2o9ual4V7y9tQNMgLdWNZmrtTtw1J9/mpPzPTSXA8XeDMw0n6U4cm0UkarF9XzWbiyvZeaCGSMxmc9MWR4rHZPL4bKYU5jJueFpSlNGJxNG6DhuiWED7STNoTiWIMkzqg9GEOfHYZUL+5JNP+iuOIetYLyTpSEMwgqFDfraPgixffI+t5arYVgplQ1VdGE0L4zJ1Ur1mPIEfbQXdG9wugymFuUwpzCUYivHx3io27apgT0kdjeEYGz4pY8MnZaQ3l9FNyGFEbqqU0Ynj1jZhR2IWgdBnJ8a7c+Kx+eKZvtjH7tYeciQSYf369QQCTnG6ZVns37+fW2+9tfciGYJyMzzUNkSJHsMSuaO/A7aC2oYIlqWYf+mpZPk91AUjrcrZdP2zNqBKgWUpahsi8Wb43Sl9600+r8lZE/M5a2I+dYFI00q5goPlAeoCEf61+TD/2nyY3AxvvMa5beN+IXrLQJ947FZCvvXWWzlw4ADl5eWceuqpbNq0ic997nPd/yqiQzFLMWFkOvtLGwhFnEujj+W3bfNvdaWcpOx2GUwalwNAdroXwzTweRQxS3VwMla1+rOnpW+9KT3VzblFwzm3aDiVtSE2FVewaVcF5TUhKmpD/PWDQ/z1g0OMyE1lyoQciibkkOHv+yGuQrTUrROPbQbwtjzxGI5aeDo5id2tbe3t27fz8ssvc+GFF3L77bfz7LPPUltbe8zf0Jo1a7jkkku4+OKLWb58eYdfb+HChcyePZs77riDWMyZqFFSUsKiRYuYM2cON998c3zFnqwCoRhzpo/BNJwfnjNDT6erD+Z5mV6G5fiYOCYTV9O+l8vUcRnOcw1do7bhs+Spaxopbh2/z+1cig2tXl83dDJSXUSarjLMSvO0W7H3tJ65N+RkePniGaP4/lVT+P8WnsbMKcPJaKphLqkI8Oq7+1m24kN+u+Zj3t1WSjDU+YxCIfqTszBq6ocda98POxDs/O9qtxJyfn4+pmkybtw4du7cyYknnkh9/bFd/11aWsrDDz/MihUrWLVqFc8//zy7du1qdczixYu58847ee2111BK8cILLwBw9913c80117Bu3TomT57M448/fkwxHCsNpxubaRz/XqZpaKS4DYom5DIiN5X0VA+moWErhWnqzqWpLb6OBmSkuvB5XbhNnYraEDT9Fm4VY/NyuYUvnjGKiuogtlKMyPXj9RgYOmSne8jLSMEwjXjC/cKUEViWTSRmoZQiErN6XPrWmzRNY0RuKnOmj2XxNVP51qWnMv3UAnxeEwXsPVzP6n/t4Wd/+ID/XfcJH31aQTjadatTIRJVtxKyz+djzZo1TJw4kVdffZUdO3YQDHbecL0rb7/9NjNmzCAzMxOfz8fs2bNZt25d/PFDhw4RCoU4/fTTAVi4cCHr1q0jGo2yYcMGZs+e3er+/qRwGgBlpXm4/Nxx6MeYlzNSnSvtRuSmAsRXyZlpHkbkppKT4SXD72FYVkp85ZuXlUJmmtNHJBKzyc3wUpDljf82Vko1/QkFWa37jRRNyOWaWSdhaBrBcJQsv5eReWl43Ua7hHvymCwuPWc86SkuQuEY6SmuHjco6iu6pjFuWDqXnTuepV89g6/NOZmpJ+bidjknK3fsr+GFv+3iZ//vfZ79y0627a0idiwlLEIMkG7tId9555288MILLF68mBdffJGvfvWr/OAHPzimL1hWVkZeXl78dn5+Pps3b+708by8PEpLS6mursbv92OaZqv7+0vz5rxSsGjWSRRNyOWT/TV8erDG6T3R5ni3qZPhd3PO5GFs3FFOSUUA09BJT3VhmgaWZTNn+hjASZgA697dT0VtiNwML1+5YiJjc31sLq5g+Rs70XWtacVqt3ru//x5O6GI04TI0DW8HoMrLyhsF79z6bLzdTYXV/CPTSUEQjF04HOnFrRKuF216EwUhq7H47wsZvHJvho2F1ewY38NUctmy+4qtuyuwus24mV044enox/rb1Eh+kG3EvK4ceO47bbbOHToEEuXLqWgoOCYv6Bt263Kl5RSrW539njb44B+K4PSNThhZAZKKRoaY1w4YzwAX7l4Io8+94FTx6g5heq2co4fkefna/NOZdopBdwAbNxeyst/30VZVZC8LB8Lzy9k2imfvY8X5qXFX7elC/PSyMjwxZ+bn936uV091pmWXyscjlEbiBxTpUdb2dmpx/0ax2pYfjrnnzWGYCjKRzvL2bCtlE/2VRGKWGzcUc7GHeVk+N2cObGAs04tYNzw9KQqoxvI97ankilWSKx4u0zIP/jBD3jooYc4fPgw3/nOd9i/fz8AEydO5MEHH2TYsGE9/oLDhg1r1dy+vLyc/Pz8Vo+Xl5fHb1dUVJCfn092djb19fVYloVhGO2e15f8KS6iMZtw1CLL7473Tx2b6+Prcyfy4t92UVodwtBhRLaPK8+fEF+Ntjz21iuLWr1uV31YW/Zp7eq5PX3dzkTDMQLBaKtexz2RSD2GJ47KYOKoDOqDEbY0daM7UNZAbUOEv248wF83HiA73cOUCbkUFeZQkOUb6JC7lEjv7dEkU6wwMPGOHJbR6WNdJuQ9e/YAcP/99zN//nxuuOEGlFKsWLGCO+64g6eeeqrHwZx99tn86le/oqqqipSUFF5//XXuvffez4IdORKPx8P777/PmWeeyerVq5k5cyYul4tp06axdu1aFixYwKpVq5g5c2aPv35P+TwGmWkewlGr1VZBs5ZbAclKA1I9Jh6XQV0gQmSQnBRL87k5e/Jwzp48nKq6EJuLK9m0q4LS6kaq6sL87cND/O3DQwzP8cWTc6aU0YkB1K0tiz179vDLX/4ScLYJFi1axIoVK47pCxYUFHDrrbdy3XXXEY1GufLKKykqKuLGG2/ke9/7HqeddhoPPPAAP/7xj2loaGDSpElcd911ANx1110sWbKEJ554guHDh/PQQw8dUwwdcZtaUw2wRkGWN74P23Jfd870MUmffLti6hpZaR6Cx7laTkTZ6V7OnzqS86eOpDGm+McHB9i0q4KahgiHK4McrtzPuvf2M3ZYGlMm5DD5hBz8Ka6BDlsMMV2OcLrgggt4+OGHefLJJ7n11luZMGEC4FRCfPvb32bNmjX9FmhvqaxswO7N65X7yECPwolZirpg91fLyfRRtTlWpRQHyhr4aFcFW3ZXtZrEomswYWQGUwpzOXVcFl73wLUKTcb3NlkM1JZFdkbHk5e6/Fs2e/Zs/uu//ott27YRjUb5zW9+w+rVq/nZz37Gj3/84z4JViQG0xi8q+VmmqYxpiCNMQVpzPv8OHaX1LJpVyUf76kiHLX49GAtnx6sxTQ0Jo7Joqgwl5NHZ7abuC1Eb+nWkNNIJEJZWRmjRo1ix44d6Lreqj9yMpEVcs/FbHXUveVkWhkdLdZozGbHgRo27apgx/7qeO8PAI/LYNL4bKYU5nDCiIz4dOaBjDeRJFOskGQrZIB3332XkpISvvCFLwBw8sknA7By5UquuOKKXgxTJKrBvLfcEZepM3l8NpPHZxOKxNi2t5pNuyooPlRLOGrxwc5yPthZTmqKi6ITcphSmMPofH9SldGJxNTlCvn3v/89zz33HGPHjuWjjz7il7/8Zbwx/RVXXMHKlSv7LdDeIivk4xOzFfWBSLvLk5NpZXSssTY0Rtmyu5LNuyrZV9r6Z5OV5ol3oxuW3btldEPhvR0oSbVCfumll3jppZfw+/384x//4NZbb+Xpp5/m5JNPlhFOQ1R8tdzUMCUZfrn1Fn+Ki89PGsbnJw2juj7M5uIKNhdXcrgySHV9mPUflbD+oxIKslKYUugMdc1O7/gfnhAd6TIhm6aJ3+8HYObMmfzoRz/i5ptv5qWXXpKPZ0Ocz23iMXXqg1FCkcFRt9wTWWkezjt9JOedPpLS6iCbd1WyqbiCqrowpdWNvL7hAK9vOMDofD9TCnM57YRs0nwycVt0rcuEnJ2dzUsvvcT8+fPxeDxcfvnlfPLJJ3zrW9865m5vYvAwdJ1Mv5vGiNWu69xQUpDlY9ZZPi6aNoqD5QE27apgS3El9Y1RDpQ1cKCsgT+/s5cJIzKYUpjDpPHZA1pGJxJXl3vI+/fvZ8mSJVx11VWtTuD9/Oc/53//93/Ztm1bvwTZm2QPuW9kZPrYd6iaUDjxV8v9sW9o24rdh+vYtKuCj/dUtfoUYRoaJ43OZEphLhPHZB21jC6Z9mWTKVZIvD3kbpW9daS0tPS4mgwNFEnIfcOJt45Q1KY+EMFK4Pe4v/8RxiybnU1ldJ/sq2nVyMnjMjh1XBZTCnOZMLLjMrpkSnLJFCskXkLu1uempUuXtrqtaRopKSmceOKJXHXVVRhGx+NIxFCj4XUZuNK9NDRGCYVjHc7+G2pMQ+fUcdmcOi6bcMSZuL1pVwWfHnTK6D78tIIPP60g1Wsy+YQcTi/MZXSBXyZuD0Hd3sjatm0bV1xxBbqu88orr5CTk0NpaSmffvopP/nJT/oyRpFkDF1z9pbdBg2BCLEEXi33N4/b4PQTczn9xFwCoShbd1exqbiCvYfrCYRivLutlHe3lZLpd8fL6LISvBud6D3dSsjFxcUsX748XnFx1VVXcf3117NixQrmz5/fpwGK5KQUeF0G7gyvU4khq+V2Ur0upp9awPRTC6hpCLOluJJNxZWUVASoaYjwj02H+cemwwzL8cWb7OdIGd2g1q2EXFtbG0/GAF6vl4aGBjRNw+WSjliic7omq+XuyPR7+MKUEXxhygjKaxrZtKuCTcWVVNaGOFIZ5EhlkL9sPMiovFSnjG5CDulSRjfodCshT5kyhR/+8IdceeWVKKV46aWXKCoqYv369aSkpPR1jCLJyWq5Z/IyU7ho2mguPHMUhyoC7DxYx3sfH6YuGOVgeYCD5QHWvrOP8SPSmVKYy+Tx2aR4pIxuMOhWlUUoFOLxxx/nH//4B4ZhcP7553PTTTfxl7/8hRNPPDGpGg1JlUXf6G68mgaNEWtAV8vJWAlQUdHA3iN1bNpVydY9lTS2KC80dKeMrmhCDqeMzcLtGriT7Mn43iZSlUW3EvKCBQs6vF/6IfedwZqQm9lKNV3lF6O/r8JP9qQRs2x2Hazlo10VbN9XTTT2WRmd29Q5pamM7sRRGRh6/7YKTfb3tj8cd9lbyyqKaDTKn//8Z0aPHt070YkhqXlvORQxqJe95R4xDZ2JY7OYODaLSLS5jK6STw/WEInZbNrlzBH0eUwmn5BN0YRcxg1PkzK6JHBMF4Yopbj66qt5/vnn+yKmPiUr5L5xPPHaStHQGKUx3D+r5cG6iguGYny8x6nU2FNS12qfPiPVzWlNZXQjcnx91otmsL63vem4V8htVVdXU1ZWdlxBCdFM1zQyUt14XSZ1jWFiscT/hZmIfF6Ts04p4KxTCqgNRJrK6Co4VB6gNhDhX5sP86/Nh8nN8DKlMJcpE3LIzZST8omkWwm57R5ySUkJX/nKV/okIDE0KQVul06Oy0tDY4xgKNrve8uDSUaqm3OLhnNu0XAqahubtjEqqKgNUVEb4s33D/Lm+wcZmZtKUWEORSfkkCETtwdct7Ys3nvvvc+eoGlkZ2fHB54mG9my6Bu9HW8kZlMfjLQ6YdVbhurHaqUUhyuDbNrl9HGuDUTij2nAuOFpFE1wWoX6vMd2fcFQfW97ok+aCyUrSch9o6/ibQhFCTT27mpZkoazb7/vSD2biyvZUlxJMByLP6ZrGieOdiZunzI2C08PyujkvT26Xt9DFqK/+L0uZ2852PWQVdEzuqYxfng644enM//ssew66Ezc3ravikjUZsf+Gnbsr8Fl6pwy9rMyOtOQidt9SRKySHimMbSGrPY3Q9c5eUwWJ4/JIhKz2LG/eeJ2DdGYzebiSjYXV5LiMZg0PocpE3IYPzwdvR8mbg81kpBFUtCAVI+Jx2V0OGRV9A63aXDaCTmcdkIOjeEYH++pYnNxJcUltTSGLTZ+UsbGT8pI8zVP3M5lZF6qjHTrJZKQRVIxdY3MNA+N4RgNjUNryGp/S/GYTJuYz7SJ+dQHI2zZ7VxwcqCsgfpglLe2HuGtrUfISfdSVJjDlAm5ZGenDnTYSa3fE3JJSQmLFy+msrKS8ePH88ADD5Ca2vqHWFZWxtKlS6moqEDXdW677TY+//nPE41GmT59equrBF9++WVpkD/EaIDPY+Jx6dQHooRktdzn0nxuzp48nLMnD6eyLhQf6lpW3UhlXYi/fXCIv31wiNH5fiaNy6aoMIdMKaPrsX6vsrjpppu49NJLmTdvHo899hjBYJDFixe3OuaHP/whU6dOZdGiRezevZtrr72Wf/zjH2zfvp2HH36Yp5566pi/vlRZ9I2Bi1fRGLGoD3Z/tSyVAL1DKcWRqiCbi50a55qGSKvHxxakMaUwh8kn5OBPScw2vYlWZdGvCbl5hfvee+9hmiaHDx/mq1/9Km+++War41577TVmzJhBRkYGlmUxffp0/va3v/HKK6/wwgsvoGkapmnywx/+kM997nM9ikESct8Y6Hht26YuGG01TLQziZrgOpMM8Sql2F/awI6DtWzYdoRAqGUZHRSOymTKhBxOHZeNx504n2gTLSH365ZFdXU1fr8f03S+bF5eHqWlpe2Omz17dvy/n3rqKU455RTS0tLQNI0LL7yQm266iU8//ZQbb7yRNWvWkJ2d3W/fg0hMuq47zYqSYMjqYKRpGmOHpTH11GFceMZIig/Vsrm4go/3VBOOWuw8UMPOAzWYxm4mjs3i9MJcThqdKWV0bfTZCvnVV1/l/vvvb3Xf2LFj2b9/P+vXrwcgFosxdepUtmzZ0uFrPP300/zhD3/gmWeeYfjw4e0ev/nmm/nSl77ERRdd1PvfgEha0ZhFbYNUYiSCSNRia3ElG7YdYUtxJbEWE7dTPCZTT8rjrEnDOHlM1pApo0txm/2/Qp47dy5z585tdV/zloVlWRiGQXl5Ofn5+R0+f9myZaxfv57ly5czbNgwAFatWsUZZ5zBmDFjAOdjUk9HSMmWRd9IvHgVdtTusBF+MmwBtJRM8XYU67j8VMblT2DB2WPZtteZuL3rUC2N4RhvbznM21sO409xOd3oJuQwOt/fb2V0A7Vl0Zl+3bJwuVxMmzaNtWvXsmDBAlatWsXMmTPbHff000/z7rvv8uyzz5Kenh6/f8eOHXz00Uf89Kc/Zffu3Wzfvp0zzzyzP78FkTS01mOjBqARvmjN6zY546Q8zjgpj/pgJD5xe39pAw2NUd7ZeoR3th4hK83DlKZWoQXZQ2vidr9XWRw6dIglS5ZQWVnJ8OHDeeihh8jIyODZZ5+lrKyM733ve3zuc5/D7/e3Ssa//e1vSU1N5fbbb2f37t1omsYdd9zBjBkzevT1ZYXcNxI5Xk2DUMSKN8JPphUnJP8K+Wiq60NNlRqVHKkKtnpsWLaPKYU5FE3IISut9yduJ9pJPWkulKASOcF1JBnibW6E7/V5kibBweBPyC2VVgXjE7er68OtHhtT4GfKhFwmn5BNWi9N3E60hCxX6okho7kRforfQ11tkJiV+L+Yh5qCbB8Xf24Ms84azcHyBj7a5XSja2iMsr+0gf2lDbzyzl4mjMhgSmEOk8Zn43UPnjQ2eL4TIbpBKUj1ushp2lvur7FRomc0TWN0fhqj89OYN2Msuw/XsXlXBVv3VBGKWOw6VMuuQ7Ws/tceTh6dRVFhDhPHZOEyk7uMThKyGJI0NNJ9brxuk7qgjI1KZLquUTgyg8KRGVx67nh2HnC60W3fV03MUny8t4qP91bhcRmc2jRxe8LIDIwkLKOThCyGNLepk5MuY6OShWnonDoum1PHZROOWGzbV8WmXZXsOlhDOGrx4acVfPhpBalek9OautGNLvAnzcRtSchiyNPQSEtxOa09+2hslOh9HrfB1BPzmHpiHoFQlC27K9m8q5K9R+oJhGL8e1sp/95WSqbfTdGEXKYU5jAsu+8mbvcGSchCNImvlvtgbJToW6leFzNOHcaMU4dR0xB2murvqqCkMkhNQ4R/bCrhH5tKyM9KoaipxjknvffL6I6XJGQh2vB7m1fLURkblYQy/R5mThnBzCkjKKtpZHNTGV1lbYiy6kb+svEgf9l4kFF5qXy+aAQThqeR3ktldMdLErIQHXAZOll+D8FwlEBjTMZGJan8zBQumjaaC88cxaGKAJt2VbCluJK6YJSD5QH++OanaBqcMCKdKRNymTQ+mxTPwKVFSchCdELTnI/CHhmymvQ0TWNUnp9ReX7mTh/LniN1zlDXvVUEQzGKD9VRfKiO1f/aw0mjM5lSmMPEsVm4zf5tFSoJWYijkCGrg4uua0wYkcGEERl8fcEk3t1cEi+ji8Zstu+rZvu+atwunVPHZjOlMIfCURkYet/XOEtCFqIbZMjq4GQaOqeMzeKUsVlEohbb9znd6HYeqCUStfloVwUf7arA5zGZfEI2UwpzGTssrc/K6CQhC9EDMmR18HK7DKYU5jKlMJdgKMrWPVVs2lXB3sP1BMMx3ttexnvby8hIdccrNYbn9G4ZnSRkIXpIhqwOfj6vi8+dUsDnTimgNhBhS9PcwEMVAWoDEf65+TD/3HyY3AxvUxLPITcj5bi/riRkIY6Roetkprl7PGRVJJeMVDfnFg3n3KLhVNQ0sqkpOVfUhqioDfHm+wd58/2DjMxNpagwh6IJuWSkHlsZnSRkIY6LRorbxGPq3R6yKpJXbmYKF545ii+eMZLDlU6r0M3FldQGIhyqCHCoIsC6f+9n3PB0Z+L2+Bx83u6nWUnIQvQCGbI6tGiaxojcVEbkpjJ7+hj2Haln064Ktu6uIhiOsedwHXsO17Hmrb2cOCqDKYW5nDI2C7er6zI6SchC9BpnbJQr3Ut9Y4Rw2ELS8uCnaxrjh6czfng6C84Zx66DtfEa50jM5pP9NXyyvwaX6VR0/OeiaZ2+liRkIXqZoWtkproJuTsesioGL0PXOXlMFiePySISs/hkXzWbdlWy80AN0ZjN5uLKLp8vCVmIPiFDVoc6t2lQNCGXogm5NIZjfLynShKyEANJ1zRnbzlixIesiqEnxWMybWI+0ybmd3lccs87ESIJKAUel0F2hhef1ySB2/GKASYJWYh+0jxkNcvvwTQkK4v2JCEL0Y+Uci7RzcnwkpriktWyaEUSshADoHlsVFaaF9OUrCwckpCFGEDNY6NktSxAErIQA67latllyj/JoUx++kIkiObVst8nq+Whqt/rkEtKSli8eDGVlZWMHz+eBx54gNTU1FbHHDp0iPnz5zNmzBgAcnNzeeqpp1BKsWzZMv72t7+h6zr33nsvZ555Zn9/C0L0KRmyOnT1+wr57rvv5pprrmHdunVMnjyZxx9/vN0xW7duZcGCBaxevZrVq1fz1FNPAfDaa69RXFzM2rVreeyxx1i6dCmxWKy/vwUh+lzzkNU0nwtZLA8d/ZqQo9EoGzZsYPbs2QAsXLiQdevWtTtuy5Yt7Ny5k8suu4zrrruOHTt2ALB+/XouueQSdF1n/PjxDB8+nA8//LA/vwUh+k3zkNWcdO9Ru4SJwaFftyyqq6vx+/2YpvNl8/LyKC0tbXecx+Ph0ksv5eqrr+af//wn3/3ud1m7di1lZWXk53926WFeXh5HjhzpUQw5Of7j+yb6UV5e2kCH0CPJFG8yxQpw4rgc6oMRAqFowvfEyM5OPfpBCSSR4u2zhPzqq69y//33t7pv7Nix7eZPdTSP6j/+4z/i/33eeefx4IMPsnv3bmzbbnW8Ugq9h5NgKysbkmKyQ15eGuXl9QMdRrclU7zJFCs48VZWNgCg2YqGBB6ymp2dSlVVYKDD6LaBiHfksIxOH+uzhDx37lzmzp3b6r5oNMr06dOxLAvDMCgvL2+14m32hz/8gfnz55OVlQU4idc0TYYNG0ZZWVn8uIqKig6fL8RgJUNWB7d+3UN2uVxMmzaNtWvXArBq1SpmzpzZ7rgNGzbw4osvAvDee+9h2zYnnHACM2fOZM2aNViWxb59+9i7dy+nnXZaf34LQgy45iGrOekevLK3PKj0e9nbXXfdxZIlS3jiiScYPnw4Dz30EADPPvssZWVl3HLLLdxxxx0sWbKE1atX4/F4ePDBB9F1nTlz5rB582YuvfRSAO677z68Xm9/fwtCJAQZsjr4aEol+imC3iV7yH0jmeJNplihe/Hatp0QQ1ZlD/noRg7LIDuj44WkXKknxCDQPGQ1M82DoUvlcrKShCzEoOGMjcpO9+L1GHJBSRKShCzEINM8ZDUjzYMpq+WkIglZiEGpabWc4SXFI2OjkoUkZCEGseYhq5l+WS0nA0nIQgxyMmQ1eUhCFmKIkCGriU8SshBDiAxZTWySkIUYgmTIamKShCzEECZDVhOLJGQhhjgZspo45N0XQgAyZDURSEIWQrTi97rIlrFRA0ISshCinZZDVnVZLvcbSchCiA41D1mV1XL/kYQshOiSaWhkpXlIS3XLarmPSUIWQhyVBqR6TLIzvHhktdxnJCELIbqtechqeqobXZoV9TpJyEKIHpEhq31HErIQ4pg0D1nN8Mtqubf0+9RpIcRgopHiNvGYOnXB6EAHk/RkhSyEOG7NQ1bTU90yZPU4SEIWQvQSjTSfW4asHgdJyEKIXiVDVo+dJGQhRB+QIavHQhKyEKLPyJDVnun3KouSkhIWL15MZWUl48eP54EHHiA1NbXVMd/+9rc5fPgwALZts3PnTl588UUmTpzI9OnTGT16dPzYl19+GcOQWkghElXzkFVXhpeGxiiN4RhKDXRUianfE/Ldd9/NNddcw7x583jsscd4/PHHWbx4catj/vu//zv+348++iinn346p512Glu3bmXq1Kk89dRT/R22EOI4NQ9Z9boM6oIRYpZk5bb6dcsiGo2yYcMGZs+eDcDChQtZt25dp8fv3r2bVatW8aMf/QiALVu2UFVVxcKFC/nyl7/Me++91y9xCyF6hwxZ7Vq/rpCrq6vx+/2YpvNl8/LyKC0t7fT4xx9/nBtuuAG/3w+ApmlceOGF3HTTTXz66afceOONrFmzhuzs7G7HkJPjP75voh/l5aUNdAg9kkzxJlOskFzxdjfWfKAxHKW2IYxl921MXcnOTj36Qf2kzxLyq6++yv3339/qvrFjx6K1+ZXY9naz2tpa3nrrLe677774fVdffXX8v0899VSKior44IMPuOiii7odV2VlA7ad+B+V8vLSKC+vH+gwui2Z4k2mWCG54j2WWHUUwcYYwVC03/eWs7NTqaoK9OvXHDkso9PH+iwhz507l7lz57a6LxqNMn36dCzLwjAMysvLyc/P7/D569evZ+bMmXg8nvh9q1at4owzzmDMmDEAKKVwuVx99S0IIfpB85BVj8ugPhghGhvA5fIA69c9ZJfLxbRp01i7di3gJNiZM2d2eOxHH33EtGnTWt23Y8cO/ud//gdw9pe3b9/OmWee2bdBCyH6hQxZHYA65LvuuosXXniBSy65hI0bN/L9738fgGeffZZHH300ftyBAwcoKCho9dzvfve7VFVVMX/+fG655RZ+/vOfx/eXhRCDw1AesqopNbQqAmUPuW8kU7zJFCskV7y9GatSEAxHCTTGsPsoTQ3UHnJ2hrfDx+RKPSFEQhqKQ1YlIQshEtpQGrIqCVkIkfCGypBVSchCiKQx2IesSkIWQiSVlkNWB9tqWRKyECIpGbpOVpqH9EE0ZFUSshAiqfncJrnpHrzu5F8tS0IWQiS95iGrmWmepB6yKglZCDFINI2NSuIhq5KQhRCDSjIPWZWELIQYhJJzyKokZCHEoJVsQ1YlIQshBrXmIavZGV583sReLUtCFkIMCc1DVrP8XkwzMbNyv0+dFkKIgeIMWdXJcXlpaIwNdDjtyApZCDHkNI+Nyk734Eqg1bIkZCHEkJXicfotp6YkxtgoSchCiCGtebWclebFZQ5sSpSELIQQOENWs9O8pPlcA9YIXxKyEEI0GeixUZKQhRCijYEaGyUJWQghOjAQY6MkIQshRBf6c2yUJGQhhDiKlmOjvH24WpaELIQQ3WToOplpbjL6aGyUJGQhhOgRjZTmsVGe3l0tD1hCfuSRR/jVr37V4WORSITFixczd+5crrjiCoqLiwFQSvHzn/+cOXPmcMkll/D+++/3Z8hCCBGn6zqZqb07NqrfE3J9fT233347v//97zs95g9/+AMpKSm8+uqr3H777SxduhSA1157jeLiYtauXctjjz3G0qVLicUSr0GIEGKo+GxsVIrHPO6xUf2ekN98803GjRvHN77xjU6P+fvf/86ll14KwFlnnUVVVRUlJSWsX7+eSy65BF3XGT9+PMOHD+fDDz/sr9CFEKJDhq7Fh6weTyP8fk/Il19+Od/61rcwjM73XsrKysjLy4vfzsvL48iRI5SVlZGfn9/ufiGEGGi90Qi/z/ohv/rqq9x///2t7jvhhBN4+umnj/pcpRRai+9GKYWu69i23eH9PZGT4+/R8QMpLy9toEPokWSKN5liheSKN5lihb6JtwAIhKLUB8JYdvef12cJee7cucydO/eYnltQUEBZWRljxowBoKKigvz8fIYNG0ZZWVn8uOb7e6KysgHbVscUV3/Ky0ujvLx+oMPotmSKN5liheSKN5lihb6PV1OKUGOUxnAM1ZR2Rg7L6PT4hCx7O++881i9ejUAGzduxOPxMGLECGbOnMmaNWuwLIt9+/axd+9eTjvttAGOVgghOtbTsVEJM8Lp2WefpaysjFtuuYVrr72WO++8k3nz5uF2u1m2bBkAc+bMYfPmzfETfvfddx9er3cgwxZCiC61GxvVRV7WlFKJ//m9F8mWRd9IpniTKVZIrniTKVYYmHgzs3y4zI6LGhJyy0IIIQarzpIxSEIWQoiEIQlZCCEShCRkIYRIEJKQhRAiQUhCFkKIBCEJWQghEoQkZCGESBCSkIUQIkFIQhZCiAQhCVkIIRKEJGQhhEgQCdPtrb/0xejuvpJMsUJyxZtMsUJyxZtMsUJixTvkur0JIUSiki0LIYRIEJKQhRAiQUhCFkKIBCEJWQghEoQkZCGESBCSkIUQIkFIQhZCiAQhCVkIIRKEJGQhhEgQkpD72K9//WvmzZvHvHnzWLZsGQBLly7l4osv5rLLLuOyyy7jjTfeAGD79u0sXLiQ2bNnc8cddxCLxQAoKSlh0aJFzJkzh5tvvplAIABAXV0d3/rWt5g7dy6LFi2ivLz8uOO99tprmTdvXjy2TZs28fbbb7NgwQIuvvhiHn744fixAxnvH//4x3iMl112GWeeeSb33HNPQr63DQ0NzJ8/n4MHDwL0+fsZiURYvHgxc+fO5YorrqC4uPiYY33++eeZP38+CxYsYOnSpUQiEcD5e33BBRfE3+fly5f3e6wdxdvXP//jjfeolOgzb731lvrKV76iwuGwikQi6rrrrlOvv/66mj9/viotLW13/Lx589SHH36olFJq6dKlavny5Uoppb71rW+pV155RSml1K9//Wu1bNkypZRSd999t/rNb36jlFJq5cqV6pZbbjmueG3bVueee66KRqPx+xobG9V5552n9u/fr6LRqLr++uvV3//+94SIt9nOnTvVrFmzVGVlZcK9tx999JGaP3++mjRpkjpw4EC/vJ9PPvmk+slPfqKUUuq9995TV1111THFunv3bjVr1ixVX1+vbNtWt912m/r973+vlFLqpptuUh988EG71+ivWDuKVynV5z//44m3OyQh96GdO3fG/xIo5fyQn376aXXGGWeoG264Qc2fP189+uijyrIsdfDgQXXhhRfGj92wYYO69tprVSQSUVOnTo0nyZKSEvXFL35RKaXUBRdcoEpKSpRSSkWjUTV16lQViUSOOd5du3apc889V1177bVqwYIF6g9/+IN699131XXXXRc/ZuXKlWrJkiUJEW+za665Rq1du1YFg8GEe29vv/12tWHDBnXBBReoAwcO9Mv7+dWvflVt2LAh/loXXnihOnToUI9jPXjwoPrXv/4Vf/zJJ59U9913n1JKqXPOOUfddNNNav78+eruu+9WoVCoX2PtKN7++PkfT7zdIVsWfejEE0/k9NNPB2Dv3r28+uqrfOELX2DGjBn87Gc/44UXXmDjxo28+OKLlJWVkZeXF39uXl4epaWlVFdX4/f7MU2z1f1Aq+eYponf76eqquqY462rq+Pzn/88jz32GE8//TTPPfccJSUlreLKz8+ntLQ0IeIF5+N/KBRi7ty5VFRUJNx7e9999zFt2rT47bax9MX72dFrHTlypMexjhw5knPOOQeAqqoqli9fzoUXXkggEOCUU05h8eLFrFy5krq6Oh5//PF+jbWjePvj53888XaHJOR+8Omnn3L99ddz2223ccIJJ/DYY4+Rn59PSkoK1157LevXr8e2bTTtszaASik0TYv/2VLb2y2fo+vH/iOdOnUqy5YtIy0tjezsbK688kp++ctfdhhXIsQL8Nxzz/GNb3wDgNGjRyfse9uss1h6M8a2zzne2EtLS/na177Gl770JaZPn05qaiq/+93vmDBhAqZpcv3117N+/foBj7U/fv69/d62JQm5j73//vt8/etf5z//8z+54oor2LFjB6+99lr8caUUpmkybNiwVieOKioqyM/PJzs7m/r6eizLAqC8vJz8/HzAWV1VVFQAEIvFCAQCZGZmHnOsGzdu5J133mkV28iRI1vF1fz1EyHeSCTChg0b+OIXvwiQ0O9ts7ax9MX7WVBQQFlZWbvXOhbFxcVcffXVXHHFFXz3u98FnBNhL774YvyY5vd5oGPtj59/b8bbEUnIfejw4cN897vf5YEHHmDevHmA85fkZz/7GbW1tUSjUZ5//nlmzZrFyJEj8Xg8vP/++wCsXr2amTNn4nK5mDZtGmvXrgVg1apVzJw5E4DzzjuPVatWAbB27VqmTZuGy+U65njr6+tZtmwZ4XCYhoYGVq5cyQ9+8AP27NnDvn37sCyLV155hZkzZyZEvDt27GDcuHH4fD4gsd/bZlOmTOnz9/O8885j9erVgPNL1uPxMGLEiB7H2tDQwA033MAtt9zC9ddfH7/f6/Xyi1/8ggMHDqCUYvny5cyaNWtAY4X++fn3ZrydfROij9x7773q9NNPV5deemn8fytWrFDPPPOMmjt3rpo1a5b6xS9+ET9++/bt6ktf+pKaPXu2+sEPfqDC4bBSSqmDBw+qr371q2ru3Lnq+uuvVzU1NUoppaqrq9VNN92kLrnkEvWVr3wlfqb5eDz88MNqzpw56uKLL1ZPP/20Ukqpt99+Wy1YsEBdfPHF6r777lO2bSdEvH/+85/V97///Vb3Jep723ziSam+fz9DoZC67bbb1CWXXKIuv/xytXXr1mOK9fe//72aNGlSq7+/jzzyiFJKqXXr1ql58+apiy++WC1ZsmTAYm373vb1z7834u2KTAwRQogEIVsWQgiRICQhCyFEgpCELIQQCUISshBCJAhJyEIIkSAkIQtxnKZOnRrvNibE8ZCELIQQCcIc6ACEGEjvvvsu9957L6+88kqr23/84x+566672LRpE2lpaRQWFgLwX//1X2zcuJF7770XTdM47bTTsG0bgEAgwNKlS9m3bx+6rjNp0iTuueeeXu11IAY3+ZsiRAcef/xxLMvi1Vdf5emnn2bbtm2A0z/jlltuYcmSJaxatYrp06cTCoUAeOONNwgEAqxevTre6+HAgQMD9j2I5CMJWYgOrF+/niuvvBJd1/H7/VxxxRUA7Ny5E9M0+fznPw/A/PnzSU1NBeDMM89k165dXHvttfz2t7/la1/7GmPHjh2w70EkH0nIYkhrbsPYLBqNAk4P3Jb3t9x2aNttoLmf7ujRo3njjTf41re+RUNDA9/4xjf461//2pfhi0FGErIY0rKzsykpKaGyshKlFH/+858Bp9vXSy+9hG3bNDY28sorr6BpGieffDJKKdavXw/Am2++SW1tLQArVqxg6dKlnHvuuSxevJhzzz03vtUhRHdIQhZDWmFhIVdffTVf+tKX+PKXv8yoUaMAuOmmm/B4PCxYsIBvfOMb5OTk4PV6cblcPPbYYzz66KPxIZo5OTkAXH755ViWxSWXXMLChQupr6/n2muvHchvTyQZ6fYmRAf+/Oc/4/f7Oe+887Btm//4j//gnHPO4Zprrhno0MQgJglZiA7s3LmTO++8k8bGRqLRKNOnT+f222/vlSb1QnRGErIQQiQI2UMWQogEIQlZCCEShCRkIYRIEJKQhRAiQUhCFkKIBCEJWQghEsT/DxOm1pc54ZLWAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.lmplot(data = ug1, x = 'ugds', y = 'ug25abv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "The regression line suggests an inverse relationship; however the clustering of dots suggests (potentially) a non-linear relationship, and there are several outliers.  For example:  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>instnm</th>\n",
       "      <th>state</th>\n",
       "      <th>preddeg</th>\n",
       "      <th>control</th>\n",
       "      <th>ugds</th>\n",
       "      <th>ug25abv</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>4880</th>\n",
       "      <td>University of Phoenix-Online Campus</td>\n",
       "      <td>AZ</td>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "      <td>166816.0</td>\n",
       "      <td>0.8368</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                   instnm state  preddeg  control      ugds  \\\n",
       "4880  University of Phoenix-Online Campus    AZ        3        3  166816.0   \n",
       "\n",
       "      ug25abv  \n",
       "4880   0.8368  "
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ug1.query('ugds > 150000')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "The University of Phoenix-Online, with over 160,000 students, stands out as a distinct outlier.  Additionally, a scatter plot is not the only way that we can assess this relationship.  We can convert our data from quantitative to categorical, and in turn assess variations amongst the categories.  \n",
    "\n",
    "### Creating new columns\n",
    "\n",
    "Creating a categorical column from a quantitative column requires computing a new column in our data frame, which will be a common part of your workflow in __pandas__.  As you learned in class, you can use basic mathematical operations to create new columns from existing ones.  In this example, we are going to organize our universities into __bins__ based on their size, and then compare universities across those bins.  \n",
    "\n",
    "__pandas__ includes a lot of methods for creating new columns; we're going to discuss here options for organizing our data into bins.  The three options are as follows: \n",
    "\n",
    "* Equal interval: all bins have the same width, regardless of the number of observations in each bin.  \n",
    "* Manual breaks: the analyst specifies where the bin breaks should be located\n",
    "* Quantile: An equal number of observations are organized into a specified number of bins; bins widths will in turn be irregular.  \n",
    "\n",
    "Equal interval and manual breaks are available via the `cut()` function in __pandas__; quantiles are available from the `qcut()` function.  Let's create a new column with the `.assign()` method that organizes the values in `ugds` into five quantiles, labeled with 1 through 5 which we accomplish with `range()`: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>instnm</th>\n",
       "      <th>state</th>\n",
       "      <th>preddeg</th>\n",
       "      <th>control</th>\n",
       "      <th>ugds</th>\n",
       "      <th>ug25abv</th>\n",
       "      <th>quant5</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Alabama A &amp; M University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4051.0</td>\n",
       "      <td>0.1049</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>University of Alabama at Birmingham</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>11200.0</td>\n",
       "      <td>0.2422</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Amridge University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>322.0</td>\n",
       "      <td>0.8540</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>University of Alabama in Huntsville</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5525.0</td>\n",
       "      <td>0.2640</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Alabama State University</td>\n",
       "      <td>AL</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5354.0</td>\n",
       "      <td>0.1270</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                instnm state  preddeg  control     ugds  \\\n",
       "0             Alabama A & M University    AL        3        1   4051.0   \n",
       "1  University of Alabama at Birmingham    AL        3        1  11200.0   \n",
       "2                   Amridge University    AL        3        2    322.0   \n",
       "3  University of Alabama in Huntsville    AL        3        1   5525.0   \n",
       "4             Alabama State University    AL        3        1   5354.0   \n",
       "\n",
       "   ug25abv quant5  \n",
       "0   0.1049      4  \n",
       "1   0.2422      5  \n",
       "2   0.8540      1  \n",
       "3   0.2640      4  \n",
       "4   0.1270      4  "
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ug_quantiles = ug1.assign(quant5 = pd.qcut(ug.ugds, 5, labels = range(1, 6)))\n",
    "\n",
    "ug_quantiles.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "We can now make comparisons by quantile if we want.  Let's try a bar chart; `seaborn`'s `barplot` function will calculate group means and plot them, and give us an indication of the uncertainty around that mean with an error bar: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='ug25abv', ylabel='quant5'>"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEJCAYAAACaFuz/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAXGUlEQVR4nO3df3BU5b3H8c+GUAgSfxA3YNWKIgYcLUoVjWBJ2poASSDkBwJOMwNaQCgMME2aAgKlQ0XpIDJe5xZsZaigQwsBxAKhFXG4Iiog/uTHRZD8AtYEkgAhWTbn/kG7vbSS3WT37G6efb9mnGFz9pzn+w2ZD49Pzj7HYVmWJQCAsWLCXQAAwF4EPQAYjqAHAMMR9ABgOIIeAAxH0AOA4Qh6ADBcbLgLuJozZ86rudn8W/wTErqquvpcuMsICXo1E71GhpgYh2644ZpvPRaxQd/cbEVF0EuKmj4lejUVvUY2lm4AwHAOtkAAgMBcbHSrvu5iWGuIiXEoIaHrtx6L2KWbac9u0Ddnzoe7DADwac3zT6he4Q36lrB0AwCGI+gBwHAEPQAYjqAHAMMR9ABgOIIeAAxH0AOA4Qh6ADAcQQ8AhiPoAcBwBD0AGI6gBwDD2Rr0586dU2ZmpsrLy+0cBgDQAtuC/sCBAxozZoyOHz9u1xAAAD/YFvRr167VvHnzlJiYaNcQAAA/2LYf/cKFC+26NAAERWNtmc5V7pPlcQd0nYKCbfJ4An+GU1xcZ+Xnj1X//g8EfK3/L2IfPAIAdjt/8lNdulAd8HUqKuqCUM1lmzdvIOgBIFiu6XGvzlW6A57R97gxPmgz+szM7ICv8+8IegBRq9N1t6rTdbcGfJ1Vzz8hl6s+CBXZg/voAcBwts/o3377bbuHAAC0gBk9ABiOoAcAwxH0AGA4gh4ADEfQA4DhCHoAMBxBDwCGI+gBwHAEPQAYjqAHAMMR9ABgOIIeAAznsCwr8E2UASCKXWx0q77uYlhriIlxKCGh67cei9j96Kurz6m52fx/g5zO+IjexzqY6NVM9Br5WLoBAMMR9ABgOIIeAAxH0AOA4Qh6ADAcQQ8AhiPoAcBwEXsf/dVu/DeR0xkf7hJChl4vu9TUqDO1TSGsBtEsYoP+0//+pZrqqsNdBmCLHxS9IomgR2iwdAMAhiPoAcBwBD0AGI6gBwDDEfQAYDiCHgAMR9ADgOEIegAwHEEPAIYj6AHAcAQ9ABjO1r1uXnrpJW3ZskWSNHjwYBUVFdk5HADgW9g2o3/vvfe0a9culZSUaMOGDfr888+1fft2u4YDAFyFbTN6p9Op4uJifec735Ek9erVS5WVlXYNBwC4CtuCvnfv3t4/Hz9+XFu2bNHrr79u13BASH3puqDS/z2jxkvNbTq/02cF8nisNp0bF9dZ+flj1b//A206H9HH9v3ojxw5ookTJ6qoqEg9e/a0ezggJHYeq1VFXQD7yV+oCGj8zZs3EPTwm61Bv3fvXk2bNk2zZs1SRkaGnUMBITX49uvU6Glu+4z+hu4BzegzM7PbdC6ik21BX1VVpSlTpuiFF15QcnKyXcMAYdHX2UV9nV3afP4Pil6Ry1UfxIqAq7Mt6P/whz+osbFRixYt8n5t9OjRGjNmjF1DAgC+hW1BP2fOHM2ZM8euywMA/MQnYwHAcAQ9ABiOoAcAwxH0AGA4gh4ADEfQA4DhCHoAMBxBDwCGI+gBwHAEPQAYjqAHAMMR9ABgONsfPNJW9056LtwlALa51NQY7hIQRSI26Kurz6m5uW0PZmhPnM74qNmXnF6B8GDpBgAMR9ADgOEIegAwHEEPAIYj6AHAcAQ9ABiOoAcAw0XsffQJCV3DXULIOJ3x4S4hZHz12tjUpLpaPkwEBFOrgr6iokKxsbHq3r27XfV4/eLPv9Y352psHweRZeW4FyUR9EAwtbh0M3PmTElSVVWVRo4cqeHDh2vo0KEaO3asTp48GZICAQCBaTHojx07Jkl69tlnlZmZqb1792rfvn3KyMjQ7NmzQ1IgACAwfv0y9tixY3ryySe9r5944glm9ADQTrQY9GfPntXHH3+s2267TUePHvV+vaKiQjEx3LADAO1Bi2mdnp6uRYsW6d1339Xzzz8vSdq4caNycnI0YcKEkBQIAAhMi3fdFBcXS5Kampp0+vRpSVKfPn302muvqXfv3vZXBwAImF/rL/Pnz9ctt9wiSUpKSlLv3r01bdo0WwsDAARHizP6efPm6dSpU9q7d69qav51T/ulS5dUVlZme3EAgMC1GPR5eXk6cuSIDh06pPT0dO/XO3TooPvuu8/u2gAAQdBi0N97772699579cgjj6hHjx6hqgkAEER+bYFQVVWlwsJC1dbWyrL+9RzXN99807bCAADB4VfQz507Vzk5Obr77rvlcDjsrgkAEER+BX1sbKzGjRvX6ou/+OKL2rZtmxwOh/Ly8tp0DQBAYPy6vbJ37946dOhQqy78wQcf6P3339emTZu0bt06/elPf9JXX33VpiIBAG3n14y+rKxMubm5+u53v6tOnTp5v97SGv2AAQO0atUqxcbG6tSpU/J4POrSpUvgFQMAWsWvoJ8xY0abLt6xY0ctW7ZMf/zjHzVkyJCQ7GOPyNJQXq/aT07Jcjf79f6CHQXyeCzfb/yHuLjOys8fq/79H2hriYDx/Ar6AQMG6OzZs2poaJBlWfJ4PDpx4oRfA0ybNk0/+9nPNGnSJK1du1aPP/54QAWjfan/wiV3zUW/319RX9HqMTZv3kDQAy3wK+hffPFFLV++XNLlD0u53W7deeedLS7dHD16VE1NTerbt6/i4uKUlpbW6nV+tH/xdzvVfMn/GX33a52tntFnZma3sTogOvgV9Bs3btSOHTu0aNEiFRUV6f3339fOnTtbPKe8vFzLli3T66+/Lkn6+9//rtzc3MArRrsSd0u84m7x/5m4K8e9KJer3saKgOjjV9B369ZNiYmJuuOOO3Tw4EFlZ2drxYoVLZ4zePBgffLJJ8rOzlaHDh2UlpamjIyMoBQNAPCf3/fRnzhxQnfccYc++ugjDRo0SI2Nvh/gPHXqVE2dOjXgIgEAbefXffQTJ07UM888o5SUFJWWliolJUUPPfSQ3bUBAILArxl9amqqUlNTJV1er//666/Vp08fWwsDAASHX0H/6quv/sfXdu/ezZYGANAO+BX0hw8f9v65qalJH374oZKTk20rCgAQPH4F/bPPPnvF61OnTmn27Nm2FAQACC6/fhn777p3766KitZ/ghEAEHqtXqO3LEuffvqpEhISbCsKABA8rV6jdzgcuvnmm1VcXGxbUQCA4PEr6MvLy+VwOLyPESwvL1dhYaH3+KpVq+ypDgAQML+C/p577tHRo0c1atQodezYURs3btSlS5fY0gAA2gG/gn7fvn1as2aNOnToIEl69NFHNWrUKKWnp9taHAAgcH7ddVNTU3PF3jbnz5/XxYv+7zEOAAgfv2b0mZmZevzxx/XYY4/Jsixt2bJFBQUFthb2u/x5tl4fkamxqSncJQDGcVj//A2rD6WlpdqzZ486deqklJQUDRgwwNbCqqvPqbnZ/wdQtFdOZ3zU7L9Or2ai18gQE+NQQkLXbz3m14xektLS0pSWlha0ogAAodGmT8YCANoPgh4ADEfQA4DhCHoAMBxBDwCGI+gBwHB+314Zale7H9RETmd8uEsIGV+9ui826mw9H5oCgilig37HzF+o4ZvqcJeBEBu26lWJoAeCiqUbADAcQQ8AhiPoAcBwBD0AGI6gBwDDEfQAYDiCHgAMR9ADgOEIegAwHEEPAIYj6AHAcLYH/XPPPafi4mK7hwEAXIWtQb97926VlJTYOQQAwAfbgv7s2bN64YUXNGnSJLuGAAD4wbZtiufOnasZM2aoqqrKriEQ4b662KD/qTsrd7Pl9zlvFBTI4/H//XFxnZWfP1b9+z/QlhKBqGBL0P/5z3/WTTfdpOTkZK1fv96OIdAOfFhfp9Nud6vOOVNR0epxNm/eQNADLbAl6P/617/K5XJpxIgRqq2t1YULF/Tb3/5Ws2bNsmM4RKgH469VUytn9F16dG/1jD4zM7sN1QHRw5agf/XVV71/Xr9+vT744ANCPgrd0TlOd3SOa9U5w1a9Kper3qaKgOjEffQAYDjbnxmbk5OjnJwcu4cBAFwFM3oAMBxBDwCGI+gBwHAEPQAYjqAHAMMR9ABgOIIeAAxH0AOA4Qh6ADAcQQ8AhiPoAcBwBD0AGM72Tc3aKnXJ78JdAsLAfbEx3CUAxonYoK+uPqfmVjywor1yOuOjZv/1aOoViCQs3QCA4Qh6ADAcQQ8AhiPoAcBwBD0AGI6gBwDDEfQAYLiIvY8+IaFruEsIGaczPtwlhIyvXpsaL6m2riFE1QDRIWKD/r8W/1W1Zy+EuwyE2KyFeeEuATAOSzcAYDiCHgAMR9ADgOEIegAwHEEPAIYj6AHAcAQ9ABiOoAcAwxH0AGA4gh4ADEfQA4DhCHoAMJytm5r99Kc/VU1NjWJjLw+zYMEC9evXz84hAQD/xragtyxLx48f144dO7xBDwAIPduWbr766itJ0vjx4zV8+HC99tprdg0FAGiBbVPturo6JScn65lnnpHb7VZBQYFuv/12DRw40K4hEaG+OfO1vir/QJc8bp/vLSjYJI/H8uu6cXGdlZ8/Vv37PxBoiYDRbAv6+++/X/fff7/3dV5ennbu3EnQR6GvKz9W/flv/HpvRUVtq669efMGgh7wwbag/+ijj+R2u5WcnCzp8po9a/XR6bbv3idPeZNfM/puCV1bNaPPzMwOsDrAfLYlb319vZYtW6Y33nhDbrdbJSUl+vWvf23XcIhgN95wm2684Ta/3jtrYZ5crnqbKwKii21Bn5qaqgMHDig7O1vNzc0aO3bsFUs5AIDQsHUtZfr06Zo+fbqdQwAAfOCTsQBgOIIeAAxH0AOA4Qh6ADAcQQ8AhiPoAcBwBD0AGI6gBwDDEfQAYDiCHgAMR9ADgOEIegAwXMRuED+lcFi4S0AYNDVeCncJgHEiNuirq8+pudm/B1C0Z05nfNTsvx5NvQKRhKUbADBcxM7oY2Ic4S4hZOjVTPRqpkjttaW6HJZlmb8+AgBRjKUbADAcQQ8AhiPoAcBwBD0AGI6gBwDDEfQAYDiCHgAMR9ADgOEIegAwXFiD/s0339SwYcOUlpam1atX/8fxL7/8Ujk5OUpPT9fs2bN16VL73dnQV6//VFRUpPXr14ewsuDz1evf/vY3jRgxQsOHD9fkyZNVW1sbhiqDw1ev27dvV1ZWljIyMlRcXKympqYwVBkc/v4Mv/POO/rRj34UwsqCz1evL730klJTUzVixAiNGDGixe9HRLDC5OTJk1Zqaqp15swZ6/z581ZWVpZ15MiRK96TkZFh7d+/37Isy/rVr35lrV69OgyVBs6fXk+ePGlNnDjR+v73v2+tW7cuTJUGzlev9fX11sCBA62TJ09almVZS5cutX7zm9+Eq9yA+Or1/Pnz1qBBgyyXy2VZlmVNnz7deuONN8JVbkD8+Rm2LMtyuVzWkCFDrNTU1DBUGRz+9Dpx4kRr3759Yaqw9cI2o3/vvff08MMP6/rrr1eXLl2Unp6urVu3eo9XVFTo4sWLuu+++yRJOTk5VxxvT3z1Kl2eQfz4xz/W0KFDw1RlcPjq1e12a968eerevbskKSkpSVVVVeEqNyC+eu3SpYvefvtt3XjjjWpoaFB1dbWuvfbaMFbcdv78DEvSnDlz9POf/zwMFQaPP71+9tln+v3vf6+srCwtWLBAjY2NYarWP2EL+tOnT8vpdHpfJyYm6tSpU1c97nQ6rzjenvjqVZKeeuop5efnh7q0oPPV6w033KDHHntMknTx4kUtX75cP/nJT0JeZzD48/fasWNH7dy5UykpKTpz5owGDRoU6jKDwp9eV61apbvvvlv9+vULdXlB5avX8+fPq2/fviosLFRJSYnq6ur08ssvh6NUv4Ut6Jubm+Vw/GtbTcuyrnjt63h7YlIvvvjba319vSZMmKA+ffpo5MiRoSwxaPztdfDgwdqzZ49SU1M1f/78EFYYPL56PXz4sEpLSzV58uRwlBdUvnq95pprtGLFCvXq1UuxsbEaP368du7cGY5S/Ra2oO/Ro4dcLpf3tcvlUmJi4lWPf/PNN1ccb0989WoSf3o9ffq0xo4dq6SkJC1cuDDUJQaNr17Pnj2rXbt2eV9nZWXp0KFDIa0xWHz1unXrVrlcLuXm5mrChAnev+P2yFevlZWV+stf/uJ9bVmWYmMj9tEeksIY9I888oh2796tmpoaNTQ0qLS0VD/84Q+9x2+++WZ16tRJe/fulSRt3LjxiuPtia9eTeKrV4/Ho0mTJmno0KGaPXt2u/4/G1+9WpalwsJCVVZWSrochv379w9XuQHx1eu0adO0bds2bdy4UcuXL1diYqLWrFkTxorbzlevnTt31uLFi1VWVibLsrR69WrvcmTECtdvgS3LsjZt2mRlZGRYaWlp1vLlyy3LsqynnnrK+uSTTyzLsqwvv/zSys3NtdLT062ZM2dajY2N4Sw3IL56/adf/vKX7fquG8tqudfS0lIrKSnJGj58uPe/WbNmhbnitvP197p9+3YrMzPTysrKsmbMmGHV1dWFs9yA+PszXFZW1q7vurEs371u3brVe7y4uDjis4knTAGA4fhkLAAYjqAHAMMR9ABgOIIeAAxH0AOA4SL7Ln8ghL7++mvNnz9fNTU1crvdysvL0/jx4yVJixYt0tatW3XddddJkm6//XYtXbq0xeslJSVp9+7d6tatm92lAy0i6IF/KC4uVk5OjvLz81VfX6+8vDz17dtXycnJ2r9/v5YsWdJuP/CE6MbSDYy0Z88eZWZm/sfrhoYGFRUVKT09XXl5eSouLlZxcbEkKS8vz3tOfHy8vve976myslJNTU364osv9MorrygrK0tTp071ftr12LFjGjdunEaNGqXU1FQ9/fTTV+xkuHTpUo0cOVIjRozQjh07JEmjR4/Wtm3bvO9ZvHixFi9ebPv3BNGLoEdUefnll+XxeLRlyxatXLlSX3zxhfdYbm6u4uLiJEnvvvuu9u/fr0cffVSnTp3Sww8/rOnTp2vTpk3q16+fJk+eLMuytHbtWmVnZ2vt2rUqLS1VeXm53nnnHe81b7nlFpWUlGjx4sUqLi5WTU2N8vPzvQ+X8Xg82rRpkxE7lyJyEfSIKjt37lReXp5iYmLUtWvXb905c8OGDSosLNSyZcuUmJioW2+9VStWrNBdd90lh8OhJ598UidOnFB5ebkKCwvVrVs3rVixQvPnz9fp06d14cIF77XGjBkjSbrrrrvUq1cv7d+/X8OGDdPHH38sl8ulXbt2qWfPnurZs2eovgWIQqzRw0gOh0P/f3cPt9stSYqNjb3i6zEx/5rrWJal5557Ttu2bdPKlSvVt29fSdLBgwd18OBBZWdnX/Hejh07aubMmfJ4PBo6dKhSUlJUVVV11es3NzcrNjZWcXFxSk9P1+bNm7V//35m87AdM3oYqVu3bqqsrFR1dbUsy9Jbb70l6fLe8OvWrVNzc7MaGhq0efNm7w6azz//vD788EOtW7fOG/LS5bBeuHChysrKJElr1qxRUlKSevTooV27dmnKlCkaNmyYJOnAgQPyeDzec0tKSiRJn3/+uU6cOOF9KMeoUaNUUlKiffv2KT093f5vCKIaM3oY6c4779To0aOVm5srp9OplJQUffrpp5o4caIWLFigrKwsxcfHKyEhQZ07d9bJkye1cuVK3XTTTRo3bpz3OgUFBcrNzdWcOXP09NNPy+PxqEePHlqyZIkkacaMGZoyZYq6dOmirl276sEHH9SJEye855eVlSk7O1sOh0NLlizR9ddfL0m655571KFDBw0ZMkSdOnUK6fcG0YfdKxFV3nrrLXXt2lWDBw9Wc3Ozpk6dqoEDB7bbh2QA/iDoEVUOHz6suXPnqqGhQW63Ww899JBmzZqljh07hrs0wDYEPQAYjl/GAoDhCHoAMBxBDwCGI+gBwHAEPQAYjqAHAMP9H3GtTxSdbg5TAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.barplot(x = 'ug25abv', y = 'quant5', data = ug_quantiles)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can analyze this further by breaking out the bar charts by college type with `sns.catplot()`: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<seaborn.axisgrid.FacetGrid at 0x7f34141fe9a0>"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABDoAAAFqCAYAAAD7mzrpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAp6klEQVR4nO3de5SVBb0//s/AEJekFBzRUx1vKVKaSZlXDtJRUBgEualolJeTpEePer7gqChqmQotU5fVV82jy5Pm4agoYSiZivnT1LxkmreTFnccQS4qzAwzz+8Pv8zRENgzs595nnl4vdZqLfae2c/zZpa8997v9t5TkSRJEgAAAAAF0CnrAAAAAADlYugAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB3w/1x//fXx0EMPtfh2ffv2jRUrVrT5/EuWLIkBAwaU5VgAHU1WHbxu3bo4//zzo7q6OoYNGxbnn39+rFu3rtXHA+iIsurgNWvWxFlnnRXV1dUxdOjQuPHGG1t9LPgoQwf8P0899VSsX78+k3Pfe++9ccIJJ8Tbb7+dyfkBspZVB//sZz+LxsbGmDVrVsyaNSvq6urihhtuaPccAFnKqoOvvfba6NOnT8yePTvuuuuuuPPOO+P5559v9xwUT2XWAaAt7rrrrrjllluiU6dOsd1228VVV10VO+20U/zXf/1X/Od//md06tQptt9++7joooti1113jZqamthmm23itddei6VLl0bfvn3jqquuinvvvTdeeumlmDZtWnTu3Dl++9vfxsqVK2PBggVx2GGHxcSJE+PSSy+NV199NSoqKmLAgAFx7rnnRmXlJ/8TWr16dXzrW9/a6Pojjzwyvve9733sumXLlsVDDz0UN998cxx55JGp/JwA0lCEDt5///3jc5/7XHTq9OH/99OvX7/4n//5n/L/sADKrAgdfOGFF0ZjY2NERNTW1kZ9fX307Nmz/D8stj4JdFCvvPJKcsABBySLFy9OkiRJbrnlluSiiy5KnnjiieTwww9Pli9fniRJktx9993JUUcdlTQ1NSXnnXdecuyxxyZ1dXVJfX19MnLkyOSuu+5KkiRJTjzxxGTOnDlJkiTJeeedl3z7299uPtfkyZOT73//+0lTU1NSV1eXnHzyyckNN9yQJEmS7Lnnns3naqtyHgsgTUXs4IULFyaHHHJI8vDDD5fleABpKVoH//u//3uy9957J+ecc06yfv36Nh8PvHWFDuvJJ5+MQw89NHbaaaeIiPjOd74Tl112Wfzud7+LoUOHRq9evSIiYtSoUbFs2bJYuHBhREQMGDAgPvWpT0WXLl1izz33jFWrVn3i8b/2ta81//mxxx6LE088MSoqKuJTn/pUHHfccfHYY49tMtvq1atjxIgRG/3vZz/7Wbn++gCZKloHv/TSS3HCCSfEiSeeGIMGDWrxzwOgPRWtg3/0ox/F73//+1i1alX85Cc/afHPA/6et67QYXXu3DkqKiqaL69bty4WLVoUTU1NG31vkiTN7zvs1q1b8/UVFRWRJMknHr9Hjx7Nf25qavrYuZqamjb7PsbPfOYzcd9995X+lwHoYIrUwffff39ceumlcdFFF8Xw4cNLvh1AVorSwb/73e9izz33jD59+sSnP/3pGDZsWMydO7ek28LmeEUHHdYBBxwQTz75ZPMHeN55550xffr0GDBgQPz6179u/gTou+++O7bddtvYeeedN3u8zp07b7K0Dz300PjFL34RSZJEfX19zJgxIw4++ODy/oUAOpCidPDDDz8cP/jBD+Lmm282cgAdRlE6eM6cOfGTn/yk+dhz5syJAw88sCzHZuvmFR10WH379o1JkybFqaeeGhERVVVV8cMf/jD69OkT3/nOd+Lb3/52NDU1Ra9eveKGG25o/qC5TfnmN78ZV199dTQ0NGz0tSlTpsQPfvCDGD58eDQ0NMSAAQNi4sSJqfy9ADqConTwVVddFUmSxJQpU5qv69+/f0ydOrUsxwdIQ1E6uKamJqZOndo8NB9++OExYcKEshybrVtFsqnXKwEAAAB0MN66AgAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYeT218suX/5eNDXl5xfCbLddj3j33Q+yjrGRPOaSqXR5zJXHTBH5zJWHTFVVPVM5rg4uTR5z5TFTRD5zyVS6PObKQyYdnK085spjpoh85pKpdHnMlYdMm+tgr+goUWVl56wjfKI85pKpdHnMlcdMEfnMlcdMRZXXn3Uec+UxU0Q+c8lUujzmymOmosrrzzqPufKYKSKfuWQqXR5z5THTRxk6AAAAgMIwdAAAAACFYegAAAAACsPQAQAAABRGbn/rSu/e22QdYSNpfbJ2W+Uxl0xbtq6uIesIsEk6uHR5zJXHTBH/m2tdXUOsWb0u4zSQXzq4dHnMlcdMEfnMJVNp6hsas47Q4eR26DjrinvjnXffzzoGpOaOaSdkHQE2SQeTpjumnRBrwtABm6KDgY/yvKHlvHUFAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYqQ4d7733XlRXV8fChQvTPA0AAABARKQ4dPzxj3+M448/Pv7617+mdQoAAACAj0lt6JgxY0ZMnTo1dthhh7ROAQAAAPAxlWkd+PLLL0/r0FCyulUL4r3Fz0XS2JB1lI1MmPBgREQ0NiYZJ/lf3bt3i1NPPSV23/3LWUcBOoDWduyECQ+2uvu6d+8WY8eOj/79v96q2wNAFtryvKS195tb831makMH5MH7S/8U6z9YnnWMT7Ro0eqsI3yiGTNmxPnnX5p1DKADaG3HtrX/Zs++d6t80AZAx9WW5yVtud/cWu8zDR0U2qd33CfeW9yQy1d07Lh9z4jI3ys6xo0bl3UMoINobcfuuH3PNr2io7p6ZKtuCwBZacvzktbeb27N95mGDgqt62e/EF0/+4WsY3yi26adEBERtbVrMk7ycVVVPXOXCcin1nbsbdNO0DMAbFXa8rzE/WbLpfrrZQEAAADaU+qv6Hj44YfTPgUAAABARHhFBwAAAFAghg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGFUZh1gU647f2TWESBV6+oaolvXLlnHgE+kg0nTurqGrCNArulg4KPqGxqzjtDh5HboWL78vWhqSrKO0ayqqmfU1q7JOsZG8phLptJ1qzJ0kE86uDR5zJXHTBH5zQV5pINLk8dcecwUkc9cMpWuqqpn1hE6HG9dAQAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAURmXWATald+9tso6wkaqqnllH+ER5y9W0viHrCEAb6eDSbSnX+vq6eHdVfTulAYpAB5cuj7nymCkin7mKlsl9fn7kduj40/89L+pXL886Bq3wtck/j4h1WccA2kAHl8+HnehBD1A6HQwdk/v8/PDWFQAAAKAwDB0AAABAYRg6AAAAgMIwdAAAAACFYegAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYRg6AAAAgMIwdAAAAACFYegAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYVSmefDrr78+5syZExERAwcOjMmTJ6d5OgAAAGArl9orOp544ol4/PHHY+bMmXHvvffGyy+/HL/5zW/SOh0AAABAeq/oqKqqipqamvjUpz4VERG77757LF68OK3TAQAAAKQ3dOyxxx7Nf/7rX/8ac+bMiV/+8pdpnY4WeqX2g5j7P+9G3fqmsh+760sTorExKftxN+jevVuMHTs++vf/emrnANiSUnu0lE7UawDQvtJ4PlTO50EeG7RNqp/RERHxxhtvxGmnnRaTJ0+OXXbZJe3TUaJ5b62KRavr0zn4B4vSOe5HzJ59r3/0QKZK7tESO1GvAUD7SeX5UJmfB3ls0HqpDh3PPvtsnHXWWXHBBRfEsGHD0jwVLTRw189GXWNTOq/o2K5P6q/oqK4emdrxAUpRao+W0ol6DQDaVxrPh8r5PMhjg7ZJbehYsmRJnHHGGfHjH/84DjrooLROQyv1q+oR/ap6pHLsr03+edTWrknl2AB5UWqP6kQAyJ80ng+5z8+P1IaOm2++Oerq6uLKK69svu64446L448/Pq1TAgAAAFu51IaOKVOmxJQpU9I6PAAAAMBGOmUdAAAAAKBcDB0AAABAYRg6AAAAgMIwdAAAAACFYegAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYRg6AAAAgMIwdAAAAACFYegAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYRg6AAAAgMKozDrApuwz8aqsI9BKTesbso4AtJEOLp/19XVZRwA6GB0MHZP7/PzI7dCxfPl70dSUZB2jWVVVz6itXZN1jI3kMVdVVc+IWJd1DKANdHBp8poL6Nh0cGnymCuPmSLymUsm0uStKwAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwqjMOsCm9O69TdYRNlJV1TPrCJ8ozVx19fWxelVdascH8mlr62BdB+TJ1tbBbZHHXHnMFJHPXHnP5PFBx5XboeP//Pel8c57K7KOsdW79aRrI8I/btjabG0drOuAPNnaOhjyyuODjstbVwAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYLRo6Fi1aFMuWLUsrCwAAAECbbHboOPfccyMiYsmSJXHMMcfE0UcfHUcddVSMHz8+li5d2i4BAQAAAEq12aHjrbfeioiIK664Iqqrq+PZZ5+N5557LoYNGxYXXnhhuwQEAAAAKFVJb11566234pRTTmm+fMIJJ3hFBwAAAJA7mx06Vq5cGS+88ELsvPPO8Ze//KX5+kWLFkWnTj7HFAAAAMiXza4VQ4YMiSuvvDIee+yxmDZtWkRE3HfffTFq1Kj47ne/2y4BAQAAAEpVubkv1tTUREREfX19vP322xERsddee8UvfvGL2GOPPdJPBwAAANACJb3/5JJLLonPf/7zERHRt2/f2GOPPeKss85KNRgAAABAS232FR1Tp06NZcuWxbPPPhsrVqxovn79+vWxYMGC1MMBAAAAtMRmh44xY8bEG2+8Ea+99loMGTKk+frOnTvHV7/61bSzAQAAALTIZoeOffbZJ/bZZ584+OCDY8cdd2yvTAAAAACtstmhY4MlS5bEpEmTYtWqVZEkSfP1v/rVr1ILBgAAANBSJQ0dF198cYwaNSq+9KUvRUVFRckHv/baa+PBBx+MioqKGDNmTJx00kmtDgoAAACwJSUNHZWVlS0eKZ5++un4/e9/H7NmzYr169fH0KFDY+DAgbHbbru1KigAAADAlpT062X32GOPeO2111p04G984xtx2223RWVlZSxfvjwaGxujR48erQoJAAAAUIqSXtGxYMGCGD16dPzDP/xDdO3atfn6LX1GR5cuXeK6666L//iP/4gjjzwy+vTp07a0tNjahWti1YvLImloatXtJzwyIRobky1/40d07lzxibfp3r1bjB07Pvr3/3qrsgC01pa6cFNdp7cAoP219TlMubTmuVA5ePzRdiUNHeecc06rT3DWWWfFv/zLv8TEiRNjxowZceyxx7b6WLTcmj/XRsOKda2+/aI1i8qYJmL27Hv9gwXa3Za6cHNdp7cAoH219TlMuZT7uVBLePzRNiUNHd/4xjdi5cqVsXbt2kiSJBobG2P+/Pmbvc1f/vKXqK+vj379+kX37t1j8ODBLX77C23X80tV0bS+9Wton89UlfUVHdXVI1uVA6AtttSFm+o6vQUA7a+tz2HKpTXPhcrB44+2K2nouPbaa+PGG2+MiIjOnTtHQ0NDfPGLX9zsW1cWLlwY1113Xfzyl7+MiIjf/va3MXr06DJEpiW6f75ndP98z1bf/taTro3a2jUtuk1VVc8W3wYgTVvqwtZ0HQCQjrY+hykXjw86rpKGjvvuuy8eeeSRuPLKK2Py5Mnx+9//PubNm7fZ2wwcODBefPHFGDlyZHTu3DkGDx4cw4YNK0toAAAAgE9S0tDRq1ev2GGHHWK33XaLV199NUaOHBk33XTTFm935plnxplnntnmkAAAAAClKOnXy1ZWVsb8+fNjt912iz/84Q+xfv36qKurSzsbAAAAQIuUNHScdtppcdFFF8Vhhx0Wc+fOjcMOOywOOOCAtLMBAAAAtEhJb10ZNGhQDBo0KCI+/LyOv/3tb7HXXnulGgwAAACgpUoaOm655ZaNrnvyySfjpJNOKnsgAAAAgNYqaeh4/fXXm/9cX18fzzzzTBx00EGphQIAAABojZKGjiuuuOJjl5ctWxYXXnhhKoEAAAAAWqukDyP9e3369IlFixaVOwsAAABAm7T4MzqSJIk//elP0bt379RCAQAAALRGiz+jo6KiIj73uc9FTU1NaqEAAAAAWqOkoWPhwoVRUVERSZI0X540aVLz12+77bZ00gEAAAC0QElDx9577x1/+ctfYty4cdGlS5e47777Yv369TFs2LC08wEAAACUrKSh47nnnos77rgjOnfuHBERAwYMiHHjxsWQIUNSDQcAAADQEiX91pUVK1ZEXV1d8+X3338/1q1bl1ooAAAAgNYo6RUd1dXVceyxx8YRRxwRSZLEnDlzYsKECWlnAwAAAGiRkoaOf/u3f4t+/frFU089FV27do3LLrssvvGNb6SdDQAAAKBFSho6IiIGDx4cgwcPTjMLAAAAQJuU9BkdAAAAAB2BoQMAAAAoDEMHAAAAUBglf0ZHe/vR2KlZRyAi6urrs44AZGBr62BdB+TJ1tbBkFceH3RcuR06li9/L5qakqxjNKuq6hm1tWuyjrGRvOYCOjYdDJAdHVyaPObKY6aIfOaSiTR56woAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAqsw6wKb17b5N1hI1UVfXMOsIn2lKuhnV1sXJNfTulAYpAB+tOIDs6uHR5zJXHTBH5zJWXTO7ziye3Q8cj5/6fWPvO8qxjFMLQ226J8A8XaAEdrDuB7OhgaF/u84vHW1cAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKI/Wh46qrroqampq0TwMAAACQ7tDx5JNPxsyZM9M8BQAAAECzyrQOvHLlyvjxj38cEydOjFdffTWt02zV3ly3Nv6/1SujoSnZ7PfdOWFCNDZu/nsiIrp37xZjx46P/v2/Xq6IALlUSn9uqTu7d+8Wp556Suy++5fTiAgAHVKpz1HyZMN9fufOFSU9b2pvG3J5vla61IaOiy++OM4555xYsmRJWqfY6j2zZnW83dCwxe97d9Giko85e/a9/uEAhVdKf5bSnTNmzIjzz7+0XLEAoMMr9TlKnrTk+VLWPF8rTSpDx3//93/HTjvtFAcddFDcc889aZyCiNi/52eivoS1tMeOfUp+RUd19cgypQPIr1L6c0vd2b17txg3blwa8QCgwyr1OUqebLjP7wiv6PB8rTSpDB2//vWvo7a2NkaMGBGrVq2KDz74IH74wx/GBRdckMbptlq7deseu3XrvsXvG3rbLVFbu6YdEgF0DKX0ZyndWVXVU78CwEeU+hwlTzbc5+f1fj2vufIslaHjlltuaf7zPffcE08//bSRAwAAAEhd6r9eFgAAAKC9pPZhpBuMGjUqRo0alfZpAAAAALyiAwAAACgOQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGJVZB9iUQVf/KOsIhdGwri7rCEAHo4N1J5AdHQzty31+8eR26Fi+/L1oakqyjtGsqqpn1NauyTrGRvKaC+jYdDBAdnRwafKYK4+ZIvKZK4+ZKA5vXQEAAAAKw9ABAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFEZl1gE2pXfvbbKOsJGqqp6Znbu+bn2sWr02s/MDWxcdXLqP5tLVQDno4NLlMVceM0XkM5dMpcs6V0d7jJPboeMn038dq1Z+kHWM3Ljg8jFZRwC2Ijq4dXQ1UA46GMibjvYYx1tXAAAAgMIwdAAAAACFYegAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYRg6AAAAgMIwdAAAAACFYegAAAAACsPQAQAAABSGoQMAAAAoDEMHAAAAUBiGDgAAAKAwDB0AAABAYRg6AAAAgMIwdAAAAACFYegAAAAACqMyzYN/61vfihUrVkRl5Yenueyyy2LfffdN85QAAADAViy1oSNJkvjrX/8ajzzySPPQAQAAAJCm1BaIN998MyIiTj755Fi5cmWMGzcuTjzxxLROVxjvvPu3eHPh07G+seFj10+YMCsaG5OPXde9e7c49dRTYvfdv9yeEQE6vE11bVt9UleXQ/fu3WLs2PHRv//Xy35sAKDjSeuxzKb8/WOczp0r2vyYJ83HN6kNHatXr46DDjooLrroomhoaIgJEybErrvuGoccckhapyyEvy1+Ida8/85G1y9atOoTv3/GjBlx/vmXph0LoFA21bVttamuLofZs+81dAAAEZHeY5lNSesxTlqPb1IbOvbbb7/Yb7/9mi+PGTMm5s2bZ+jYgp3/4avRuLB+o2WuV+9tPvEVHePGjWvPeACFsKmubatP6upy6N69W1RXjyz7cQGAjimtxzKb8vePccr1io60Ht+kNnT84Q9/iIaGhjjooIMi4sPP7PBZHVu2/XY7x/bb7bzR9RdcPiZqa9dsdH1VVc9PvB6ATdtU17bVproaAKCc0nossyl//xgn789DU/v1smvWrIlp06ZFXV1dvPfeezFz5sw44ogj0jodAAAAQHqv6Bg0aFD88Y9/jJEjR0ZTU1OMHz/+Y29lAQAAACi3VN9LcvbZZ8fZZ5+d5ikAAAAAmqX21hUAAACA9mboAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAArD0AEAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDCMHQAAAAAhWHoAAAAAAqjMusAm3LGpKFZR8iV+rr1WUcAtiI6uHV0NVAOOhjIm472GCe3Q8fy5e9FU1OSdYxmVVU9o7Z2TdYxANqFDi5NXnMBHZsOLk0ec+UxU0Q+c8lUurzmyjNvXQEAAAAKw9ABAAAAFIahAwAAACiM3H5GR6dOFVlH2EgeM0XkM5dMpctjrjxmishnrjxmKoc8/r3ymCkin7nymCkin7lkKl0ec+UxUznk8e+Vx0wR+cyVx0wR+cwlU+nymCuPmTaoSJIkP590BAAAANAG3roCAAAAFIahAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhGDoAAACAwjB0AAAAAIVh6AAAAAAKw9ABAAAAFEamQ8evfvWrGDp0aAwePDhuv/32jb7+yiuvxKhRo2LIkCFx4YUXxvr16zPPtMHkyZPjnnvuST1PKZkeeuihGDFiRBx99NFx+umnx6pVq3KR6ze/+U0MHz48hg0bFjU1NVFfX595pg0effTR+OY3v5l6nlJzXX/99TFo0KAYMWJEjBgxYrPZ2yvTm2++Gd/61rfi6KOPjlNOOSUX/1298sorzT+jESNGxIABA6K6ujrTTBERL7/8cowePTqOPvroOO2002L16tWpZyoHHVy+XFn0sA4uXy4dXFouHVxeOrh8uXRwaZk20ME6uFyZInLcwUlGli5dmgwaNCh59913k/fffz8ZPnx48sYbb3zse4YNG5Y8//zzSZIkyfnnn5/cfvvtmWdaunRpctpppyVf+cpXkrvvvjvVPKVkWrNmTXLIIYckS5cuTZIkSa655prk+9//fua53n///eTQQw9NamtrkyRJkrPPPju58847M820QW1tbXLkkUcmgwYNSjVPS3KddtppyXPPPdcueUrJ1NTUlAwePDiZN29ekiRJMn369GTatGmZ5/qoDz74IBk2bFjyzDPPZJ7p+OOPTx599NEkSZLkiiuuSK6++upUM5WDDi5frix6WAeXN5cOLi3XR+ngttHB5culg0vLtIEO1sHlzpTXDs7sFR1PPPFEHHjggbHttttGjx49YsiQIfHAAw80f33RokWxbt26+OpXvxoREaNGjfrY17PIFPHhovXP//zPcdRRR6WapdRMDQ0NMXXq1OjTp09ERPTt2zeWLFmSea4ePXrEww8/HNtvv32sXbs2li9fHp/5zGcyzbTBlClT4l//9V9TzdLSXC+99FLccMMNMXz48Ljsssuirq4u00wvv/xy9OjRI/7pn/4pIiImTpwYJ5xwQqqZSsn1UTfccEPsv//+8fWvfz3zTE1NTfH+++9HRMTatWujW7duqWYqBx1cvlxZ9LAOLm8uHVxaro/SwW2jg8uXSweXlmkDHayDy50prx2c2dDx9ttvR1VVVfPlHXbYIZYtW7bJr1dVVX3s61lkiog49dRTY+zYsanmaEmm7bbbLo444oiIiFi3bl3ceOONcfjhh2eeKyKiS5cuMW/evDjssMPi3XffjUMPPTTzTLfddlt86Utfin333TfVLC3J9f7770e/fv1i0qRJMXPmzFi9enX89Kc/zTTT/PnzY/vtt48LLrggjjnmmJg6dWr06NEj1Uyl5NpgzZo1MWPGjHa5oy4lU01NTUyZMiUOPfTQeOKJJ+K4445LPVdb6eDy5cqih3Vw+XLp4NJzbaCD204Hly+XDi49kw4uLZMOblmmvHZwZkNHU1NTVFRUNF9OkuRjl7f09SwyZaHUTGvWrInvfve7sddee8UxxxyTm1wDBw6Mp556KgYNGhSXXHJJpplef/31mDt3bpx++ump5mhprk9/+tNx0003xe677x6VlZVx8sknx7x58zLNtH79+nj66afj+OOPj5kzZ8YXvvCFuPLKK1PNVEquDWbNmhWHH3549O7dO/NM69atiwsvvDBuvfXWePzxx2P8+PFx3nnnpZ6rrXRw6fLYwzq4fLl0cOm5NtDBbaeDS6eDy5NJB5eeSQeXninPHZzZ0LHjjjtGbW1t8+Xa2trYYYcdNvn1d95552NfzyJTFkrJ9Pbbb8f48eOjb9++cfnll+ci18qVK+Pxxx9vvjx8+PB47bXXMs30wAMPRG1tbYwePTq++93vNv/c0ralXIsXL4677rqr+XKSJFFZWZlppqqqqth5551jn332iYiI6urqePHFF1PNVEquDR566KEYOnRo6nlKyfT6669H165d4ytf+UpERBx77LHx9NNPt0u2ttDBpctjD+vg8uXSwaXn2kAHt50OLp0OLk8mHVx6Jh1ceqY8d3BmQ8fBBx8cTz75ZKxYsSLWrl0bc+fObX4fVETE5z73uejatWs8++yzERFx3333fezrWWTKwpYyNTY2xsSJE+Ooo46KCy+8sN2W9y3lSpIkJk2aFIsXL46ID8u1f//+mWY666yz4sEHH4z77rsvbrzxxthhhx3ijjvuSDVTKbm6desW06dPjwULFkSSJHH77bc3vwQzq0z77bdfrFixIl599dWIiHj44Yfjy1/+cqqZSskV8eF/Wy+//HLst99+qecpJdPOO+8cS5cujTfffDMiIn7729823zHmmQ4uX64selgHly+XDi49V4QOLhcdXL5cOri0TDq49Ew6uPRMue7gtD7ltBSzZs1Khg0blgwePDi58cYbkyRJklNPPTV58cUXkyRJkldeeSUZPXp0MmTIkOTcc89N6urqMs+0wXnnnddunza9uUxz585N+vbtmxx99NHN/7vgggsyz5UkSfKb3/wmqa6uToYPH56cc845yerVqzPPtMGCBQva7dOmS8n1wAMPNH+9pqYmF/+tv/DCC8no0aOToUOHJieffHLyzjvvpJ6plFzvvPNOcvDBB7dLllIzPfroo8nw4cOT6urq5Nvf/nYyf/78ds3XWjq4PLmy6mEdXL5cOrj0XDq4fHRweXLp4NIzbaCDdXA5M+W1gyuSJEmyHlsAAAAAyiGzt64AAAAAlJuhAwAAACgMQwcAAABQGIYOAAAAoDAMHQAAAEBhVGYdAPLgb3/7W1xyySWxYsWKaGhoiDFjxsTJJ58cERFXXnllPPDAA/HZz342IiJ23XXXuOaaazZ7vL59+8aTTz4ZvXr1Sjs6QIengwGyo4MpIkMHRERNTU2MGjUqxo4dG2vWrIkxY8ZEv3794qCDDornn38+rr766ujfv3/WMQEKSQcDZEcHU0TeukLhPPXUU1FdXb3R5bVr18bkyZNjyJAhMWbMmKipqYmampqIiBgzZkzzbXr27Bn/+I//GIsXL476+vr485//HD//+c9j+PDhceaZZ8bixYsjIuKtt96Kk046KcaNGxeDBg2K733ve1FXV9d83muuuSaOOeaYGDFiRDzyyCMREXHcccfFgw8+2Pw906dPj+nTp6f+MwFoLzoYIDs6GD5k6GCr8dOf/jQaGxtjzpw5ceutt8af//zn5q+NHj06unfvHhERjz32WDz//PMxYMCAWLZsWRx44IFx9tlnx6xZs2LfffeN008/PZIkiRkzZsTIkSNjxowZMXfu3Fi4cGE8+uijzcf8/Oc/HzNnzozp06dHTU1NrFixIsaOHRv33HNPREQ0NjbGrFmzYuzYse36cwDIgg4GyI4OZmtj6GCrMW/evBgzZkx06tQpttlmmzjmmGM2+p577703Jk2aFNddd13ssMMO8YUvfCFuuumm2HPPPaOioiJOOeWUmD9/fixcuDAmTZoUvXr1iptuuikuueSSePvtt+ODDz5oPtbxxx8fERF77rln7L777vH888/H0KFD44UXXoja2tp4/PHHY5dddolddtmlvX4EAJnRwQDZ0cFsbXxGB4VTUVERSZI0X25oaIiIiMrKyo9d36nT/+58SZLEVVddFQ8++GDceuut0a9fv4iIePXVV+PVV1+NkSNHfux7u3TpEueee240NjbGUUcdFYcddlgsWbJkk8dvamqKysrK6N69ewwZMiRmz54dzz//vBUbKBwdDJAdHQwf8ooOCqdXr16xePHiWL58eSRJEvfff39ERAwcODDuvvvuaGpqirVr18bs2bOjoqIiIiKmTZsWzzzzTNx9993N5R7xYUlffvnlsWDBgoiIuOOOO6Jv376x4447xuOPPx5nnHFGDB06NCIi/vjHP0ZjY2PzbWfOnBkRES+//HLMnz8/9t1334iIGDduXMycOTOee+65GDJkSPo/EIB2pIMBsqOD4UNe0UHhfPGLX4zjjjsuRo8eHVVVVXHYYYfFn/70pzjttNPisssui+HDh0fPnj2jd+/e0a1bt1i6dGnceuutsdNOO8VJJ53UfJwJEybE6NGjY8qUKfG9730vGhsbY8cdd4yrr746IiLOOeecOOOMM6JHjx6xzTbbxP777x/z589vvv2CBQti5MiRUVFREVdffXVsu+22ERGx9957R+fOnePII4+Mrl27tuvPBiBtOhggOzoYPlSRfPQ1RlBg999/f2yzzTYxcODAaGpqijPPPDMOOeSQGD9+fNbRAApPBwNkRweztTF0sNV4/fXX4+KLL461a9dGQ0NDHHDAAXHBBRdEly5dso4GUHg6GCA7OpitjaEDAAAAKAwfRgoAAAAUhqEDAAAAKAxDBwAAAFAYhg4AAACgMAwdAAAAQGEYOgAAAIDC+P8BKYDz5tJlyR0AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1080x360 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.catplot(data = ug_quantiles, x = 'ug25abv', y = 'quant5', col = 'control', kind = 'bar')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "Our results appear quite different when we break them out by college type!  This starts to reveal the importance of detailed investigation of your data through visualization; we'll be exploring this much more in the next few weeks.  We can also take a look at the numbers behind these plots by creating a `groupby` object that groups by `control` and the quantiles: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "control  quant5\n",
       "1        1         0.498289\n",
       "         2         0.298886\n",
       "         3         0.265841\n",
       "         4         0.224723\n",
       "         5         0.183518\n",
       "2        1         0.410560\n",
       "         2         0.225184\n",
       "         3         0.162080\n",
       "         4         0.176742\n",
       "         5         0.210676\n",
       "3        1         0.701017\n",
       "         2         0.641752\n",
       "         3         0.610130\n",
       "         4         0.594723\n",
       "         5         0.787532\n",
       "Name: ug25abv, dtype: float64"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "groups2 = ug_quantiles.groupby(['control', 'quant5'])\n",
    "\n",
    "groups2.ug25abv.mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "## Exercises\n",
    "\n",
    "You'll now complete some exercises so you can practice what you've learned.  Your job, in general, is to replicate the analysis above, but for a different column in the dataset.  You'll be analyzing the column `PCTPELL`, which denotes the percentages of students receiving Pell Grants, a federal grant program to help students pay for college.  Unlike loans, Pell Grants do not need to be repaid.  \n",
    "\n",
    "While you'll be using a lot of the code I provided for you and modifying it slightly, make sure that you know what the code is doing at every step! \n",
    "\n",
    "---\n",
    "\n",
    "__Exercise 1__:  Similar to what you did earlier in the notebook, create a data frame that is subsetted for the columns `'INSTNM', 'STABBR', 'PREDDEG', 'CONTROL', 'UGDS',` and `'PCTPELL'`, and that retains only those rows where `PREDDEG` is equal to `3`, representing primarily bachelor's-granting universities.  What is the mean, median, maximum value, and minimum value of the `PCTPELL` column?  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "__Exercise 2:__ How many null values are contained in the `PCTPELL` column in your subsetted data frame?  Once you've found this out, drop them from your data frame.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "__Exercise 3:__ Compare the means of `PCTPELL` by the different groups of `CONTROL` (the college/university type).  Draw a visualization that shows how the distributions vary.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "__Exercise 4:__ Draw a visualization that breaks your visualization from Exercise 3 down by state with `catplot`, similar to what you did earlier in the notebook; in this instance, however, compare Texas with Florida and Illinois.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "__Exercise 5__: Show how the percent of students receiving Pell Grants varies by institution size and institution type (public, private, for-profit).  Break up the `UGDS` column into five quantiles in your data frame as you did before, and then compare the means of `PCTPELL` by institution type by institution size.  Draw a visualization with `catplot` to show this graphically.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "bb36249fbdba1d6ade34ba4fc02bf8ccd33de44ac7660b9968452a14add822d2"
  },
  "kernelspec": {
   "display_name": "Python 3.9.5 64-bit ('python39': conda)",
   "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.9.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}