[d869a2]: / Notebooks / 04_first_subject_data_analysis_test.ipynb

Download this file

189 lines (188 with data), 52.0 kB

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "#%cd ..\n",
    "\n",
    "from wearsed.dataset.WearSEDDataset import WearSEDDataset\n",
    "from collections import Counter\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "dataset = WearSEDDataset()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "subject_infos = {\n",
    "    'age':         {i: 0 for i in range(-1, 100)}, \n",
    "    'bmi':         {i: 0 for i in range(-1, 100)},\n",
    "    'sex':         [],\n",
    "    'race':        [],\n",
    "    'cur_smoker':  [],\n",
    "    'ever_smoked': []\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "def int_or_NaN(n):\n",
    "    try:\n",
    "        n = int(n)\n",
    "        if n >= 100:\n",
    "            print('Large')\n",
    "        return n\n",
    "    except:\n",
    "        return -1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "try:\n",
    "    for idx, recording in enumerate(dataset):\n",
    "        subject_info = recording.get_subject_info()\n",
    "        subject_infos['age'][int_or_NaN(subject_info['age'])] += 1\n",
    "        subject_infos['bmi'][int_or_NaN(subject_info['bmi'])] += 1\n",
    "        subject_infos['sex'].append(subject_info['sex'])\n",
    "        subject_infos['race'].append(subject_info['race'])\n",
    "        subject_infos['cur_smoker'].append(subject_info['cur_smoker'])\n",
    "        subject_infos['ever_smoked'].append(subject_info['ever_smoked'])\n",
    "except:\n",
    "    pass"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Axes: >"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGdCAYAAAA44ojeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABKlklEQVR4nO3deXSb5bU/+u8ryZI8yY7t2LITxwlDJhJCmgkXSmmTH0ng0HLI7So0pYGTCxeO01OacwpNy1Dg0PS0XaUHfinc9rbQ/kpKy+8ylFyaNiQQoHUGDAFCIBMBZ7Ds2I4t27Lm5/4hPa8kW6MtWcP7/azl1Vh6bT15i+ydvfezH0UIIUBERESUQ3TZXgARERHRSAxQiIiIKOcwQCEiIqKcwwCFiIiIcg4DFCIiIso5DFCIiIgo5zBAISIiopzDAIWIiIhyjiHbCxgLv9+PM2fOoLy8HIqiZHs5RERElAQhBAYGBtDQ0ACdLn6OJC8DlDNnzqCxsTHbyyAiIqIxOHnyJKZOnRr3mrwMUMrLywEE/oIWiyXLqyEiIqJk2O12NDY2qr/H48nLAEWWdSwWCwMUIiKiPJNMewabZImIiCjnMEAhIiKinMMAhYiIiHIOAxQiIiLKOQxQiIiIKOcwQCEiIqKcwwCFiIiIcg4DFCIiIso5DFCIiIgo5zBAISIiopzDAIWIiIhyDgMUIiIiyjkMUIiIcsQHZ/qxdW87hBDZXgpR1uXlacZERIVGCIE7fv822nsdmFlXhsXTq7K9JKKsYgaFiCgHHO4cQHuvAwBwum84y6shyj4GKEREOWDnh13qn3uH3FlcCVFuYIBCRJQDdhzqVP/MAIWIAQoRUdadHXDh3VN96ucMUIgYoBARZd2rH3UhfOMOAxQiBihERFm348NAeefC2jIAQA8DFCIGKERE2eT0+PDm0W4AwFeXNAIAzjFAIWKAQkSUTa3HezDs8aG+wozLLqgBwBIPEcAAhYgoq2R5Z/mcWlSXGgEA5xxu+P2cJkvaxgCFiChLhBDYFZx/snxOHSYFAxS/APqHPdlcGlHWMUAhIsqSD87YYbM7UWLUo/m8ahTpdSg3B04gYaMsaR0DFCKiLHklWN753IU1MBfpASCizEOkZQxQiIiy5BW1/6ROfawqGKD0DDJAIW1jgEJEFMdrh7uw5dVjECK9Taunzjlw8LQdigJ8cXat+rgMULiTh7TOkO0FEBHlsgdeOoQT3UP43IU1uHhqZVq+55HOAaz/7X4AwOKmSagpM6nPVbHEQwSAGRQiorj6goHCmT5nWr7fro86cf0v/oGTvcNorCrG5usvjnh+Eks8RACYQSEiisvh9gEAzg66xvV9hBD41RsfY/NfPoIQwLIZVXj864vUjInEJlmiAAYoREQx+P0CLq8fANA9ML4A5ZEdR/DormMAgBuXNuKBL82D0TA6iV1VGij3cJsxaR0DFCKiGIY9PvXP482g/OmtUwCAu1bNwh2fPx+KokS9rqq0CADQOzS+1yPKd+xBISKKQZZ3AODsODIobq8fnQOBHpb/Y9HUmMEJEMqgnBviJFnStpQClM2bN2PJkiUoLy9HbW0trrvuOhw+fDjimiuvvBKKokR83H777RHXtLe345prrkFJSQlqa2vxne98B16vd/x/GyKiNHJ60hOg2PqdEAIwGnSoKTXFvbaqJNgkywwKaVxKJZ7du3ejpaUFS5Ysgdfrxfe+9z1cddVVOHToEEpLS9Xrbr31Vjz44IPq5yUlJeqffT4frrnmGlitVvzjH/9AR0cHvvGNb6CoqAg//OEP0/BXIiJKj3RlUE71OQAAUyqLodPFzp4AQFVZIEBxevwYdvtQbNSP+XWJ8llKAcr27dsjPn/qqadQW1uLtrY2XHHFFerjJSUlsFqtUb/H3/72Nxw6dAivvPIK6urqcMkll+Chhx7C3XffjR/84AcwGo1Rv46IaKI53KHMbvegC0KIuOWZWOQW5SmVxQmvLTXqYTTo4Pb60TPkwlRjScKvISpE4+pB6e/vBwBUVVVFPP7000+jpqYG8+bNw6ZNm+BwONTnWltbMX/+fNTVhUY7r1y5Ena7HR988EHU13G5XLDb7REfRESZFt4k6/L6MeAaWyn69LlhAMkFKIqiqGUeTpMlLRvzLh6/348777wTl112GebNm6c+/rWvfQ1NTU1oaGjAe++9h7vvvhuHDx/Gc889BwCw2WwRwQkA9XObzRb1tTZv3owHHnhgrEslIhqT4bASDxAo81jMRSl/n9OyxDMpcYACBKbJ2uxOBiikaWMOUFpaWnDw4EG8+eabEY/fdttt6p/nz5+P+vp6LF++HMePH8f5558/ptfatGkTNm7cqH5ut9vR2Ng4toUTESXJESVAOX9yWcrf53Rf8hkUgOfxEAFjLPFs2LAB27Ztw6uvvoqpU6fGvXbZsmUAgGPHAgOKrFYrOjs7I66Rn8fqWzGZTLBYLBEfRESZFl7iAQJ9KGOhlnhSyKAADFBI21IKUIQQ2LBhA55//nns2rULM2bMSPg1Bw4cAADU19cDAJqbm/H++++jq6tLvWbHjh2wWCyYO3duKsshIsqoaCWeVPn9IqUmWYABChGQYoDS0tKC3//+99i6dSvKy8ths9lgs9kwPBz418Hx48fx0EMPoa2tDZ988gn+/Oc/4xvf+AauuOIKXHxx4ECsq666CnPnzsVNN92Ed999F3/9619xzz33oKWlBSZT/PkAREQTaWQGZSwBSvegC26fHzoFsFaYk/qa8QYorx85i/tePBgxx4Uo36QUoDz++OPo7+/HlVdeifr6evXjj3/8IwDAaDTilVdewVVXXYXZs2fj3//937FmzRq89NJL6vfQ6/XYtm0b9Ho9mpub8fWvfx3f+MY3IuamEBHlgmg9KKk6Few/sVrMKNIn9yN3vAHKI68cwe9aP8WrH3UlvpgoR6XUJCuEiPt8Y2Mjdu/enfD7NDU14eWXX07lpYmIJtxwcA5KZUkR+hyeMZ3Hk2r/CTD+AKV/ODAm/5MeR4IriXIXz+IhIopBlniaqgLD0sbSJCt38DQk2X8CjD9AcbgC627vZYBC+YsBChFRDLLE0xgMUMZS4jmT4hZjAKiWAYpjbAHKUDDz0947NKavJ8oFDFCIiGKQu3imqRkUN/z++KXukcZS4pkUDFD6HB54ff6UXk8IgSGXDFCYQaH8xQCFiCgGWeKRGRSfX+BcilmNVIe0AcCkEiPkkT99wX6SZLm8fsgY6kyfE54UAxyiXMEAhYgoBlniKTcb1L6QVBtlZQZlagoZFL1OQWVxYKR+qn0oQ2HnBfn8Qi0xEeUbBihERDHIEk+JUY+askCA0j2QfMDQP+xRDxhMpUkWCJV5egZTDVAit0azzEP5igEKEVEMssRTXGTA5PLAIMmzg86kv15mT6pKjSgxpnb0mWyUTbWkJBtkpU+51ZjyFAMUIqIYZAal2KjH5LJggJLCTp6x9J9Ik0qCGZRxlHgA4CQzKJSnGKAQEcXgCGYjSoz6UAYllQDlXCA4GEuAUh0sKZ1LNUBxs8RDhYEBChFRDKESzxgDlL7UtxhLYx3W5nCxxEOFgQEKEVEUfr+A0xPYolts1KMmWOLpTqFpNRslnsFggCJf82SvI+ExJUS5iAEKEVEU4ScZj73EM/YMylhLPHJr9GxrOQBgwOXFOUdqs1SIcgEDFCKiKMIDFLMhLEBJYQ7K6b7Ajp8JbZIN9s1UlRpRZwmsOdt9KK8fOYuWp98e89lCpE0MUIiIopA7eMxFOuh0irqLp3fIndR0VqfHpx4umMqQNqm6VL5eaoPh5C6eUpMBTVWlALIfoDyx+zj+v/c78OcDp7O6DsovDFCIiKJwqEPaAvNLJpUYodcF5s8nMzxNTnAtNepREZwKm4oqtcTjSamHRA5qKzHq1RH97T3ZPTSwoz+QSTrcOZjVdVB+YYBCRBRF+A4eANDpFHV4WncSZR7ZINtQWQxFHqyTgqpgicft86uNr8mQW6NLTQb1kMNsZlCEEOjoD9yLwzZ71tZB+YcBChFRFPIXfbFRrz6WSqPseBpk5evK4OjcUPJNrjKDUmrUo6k6+wFK/7BH3Q11pHOQO4ooaQxQiIiiCD+HR0opQBnHFmNJzkLpSaEPRTbJlpgMYSWe7AUosrwDBLZAn+lP/qgA0jYGKEREUYws8QAIjbtPpsQzzgwKMLZhbbJJtiysxNNhd8Ll9cX7soyR5R3piG0gK+ug/MMAhYgoCkfYOTxSTQoZlFNpzKCkFqBEnsBcYtRDiFDANNE6RmRMDncyQKHkMEAhIooiaolnDBmUsWwxlsYSoIQ3ySqKomZRPs1SH4otGKDIPmFmUChZDFCIiKIIlXgM6mPJ9qB4fX7Y7HJIW8mY16AGKI7kA5RBtUk2sG7Zh5KtU41lBuXiqZUAmEGh5DFAISKKIlTiCf2YlAFKd4IApXPABZ9foEivoDb4NWOhBigpnP8TyqAEMj9NMoOSpUZZmUH5/MzJAICjXYPwJjHojogBChFRFMNyN4wxlEGRBwYmyqDIIW31FcXQ6VKfgSKlWuLx+4UaWJWaAuueluWtxrJJdun0KhQX6eH2+rNWbqL8wgCFiCiKqLt4gtmQAZcXTk/sXTGHg30WjVVj7z8BQgFKd5IBiiNsTblQ4gkMaQtkUBoqzZhZVwaAfSiUHAYoRERRRNvFYzEbYDQEfmzGy6Ls+qgLAPDZ82vGtQarxQwA6ExydogjuMVYpwTOEAJCJZ72XseED0kbcHnV+2itMGNmXeCEZfahUDIYoBARRRFtF4+iKAl38gy7ffj7sW4AwPI5teNaQ31FIEDpGnAm1bcx5A41yMrx+lMmFUNRAgFXdwq9LOkg+08qiotQYjRgljUQoBxhgEJJYIBCRBRFtBIPkHgWypvHuuHy+jF1UjFmBTMGY1VTZoJBp8AvktvaLIe0lZhCazYZ9KgPZmImug8l1IsTeH01g8ISDyWBAQoRURTRSjxA2CyUGAHKK4c6AQAr5tSN6ZDAcDqdgrpgcDFy4Fk0MkApDWvsBcIbZSf2VGOZQZEBisygfNLjiNvDQwQwQCEiiipaiQcI22ocJaPh9wvsDPafjLe8I8lf7h19SQQoYUPawqmnGvdM7DRZGVRZKwLNwrXlJlSWFMHnFzh+dnBC10L5hwEKEVEU0Qa1AfGHtb13uh/dgy6UmQxYNqM6LeuwygClP3FwET7mPlxommx2MyiKoqhlHvahUCIMUIiIohiOVeIJBijR+jlkeefzMyeru33GS/5ytyVR4nHEyqBUlwKY+FONO+wyg2JWH5ul9qEwg0LxMUAhIorCoQ5qiwxQFjdNgqIAbxztxr4TvRHPvfJhIEBJV3kHCJVH5C/7eNQx9yMClPNqAgHK8bODE7rV2NYf2SQLADO5k4eSxACFiCiKWLt45tRbcMOSRgDAfS8eVLf/njrnwEe2AegU4Auz0hegpJRBUZtkI9d8QW0ZFAU45/AktRsoXTpGlHiA8AwKAxSKjwEKEdEIfr+A0xMIPEaWeADgrpWzUVlShI9sA/ht66cAgJ0fBppjFzdVYVJwAmw6WFMIUIbc0TMo5iI9pgfLPEcmqLQy6PJiwBkImGQWCAgFKKf7hjHg9EzIWig/MUAhIhphOGwL7MgSDwBMKjXi7lWzAQCP7DiCTrszI+UdAGgI/nLvtDvh88cvzwzFyKAAocDgI5s9reuLRQZU5SYDysICpoqSInVC7pFO9qFQbAxQiIhGCA9QzIbRv+wB4KuLG3FJYyUGXV7c88JB7P040I+yYm5dWtcyudwEvU6B1y/Qk6A8I7cZl4zIoAAT3/shdx2FN8hmay2UnxigEBGNIHfwmIt0MU8j1ukU/Od186AowI5DnXD7/JhRU4rzJ5eldS16nYLa4M6hRMPaHDGaZIGw3o8Jylqo/SeVow9MnBU8NJB9KBQPAxQiohEc6pC20b/ow82bUoGvL2tSP1+R5vKOFJqFEj9AUQe1RSvxWANBwdHOAfgTlIrSQZ2BYomSQeEsFEoCAxQiohFi7eCJ5j+umoXqYFPsqnnWjKynPslhbepZPFECq6bqUhj1OjjcPpzuy/xE2dAU2dgBytEu9qBQbPH/eUBEpEFyBkq0HTwjVZQU4U+3N+OT7iEsaqrKyHqslkCZJNFOHjlJtixKiadIr8P5tWX4sMOOw7YBNAany2ZKtBkoktxRdHbABYfbmzBTRdrEDAoR0QixzuGJ5fzJZVg+J73NseHqUyzxhJ9mHE7t/ZiA0kq8DEpFSREqiosAACd7J/Z8IMofDFCIiEaQJR5zEiWeiZDsLBTZOxMtgwKEds9MRHOqzS6HtI1ukgXCDjCMcmQAEcAAhYhoFEeKGZRMUzMo9vjZhkFX9PH80qwJak4ddvvQ5wgMYYuWQQHCDjDsmdgDDCl/MEAhIhrB6cmxACW4Vbez3xVzB47H54fbG5h+Wxqjp0M2px4/OwhPcER/JsjsSYlRD4s5+lqmVQcClJPMoFAMDFCIiEZwuHOrxFNbboKiAG6fH70Od9Rr5JqB6HNQAGDqpGKUGvXw+AQ+6c5c5iJ8SJuiRJ8jwxIPJcIAhYhohFwr8RTpdZhcFhjWFqsPRW4xLtIrMBqi/2hXFCXUh5LBMk9HX2CNDTH6T4CwEg8DFIqBAQoR0QihEk/ubH9NtJNHbo1OtGa1DyWDjbKyxBOr/wQIBSineocnZHAc5R8GKEREI8hf9rlS4gHCp8lGb5SNNwMl3Ez10MAMZlDizECR6ivMMOgUuH1+NaAhCscAhYhohFwr8QCh7bqxMihDCXbwSLMm4KA+W5wZKJJBr8PUSYG/E/tQKBoGKEREI+TaLh4g8SyUIRlUJcigyADl016HOpAu3dSDAuMEKADUabYMUCialAKUzZs3Y8mSJSgvL0dtbS2uu+46HD58OOIap9OJlpYWVFdXo6ysDGvWrEFnZ2fENe3t7bjmmmtQUlKC2tpafOc734HX6x3/34aIKA1ybRcPkPg8HplBKYsxRVaqKTOhutQIIYBjGToLR82gWGI3yQJhO3l6GKDQaCkFKLt370ZLSwv27NmDHTt2wOPx4KqrrsLQUGi72re//W289NJLePbZZ7F7926cOXMG119/vfq8z+fDNddcA7fbjX/84x/47W9/i6eeegr33Xdf+v5WRETjkMslntgZlOSaZIFQH0omdvI4PT70DAW2QifKoDRVM4NCsaXUor59+/aIz5966inU1taira0NV1xxBfr7+/HrX/8aW7duxRe/+EUAwJNPPok5c+Zgz549uPTSS/G3v/0Nhw4dwiuvvIK6ujpccskleOihh3D33XfjBz/4AYxGY/r+dkREY5CLJZ7wXTxCiFHzRRxJNskCgTJP68c9GelD6bK7AAAmgw6VJUVxr+UsFIpnXD0o/f39AICqqsAJnm1tbfB4PFixYoV6zezZszFt2jS0trYCAFpbWzF//nzU1YUO1lq5ciXsdjs++OCD8SyHiCgtcrHEU2sJzEFxef3qGPlwicbch1MzKBnYyXO6L1CCaqgsjjmkTWIPCsUz5k3+fr8fd955Jy677DLMmzcPAGCz2WA0GlFZWRlxbV1dHWw2m3pNeHAin5fPReNyueByudTP7Xb7WJdNRJRQ6DTj3JmDYjLoUVNmRPegGx39Tkwqjcw2y63RsabIhpuVwUMD23sDJX8ZfMQjMyi9Q24MOD0oN8fPuJC2jDmD0tLSgoMHD+KZZ55J53qi2rx5MyoqKtSPxsbGjL8mEWnXcA6WeICwnTxRDg2Uu3hincMTbmZdWfD7ONEfJRszHjIbMq0qfoMsAJSbi1AVDLSYRaGRxhSgbNiwAdu2bcOrr76KqVOnqo9brVa43W709fVFXN/Z2Qmr1apeM3JXj/xcXjPSpk2b0N/fr36cPHlyLMsmIkqKzEYU51CJBwjtijnTN7pRVu7iKU2wiwcIBAZTggcQxmuUffLvJ/CFn76mlm2S0d4buHZaEhmU8Ot4aCCNlFKAIoTAhg0b8Pzzz2PXrl2YMWNGxPOLFi1CUVERdu7cqT52+PBhtLe3o7m5GQDQ3NyM999/H11dXeo1O3bsgMViwdy5c6O+rslkgsViifggIsoEv1/A6Qmc9FucYxmU+jizUOQk2WTLUtNrAoHB6b7YgcEL75zGie4h7DvRk/Qa23sCJZ5pVaVJXc9GWYolpQJrS0sLtm7dihdffBHl5eVqz0hFRQWKi4tRUVGB9evXY+PGjaiqqoLFYsE3v/lNNDc349JLLwUAXHXVVZg7dy5uuukm/PjHP4bNZsM999yDlpYWmEym9P8NiYhS4PSGhpflaokn2jTZUA9KcmuuCR4+2DMY/XRkAOgOPieDn2SESjypZVA+5SwUGiGlAOXxxx8HAFx55ZURjz/55JO4+eabAQCPPPIIdDod1qxZA5fLhZUrV+IXv/iFeq1er8e2bdtwxx13oLm5GaWlpVi3bh0efPDB8f1NiIjSwBE2XdVsyK0ApT5eD4os8SSZQakuDQQo3TECFCEEeoYCmxNk8JOI3enBuWBPy7TqJAMUzkKhGFIKUIRIfOKk2WzGli1bsGXLlpjXNDU14eWXX07lpYmIJsSwusVYB50u/jbZiRbvPJ7QqPvkgqrqskBzas+gK+rzDrdPLXUNJplBkRNhq0uNSc1jAVjiodh4Fg8RUZjQDp7c2WIshfegjPwHY2jUfXLrrpEBylD0DEp46cfhSi6DIhtdk9liLMkA5fS5YXh9/qS/jgofAxQiojCyxJNrO3iAUA+Kw+2D3RkZNIROM06txBMrg9I9FHp8KMlDBT8NBihNSZZ3AMBqMcOo18HrFzFPaiZtYoBCRBRG3WKcYw2yQGCy7aTg+PjwnTxCCDWwSjaDIks8sXpQwjMoQ0lmUFJtkAUAnU7B1ODMFJZ5KBwDFCKiMLl4Dk842Ydy6lzol7nL64fXHyj5JNuDou7iGXJF7S8Mz6wk2yQ7lhIPwD4Uio4BChFRmFw8hyfc+bWBKbBHuwbVx8J3HpUkuW6ZQXF6/FFLOOG9KYNJZlDkVuGmFAOUJgYoFAUDFCKiMA53bmdQZgXH1B8JO0dHlmDMRToY9Mn9WC8xGtS/Y7Q+lO6IDEriHhSvz69OnE12i7GkHhrIWSgUhgEKEVGYXC/xzLIGJml/FB6guFObgSLF60PpTrEHpaPfCZ9fwGjQoa7cnNI6mqoDU2eZQaFwDFCIiMLkeolnVl3gJOJjZwfVbbnqmPsk+0+keDt5wh9LZpKsLO80TipOeX5MaJrsUEpfR4WNAQoRUZhcL/FMnVSM4iI93F6/uq3XMcYMSrxZKBG7eJJokh3LDh6pMbiLx+70pv10ZcpfDFCIiMI4c3hQGxDYljtzRB9K6CTjFEs88TIoQ5E9KIkmiberM1CSOyQwXInRgMnlpojvQ8QAhYgojMxG5GqJBwBmBss8hztlgBIIqlIOUGL0oPj8Ar1hWRWfX8DljT/ltb03UJ5JdYuxJDMvn7DMQ0EMUIiIwuR6iQcAZlmDAYrMoKglnhR7UNRZKJEBSp/DDf+IhEmiRtnxlHgAqFmh90/3j+nrqfAwQCEiCpPru3iAsABlRAYl1bJUTYwDA2XAUllSpN6HRFuN5RbhVMbch1syvQoAsO9E75i+ngoPAxQiojC5vosHCO3k+aR7CE6PTy1LlY15F09kBkXOQKkpM6lBT7xhbX0Ot3o2UOOksQUoS2cEApSDp/uTHq1PhY0BChFRmHwo8UwuN6GypAh+ARw/O6gGDyVj7kEZkUEJBizVpUaUmmQGJXbQIMs7k8tNYz7DaOqkEkypLIbXL/BOe9+YvgcVFgYoRERh8qHEoyiK2ih7pHMADtkkm3IPSiBA6XW44QtrOukJy6CUqhmU2CUedQfPGPtPJJlF2XeiZ1zfhwoDAxQiojD5UOIBQmWej2wDoSbZFDMoVSWBAEUI4JwjVOaRPSjVZWEZlDhlFzmkbawNspIMUPayD4XAAIWIKMKwO7fnoEiyUfaIbSA0ByXFNRv0OkwqKQIQ2YfSrZZ4Qj0o0Q4UlMZ6ivFIMkB552QfXN7E02upsDFAISIKMxws8RTnegZFBiidg2rwkOqoeyBQxgEid/LIP1eXGVEWzMrEa1z9dJw7eKTzakpRU2aE2+vH+6e43VjrGKAQEYWRzaC53IMCADNrAwHK6b5hdNmdAFIv8QBhjbJDo0s8NWVG9T7EG3c/3hkokqIo6nZjlnmIAQoRUZDfL+D0BCamjnU3ykSpKCmC1RI4NfiTYAYj1RIPEDasLWoGxaQGPY4YTbJurx8d/cMAgGnjzKAA4Y2yDFC0jgEKEVGQM6zvIddLPAAwM1jmkUrHUuIplcPaovWghJpkY81BOd03DL8AzEU6TA4GO+MhA5S2T8+ppzWTNjFAISIKCp+Wmg8ByqzgeHhpXBmU4OGATo9PDUaqwwa1xZqDEl7eURQl5dcfabbVgnKzAYMuLz7sGBj396P8xQCFiChoWN1irINON/5ftpk2y2qJ+HwsTbIjDwyU/SdGvQ4WsyGsSTZ6iScUoKR+inE0el14HwrnoWgZAxQioqB82cEjyVkoUtlYmmRLI3tQwnfwKIqSsEm2PXj68HgbZMPJAGX/J+xD0TIGKEREQY48mYEiXVBbBllVURTAbBjLNuMRGZTB0JA2AAmbZEMZlOKUXzuW8EZZIUSCq6lQMUAhIgqSfRa5voNHKjbq1fHyJUX6MZWlRu7ikefyyMyKDFBiNcmeHQhcb60wp/zascyfUgFzkQ7nHB4c6xpM2/el/MIAhYgoyJlnJR4A6pk8Y5mBAoQyJUNuH4bdvogx90DofJ9YTbJ9wx4AQGVwbH46GA06fGbaJACch6JlDFCIiIJkiSdfMigAMNs6vgCl3GSAUR/4VdAz5Io4KBBAwlH3/Q4ZoBSN6fVj4TwUYoBCRBQU6kHJnwBF7uSxmMcWoCiKomZLegbdoR6U4HyUeKPuhRChDEpx+jIoALA02Cj7dvu5tH5fyh/50QlGRDQBBp2BX8Jj2Q2TLcvn1OLGpY1YMaduzN+jpsyEjn4neoZc6sh72Zsity473D74/SKiz2XQ5YXPH2hiTXcGZeqkQG9N+AA50pb8eRcSEWWY3RnIBliK0/vLNpPMRXpsvv7icX2P8Fko4duMgchgzeHxRXzeFyzvmIt0MKe5b6ciGPAMe3xweX0wjWGHEuU3lniIiIIGghmU8jGWS/JVaBZKqMRTE3zMZNBBJk0cI8o8MkBJd3kHCPTGyC3U/cEyEmkLAxQioiB78BehxZw/GZR0CM1Ccakj72UGRVEUdYT+yEbZvuFAMJPu8g4A6HSK+v+DnQGKJjFAISIKkhmUsTac5isZjHzSPQSPT0Q8BoR2CI1slJUZlIoMlcTk92UGRZsYoBARBckelHKNZVBkiecjW+BwvnKzIaLnQzbKjgpQhjOzxVhigKJtDFCIiILUDEqxNjMop/uGAYRmoEil6onGkSWefkewxJOBHhQgFPgwQNEmBihEREFazaCMDEjkDBSpNJhBGTnu/pxski3NzP2Su6lkKYm0hQEKEVFQqAdFWwFKeL9JtM9DGZSJ28UDsMSjdQxQiIgQmIoqd4tobZtxVenIACUyo1KiHhg4osSTwV08AAMUrWOAQkQEwOnxwxucippPg9rSwWTQRwRlNSMCljI5TTbmHBQGKJR+DFCIiAAMBPtPdEroBF8tCe9DGZVBiTkHJbjNOMMZFM5B0SYGKEREiGyQVRQlwdWFJ7wxdlQPSsw5KIESz6SSDO3iYQZF0xigEBEBsGt0zL0UkUEpHbnNODgHJaxJVggRKvFkOIPCXTzaxACFiAjaHXMvhWdNJpdHZkRkk6wjrEl2yO1Te3YytYvHwgyKpjFAISKCdg8KlKrjZFBkk2x4BkWWd4wGHcxFmflVwiZZbWOAQkSEUA+K1nbwSPLAQL1OGXW2jtok6woPUEI7eDLVsyObb11eP5weX4KrqdAwQCEiAjMoMmtSVWqEThcZcEQbdS+zGplqkAWAcpMB+uBauJNHexigEBGBPSgzakoj/jdctFH354IlnkxtMQYARVHUk6X7GKBojjb/qUBENEJozL02fyzObbBg663LYgQoozMomR7SJlUUF+Gcw8M+FA3S5juRiGgErfegAMBnz6+J+niJcXQGRQYMmdpiLKmNstxqrDks8RARgT0o8ZQFMyhurx8enx9AaBdPZQZ7UABuNdaylAOU119/Hddeey0aGhqgKApeeOGFiOdvvvlmKIoS8bFq1aqIa3p7e7F27VpYLBZUVlZi/fr1GBwcHNdfhIhoPEIHBWo3gxKL3MUDhMo8mR7SJnGrsXalHKAMDQ1hwYIF2LJlS8xrVq1ahY6ODvXjD3/4Q8Tza9euxQcffIAdO3Zg27ZteP3113HbbbelvnoiojQJ9aAwQBnJaNDBqA/8upBbjc+pPSiZzaDIAIhNstqTci5z9erVWL16ddxrTCYTrFZr1Oc+/PBDbN++Hfv378fixYsBAI899hiuvvpq/PSnP0VDQ0OqSyIiGrfQWTws8URTYtLD7fDDERzW1j8sSzwTk0GJtc24y+5EqcmgNvJS4chID8prr72G2tpazJo1C3fccQd6enrU51pbW1FZWakGJwCwYsUK6HQ67N27N+r3c7lcsNvtER9EROmkZlA03CQbT6k6rG1EiWcCdvEA0Us8vUNuXPGTV3Hjr/ZkdA2UHWkPUFatWoXf/e532LlzJ/7rv/4Lu3fvxurVq+HzBf6jttlsqK2tjfgag8GAqqoq2Gy2qN9z8+bNqKioUD8aGxvTvWwi0jCfX6g7VJhBiU7OQpElHllyyeQcFCB+gHKkcwBOjx8fnLHDFzwXiApH2t+JN9xwg/rn+fPn4+KLL8b555+P1157DcuXLx/T99y0aRM2btyofm632xmkEFHaDDpD22cZoESnjrt3+yCEULf9ZnKSLBA/QLH1OwEEAsyeIRdqy80ZXQtNrIxvMz7vvPNQU1ODY8eOAQCsViu6uroirvF6vejt7Y3Zt2IymWCxWCI+iIjSRfafmAw6mAz6LK8mN4VnUIY9PriD240z34MSCIDktuZwHcEABQC67K6MroMmXsYDlFOnTqGnpwf19fUAgObmZvT19aGtrU29ZteuXfD7/Vi2bFmml0NENAqHtCWm9qC4veoOHqNeh+KizAZ0oQyKd9Rztv5h9c9dA85Rz1N+SzmXOTg4qGZDAODEiRM4cOAAqqqqUFVVhQceeABr1qyB1WrF8ePHcdddd+GCCy7AypUrAQBz5szBqlWrcOutt+KJJ56Ax+PBhg0bcMMNN3AHDxFlBYe0JaaOu3f51GxGRUnmTjKWZI+LfdgDIUTE64VnUDqZQSk4KWdQ3nrrLSxcuBALFy4EAGzcuBELFy7EfffdB71ej/feew9f+tKXMHPmTKxfvx6LFi3CG2+8AZPJpH6Pp59+GrNnz8by5ctx9dVX4/LLL8cvf/nL9P2tiIhSoPWDApMRfmBg/wTt4AFCGRS3zw+nxx/xXGSAwgxKoUn5nwtXXnklhIjdLf3Xv/414feoqqrC1q1bU31pIqKMYAYlMVnicbi96g6eTPefBF5XD71Ogc8v0D/sQbExVFJiBqWw8SweItI89qAkFr6LJzTmPrM7eABAUZSoO3ncXj+6B0NBSRczKAWHAQoRaV5ozD0zKLGE7+I5Jw8KnKCATr5O+E6ekSWdrgFmUAoNAxQi0jz2oCQmm2SHXD41kzERJR4g+onGthEBCntQCg8DFCLSPPagJFZiDGVQZCZjIko8QPRhbbL/5LyaUgBA96ALXp9/9BdT3mKAQkSaxx6UxMpMYU2ywR6Uigm6X9ECFDkDZW6DBXqdAr8AeoZGD3Oj/MUAhYg0jxmUxCKaZIcnZsy9FO1EY5lBmTKpGJPLAmMsOE22sDBAISLNkxmUchMzKLGEN8mqc1AmqAdFvk5fRAYlEKDUW8yoswQCFPahFBb+c4GINE/dxcMST0yhJlmvenJwNks8Z2SAUlmMyeVmAP3o5Lj7gsIAhYg0b0BmUFjiiak0rMTj9E7MQYFS1F08wR6U+orwDApLPIWE70Yi0jz7MDMoicgSj88v1AxKtnbxeHx+de6JtcKMOosZAIe1FRr2oBCRpjk9PriD21OZQYlNNslKRXoFpcbMnmQsjQxQzg64IARg0CmoKTWpGRQOayssDFCISNNkg6yiAGVGBiix6HUKzEWhXxkVxcaMn2Qceq1ggBJszpU7eOosZuh0CmrLAxkUNskWFgYoRKRpskG2zGSATjcxv3DzlZyFAkxc/0n4a/UPeyCECO3gqQgEJrXsQSlIDFCISNM45j554WWeiTqHBwhlULx+AYfbh45gg6w1GKDIHpSeIRc8nCZbMBigEJGmcUhb8krCek4mMoNSXKRHkT6Q3eof9ozKoFSVGGHQKRACESccU35jgEJEmsYx98mLLPFMzA4eAFAUJaJRtiPYa2KtKAYA6HQKJpdzmmyhYYBCRJqmDmljBiWhElN2SjxA5CyUjr5AiachmEEBgFoLG2ULDQMUItI09qAkrzRLJR4g1IfS5wiVeKxhAUpdMIPSya3GBYMBChFpGntQklcalkGpmMASDxDK2JxzuNUgpD5Y4gHAYW0FiAEKEWkae1CSF5FBmeD7JTMox7sG4fML6MP6TgCglj0oBYcBChFpGjMoyQvPoEya4AyKDFAOdw4ACAQk+rC5NTKDwgMDCwcDFCLSNNmDUs4elIRKszSoDQgFKB/ZAgFKeP8JwGFthYgBChFpWmgXDwOURMLnoFRkaRfPWbX/JDJAYQ9K4WGAQkSaJntQWOJJLJsZlJFzV6yW4ojPZQ9Kz5Abbi+nyRYCBihEpGlqBoVNsgmVBkfd63VKxNC2iTAyY9NQGZlBmVRiVKfNcppsYWCAQkSaxgxK8kpMgRJPZXHRhJ1kLI0MUEb2oPBU48LDAIWINMvvFxh0sQclWbLvY+qk4gRXpt/IAGVkDwrARtlCw38yEJFmDbq9ECLwZ2ZQEpttteD/vmkRLqwtm/DXHp1BGR0kqbNQuNW4IPAdSUSaJbcYGw06mIv0Ca4mAFh5kTUrrxseoChKKBgJV8fzeAoKSzxEpFk8KDB/mIt0MBoCv7Iml5lQpB/96yu01ZglnkLAAIWINIsHBeYPRVHULEq0/hMglFXhgYGFgQEKEWkWx9znFxmgjNzBI9VyWFtBYYBCRJrFgwLzSyiDEn0XUZ26iycyQBl0eTm8LQ8xQCEizWIGJb9UJsig1AXnoJxzeODy+gAAL7/fgaUPv4Kv/rJ1YhZJacN3JRFpFntQ8st1C6fgdN8w/sfcuqjPV5YUwajXwe3zo8vuwnNvn8YjrxwBALzT3geX1weTgbu18gUDFCLSrAEXMyj55NoFDbh2QUPM5xVFweRyE073DeNfn34b75/uj3i+o8+J6TWlmV4mpQlLPESkWcygFB7Zh/L+6X4U6RX815r5uCA4WO5033A2l0YpYoBCRJrFHpTC01AZaKCtKjXi6f/zUnx1yTRMCT52+hwDlHzCdyURaVbooEBmUApFyxcuQJ3FjJs/Ox2NVSUAQkHLKWZQ8goDFCLSLLucJMttxgVjTr0F9/7T3IjH5OGGzKDkF5Z4iEizBtQMCv+tVsjUEk+fI8sroVQwQCEizRoMZlDKTAxQCtkUmUFhiSevMEAhIs0a9gSGeRUbORujkMkMSkefEz6/yPJqKFkMUIhIs1zB8efmIgYohazOYoZBp8DrF+ga4Dk9+YIBChFpkt8v1PNZzAb+KCxkep2ijsdno2z+4LuSiDTJFXZ4nIkZlIIXapRlgJIvGKAQkSY5g/0nADMoWsBG2fzDdyURaZIzeNqtQafAoOePwkI3ldNk8w7flUSkSU4PG2S1pIElnrzDAIWINMkVzKCYi/hjUAumcJps3uE7k4g0SWZQTAZmULQgvElWCM5CyQcpByivv/46rr32WjQ0NEBRFLzwwgsRzwshcN9996G+vh7FxcVYsWIFjh49GnFNb28v1q5dC4vFgsrKSqxfvx6Dg4Pj+osQEaVCNsmamEHRBFnicbh96HN4srwaSkbK78yhoSEsWLAAW7Zsifr8j3/8Yzz66KN44oknsHfvXpSWlmLlypVwOkPDcdauXYsPPvgAO3bswLZt2/D666/jtttuG/vfgogoRTJAMTODognmIj1qykwA2IeSL1I+gGL16tVYvXp11OeEEPj5z3+Oe+65B1/+8pcBAL/73e9QV1eHF154ATfccAM+/PBDbN++Hfv378fixYsBAI899hiuvvpq/PSnP0VDQ8M4/jpERMkJNckyg6IVUyYVo3vQhVPnhjFvSkW2l0MJpPWdeeLECdhsNqxYsUJ9rKKiAsuWLUNraysAoLW1FZWVlWpwAgArVqyATqfD3r17o35fl8sFu90e8UFENB6hJllmULRiKnfy5JW0Big2mw0AUFdXF/F4XV2d+pzNZkNtbW3E8waDAVVVVeo1I23evBkVFRXqR2NjYzqXTUQa5OI2Y83hTp78khe5zU2bNqG/v1/9OHnyZLaXRER5zsltxpoT2snjyPJKKBlpfWdarVYAQGdnZ8TjnZ2d6nNWqxVdXV0Rz3u9XvT29qrXjGQymWCxWCI+iIjGQ93FwyZZzeB5PPklrQHKjBkzYLVasXPnTvUxu92OvXv3orm5GQDQ3NyMvr4+tLW1qdfs2rULfr8fy5YtS+dyiIhiYpOs9jRw3H1eSXkXz+DgII4dO6Z+fuLECRw4cABVVVWYNm0a7rzzTvznf/4nLrzwQsyYMQP33nsvGhoacN111wEA5syZg1WrVuHWW2/FE088AY/Hgw0bNuCGG27gDh4imjDMoGiP7EE55/DA4faixJjyr0CaQCn/v/PWW2/hC1/4gvr5xo0bAQDr1q3DU089hbvuugtDQ0O47bbb0NfXh8svvxzbt2+H2WxWv+bpp5/Ghg0bsHz5cuh0OqxZswaPPvpoGv46RETJ4Vk82lNRXIRykwEDLi9OnxvGhXXl2V4SxZFygHLllVfGHROsKAoefPBBPPjggzGvqaqqwtatW1N9aSKitOFZPNo0ZVIxPrIN4FQfA5Rcx3cmEWkSz+LRpinsQ8kbDFCISJO4zVibZB/KGe7kyXl8ZxKRJrk8nCSrRdxqnD8YoBCRJnGbsTZxmmz+4DuTiDSJpxlrEzMo+YMBChFpksvLbcZaJDMonXYnPD5/lldD8TBAISJNCg1q449BLakpNcGo18EvAFu/M9vLoTj4ziQiTZK7eEzMoGiKTqegoTIwOPQU+1ByGgMUItIkNslql9ooyz6UnMZ3JhFpkpPbjDXrwtrABNnXj5zN8kooHgYoRKRJLp7Fo1lrPjMVALD9oA29Q+4sr4ZiYYBCRJrj9wu4gzs4zGyS1Zz5Uyswf0oF3D4//nfbyWwvh2LgO5OINEduMQbYJKtVX1s2DQDwh30n4x6AS9nDAIWINEf2nwDMoGjVlxY0oMxkwInuIbQe78n2cigKvjOJSHPkFmODToFBzx+DWlRqMuDLlzQAAJ7e157l1VA0fGcSkeY42SBLCJV5/vaBDd2DriyvhkZigEJEmhPaYswfgVp2UUMFFjRWwuMTePatU9leDo3AdycRaY5skjXxoEDNW7s0kEV5Zn87/H42y+YSBihEpDnqOTzMoGjePy2oR7nJgE97HPgHm2VzCt+dRKQ5aomHGRTNKzEa8M+fmQIA+G3rJ9xynEMYoBCR5vAcHgonm2V3HOrEvz/7bsQ2dMoevjuJSHNcXp7DQyGzrRY89OWLoNcpeO7t0/jar/aga8CZ7WVpHgMUItIcHhRII93UPB1P3bIEFrMBb7f34br/+XccPN2f7WVpGgMUItIcuYuHJR4K97kLJ+OFlstw3uRSnOl34itPtDJIySK+O4lIc9RdPGySpRHOm1yG5//1MiydXoVhjw9/3M/DBLOFAQoRaQ6bZCmeiuIi3HLZdADAvhO92V2MhvHdSUSawwwKJbJkRhUA4HDnAPoc7iyvRpsYoBCR5vAsHkqkpsyE8yeXAgD2f3Iuy6vRJgYoRKQ5Ti/P4qHEls6oBgDsO8EJs9nAdycRaY6LGRRKwrJgmYd9KNnBAIWINEdmUEwG/gik2JYGA5SDZ+wYdHlHPT/s9uEv73fAHdy2TunFdycRaY6Lg9ooCQ2VxZg6qRg+v8Dbn47uQ7n3xYO44+m3sXXvp1lYXeFjgEJEmsNtxpSspdOjl3nODbnx53fPAAAOddgnfF1awHcnEWkOTzOmZMkyz75PIgOU//ftU2ppp73XMeHr0gIGKESkOU4eFkhJkgHKgZN9amArhMDWfe3qNe09DFAygQEKEWmO3MVjYomHEphRU4qaMhPcXj/eOxU4l2fPx734+OwQjPrAfz8ddqd6QjalD9+dRKQ5oV08zKBQfIqihG03DsxDkdmTNYumosSohxDA6XPDWVtjoWKAQkSawyZZSoUs8+w90YueQRe2H+wAAKxdNg3TqkoAAJ+yDyXt+O4kIs1xcpsxpUAGKG2fnsMz+0/C4xNYMLUC86ZUoDEYoJxMMkDx+QVu+vVe3P6/2jK23kJhyPYCiIgmGifJUipm1ZXDYjbA7vTiF68eAwB8bdk0AFAzKMk2yp7oHsQbR7sBBAa9FRv532AszKAQkab4/AJuXzBA4SRZSoJOp2BJcB7KkNuHcpMB1y5oAAA0VadW4jlsG1T/POD0pHmlhYXvTiLSlPCx5MygULJkmQcArls4BSXGQAEi1RLP4c4B9c8DUcbnUwgDFCLSFNl/AvAsHkpeeIAiyztAWImn1wEhRMLvc9gWmjo74GSAEg97UIhIU+QWY4NOgUHPAIWSc/HUSnxl0VRUlRkxp96iPj51UjEUBXC4fegedGNyuSnu9znSyRJPshigEJGmONkgS2Og1yn4yVcWjHrcZNCj3mLGmX4n2nsdcQMUp8eHT3qG1M8HmUGJi/98ICJNCW0x5o8/So9k+1COdQ0ivArEEk98fIcSkabIAIVTZCld1J08CbYaH7YNRHxuZ4knLgYoRKQpLi+nyFJ6hTfKxhO+gwdgBiURvkOJSFOYQaF0a1QDlKG418kMSkVxEQBgkNuM42KAQkSawnN4KN2aqksBJM6gHAlmUBY1TQLAXTyJ8B1KRJri8vIcHkovWeLptLsi5uyE6x/2oKPfCSA8QGEGJZ60Byg/+MEPoChKxMfs2bPV551OJ1paWlBdXY2ysjKsWbMGnZ2d6V4GEVFUPCiQ0m1SSRHKTIGpHbF28sjsSUOFGQ2VZgAs8SSSkQzKRRddhI6ODvXjzTffVJ/79re/jZdeegnPPvssdu/ejTNnzuD666/PxDKIiEZhiYfSTVGUhI2ysv9kprUcZaZAD4qdGZS4MjKozWAwwGq1jnq8v78fv/71r7F161Z88YtfBAA8+eSTmDNnDvbs2YNLL700E8shIlKpJR42yVIaTasqwaEOe8wARWZQZtWVo9wc+NXLHpT4MvJPiKNHj6KhoQHnnXce1q5di/b2dgBAW1sbPB4PVqxYoV47e/ZsTJs2Da2trTG/n8vlgt1uj/ggIhoLmUExMYNCaZRoFoqaQYkIUJhBiSft79Bly5bhqaeewvbt2/H444/jxIkT+NznPoeBgQHYbDYYjUZUVlZGfE1dXR1sNlvM77l582ZUVFSoH42NjeleNhFpBLcZUybEmyYrhAhlUKzlKA+WeDjqPr60l3hWr16t/vniiy/GsmXL0NTUhD/96U8oLi4e0/fctGkTNm7cqH5ut9sZpBDRmPAsHsqEeD0oZwdcOOfwQKcAF9SWYdgdCJKHPT54fH4U8dDKqDJ+VyorKzFz5kwcO3YMVqsVbrcbfX19Edd0dnZG7VmRTCYTLBZLxAcR0Vg4vTyLh9JPlnjaex3w+0XEc3KC7PTqUpiL9Cgzh3IDzKLElvF36ODgII4fP476+nosWrQIRUVF2Llzp/r84cOH0d7ejubm5kwvhYiI24wpIxoqi6FTAkcpnB10RTwX3n8CAEV6HYqD//1xq3FsaS/x/Md//AeuvfZaNDU14cyZM7j//vuh1+tx4403oqKiAuvXr8fGjRtRVVUFi8WCb37zm2hubuYOHiKaEOpZPAZmUCh9ivQ6NFQW49S5YbT3OlBnMavPyf6TmdZy9bEyswHDHh8PDIwj7QHKqVOncOONN6KnpweTJ0/G5Zdfjj179mDy5MkAgEceeQQ6nQ5r1qyBy+XCypUr8Ytf/CLdyyAiisolm2SZQaE0a6ouwalzw/i0x4El06vUxw93DgIAZocFKOVmA84OuLiTJ460ByjPPPNM3OfNZjO2bNmCLVu2pPuliYgS4qA2ypRpVSX4O3oiGmX9foGjnZElHgAoNwd28jBAiY3vUCLSFLUHhduMKc2ibTU+dW4YDrcPRr0O04ONtABQHhyNP+hiiScWBihEpClOHhZIGdJUFTjV+NOeIfUxuYPn/NoyGMK2E3NYW2IZGXVPRJSrOEmWMkXOQnn/dD/+x892AwD6hgMZkll1ZRHXMkBJjAEKEWmKixkUypDzJpfCYjbA7vTiaNdgxHPLzquO+Jw9KIkxQCEiTVEzKNxmTGlWajJgx8bP4/jZyOCkzGTAvIaKUY8BPDAwHgYoRKQpHNRGmVRnMUfMQImFJZ7E+E8IItIUF8/ioRxgUUs8zKDEwgCFiDTD5xdw+zhJlrJPnsfDUfex8R1KRJohG2QBZlAou1jiSYwBChFphizvAGySpeziLp7E+A4lIs2QQ9oMOiViaBbRRAtlUNiDEgvfoUSkGU42yFKOCI2690IIkeXV5CYGKESkGaEtxvzRR9klSzx+AQy5fQmu1ia+S4lIM2SAYuJBgZRl5iIdDDoFAMs8sTBAISLNCJV4+KOPsktRlNBWYzbKRsV3KRFpBs/hoVwiG2XtDFCiYoBCRJrBc3gol5SbOE02Hr5LiUgzmEGhXFLOabJxMUAhIs3gQYGUSzhNNj4GKESkGWySpVxSzgMD4+K7lIg0Q82gcJsx5QBmUOJjgEJEmuHyBptkWeKhHFBmYoASDwMUItKM0KA2/uij7OOBgfHxXUpEmsGzeCiX8MDA+BigEJFmOL08i4dyB7cZx8d3KRFpBrcZUy5hk2x8DFCISDNcssTDHhTKAdxmHB/fpUSkGZwkS7mEGZT4GKAQkWaoZ/GwB4VygLrNmD0oUfFdSkSawUFtlEtkicft9avZPQphgEJEmuFkiYdyiMygACzzRMMAhYg0gyUeyiV6nYJSYyBYHmSAMgrfpUSkGdxmTLmG02RjY4BCRJohz+JhDwrlCk6TjY0BChFphnoWD0s8lCNkgGJnBmUUvkuJSDNcPIuHckxZsMTDcfejMUAhIk3w+QXcPk6SpdzCEk9sfJcSkSaEz5lgBoVyhYXTZGNigEJEmiC3GAMMUCh3yFkoLPGMxgCFiDRBZlCK9Ar0OiXLqyEK4IGBsTFAISJNUIe0cYsx5RDu4omNAQoRaUJoSBt/7FHu4KC22PhOJSJNUGegMINCOUTtQWGJZxQGKESkCU51Bgp/7FHu4C6e2PhOJSJN4EnGlItY4omNAQoRaQKnyFIuKjNzm3EsDFCISBPkNmMTp8hSDikPC1B8fpHl1eQWvlOJSBPkNk5mUCiXyAAFYBZlJAYoRFTwXF4ffvPmCQDAnPryLK+GKMRk0MMYzOpxWFskBihEVPD+nzdO4ET3ECaXm/B/ff78bC+HKEI5x91HxQCFiArayV4HHtt1FADw/avnwBLcNUGUK8q51TiqrAYoW7ZswfTp02E2m7Fs2TLs27cvm8shogL00LZDcHr8WDajCl++pCHbyyEahefxRJe1AOWPf/wjNm7ciPvvvx9vv/02FixYgJUrV6KrqytbSyKiAvPqR13426FOGHQKHrpuHhSFhwRS7pHTZJlBiZS1AOVnP/sZbr31Vtxyyy2YO3cunnjiCZSUlOA3v/lNtpZERAXE6fHh/j9/AAD4l8tnYGYdm2MpN7HEE50h8SXp53a70dbWhk2bNqmP6XQ6rFixAq2traOud7lccLlc6ud2uz0j62r7tBfb3uvIyPcmoon1SfcQ2nsdsFrM+LflF2Z7OUQxyRLPiwdO4/jZwSyvJmRR0yT808XZK4tmJUDp7u6Gz+dDXV1dxON1dXX46KOPRl2/efNmPPDAAxlf12HbIJ78+ycZfx0imjj3/NMcNYVOlIvqLCYAwP5PzmH/J+eyvJoQl9evvQAlVZs2bcLGjRvVz+12OxobG9P+Ohc1WNDyBW5BJCoUTdWluGZ+fbaXQRTX+stnoNRkgMOdWyWeBVMrs/r6WQlQampqoNfr0dnZGfF4Z2cnrFbrqOtNJhNMJlPG17WgsRILGisz/jpERERSdZkJLV+4INvLyDlZaZI1Go1YtGgRdu7cqT7m9/uxc+dONDc3Z2NJRERElEOyVuLZuHEj1q1bh8WLF2Pp0qX4+c9/jqGhIdxyyy3ZWhIRERHliKwFKF/96ldx9uxZ3HfffbDZbLjkkkuwffv2UY2zREREpD2KECLvzne22+2oqKhAf38/LBZLtpdDRERESUjl9zfP4iEiIqKcwwCFiIiIcg4DFCIiIso5DFCIiIgo5zBAISIiopzDAIWIiIhyDgMUIiIiyjkMUIiIiCjnMEAhIiKinJO1UffjIYff2u32LK+EiIiIkiV/byczxD4vA5SBgQEAQGNjY5ZXQkRERKkaGBhARUVF3Gvy8iwev9+PM2fOoLy8HIqiZPS17HY7GhsbcfLkSZ77k0G8zxOH93pi8D5PDN7niZGu+yyEwMDAABoaGqDTxe8yycsMik6nw9SpUyf0NS0WC//jnwC8zxOH93pi8D5PDN7niZGO+5wocyKxSZaIiIhyDgMUIiIiyjkMUBIwmUy4//77YTKZsr2Ugsb7PHF4rycG7/PE4H2eGNm4z3nZJEtERESFjRkUIiIiyjkMUIiIiCjnMEAhIiKinMMAhYiIiHIOA5QEtmzZgunTp8NsNmPZsmXYt29ftpeU1zZv3owlS5agvLwctbW1uO6663D48OGIa5xOJ1paWlBdXY2ysjKsWbMGnZ2dWVpx/vvRj34ERVFw5513qo/xHqfP6dOn8fWvfx3V1dUoLi7G/Pnz8dZbb6nPCyFw3333ob6+HsXFxVixYgWOHj2axRXnH5/Ph3vvvRczZsxAcXExzj//fDz00EMR57nwPqfu9ddfx7XXXouGhgYoioIXXngh4vlk7mlvby/Wrl0Li8WCyspKrF+/HoODg+lZoKCYnnnmGWE0GsVvfvMb8cEHH4hbb71VVFZWis7OzmwvLW+tXLlSPPnkk+LgwYPiwIED4uqrrxbTpk0Tg4OD6jW33367aGxsFDt37hRvvfWWuPTSS8VnP/vZLK46f+3bt09Mnz5dXHzxxeJb3/qW+jjvcXr09vaKpqYmcfPNN4u9e/eKjz/+WPz1r38Vx44dU6/50Y9+JCoqKsQLL7wg3n33XfGlL31JzJgxQwwPD2dx5fnl4YcfFtXV1WLbtm3ixIkT4tlnnxVlZWXiv//7v9VreJ9T9/LLL4vvf//74rnnnhMAxPPPPx/xfDL3dNWqVWLBggViz5494o033hAXXHCBuPHGG9OyPgYocSxdulS0tLSon/t8PtHQ0CA2b96cxVUVlq6uLgFA7N69WwghRF9fnygqKhLPPvuses2HH34oAIjW1tZsLTMvDQwMiAsvvFDs2LFDfP7zn1cDFN7j9Ln77rvF5ZdfHvN5v98vrFar+MlPfqI+1tfXJ0wmk/jDH/4wEUssCNdcc434l3/5l4jHrr/+erF27VohBO9zOowMUJK5p4cOHRIAxP79+9Vr/vKXvwhFUcTp06fHvSaWeGJwu91oa2vDihUr1Md0Oh1WrFiB1tbWLK6ssPT39wMAqqqqAABtbW3weDwR93327NmYNm0a73uKWlpacM0110TcS4D3OJ3+/Oc/Y/HixfjKV76C2tpaLFy4EL/61a/U50+cOAGbzRZxrysqKrBs2TLe6xR89rOfxc6dO3HkyBEAwLvvvos333wTq1evBsD7nAnJ3NPW1lZUVlZi8eLF6jUrVqyATqfD3r17x72GvDwscCJ0d3fD5/Ohrq4u4vG6ujp89NFHWVpVYfH7/bjzzjtx2WWXYd68eQAAm80Go9GIysrKiGvr6upgs9mysMr89Mwzz+Dtt9/G/v37Rz3He5w+H3/8MR5//HFs3LgR3/ve97B//37827/9G4xGI9atW6fez2g/R3ivk/fd734Xdrsds2fPhl6vh8/nw8MPP4y1a9cCAO9zBiRzT202G2prayOeNxgMqKqqSst9Z4BCWdPS0oKDBw/izTffzPZSCsrJkyfxrW99Czt27IDZbM72cgqa3+/H4sWL8cMf/hAAsHDhQhw8eBBPPPEE1q1bl+XVFY4//elPePrpp7F161ZcdNFFOHDgAO688040NDTwPhcwlnhiqKmpgV6vH7WzobOzE1arNUurKhwbNmzAtm3b8Oqrr2Lq1Knq41arFW63G319fRHX874nr62tDV1dXfjMZz4Dg8EAg8GA3bt349FHH4XBYEBdXR3vcZrU19dj7ty5EY/NmTMH7e3tAKDeT/4cGZ/vfOc7+O53v4sbbrgB8+fPx0033YRvf/vb2Lx5MwDe50xI5p5arVZ0dXVFPO/1etHb25uW+84AJQaj0YhFixZh586d6mN+vx87d+5Ec3NzFleW34QQ2LBhA55//nns2rULM2bMiHh+0aJFKCoqirjvhw8fRnt7O+97kpYvX473338fBw4cUD8WL16MtWvXqn/mPU6Pyy67bNQ2+SNHjqCpqQkAMGPGDFit1oh7bbfbsXfvXt7rFDgcDuh0kb+u9Ho9/H4/AN7nTEjmnjY3N6Ovrw9tbW3qNbt27YLf78eyZcvGv4hxt9kWsGeeeUaYTCbx1FNPiUOHDonbbrtNVFZWCpvNlu2l5a077rhDVFRUiNdee010dHSoHw6HQ73m9ttvF9OmTRO7du0Sb731lmhubhbNzc1ZXHX+C9/FIwTvcbrs27dPGAwG8fDDD4ujR4+Kp59+WpSUlIjf//736jU/+tGPRGVlpXjxxRfFe++9J7785S9z+2uK1q1bJ6ZMmaJuM37uuedETU2NuOuuu9RreJ9TNzAwIN555x3xzjvvCADiZz/7mXjnnXfEp59+KoRI7p6uWrVKLFy4UOzdu1e8+eab4sILL+Q244ny2GOPiWnTpgmj0SiWLl0q9uzZk+0l5TUAUT+efPJJ9Zrh4WHxr//6r2LSpEmipKRE/PM//7Po6OjI3qILwMgAhfc4fV566SUxb948YTKZxOzZs8Uvf/nLiOf9fr+49957RV1dnTCZTGL58uXi8OHDWVptfrLb7eJb3/qWmDZtmjCbzeK8884T3//+94XL5VKv4X1O3auvvhr15/G6deuEEMnd056eHnHjjTeKsrIyYbFYxC233CIGBgbSsj5FiLBRfEREREQ5gD0oRERElHMYoBAREVHOYYBCREREOYcBChEREeUcBihERESUcxigEBERUc5hgEJEREQ5hwEKERER5RwGKERERJRzGKAQERFRzmGAQkRERDmHAQoRERHlnP8f+0u+8RwaBFkAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.Series(subject_infos['age']).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Axes: >"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGdCAYAAAA44ojeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABAaElEQVR4nO3de3iU9Z338c/MJJmcMySQhEMSwBOgoAgK8dB6oKCytlaePtWlSlu3XrrRVXlWrVu1VtfiY7fa2kXddi3YR61buh4q6wlRsUpACKAIFlGQhMMkQjKZHMhp5n7+SO47CQRIyGTue2ber+ua6yIzdzK/uS8kH3+/7+/7cxmGYQgAAMBB3HYPAAAA4FAEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DhJdg/geITDYe3du1dZWVlyuVx2DwcAAPSDYRhqaGjQqFGj5HYffY4kJgPK3r17VVRUZPcwAADAcaiqqtKYMWOOek1MBpSsrCxJnR8wOzvb5tEAAID+CAaDKioqsn6PH01MBhRzWSc7O5uAAgBAjOlPeQZFsgAAwHEIKAAAwHEIKAAAwHEIKAAAwHEIKAAAwHEIKAAAwHEGFFDuu+8+uVyuXo8JEyZYr7e0tKisrEx5eXnKzMzUvHnzVF1d3etnVFZWau7cuUpPT1d+fr5uv/12dXR0RObTAACAuDDgPiinnnqq3nrrre4fkNT9I2677Tb9z//8j5YtW6acnBzddNNNuvLKK/XBBx9IkkKhkObOnavCwkKtXr1a+/bt07XXXqvk5GT9/Oc/j8DHAQAA8WDAASUpKUmFhYWHPV9fX6+nnnpKzz33nC666CJJ0pIlSzRx4kStWbNGM2fO1JtvvqmtW7fqrbfeUkFBgc444ww98MADuvPOO3XfffcpJSVl8J8IAADEvAHXoGzfvl2jRo3S+PHjNX/+fFVWVkqSKioq1N7erlmzZlnXTpgwQcXFxSovL5cklZeXa/LkySooKLCumTNnjoLBoLZs2XLE92xtbVUwGOz1AAAA8WtAAWXGjBlaunSpXn/9dT3xxBPauXOnzj//fDU0NMjv9yslJUU+n6/X9xQUFMjv90uS/H5/r3Bivm6+diSLFi1STk6O9eCgQAAA4tuAlnguvfRS689TpkzRjBkzVFJSoj/96U9KS0uL+OBMd911lxYuXGh9bR42BAAA4tOgthn7fD6dfPLJ+vzzz1VYWKi2tjYFAoFe11RXV1s1K4WFhYft6jG/7quuxeT1eq2DATkgEEPhlY/2as2OA3YPAwDQZVABpbGxUV988YVGjhypadOmKTk5WStXrrRe37ZtmyorK1VaWipJKi0t1ebNm1VTU2Nds2LFCmVnZ2vSpEmDGQpw3LZXN+jmP27Ujc9UyDAMu4cDANAAl3j++Z//WZdffrlKSkq0d+9e/fSnP5XH49HVV1+tnJwcXXfddVq4cKFyc3OVnZ2tm2++WaWlpZo5c6Ykafbs2Zo0aZKuueYaPfzww/L7/br77rtVVlYmr9c7JB8QOJZ1X9ZJkuqa27W/sU0jsvi7CAB2G1BA2b17t66++modOHBAI0aM0Hnnnac1a9ZoxIgRkqRHH31Ubrdb8+bNU2trq+bMmaPHH3/c+n6Px6Ply5frxhtvVGlpqTIyMrRgwQLdf//9kf1UwABsqKyz/lxZ20RAAQAHcBkxOKcdDAaVk5Oj+vp66lEwaBf98l3t+KpJkvTL75yuedPG2DwiAIhPA/n9zVk8SGiB5jYrnEjSrtpmG0cDADARUJDQNlYFen2960BT3xcCAKKKgIKEtnFXZ/1JTlqyJGnXAWZQAMAJCChIaOYMytwpIyUxgwIATkFAQcIKhw1tqgxIkr49dbSkzq3GwZZ2G0cFAJAIKEhg22sa1dDaofQUj6YW+TQ8s/M07UqWeQDAdgQUJKyNXf1PpozJUZLHrZK8DEnSlyzzAIDtCChIWGaDtjOLh0mSSnLTJVEoCwBOQEBBwtrQVX9iBZSuGRQKZQHAfgQUJKT6g+36vKZRknRGsU+SVJLHDAoAOAUBBQlpU9f24pK8dA3P7Dx7p5iAAgCOQUBBQtqwq3f9iSSN7Vri8Qdb1NIesmVcAIBOBBQkJLNAdmrX8o4kDUtPVpa384DvKs7kAQBbEVCQcMJhw1ri6TmD4nK5VDK8c5nnS5Z5AMBWBBQknC++alRDS4dSk92aUJjV67WSXHbyAIATEFCQcDZYDdp8SvL0/k+AnTwA4AwEFCScjYf0P+nJCijUoACArQgoSDh9FciailniAQBHIKAgoQRb2rW9q0FbXzMoY7uKZPfUHVRHKBzVsQEAuhFQkFA+qgrIMKSi3DSNyPIe9npBVqpSktzqCBvaG2ixYYQAAImAggSzYVdAUt+zJ5LkdrusQwM51RgA7ENAQUI59ATjvlAoCwD2I6AgYfRs0NZXgazJOtV4PzMoAGAXAgoSxo79Tao/2K7UZLcmjsw+4nXMoACA/QgoSBhWg7bRPiV7jvxXvzjXbNbGDAoA2IWAgoSx8Sj9T3oyTzWurG2WYRhDPSwAQB8IKEgYZgfZqUcpkJWk0cPS5HG71NIeVk1DaxRGBgA4FAEFCaGhpV3bqhskSWeW+I56bbLHrdG+NEnSlxTKAoAtCChICB9V1cswpDHD0pSflXrM6zk0EADsRUBBQuiuPzn68o5pzLDOGZS99QeHbEwAgCMjoCAhdDdo8/Xr+hFdsyxfUYMCALYgoCDuGYahjV0N2o7WQbYn85weimQBwB4EFMS9nfubFGhulzfp6A3aesonoACArQgoiHsburYXTx6do5Sk/v2VNwPKfgIKANiCgIK4Z9WflPRveUfqXuL5qqGVZm0AYAMCCuLehl0DK5CVugNKWyis+oPtQzEsAMBREFAQ1xpbO/RZV4O2/m4xliRvkke+9GRJ1KEAgB0IKIhrH3y+X2FDKspNU0H2sRu09TQis6tQNkhAAYBoI6Agrr2xxS9JmjWxYMDfm5/dVYfS2BLRMQEAjo2AgrjVHgpr5ac1kqQ5pxYO+PvNlvjMoABA9BFQELc+3Fmr+oPtys1I0Vljcwf8/TRrAwD7EFAQt7qXd/LlcbsG/P35PbYaAwCii4CCuBQOG3pzS7Wk41vekXrOoFCDAgDRRkBBXPp4T738wRZlpHh07onDj+tnsMQDAPYhoCAumcs7F5ySr9Rkz3H9jHxONAYA2xBQEJfMgDL71IFvLzaZMygNLR1qaQ9FZFwAgP4hoCDufF7ToB1fNSnZ49KFE/KP++dkpybJ23W4ILMoABBdBBTEnTe6imPPOWG4slOTj/vnuFwuq1kbhbIAEF0EFMQdc3nneHfv9ES7ewCwBwEFcWVv4KA+3l0vl0v6xqTjrz8xWYWyjQQUAIgmAgriyptdsyfTiodZRa6DYS3xMIMCAFFFQEFcWf3FAUnSrAjMnkg9lnioQQGAqCKgIK5U1R2UJJ1SkBWRn2edaMwuHgCIKgIK4sqeumZJ0uhhaRH5edaJxgQUAIgqAgriRrClXcGWDknSaF9kAgrt7gHAHgQUxI09Xcs7w9KTleFNisjPNE80PtDYqlDYiMjPBAAcGwEFcWN3V0AZMyw9Yj8zNyNFLpcUNqQDTcyiAEC0DCqgPPTQQ3K5XLr11lut51paWlRWVqa8vDxlZmZq3rx5qq6u7vV9lZWVmjt3rtLT05Wfn6/bb79dHR0dgxkK0F1/EqHlHUlK8riVl0GhLABE23EHlHXr1uk//uM/NGXKlF7P33bbbXrllVe0bNkyrVq1Snv37tWVV15pvR4KhTR37ly1tbVp9erVevrpp7V06VLde++9x/8pAPWcQYlcQJG6l3moQwGA6DmugNLY2Kj58+frd7/7nYYNG2Y9X19fr6eeekqPPPKILrroIk2bNk1LlizR6tWrtWbNGknSm2++qa1bt+qZZ57RGWecoUsvvVQPPPCAFi9erLa2tsh8KiSkPYHOgBKpHTwms1D2K5q1AUDUHFdAKSsr09y5czVr1qxez1dUVKi9vb3X8xMmTFBxcbHKy8slSeXl5Zo8ebIKCrobac2ZM0fBYFBbtmzp8/1aW1sVDAZ7PYBDDUUNitQ9g0K7ewCIngFvdXj++ee1YcMGrVu37rDX/H6/UlJS5PP5ej1fUFAgv99vXdMznJivm6/1ZdGiRfrZz3420KEiwVgzKBGsQZF6trunmywARMuAZlCqqqp0yy236Nlnn1VqaupQjekwd911l+rr661HVVVV1N4bsaG5rUO1TZ1LhBFf4smkBgUAom1AAaWiokI1NTU688wzlZSUpKSkJK1atUqPPfaYkpKSVFBQoLa2NgUCgV7fV11drcLCQklSYWHhYbt6zK/Naw7l9XqVnZ3d6wH0ZPZAyUpNUk5ackR/dn5214nGBBQAiJoBBZSLL75Ymzdv1qZNm6zH9OnTNX/+fOvPycnJWrlypfU927ZtU2VlpUpLSyVJpaWl2rx5s2pqaqxrVqxYoezsbE2aNClCHwuJZqjqTyR28QCAHQZUg5KVlaXTTjut13MZGRnKy8uznr/uuuu0cOFC5ebmKjs7WzfffLNKS0s1c+ZMSdLs2bM1adIkXXPNNXr44Yfl9/t19913q6ysTF6vN0IfC4lm9xDVn0g92923yDAMuVyuiL8HAKC3yPQD7+HRRx+V2+3WvHnz1Nraqjlz5ujxxx+3Xvd4PFq+fLluvPFGlZaWKiMjQwsWLND9998f6aEggezuatIW6R4oUveBgS3tYTW2digrNbJLSACAww06oLz77ru9vk5NTdXixYu1ePHiI35PSUmJXn311cG+NWDZM0RN2iQpLcWjLG+SGlo7VNPQSkABgCjgLB7EhaHqImuymrVRhwIAUUFAQVzo7oES+SJZqWcdCgEFAKKBgIKY19IesmY2hnoGhWZtABAdBBTEvL1dsyfpKR750oemPsQslKXdPQBEBwEFMa9n/clQbQE2291zYCAARAcBBTFvqM7g6Yl29wAQXQQUxLw9Q9hF1mTNoBBQACAqCCiIeWaTtkgfEtiTWYNS00CRLABEAwEFMc9c4hmqHTySVJjTGVDqmtt1sC00ZO8DAOhEQEHMM4tkh7IGJSctWVnezsbLewLNQ/Y+AIBOBBTEtLaOsKq7epMMZQ2K1L2EZAYiAMDQIaAgpvnrWxQ2JG+SW8MzU4b0vcYQUAAgaggoiGm7u5ZbRvuGrgeKyZyhMWteAABDh4CCmGbVnwxhgazJrHFhBgUAhh4BBTFtzxCfYtyT+R576iiSBYChRkBBTNsdhSZtJopkASB6CCiIaXt61KAMNTME1TS0qqWdXigAMJQIKIhpu6O4xDMsPVlpyR5J0r56OsoCwFAioCBmdYTC8ncFhWgUybpcrh5bjalDAYChREBBzKpuaFVH2FCyx2WdlTPURluFstShAMBQIqAgZv1tX1CSVJybLo97aHugmGjWBgDRQUBBzNpYGZAkTS0eFrX3HO2jWRsARAMBBTFrQ2WdJOnMKAYUalAAIDoIKIhJobChj6oCkqSpxb6ovS81KAAQHQQUxKTPqhvU1BZSpjdJJxdkRe19zRkUf7BF7aFw1N4XABINAQUxyVzeOb0oJ2oFspI0PMOrlCS3woasLc4AgMgjoCAmbdgVkCRNLYpe/Ykkud0ujenqWltFHQoADBkCCmLSxqquAtkSX9TfmzoUABh6BBTEnEBzm3Z81SQp+jMoEr1QACAaCCiIOWb/k3HDMzQsIyXq728eTEgvFAAYOgQUxJyNXQWy0dxe3JN5qjG9UABg6BBQEHM2dM2gRLNBW09WDQozKAAwZAgoiCmhsKFNNjRo68msQdkXaFEobNgyBgCIdwQUxJTPaxrV2Nqh9BSPTolig7ae8rNSleR2qSNsqDpILxQAGAoEFMQUq0HbGJ+SPPb89fW4XRrlYycPAAwlAgpiyoZd9hbImrp38lAoCwBDgYCCmLKxq/7ErgJZk9ULpZYZFAAYCgQUxIz65nZ9XtMoyQEzKOzkAYAhRUBBzNi0OyBJKslLV16m19axdPdCIaAAwFAgoCBmmPUndi/vSHSTBYChRkBBzNi6LyhJOn1Mjs0j6a5B2VN3UGF6oQBAxBFQEDMONLZKkgpzUm0eSecY3C6pLRTW/q5xAQAih4CCmBFobpck5aRF/4DAQyV73BqZ0zmLUkUdCgBEHAEFMSNwsDOgDMtItnkknahDAYChQ0BBTAiHDQWa2yRJw9Ltn0GRpDG5nQFl696gzSMBgPhDQEFMaGjpkFmL6kt3xgzKhafkS5Je+WgvhbIAEGEEFMSEuq7Zk/QUj7xJHptH0+kbkwqUlZqkPYGDWrPzgN3DAYC4QkBBTKhz2PKOJKUme/R3U0ZKkv67Yo/NowGA+EJAQUwwd/A4ZXnHNO/MMZKk1z7Zp+a2DptHAwDxg4CCmODEGRRJmlYyTCV56WpuC+mNLX67hwMAcYOAgphQ59AZFJfLpSunds6isMwDAJFDQEFMMLcYOy2gSNKVZ46WJH3wxX7tq6cnCgBEAgEFMcGsQXHaEo8kFeWm6+xxuTIM6cWNzKIAQCQQUBAT6qwZFOcFFEma1zWL8sKGPTIMeqIAwGARUBATumdQnLfEI0mXTR4pb5Jbn9c06uPd9XYPBwBiHgEFMcGpu3hMWanJmnNqoSTphQ27bR4NAMQ+AgpiglP7oPQ0b1rnbp5XPt7HMg8ADNKAAsoTTzyhKVOmKDs7W9nZ2SotLdVrr71mvd7S0qKysjLl5eUpMzNT8+bNU3V1da+fUVlZqblz5yo9PV35+fm6/fbb1dFBgyscndNnUCTpnBPylOJxq7apTbvr2M0DAIMxoIAyZswYPfTQQ6qoqND69et10UUX6Vvf+pa2bNkiSbrtttv0yiuvaNmyZVq1apX27t2rK6+80vr+UCikuXPnqq2tTatXr9bTTz+tpUuX6t57743sp0Jcae0IqbktJMnZASXZ49aJ+ZmSpK37OOEYAAZjQAHl8ssv12WXXaaTTjpJJ598sh588EFlZmZqzZo1qq+v11NPPaVHHnlEF110kaZNm6YlS5Zo9erVWrNmjSTpzTff1NatW/XMM8/ojDPO0KWXXqoHHnhAixcvVltb25B8QMQ+c3nH7ZKyUpNsHs3RTRyZLUn6lIACAINy3DUooVBIzz//vJqamlRaWqqKigq1t7dr1qxZ1jUTJkxQcXGxysvLJUnl5eWaPHmyCgoKrGvmzJmjYDBozcL0pbW1VcFgsNcDiaPnFmO322XzaI5u4sgsSQQUABisAQeUzZs3KzMzU16vVzfccINefPFFTZo0SX6/XykpKfL5fL2uLygokN/feUaJ3+/vFU7M183XjmTRokXKycmxHkVFRQMdNmKYVSCb5twCWdMkawalweaRAEBsG3BAOeWUU7Rp0yatXbtWN954oxYsWKCtW7cOxdgsd911l+rr661HVVXVkL4fnMXJbe4PZS7xVNY2q6Gl3ebRAEDsGvCCfkpKik488URJ0rRp07Ru3Tr9+te/1ne/+121tbUpEAj0mkWprq5WYWFnf4jCwkJ9+OGHvX6eucvHvKYvXq9XXq93oENFnKhzcJv7Qw3LSNHInFTtq2/R3/wNOmtsrt1DAoCYNOg+KOFwWK2trZo2bZqSk5O1cuVK67Vt27apsrJSpaWlkqTS0lJt3rxZNTU11jUrVqxQdna2Jk2aNNihIE45vc39oSiUBYDBG9AMyl133aVLL71UxcXFamho0HPPPad3331Xb7zxhnJycnTddddp4cKFys3NVXZ2tm6++WaVlpZq5syZkqTZs2dr0qRJuuaaa/Twww/L7/fr7rvvVllZGTMkOCKnt7k/1MSRWXr7bzUEFAAYhAEFlJqaGl177bXat2+fcnJyNGXKFL3xxhv6xje+IUl69NFH5Xa7NW/ePLW2tmrOnDl6/PHHre/3eDxavny5brzxRpWWliojI0MLFizQ/fffH9lPhbhS19TVpC0jtmZQtu4loADA8RpQQHnqqaeO+npqaqoWL16sxYsXH/GakpISvfrqqwN5WyS4uhhoc9+TuZNnW3WDQmFDHodvjQYAJ+IsHjheIAba3PdUkpehtGSPWtrD2rm/ye7hAEBMIqDA8epiaJuxJHncLp1SSMM2ABgMAgocr/6g2agtNmZQpB51KAQUADguBBQ4mmEY3bt4MmJjBkWSJtHyHgAGhYACR2to7VBH2JAUOzUokjRpFL1QAGAwCChwtEBT5+xJarJbqckem0fTf6cUdgaU6mCrDjS22jwaAIg9BBQ4Wl2M7eAxZXqTVJKXLomDAwHgeBBQ4Gix1ua+p4mFLPMAwPEioMDRYq3NfU/UoQDA8SOgwNFidYlHYqsxAAwGAQWOFmtt7nua2LXV+IuvGtXWEbZ5NAAQWwgocLRYa3Pf02hfmrJTk9QeMrS9hkJZABgIAgocLRDDMygul0sTRpp1KAQUABgIAgocLZZ38UjSxK4zeT6rJqAAwEAQUOBosbyLR5JOzM+UJH1R02jzSAAgthBQ4GixPoNyQldA+fwrAgoADAQBBY4WLzMoVbXNamkP2TwaAIgdBBQ4VltHWI2tHZJicxePJI3I9CorNUlhQ/ryQJPdwwGAmEFAgWMFDnYu77hcUnZabM6guFyuHnUoBBQA6C8CChzLXN7JSUuWx+2yeTTH74QRXXUoFMoCQL8RUOBYdU2x26StJ2sGhUJZAOg3AgocK3Awdpu09cQMCgAMHAEFjmW2uffFaP2JyZxB2bG/UeGwYfNoACA2EFDgWHXWFuPYXuIpGpamFI9bLe1h7QkctHs4ABATCChwrFhv0mZK8rg1dni6JOpQAKC/CChwrEBTbDdp64k6FAAYGAIKHMuaQcmI7RkUqedOHnqhAEB/EFDgWLHe5r4ncwaFQwMBoH8IKHAscwYl1otkJXqhAMBAEVDgWOYunljvgyJJ40dkSJIONLVZDegAAEdGQIEjGYah+oPxM4OSnpKk0b40ScyiAEB/EFDgSE1tIbWHOpuaxcMMitQ9i8JOHgA4NgIKHMlcBklJcist2WPzaCKDOhQA6D8CChyp5w4elyt2TzLuiV4oANB/BBQ4Ujzt4DHRCwUA+o+AAkeKx4BizqBU1TWrpT1k82gAwNkIKHCk2q4alNw46CJrGp6Zopy0ZBmGtHM/sygAcDQEFDiSdZJxRnzs4JEkl8ulE7p28lAoCwBHR0CBI5m7eOJpiUfqrkOhUBYAjo6AAkeKxxoUqceZPBTKAsBREVDgSFZAiaMlHokZFADoLwIKHKmuyeyDEp8zKDu+apRhGDaPBgCci4ACR4rXJZ6RvlRJUmtH2GpGBwA4HAEFjmQGlHjaZixJ3iSP9ZmqG1psHg0AOBcBBY5zsC2klvawJGlYnAUUScrP8kqSqoOtNo8EAJyLgALHMWdPkj0uZaTEx0GBPRVkdy7zVAeZQQGAIyGgwHFqe/RAiZeDAnsyZ1C+amAGBQCOhIACx+k+yTj+lnckZlAAoD8IKHCc2jjtgWIqyDZrUAgoAHAkBBQ4Try2uTflWzMoLPEAwJEQUOA43V1k4zSgUIMCAMdEQIHjdM+gxOsST+cMSk1Di8JhuskCQF8IKHCcujgvkh3RNYPSHjKs2SIAQG8EFDhOvHaRNSV73Bqe2dVNljoUAOgTAQWOE6/n8PSUn9W9zAMAOBwBBY5jnWQcpzMokpTftdW4hhkUAOgTAQWO0z2DEp9FspJUkEWzNgA4mgEFlEWLFumss85SVlaW8vPzdcUVV2jbtm29rmlpaVFZWZny8vKUmZmpefPmqbq6utc1lZWVmjt3rtLT05Wfn6/bb79dHR0dg/80iHkt7SE1t4UkxfcMitWsjSUeAOjTgALKqlWrVFZWpjVr1mjFihVqb2/X7Nmz1dTUZF1z22236ZVXXtGyZcu0atUq7d27V1deeaX1eigU0ty5c9XW1qbVq1fr6aef1tKlS3XvvfdG7lMhZplt7pPcLmV5k2wezdAxm7WxxAMAfRvQb4DXX3+919dLly5Vfn6+Kioq9LWvfU319fV66qmn9Nxzz+miiy6SJC1ZskQTJ07UmjVrNHPmTL355pvaunWr3nrrLRUUFOiMM87QAw88oDvvvFP33XefUlLi9/+acWzmQYG+OD0o0GQ2a6umWRsA9GlQNSj19fWSpNzcXElSRUWF2tvbNWvWLOuaCRMmqLi4WOXl5ZKk8vJyTZ48WQUFBdY1c+bMUTAY1JYtW/p8n9bWVgWDwV4PxKdEqD+RejRrowYFAPp03AElHA7r1ltv1bnnnqvTTjtNkuT3+5WSkiKfz9fr2oKCAvn9fuuanuHEfN18rS+LFi1STk6O9SgqKjreYcPh4r3Nvam7m2wr3WQBoA/HHVDKysr0ySef6Pnnn4/kePp01113qb6+3npUVVUN+XvCHmab+9w47oEiScMzU+RySaGwoQNNdJMFgEMdV0C56aabtHz5cr3zzjsaM2aM9XxhYaHa2toUCAR6XV9dXa3CwkLrmkN39Zhfm9ccyuv1Kjs7u9cD8clqc58R30s8SR63hmd29UJhJw8AHGZAAcUwDN1000168cUX9fbbb2vcuHG9Xp82bZqSk5O1cuVK67lt27apsrJSpaWlkqTS0lJt3rxZNTU11jUrVqxQdna2Jk2aNJjPgjhQ2xT/XWRNZqEsO3kA4HAD2sVTVlam5557Ti+//LKysrKsmpGcnBylpaUpJydH1113nRYuXKjc3FxlZ2fr5ptvVmlpqWbOnClJmj17tiZNmqRrrrlGDz/8sPx+v+6++26VlZXJ6/VG/hMipgQSoM29qSA7VVv2BmnWBgB9GFBAeeKJJyRJF1xwQa/nlyxZou9///uSpEcffVRut1vz5s1Ta2ur5syZo8cff9y61uPxaPny5brxxhtVWlqqjIwMLViwQPfff//gPgniQm1z/Le5N1nN2phBAYDDDCigGMaxdxukpqZq8eLFWrx48RGvKSkp0auvvjqQt0aCCCTINmOJAwMB4Gg4iweOYtWgJMAMSj4zKABwRAQUOIrZ6j4halCYQQGAIyKgwDHaOsJqbO08NDLe+6BI3c3aKJIFgMMRUOAYZv2Jx+1SVmr8HhRoMotk9ze2KUQ3WQDohYACx6jtCii+tGS53fF7UKApL9Mrt9VNljoUAOiJgALHSKQCWalzpsjqJkuhLAD0QkCBY3QXyMb/FmMTdSgA0DcCChwjkdrcm8w6lJoGZlAAoCcCChwjkdrcm/KZQQGAPhFQ4Bi1TYnT5t5kHhhIszYA6I2AAsdIpDb3JrMGpYYZFADohYACxzC3GSfSDIp1YCDdZAGgFwIKHKOuaxdPInSRNVkHBrLEAwC9EFDgGHVWH5TEWeLJt7rJtqojFLZ5NADgHAQUOEZdAu7iycvwyuN2KWxIB7oCGgCAgAKHaA+F1dDSeVBgIgUUj9ulEZnmTh7qUADARECBI5izJ26XlJ2WOEs8Uo9mbdShAICFgAJHMNvc56Qly5MABwX2ZDVrYycPAFgIKHCERDsosKfCroBSeaDZ5pEAgHMQUOAIZpO2RNpibDq9yCdJWvdlrb0DAQAHIaDAEcw2974EDCgzxuVKkj7eXa/mtg6bRwMAzkBAgSOYRbK5CdQDxTRmWJpG5aSqI2xow65Av7/vw521Ws+sC4A4RUCBI1hN2hJwBsXlcmnG+DxJ0tqdB/r1PQ0t7freU2t17e8/VFsHDd4AxB8CChzBbHOfiEWyUvcyz9od/ZsR+ay6QW0dYTW3hRQ4SIM3APGHgAJHqEvAk4x7MmdQNlUF1NIeOub12/yN1p/ru8IdAMQTAgocIRHb3Pc0Ni9d+VletYXC2lBZd8zrP6tusP4cOEhAARB/CChwhLoE7oMiHVKH0o9lnp4BhRkUAPGIgAJHqE3gIlmTVYfSj0JZZlAAxDsCCmzXEQoraB0UmJg1KJI0c3xnQNlYGVBrx5HrUA40tmp/Y3dhrNnkDgDiCQEFtjNnAFyuzrN4EtUJIzI1PDNFrR1hfVRVf8TrPqtu7PV1kBkUAHGIgALbmTMAOWnJSvIk7l9Jl8uls63txkde5tle09Dra5Z4AMSjxP1tAMcw29wncv2JacY4s2HbkQtlt/k7A0p6ikdS90nQABBPCCiwXaL3QOlpRlcdSsWuOrWH+u4QaxbITisZJokZFADxiYAC2yVym/tDnZyfJV96sg62h/Tx7sPrUAzDsGpQzhrbGWbqCSgA4hABBbZL9Db3PbndLp099sjbjWsaWlV/sF1ul3RmcecMSj27eADEIQIKbMcST29Ha9hm1p+MHZ6hgmyvJJZ4AMQnAgpsV5vgXWQPZTZsW/9l7WH9UMz6k5Pzs6wt2cGD7QqHjegOEgCGGAEFtgsk+Dk8h5o4MluF2alqagvp7U9rer1mBZTCLGV3BZSwITW0dkR9nAAwlAgosB1t7nvzuF26YupoSdJ/b9jT6zWzQPbkgkylJnuUlty51ZjzeADEGwIKbGf28aAGpdu8MzsDyrvbanSgsVWSFA4b2t41g3JKQZYkydd1zwIHKZQFEF8IKLBdbdcSTy41KJaTCrI0ZUyOOsKG/vLRXknSnsBBNbWFlOxxaezwDEndRwOw1RhAvCGgwFahsGH9cqVItrcru5Z5Xuha5jFb3I8fnqnkriMBzIBCN1kA8YaAAlvVH2yX0bUBxZfABwX25ZtnjFaS26XNe+r1WXWDtvm76k8Ks6xrupd4CCgA4gsBBbYye6BkpyYl9EGBfcnNSNGFE/IlSf+9YbdVf3JyfqZ1Tc+txgAQT/iNAFvV0QPlqMxi2Zc27tHWfUFJh86gdN63AN1kAcSZJLsHgMRmtblni3GfLpyQL196sqqDraoOdu7mMXfwSNSgAIhfzKDAVt0HBVJ/0hdvkkeXTxnV42u3inLTra+tgMISD4A4Q0CBrcwtxizxHNm8aWOsP59UkCmP22V9bRbJss0YQLwhoMBWdbS5P6bTx+Ro/IjOvicn52f1es2X1nnf6CQLIN4QUGArc4mHJm1H5nK5VHbBiXK7pEtOK+z1Gp1kAcQrimRhK7NI1kcNylHNmzZGV0wd3Wt5R6KTLID4xQwKbGXNoLDEc0yHhhNJyukKdi3tYbW0h6I9JAAYMgQU2MqsQfERUI5LljfJCi7MogCIJwQU2Mpc4qEG5fi4XC5lp3au1NILBUA8IaDANuGwYXVAHZZBDcrxMmefmEEBEE8IKLBNsKVdYeugQGZQjld3N1l28gCIHwMOKO+9954uv/xyjRo1Si6XSy+99FKv1w3D0L333quRI0cqLS1Ns2bN0vbt23tdU1tbq/nz5ys7O1s+n0/XXXedGhsbB/VBEHtquwpks7xJSkkiKx8vTjQGEI8G/FuhqalJp59+uhYvXtzn6w8//LAee+wxPfnkk1q7dq0yMjI0Z84ctbS0WNfMnz9fW7Zs0YoVK7R8+XK99957uv7664//UyAmWVuMWd4ZFE40BhCPBtwH5dJLL9Wll17a52uGYehXv/qV7r77bn3rW9+SJP3hD39QQUGBXnrpJV111VX69NNP9frrr2vdunWaPn26JOk3v/mNLrvsMv3bv/2bRo0a1efPRvxhi3Fk+DgwEEAciui8+s6dO+X3+zVr1izruZycHM2YMUPl5eWSpPLycvl8PiucSNKsWbPkdru1du3aSA4HDscW48jI6bp/dJMFEE8i2knW7/dLkgoKCno9X1BQYL3m9/uVn5/fexBJScrNzbWuOVRra6taW1utr4PBYCSHDZuYAYUtxoOTwwwKgDgUE5WJixYtUk5OjvUoKiqye0iIANrcR4aPdvcA4lBEA0phYedBZtXV1b2er66utl4rLCxUTU1Nr9c7OjpUW1trXXOou+66S/X19dajqqoqksOGTahBiQwz4BFQAMSTiAaUcePGqbCwUCtXrrSeCwaDWrt2rUpLSyVJpaWlCgQCqqiosK55++23FQ6HNWPGjD5/rtfrVXZ2dq8HYp9Vg8ISz6BY24xZ4gEQRwZcg9LY2KjPP//c+nrnzp3atGmTcnNzVVxcrFtvvVX/+q//qpNOOknjxo3TPffco1GjRumKK66QJE2cOFGXXHKJfvSjH+nJJ59Ue3u7brrpJl111VXs4EkwdU1dbe6ZQRkUTjQGEI8GHFDWr1+vCy+80Pp64cKFkqQFCxZo6dKluuOOO9TU1KTrr79egUBA5513nl5//XWlpqZa3/Pss8/qpptu0sUXXyy326158+bpsccei8DHQSypNdvcU4MyKDldXXiDLe0KhY0+Tz0GgFgz4IBywQUXyDCMI77ucrl0//336/777z/iNbm5uXruuecG+taIM93n8DCDMhjmDIphSA0t7WzbBhAXYmIXD+KPYRicZBwhKUlupad4JFGHAiB+EFBgi2BLh0JdJwWyzXjw2GoMIN4QUGALc4txRopH3iSPzaOJfd3dZAkoAOIDAQW2oM19ZHWfx0O7ewDxgYACW9DmPrI40RhAvCGgwBZmDxTqTyKDZm0A4g0BBbZgBiWycsyAwgwKgDhBQIEtapvMJm0ElEjgRGMA8YaAAluYPVAIKJHh6+omyzZjAPGCgAJbmNuMh2VQgxIJ3Scas4sHQHwgoMAWdc0s8USSjyUeAHGGgAJbUCQbWdlH6CT7UVVAm6oCNowIAAaHgAJbmDUobDOODF+PXTzmYZ47vmrU/3pyta76bbnqmVkBEGMIKIg6wzCsGhRmUCLD7Mjb1hFWS3tYkvTQa39Te8hQS3tY676stXN4ADBgBBREXWNrhzq6DgqkBiUyMlI88rhdkqTAwTZ9uLNWb26ttl5fu/OAXUMDgONCQEHUmV1k05I9Sk3moMBIcLlcVqFsXVO7HvyfrZKk0b40SdLancygAIgtBBREXU1DiyRpGPUnEWV2k3127S59tLteGSkePfG9MyVJn+ypV0MLdSgAYgcBBVG3obJOkjRpVLbNI4kvZjfZ5z6slCTd8PUTNGWMT0W5aQob0vpddXYODwAGhICCqFu7o3O5Yca4PJtHEl/MJR7DkAqyvfqH88dL6r7P5n0HgFhAQEFUhcKGPuzaUTJjfK7No4kvvh4Fx/9n9ilKS+ms75kxrvM+UygLIJYQUBBVn+4LqqGlQ5neJE0ayRJPJJm9UCYUZmnemWOs52eO75xB2by7Xs1tHbaMDQAGioCCqDJ3k0wfO0xJHv76RdJVZxVr9qQCPfK/z7C2HEvSmGFpGpWTqo6woQrqUADECH5DIKrW7uhcZqD+JPJOKczSb6+dfljxscvl0ozx1KEAiC0EFERNmPoT21CHAiDWEFAQNZ/VNCjQ3K70FI8mj86xezgJxZxB+aiqXi3tIZtHAwDHRkBB1JjLC9NKhimZ+pOoGpuXrvwsr9pCYasPDQA4Gb8lEDXm8oK53IDooQ4FQKwhoCAqDMPQhzvN+hMKZO1AHQqAWEJAQVR88VWj9je2yZvk1pQx1J/YYWZXYfLGyoBaO6hDAeBsBBRExZquZYUzi4fJm8QJxnY4YUSmhmemqLUjrI9319s9HAA4KgIKomLtTrYX283lcunsrmWe+/6yRa9/4lcobNg8KgDoGwEFQ84wDBq0OcSC0rFK8bi1ZW9QNzxToYt/+a7+X/mXOtjGkg8AZyGgYMh9eaBZNQ2tSvG4NbXYZ/dwEtqM8Xl6/84LVXbhCcpJS9aXB5p1z8tbdPEv39U2f4PdwwMACwEFQ86cPTmjyKfUZOpP7Jafnarb50zQ6h9fpPsun6TRvjTtrW/Rd55crXVfsgUZgDMQUDDkqD9xpgxvkr5/7jj9zz+dp2klwxRs6dD3/nOt3tzit3toAEBAwdAyDEPvf75fkjST/ieO5EtP0TPXzdCsiflq7Qjrhmcq9McPK+0eFoAER0DBkPqbv0FfNbQqLdmj6WOH2T0cHEFaikdPfm+avju9SGFDuuuFzXpp4x67hwUggRFQMKTe++wrSZ1Nwuh/4mxJHrcemjdZ3z9nrCTpD+Vf2joeAImNgIIh9d72zoDytZNH2DwS9IfL5dINXz9BkrShMqCaYIvNIwKQqAgoGDLNbR1at7Pz5FwCSuwozEnV6UU+SdKKT6vtHQyAhEVAwZBZu6NWbaGwRvvSNH54ht3DwQDMObVAkvTGFgIKAHsQUDBkVn3WvbzjcrlsHg0GYs6phZKk8i/2K9jSbvNoACQiAgqGjFl/8vWTh9s8EgzUCSMydWJ+ptpDht75W81hr2+qCujKxz/Qn9ZX2TA6AImAgIIhsbuuWTu+apLH7VLpCQSUWNS9zNO7cVtHKKw7/vyRNlQGdMefP9av39ouw+DQQQCRRUDBkPjr9s7mbGcU+ZSTlmzzaHA8zGWed7d9pZb27sME/1yxW59VNyrF0/nPx6Nvfaa7X/qEk5EBRBQBBUPC7H/ytZPYvROrJo/O0aicVDW3hfR+V+Bsau3QIys+kyTdcckpeuBbp8rlkp5dW6myZzf0CjIAMBgEFERcRyhstbf/GvUnMcvlcml21yyKuczzu7/uUE1Dq4pz03VNaYmuKR2rxX9/plI8br2+xa8fLFmnjlDYzmEDiBMEFETcR7sDamjpkC89WVPG+OweDgZhdlcdylufVmtf/UH9x6odkqQ7L5lgdQa+bPJIPf3Ds5XpTVL5jgP6LwpnAUQAAQURt+qzztmTc08cLo+b7cWx7OyxuRqWnqy65nb9w9PrdbA9pKnFPl02ubDXdaUn5OmfZ58sSXp0xWdqbO2wY7gA4ggBBRFn1p98nfqTmJfkceviiZ2zKFv2BiVJP7lsYp99bf5+RonGDc/Q/sY2/ceqL6I6TgDxh4CCiKpratPHuwOSpPOpP4kL5m4eSbr0tEJNH5vb53UpSW7deckESZ21Kv56zvEBcPwIKIcwDENhtkseF8Mw9LNXtihsSBNHZmtkTprdQ0IEnH/ScOVlpMib5NYdXQHkSOacWqDpJcPU0h7WL9/cFqURAohHBJQeKnbV6conVuuFjXvsHkpMWlaxWy9t2iuP26V/veJUu4eDCElN9uilsnP12i3na9wxzlRyuVz6ydyJkqQ/b9itrV3LQgAwUASUHtZ/WauNlQH92xvbdLCNfg4D8XlNg3768hZJ0sJvnKxpJX0vAyA2FeWma/yIzH5dO7V4mP5uykgZhrTotU+HeGQA4lWS3QNwkgXnjNUfyndpT+Cgfv/BTpVdeKLdQ4oJLe0hlT27UQfbQzrvxOG68esn2D0k2OzOSybozS3V+uv2/frh0nXyJnX/v9DJBVn6xwtPsLYpA0BfmEHpITXZozsuOUWS9Pg7n+urhlabRxQbfvbKVm2rbtDwzBQ98t3T5WZrccIryk3XgnNKJElv/61Gr33itx6/Xrld1y1dz1ZkAEflMmLwlK9gMKicnBzV19crOzs7oj87HDZ0xeMf6OPd9frezGL96xWTI/rz483Lm/boluc3yeWS/vDDs3U+W4vRpbUjpOUf7VNzj/b3za0d+vXK7WpuC+m00dla8v2zNSLLa71e39yuZRVVamzt0P+aNkZjhqXbMXQAQ2Qgv79tDSiLFy/WL37xC/n9fp1++un6zW9+o7PPPvuY3zeUAUWS1uw4oKt+u0Yet0tv3Po1nZjfv7X3RPOndVW668XNCoUN/eMFJxxzhwcgSR/vDugHS9bpQFObSvLS9Ycfni2XXPr9Bzv1p/VVau6q//K4Xbps8kj96PxxdCQG4kRMBJT/+q//0rXXXqsnn3xSM2bM0K9+9SstW7ZM27ZtU35+/lG/d6gDiiT96A/rtWJrtWZNLNB/Lpg+JO8RqwzD0OPvfqFfvNG5jXTemWP0f+dNVpKHFUP0z879Tbr292tVVXtQWd4kNbV1yNzdP6EwS3mZKfrg8wPW9WePy9X154/XRRPy+7WE2NzWoZc27lXgYJu+PXV0v7e8b90b1Btb/JpQmKXZpxbSCRmIsJgIKDNmzNBZZ52lf//3f5ckhcNhFRUV6eabb9aPf/zjo35vNALKF181avaj7ykUNvTHH81U6Ql5Q/I+sSYUNnT/K1v0dPkuSdKNF5ygO+ac0mdnUeBoahpa9P3fr9PWfZ1bkb928ghdf/54nXtinlwul7buDeo//7pDf/lorzq60sv4ERn60fnj9e2po5WafHiRbU1Di/6wepeeWbtLgeZ2SVKS26Vvnj5K/3D+eE0adfi/F4Zh6L3t+/W793ZYh1xKUnFuuq47b5y+M32M0lPYTwBEguMDSltbm9LT0/XnP/9ZV1xxhfX8ggULFAgE9PLLL/e6vrW1Va2t3QWrwWBQRUVFQxpQJOmelz7R/1uzS+OHZ+jrp1BbIUmf1zTqr9s7/xG/9+8m6YfnjbN5RIhlDS3tennTXk0fO0wTCvv+b3lf/UEt/eBLPbe2Ug1dhbV5GSm65LRCpfTYHXSgsU2vf+JXW9dpysW56SrI9mrdl3XWNeeemKeTC7Ksrw1DKv/igLZVN0iS3K7OoLSpKmAFnJy0ZF02ubDPQATEs2klw/R3U0ZF9Gc6PqDs3btXo0eP1urVq1VaWmo9f8cdd2jVqlVau3Ztr+vvu+8+/exnPzvs5wx1QNnf2KoLfvEuuw0Okexx6Zf/+wx98/TI/sUFjqaxtUP/ta5Kv39/p/YEDh7xujOLfbr+a+P1jUmdSzQf7w7od3/dqVc371PoCF2iM1I8+u5ZxfrBuWNVlJuu5rYO/XfFbv3n+zu160DzUH0kwNH+fkaxfv7tyG4UibuAYtcMitTZvO2dbTVD+h6xxO1y6RuTCihahG06QmG9saVaW/fV93re43Lp66eMOGKTwN11zXp50141t/X+H478rFRdMXW0ctKSD/ueUNjQyk+r9VHX+VJAIjl9jE+zTy089oUDMJCAYsvC6vDhw+XxeFRdXd3r+erqahUWHn4zvF6vvF7vYc9Hw/SxuUc8HA1A9CV53Jo7ZaTmThk5oO8bMyx9wM0XPW6XZp9aGPF/pAEcmy3bLlJSUjRt2jStXLnSei4cDmvlypW9ZlQAAEBisq00feHChVqwYIGmT5+us88+W7/61a/U1NSkH/zgB3YNCQAAOIRtAeW73/2uvvrqK917773y+/0644wz9Prrr6ugoMCuIQEAAIeg1T0AAIiKgfz+pvUnAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHNta3Q+G2fw2GAzaPBIAANBf5u/t/jSxj8mA0tDQIEkqKiqyeSQAAGCgGhoalJOTc9RrYvIsnnA4rL179yorK0sul2tI3ysYDKqoqEhVVVWc+zOEuM/Rw72ODu5zdHCfoyNS99kwDDU0NGjUqFFyu49eZRKTMyhut1tjxoyJ6ntmZ2fzlz8KuM/Rw72ODu5zdHCfoyMS9/lYMycmimQBAIDjEFAAAIDjEFCOwev16qc//am8Xq/dQ4lr3Ofo4V5HB/c5OrjP0WHHfY7JIlkAABDfmEEBAACOQ0ABAACOQ0ABAACOQ0ABAACOQ0A5hsWLF2vs2LFKTU3VjBkz9OGHH9o9pJi2aNEinXXWWcrKylJ+fr6uuOIKbdu2rdc1LS0tKisrU15enjIzMzVv3jxVV1fbNOLY99BDD8nlcunWW2+1nuMeR86ePXv0ve99T3l5eUpLS9PkyZO1fv1663XDMHTvvfdq5MiRSktL06xZs7R9+3YbRxx7QqGQ7rnnHo0bN05paWk64YQT9MADD/Q6z4X7PHDvvfeeLr/8co0aNUoul0svvfRSr9f7c09ra2s1f/58ZWdny+fz6brrrlNjY2NkBmjgiJ5//nkjJSXF+P3vf29s2bLF+NGPfmT4fD6jurra7qHFrDlz5hhLliwxPvnkE2PTpk3GZZddZhQXFxuNjY3WNTfccINRVFRkrFy50li/fr0xc+ZM45xzzrFx1LHrww8/NMaOHWtMmTLFuOWWW6znuceRUVtba5SUlBjf//73jbVr1xo7duww3njjDePzzz+3rnnooYeMnJwc46WXXjI++ugj45vf/KYxbtw44+DBgzaOPLY8+OCDRl5enrF8+XJj586dxrJly4zMzEzj17/+tXUN93ngXn31VeMnP/mJ8cILLxiSjBdffLHX6/25p5dccolx+umnG2vWrDH++te/GieeeKJx9dVXR2R8BJSjOPvss42ysjLr61AoZIwaNcpYtGiRjaOKLzU1NYYkY9WqVYZhGEYgEDCSk5ONZcuWWdd8+umnhiSjvLzcrmHGpIaGBuOkk04yVqxYYXz961+3Agr3OHLuvPNO47zzzjvi6+Fw2CgsLDR+8YtfWM8FAgHD6/Uaf/zjH6MxxLgwd+5c44c//GGv56688kpj/vz5hmFwnyPh0IDSn3u6detWQ5Kxbt0665rXXnvNcLlcxp49ewY9JpZ4jqCtrU0VFRWaNWuW9Zzb7dasWbNUXl5u48jiS319vSQpNzdXklRRUaH29vZe933ChAkqLi7mvg9QWVmZ5s6d2+teStzjSPrLX/6i6dOn6zvf+Y7y8/M1depU/e53v7Ne37lzp/x+f697nZOToxkzZnCvB+Ccc87RypUr9dlnn0mSPvroI73//vu69NJLJXGfh0J/7ml5ebl8Pp+mT59uXTNr1iy53W6tXbt20GOIycMCo2H//v0KhUIqKCjo9XxBQYH+9re/2TSq+BIOh3Xrrbfq3HPP1WmnnSZJ8vv9SklJkc/n63VtQUGB/H6/DaOMTc8//7w2bNigdevWHfYa9zhyduzYoSeeeEILFy7Uv/zLv2jdunX6p3/6J6WkpGjBggXW/ezr3xHudf/9+Mc/VjAY1IQJE+TxeBQKhfTggw9q/vz5ksR9HgL9uad+v1/5+fm9Xk9KSlJubm5E7jsBBbYpKyvTJ598ovfff9/uocSVqqoq3XLLLVqxYoVSU1PtHk5cC4fDmj59un7+859LkqZOnapPPvlETz75pBYsWGDz6OLHn/70Jz377LN67rnndOqpp2rTpk269dZbNWrUKO5zHGOJ5wiGDx8uj8dz2M6G6upqFRYW2jSq+HHTTTdp+fLleueddzRmzBjr+cLCQrW1tSkQCPS6nvvefxUVFaqpqdGZZ56ppKQkJSUladWqVXrssceUlJSkgoIC7nGEjBw5UpMmTer13MSJE1VZWSlJ1v3k35HBuf322/XjH/9YV111lSZPnqxrrrlGt912mxYtWiSJ+zwU+nNPCwsLVVNT0+v1jo4O1dbWRuS+E1COICUlRdOmTdPKlSut58LhsFauXKnS0lIbRxbbDMPQTTfdpBdffFFvv/22xo0b1+v1adOmKTk5udd937ZtmyorK7nv/XTxxRdr8+bN2rRpk/WYPn265s+fb/2ZexwZ55577mHb5D/77DOVlJRIksaNG6fCwsJe9zoYDGrt2rXc6wFobm6W293715XH41E4HJbEfR4K/bmnpaWlCgQCqqiosK55++23FQ6HNWPGjMEPYtBltnHs+eefN7xer7F06VJj69atxvXXX2/4fD7D7/fbPbSYdeONNxo5OTnGu+++a+zbt896NDc3W9fccMMNRnFxsfH2228b69evN0pLS43S0lIbRx37eu7iMQzucaR8+OGHRlJSkvHggw8a27dvN5599lkjPT3deOaZZ6xrHnroIcPn8xkvv/yy8fHHHxvf+ta32P46QAsWLDBGjx5tbTN+4YUXjOHDhxt33HGHdQ33eeAaGhqMjRs3Ghs3bjQkGY888oixceNGY9euXYZh9O+eXnLJJcbUqVONtWvXGu+//75x0kknsc04Wn7zm98YxcXFRkpKinH22Wcba9assXtIMU1Sn48lS5ZY1xw8eND4x3/8R2PYsGFGenq68e1vf9vYt2+ffYOOA4cGFO5x5LzyyivGaaedZni9XmPChAnGb3/7216vh8Nh45577jEKCgoMr9drXHzxxca2bdtsGm1sCgaDxi233GIUFxcbqampxvjx442f/OQnRmtrq3UN93ng3nnnnT7/PV6wYIFhGP27pwcOHDCuvvpqIzMz08jOzjZ+8IMfGA0NDREZn8swerTiAwAAcABqUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOP8f/IuqsYlkXL1AAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.Series(subject_infos['bmi']).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Counter({'female': 3183, 'male': 2724})\n",
      "Counter({'white': 2136, 'black or african american': 1647, 'hispanic': 1413, 'asian': 711})\n",
      "Counter({'not reported': 3246, 'no': 2271, 'yes': 390})\n",
      "Counter({'no': 3216, 'yes': 2625, 'not reported': 66})\n"
     ]
    }
   ],
   "source": [
    "print(Counter(subject_infos['sex']))\n",
    "print(Counter(subject_infos['race']))\n",
    "print(Counter(subject_infos['cur_smoker']))\n",
    "print(Counter(subject_infos['ever_smoked']))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "master",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}