{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.05656,
"end_time": "2020-05-08T03:34:32.329300",
"exception": false,
"start_time": "2020-05-08T03:34:32.272740",
"status": "completed"
},
"tags": []
},
"source": [
"# User Profile Data Prep\n",
"\n",
"This notebook does the data preparation for Bayesian inference for the model analysis."
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.05407,
"end_time": "2020-05-08T03:34:32.438247",
"exception": false,
"start_time": "2020-05-08T03:34:32.384177",
"status": "completed"
},
"tags": []
},
"source": [
"## Setup"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"papermill": {
"duration": 0.06991,
"end_time": "2020-05-08T03:34:32.565188",
"exception": false,
"start_time": "2020-05-08T03:34:32.495278",
"status": "completed"
},
"tags": []
},
"outputs": [],
"source": [
"from pathlib import Path"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"papermill": {
"duration": 2.090999,
"end_time": "2020-05-08T03:34:34.711129",
"exception": false,
"start_time": "2020-05-08T03:34:32.620130",
"status": "completed"
},
"tags": []
},
"outputs": [],
"source": [
"import pandas as pd\n",
"import numpy as np\n",
"import seaborn as sns\n",
"import matplotlib.pyplot as plt\n",
"from plotnine import *\n",
"import ujson"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"papermill": {
"duration": 0.194735,
"end_time": "2020-05-08T03:34:34.960602",
"exception": false,
"start_time": "2020-05-08T03:34:34.765867",
"status": "completed"
},
"tags": []
},
"outputs": [],
"source": [
"from bookgender.config import data_dir\n",
"from bookgender.nbutils import *"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"papermill": {
"duration": 0.076641,
"end_time": "2020-05-08T03:34:35.092662",
"exception": false,
"start_time": "2020-05-08T03:34:35.016021",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"using figure dir figures/ProfileData\n"
]
}
],
"source": [
"fig_dir = init_figs('ProfileData')"
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.055835,
"end_time": "2020-05-08T03:34:35.204268",
"exception": false,
"start_time": "2020-05-08T03:34:35.148433",
"status": "completed"
},
"tags": []
},
"source": [
"## Load Data\n",
"\n",
"We need to load the author-gender information:"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"papermill": {
"duration": 7.393277,
"end_time": "2020-05-08T03:34:42.657628",
"exception": false,
"start_time": "2020-05-08T03:34:35.264351",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/plain": [
"count 12234574\n",
"unique 6\n",
"top male\n",
"freq 3645216\n",
"Name: gender, dtype: object"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"book_gender = pd.read_csv(data_dir / 'author-gender.csv.gz')\n",
"book_gender = book_gender.set_index('item')['gender']\n",
"book_gender.describe()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"papermill": {
"duration": 5.524567,
"end_time": "2020-05-08T03:34:48.236303",
"exception": false,
"start_time": "2020-05-08T03:34:42.711736",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/plain": [
"[male, unlinked, female, unknown, ambiguous]\n",
"Categories (5, object): [male, unlinked, female, unknown, ambiguous]"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"book_gender[book_gender == 'no-viaf-author'] = 'unlinked'\n",
"book_gender[book_gender == 'no-loc-author'] = 'unlinked'\n",
"book_gender[book_gender == 'no-loc-book'] = 'unlinked'\n",
"book_gender = book_gender.astype('category')\n",
"book_gender.unique()"
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.057075,
"end_time": "2020-05-08T03:34:48.350561",
"exception": false,
"start_time": "2020-05-08T03:34:48.293486",
"status": "completed"
},
"tags": []
},
"source": [
"And we load book hashes, to set up our dummy bias:"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"papermill": {
"duration": 18.899345,
"end_time": "2020-05-08T03:35:07.307688",
"exception": false,
"start_time": "2020-05-08T03:34:48.408343",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" nisbns | \n",
" md5 | \n",
" dcode | \n",
"
\n",
" \n",
" item | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" 17 | \n",
" 3781b82fabd530590c70cac955b52bb0 | \n",
" 0 | \n",
"
\n",
" \n",
" 1 | \n",
" 2 | \n",
" 4c6606ab43bfbe946a436c0ce7633a7a | \n",
" 0 | \n",
"
\n",
" \n",
" 2 | \n",
" 38 | \n",
" e16249d40bf94b35d8a784d73d0511c5 | \n",
" 1 | \n",
"
\n",
" \n",
" 3 | \n",
" 2 | \n",
" 289071ab1041c090ac252616a76fe079 | \n",
" 1 | \n",
"
\n",
" \n",
" 4 | \n",
" 4 | \n",
" 7308735b39347b616ee6be0ab093541e | \n",
" 0 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" nisbns md5 dcode\n",
"item \n",
"0 17 3781b82fabd530590c70cac955b52bb0 0\n",
"1 2 4c6606ab43bfbe946a436c0ce7633a7a 0\n",
"2 38 e16249d40bf94b35d8a784d73d0511c5 1\n",
"3 2 289071ab1041c090ac252616a76fe079 1\n",
"4 4 7308735b39347b616ee6be0ab093541e 0"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"book_hash = pd.read_parquet(data_dir / 'book-hash.parquet').rename(columns={'cluster': 'item'})\n",
"book_hash['dcode'] = book_hash['md5'].apply(lambda x: int(x[-1], 16) % 2)\n",
"book_hash = book_hash.set_index('item')\n",
"book_hash.head()"
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.036231,
"end_time": "2020-05-08T03:35:07.380090",
"exception": false,
"start_time": "2020-05-08T03:35:07.343859",
"status": "completed"
},
"tags": []
},
"source": [
"Load the sample user ratings for each data set:"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"papermill": {
"duration": 1.922147,
"end_time": "2020-05-08T03:35:09.338882",
"exception": false,
"start_time": "2020-05-08T03:35:07.416735",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Set | \n",
" user | \n",
" item | \n",
" timestamp | \n",
" nactions | \n",
" first_time | \n",
" last_time | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" AZ | \n",
" 4975592 | \n",
" 462114 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 1 | \n",
" AZ | \n",
" 4975592 | \n",
" 1662785 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 2 | \n",
" AZ | \n",
" 4975592 | \n",
" 7287509 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 3 | \n",
" AZ | \n",
" 4975592 | \n",
" 8866889 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 4 | \n",
" AZ | \n",
" 4975592 | \n",
" 10031188 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Set user item timestamp nactions first_time last_time\n",
"0 AZ 4975592 462114 1.233101e+09 NaN NaN NaN\n",
"1 AZ 4975592 1662785 1.233101e+09 NaN NaN NaN\n",
"2 AZ 4975592 7287509 1.233101e+09 NaN NaN NaN\n",
"3 AZ 4975592 8866889 1.233101e+09 NaN NaN NaN\n",
"4 AZ 4975592 10031188 1.233101e+09 NaN NaN NaN"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"user_ratings = pd.read_csv(data_dir / 'study-ratings.csv')\n",
"user_ratings.drop(columns=['rating'], inplace=True)\n",
"user_ratings.rename(columns={'dataset': 'Set'}, inplace=True)\n",
"user_ratings.head()"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"papermill": {
"duration": 6.717506,
"end_time": "2020-05-08T03:35:16.116636",
"exception": false,
"start_time": "2020-05-08T03:35:09.399130",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Set | \n",
" user | \n",
" item | \n",
" timestamp | \n",
" nactions | \n",
" first_time | \n",
" last_time | \n",
" gender | \n",
" dcode | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" AZ | \n",
" 4975592 | \n",
" 462114 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 0.0 | \n",
"
\n",
" \n",
" 1 | \n",
" AZ | \n",
" 4975592 | \n",
" 1662785 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 1.0 | \n",
"
\n",
" \n",
" 2 | \n",
" AZ | \n",
" 4975592 | \n",
" 7287509 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 1.0 | \n",
"
\n",
" \n",
" 3 | \n",
" AZ | \n",
" 4975592 | \n",
" 8866889 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" unknown | \n",
" 0.0 | \n",
"
\n",
" \n",
" 4 | \n",
" AZ | \n",
" 4975592 | \n",
" 10031188 | \n",
" 1.233101e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 0.0 | \n",
"
\n",
" \n",
" 5 | \n",
" AZ | \n",
" 4975592 | \n",
" 43930 | \n",
" 1.233274e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 1.0 | \n",
"
\n",
" \n",
" 6 | \n",
" AZ | \n",
" 4975592 | \n",
" 11005013 | \n",
" 1.233274e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" unknown | \n",
" 1.0 | \n",
"
\n",
" \n",
" 7 | \n",
" AZ | \n",
" 4975592 | \n",
" 3649898 | \n",
" 1.233619e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 1.0 | \n",
"
\n",
" \n",
" 8 | \n",
" AZ | \n",
" 4975592 | \n",
" 8462465 | \n",
" 1.233619e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 0.0 | \n",
"
\n",
" \n",
" 9 | \n",
" AZ | \n",
" 4975592 | \n",
" 3727790 | \n",
" 1.234656e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 0.0 | \n",
"
\n",
" \n",
" 10 | \n",
" AZ | \n",
" 4975592 | \n",
" 856425 | \n",
" 1.234829e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 0.0 | \n",
"
\n",
" \n",
" 11 | \n",
" AZ | \n",
" 4975592 | \n",
" 4457778 | \n",
" 1.234829e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" unknown | \n",
" 1.0 | \n",
"
\n",
" \n",
" 12 | \n",
" AZ | \n",
" 4975592 | \n",
" 1551960 | \n",
" 1.235002e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" male | \n",
" 0.0 | \n",
"
\n",
" \n",
" 13 | \n",
" AZ | \n",
" 4975592 | \n",
" 2835913 | \n",
" 1.238026e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" unknown | \n",
" 0.0 | \n",
"
\n",
" \n",
" 14 | \n",
" AZ | \n",
" 4975592 | \n",
" 8845927 | \n",
" 1.238544e+09 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" unknown | \n",
" 1.0 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Set user item timestamp nactions first_time last_time \\\n",
"0 AZ 4975592 462114 1.233101e+09 NaN NaN NaN \n",
"1 AZ 4975592 1662785 1.233101e+09 NaN NaN NaN \n",
"2 AZ 4975592 7287509 1.233101e+09 NaN NaN NaN \n",
"3 AZ 4975592 8866889 1.233101e+09 NaN NaN NaN \n",
"4 AZ 4975592 10031188 1.233101e+09 NaN NaN NaN \n",
"5 AZ 4975592 43930 1.233274e+09 NaN NaN NaN \n",
"6 AZ 4975592 11005013 1.233274e+09 NaN NaN NaN \n",
"7 AZ 4975592 3649898 1.233619e+09 NaN NaN NaN \n",
"8 AZ 4975592 8462465 1.233619e+09 NaN NaN NaN \n",
"9 AZ 4975592 3727790 1.234656e+09 NaN NaN NaN \n",
"10 AZ 4975592 856425 1.234829e+09 NaN NaN NaN \n",
"11 AZ 4975592 4457778 1.234829e+09 NaN NaN NaN \n",
"12 AZ 4975592 1551960 1.235002e+09 NaN NaN NaN \n",
"13 AZ 4975592 2835913 1.238026e+09 NaN NaN NaN \n",
"14 AZ 4975592 8845927 1.238544e+09 NaN NaN NaN \n",
"\n",
" gender dcode \n",
"0 male 0.0 \n",
"1 male 1.0 \n",
"2 male 1.0 \n",
"3 unknown 0.0 \n",
"4 male 0.0 \n",
"5 male 1.0 \n",
"6 unknown 1.0 \n",
"7 male 1.0 \n",
"8 male 0.0 \n",
"9 male 0.0 \n",
"10 male 0.0 \n",
"11 unknown 1.0 \n",
"12 male 0.0 \n",
"13 unknown 0.0 \n",
"14 unknown 1.0 "
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"user_ratings = user_ratings.join(book_gender, on='item', how='left')\n",
"user_ratings['gender'].fillna('unlinked', inplace=True)\n",
"user_ratings = user_ratings.join(book_hash['dcode'], on='item', how='left')\n",
"user_ratings.head(15)"
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.060949,
"end_time": "2020-05-08T03:35:16.238148",
"exception": false,
"start_time": "2020-05-08T03:35:16.177199",
"status": "completed"
},
"tags": []
},
"source": [
"Now we will summarize user profiles:"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"papermill": {
"duration": 0.089575,
"end_time": "2020-05-08T03:35:16.389783",
"exception": false,
"start_time": "2020-05-08T03:35:16.300208",
"status": "completed"
},
"tags": []
},
"outputs": [],
"source": [
"def summarize_profile(df):\n",
" gender = df['gender']\n",
" dc = df['dcode']\n",
" data = {\n",
" 'count': len(df),\n",
" 'linked': np.sum(gender != 'unlinked'),\n",
" 'ambiguous': np.sum(gender == 'ambiguous'),\n",
" 'male': np.sum(gender == 'male'),\n",
" 'female': np.sum(gender == 'female'),\n",
" 'dcknown': dc.count(),\n",
" 'dcyes': dc.sum(skipna=True),\n",
" 'PropDC': dc.mean()\n",
" }\n",
" data['Known'] = data['male'] + data['female']\n",
" data['PropFemale'] = data['female'] / data['Known']\n",
" data['PropKnown'] = data['Known'] / data['count']\n",
" return pd.Series(data)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"papermill": {
"duration": 79.232407,
"end_time": "2020-05-08T03:36:35.684635",
"exception": false,
"start_time": "2020-05-08T03:35:16.452228",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" | \n",
" count | \n",
" linked | \n",
" ambiguous | \n",
" male | \n",
" female | \n",
" dcknown | \n",
" dcyes | \n",
" PropDC | \n",
" Known | \n",
" PropFemale | \n",
" PropKnown | \n",
"
\n",
" \n",
" Set | \n",
" user | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" AZ | \n",
" 529 | \n",
" 8 | \n",
" 8 | \n",
" 2 | \n",
" 1 | \n",
" 4 | \n",
" 8 | \n",
" 3 | \n",
" 0.375000 | \n",
" 5 | \n",
" 0.800000 | \n",
" 0.625000 | \n",
"
\n",
" \n",
" 1723 | \n",
" 25 | \n",
" 24 | \n",
" 3 | \n",
" 15 | \n",
" 6 | \n",
" 25 | \n",
" 14 | \n",
" 0.560000 | \n",
" 21 | \n",
" 0.285714 | \n",
" 0.840000 | \n",
"
\n",
" \n",
" 1810 | \n",
" 14 | \n",
" 6 | \n",
" 0 | \n",
" 6 | \n",
" 0 | \n",
" 8 | \n",
" 1 | \n",
" 0.125000 | \n",
" 6 | \n",
" 0.000000 | \n",
" 0.428571 | \n",
"
\n",
" \n",
" 2781 | \n",
" 8 | \n",
" 8 | \n",
" 1 | \n",
" 5 | \n",
" 1 | \n",
" 8 | \n",
" 5 | \n",
" 0.625000 | \n",
" 6 | \n",
" 0.166667 | \n",
" 0.750000 | \n",
"
\n",
" \n",
" 2863 | \n",
" 6 | \n",
" 6 | \n",
" 0 | \n",
" 6 | \n",
" 0 | \n",
" 6 | \n",
" 4 | \n",
" 0.666667 | \n",
" 6 | \n",
" 0.000000 | \n",
" 1.000000 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" count linked ambiguous male female dcknown dcyes PropDC \\\n",
"Set user \n",
"AZ 529 8 8 2 1 4 8 3 0.375000 \n",
" 1723 25 24 3 15 6 25 14 0.560000 \n",
" 1810 14 6 0 6 0 8 1 0.125000 \n",
" 2781 8 8 1 5 1 8 5 0.625000 \n",
" 2863 6 6 0 6 0 6 4 0.666667 \n",
"\n",
" Known PropFemale PropKnown \n",
"Set user \n",
"AZ 529 5 0.800000 0.625000 \n",
" 1723 21 0.285714 0.840000 \n",
" 1810 6 0.000000 0.428571 \n",
" 2781 6 0.166667 0.750000 \n",
" 2863 6 0.000000 1.000000 "
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"profiles = user_ratings.groupby(['Set', 'user']).apply(summarize_profile)\n",
"profiles = profiles.apply(lambda s: s if s.name.startswith('Prop') else s.astype('i4'))\n",
"profiles.head()"
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.061786,
"end_time": "2020-05-08T03:36:35.809283",
"exception": false,
"start_time": "2020-05-08T03:36:35.747497",
"status": "completed"
},
"tags": []
},
"source": [
"How are profile sizes distributed?"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"papermill": {
"duration": 2.877552,
"end_time": "2020-05-08T03:36:38.750065",
"exception": false,
"start_time": "2020-05-08T03:36:35.872513",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAsIAAACICAYAAADtTo36AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAesUlEQVR4nO3de5RcZZnv8e+vu9O5QyBBSAKYINEzcEYRIw4DOoIzEtFDdA5K1HHAw4jrCONtxlmgR4flGmZEZoDhiLpQOCiCwCCMOSoil3AEQSBczAUMBBIhEHOBEALk1l3P+WO/lVTavlRXV9Xe1fX7rLVX79q1L09VPb3rqb3fvV9FBGZmZmZm7aYj7wDMzMzMzPLgQtjMzMzM2pILYTMzMzNrSy6EzczMzKwtuRA2MzMzs7bkQtjMzMzM2lLbFcKSviRpuaQlkh6R9LYh5j9N0ow6bv9zkrZJ2rti2nkplvLwuKReSZPqtV0bnrzyRNIsSVvTNn8j6R5Jb0jP/aWk2yvmPTbN19XPeu6UtKIip24YaWw2PKMkh+aONB6rXp7fT5LmSPqJpCclPShpkaR3VGxnQ4rpt5I+N8h6VktaWrHvuaQe8Vl1RlEOTatHTFWJiLYZgKOBe4Gx6fE0YMYQy9wJzK1jDPcDdwGnDTLP1cA/5f1+teuQZ54As4BlFY8/CXyv4vFPgY8AXcAS4E8bGY8H51De72W7DDnnzDjgceCkimn/tfw9BZwGfCONTwU2AgcNsK7VwLS83892HJxDtQ1/cBRglJsObIyI7QARsbH8hKS3ABcCk8g+oNOAY4C5wNWStgJHR8TWWjcu6XVp/V8Avghc2c88fwUcCpxa63ZsxHLNkz72AjZVPP5b4DbgcOCBiLinTtux+nIO2XDlmTMfBe6NiIXlCRGxDFjWd8aIeF7SyhTvMzVuzxrDOVSLvH/BNHMgS4BHyH61fBP4szR9DHAPsF96fApwRQzxa4msoH2kn+GSAeb/X8CXyZqkrAZe0+f5WcAa4NC836t2HvLMk5QDW9PzTwJrgYP7zPMvwGYG+cWc4llRsa0L8n5f22kYRTnkI8LtkTMXAp8ZJLbT2H007+C0nnEDzLsaWFqxvc/l/d62yzDKcshHhBshIl5Ov4reDhwHXCfpbGAx2SH8WyUBdJJ9eQy1vguAC4YRwgLgAxFRknQj8EHgUgBJncAPgC9HxMphrNPqrAB58mREHAEg6RTgMmBeetwB/DnwMvBasl/2A/loRCwexnatTkZRDlmTFCBndpF0EzAHeDwi/jJNPkXSccAbgE9ExLZBVnFcVByNtOYYZTnUNG1VCANERC/ZL6A7JS0la4LwILA8Io4ezrokfYHsdEBfv4yIT/eZ941kSVFOxG7gKVIhTHa0eG1E/J/hxGCNkVee9GMhUJkTZ5KdavoycKmkoyMiJN0C7A8sjoi/GU581hjOIRuuHHNmOfCOijg+kC6U/NeKea6LiLMkHQ38VNLNwIYUH8DCiPjKcGK0+nMODV9bFcLpyulSRDyRJh0B/I7sFPJ+6QvhXkljgNdHxHJgCzC5v/UN89fSh4FzI+JfKuJZJem1ZO1kTgOOrOFlWZ3lnCd9HUt2ehtJBwCfB46KiA2SPgH8DfCdiDihxvVbAziHbLhyzplrgHMknRS723hOGGC990q6iuw0+DkpTisA51Bt2qoQJms/878lTQF6gJXAGRGxQ9LJwCXKbmvWBVxM9gvnSuDbdWhIvgB4T59pN6Xpx5MlzKJ0tLjsv0fEkzVuz2qXZ54AvE7SI4CAHWSFCmRtsL4eERvS488Cd0n6UUS80M96yhdAQHYBxZ+PICYbntGSQ9Y8ueVMRGyV9D7gQkkXA+vICqR/GmCR84GHJP1zRGzp5/lFknrT+JKI+Ota4rJhG0051DRKDZPNzMzMzNpK23WoYWZmZmYGLoTNzMzMrE25EDYzMzOztuRC2MzMzMzaUuEK4Xnz5gXgoTnDqOLcce6MhPPHuVMr545zp1bOnfxzp3CF8MaN7ozGauPcsZFw/litnDtWK+dO/gpXCJtZe5E0T9IKSStTd6B9nx8r6br0/H2SZlU8d06avkLSCX2W65T0sKSfNP5VmJlZK3IhbGa5kdRJ1s34e4DDgA9LOqzPbKcDmyLiUOAishuxk+ZbABwOzAO+mdZX9hngsca+AjMza2UtUQhfc9/TeYdgLcq5U3hHASsj4qmI2AFcC8zvM8984Htp/AbgXcq6YJwPXBsR2yNiFVkvSkcBSDoQeC/w3VoDc+7YSDh/rBbX3Pe0c6fJWqIQNrNRaybwTMXjNWlav/NERA+wGZg6xLIXA/8AlAbbuKQzJC2WtHjDhg2DzWpmZqOQC2Ezy5P6mdb3Ct+B5ul3eurvfn1EPDjUxiPisoiYGxFz99tvv6GjNTOzUcWFsJnlaQ1wUMXjA4HnBppHUhewN/DCIMseA5wkaTVZU4vjJf2gEcGbmVlrcyFsZnl6AJgjabakbrKL3xb2mWchcGoaPxm4IyIiTV+Q7ioxG5gD3B8R50TEgRExK63vjoj4q2a8GDMzay1deQdgZu0rInoknQXcAnQCV0TEcklfBRZHxELgcuAqSSvJjgQvSMsul3Q98CjQA5wZEb25vBAzM2tJLoTNLFcR8TPgZ32mfaVifBvwwQGWPQ84b5B13wncWY84zcxs9HHTCDMzMzNrSy6EzczMzKwtuRA2MzMzs7bkQtjMzMzM2pILYTMzMzNrSy6EzczMzKwtuRA2MzMzs7bkQtjMzMzM2pILYTMzMzNrSy6EzczMzKwtuRA2MzMzs7bkQtjMzMzM2pILYTMzKzRJ8yStkLRS0tn9PD9W0nXp+fskzap47pw0fYWkEyqmr5a0VNIjkhZXTN9X0q2Snkh/92n06zOz/LgQNjOzwpLUCVwKvAc4DPiwpMP6zHY6sCkiDgUuAs5Pyx4GLAAOB+YB30zrKzsuIo6IiLkV084Gbo+IOcDt6bGZjVIuhM3MrMiOAlZGxFMRsQO4FpjfZ575wPfS+A3AuyQpTb82IrZHxCpgZVrfYCrX9T3g/XV4DVYwDTrLcIWk9ZKWNedVWD24EDYzsyKbCTxT8XhNmtbvPBHRA2wGpg6xbAC/kPSgpDMq5tk/Itamda0FXtNfUJLOkLRY0uINGzbU9MIsHw08y3BlmmYtxIWwmZkVmfqZFlXOM9iyx0TEkWTF0JmS3jGcoCLisoiYGxFz99tvv+EsavlryFmGiPgl8EIzXoDVjwthMzMrsjXAQRWPDwSeG2geSV3A3mQFyYDLRkT573rgJnY3mVgnaXpa13Rg/XADfvHVHXzxpqVcv/iZoWe2PDTqLENVfDahWFwIm5lZkT0AzJE0W1I32WnphX3mWQicmsZPBu6IiEjTF6T2nrOBOcD9kiZKmgwgaSLwbmBZP+s6FfjxcAN+auMrAFx939PDXdSao1FnGariswnF0pV3AGZmZgOJiB5JZwG3AJ3AFRGxXNJXgcURsRC4HLhK0kqyI8EL0rLLJV0PPAr0AGdGRK+k/YGbsjPddAHXRMTP0ya/Blwv6XTgaeCDw415R08JgLGdPtZUUMM5y7Cm2rMM1ppcCFtdSJoH/DvZF9V3I+JrfZ4fC3wfeAvwPHBKRKxOz51DdmFCL/DpiLglTV8NbEnTe/rc4sjM2kRE/Az4WZ9pX6kY38YABWtEnAec12faU8CbBpj/eeBdI4l3285eALq7XAgX1K6zDMCzZD+cPtJnnvKZgXupOMsgaSFwjaQLgRmkswxNi9zqzoWwjVjFFbh/QfZr+QFJCyPi0YrZdl2BK2kB2RW4p/S5AncGcJuk10dEb1ruuIjY2LQXY2Y2Qi9v7wFcCBdVI84yAEj6IfBOYJqkNcA/RsTlTX55NkxVFcK1Hu2TNJXsasu3AldGxFn1DN4KY9cVuACSylfgVhbC84Fz0/gNwDf6XoELrEo7naPIfoWbmbWcl7amQthNIwqr3mcZ0vQP1zlMa4Ih/0tHcr89YBvwZeDv6xaxFVGz7/O5i6++NbOieWnbTgDGjnEhbFZ01fyX1ny/vYh4JSLuJiuIbfTK7T6fvvrWzIrmpa1ZIdyh/nZvZlYk1RTCIznaVxUf1Wt5zb7Pp5lZYZWPCJcvmjOz4qqmEB7J0b6q+Khey2v2fT7NzAqr3EZ4e7qNmpkVVzWF8EiO9lkbSGcBylfgPgZcX74CV9JJabbLganpYrjPA2enZZcD5Stwf87uK3D3B+6W9BuyW9P8tOI+nzaKSJonaYWklZLO7uf5sZKuS8/fJ2lWxXPnpOkrJJ2Qph0kaZGkxyQtl/SZ5r0aMx8RNmsl1dw1oub77dUzUCu2Zt7n00aPRtx6j+yWRn8XEQ+lswoPSrq1zzrNGubVHVkB7ELYrPiGPCI8kqN9sKtThAuB0ySt6eeOE2bWvmq+GJeKW+9FxCpgJXBURKyNiIcAImIL2X6r73UNZg3TW8qOA7lphFnxVXUf4REe7Zs1gvjMbHTr72Lctw00T7oRfuWt937dZ9k9Ct7UjOLNwH39bTzdlu8MgIMPPrjGl2C2p55UCPuIsFnx+SaHZpanRt16D0mTgB8Bn42Il/rbuC/UtUboLWVHgn1E2Kz4XAibWZ4acus9SWPIiuCrI+LGhkRuNoCe3vIRYRfCZkXnQtjM8tSIW++J7LqFxyLiwqa8CrMK5aYR2900wqzwqmojbGbWCKnNb/li3E7givLFuMDiiFhIVtRelS7GfYGsWCbNV771Xg/p1nuSjgU+BiyV9Eja1BfTtQ5mDeeL5cxahwthM8tVA269dzf9tx82a4qdvVkBvKO3RG8p6OxwOpoVlZtGmJmZ1VH5iDDA9h43jzArMhfCZmZmddRTUQi/vK0nx0jMbCguhM3MrLCa2QW3pHMlPSvpkTScWEvMvaVgQncnAJu37qxlFWbWJG4jbGZmhZRTF9wXRcS/jiTunb0lJnZ38eqOXl7a5kLYrMh8RNjMzIqqJbvg7i0FE8f6iLBZK3AhbGZmRdVfF9x9i9Y9uuAGKrvgHnTZAbrgPkvSEklXSNpnoMAknSFpsaTFGzZs2OO5nlIwcWx2wtWFsFmxuRA2M7OianYX3N8CXgccAawF/m2gwAbrnrunt8SE7qwQfmmrL5YzKzIXwmZmVlRN7YI7ItZFRG9ElIDvkDXNGJZSKSgFbhph1iJcCJuZWVE1tQtuSdMrHn4AWDbcgHsjO+g8prODCd2dvORC2KzQfNcIMzMrpBy64P66pCPImlCsBj453JjLnWl0SOw9foyPCJsVnAthMzMrrGZ2wR0RHxtpvOXulTsEe41zIWxWdG4aYWZmVieVR4R39JZY8fstOUdkZoNxIWxmZlYn5e6VOzvEpLFdvLzdd40wKzIXwmZmZnXS07v7iPCkcS6EzYrOhbCZmVmd9JR2txGeNLaL7T0ltu3szTkqMxuIC2EzM7M66e3TNALg+Vd25BmSmQ3ChbCZmVmd9FRcLFcuhDdu2Z5nSGY2CBfCZmZmdbKrjfAeR4RdCJsVlQthMzOzOim3Ee4U7D1hDAC/e/7VPEMys0G4EDYzM6uTyvsITx7bxcSxXSx/7qWcozKzgbhnOTMzszrZWdE0QhIzp4xzIWxVe+SZTdyyfB1TJozhxD+ennc4bcFHhM3MzOqk8ogwwIwp43li3Ra2bHNXyza0pWs2s3nrTu5f9ULeobQNF8JmZmZ1sus+wunbdc5rJtNTCu56YmOOUVmr2JruOb1+y7acI2kfLoTNzMzqZNd9hNMR4YP3ncCE7k5ue3RdnmFZiygXwute8p1GmqUlCuGnNrzMxy6/j1fcVaWZmRVYZRfLkHWs8Yb9J3PHivX09JbyDM1awNYdPiLcbIUvhHt6S3z37lXc9cRGHn76xbzDMTMzG9CuDjU6tGvaH03fixdf3cntv12fV1jWIiqPCEdEztG0h8IXwms37/5VtOy5zTlGYq1o0ys7uOyXT7K9pzfvUMysDfSW2wjvroP5LwdM5vX7T+LL/7mM32/2kT7r346eEjt7gwndnezoKfHSVp8Fb4bCF8JrNm3dNb70WRfCVr2I4OLbH+eff/ZbfrHc7fNseH54/9N84Ju/Yu3mrUPPbJbs7NM0AqCrs4N3H3YAz7+8g7d//Q4W+ciw9WPz1uzOIgftMwGAOx93njRDVYWwpHmSVkhaKensfp4fK+m69Px9kmZVPHdOmr5C0gnDDfCZTVmPPLOmTuQXy3/Psy/6S6mIGpEjQ61zKGs2bd31pXS3r9gurCLmzradvSx9djMPP/0iF9yyotaXZnXSzByRNDut44m0zu7hxPrEui3A7ovlymZMGc8n3j6byePG8PErH+DMax5iyZoX2dHjdsN5KOJ+p9wV9xsP3Jv99xrL127+Letf8hmERhuyQw1JncClwF8Aa4AHJC2MiEcrZjsd2BQRh0paAJwPnCLpMGABcDgwA7hN0usjoqrz1K/u6OFLNy0FYP4RM7h00UreecEiTjj8AN5/xEwkmDVtInuPH8OksV10SPSWgnFjOna10xLZr/GBRASlgFIEpQgiICK79c3Yrs4Bl0nvTTUvIzcRQU8p6C0FO3tLTB43piHbaUSOpGWGWueASqXg0kUrAdhnwhhufHgNM6aM59g5U5k5ZQIvbt3BgftMYMKYzj3a8llzFTF3AC689XEApk3q5saHnmVidxcz9xnPmw6cwuxpE1m7eSuz074HYHtPiXFjdu8vdvSU6NDg+57RrLxf3d7Ty4TukfXblEOOnA9cFBHXSvp2Wve3qon12Re3cskd2X6nv/3KwVMn8sl3HMK9Tz7Pz5f+np8uWcu0SWN53xunM3ViN3P2n8x+k7sZN6aTA/YahyQigiC7G8WmV3fw2n0n0hvB2K4OxoyC/CqVgp2lEhHZ/9GOnhLPvbiVJWte5ENvPWjA7+GRKOJ+Z/OrO5l38V0ATOju4oNvOYjL717FiZfcxfwjZnLiHx/AayaPY3tPiel7j6NDoqMj+8GVjbff91i5FitFdtvCUgm2bN/Jmk1beePMvave/1azhzoKWBkRTwFIuhaYD1R+uPOBc9P4DcA3lFWJ84FrI2I7sErSyrS+e6sJ7tZH1+06orf/XuM46/hDuX/VC9z+2Hp+smTtgMtJWTErZUkyprNjd5HL7sJ3qHbo3Z0du+4FKbRr3T2psOzu7KBvLVyerzKWRii/luwvUPG4/A+yPR1pyKbB3737DZx53KGNCKcROUIV6xzQr57cyLUPPAPA6ccewo0Pr+Gi2x7notv2nG9Mp+geBV8mjfahtx7EP/63wxux6sLlzvLnNnPZL58C4LQ/nc1PlzzHD+9/eteP60odyu4KsLM3+wHeW8r2L+VCcPyYzobtA6pV3s9l5VTfafQZ2T3fH8zD7i+ePaf94Ta7OkRv2sfOP2IG/77gzTXHTxNzRNJjwPHAR9I830vrraoQnjllPJPGdvHy9h46ByhMJo8bw7sPP4A3HTSF3z3/Kr96ciNX3rO6mtX/gfFjOmn1+mdnb7BjgLtp3Lzs9/zg9Lc1osgr3H5nr/FdTN97HGs3b2N8dyczpozn9GNnc9tj67jyntVcfveqIdfR2SE6lNUAA+VfK9mjzond+6YOCSk74KB0ALRc95X90fS9+PGZx9DdNfT3ezWF8EzgmYrHa4C3DTRPRPRI2gxMTdN/3WfZmX03IOkM4Iz08GVJlecipwEbP3p+FZHW3zQgr3Pqdd/2WV+Ds/ac9POImFeHVTcqR4Za51C5AzDtb8/P7TPsK8986k/V8ZzL7m+EpOVzB4be9xQkd4qUN8OO5RLgkg/vMWm4udPMHJkKvBgRPf3Mv4ehcudTxcidspbMod8BPzxjj0ltsd85p1i5M1K55N7vgLGf3WPSgLlTTSHc38+KvscABpqnmmWJiMuAy/rduLQ4IuYOFWQjtOu2a9CIHOnvZ9ywcgeK9T4WKRYoTDy55Q4Ud99TxDggt1iamSNVfWdBa+ROWZHiKUgs3u80SSu8nmrOCa8BDqp4fCDw3EDzSOoC9gZeqHJZa32NyBHnTntw7thQmpkjG4EpaR0Dbctan/c7tks1hfADwJx0JW03WSPxhX3mWQicmsZPBu6IrDHZQmBBuvpyNjAHuL8+oVuBNCJHqlmntT7njg2laTmSllmU1kFa548b+NosH97v2G4RMeQAnAg8DjwJfClN+ypwUhofB/wHsJIsIQ6pWPZLabkVwHuq2V6fbZ8x3GXqNbTrtmuMt+450t86W/l9LFIsRYrHudMaceQZSzNzBDgkrWNlWufYVnmfWiGeosTi/U57fd6DDUqBmpmZmZm1Fd83yszMzMzakgthMzMzM2tLhS6ENcLuCgdY5xWS1ktaVjFtX0m3KutS81ZJ+6TpknRJ2v4SSUdWLHNqmv8JSaf2t60+2z1I0iJJj0laLukzzdp2O2pE7vSzjYE+03MlPSvpkTScWLFMQ7rmrFjPaklL03YXp2nOsWFoUu7k9jk1eh8o6S3pta1My7b+nf2HIed9T27/65I6JT0s6Sfp8Wz101W1aui6uF00I3dq0Yx8y3W/kXcj5UEaWHeSNTg/BOgGfgMcVof1vgM4ElhWMe3rwNlp/Gzg/DR+InAz2X0D/wS4L03fF3gq/d0nje8zxHanA0em8clkDeoPa8a2221oVO4M4zM9F/j7fuY/LMUyFpidYuysZ7zAamBan2nOseLlTm6fEw3eB5JdWHR0WuZmarhIulWHAux7cvtfBz4PXAP8JD2+HliQxr8N/M80/ing22l8AXBdGu93/5j3Zzracqeo+ZbnfqPIR4R3dYEYETuAcneFIxIRvyS7F2Cl+WRdaZL+vr9i+vcj82uy+0tOB04Abo2IFyJiE3ArMGhvNxGxNiIeSuNbgMfIeqNp+LbbUENyp69BPtOB7OqaMyJWkV2NfFQT4nWOVa8puTOApnxOjdwHpuf2ioh7I/t2+37FutpB3vueXP7XJR0IvBf4bnossq6qbxgglnKMNwDvSvMPtH9sF3nuewbV6HzLe79R5EK4vy4QBysyRmL/iFgL2QcOvGaIGEYUWzoV9GbgvmZvu000/T3q85kCnJVOCV1RPl00SFz1jDeAX0h6UFk3nuAcG45mvfaifU712vbMNF6PmFpR3vuevHLoYuAfgFJ6PFhX1Xt0XQxUdl3crvsdaJHX36B8y3W/UeRCuOquLnOIoebYJE0CfgR8NiJeaua220hT36N+PtNvAa8DjgDWAv82RFz1jPeYiDgSeA9wpqR3DBZ6E+JpNc167a3yOQ132+2cO5D/vmfAWfuZVpfPS9L7gPUR8WAV22toLC2u8K+/gfmW62svciHczO4K16VD86S/64eIoabYJI0hS6KrI+LGZm67zTTtPervM42IdRHRGxEl4DvsPr3X8M80Ip5Lf9cDN6VtO8eq15TXXsDPqV7bXpPG6xFTK8p130M+OXQMcJKk1WSn848nO0I8UFfVw+26uF0U+vU3ON/y3W9U05A4jwHoImtIPZvdDccPr9O6Z7HnhSIXsGeD76+n8feyZ4Pv+2N3g+9VZI2990nj+w6xTZG1e7m4z/SGb7vdhkbmTpWf6fSK8c+RtXsDOJw9LwZ5iuwCibrEC0wEJleM30PW3s85VqDcKcLnRAP3gWRdzf4Juy96OTHvz3U05U/aTiG/T4B3svtiuf9gz4vlPpXGz2TPi+WuT+P97h/z/kxHW+4UNd/y3G/k/gYP8eaPuLvCftb5Q7LT1TvJfoWcTtY+6XbgifS3/MEIuDRtfykwt2I9/4OsMf9K4ONVbPdYskP9S4BH0nBiM7bdjkMjcmcYn+lV6TNbQtbXfGVh3MiuOQ9JO8/fAMvZ3W2oc6xAuZP359TofSAwF1iWlvkGZD2YtsuQ874n1/919iyE++2qmhq6Lm6XoRm5U9R8y3O/4S6WzczMzKwtFbmNsJmZmZlZw7gQNjMzM7O25ELYzMzMzNqSC2EzMzMza0suhM3MzMysLbkQbjJJn5U0Ie84rPU4d6xWzh2rlXPHatUquePbpzVZ6n1nbkRszDsWay3OHauVc8dq5dyxWrVK7viIcD8k/bWkJZJ+I+kqSa+VdHuadrukg9N8V0o6uWK5l9Pfd0q6U9INkn4r6WplPg3MABZJWpTPq7NGcu5YrZw7VivnjtXKuUOxe5bLqQeVw8l6tZmWHu8L/F/g1IpeUf4zjV8JnFyx7Mvp7zuBzWT9ZXcA9wLHpudWl9ftYXQNzh0Pzh0Pzh0PrTI4d7LBR4T/0PHADZEO5UfEC8DRwDXp+avIuhscyv0RsSYiSmTdEc5qQKxWLM4dq5Vzx2rl3LFaOXdw04j+iKxP7cGUn+8hvYeSBHRXzLO9YrwX6KpXgFZYzh2rlXPHauXcsVo5d3Ah3J/bgQ9JmgogaV/gHmBBev6jwN1pfDXwljQ+HxhTxfq3AJPrFawVinPHauXcsVo5d6xWzh1arGpvhohYLuk84P9J6gUeBj4NXCHpC8AG4ONp9u8AP5Z0P1lCvVLFJi4Dbpa0NiKOq/8rsLw4d6xWzh2rlXPHauXcyfj2aWZmZmbWltw0wszMzMzakgthMzMzM2tLLoTNzMzMrC25EDYzMzOztuRC2MzMzMzakgthMzMzM2tLLoTNzMzMrC39f0xthaMduTcBAAAAAElFTkSuQmCC\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"g = sns.FacetGrid(profiles.reset_index(), col='Set', sharex=False, sharey=False, height=2)\n",
"g.map(sns.distplot, 'count')\n",
"plt.savefig(fig_dir / 'profile-size-all.pdf')"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"papermill": {
"duration": 2.759299,
"end_time": "2020-05-08T03:36:41.571838",
"exception": false,
"start_time": "2020-05-08T03:36:38.812539",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAsgAAACICAYAAAD6bB0zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAe+0lEQVR4nO3debxcZZ3n8c/3bgmQsGRhkS0IgRFsRYwIzdLSKgR0wLXBtgUcFF/zgh6XVhu08eXQ0C3NDDooOoMNLaZVcFTatAsoLQwgYQnKFmLgEgLGAElIzEa2e+s3f5ynQqWoe2/l3lP3nLr3+369zivnnq2ec+qXOr966nmeo4jAzMzMzMwyHUUXwMzMzMysTJwgm5mZmZnVcIJsZmZmZlbDCbKZmZmZWQ0nyGZmZmZmNZwgm5mZmZnVGPMJsqTPS1og6RFJD0l68xDbnyvpVTm+/iclbZK0W82yy1NZqtMTkvolTcrrdW1kioobSTMkbUyv+bCkeyQdlta9R9J/1Gx7fNquq8Fx7pC0qCbGfjDSstnwjZF4mjXS8tjwFXkvkzRT0k8kPSXpQUm3Szqx5nVWpDL9TtInBznOEkmP1nwuXZ1H+WzHjaF4mpZHmRqKiDE7AccC84AJ6e9pwKuG2OcOYFaOZbgfuAs4d5BtvgNcVvT18lR83AAzgMdq/v4YcEPN3z8F/hLoAh4B/rSV5fHkeMqzPJ7aMoYmAk8Ap9cse231ngacC3wtzU8FVgL7D3CsJcC0oq/neJ8cT81Nr6gpGGP2AVZGxGaAiFhZXSHpjcBVwCSyN+Bc4DhgFvAdSRuBYyNi43BfXNLB6fifAT4HfKvBNn8FHAKcM9zXsdwVGjd1dgVW1/z918BtwBHAAxFxT06vY63jeLKRKjKGPgjMi4i51QUR8RjwWP2GEfGipN5U3t8P8/Ws9RxPzSj6m0wrJ7I3+CGybytfB/4sLe8G7gGmp7/PBK6PIb4lkSW6DzWYrh5g+78DLiFryrIE2LNu/QxgKXBI0dfKUzniJsXExrT+KeA54IC6bf4RWMMg35xTeRbVvNaVRV/X8TqNoXhyDfL4jKGrgI8PUrZzebnG74B0nIkDbLsEeLTm9T5Z9LUdj9MYiyfXIA9HRKxP34ZOAE4CbpJ0ETCfrEr/l5IAOsluHEMd70rgyh0owlnAuyOiIulHwPuBawAkdQL/ClwSEb07cExrsRLEzVMRcSSApDOBa4HZ6e8O4G3AeuBAsm/4A/lgRMzfgde1FhhD8WQFKUEMbSPpZmAm8EREvCctPlPSScBhwEcjYtMghzgpamosbfSNsXhqmTGdIANERD/ZN587JD1K1pThQWBBRBy7I8eS9Bmynwfq3RkR/61u29eRvenVQOsBFpMSZLLa5eci4l92pAw2OoqKmwbmArUxcgHZT1GXANdIOjYiQtKtwF7A/Ij4yI6Uz1rP8WQjVWAMLQBOrCnHu1OHzf9Rs81NEXGhpGOBn0r6ObAilQ9gbkR8YUfKaK3leBramE6QU2/tSkQ8mRYdCTxD9tPz9HQzmCepGzg0IhYA64DJjY63g9+SPgB8MSL+saY8T0s6kKw9zbnAUcM4LWuxguOm3vFkP40jaW/gU8DREbFC0keBjwDfjIhThnl8azHHk41UwTH0XeBiSafHy+1Gdx7guPMkzSH7Cf3iVE4rGcdTc8Z0gkzWzuarknYH+oBe4PyI2CLpfcDVyoZf6wK+QvbN5lvA/86hIfpZwKl1y25Oy/+cLCBuT7XLVe+NiKeG+XqWnyLjBuBgSQ8BAraQJS2Qtd36p4hYkf7+BHCXpB9GxKoGx6l2qICsQ8bbRlAmG76xEk9WnMJiKCI2SnoncJWkrwAvkCVLlw2wyxXAbyT9Q0Ssa7D+dkn9af6RiDh7OOWyERlL8dQySg2dzczMzMyMcfCgEDMzMzOzHeEE2czMzMyshhNkMzMzM7MaTpDNzMzMzGqULkGePXt2AJ7ymcY1x5JjKS+OJcdRHhxHjqW8OJZaH0ulS5BXrvQDdiwfjiXLi2PJ8uA4srw4llqvdAmymZmZmVmRnCCbmZmZmdUofYL83fue5bv3PVt0MWwMcBxZHvyZZHlyLFke/LmUv9InyGZmZmZmo8kJspmZmZlZDSfIZmZmZmY1nCBb7iTNlrRIUq+kixqsnyDpprT+Pkkz0vKpkm6XtF7S1+r2eaOkR9M+V0vS6JyNmZmZjTdOkC1XkjqBa4BTgcOBD0g6vG6z84DVEXEI8GXgirR8E3AJ8OkGh/4GcD4wM02z8y+9mZmZmRNky9/RQG9ELI6ILcCNwBl125wB3JDmfwC8VZIiYkNE3E2WKG8jaR9g14iYFxEBfBt4V0vPwszMzMYtJ8iWt32B39f8vTQta7hNRPQBa4CpQxxz6RDHBEDS+ZLmS5q/YsWKHSy6mZmZmRNky1+jtsH1zzpvZpthbR8R10bErIiYNX369EEOaWZmZtaYE2TL21Jg/5q/9wOWDbSNpC5gN2DVEMfcb4hjmpmZmeXCCbLl7QFgpqSDJPUAZwFz67aZC5yT5t8H/Cq1LW4oIp4D1kk6Jo1ecTbw4/yLbmZmZgZdRRfAxpaI6JN0IXAr0AlcHxELJF0KzI+IucB1wBxJvWQ1x2dV95e0BNgV6JH0LuDkiHgc+K/At4CdgJ+nyczMzCx3TSXIkmYD/4ss4fnniPhS3foJZCMLvBF4ETgzIpZIejvwJaAH2AJ8JiJ+lWP5rYQi4mfAz+qWfaFmfhPw/gH2nTHA8vnAa/MrpZmZmVljQzaxGOG4tiuB/xwRf0L2k/qcvApuZmZmZtYKzbRBHsm4tr+NiGpnqgXAxFTbbGZmZmZWSs0kyHmNa/te4LcRsbn+BTx2rZmZmZmVRTMJ8ojHtZV0BFmzi481egGPXWtmZmZmZdFMgjyicW0l7QfcDJwdEU+NtMBmZmZmRZA0W9IiSb2SLmqwfoKkm9L6+yTNqFl3cVq+SNIpadn+km6XtFDSAkkfH72zscE0kyAPe1xbSbsDPwUujohf51VoMzMzs9E0kkEL0nZnAUcAs4Gvp+P1AX8TEa8BjgEuaHBMK8CQCXJqU1wd13Yh8P3quLaSTk+bXQdMTePafgqofqu6EDgEuETSQ2naM/ezMLMxKe/ampp1nZJ+K+knrT8LMxsjhj1oQVp+Y0RsjoingV7g6Ih4LiJ+AxAR68jyrPp+XlaApsZBHu64thFxGXDZCMtoZuNQTW3N28macT0gaW56cEzVttoaSWeR1dacWVdb8yrgNkmHRkR/2u/jZDeiXUfpdMys/TUatODNA22THpxVHbRgX+Deun23S4TTF/w3APc1enFJ5wPnAxxwwAHDPAVrlh81bWZllXttDWzrF/EO4J9H4RzMbOwYyaAFQw1mMAn4IfCJiFjb6MU9oMHocoJsZmU1kiEmB9v3K8BngcpgL+7hJ82szkgGLRhwX0ndZMnxdyLiRy0pue0wJ8hmVla519ZIeiewPCIeHOrFXVtjZnWGPWhBWn5W6jdxEDATuD/94nUdsDAirhqVs7CmOEE2s7JqRW3NccDpkpaQNdn4c0n/2orCm9nYMpJBCyJiAfB94HHgFuCC1CfiOOBDZJ9F1cEMThvVE7OGmuqkZ2ZWgG21NcAfyGpr/rJum2ptzTy2H2JyLvBdSVeRddKbCdwfEfOAiwEkvQX4dET81WicjJm1v+EOWpDWXQ5cXrfsbhr/4mUFc4JsZqWUeoBXa2s6geurtTXA/IiYS1ZbMyfV1qwiS6JJ21Vra/p4ubbGrFQqEUQE2S/tZlYWTpDNrLTyrq2pW38HcEce5TQbrr/7t8e4/+lVXP2BNxRdFDOr4TbIZmZmBahUsj6ncx+ub1pvZkVzgmxmZlaANRu3Fl0EMxuAE2QzM7MCvLhhc9FFMLMBOEE2MzMrwMr1WwDocP88s9JxgmxmZmOepNmSFknqlXRRg/UTJN2U1t8naUbNuovT8kWSTqlZvkTSo2ns2vk7WqYXU4K8ywT3lzcrG/+vNDOzMU1SJ3AN8Hayh8g8IGluRDxes9l5wOqIOETSWcAVwJmSDicbPvAIsjG1b5N0aM2wgSdFxMrhlKvaxGKSE2Sz0nENspmZjXVHA70RsTgitpA9RfGMum3OAG5I8z8A3poeA3wGcGNEbI6Ip4HedLwRW7Uhq0Heqbszj8OZWY6cIJuZ2Vi3L/D7mr+XpmUNt0mPFF4DTB1i3wB+IelBSec3emFJ50uaL2n+ihUrtlu3tb8CwJb0r5mVhxNkMzMb6xp1g4smtxls3+Mi4ijgVOACSSe+YsOIayNiVkTMmj59+nbrqnnxpq1OkM3Kxgmy5a6MnWHMbFxbCuxf8/d+QP3TObZtI6kL2I3s8eUD7hsR1X+XAzezg00v+itZYrxpq5+CblY2TpAtVzWdYU4FDgc+kDq51NrWGQb4MllnGOo6w8wGvp6OV3VSRBwZEbNafBpmNrY8AMyUdJCkHrLPmbl128wFzknz7wN+FRGRlp+VvtgfBMwE7pe0i6TJAJJ2AU4GHtuRQvWlJ+k5QTYrHyfIlrdSdoYxs/ErtSm+ELgVWAh8PyIWSLpU0ulps+uAqZJ6gU8BF6V9FwDfBx4HbgEuSCNY7AXcLelh4H7gpxFxy46Uqz8lyH2VoM/tkM1KxWPLWN4adWh580DbRESfpNrOMPfW7VvfGSaA/xMR17ag7GY2RkXEz4Cf1S37Qs38JuD9A+x7OXB53bLFwOtHUqZqDTLApr4KkzpdZ2VWFk39bxxum1JJUyXdLmm9pK/lW3QrqcI6w8DgPcbNzMqkv78mQXYzC7NSGTJBHkmbUmATcAnw6dxKbGVXaGeYwXqMm5mVyXY1yE6QzUqlmRrkYbcpjYgNEXE3WaJs40MpO8OYmZVNJZwgm5VVMwnySAZYb4p/Fh87ytoZxsysbGprkDdsdoLcDlo0jOn1kpZLcsVPiTSTII+kTWlT/LP42BIRP4uIQyPi4NS5hYj4QkTMTfObIuL9EXFIRBydOrtU97087XdYRPw8LVscEa9P0xHVY9rYl/fNSNL+qV/EQkkLJH189M7GbHvVcZABNmzuK7Ak1owWDmP6rbTMSqSZBHkkbUrNzIalRTejPuBvIuI1wDFkHT7rj2k2KvpqOumtc4LcDloyjGlE3IlzptJpJkEeSZtSM7Phyv1mFBHPRcRvACJiHVkzoPomY2ajor8STOzObsPrNzlBbgMjaXLazL6DcnPU0TVkgjySNqWQPSIYuAo4V9JS19aYWZNaejNKzTHeANzX6MV9M7JW66sEO/dkjyNY7xrkdtCqYUyb4uaoo6upB4WMcID1GSMon5mNXy27GUmaBPwQ+ERErG304ulhNNcCzJo1y7+IWe76K8FO3Vkz1HWbthZcGmvCjjQ5XdrsMKZWTn5sj5mVVUvG1JbUTZYcfyciftSSkps1oa9SobtTTOjqcBvk9pD7MKajVG4bBifIZlZWrRhTW2RNwhZGxFWjchZmA6hUoENi8sQut0FuAy0axhRJ3wPmAYelpqjnjeZ5WWNNNbEwMxttEdEnqXoz6gSur96MgPlp2MDrgDnpZrSKLIkmbVe9GfWRbkaSjgc+BDwq6aH0Up9LzcjMRlVfpUJHh5g0octtkNvECJucXg68YpjSiPhAzsW0HDhBNrPSyvtmlJ7s2ah9stmo668EHYLJE7tdg2xWMk6QzczMCtBXCTokNmzuo9ed9GyEnluzkYgga0lmI+U2yGZmZgXoTwnyLhO6/CQ9G5HfPruar/6ql9sXLS+6KGOGE2QzM7MC9KUmFrtO7GKtm1jYCCx8LhutctUG/xKRFyfIZmZmBeivBB0dYvLEbrb0VdxRz4Zt5fotgDtY5MkJspmZWQGqTSx23SnrDvTC2k0Fl8ja1drUht1fsvLjBNnMzKwA1QR58sRuwAmyDV+15thPZMyPE2QzM7MC9FUqaZi3rAb5+TVOkG14KpH96ycy5scJspmZWQGqbZCn7NxDh2Dxig1FF8naVF+lAsA6d/bMjRNkMzOzAlTHQe7q7GDKLhN4cvm6ootkbahSCbb2Z1XIfuBMfpwgm5mZFaC/P+hMjUf3nDyBJ5evL7ZA1pY291W2zbuTXn6cIJuZmRWgr/LyU8/23m0iS1Zu2DYagVmzNm7t3zbvTnr5cYJsZmZWgOooFgCvnrYLlYD7F68quFTWbrZPkF2DnBcnyGZmZgXoj6Aj3YX3n7IzXR3inqdeLLZQ1nY2pQR50oQuepevZ0tNkwsbPifIZmZmoywitqtB7u7s4MCpO3PPUysLLpm1m41bsgR55p6T6KsEC5atKbhEY4MTZDMzs1HWnwaurSbIAAdPn8Tvnl/Hsj9uLKpY1oaqNcgH7zkJgPlLVhdZnDGj9AlyfyV4fs0mtvb7JwMbmfWb+/jGHU9t+7ZtNlyrX9rCPb2u6bPh60sJcufL+TGv2293ero6uOKW3xVUKmtH1TbIU3buYcouPTywxO3Y89BVdAGG8uveldyy4Hn+uHELl57x2qKLY23sX379NM+t2cTUXXr4izftX3RxrE1VIrjy1kUALLx0Njv1dBZcImtH22qQO17OkKfs0sNxB0/jxw8tY8GytXz+tNdw0n/as6giWpuoVvp0d3UwY+rOzH9mNZv7+pnQ5c+mkWiqBlnSbEmLJPVKuqjB+gmSbkrr75M0o2bdxWn5Ikmn7GgBf7/6JQBue/wFImJHd7cCtCJehjrmULb2V1i+djMA8xa7E0y7KGMsPfPiS9vm73xyxQ6fkxWjbLH04vot1WNst/wth03nwKk707t8PR+b8yDX3N5Lrx8gUhpliyNg2/jZ3Z3idfvtzqoNW/jM/32EzX3+tXQkhqxBltQJXAO8HVgKPCBpbkQ8XrPZecDqiDhE0lnAFcCZkg4HzgKOAF4F3Cbp0Iho6l17fNlaFixbC8CyNZt4x9V3c8oRe/O6/XZj1526WLuxj9fvvzuTJ3bRKdHRIforQSWC7s7tc/9Kg2/rVRHxig+pdrNxSz8TuzsKP49WxEvaZ6hjDmrOvGfoj6BT4t8fXsZr992NNxywO/vtsRMTujqZ2N1BT2fj61f7xazo69tqff0V+iNKUfNQxljasLmP3zz7cvu+j815kHOOPZCD95zErAOnsM9uE+nqFF0dHelfbYuZrEMWbO0PerpK37ptTCljLJ0/Z37D5d2dHXz0hFezesMW5j68jCtvXcRVv3yCPzt0OgdN24V9d9+JDsFbX7MX0ydPYPVLW5g0oYtKwNqNW+mvBNMnT0DK2jf3dHY0vO+VTURQCbY1p+yvBBs297FszSbuXfwi7zpyX/bebWKhZSxjHC1fu2nbL1rdnR0cutdkTj58L+Y+vIy7e1fyjj/Zh9OPfBV/fGkr++6+E3vvNpGIYEt/hb13nTim72mVStCfOsNu2trPH/64kaWrN3Ly4Xs1dd7NNLE4GuiNiMUAkm4EzgBq37wzgC+m+R8AX1P26mcAN0bEZuBpSb3pePOaObnLfpq9xIkzp9PdJR54ehVfvu2JAbfv7BACguybVKWSPZ+8EtDVIfoqkSWRZBemen02bu2nu7ODzgYXrH5RGUMpgJe29LP3rhO5629PesWXg1HWinihiWMOaOX6zVz6k2zTj5xwELcueIG//8krd+1OiU1V9b3vrwR9lew/WU9XB9V7TX0ctbtKBJu2VpDgQ8ccWIYmTaWLpdsWvsCDz6xmyi49fPDNBzBn3jPcMO+ZQffp6hBdndp2bSNgQlcHXW2QtOThm2fP4k8PmVZ0MUoXSx8+bgZ/+8NH2dDgyWcdElMnTeDDxx3E6pe2cOuC57nziRXc+cSKbW2Xv/jvzdUP9HR20NkhgiAiu19Q82NsR8f2HQUhi9HB1H/mVbePdODa/aXss1LKPmMqkSUuHR2iU2Jrf4UgS5BjkNe+8f5nmfvXx7PrxO6hTrmVShdHe+46kZ6uDrb0VehJ9/63HLYn++2xMw8sWcV373+WOfc2/ozq6Wqc92TlGPg1Kw3eJDXIjnb03jicuKtsi5usorM2pqrlrD/uuX86gy+efsSQ5WkmQd4X+H3N30uBNw+0TUT0SVoDTE3L763bd9/6F5B0PnB++nO9pEU1q6fNAfeGacIzQM/nt/05DbglImaPcjFaFS9DHRMYOpY+d0XbxNI0Coz7v09TDcfSy7E07RlY+dsdO5eiFRZPx20fSEXEERQYS0N8JvFtmPbt9rnHFRZHzwC7fXa7Rf5Mqru/nd8+97dGRi22/nuaajSMpWYS5EbfAerz/IG2aWZfIuJa4NqGLy7Nj4hZQxXStpeuWxE3olbES6Mq8YbfNcdKLLVTWVuolLHUju9NO5Y5Z4XF0mCfSdBe7007lbVFSvmZBO3/3pSx/M38Fr8UqO3yvx+wbKBtJHUBuwGrmtzXxpZWxIvjaHxyLFleHEuWB8fRONJMgvwAMFPSQZJ6yBqZz63bZi5wTpp/H/CryHo2zQXOSr06DwJmAvfnU3QrqVbESzPHtLHHsWR5cSxZHhxH48iQTSxSG5oLgVuBTuD6iFgg6VJgfkTMBa4D5qRG56vI3mDSdt8na2zeB1zQ7AgWNQb8acoGVch1a1W8NDrmMIrXTrHUTmVtiRLHUju+N+1Y5tyUOJagvd6bdipr7hxHLVW68stjC5uZmZmZvcyDcZqZmZmZ1XCCbGZmZmZWo7QJskb46MWxRtL+km6XtFDSAkkfT8unSPqlpCfTv3uk5ZJ0dbp+j0g6quZY56Ttn5R0zkCvOVaUIZYGef++KOkPkh5K02k1+7TksaQ2fGW69pKWSHo0xc38tMyfB22iDLHk+8rYUIZYSuW4XtJySY/VLMstliS9MX3m9aZ9W/u0pYgo3UTWUP0p4NVAD/AwcHjR5Sr4muwDHJXmJwNPAIcD/wRclJZfBFyR5k8Dfk429uIxwH1p+RRgcfp3jzS/R9HnN9ZjaZD374vApxtsf3gq6wTgoHQOnWU5n/E4le3aA0uAaXXL/HnQBlNZYsn3lfafyhJLqSwnAkcBj9Usyy2WyEb9ODbt83Pg1FaeT1lrkLc9zjEitgDVRy+OWxHxXET8Js2vAxaSPYXnDOCGtNkNwLvS/BnAtyNzL7C7pH2AU4BfRsSqiFgN/BIo4oEio6UUsTTI+zeQbY8ljYingepjSUtxPuNUO1x7fx60h1LEku8rY0IpYgkgIu4kG7mjVi6xlNbtGhHzIsuWv11zrJYoa4Lc6HGOgyUT44qkGcAbgPuAvSLiOcg+7IA902YDXcPxdm1Ld7517x/AheknpuurPz/h96+MynbtA/iFpAeVPYIW/HnQLkp33X1faVtlv/Z5xdK+ab5+ecuUNUFu6hHV45GkScAPgU9ExNrBNm2wrOnHf48hpTrfBu/fN4CDgSOB54D/Wd20we7j8f0rk7Jd++Mi4ijgVOACSScOsq3jqVxKdd19X2lr7XrtdzSWRv08y5og+9GLDUjqJvsQ+05E/CgtfiH99ED6d3la7sdaZkpzvo3ev4h4ISL6I6ICfJPs5zLw+1dGpbr2EbEs/bscuJksdvx50B5Kc919X2l7Zb/2ecXS0jRfv7xlypog+9GLdVJvzeuAhRFxVc2q2sdangP8uGb52amn6DHAmvTzxq3AyZL2SD/nn5yWjVWliKWB3r/qB0fybqDa+9ePJS2f0lx7SbtImlydJ/t//Bj+PGgXpYgl31fGhFLE0iByiaW0bp2kY1Lcnl1zrNZoZQ/AkUxkPRyfIOud+fmiy1P0BBxP9nPCI8BDaToNmAr8B/Bk+ndK2l7ANen6PQrMqjnWfyHr9NULfLjocxsPsTTI+zcnvT+PkH1g7FOzz+dTmRdR01u3DOczXqeyXHuyHusPp2lBtSz+PGifqQyx5PvK2JjKEEupHN8jayq4lazG97w8YwmYRVYR8BTwNdLToFs1+VHTZmZmZmY1ytrEwszMzMysEE6QzczMzMxqOEE2MzMzM6vhBNnMzMzMrIYTZDMzMzOzGuM2QZa0vmb+NElPSjqgyDJZe3IsWR4cR5YXx5LlZTzHUlfRBSiapLcCXwVOjohniy6PtS/HkuXBcWR5cSxZXsZjLI3bGmQASSeQPd73HRHxVFr2LUlXS7pH0mJJ70vLJelKSY9JelTSmWn51yWdnuZvlnR9mj9P0mWSZkhaKOmbkhZI+oWknYo5Y2sVx5LlwXFkeXEsWV7GayyN5wR5AtljCt8VEb+rW7cP2ROG3gl8KS17D3Ak8HrgbcCVyh4TfCdwQtpmX+DwNH88cFeanwlcExFHAH8E3pv72ViRHEuWB8eR5cWxZHkZt7E0nhPkrcA9ZI9CrPdvEVGJiMeBvdKy44HvRUR/RLwA/D/gTWRv7AmSDgceB15IwXBsOj7A0xHxUJp/EJjRihOywjiWLA+OI8uLY8nyMm5jaTwnyBXgL4A3Sfpc3brNNfOq+3c7EfEHYA9gNtk3pLvScddHxLoGx+vHbb/HGseS5cFxZHlxLFlexm0sjecEmYh4ieyngQ9KavTtqNadwJmSOiVNB04E7k/r5gGf4OU3/dO8/JOBjQOOJcuD48jy4liyvIzXWBr33/QiYpWk2cCdklYOsunNZD8FPAwE8NmIeD6tu4usZ2evpGeAKZT4TbfWcCxZHhxHlhfHkuVlPMaSIqLoMpiZmZmZlca4bmJhZmZmZlbPCbKZmZmZWQ0nyGZmZmZmNZwgm5mZmZnVcIJsZmZmZlbDCbKZmZmZWQ0nyGZmZmZmNf4/uDL92Qy4ysIAAAAASUVORK5CYII=\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"g = sns.FacetGrid(profiles.reset_index(), col='Set', sharex=False, sharey=False, height=2)\n",
"g.map(sns.distplot, 'Known')\n",
"plt.savefig(fig_dir / 'profile-size-known.pdf')"
]
},
{
"cell_type": "markdown",
"metadata": {
"papermill": {
"duration": 0.065854,
"end_time": "2020-05-08T03:36:41.701767",
"exception": false,
"start_time": "2020-05-08T03:36:41.635913",
"status": "completed"
},
"tags": []
},
"source": [
"For the paper, we want to make some changes - we're going to show this on a scatter plot."
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"papermill": {
"duration": 0.143,
"end_time": "2020-05-08T03:36:41.907513",
"exception": false,
"start_time": "2020-05-08T03:36:41.764513",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Set | \n",
" user | \n",
" Type | \n",
" Size | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" AZ | \n",
" 529 | \n",
" All | \n",
" 8 | \n",
"
\n",
" \n",
" 1 | \n",
" AZ | \n",
" 1723 | \n",
" All | \n",
" 25 | \n",
"
\n",
" \n",
" 2 | \n",
" AZ | \n",
" 1810 | \n",
" All | \n",
" 14 | \n",
"
\n",
" \n",
" 3 | \n",
" AZ | \n",
" 2781 | \n",
" All | \n",
" 8 | \n",
"
\n",
" \n",
" 4 | \n",
" AZ | \n",
" 2863 | \n",
" All | \n",
" 6 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Set user Type Size\n",
"0 AZ 529 All 8\n",
"1 AZ 1723 All 25\n",
"2 AZ 1810 All 14\n",
"3 AZ 2781 All 8\n",
"4 AZ 2863 All 6"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"up_sizes = profiles[['count', 'Known']].reset_index().melt(id_vars=['Set', 'user'], var_name='Type', value_name='Size')\n",
"up_sizes['Type'] = up_sizes['Type'].astype('category').cat.rename_categories({\n",
" 'count': 'All',\n",
" 'Known': 'Known-Gender'\n",
"})\n",
"up_sizes.head()"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"papermill": {
"duration": 0.207822,
"end_time": "2020-05-08T03:36:42.154324",
"exception": false,
"start_time": "2020-05-08T03:36:41.946502",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Set | \n",
" Type | \n",
" Size | \n",
" Users | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" AZ | \n",
" Known-Gender | \n",
" 5 | \n",
" 1300 | \n",
"
\n",
" \n",
" 1 | \n",
" AZ | \n",
" Known-Gender | \n",
" 6 | \n",
" 829 | \n",
"
\n",
" \n",
" 2 | \n",
" AZ | \n",
" Known-Gender | \n",
" 7 | \n",
" 572 | \n",
"
\n",
" \n",
" 3 | \n",
" AZ | \n",
" Known-Gender | \n",
" 8 | \n",
" 416 | \n",
"
\n",
" \n",
" 4 | \n",
" AZ | \n",
" Known-Gender | \n",
" 9 | \n",
" 284 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Set Type Size Users\n",
"0 AZ Known-Gender 5 1300\n",
"1 AZ Known-Gender 6 829\n",
"2 AZ Known-Gender 7 572\n",
"3 AZ Known-Gender 8 416\n",
"4 AZ Known-Gender 9 284"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"size_counts = up_sizes.groupby(['Set', 'Type', 'Size'])['user'].count().reset_index(name='Users')\n",
"size_counts = size_counts[size_counts['Users'] > 0]\n",
"size_counts.head()"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {
"papermill": {
"duration": 0.099594,
"end_time": "2020-05-08T03:36:42.324672",
"exception": false,
"start_time": "2020-05-08T03:36:42.225078",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"data": {
"text/plain": [
"count 4202.000000\n",
"mean 11.899096\n",
"std 46.498697\n",
"min 1.000000\n",
"25% 1.000000\n",
"50% 2.000000\n",
"75% 7.000000\n",
"max 1300.000000\n",
"Name: Users, dtype: float64"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"size_counts['Users'].describe()"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"papermill": {
"duration": 17.040578,
"end_time": "2020-05-08T03:36:59.435841",
"exception": false,
"start_time": "2020-05-08T03:36:42.395263",
"status": "completed"
},
"tags": []
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/MICHAELEKSTRAND/anaconda3/envs/bookfair/lib/python3.7/site-packages/plotnine/ggplot.py:729: PlotnineWarning: Saving 7 x 3.2 in image.\n",
" from_inches(height, units), units), PlotnineWarning)\n",
"/home/MICHAELEKSTRAND/anaconda3/envs/bookfair/lib/python3.7/site-packages/plotnine/ggplot.py:730: PlotnineWarning: Filename: figures/ProfileData/profile-size.pdf\n",
" warn('Filename: {}'.format(filename), PlotnineWarning)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmgAAAHOCAYAAAAsQvUnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdeZgcVfXw8e+ZmSSQQBI2WS6LEBHZd2URBBUEwVYRBBEhoLggm/ADXgQEFUQEAZFNdhAUQbZmB0GQfZEtyCYBWYo17AkkmSTn/eNUMz2d7pnunuqqnq7zeZ55kqnqrjnTd6rq1l3OFVXFOeecc861j66sA3DOOeecc/15Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs14Bc0555xzrs30ZB2Aaz8hhKuALYDFoyh6O942tcpLRwI9URR5RX8IQgi3AesBswAFngb2j6Lo9hDCb4ENgC9EUaTx6w8AdgLWiaKot8rxjgAOAWZU7FouiqI3WvV7dKIWlc1noijaIZ3foHOEELYF9gNWw/62XwAuBk4CdgDOBj4C5gARcHQURRcMcLyJZe8pt2kURQ8kHX8na1HZ/DiKovVaG3l78xur6yeE8AlgS2AasGNpexRF85V/AUsCLwK/zCbSjrNv/LmOA/4EXBlCGAEcDiwM7AUQQlgB+AWwa7UKQJnLKsvMK2dNS7psXINCCPsCp2E3/CWiKFoIqwgvAywVv+yBsnI6EDg7LpOBPFDlPPHKWQNaWDa55xU0V2kn4L/ACcCu1V4QQhDgQqw14Vfphdb5oiiaA1wEjMdaMGcAuwFHhhCWB84F/hBF0UMZhplLXjbZCCGMA44E9oii6OIoit4HiKLoiSiKfhpF0X/LXx9FkUZRdDXwLtai41rEy6a1vILmKu2K3YQuBNYKIVQ7iX4BrAjsVOracckIIXQDu2Ctk68ARFF0L9ZyczcwFq8UZ8LLJjPrA/MAV9bz4hBCVwhhG2Ah4NlWBua8bFrJx6C5j4UQ1gFWBi6KouiFEMJdWIXtZ2Wv2QI4ANgwiqJ3som0Ix0fj2maFxunsUsURbPK9v8L+D/gj1EUzazjeNuEEN4t+/71KIq8S6E5SZeNa8zCwJTybuMQwvVY5WAk8KN48zrx3/xorPFh3yiKHh7k2OtUnCcAi8ato25wrSyb3PMWNFduV+DOKIpeiL//M7BTPN6GEMKyWOvaHlEUPZpRjJ1qvyiKxmOVgM2AM0MIawOEEMYCp2DdzvuHEJYpvSmEcHoIYWr8dXrZ8S6Pomh82ZdXzpqXdNm4xkwBFi5dhwCiKNoyLpPHgO5484PxtrHA6cCmpdeHEH5eVhbXlx37wYrzZLxXzhrSyrLJPW9BcwCEEEYB3wFGhRBeizf3YE3RW8cnzmXAxQPNvnFDE49zujOE8Ax2Efs38HvgoSiK9gsh9ABnAF+JX/9j4MdZxZsnXjaZuQebGfgN4NLBXhxF0fQQwv7AMyGEr0dRdFUURb8BftPiOPPIy6aFvAXNlXwDq5CtAqwRf60C/A1rWTsNOxF/VusALhkhhPWAlYDHQwibAd8CfhLvPhj4dAih6gQO11peNumLoug9bNzrKSGE7UMIY0MIEs8CXKTGe2YAxwNHxJOaXAt42bSWt6C5kl2BP0dR9Hz5xhDCCcC98bczgLdDCJXvXSmKohdbH2JHOzGEcFz8/9eAQ4G7gMexLrZXAaIomhZC2B24JIRwQ2l7Fd+qkrtu/SiKJrUi+A6XdNm4BkVR9PsQQoQ9IJ6D5dR6ATgVa7nZrsrbzgIOA75O7UHs61Y5T3aIouiaRALPgRaWTe6Jqk/Cc84555xrJ97F6ZxzzjnXZryC5pxzzjnXZryC5pxzzjnXZryC5pxzzjnXZryC5pxzzjnXZjzNRoomT568NLY0hkvXlAkTJlRNA+Jlkgkvj/bjZdJevDzaT80yaRVPs5GS+KR6EluLzKXrQ2DFypPLyyQzXh7tx8ukvXh5tJ+qZdJK3oKWnoWxk2on7ARz6VgRuBD7/CtPLC+T9Hl5tB8vk/bi5dF+BiqTlvEKWvqenDBhwkNZB5EXkydPrudlXiYp8fJoP14m7cXLo/3UWSaJ80kCzjnnnHNtxlvQMiAigi24vCDwhKq+lXFIzjnnnGsj3oKWspdffnkEUMQWWv4X8JKIbJ1tVM4555xrJ15BS9nEiRN/BGxetmle4FIRCRmF5Jxzzrk24xW0lL3++uvrAyMrNvcAa2cQjnPOOefakFfQUtbT0zMVqEw+1wVMzSAc55xzzrWhYTNJoFAo7AlMBFYFrigWizuU7VsFOAtYDXgO+EmxWLyj4r0HA2OB64Ddi8Xi+/G+8cAZwJbA+8BRxWLx1Fb9HhtssMEFN95441qAxJtmAk8Dd7bqZzrnnHNueBlOLWivAEcCZ5ZvLBQKpUH3VwALAL8FrioUCgvE+zcDDge2AgIwCvhj2SFOxiqqSwBbA78uFAqbtuqXOOWUU+4Cvgk8giW8uwzYVFVntupnOuecc254GTYtaMVi8XKAQqGwBv3XIdsEy6x8bLFYnANcWCgUfgZsA5yNtbqdWywWH4nffxjwQKFQ+AnWirUdsGaxWPwAeLhQKJwH7Ab8s1W/i6oWsUqlc84559xchlMLWi2rAJPiylnJI/H20v5Hy/Y9DnQDywOfBqRYLD5R473OOeecc6nrhArafMB7FdveBeavtr9YLCo21mz+eN/7A7zXOeeccy51w6aLcwBTscH/5cYBHwywf2y8v2uQ97qcEJEuYHHgI1V9O+t4nHPO5VsntKA9DqxaKBTKf5c14u2l/auXdsQzPmcD/wWeAbRQKKxY470tJSKLikhRRKaIyLMismMaP9f1d+aZZy4JPAW8DLwlIpeLyOiMw3LOOZdjw6YFrVAo9GDx9gBdhUJhHqyidRswHdi/UCicBHwLWA6b1QlwHnBRoVC4CHge+DXwt2Kx+GF83L9jMzd3jd83Efh2q38fEZknjn05LHHtQsCFIjJdVS9v9c93pre3lxNOOOFUYNGyzVsBJwI/zCYq55xzeTecWtAOBT4CDsFmXn4EnFksFnuBArAtNn7sEOAbxWLxbYBisXgzVim7AXgV6AX2KjvuT7HEsa8C1wO/KBaLt6bw+2yATVIoX1VAgH1T+Nku9sILLzBz5sxA/4eVkdgsYOeccy4Tw6YFrVgsHgEcUWPfJOBzA7z3j/TPfVa+712swpe2eYE5zF1Jni+DWHKrp6fmKdCbZhzOOedcueHUgtZpHsRWESjXC1yTQSy5tfTSSzNmzJjH6V8WvdjqEi5DIjKfiCzv4wGdc3nkFbSMqOrrWNdseZqPy7HuWJeSrq4ujj/++H2Af2CVtA+A44BfZRpYzonIT7EhC88A74jIzhmH5JxzqRo2XZydSFVvEZGAJc19R1X/l3FIufSlL33pXVXdKus4nNlnn33Ww4YklNarHQmcKyLPqOq92UXmnHPp8Ra0jKnqVFV92CtnzpnHHntsQ2yGdrle4CsZhOOcc5nwFrQ2ISLLAesD04CbVXVaxiHlkoisAWwPjACuUdXbso0of7q6unqxmdWVKsdsOudcx/IWtDYgIttgiVLPBi4FHhORJbKNKn9EZHPgAWA/YG/gVhH5frZR5c9XvvKV6+P/llfSRgG7i8gGGYTknHOp8wpaxkRkQeAvWIvNKKxVcyng9CzjyqmzgG5szNMIbAzUqXFSYZeSAw888L/AFtjKDuWWAW4RkU+lH5VzzqXLK2jZ+zRWMSs3Alg3g1hyS0QECPQNTC8ZCXwi/YjyTVVvBf5G/3x0XVj5bJtJUM45lyKvoGXvzSrbFHgj7UDyTFUVeBFLHlxuBvB6+hE5+loxq213zrmO5hW07D0H/Jm+lgKNvw7OLKL8mgjMwiplM7DK2m6qOiPLoHLsWqzLudwI4MYMYnHOuVR5BS1jccvNrtgaorcCVwFfVtXrMg0sh1T1dmB14PdYaxrAeSJyhohUdkO7FlPVm7G1cmfFm2YAO6vq/dlF5Zxz6fA0G21AVWcDx8ZfpSVuulS1srvNtd7T2EL2y2APMF3ALlherp9kGFcuqeppInI+sBjwqqp+lHVMzjmXBm9BayMisr6IvIwtNzQ1Xu7GpWsRYBNsckDJSMCXGsqIqn6oqs8B40RkORHxB0vnXMfzClqbEJElgZuBxeNN8wInxTnSXHoqxzyV+LmSEREZLSJF4FVgMjBZRFbKOCznnGspv+m0j82wLufyMukCdsgmnNx6DXiI/ukdZgJXZBOOA06k/zJPAbjRxwU65zqZV9CcKxNP2tgaeKRs8/XAD7OJyAEF+nc5dwNLAitkE45zzrWej+VoHzdjs9VG0FdxngNcnFlEOaWqrwKfFZEFgFmq+kHWMeVcb43tvjanc65jeQtam1DVl7FuzlfjTR8Be6nq5dlFlXs9wDEicruInBWPE3TpO5m+VBtgM2qfBZ7PJhznnGs9b0FrI6p6D7CkiMwHTIu721wGRGQctnD64lj32vpAQURWU9XXMg0uf47FrlWHAqV1UT8J3CEim6jqh1kF5pxzreItaG1IVad65SxzO9NXOQPreh4H/DiziHIqzgd4P31LPHVjFbbVgf/LKi7nnGslb0FrQyLSDayFVQgeU1VflzN9i2BLbpXrjre79K2MjTmbt2zbSGC1bMJxzrnW8ha0NhN3b96KtRjcCLwgIltnG1UuPcTcDzBzgIcziMVBxNyLpPcCL2UQi3POtZxX0NrPMcB68f+7sDE3l4rIotmFlEtXAWdgrWjTscrZlcA5WQaVY1diYwJLMzd7gfeIl0dzzrlO412c7eeL9M/5RPz96sBN6YeTT/EYwD1E5EJsdu1yWBb7FYH/ZBlbHqlqr4h8Edgf69aMgGPjlCjOOddxvAWt/bxbZVsX1lrg0heAXwDbAz8HHhaRL2cbUj6p6nTgt8BjwMbA30Vkx2yjcs651vAWtPZzFNa9Vqo8z8S6dh7MLKKciidrnIeVRWlZIQXOxypuLn2nARPpG4+2nojMp6pnZBeSc84lz1vQ2oyqXgN8HbgPeAYb87SFqs7ONLB8WgQYXbFNgCVExB9uUiYiCwO703+yQBfwy2wics651vGbTBuKK2nXZB2HYwo2QWCeKvsOBn6dbji5N77G9gVSjcI556oIIdyKPcTXFEXRpvUez1vQnKtBVWdRe5H0X4nIymnG43gReIu589ONEpH9MojHOefKXQz8Nf66ElgTGz/+17KvunkFrU2JyBIiUhSRV0TkYRHZLOuY8khV/wy8XGP3ZSJSmZvLtYiqzgQKNXb/TkRG1djnnHMtF0XRGVEUnQFcAGwNPAd8AXi+bF/dvILWhkRkfuAuYAtsuaHVgRtEZKNMA8uv/2B50CpNAL6Rcix5d2+N7d3At9IMxDnnKoUQRgPXYUNjNgK+C1wSQtiw0WN5Ba09bYHNEiy1zghWVseLyID9264l9qN6Ba0Xq0C7lMTrcv63xu7TvRXNOZeVEMIYrHLWDWwRRdG0KIquwyY3FUMIazVyPK+gtacxQLVZm2tjiTpdilT1CeBk5i6TUcCk9CPKvVq5z+YHVkozEOecK3Md9jC/RRRF00oboyj6O/agf30jB/MKWnu6h+ozbAU4POVYnPk5lotuFvARdhKepKr/zDSqHFLVfwO1VhD4fZqxOOdcmVnAVlEUfVS5I4qi82lw5r+n2WhDqvq0iPwU+FOV3aNFpNvzoqVLVT+KxwBuC3wN6968X0S64m43l659gb9V2b6piIxX1WorcjjnXCutDDwRwtx5zKMoWjaKopMbOZhX0NrXmVhr2WL0tXTOAh71yllmRgAHYCehYl1t24vINl5JS92lwEVUv4Z5z4BzLgs7V3y/KLA3DabXKPEKWptSVRWRrYGbgQWx7s0eYKqILKyqUzINMJ/2wSpn5YvZb4XNHrw0k4hyKj4/rsTSbpTKoxd4UFXfzi4y51xeRVF0U+W2EMItQBE4vtHj+ZNmG1PVh4Gz6T+DcH3gGhHxskvfp+m/zBBYq+anM4jFwW7AjWXfPwB8M6NYnHOumneB5Zt5o7egtb+dsSm7JSOBzwHLAc9mElF+vYC10pS3oPXE213KVPUDEZkI/AX4MrAu8GsR2UtVZ2QanHMud0IIP6rYNBobt3xHM8fzClr7q5Wp3ssufScCuwDL0LfeWg+Wn+4uVX0+s8hyKM4JeDnWqtwdf+2CtWrukWFozrl82r7i+6nYMKUTmjmY3+Tb3+XAROauqP0EGxPlUqKq74rIxcCB9F8QdxHgPhGZoKofZBNdLi2CLaNSbiTW6uwVtISJyCeBr2P3jZtU1XMAOlcmiqIvJnk8r6C1v32x5YQWqdi+t4icoqrPZBBTni1H9bGbiwA/Ao5LN5xc625wu2vSEUccsTJwGnbPUOC3IrKdql6ZbWTOtZcQwvLAt7HVgCLgkiiKaq1+MiAfaN7mVPVDrJm00hxg6ZTDcTbeTGvs+52IfDbNYHLuNeAhbFxgSS/W6uwSdOmllx4JzIutnjEPVlG7SERGDvhG53IkhLAVdk1aE0tovhbwSAjhq80cr2Na0AqFwp5YV+CqwBXFYnGHsn2rAGcBq2Gry/+kWCzeUfHeg4Gx2FINuxeLxffTi35QTwBL0b+8BJicTTi5dizWUja2yj4BLsRndaaiLBXNbfR95iOAICLzqupc2bxd41SVGTNmLMncD/SjsTyNL6YflXNt6WhgmyiKbi5tCCFsARyD1S0a0kktaK8AR2IJXj9WKBRGYDlIrgAWAH4LXFUoFBaI92+GJYTdCmuSHAX8Mb2w67I3NlV3Jn0pN/6BNZ+6FKnqm8CnsIHo1XirZrq6gE9WbPsCcHX6oXQmEaGnp6dabrnZwJtpx+NcG1sKuLVi2000eV/omApasVi8vFgsXglUJnDdBHvSO7ZYLM4oFosXAs8D28T7JwLnFovFR+JWs8OA7QuFwuh0Ih+cqj6HtdqMoG9w+ibApfFMNpeiuJJ2So3dXmlO19pUH3P2JRHZNO1gOtXGG2/8G+zhcHb8rwL71WqlFJEuEfk/EXlMRB4Vkf38WuVy4GlsJnm53YAnmzlYx3RxDmAVYFKxWCxP9vpIvL20v7zp8XHsgr888GgqEdZnv/jf0kVuBJZFfXXs93HpOgAbX7BR/L1iN64fZxZRPr1P9QqaApsBvph9As4444zbP/WpT20E7ITdN65W1YFaKX8D7E/fPeYYrAfjsJYG6ly29gauCyHsBfwPa91fCtiymYN1TAvaAOYD3qvY9i4wf7X9xWJRsYv+/LSXRemf2gGsQrBwBrHknqr2qurGWAvsC8A07OlpZpZx5dBdwKtVtgtzp+BwQ6Cqd6vqHqr6w4EqZyLSA/wf/RsAeoCDROQ7IrJSq2N1LgtRFN2PDYE5Frg7/ndCFEUPNHO8PLSgTWXuAd3jgA8G2D+2bH+7uA/rxy6fNTUbm0DgMiAi47Bxj4tiLZqfAW4RkU1U9c5Mg8sJVe0Vka9hyzxVPsBsEOem88k0LSQiGwJfxcZlXoJdP6u1ao7AFrhHRH6mqn9ILUjnUhJF0bvY6iYAhBBGhBA+H0VRw/eEPLSgPQ6sWigUyn/XNeLtpf2rl3bEMz5nA03lLWmhfbClncpTPPyX/ut0unRtTV/lDOx86gZuEJHlMosqZ1T138AbNXZ/Kc1Y8kJEFhSRFURkN+BfWIvZQcDDDJzmROKvE0RkjdZH6ly6ShWyEMKhIYR/AO8A5zRzrI5pQSsUCj3Y79MDdBUKhXmwitZtwHRg/0KhcBLwLSzZ6BXxW88DLioUChdhkwd+DfytWCx+mOovMAhVfUtELgN+Tt/T6fLAjSKylqrOzi663Jof+xurXOVhNPB3bIyaS8fjWGW50l7AGSnH0rHigf7HYGMwy5Va9hVLvTHooYDbRWRHVb02wRCdy0wI4WZgA2zYxW1YxWznKIpeaeZ4ndSCdiiWGO4QYLv4/2cWi8VebDD9ttjYs0OAbxSLxbcBisXizVil7AbsQ+3FLurtaCL9uw5GYLndVsgkGncX/bucSwRYw5N4purEGtuXSjWKzrc7fROWqmlkpuZY4BoR+c7QQnKubSwIzMCGJN0PPNhs5Qw6qAWtWCweARxRY98k4HMDvPePtF/us2pqlZcvbZMBVZ0kIrsDZ1fZPZP+Ge5da92DdfdXPnSOFZEvqOrtGcTUiQokf705Q0Q+wCrTTwG3qWqt1TrcIERkFDb57W3/HNMVRdHaIYTxwMZYKqyfhBAWAe6MomiHAd9cRSe1oOXBJfSfJTgL65Z9OptwnKqeA/wJ6+osmQX81i+O6VHVt6iewkGAK0WkshvaNWc6tZc6a9Z8wJXA77EE3Kd5zrTGxbnnjsN6j6YAz4nIqhmHlTvxJIE7gNuxbs4ZWLdnw7yCNrwchFXSSiYDm6uqp3bI1mnYbNrp2Ozf84BfZRmQ62c8tkqIG7ozB39JU7qxtT67gB9gs0LdIERkTRHZVUS2wsa97kdfN/NS2KzycZkFmDMhhBNCCI9gOVS3AyYBm0VRlO+VBPJAVWeo6vewwemfAFZU1WczDivX4ifU+7CEx/NgZfMD4HpvtUndIgPsq7ZUkWuQqt4IbE9rl3jqpWxmvatORA4C/g2cii1t9k36jwHsBhYC1k0/utz6PrAMcDNwI3BTFEXPNXswr6ANQ6o6VVXf9C60tnAI/ZfgKvki7TvZpFM9RP+u5pKrVfX9tIPpVKp6qap+AtgCuJfBc0Yq1qKwB9aNOZgu4LUhBdnhjj766BWwhbkFezCs1SUseCqmNC2Atf4+i6268WQI4dkQwlnNHKxjJgnkjYgsA6yPjTe4RVWnZhxSXi1O9QedHjzNRtouBL6MXRjByuVu4OuZRdTB4ta0G0VkX+B3zJ1uBmw85gOqugGAiIzBlsOp1bo8E7u5XZx8xJ3jscceWwEb2zTPIC99GbhPRFYG1sTGpt2iqj6BqQWiKJqNTVi6Bzg6hNCNtWA2tS6wV9CGIRHZCzihbNNLIrKRqr6cVUw5dg82ALTyXJpF7eSprgVUVUVkIjardkesu/kJYAUReUZVvSWhNU7GZqx9Dfu7H4m12szBWsJKFWZUdZqI7IKtKFBt6borgd1Vta3yULabBRZY4G1qV3JLpmMt+bsBf8Aqvz3A/SKymapOAxCRbuxB8wNVrVwW0Q1BXGG7N/5qmFfQhhkR2QY4qWLzUsBZWJeDS9evsJtTZRqXXuyi6NK3B5aQuiv+OhJ4QUS2VNUnM42sA6nqLBH5JrARsAQ2Pm0cdg7cpqofVLz+ryJyLHNP3OjC8lWehY3hcTX88pe/vPemm256AFiH6vfx97AHx1HYdUji/xO/50wR+TGwLHANsCSAiPwV2FlVZ7X2N+hMIYRbqdHdHEVRw61oXkEbfqplRe8G1ks7EAeq+qGIfB5b4eGXZbtGYOPTfphJYPm1GXaTr8zVtTRwk4gsr6rT0w+rs8XjYf/VwFsq1z8u6cLK6T/AHqrayDFzY5FFFpmFLWN2F7Z0YbmZwKmq+oSIfDf+flTZ/hHAd+h7oB9ftu872DiqLVsRdw5Uds0vhiWYv7CZg3kFbRiJZwUuVGP3WBFZTVUfSzMm93ELwtexwdClp6ceYHcROcYX607VBOyGNG/FdsFaCVbCJhO4bE3FuqBrWRn4Z5xkuOFFpvMgfji8HvubLl+1RIDvicj3gVeo3RU6nuqtPVuIyIaqeleiAedAFEVzNaCEEP6GpWJqmM/iHEbigZ21prcrPmswS4tT/WLXVIJC17Tn6d9aUMm7bjImIutQ33qdXXg+wcGcAkyj7++6F6uQBSwV0yrUTiw8UDLg7ycVoGMy1q3cMG9BG352BYrMXbnuonbrmmu9p7FKWqXTReQyH/ScmhuB9+nfbVPyITZpwGVrLeqbgQiwsYgsoqqtzLs2bKlqJCJrYRXZ5bC1IFeg7/7QQ3MrP0wUkbO9Fa0xIYSvVGwajU1YerSZ43kL2jCjqtcCO9fYvaAv0J2ZP9fYPpr+M25dC8VjoWqlEPiQwWe+udZ7i/obB7qx2aGuBlX9n6rurKqfxyoClff1ZpbNEuB4X3KrYadVfP0Sa+HcsZmDeQvaMKSqF4nIAsy9wPv62AyoWhU41zr/GWDf7iJyos8gTE0ELMzcN6bxwDEiso8nec7UtVhL5mfoP3aqloVbG05H+Qc2SSaJe/u62OSnoxI4Vi5EUbRcksfzFrTh6xQsSW25kdjg0IEG37rWuB+oVQETmnyCck35P6p36/Rg4zT/IiL+cJqReBbtxsCfsGXSbh3kLTe0PKjOcQ52b6im0TyAAuw/tHDyK4Rw8FCP4RW04Uuo/fS5apqBuI+71tbAxj9Vs2uK4eSaqt6CpTip1Uq2PT74PFOq+p6q7q2q6zHwSg+T8YTPdVOzL1bxrdTM/X5+7+Zs2s+HegCvoA1TcVb0WgMP/aknA6o6E/h/NXYHEdk2zXhy7vfAgzX2CbCn33jaQ7xMXa1VUJYBnhSR1QBEZB4R+XQ8xMNVISJfBj6b0OF6gH+JiNcVGjfkYRT+oQ9vp1P9j2CZtANxHytSfcFugF+IyAppBpNXcUqa/QZ4yfzAoimF4wa3MdVbn3uwVQluFJEdsbUknwbeEpGjvJJd1bdJoHJQZkMs2aprzB1DPYBX0Ia3aut7zaTJKb1u6FQ1Ak6tsXtl4DER+VKKIeXZFwbZP1C+NJciVX0eG6dZbZyUYHnTLgLGlG07CB86kAYBjhWRcVkHMpxEUbTVUI/hFbTh7QlsCYnyp6UebMFhl51jsIpypS4szcNf0g0nt+Zn4IHRk0XkEG+FaRvz0VjLj6fgqO4ymkutMZBxwGEJH7MjhRAmhBAuDCG8GkLojf/9cwhh2UaP5RW0YUpERgE3ATvR/2QU4GIRWSqTwBxYEs5aswQF+ISIjKmx3yXnfwx8jevG8hR5Wpr28G8auycp4OuqVlDVG4EfYHn/ktINbC8ie/kM6NpCCBOwGf3jgIOBrbBxyZEhJNYAACAASURBVOOBB0IIDaXh8Ara8HUg8HnmflIS7GTaNPWIHHHF624GPrdmkuzF01W3G4O3yHRjY3ZcxlT1WeC/DbxFgFVEZOMWhTRsqeo5WAvyhtjygLXGxTZiSWzyzdU+aaCmXwN/jqLoa1EUnRdF0U1RFJ0fRdHXsGTmDc0e95rw8LUetdNsCI3nvHFDFF+0rgE+NchLn/BEqalYg/q6ehZpdSCubi8Anx7kNUpfua4M3CoiG6nqPS2NbJiJZ/rfLSLLA5/DltZS7AFxe2AXGm+kGQFsBmyJJRx2/X0JWLvGvoFmllflteDh61WqPxXNwVpnbk43HAesDmzC4JWCWikFXLLeqfN1PqmmDcQPOINN7IC5h3QI8LOWBNUB4pxzN6lqUVWvjrtAhzoO1jMFVDdvFEUfX99DCDuV/h9vn7eRg3kFbfg6FhvrNKti+7PApqr6evoh5d5Y6hvkfF2rA3GAzfIbrDwUuCSFWNzgRlDf0k+VurBFwl39xtJ8t2c3cICv+1zVqyGE8tmuHy/HGEIYT4NJl72Lc5hS1adFZF3gF9jYgAeAx4AlgJVF5HFVray8udZ6HFt+a/Qgr7s6hVhyT1XPE5H3gB9iN6T5gJWw694srOXlz9j6hS5jqjpDRB4BVqOxxoNeBl8uyvX37yG+fxlsHPSRCcTSSe4Ajg4hlFrlu8v2bQPc3sjBvAVtGFPVJ1R1B6xb4NPAGcChwLnAzSIyIsv48kZV38JOwg+pPQbwemA3EdkstcByTFWvUNUtgb2BVeh7KO0B3gMu9PGAbWVbGl/a6WqsR8HVSVVfoHbeuXoIsFFyEXWM3wGfwcb4bY9NGCsZHe+vm7egdYbvYgM3e+gr0w2wloNaC+e6FlDVG0VkGawVYH1sPbbRWDqAGcAXgS8DPSJyMfBdryCk4hbmfiBdEHuQOVRVf5NBTK6Cqk4WkQnAAViagsGSCc8BDo9XjnANUNW/i8gONNfFr9ikA1cmiqJnsGt8tX0nN3o8b0HrDJ9h7ieh7ni7S5mqTlHVW1X1KKxrbRFsCajR2A1nBPYE+h3gEk+U2loi8kksL1HV3cCRcaXAtQFV/VBVf0nfqgED6QLuimcqusZdDtxJ9cTaAxFgYxE5PfmQhrdgDgohnBJCOCCEsFizx/IKWmd4mblnDs4BxojIKhnE42KqOht4C1gHq5hV2gZrUXOtM1jaEwW28Ipy25mvzteNxtYldg2Kr09bAKfRXG7GH4nIrp4XzYQQ1sJSaayA5WFcC3g8hLBuM8fzD7UznAs8V7FtBJbnZpKInOI3n2zE4wAvA2plkO7FWzpb7cVB9ncBJwNX+rjNtvI+NixgMD3A6n6Na46qTsPWOW0oBUSZc4Ciz+oEbIzZ7lEU7QZ8FEXRd7D1Yn/fzMG8gtYZxmMzOSuVyvdH2Dg1l779GXi9wFHAjiIy2MxP1yRVfQY4kcHTCmyJzUxz7aPefI4LAVNEpNp10A1uVeqrDNeyOTZmMO/WwCaClbsGm6DUMK+gDXPxU+NVDN4dsGEK4bi5bc7gk3HWoixfjmuJ/YCHGDgv2gh8ibR28lOs+61eCwIPektaU16n+hCMeo2gxuD4nFH6UmuU/g7XwRLLN8wraMPfgsBg/duzgXdTiMXNrZ6L3khsLJprnW7sQjnQzXsONl7QtYdv03imgUWpbzUC19+NWEqIRicLlFtMRPL+2T+KtaIBdIcQzgJuwNJfNcwraMNfPcloZwPntzoQV1W9SWmTWMzY1bYugy/B1QXc5y0wbaPZ1BkHJBpFDsRJzTcH/k59q6FUMwH4p4jsnlhgw89+WH5FgAuA/wEbRFF0RTMH8wraMKeq72F93gM9+cwL3Cki66cTlStzFoNXonvj17nWmY/6bjy/xwY815PiwbXWWTSXSDXvrThNUdXp2HipZlvRurGHoFNFpN4ZuB0liqLHoih6Ov7/nlEUHVn6vhmeqLYz7ACcjWXhrmVB4HoRmRBnvHcpUNW3ReRVYKlqu7HK2f+AxUVkfVW9J834cuQR6luGC6wl4QQs0bPLiKr+VUQWwJYTGs/gLaAlY0TkIuAcVb2lZQF2pvsYesNND7bk4DNDD2d4CSGcW2tfFEW7Nno8b0HrAKr6PjZLc6AnH8FaEdZOJSgHgIiMovoM25JubJmunYA7fAmo1lDVN4FvUF/rwEig0NqIXD1U9VRVXRDL09WIHbFVIrZrQVgdS1WfwxJoD2VG5yzglWQiGnZerfhS7FrSTI45b0HrIN9k8Ap3F0M78VzjZmPdNN1V9knZ9i7sZD4em/LukvcvqpdDNX6etAkRWQjYo5m3An8ALk02os6mqpeJyO3AU1j6kkZNxXI7PphoYMNAFEU/r9wWQvgcTabv8Ra0zrE4Aw+qnQ28hDVhu5TEg2/Pob6xNMLArW1uaBakdgWtfJLGLOAGEfm+iHzBJw1kbtEhvHexuILnGqCqU4Bf09zkpbHATf65f+x+bK3shnkLWueYxMALCws2DuoX2ALeLj17YklQ66l8DWWauxvYG1i6mXHMPZ5pNn2Vtx5gd2BnrLvzPBH5gS9qn5kXsPGDzWS6F+BJEdlYVZ9KNqyOdxL29/+7Bt/XBcwPbED9s9g7QghhhYpNo7HrSOVKP3XxFrQOEQ+GPZn+LTXlN5Qu7GJ1kIhskmJozvLihDpfu4iILNPKYPIqXndwO6q3ZlYuUyPAPNh5szMDT8BxLRQvRfRdrGWznrRClRYA/ppoUDkQP5D8bQiHaGYG7nD3BPBk/O8TwO1YIvKdmzmYt6B1EFXdR0SuBNbHZtGsj/1xlJsBrAnclm50ubYG1v1cz1p1c4CVsVYDlzBV/YeIPA2s1MDbZmPnjI9lyoiqXiEiq2J5ulZu8O09NLnUjmMl7EG/0W7+KcCdyYfT3qIoqneMa11yUUErFAp7AhOxwddXFIvFHcr2rYLl21kNa4b8SbFYvCOLOJOgqv8E/gkgIldglYPyltIerKvHpecN6m+t7sbLp9WmNPh6wcskc6r6lIjcB6xA4/euj1oQUh68QeOVM7CJAvuLyPGqmqtVbEIIPfQ1krwK3B1FUTMtv7np4nwFy6VzZvnGQqEwAigCV2DN4L8FrioUCgukHmHCRGQ/4OvMXcZdwL/TjyjXrgMeruN1vdjfo5dPa13XwGtnAa8B57UmFNegY7BegEZvePOLyNdbEE+nexi4ksYnC3wKOAx4PE+TBUIIq2H53y7FZm5eAjwTb29YLipoxWLx8mKxeCVzPzlvgg3iO7ZYLM4oFosXAs8zzNdFFJFFgeOo/uSjwNHpRpRvqjoTW4T7/YFehj0gbOuD0VuukYHLM4C189YK0K5U9Rls2MZF2Oy456l/aaILfEZuY+Jr0bexyWUfNHGIJYBDEg2qvZ0NnB5F0WJRFK0dRdFiwJ9ocqWYXFTQBrAKMKlYLJYPZnyE4T9eYRlqN0v3YF0ELkXxQOetBniJAJNUtdn1B12dVPUJ6h+DOQab+enahKo+o6oTsZvestTfBTcWS9btGqCqvar6G5pb8FuwnGh5sQKWy7LccTR5z817BW0++hY2LXkXmyI8nL1A7afKWeRwCY52oKp3ApMHeMkJ/oSfmi8C/6njdUpOxuoOJyKyMNYy0Yj3sbFRrkEiMgJbp7YZ9Syv1ilux3rmyn0RS5LdsLxX0KZiT1XlxtFcU27bUNXXgYOwm0tlRa0H+JyI7JZ6YA7g4gH2BeAtEflsWsHkkYiMxJ5y63mqFaxF2rWXVWl88PoY4CMR+Y+IfL4FMXWypWj+QWU9EfmFiOShvvEycFUI4bIQwgkhhMuxMXxRCOHo0le9B8vDBzaQx4FVC4VC+eewRrx9WFPVY7HB0NVa0hYDzhSRHarsc631EAMPcF4AuEtEFkkpnjz6A/BT6r/h+ODy9tPM0k/dWDLvFYFbRKSRVCt59/YQ3jsKmzDQTBfpcDMKyx33HtbY8278/UhsRYzSV11y0XRfKBR6sN+1B+gqFArzYLNSbgOmA/sXCoWTgG8By2GzOoc1EZkH2ILalfAu4GAReQR42gemp+ZarGttFWovO9QD/ACfzJE4EenGPttGrn0/itcmfB94ys+VbMXdbUOZyCXx1+7AzxIJqsOp6rsicjNNLlmEnW97Ab9KLqr2E0VRoj1TeWlBOxTLg3MIlkn8I+DMYrHYi600vy1W0z0E+EaxWBzK00LmRGRZrBVwsKR5q2FZj6+JK3Su9WYCdzF42XwuhVjyqJv6F0wvWQ64F8sMfr2I5GlMTTvqYej3rhHA3iLSVIb3nNp1iO9fUEQ+mUAcbS2Yg0IIp4QQDgghLN7ssXLRglYsFo8AjqixbxKddzO8gsbGzXwZOArYvzXhuDITgR/V8bqviciqqjqpxfHkiqrOFJE7sESSI5o4xKZYLq69Eg3M1U1VPxKRe4B1aK4MS7qAc0XkEVV9LJnoOtqrwFM0PytTgMuZe3WbjhFCWAvrJbke+A42/uygEMKWURQ90Ojx8tKClhsiMgZYncYq3yOx7lDXel+kvvNuJrCxiCwrIkv57M5EbU/z40xHAl9JMBbXnG2x1v+hmgH4hIE6qOochtaKJsCaHd4C/Ttg97ir86Moir6DfWZNzYD1ClrnmUlzi9QO65mrw8hU6iufkVjX/HPAi8C/RGTBVgaWF6r6Gtb68k6Th5iWYDiuCar6CrY+6nNDPFQ3nnqjEUOtFM/G7lGdag2s9azcNTSZW9UraB0mTnR6GrZsUN1vw5LpudY7o87XddF/ts9ngfOTDyef4taAY2j8YWYOcGzyEblGxWV4FPWvJPDxW+N/e4G3aGxliVxT1fewVRyamSijwKmq2tS6lMOE0jfGtdTrsQ7WPdwwr6B1pn2xm8/L2ISIaVhrQa2lhv6nqn9PKbZcU9WHsQHn9Sjv1hwJfMW7OhP1Oyypc72eAXZR1b+0KB7XIFU9B/gh9fcAzAQi4HXgVmB9VW22JTWvDqO5BdRnA/slHEu7eRRrRQPoDiGcBdxAkylGcjFJIG9UdZaIHI897S+JVcxGAJ+k+nJDTdXuXdOaTdPgy0AlSFVVRJ7Dlguqx0qq2uii0a7FVPUsEXkXW6B6MD3A/4ALgLPjVjjXmGaHwwjw7zhlzVFxQvVOsx82rhHsb+w1YIMoip5u5mBeQetA8Vilh7GEtCOwE2OgSsEjacTlQET2BzZo4q29wJ88B1fiTgW+VMfrZmKJJ4d1Cp5OJCILASfU+fIuYENgPWAjEdnFz6mGzcJ6ZuZt8H3dWGqn1YAd41nqHdU4EEXRY2X/33Oox/MKWmf6GVY5G1m2baAm6YkisqdfqFpLRMZi3WqNDi2YDpyE5elzyWqkNWA/8pENfbjZD/hEA68X7N73PeBk4P5WBNXB9sYy5g/FAsDP6bB0NSGEc2vti6Ko4RmwPgatMy1N/8rZYEZjC8e71lqc5s65oqoe1OGDa7OyXJ2vG4mdV679LENj17uSXmyNSdeYpRl63aGL+s+94eTVsq/XsYrst+jr9myIt6B1pqewLpl6L1rv4lPN0/ASdqI28vTZC8wnIt8ArvFKWuL+U+fr5mDlIMDmwPLYBINrfRxT5p7EzpNGk9b2AFvHuSOvUVXvvq7PU9j5MJRK2hw6YM3rSlEU/bxyWwhhY5pMAu8taJ3pRGASVkkr1dznYLNoqpkXG5fhWkhVPwR2w8qi3gH/I7AxUpcC/xCRoXYtuJiIdGHdNfV07XcB38RmAF6LdVVfAVwtIv6gm63naK6xQbCVPc4BnhCRCUkG1cFOxsY4D8VHwG8SiKXtRVH0L2wFkoZ5Ba0DqepHWIXrx1ieoD2AA4Ai1fM+jQQu9RQOrRenaPgqja0HOQq7Aa0P7NOKuHJqR6zS1cjf/eJY2c0b/7sZtui2y4CIzA+cS3NpH0q6gYWBsxMJqsOp6nRsotObQzhMacWbjhdCGAX8IITQ8EOEP/l1KFWdgV24PhYvol6ttUCwSQXzUztXmkvODJq7oYzEZkC5ZKxKc6tulOuKj+OysSxDH7AOVknzcqzfHGCRIR5jU+BfCcTS1qIomgFc0sx7vQUtX16jdjfndHwcWlpeo7kKWi+esy5JryVwjNkJHcc1J6lcWprgsTpePBb23SEe5tkkYulkXkHLl9uxtAKVrQaKDWJcWkQOFJF9fd3H1lHVp4DzsHxCUF8rzmzgPeD4FoWVR2djq200sjag0lduvcAU4JSE43J1ipOdHkfzyZ8/PhRwr4gsMPSocmNfmv/cpwO7xPebeRKMqaN4BS0HxJwF3AEsxNzlfg3WevYctkTUCcDrItJMQlVXn+8Dl8X/r+c87Ma6cjyTfUJU9X1gXSzH3KN1vq2UQ2s68GdgTVV9qzURujodiI21HUp3dRewK/CmiHwukag6nKqeD2xNc2PR5sHGbx4DTBGRJZKMrVN4BS0fdgN2GWD/V7EWnfJutx7gxhbGlHerAt9u8D3zAze1IJbcilMrnEbjY/vmAVZR1TeSj8o1SIAfkcz9rBs/xxqxOEMfizYGW6/SVfAKWj5swMCzBpXqY6Lmi2dJueR9luae+D+TdCCOtWluTODaPvO5LSRRSSg31tPZ1K2Q0HFWSug4HcUraPnwLn3jZqoZqKLwUcKxOPMezVUKpicdiOO9Jt/3oS+P1hZaMfO8kXGJeZbUxIqmMu13Oq+g5cNp2IDmajeTOVglbHKVfUXPXN8y1wIvNvG+Y5IOxHEb8HyD71HgV8mH4hqlqh8AZyZ4yMu94l23oY79KzksgWN0HK+g5YCqPgusB9yMzTr7AHvqfA24HhsovTJwC9ZC8yFwAZbE07WAqk7DujmvBqZhT5DTsZZOZe7K9FTg56p6dJpx5oGqzgTWBP6BPazMjP8t3XgUeAd4Ayuj14CfAr9PPVhXy0+A32LXtXoqVzOxGbzlrW+zsevetolH16FU9QVgLewa1ozZwCGq6rPTq/AKWg6IyKLA4cA6wFjspLgRS7uxBpYs8AVsfNO1wKKquktpjUERWUtEzhWRx0XkMRE5U0SOEJHLROQ8EVkvi99ruFPVN4FvAAcBD9HXnVyqoM2kr5IwBthPRB4VkVtF5HgRSXLcTa6p6nvY7OUrsdnM5es6CrAA8AksWfBiwKnAHBFREXlTRO4SkXdE5EMReUJEVkn5V8g1VZ2tqger6jhsMs1RWAWs9MBTaSSwJHY9LOkCdgZmicidIrJ2i8PuCKr6KHZ+nEPjQzC6gV+LyBQRudZn0PbnKwl0uPgm/iR2ApWMBLar8ZZvAeuLyARVnR6n2riN/jes8ozbc4CdRGRzVb01uchz43Rslm21SRzli90LthzNwvH3nwe2E5HVfZHnoRORidgNBgYeG1jtoba8XABWBB4VkTVUdVIyEbp6xIP778AePBsd41l6fRe2VN79IrKxqt6VYIgdJ17T9l6sJa0ZXVj6p68CW4jIZn4vMd6C1vn2pH/lrB5LAN+J//87Bp4B2hV/eRN1g0TkU9g6jo2sy1kyApu59pNEg8qvE7EbdFKzMruAPyZ0LFe/bbEHyCTKsQs4NoHjdLoNab5yVsnvJWUSqaCJyPgkjuNaYrEm3/fJ+N/A4H8npbU8XWOG+pn1JHCM3BOREcC4Fhw6tOCYbmCLkWwyZ0+gOrikr0F+TYs1XEETkYNEZPuy7y8B3hKRSERysTr9MDOJ5pbjuDf+9yEGTtEBNkO03kzsrs9/GfyzHcgcrHzdEKhqLzaLOYnZaOUeTPh4bnCT6D8cY6geTvBYneo/CR5L8XvJx5ppQfsR8BKAiGyGLdewJTYb0JuD28/pwN0NvqdIX2bnvYBX6P9UWrqRTccGsk/B/i5cA+J1BH9A4xUDxcrjJmw9STd0O2Az0XoTOt5bwB4JHcvV72YsrVASaTJKs3XdAFT1CWyCTRLeBH6Y0LGGvWYmCSxOXEHD1uG6RFVvEpH/AfclFZhLhqrOEpEvYGPKtgGWw8rvAuxmtB1WQZiD/T1cA1xcygOkqq+IyKrA17CxHXOwJ5zJ2KzQacA1qvpOmr9Xp1DV80XkYWyA7HrABGzM4BzsYjUeWBCb4XldvO0NrPXtutJMWzc0qvqgiKwIfAXr3v8mNvB/NDZZYwSWCuUFYGlgXqyM3sJmRD+APawuiqWrOUpVP0z3t3DxdWtPEbkM2B7YCCuTMdj1Teh7IOqmb6zaNOyhpws7v04ELojXa3WDUNX9ROQm4GBgfeyzLq8klz7n8lVr5mAP+G8C9wD/xOoTfi+JNVNBewdYCrvJbwEcGm8Xmhvs7FpMVWcDF8Zfla6s4/3vAxdV2eVdOAlQ1ceAx7KOI+9UNaJvJucvmjjEyQmG44ZAVf+J3fBdSlT1BnxNzUQ1U0G7HPiLiPwXmxp7fbx9DeDZpAJzzjnnnMurZipoP8OWRVkaOFBVp8bbF8eSN7phIM4X9FVs/OAKWJeNYOPKxmHdO4o1P0/D1iu8BzhXVV8Wka2x7qBXgFPjRJ9uiOLFt78MfA9bQLjU5bIYNitwKpbX7j6s1e0TWPfNi1gXzWjgCVX1wc0Nij/7bYAvYte4l4AvYOfD4lgqgTHYefJB/LbSeTMNO1+6sfUJzwLOB76PdVu/ibU4/wsr16WAp1TVW6FbQESWx4YMbI3NxByPfebz0T+dyhz6d78pNvRjDnbePYK1Ct0ar8jiBiAiC2PnzLZYl/988a5S93Jl+pPZ2NCBHmwYR2k1lcnY5/4C8Cp2/ZsG3JKnoQMNVdDi6ehnAL9W1efK96nqiUkG5lonTotyF3ajGEz5NPMtgcNF5D5gg7Lth4nIan4BG5o44ePfGXiJrYWxMVJb1tg/ExghIr9U1V8mG2HniitnN2GV43pUphYqTyo8Dssf+LuK15Ru/iOwchopIsep6oGNR+xqEZHvAedS35Cbaq8p3ReXjb++CcwWkZ1V9S/JRNl5RGQNbPzlAtSfh66Hvs97VNn2JbGK3hzsIbU02/1/IrKRqr429IjbX0OzOOPp6L4+4/B3HLasUzO66V85A2tFuL7Ka11jdsaWfhqKkdjF8XAR2WjoIeXGbtRfOWuW0Fc+o+J/9xeRzVv8c3NDRJbCxhEmPR66GzhPRDwvWm1/xx5ckkr2DH11lFJFbmksM0EuNJNm4wqGfhNx2Vqf5FeRWDbh4+XROiR3cZtOctm98+CLGf3cmYCv+ZicVUi2glCuG1i5Rcce1kRkNNaV3+rViUYCn23xz2gbzYxBexbr0toA+DcVq9ir6klJBOZa6g3q695sxEeDv8QNYkqCx+pJ+Hid7pWMfm4XXk5JmkLrsgl4WdU2Pf6aJ4Wf9WYKP6MtNFPb/QHwLvbU90Ns0kDpa9/kQnMtdBhDy5peLZnnwUM4njOnA0nkXZqJPUhdkcCx8uK32A0mDaVzbybwMnBxSj83D/6N5QtMKuFwyWzgamzSgKsQ52M8jGSX2Zrrx2Dnzs9b+DPaSsMtaKrqXVnDnKreKSKfx8airUX/wZkwdxdB6cR4B8sYfR423mAlbFbh4arqGe2HSFVfE5FVsAHOG2GDyUuJNStnnr2N3dzHYE+tpRag2djyXIflabbTUKnqW/HMvyuB5bG/6w+x2Zsj6ZuFVq+3sbyD38ISDc/EEjxfDayOJYyeBByiqh/UOohrjKrOEZFtsJv4dsCn6LvPNdP1ORub0XshcHQpgbebm6oeJyJvAD8G1qW5HrqPD4fdb57CzsNRWCaBP6jqP4Ya63DR9AcoIiOxcUeTVXUo6wm6DKjqPcCGQzhE5UQBlwBVfYnWD1Z3Vajqy9g4wCTtk/Dx3CBUdQZwePzlUqSqF2Cr1LgENLNY+mgRORur1f4Hm1WBiJwkIv8v4ficc84553KnmTFoR2NN9JvQf8zGP7C1z5xzzjnn3BA008X5DWB7Vb1XRMr745/Aptk655xzzrkhaKYFbREsTUOlMfRfvd4555xzzjWhmQraA8BWZd+XKmW7Y2s1Ouecc865IWimi/Ng4AYRWSl+/z4isjKWnf4LSQbnnHPOOZdHzeRBu1tENgT+D1txfnMs79L6qjop4fhSUSgU9gQmAqsCVxSLxR2yjcg555xzedZUHrS4IrZLwrFk6RXgSCz/1MIZx+Kcc865nGu4giYiawG9pdYyEfk6sCs2i/MIVZ2ZbIitVywWLwcoFAprkIMKmoiMwWbjboxlNF8emBfLdt4FzIdlq38Py769KLa+3Rz6xi3OKvu/Yst/TcEyb4/F1uacg2XhfhlbImUO8G1gNLYSwZ+AQtnPfw64ERgPfB5b5/VqVX2nyu8wAfhSHMf1qvrqUD+XrInIPNgQghWxPINdwHpAwDJpl1YUqFxrUOOv2VgG+3uBz8Tv+xB7AJmKLf10qqrOEJE1sQeSUqb1scBrwBWqemuT8Ut8zM8AL2Fl18qlX1oi7iFYE/sbXhLYGltRYAFsdYcuao/fVaqvxFG5rRcrx9L2GcCr2Hn0evyzlwOWiLc9ga1EcFsc01jgblUdcOmheGWKjRjgXMoDERkP/AW75o2mepkMpnSelZf9LOz61g2Mw8rqQeCa+GspYG/sHDsbK9vdscl2twAXqWrSy1JlJr4G7AjsjKXjWoDGV+Go1zvA7cD82LkzFXg4/n51rIfvYuxc+SlWRmer6u0ishx2rZqN3T+yWot3QM20oP0JW7duUvxL/g24HFtWYzS+HmdbE5GFsckcExj8pBk/wL7Kv51F469KK8f/Vl4QN8By6nWXHWs2dsEbQV9lcYqIbKCqL5T9DpsDxbJjTheRTQa7WbUzERmLXVAGekCotQi0lO1fDKt8l8xPX7lsBOwpIn/AluyqVsn4qYicoqp71Rs7fHxhPgsbKjADK8PbRWTL4XQDEpEjgF9gv0MzCz9XO6eqbRtZ8f28WIWsmk8An8bK9XDsPOkFRorIj1X1zKqBiHwXy+peOpfejs+l5wf9LTqIiCyDrU1bAw70uQAAIABJREFUfs1qpsJQ7T0jsBV1SsYBywDbAB9glYOSneh/HdwJ2FtENlTVj5qIp62ISBf2gJ3WSigLAt8c5DWH0/869z0ROQf4btm2GSKyqao+lHSAQ9XMLM5P07dg7HbA7aq6I3Zh/lZCcbnWORb4JK15ohlItZ83iv4Xze54Wxd2cxyJPWme/vFBRHqAS+LXzRP/Ox/2dDycnUM6rbfLASdS+9wXYA8R2azB434NG/bQhVU2erAK4U+bjDN1IrI2VjkTmqucpaUbi68LOE1Elqx8gYgshK3pWn4uLQSckWKc7eJ6hrYuZDOE/pWz8u3l1sDGc3eC72K9Gu2k2nVuN+x8GBV/jQH+mmZQ9WqmgiZl7/sycF38/5fIQfdgB1ib9C9WQzECu4iVLIY9pZbrBlaIW3GGqzUGf0liBvuc5mATZhqxKtaqU25EE8fJ0ipYy9lw0gWsUGX7BOzzLzcCWK3lEbWfWi2T7UBI99xvpVVJ/8G/WeVxdgPLxy2AbaWZgB4EDhWR72FpNa6Nty+L9a+79vYSdgMeLuZgY6hK3sK6eCq9parDOVFyO42BEGxMRyNeZe7rSS82rm24eJW5KzXtrlZZVfvctcZrO917WQcwiJezDiAhw/lv621Vbbv7YjMVtH2BtYCTgaNU9dl4+7bA3UkFlqZCodBTKBTmwVqWugqFwjyFQmG4XajrdQhzt3RkqdpJMafsXwX2L+2Ix2ocUvaa0sDdfVoYYxr2IJ2K80zgKmysXzWKDUa/vMHj/hV4Jj4+2N/YO8AfmogxK7dgg47bfaJT6UFkFnA+8ORcL1B9EbtGl8p5Tvy1f+Vrc+B7Gf3ceh4YpwLHtTqQlJxN9VWG2s0zzH3/aMux883kQXuM6t0WB1C9ZWM4OBQbTFiyHXbhm5hJNC2kqo/EY20OxMYIfYK+8SylAaxC3x8u1FeRL80inIP9XZUqV9OwJ9iHsDEwa2BNyo9iN4uJ8bZ5sBvNX7DZh1/GBtmepqr3V/wOx4jIC9ig6VnABap6U2OfRHtR1cdFZHVs3NCS2Gc3Apuo0ch5qthFfzT2OSv2Gc3AZjjtAryIjQ37JrA09tmPwCpUlwFHqmpDXX2q+pGIbAD8P6yr8Hngt6o6HC7YAKjqbBHZEhsTtC722X0Km2U8kuYeaBsxByur6di5NB99s6enYgPd78D+JsZhMzpPGaDleG9gEpar8gPgdFW9r4XxtyVVvUlEvoQ9mIwh+W640s2+K/7/u8A/gD9is0Z/gJ3DlwFPA3thk3fuB36mqi8lHE8mVPV9EfkMNmZ4S+yzbuU5o8D78f9nYPeZZ7HPdimssvhHrHdv1ziWy7Dze1tsIscs4M+qemML42xaYmORVHV6UsdKW7FYPAI4IuMwUqOq/6F98tjdW2P76TW2A6CqF2NTqDuGqj6OVQzScFL8lRhV/QBr3Ry24jRBv8k6jiTEFbczyOfEgH7i1DHzZ/Cj72Tuv6cBr23DWZzGZfus46jiVxXfXxJ/tbW6a7ciMkdEZlf5ekdE7hWRbVoZqHPOOedcXjTSglYr38h44LPAhSKyi6peOvSwnHPOOefyq+4KmqpeNcDu80XkCaxv1ytozjnnnHNDkGQ+rJuw9SzdwFacPHly1jHkyYr1vMbLJDVeHu3Hy6S9eHm0n3rKJHFJVtDmxWYfueqmYOsiXph1IDn0Ifb5V/IyyYaXR/vxMmkvXh7tp1aZtIwkldtTRP4ITFDVryZywA40efLkpfHVFrIwZcKECS9W2+Flkgkvj/bjZdJevDzaT80yaZW6K2gicnyNXeOAdbClRTZS1YcTis0555xzLpcaqaD9s8au94GngFNV9YWkAnPOOeecy6vEujidc84551wy/j979x0mSVn1ffz7m52ZZYElI2EQUEBAyQoo+igGBAwtiiIiAoKoIIqKAZRkTqAIiI+ioAQVXhUoQYL6ACIqKEGCBEkLO2SWhQWWDbPn/ePUuD093TPdPV3Voc7nuuba2a7uqjNdXVV333Xf53Rc9fYQQgghhKJr5SzOMIkY3Nk2MeC2s8T+6DyxTzpL7I/O07mTBMLUpAfVbXgR65Cv54BNKw+u2CdtE/uj88Q+6SyxPzpP1X2SpehBy89q+EG1N36AhXxsiucMWg2oPLBin+Qv9kfniX3SWWJ/dJ6J9klm6mqgSboeeKOZPSnpaOA4M3su29B61m0bbLDB9e0OoijqzLYd+yQnsT86T+yTzhL7o/O0q2pDvZMENgWWS38/Blg+m3BCCCGEEEK9tzhvBE6X9BdAwGckPVPtiWb25VYFF0IIIYRQRPU20PYDvgS8DTBgV2BxlecZEA20BklaFTgdeD0wH/gOfhs5ZnC0iaT18P2xCLjUzHKtwRZCCKHY6mqgmdkdwJ4Akpbg49EezTKwopC0DfBbYAjfH8sDX8cbBie0MbRCkrQ68GngMPwLB8A8Sa8zs1vbF1kxSeoDdgLWB/4DXB5fXEIIRdDwLE4zi+S2LSLpY8BJ+G3jcv3A4ZJOMbOF+UdWTJI2A64EVqlYtBLwC2DL3IMqMEnTgN/gPfeLgEHg55IOiEZaCKHXNdXYkrSBpJMk/VHSHySdKGmDVgfXyyRtCJzI+MbZqDWAv0parsby0Hq/BFas8vg04KU5xxLgAOAt+Pu/DH6++gCwezuDCiGEPDTcQJO0M/BvYDvgJuAWYHvgVkk7tTa8nvYyYGSS52wJHJdDLIUnSfhs5Wk1njInx3CC24rx56gRYIs2xBJCCLlqJlHtN4Hvmdnh5Q9K+ibwLeAPrQisAB4BBiZ5Tj/wUUnXmtnpOcRUWGZmkuYAq1dbjI9JC/l6GG+QlTea+4GZ7QknhBDy08wtzk2Bn1Z5/DTiNlAjrgEuxMfWACyZ4LmnSnph9iEV03XXXbecpF8wvnyKAc+nP4dL+p/cgyu2H+I9l+UzxvuAQyUd3J6QQgghH8000B7Dbz1U2gqImZ11Sgc5vwtP/HsxcDaeb67a4GcBm+UXXXGYGR/60Ie+i49rGh3vZ/iMwYV4L+cM/IvJn9KJBCEHZvYYsDXwBEuPC6U/J0qKgtEhhJ7VzC3OU4EfS3ox8Ff8xPka4PPA8S2MreeZ2SLgG+kPkpbB66utX/HUPuCbkq4ws/m5BtnjHnjgAebNm/eKioeFT9KApbfX+vDbbR/AP+shB2b2sKQBxk+mmQasB0R+uhBCT2qmB+0reDLaj+MpCf4MHAIcC3ytZZEVkJk9jxfBrTZ5YBPgi/lG1PsWLFhQa1G18YHCe9NCvmZRfQjAoWmetBBC6DkNn9zMfc/M1sFTEqxoZuuY2fcjN9HUmdnV+ESMygvSID5bNrTQeuutx8DAwKOMbRQvxL94DFY8fRpwQ16xhf86mOqVS94HHJhzLCGEkIspffs0s3lmNq9VwYT/uobxDbQR4ME2xNLTBgcHOfDAAw/BZwyO+gvwHuDDjG24jQAnSYqGco7M7O94+bPKRlo/Xo4rhBB6Ttwe6EyXANfhPTngDYOFeM9aaLFPf/rTdwMvxmchvwh4U/rF4zTgmbKn9uO3OH+Ze5Dhbqp/aXm6DbGEEELmooHWJpLWkLRJOjFgjHTywBuA7wJ34rPYbsJntIUMmNlCM7vNzO4ru1W/OuMrC/QBL5LUzASb0Lzz8AkB5WlpluCpOEIIoedEAy1nTz755DRJP8Fvqd0GPCTptZXPM7Pn8B6CFwMvwMefnSnpkDzjLbg5LG0QjHnczKqNiQotIreDpL0kbWtmc4FXApcCs/FqJh/FK5mEEELPiQZazvbaa6/9gH3KHloRuEjSC8qfl/asHcHYVCh9wFezjjG4tCfzMMbfWpsmads2hFQI8+fPF/BzfCzgacA1kk7EG2ZHA9PxvIA/BW6RNNSuWEMIISvNFks/WdIqrQ6mCIaHh3dkbAoH4YWgt6t46opU3z8rRGqB/JjZScCvGdtImwlcImml9kTV2w444IC34jM0hTfGhM/kfCdwEbBq2dNfAtwpaYO84wwhhCzVfaGXtE7Zf/cClk8fvznKENWvr6+vWuKtPqDy8cfw26CVvTcCDsogtFDbZow9VvqAlYBt2hNObxseHn4p4xPTLsbHZa7F+PPWssCFacH7EELoCY30xNwuaVZas3AZYLRRtj6TF/0Oqa222upcxja6FgH3AleXP8/MlgC7Ub300/clrZxZkKFStUa1WDrLNrTQMsss8xTVkzU/XOWxUZtQvdB9CCF0pUYaaCviuaGuS1/3e0l34rcgdpa0Zgbx9Zyf/exnlwEfAobxFA5XATumkwLGMLNrgHGP4wlTo9cyPz9gbINhCT6j8F/tCae3HXzwwb8B5rF0gsZC/P3+AXDOBC+NMmghhJ7RSANtwMyuNbPj8RPh1sAH8QvX/sDdku7IIMaeY2anp9UXZprZG81s9gRPv5fxtzmX4AOmQz5Ow4vaj/Zmjt7ivDDSbbTeO97xjsfx28fnAH8HzgJeYWZPAvsC/6h4yULgV5E0O4TQSxq5uDwt6Qb8VtwgsKyZXS1pMfBevMFQOdA9TN2HgSvwfTUNv7V2HdVvAYUMmJlJWh0fBzV6O38A2AF4B/CbdsXWq8zsfrwwfeXjCyS9Cp/NeTA+/uxu4Ff5RhhCCNlqpAdtbTzFwwK8sfBPSVfhjbVt8OvYX1ofYrGltzkPx9/z0UHQWwKXRu9NrtZl/FjLRUCkeMiZmY3gQwNWwvfJS4HzJB3a1sBCCKGF6m6gmdnjZvY7MzsCHxe1LXASftvnOLyH7cpswiy8fSr+P1o4PSoL5OcWxk8KGMSTDYf8/QjvUR5gac/y8ZJmtjWqEEJokank03rKzM7FexHegNcwPKUlUYVKK1R5zEhTnYRcfAvPXr8IH4M5gidK/WM7gyqwtRmfimMaEJOVQgg9odlbZFvgsxABZgGLzOxhJp5hFZp3GX6LbbDsMQE7AZe3JaKCMbNnJb0S2APPxXULcHFZ3c6Qr3vw1BqVXzLXBf6TfzghhNBaTfWgmdkDaZ4uzGwzM3ugtWGFCp8Fnq7y+BGRQT0/ZrYAuBV4C57y4XeS1m9nTAW2H+N70Az4ZVTaCCH0gjiRdYE0R1q1BtoSvJh6yIGkTfBZzK/BEzS/Ga8TuVo74yqo66o8JjxZ7QuqLAshhK4SDbTucRee5qFcH36LOeTjw/g4p2np/wfwBM57tC2igkp78J+qsTgaaCGErhcNtO7xKXxw+kJ8gPoS4BQzu7OtURXLCowft2lUn8QRsveFKo8tAX6WcxwhhNBy0UDrEmb2b2BHvNdgtAdnX0lvbFtQxXMV46s6TAci/197/L7KY33A5nkHEkIIrRYNtO5yFH5LDXzfLYcn6IzcT/k4A/hh+rulPycDUeKsPZ5gafmtcnPyDiSEEFotGmjd5dWMTbUBMBPYqA2xFI65jwObAX/FB6V/HHhA0jvbGlwBmdkzwNcZW/bMgM+1J6IQQmidKBXUXebis9SqPR7ysydeyWHUdOAcSRuZWUzayNdR+ESZd+FJhH9iZkl7QwohhKmLHrTucgxjb+ksxAt139uecIpH0gzgk4z/cjMAfFtSZW6ukKG0V/NUM9vVzErROAsh9IroQesiZvZLSQvxGZ3L4YOkvxTZ7HN1GLBsjWV74Lc+v59fOAFA0gDwfjwv4N3A2WZWmZYmhBC6RjTQuoyZ/QbvNQvtsTkT9zwfRTTQciVpdeBmYA28h3kE2E/Sm81sUVuDCyGEJsUtzi4kaaakMyQ9ImmWpEPj1lpuhvGxTrXU6l0LGUh7zv6JN87AJ27049Ue9mtTWCGEMGU904NWKpUOwU/ImwPnJUmyZ9myzYCf4EXe7wEOSpLkqorXHoEnHP09cGCSJNVKK7Vd2hBLgB1YOqPz+PTf6LnJ3neAvYFVWJqPrtxjkhS3nXOzA14gvdI0YnZzCKGL9VIP2oPAV4FTyx8slUoDeIPmPGBl4JvABaVSaeV0+U744Pu3AkP4jLyT8gu7YRviCWvL021MAz7dlmgKxsweArYETsQb+5WGgP1zDarYVmR88mDwnrR1co4lhBBapmcaaEmS/DZJkvOBxysW7YjfdvpOkiQLkiQ5C5/1+K50+X7A6UmS3Jj2mh0FvLdUKnXqraoZNR5fLW5z5sPMHjKzTwPfAhZULJ4GvCX/qArresbXqB21h6S18gwmhBBapWcaaBPYDLg5SZLyb9k3po+PLv9X2bJb6OzbI3fiGdQrzQD2zTmWonse76kpZ8CzbYilkMxsNvBeqlcUmMbS4zyEELpKERpoy+P1K8vNxTPwj1ueJIkBT5ct7zRrUX3sk/DbtCE/F+GflXICVo3ezPyY2fl4T1o18aUlhNCVitBAewYf/F9uRWDeBMtXKFveac6leuNxCTA/51gKzcyeAP7E+N6bXfBZhCE/x9d4fNdcowghhBYpQgPtFmDzUqlU/rdulT4+unzL0QXpjM8R4D+5RVinNKXAy6ndg3ZTvhEFYG3G3+ZciCdMDfmp1YO2iqS9c40khBBaoJfSbPTjf08/0FcqlZbBG1pX4GOFDiuVSicCu+MXz/PSl/4MOLtUKp2NTx74CnBOkiTP5foH1GcxPih9mSrLDC81NM/MfpRvWIV2B/BKvNTTqOl4vrSQn//gFQQ2qLLsNEm/NrPnc44phBCa1ks9aEfit/i+CLwn/f3UJEkWASXg3fjYsy8CuyVJMgcgSZI/4I2yS4CH8CSkH889+jqkubW+TfVZa314T86JkgarLA/ZOBr/IlDO8LJDISdmtgS4rMbiAeAVOYYTapDUJ2lHSXtK2nLyV4RQXD3Tg5YkybHAsTWW3QxsP8FrT6Kzc5+V+xI+MH1/4KVVlg8CqxM9OHl5ivE9mn3APpI+ZGaVjbeQnccmWHYUsHNegYTx0i+OCbAT/mVSkhYDlwIHm9n97YwvhE7TSz1ohWBmS8zseDN7GX5LtnKA+gjjZxaG7NT6ktMHnJ5nIIEzqF2G602SVs4zmDDOZ4A3sbS3H/z4eStwq6T12hVYCJ0oGmhdSlIfMIvqebiOyT+iYjKzp4G/UT0P1wck7Z5zSIVlZncDH6mxuA/YOMdwQpk07cxHqT7BCTyZ+OH5RRRC54sGWvc6hOqpHPqB1+ccS9G9i/Hj0EadKanapI6QATM7nfHVREZVS/Ac8vEK4IUTLO8DogcthDLRQOteb6L27bXNJG2aZzBFZmYPAzfXWDwD+J8cwwlwAGN7NBcB5wB3tSecgBe0r3X7eVTHpTYKoZ2igda95lK9SDT4bYRf5RhLgM/WeNzwahUhJ2aW4MmCLwH+CnwZ2DudBR3a404mn5S2vaQPRhWOEFzPzOIsoB8Ae+GNtMqG9jRgc0l9afqBkL0H8MkZlVUpFgL/yD+cYjOzyyQtwRM7P4bnpqtVVD20iKTXA9sCc4Bz0zGa4LcvnwBWm+Dl2+O3Ql8AfCvLOEPoBtGD1qXM7Br8NuecGk95Lhpn+ZD0ArynZtmKRQuA3dOC3iFHkr6Fp284FvghcEPM4syWpGPx0mdfwr9A/kvS6mmj7QImbpyNmgZ8WVKtyQQhFEY00LqYmV0BXFdj8Xk1Hg+t9078NmZ5j/QIcJaZXdSekIpL0rb4Lec+PEfdIN6D8+V2xtXLJG2OJ20WS9/ztYFvAgczfrb5RAYZ/2UnhMKJBlr3+xt+G62cASe0IZaimsH4NBt9jC3/FPKzCeOPiUFg3xjflJmNGT8JYBCvc7wc9TfQlgCzzGxeC2MLoStFA637fQu4Gj+xLcAbCp8zs1o9a6H1rsQvRpW2l3S+pB3yDqjghqm+P2bitz1D6w1T/QvJxviXyHrH/83H09aEUHjRQOtyaQHonfAyNvsDW5rZce2NqljM7Abgg4ztQRB+cSoBf5b02nbEVlBXULtBsJOkyokcYQrSEk7zqN74XR44gtpjZSstB+zaotBC6GrRQOsBZjZiZn8EfgM8HYlR82dmZwBrAtdWLBJ+nH0p96AKKp0cM1GD4MC8Yul1kl4C3IHnAdylxtNm4DMz6/VVSa+aamwhdLtooPUISXvhudHuA56SFBeh/K0FbFflcQHr5BxL0X1ngmWfzy2KHpaWm7uQ1n+2R4DXpusPDZK0kqTPSzpB0oExI7Z7xQHQAyRtD5yJz54CH3/zI0lvaF9UhTRRL8Gs3KIIAN8FTquxbKU8A+lhawAb0Vw+zVql0cBTbXwOWCzpcUl7NxNcEUlaBbgB77H/GJ7uJInGbneKndYbdmX8DKrFwFvbEEuR3Ub1oungubhCTtKqAR/Fe5XL98liapflCo2pdRuynvyLtwLfnmD5KnjP86rAGZJ2bjC2ovocnt5kOt5wHgDejI+FJc1Ld7akOyX9WVK1es6hQ0QDrTfUqnEXmdNzlNbkPL7KonuB3+ccTuGZ2SJ8RuB8PO3GIuBJ4P3tjKvbSVpD0pn4mNdaFkywbAT4m5l9Hq+bOukmgW/F5I66bMj4GcyLgPUlLQv8BXg33vP5auCK9A7Mf0naTtJ3JZ0oaac8gg7VRQOtN/w6/be8p2AAeGVkT8+XmX0W+CQ+q83w3oTZxG21tjCzy/HZtB8E9gY2MbPb2xtV95K0NnATEzdyJ8sBKOCs9Pd6EzlvDlwXjbRJ3Un1HIB345M4XszSBlwfvi+OlTRT0pslHY2nRfk4cBBwaYxnbp9ooPUAM7sDOJ/xySD/Bx9/EMk587UQz4Q+OoPzlcCf0nQEIUfpBf3DwNuALajv9luo7ShgZSZPPFvZSCi3t5n9BcDMHgGeqmO7fcC6wCfqCbLAvoWPdy3/nE/DG7grMP6uSh/ecHsSuBgfu9aH3x7tx/fzyXHuao9ooPWOt1R5TMBr8DI3IT8H4yfFUQPApsA27QmnmCQtj6c9+TzwvvTfO9JeoNCcDamvQkYf8HSVxw2vIVyu3jGB/XgjLdRgZk8Bp1RZ9FX81nKthtY0arcHBoHVpx5daFQ00HrHRCfN1+cWRYDaJ8HpuUYRDgBexNL90Y/PtL1f0svbFlV3u4XavWOjQywW471iX63yHOE9++W+18D214phG5PahvE9nM/jCYObueY/Bzw61aBC46KB1jsuovYMwuPTAaIhH1fWePx8SZ+IW865qdVTNg24PPJDNeUrwAN4I+15vFfm98A7gX/hF/KrgB2A/6vyemN8EuE/1bntPnxG4r8kRY9ObXMYfytzBt6L34jRMbT7phNuQs6igdY7DqX2rM0V8UGfIR//rPH4Snh+rmPzC6XQ/k3tc9xMInlww8xsDrA1PhHmW/iMwLeZ2flmtrWZrWFmbzCzu/B8XH9kaY+bpT+VVTWepv5SUIN4xY6jpvaX9LST8Zmbo7nmJso5NxEB+5jZryd9ZshENNB6R+W4jnLCZ++EfDw5wbJpwBGRODIXZ1G7sQyeJy00yMzmmdkPzezYtGFWtec+LblVAk4CbsR7lnc1s4srnmfAdQ2EMABs0lz0vS9tHG+L92zeDDw8hdW9rSVBhaY0kwE6dKaV8B60amPRFuO5uEI+Lsbf81rH1wBe9eG53CIqIDMbkbQbMMzYSRujPiPpJDN7MOfQCsPM5gOfqeOpjeZsjMocFSS9FU8lsyme2ucKPHn2W/Fex2Zu6a/WqvhC46KB1jv+ztJST5UGgMdyjKXQzOwZSXuyND9dpQfMLBpn+TgCH0dT7eI0DTgB2CPXiLqUpPfjKUv68c/299Nesla4EtiJ+q9Jl7Vouz1B0sfwnsrR8a1b4g2zJeljqvi9Hkbt8bQhB3GbpUeY2TXAp6k9UeCEGJyeHzP7DZ40spoBSUfGIPVcbEXtGc7CL2RhErvvvvt7gDOA1+ITAL4NHNfCTRxPY42u4yS9tIXb71qSpuMzYaud3/sqHq91Tqq6auAYSY9J+oaketKrhBaKBloPMbMTqJ1ra3k8eWrIgaRXAi+psXhN4Gh8RlzI1ixq3z4bAe7PMZaudcsttxzM2OtFP/BJSTNbsX4zW4yPd9obeKaOl6wL/FPSUCu23+VWo77cdAvxc08j+tP1Hw4MS9pEUn+Moc1HvMm951/A44ztSTPgUTN7tj0hFdI+kywfAD4VvZqZ+xLwLNUrCPThxaXDBMyMkZGR5aosEi0sYWbubHyWaD1mAKe3avtd7BHqa9QOMrVr/ur49WUB8LykU9Peu5CRaKD1mHRG1PvwadajFyUBV0uKMYf5qSdv0ABxDGbKzO7Bb2NWK94t4JB8I+o+klh22WXvYGxP5BI851nLJ1iksxCPrfPpO0nartUxdJO09/H9+P6pHOIywtLrQB9eI3gq4wZHG3kDwL542qCQkbg49KYrgbmM3b9vZ3z+oZCdy5j4RGjAX82s2RxFoU5mNtFtzv0lHZxnPN3o0EMPPQLvqRnB38tngHdk+Pn9MvXn7/q1pEZv3fUUM0vwnsejgXPxHHVvBm5nbKNtbVp33R8APtCidYUqooHWm/bGS9qU68dL34QcmNlFwJHUnrQh4Lf5RVR4D0yw7CRJP5JUaxZ04R1wwAGz8fQNJTw57YZm9veMN1tvT88Q8PuijouStIuk+4Dr8Z60483scHxm/8uoPoO5VbNvC/me5yXe3N60S43HZ+QaRcGZ2Tfwi0e122sAe+UYTtH9e4JlfcB+wP/mE0p3ShPU/t7MLjCzTNP2pEM1Emp/wSnXh/ceFS4Zt6RtgQvxSRMDwEZ4GbMXAQdN8NI+llZ2gOqF7etxU5Ov6zlDQ0MXDw0NvXdoaKhWLeaGRQOtN61Y43HFVOncPUftIumRZiMHkkZ7fiYyCHwgxml2lFOpP2cX+C29onk/3sgafZ+mpT/vor47JqvhDdvVmLgCSi3bS9q/idf1ov8APwIeGhoa+sHQ0NArprrCaKD1GEl74Qkfq5kJbJ9jOIWWznC6fIKnPByzOHPxIeq70ItoNHeSxTR2K+5kSR/JKpgONcj4z7bhvWkTfRkfAS42szmj/lRLAAAgAElEQVRmdi+efmP5Jrbfh6fgKLzh4eFPAFvg++RlwN+HhoZuGhoa+vTQ0NDqzawzGmi95zhq79cRmjsIQ3PeDmw+wfK3AH+NRlrmZjJ5w2sJcIWZ1bodHfJ3HT7ZqV4CfiCpWkqQXnUR48/3g/gkpd+ytFB9pT/hvW9IWpWJrxuTiWvKUs8DzwwPD+8IrA+cDXwQmD00NHRBoyuLBlrvWWWCZQuBG/IKJPACJk+38UpgtxxiKbIrmHwsUx8+czB0CDObi/d+NmIa8NMMwulI6WSkT7J0xuvzwF5mdj3wReA3ZU9/EHgNsLyZ7WxmcyXtjCdr3oPmeo9HgD82G38vGx4eng2cA/wKT2+yc6PriAZa77mZ6ikFngXeaWaP5BxPkd1I7fqo5WpVfwit8Uvqy9c00Zeb0B7NVNvYQ1JhJgyY2Yn4uOMNgJXN7Jz08QVmthc+vmx9YF0zu3o0Ybmk1fBetqlUmJkDfHwKr+9F/UNDQx8ZGhr6M3AP8E48/UnDqWCigdZ79sYPmkX47MEFeLbtn+MDOqM0Sk7M7K94PqLJ3JJ1LEWWZqj/HPAiJr5lVuhcWp1G0vL4WJ5GLQGmPEC7m5jZs2Z2j5k9X2XZE2Y2q0rOui2Y+sz+C8zsqSmuoycMDQ1tB5yAf9H7MvAPYIvh4eFXDA8PnzI8PNzI7XrAc2OFHmJmd6Sz1nbCB4m+Dr8HvgS/zfMpSdum2bpD9v6B3waY6PbBwznFUmhmdp+kZ6ldnugHkjCzU/KMK9T0PJMfO7VEWbvJPU1js2Sr2U/S4Wb2RCsC6nJXAZfgM2gvHB4ennIS52ig9SAzmwOcI+llwJnpw6MnuT7geOAd7YitgKYz+UXmMkkrm9lzOcVUZNfitxxqOUnSmWY2L6+AQnVmtljS/wIfa/Cl0/Dz395mdn4GofWKG4CrgW3xiQXN6Acel3S0mTVzO7qX/AQvf7YVsNXQUPWbVcPDw3VX9IlbnL3tRYwfj9ZPjHnK019Y2ntZyyBwaT7hFN5kyWj7gLXyCCTU5SdNvm45vATURLOou5qkPkn7SPqepCPS2Zh1S2957ooPf7kbuJXmKwx8WdK7mnxtr1gZ2HiSn00aWWH0oPW2x6m+j9eR9BUzOyrvgIrGzB6Q9HZ8NtUKEzz11ZJWMLNmM3qHOpjZZZIWUDt58AgZFAAPTZtF87c5pwGfBfZpaUQdYNGiRQBnAe9haUWAj0vaxszqHjKR9hR/ePT/knbDZx4206P2MQpcvm54eLjllWGigdbbdsO/EVXrKT1S0klm9mjOMRWOmf1R0hp43qFa2dEFPCTpVOBzZlYrf1GYuj2A86m+Hx7Ea07+I9eIQlVm9qSkw4HvNLmKPSTtm5aO6hlHHHHEK4D3MvbcvhpwDBOXeJqQmZ2fFp6/GS9T14hXNbvdXjA0NLRePc8bHh6eVe864xZnb9uQiffx6XkFUnTp7KpBatflBJ/ufhBwci5BFZSZJcB2wMX4mJFyQ8CVkjbIPbBQy+1TeO10ape+61qzZ89em/FJaAeAl0x13Wb2JF6btlHLjCbdTm+/bihpkwKVT7sHuDf9t9bPvY2sMBpovW2yGTpvSaeyh3y8hMnzog0C+xfopNYuNwPz8WTC5frwW2N75h5RqOW1U3jtAqDn0kCsu+66sxl/G3IRcEeLNnED9RWqLyfgEUkvAf6O16a8DbitIHnpZuDn9xkT/NSTF/O/4iLQoyRtzOQFosG/YT6TcTiFJ2kH4NA6nz4N/zZcLeFwaI2v4KW4qhmgwRNpyNRUUkFMx9Me/GayJ3aTr33ta9f/9re/PRMfNjE6Cekx4NhWrN/MnpB0OfCGBl+6Ot64K68Duj5woaSX9dqt5nLDw8MtH5YSPWi967VMfIFfjPcizMknnMJ7IxPf3hy1CPibmc3POJ6iez+1i0kL+EOOsYSJXUbjvTnl3tyqQDrFwMAAeH7LffDkqEcCW7R4TPE3ae59X5axx1Y/Pq6zsrc6TCIaaL3reSbevw8CpV7+RtNhnqe+noAHgHdnHEuhpdU0JkqlscTM/pxXPGFiZvYH4NM030jrybQpaYWMX5rZ583suHTsWCvX/we8zmezqTcqjatyECYWDbTedREw0QH7QuBESbV6EUJrncP4Qb3V/MLMIs1DtrZg4ov9NXkFEupjZicAqwKvp/EvMD03Bi0vaZ3P/8fUejAXAudGSajGRQOtR6XVBF5D7UaBgF2AwyStL+nFkprJNRTqYGb347edJyv/EWWfsjeHic99V0h6raSem/3XzczsSTO7Ak+R0kgZnS0l9UvaSNK6ozMNQ90eobkG2qN4Ls7TaG5WaOFFA623zQMemmD5APB5fOrv3cCNktbJI7AiMrMbgcMmedrZecRScP9g4l6yI4ArgbmSPpVPSKFeaQb8Wxt4ydrAnenPLOCPkiZKGh1SkrYF9qXxtsJDZraGma1uZgfFmNrmRAOtt52Dn5xqMcZmt98Y/3YaMmJm3we+McFTei7reacxsyV4ks96HC/pLVnGE5ryf9Q/NmoFYN2y/78G+FHLI+oxaQ/ypTSXR27NNOFtmIJooPUoSYP4LbWJxpiJsZ+BAeDlcWsnc1+k9q3nvfMMpKjMbBbwU+q7VbZzxuGExv2Y+m9zDjC2VNQgEI3uyW0FrNTkawVs38JYCqkQedBKpdIh+D3wzYHzkiTZs2zZZnhB3i3wTL8HJUlyVTvibLGR9KeZfRxlhrK1DbUbzltKWtXMnsgzoIL6CLAe8KZJnldPepSQk3QM2VenuJo4x02u3pnntWwPXNCiWAqpKD1oD+IH9KnlD5ZKpQEgAc7DK9F/E7igVCqtnHuELZaO0/gJnlerXgvxWYQxXiBbEyWs7QO+EL2Y2UuPkcm+wBhwZg7hhPqtiyefbWQGenlOyEVEObV6XA/cQvMzOPeUtIOkmS2MqVAK0UBLkuS3SZKcj88oKbcjnlTvO0mSLEiS5Cx8wPy7cg4xK4cCp+CTBerxBJ73JmRrNWp/M+3Hcz7NkvTq/EIqHkl7Aq+b5GnnmdnNecQT6tbMl5fRhvjT+Jf1r7QunN5kZouAnWg+TcmLgKuB+yRt17LACqQQDbQJbAbcnCRJ+WDTG9PHu56ZLTSzTwIvpb7xGqsC3842qgBczuQ9myvg5VGiVmoG0nqBZzH5LZy3SerJRKdd7C6aT3q6HPCrdKJImNx7gameg1bCz2VRPq1BRW+gLc/4bwdzgZ7qkjWz2cBNdTx1kMnH44Sp+x5+a30iwk9sm2QfTiG9ivpqnQ4AH5M0PeN4Qp3M7DngviZfvgTf96E+b2TqY9X78BqdG049nGIpegPtGcammQDvPq/3lmA3uZP6xhJE4fSMmdli4D3AZ+t4euyPbMxj7My+WvqAw4FrIndWRxlu8nUDRE3IRsylejqTZsalxbmsQUVvoN0CbF4qlcrfh63Sx3vN9/EDbaIDawkT5+gKLZLWQD2AyffHHflEVDiX4mMu6zENL/b85ezCCQ36Ds0PXj+klYH0uJPx97nyvW50dufv8CTBoQFFSbPRj/+t/UBfqVRaBh+TdQU+luGwUql0IrA78GJ8VmdPMbO/SdoJ+BawNeP3vQHfM7Mzcg+ugNKyWhsz8YmuD+/hjRp2LWZmz0r6NXAQ9X1RHQRekW1UoV5mdqmkEnASsA6NXcvWlTTdzCJ9yiTM7FpJb8QbxGvh73Uz9km/lIYGFKUH7UhgPp4g9D3p76cmSbIIKOHFd+emy3dLkmROuwLN2LX4LKZqJzMB60SdunykKR7mTvY04rZAlu6hsTQ0hfhC2y3M7EIzexGNJxJ+jsiDVjczu9LMtsPTmzR7PopjpwmFeNOSJDkWOLbGspspTsbj44H/mWD57sD+eIb1kL3D8Pe6VqP4jLQhF7JxGvAJ4IXU92V1W0kvM7NG6kCG7O2Bf5mp98vlwujNaZyZmaTDaK5M1l+Bl7Q4pJ5XlB604HbGb9XUMg2ftRNyYGanA+8Anq22GPhZrgEVjJnNxW9bPlrnSxbidRxDZ9mZxsZERRLoJpnZj2kuxckGrY6lCKKBViyTdU8Lr8U5WfLOMEWSXi7pJGA3qt9uERDFhjNmZo9T/6SgQeDtktad9JkhT43edpOk30naIpNoel8zWQ4k6XJJkeKkAdFAK5Z/1fGclwCXS9oj62CKKh10+3fgo/gt5VqlxX4hKQamZ+/BOp/Xh/fW3CQpegQ6x7dpfEbn24Dro1pHYyQN4GO4G34pXrnnaklvbWlQPSwaaAUhaTlgz0mfmD4d+GGG4RTdD/DbyZONARXw8+zDKa60SsA+DbykH89GH6WCOsflNFfUexpRk7NRe+CzOZslYoxz3aKBVhxrUF9izlGrpN+WQusNUf8FZY0sAwkMNfGafrzOYOgMzaZ+AJ+ZGOr3QuorGziRVVsRSBFEA604HsSnl9frgbRYbmi9O6n/JNerKV86xb00fsFZSH2l00I+7mEKjQZJX5a0hqTlJe0v6XOS3tDC+HrJ7Uw9+8NT6Xv8lkjrNLFooBWEmT2P38oZob4cQGdmG1GhfRCfuVlPLciNJJ2QcTyFZWZPAB/BqzY8z+THxiK8zNAXMw4t1MnMHgUOpnpJosmsAhyFF2C/BR/acSzwB0lfa1WMPeQC4Nf4uavZgvWr4lU5EuBH0UirLRpoBWJmvwG2ob7knJ+XVGvwepgCM7sJeDX133I+VNKmGYZUaGb2U7yA9jFMfmw8CmyZzv4MHSJN/7ADcGWTq1gev905CMzAr41HSNqmNRH2hjR/3F74WLQ/Ud+XzGqm4+e/A4A3tSa63hMNtOK5Cx/kPJlpxDibLE2nsYHNG2cVSPCSNmb2bWDZSZ66qpk1k2YgZMzMrqH+GbnVVB6PC/EarKGMufPwnuRGxjVXs4B4j2uKBlrxzKf+sWhTOdmFiTX63g5nEkWoNNmx8UguUYRm3d3CdQ0Qx91Ehmm+B21UvMcTiAZawaRd1AfV8dQbgGUyDqewzOwh6p9ufh/wVknfkfQBSXHcZkDSlsA/J3nampLmSPqVpB+l++SDkqbakxBa44c0X2dzCd7gsPTfu4DXS3pZi2LrNScBT9J4DrpRS/Dj7YKWRdRjClGLM4xlZmdI2hEfrF7L1sBtkrYzs5vziaxwDsRnRX0Yv508jeq3Pddn6fgo4Y2190U9wdaRtAtwIZPfspme/rw3/f8IfqHZTdI7zayZgeqhBSStg1/wm/0C05f+3Aq8FD/uDsfHou1sZpe3Is4eciCwOv75b+YLSh+wJf4+39W6sHpHfBMvIEmDwL51PHUZ4JSMwymsdCzHccDXqa/Y8wD+perdwE4Zh1c0p9PcRWYavl/eApRaGlFo1DfwWZlT7Xh4GX4sDqY//fjnI6QkrQ18E3+fptJ7PAOIWeo1RAOtmF5A/ft+oywDCYAXEm6k52Uh8OKMYikcSf1Mve7pImKftNsmeGO51QSsE+kgxlif5qo3VPOSFq2n50QDrZgepf5xA9H1nL17aOxYHMQTrIYWMLPFTH3w/wCxT9rtdupLIdQoA2bHkIIxZtH82LNK/2nRenpONNAKyMwW4skZJ7MATwAZsnUWXjx9IROf9BanP+cBl+UQV5HsT3PZ6EfwRsGlxGDndvsCPmh9KjMLDfgO3qO9MP1ZjOfrCikzG8bf79H3qdnG2nzgU62Kq9fEJIGCkTQd2A9PzPi/eOLaDfHxZtNYOlB2EZ7t/nhJx6XPWQe4DTjbzKZajy2kzGyRpKOBz+LjX4ZYOq5jId5QHsAzd9+B9xQcI+kiM/tHG0LuOWb2e0nb4pMx3sDSfGijYwP70t9HL0TP4r1u84DZwDdigkB7mdkDkrbAG1Pb48mHVyt7Sq1bclb275P4dXFfYDO8cXZOTJSq6mRgPXxC2YPA2sAW+LVkstufo9eP6cA/JJ0GHBbH0FjRQCsQScsAf8ZnzhhLe1CnMb43dRp+oL0p/VmMH1R9wHskvSMOptaQ9GG8sVy+T0aNzhoEH1D7yvRnAXBUOpvz3Lxi7XEb4AP9a11cyh9fIf0B2Bx4m6T3m9mvMowvTMLMHpH0DeB3+AzDeqjs31XxHp1FwJvN7IqWB9kDJK0IXIuPRTNgOxobk1Y+sWAF4JPAzpJeFreSl4pbnMXyUbxxNohf9AfSn3o+B/1lr9kFeFdGMRZKeqL7AUt7aeo1PX3+6ekg9zAFaR6zn9HcwOfRY+g0SVkMUg+N2Rl46xTXMYB/HkJ1n8cbZ6PXklZMGNiUiVM/FU400IplQ6ZemgO8N23DFqwneP2/qTSwlsVn5YapWY36SqBNZAawRgtiCVPTqnPTupEUuqaN8cZZq0Xt0zLx4SuW+2huIHSl/nRdYepm01iKjUoLgMdaFEuRPYGP8ZuKhfgM6dBe97VoPQ/GMI6a7qb5ig0TibF+ZaKBViw/BO5k6cykRSzNhD6Z0TFoC4GrgF9nFGOhmNmTwGGMHYBej0X4fvuImWWRWqBQ0lQbH6a52Wij5YEOSmdIh/b6PXDlFNexBM+UH6r7FvAw/gVx9PM/VbOAU1uwnp4RY1cKxMyelbQ9njpjPTzH2e+BfYBX47d5lsUHyg7gDbLbgePxUkTrpP//cXpBCy1gZidIuh0fO7Md3s0/HX//n8THdywDPAXchOcNWgT8zsymeiEKKTM7U9Is4BP4vhi9hTM6s/lZ/Bv+Gvi58wZgLjAH3xdX5B1zGM/Mlkh6I/BF4GPAivj5TCwdKzX6pdTw42xa+tgzwPnAyWZ2fZ5xdxMze0LSVnhd57WAfwM7ALsCK+HHS61xaSN4gfQ18f2yADgXOCB6LMeKBlrBmNlzwHEVD9eTEy1kyMwuAS5pdxxFZ2Z/xmc6hy6WpgH6cvoTMpD2/n+97KEftiuWXhW3OAtG7qOS7pf0tKSbJX1S0hWSHpf0kKSnJC2RZOm/cyXNkbQo/f8SSc9L+vrkWwz1ktQn6TOS/lqxDyp/bpF0u6R5ku6TtEe7Y+816b44RtL8Gvug2s+SdL/dIOkfkk6SNLPdf0vRSNpI0rmSbpT06ATH0WQ/89P9+YikH8a+rE3SmunnfmQK7/foteYrUVbLRQ9a8RwBfJWl3c+bAd+b4PnCbxFUmg4cIWmGmUUm6NY4EfgIkx+XLyv7fXngHEnTzezMzCIrnpNovIqG8JxOW6X/3wLYQdKrYmxaPiStB/wTH6ox1evbMunPCvhxub2kV8a+HEvSSvjY5qk2YEevNUfiRe8/NsX1db3oQSuQdMr4sbSuyC3AIS1cV2FJGj0hNXtR+XYLwym0NDddK0qcDeJ5B3dqwbpCfQ7B0520uvNBxL6s5f1MvXFW6aDIKRgNtKJZBh+U2UrRC9saK0/x9StM/pRQp5VauK7F+KSbkI9VyO6cNELsy2rqrdjQCNH6Rl/XiQZagaQTBO5t8WqfaPH6iuoBfDZgs9PVY8ZZ6wwDT7doXQPEvsnTP5hasfSJ9BP7spq/Z7DOp/AZ7IUWDbTi2Q0v8FxuAc01DEbwOp1hitJxLe/EUzk0mkz4UWCvlgdVUGkKmbcwtQv9aJ66j5vZLS0JLNTjx8BvWJo+o1UMOCT25XjpDPTTWrjK0Tqoha/JGbenCsbMbpK0Pl5Lc0O84O2leA6brfEGwkLgPcCL8cGf5+GD0TfHywoth+eA+pqZRQ9ai5jZnyVtBLwKzzn3brze3XN4d/8M/H3/PL7vXg7cBpxrZs+2I+ZeZWZXS1oHOBQvoL4cPoB5On6xfgbfL6uk/38cuB8v0n07vr9uMrM78o++uNIcaHvhkzzWwTsh3o/nrusDXoJPIADfb4vwL6h9+PVwBO89vQ24In3sSeDS2Je1mdkBks7Aj5c18eNlA/ycJap3ACzB3//n8LsHd+FJ0H9iZlEdhWigFZKZzQF+UvHwH9KfUZHRuQ3M7GG8QQx+kanlOuCc7CMqLjN7BPhC+hO6RNrz8teyh37VrliKJE2aHYmzWyhucYYQQgghdJhooIUQQgghdJhooIUQQgghdJhooAUkLSfpB2k5p8rSG7XKdowuG5G0sGLZ1ZL6Ja0n6W9pWahnJf2/KJdSm6TXSbpK0qz0PWu0VMrTaXmbD42WSpHbPy3D8m9JX5c0OFksRSYv8/T19DO7RF7irJnSNaM/z0q6VdKnJR0gL6m2UNKDkt6SbnN0398p6QxJkW9riuQlny6S9MQU91/lMba3pKPkJdf+JemQ9DibLuk4SbdJul7S3u1+D/ImaX1Jv2vy/DVRCag7JK3Z7r8vbzFJoODk1QUuAl4DTKtcPNFLy/6tbOjvgM82XIexST/fDWwgabs0lUFISdoB+BP+XjZb6WEmnu38f/HZUycBB6X/ju6jDYH18JltoboTGVtmZqrnyWWBl+LVHsqPsbWAiyR9HDiBpft+PeDlkl5uZs9PcduFJGkIn6E+k/HntamYCZyJz/YcXe/38ETR2+LpWUa/AP1c0qCZtTIFRceStDpwDZ64tpXVaoTPvr1T0mpFKrUVPWhha+B1tPYkBl7js1oNz63xE1kY6zP4iagVJ7ZpeD07gKMYe5wPAHtJWqMF2+k5ae9iK8o8VVPrGDuGsft+EL8gRY7B5u2Lf0lp9XltVPl6+/GZvruxtHEGftwdSXHshZ/zsyp0PhMoVK9kNNDCCjSfvb4ZRvWGW9GtQmuPx+Ur/q0U+6C60bxNeW+zct+PEPtoKlYg3/04Y4I4iiKPv7VQXyyjgRZuBrK4jTJC9RPkYuDGDLbX7f6EJwhuhUXA39Lfr07/P2oJ8BhwX4u21VPM7Cm87FYmq6/x+I2M3UfgvTL/yCiOIvgL2Q7hKd+Xi/G8hHMrHl+IJ14tiqsZ24OYhQszXn9HiQZawZnZ48DuwPwWr3pX/PZauRFg9zQZaxjrm8Al6e9T7dG8G9gn/X0/4D9l630KeHuRxnE0YWdaV4uz3DDjG3/XArsAF6f/X4Jf8PczszsziKEQzOxC4OsZrX4W/iXH0p/7gT2At+Nl9Jakz7sdODCjGDqOmf0f8EWyuyPzFTO7OaN1d6SYJBAws4slrQdshZeAehV+MbkY2Bgf6D+Idy8vi5fk+Cc+mPlBvEHwTuBF+Lf+L6cloP4g6Vf4BW8+cFGanT1UMLNFknbDB/mvjt+a3A3YAr998gR+vA7hF/A5wE34xWJrvBF2GV6C6J9mtiBd78OStgZegZcpusHM5ub4p3UdM7stHWT+NmA7/L3dDNgJ/5Jh+D5ZFh+L9AzeAzYdbwDPx/fhfLxn9I94uaB/4r0qbwZelv7/yrQ8Ufm+/7eZDefyx/YwMztKXn5oA/w9fxO+f67B9+eW+G3k8vFkoz3/S/Cyd4vwRtdsvNf5QuAC/PjcJn3eP81sPnCfpBfh59H56eOVPaM9zcy+Iekc/L05BD9/LcvSMZaVk6AqG3PP4O/pPODf+PFyB/BzM7s12+g7TzTQAgBp7bPKck+N+HWN9d6FN+jCJNISNeW3f8+r9dwG17uQsaVvwiTM7Bm8RFAWZYJ+n/6Ub69y34cWMLP/4A3sSyZ7boPmUaWsUVpG7/9avK2uYmb3APdQ45oQ6he3OEMIIYQQOkw00EIIIYQQOkw00EIIIYQQOkw00EIIIYQQOkw00EIIIYQQOkw00EIIIYQQOkw00EIIIYQQOkw00EIIIYQQOkw00EIIIYQQOkxUEgBKpdIheM3CzYHzkiTZs70RhRBCCKHIooHmHgS+itdqW63NsbSNpD5gbeA5vP7mS/H6grsBz+O16K7Ba829GbgFr5U2iNeFHMbrdW4FPAzcmy6fhhefXj197hy8tt0L8YLCT+D18u4HHk+3uzVwPV5SBeDpyWpISloeWBV4sBdq4ElaEa/ZuAivT7crS2sHLsFrofbh7/vj+PH8Qnyf/Qe4OX3+fOARYFvgX/h+ucvMFqfb6UtftxHwHzObJWk5fH892GxhdUkC1sTrGz6WljPqGpKm4cfDDODFeJ1AgNcAA3i9wWeAtfDP9WP4l7zH08cFrJz+Ows/Nubh+6cfmIsfC3fhdR+3wY+Lq/AahY/jx8896XqH0vXekG7vcTN7VlJ/GudcM5uwyLuk6Wm8j6T1Iwsh/SyuhZ8f5uN1VjfAj6E18f0y+pnvB1bA3/vF+Lnrebz+5hz88/wA8HdgJvBqvNZqAlyKH0urp6+/xcyeSD9LGwGb4OfRm81sRNKywIbp9meZ2Wih9a6UnkvWZ+kxcQP+Wf4gsD1+PlsZr8/Zj5/L+tLnLGZpm2Rx+q/h578+/P1chF+frsOPo+XwkltX4tfxacC6+LG4IH39dGDYzJ6SNAM/Dh8arVdcft0zszllx/3TZvZUS9+gBkUDDUiS5LcApVJpKwraQJO0CXARfiGq104ZhDLC2OLF/5UWPv5QZeMrPfkeDRyDXwyfkrSHmV2WQXyZS08QpwMfyHAz8yW9GW+4XYmfUEe3/yR+gZoGPCtpLzNLGlm5pLWB3+GNDoArJL3rrru6oyyrpM3w42HdJl6+UZXHtqnyWC2fr/N5I5L+F3gP8AIASScDh1a70Et6D/BzvMG5WNKnzOzkBuLqSmnh+z8Amza5irWrPPYqYI+Kx95bY/t/wxvYy5c9vEDSj/GC4qPFwx+W9PpuOUYqSdoI/3KxRg6b27js913qeYGkq/FGYj9+/vsg/oX1v9c9SZfhBd7XTP//I+Bj7donMQYtjH6rvgz/NtluVRtnqfcBX67y+AeAo1h6olsRSCS9qMWx5eVwYO+MtzED+CNe2HmtimUrs3Q/LAf8RlLdF7e0wXw+flEatQNwZtPR5ijt1fgD3hPSyaYBH8N7a0Z9FDis8omStsILv89IH+oHTpS0a9ZBtlP6WfwdzTfOWuFVjG2cgffqfJyl5yzwRsHls2fPHsgrsFaRNABcQT6Ns2a9mqWdUjOAX+Axl1/33szYv2F/4At5BFdNNNAC+MnrhS00IFcAABD3SURBVEzcOOoEA8C7qzy+G9Vjf3224WRmT8aeuLMyCKxTx/MW47f/67Uyfju1/EIzCOwyb968bjjnjH6DzmMftEJ5nP14j1qlXfBbROWWAG/PKqgOsSo+XKJbrPmLX/xi/XYH0YQNqd7T2MkW4Y2xymtH+fFU65qTi244WYbsddN4rcVVHluIjzUopxrP7QYL2h1AhUbfy1qfpyXLLLNMN4yx6abjoZpqYwYXMb7BOTq+p5d13d83ODg40u4YmtB17zP1fwFr23UkGmgBfKD+9XT+QTYC/LDK46dXed4zwCWZR5SNUxjf4MzCk/gg3omM4Bf839W7UjObB1zA2IbCQuD0gYGuuHtzE3Ar/rd3g/I4lwA/qvKc89Jl5Q3kPuDsDONqu3SQ94Xkczy1wi0HHHDAfe0Oogn34OO5usUS/Lj5N+OveyMVv1e75uQiGmhAqVTqL5VKy+C3B/pKpdIypVKpK64krWBmI/gMwT/hvTd5f2MYPXmO4DPbqi1/CjgSOGncQrNLgffjM0cX4xfX15nZo5lEm73T8YHiWTaY7wVeiY+5uKpi2b+Ah/D38k5gRzOb3eD69wZ+jc+Yexb4KfCJqQScl3QSyk74+JRO7fFbgu+jg4Fr8c/KE8AnzGzcWD8zuwd4Az5jdDE+W/TtZnZtbhG3z174mMh2NLiXAKfis6rLDQMfwWeIgp/j/ga8cebMmZ36maspnZTyRnxma3ljuFMaxvOBH+OzbxfjDco34sNg/g//AvkU8D18PyzCZ+x+Bj93tUXM4nRH4jMAR43OdtqvLdG0QdqY6doBw2b2S+CX7Y6jFdJ0FN9Jf/Lw2lav0MyewRvNY9x9992t3lQmzOwhGht31051fcM3s78zdvZbIaQ9uu9qdxw1/LjygbvvvruecaEdx8yewCdEdJtJZ4G267wVDTQgSZJjgWPbHEYIIYQQAhC3OEMIIYQQOk400EIIIYQQOkzc4szfpt0yDqdH1JOgMvZJfmJ/dJ7YJ50l9kfnaUui42ig5edxvIbYWe0OpICew9//SrFP2iP2R+eJfdJZYn90nlr7JDPqsvrFXe3uu+9el4LW+myzxzfYYIP7qy2IfdIWsT86T+yTzhL7o/PU3CdZiQZaCCGEEEKHiUkCIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdpr/dAQRXKpUOAfYDNgfOS5Jkz7JlmwE/AbYA7gEOSpLkqnZtL33tEcAKwO+BA5MkeTrvbZVKpZWAHwO7Ak8DX0uS5JTm3pHWxZz1+ut5/7vxPa83/kle93rgaGAbYH6SJGtmsI3PAh8A1geeAs4Ajk6SZKSF2/gI8FlgdeB54GLgE9WOs2a3UbGOy4EdgRlJkjxf52syP2flcZ7qxmMly+Mjq+Miq2OhXZ//PEQPWud4EPgqcGr5g6VSaQBIgPOAlYFvAheUSqWV27G9Uqm0E3AM8FZgCJgOnNSmbZ2Mf8lYG3gb8JX0BNQKWe+PrN//bnzPJ42/Ds8CpwGfznAbfcD+wKrAa/D34TMt3sZlwHZJkqwIbAgMAt9o8TYAKJVK+zXzulrbbfE5K4/zVDceK1keH1kdF1kdC+36/GcuetA6RJIkvwUolUpbAauVLdoRWBb4TpIkS4CzSqXSp4B3AT9tw/b2A05PkuTG9PVHAf8olUoHJUnyXF7bAgS8B9g6SZJ5wA2lUuln+Ani8mbflxbEnPX696OO978b3/NyE8RPqVTaFvgusBnwMPDF0ecnSXItcG2pVNoxw218q+yps0ql0tnA/wDlj091G/eWPVXAEmCjVv4d6fJVgS8A7weurbb+WvI4Z+VxnurGYyXL4yOr4yKrY6Fdn/88RA9a59sMuDk9QYy6MX28HdvbDPhX2bJbgGnUuHhkuK2XAEqS5N81XpuVrPdH1u9/N77n/1UqldYCLgGOx0/G+wE/KZVKm7ZxG6/D36eWbqNUKr21VCo9hd8ueid+oWn133E88D3gsUbWPYk8zll5nKe67ljJ8vjI6rjI6lho4+e/ZaKB1vmWxz+U5eYCM9u0vTHLkyQxfHxFM/FMZVvLp7/Xem1Wst4fWb//3fiel/sA8MckSc5PkmQkSZJr8FtQ72nHNkql0sfxsS/HtXobSZJclN7WWQ+/iNzTym2kvSgvBX7U4Honk8c5K4/zVDceK1keH1kdF1kdC+36/LdM3OLsfM/gA1DLrQjMa9P2qi1focl4prKtvklem5Ws90fW7383vufl1gfeUSqV5pY91g+cmfc2SqXS3vgA8TckSfJEFtsASJLk/lKpdCnwS+DlrdhGqVQaBH4A7JskyZJSqdRg+BPK45yVx3mqG4+V9cnu+Khr3U0cF3WtFxo+FmquN+PPf8tEA63z3QJ8rlQq9ZV1tW8F/LBN27sF2BL4Bfx3ltMI8J+ctyXASqXSpkmS3Fb22oZuNWUQc9brn+r7343vebn7gV8lSbJfO7dRKpX2wnsH3pQkye1ZbKNCPz5AuiXbKJVK6wObAJekF6fRuymzS6XSB5IkubjBbZXL45yVx3mqG4+VLI+PrI6LrI6Fdn3+W0Zm1u4YAlAqlfrxD96R+PiFffCDHeBO4BTgRGB3fHbQhkmSzMl7e+nMpbOBNwH34lOpn06SZN+8t5UOQp0OfBB4MfAnYI8kSf6v2fdlqjFnvf563/9ufM/rjH8N4DrgQ/h0+z5g6zTG20qlUh8+y+u1aezrA5YkyYIWbuN9+Hv35iRJbsjo7/ggcEmSJA+VSqUX470JDydJsnsrtoF/BtYqW806wN/wC99wUkeqgTzOWXmcp7rxWMny+MjquMjqWGjX5z8PMQatcxwJzAe+iN8jnw+cmiTJIqAEvBsfw/BFYLepNM6msr0kSf4AfAUffPkQsAj4eJu29THA0mUX4zl3WtVQyHp/ZP3+d+N7Xk/8s/G0Bp8EHknj+AZ+IQS/8MwHLsVP0POBO1q8ja8DKwFXlUqlZ9KfWt+4m93GK4DrS6XSs8CVwM3Aga3aRjomZ/boDz7DDRq7OOVxzsrjPNWNx0qWx0dWx0VWx0K7Pv+Zix60EEIIIYQOEz1oIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYQQQggdJhpoIYTMSdpE0t8lPS/pxnbH000k3Sfpk+2OI4SQr2ighRD+S9LqkhZJWlZSv6RnJa3bglV/CXgW2Bh44wTbX1PSSZLukbRA0gOSfiep5muKTtKx5Y1eST+TdH47YwohTF1/uwMIIXSUVwE3mtlzkrYH5pjZ/S1Y7wbARWY2q9YTJK0PXI3XPfwccBMwAOwM/ADYpAVxhBBCV4getBBCuR3wRhLAa8p+r0lSn6SjJc1Oe71ulLRL2XIDXg4cLckkHVtjVafgxaW3M7Nfm9mdZnarmX0XeGXZ+taVdIGkZyQ9LelcSWuULT82jeED6e3BpyT9StLMsue8W9LNkuZLekLSHyUtly67QtIJFX/j+ZJ+Vvb/+yQdKemMNI5Zkt6R9kCOxnazpFdUrGcHSX9Ot/uApBNHt5suf0HaYzhf0r2S3j/Z+1+x/mOBfYF3pO+1SdoxXTYk6RxJT6Z/8wVpo3j0tT9L/84vSHpE0lxJx6Q9qd+RNCfdx/uXvWZQ0smSHkpvX98n6YhGYg4hVBcNtBAKLm3wzJU0F/g08JH0968Du6XLTplgFYcChwGfAbYALgUSSRuly9cCbgWOT38/rkoMqwC7AD8ws2crl5vZ3PR5As4HVgFeB+yE986dU/GSDYDdgLelP68DDk/XsRbwS+A0YFNgR+C3gCb4G6v5FN6A3Rq4CDgTOAM4C9gGuAs4I40ZSZvj781v8ffpvXgj+OSydf4MWB94A/Bu4GDgBQ3EdBxwLnAJ/l6vxf9v595CrKriOI5/f6NGBZWBUUSlUmpqlGg13YNIKQXB7PbQQwU+VCRWD0IvDV2gC0IagWaEFT5UIApRQXYzupBdLVMzalK0vKSkQ8Vk/Xv4r6O705lxzjgPB+b3gc3svc7a67JnmPNnXTZ8KOlY4B2gC7ii1NsFvCHpqMr9VwGnljz3Ah3Aq8BeoB1YDCyWdHrJPxeYCdxITl/fAnQ20V4z60lE+PDhYxAf5FKHUWTQ0A2cRwY4+8kv6lHAiF7u3wbcX5f2CRls1a6/BDp6KeNCcvRs1mHaOhU4AJxeSZtQ7r2gXHeQ692Oq+R5HPi4nE8u+Uf2UMe7wJN1aSuBZZXrTuDFyvUppcwHK2kXlbRTyvULwJK6ci8D/gaOBsaW/O2Vz88uafN6eSYd5LR07XoZsLIuz+3ARkCVtKOA34Fplfs6gbZKno3Amsr1EDKwu7lcLwLeqpbrw4ePgTk8gmY2yEXEgYjoJIOBtRHxFRlw7IiINRHRGRG7G90r6XhyxKV+KvQDcnSqr2qjV3GYfOOBrRGxtdL+b8l1a9X6OiNif+X6Zw6NRH1FBhVfS3pF0hxJJzbR1pp1lfMd5efXDdJq9U4Bbi3Tn12SusgRtTZgdGn/AeDTSt82lr4dqSnAWcD+St17yMDwzEq+9RHxT10fDvYpIv4Gfq30aRkwCdhUpmunDUBbzQxvEjAb9CStB0aSC/Lbypf3UGBoOf8pIiYeppj6wEoN0nqzueQfT45W9djcHsqtT/+rQfvaIIMMSVPJ9XbTgLuBRyS1R8SPwD/8f7pzWIM6D9YREVFmMqv11trTVvm5hBx1qreFnCKs3jeQ2oDPgEZr2nZVzhs9t96e5eeSRgPXAlcDL0taHRHXD0irzQYxj6CZ2XRyFOQXcg3RJOAbYF45n97TjRGxD9hOTtVVXQJs6GsDImIPOZp0V3XRfI2k4eX0W+CMyhooJE0ATmiyvoiIDyLiAXINWTcwq3y8i1y7VSt/CHBOX8vuxefAxIj4vsHRXdo/FDi4sUDSOGB4D+X1pJuciqyvewyws0Hdv/W/S/k3EBEvRcQccl3d7LKm0MyOgAM0s0Eu8tUXXcDJwCpyNGcCsKJ8gff4aoziCWC+pJskjZP0KBnYLWyyKXeSgcUnkmZLGiNpvKS5wEclz2pyanG5pMmSLiTXdr0XEZ82Lva/JLWXnYrnK9/xdh1wEocCvLeBGZJmSDqb3F3abJDUyGPAxZKeljSp9G+mpKcAImITubh/aWnjFOBZ4I8m6+kEzi2/ixGShgHLgd3AKkmXSxot6UpJCyWd1t8OSbpH0s3KFxGPBW4gA/2BmJY1G9QcoJkZ5E7GtRHxJ7lbb1tEbO/jvYvIHZoLyPVK1wAzI2JzMw0o04uTyd2GC8hRvDfJF9veUfIEuTtzL7CGDNh+IEdu+mofufnhNeA74GHgvoh4vXz+HPA8JfADfixtOiIRsY7cTToGeB/4AniIXB9XcxuwtdS7AngG2NlkVUuBTeRatl3ApRHxO9nnLaXcDWQ/jyGfR391AfNLXWvJDSXT69axmVk/KP/fmZmZmVmr8AiamZmZWYtxgGZmZmbWYhygmZmZmbUYB2hmZmZmLcYBmpmZmVmLcYBmZmZm1mIcoJmZmZm1GAdoZmZmZi3GAZqZmZlZi3GAZmZmZtZiHKCZmZmZtRgHaGZmZmYt5l96OshvbooFOgAAAABJRU5ErkJggg==\n",
"text/plain": [
"