Data analysis_2021.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "RangeIndex: 3663 entries, 0 to 3662\n",
      "Data columns (total 42 columns):\n",
      " #   Column                           Non-Null Count  Dtype  \n",
      "---  ------                           --------------  -----  \n",
      " 0   OSEBuildingID                    3663 non-null   int64  \n",
      " 1   DataYear                         3663 non-null   int64  \n",
      " 2   BuildingName                     3663 non-null   object \n",
      " 3   BuildingType                     3663 non-null   object \n",
      " 4   TaxParcelIdentificationNumber    3663 non-null   object \n",
      " 5   Address                          3663 non-null   object \n",
      " 6   City                             3663 non-null   object \n",
      " 7   State                            3663 non-null   object \n",
      " 8   ZipCode                          3663 non-null   int64  \n",
      " 9   Latitude                         3663 non-null   float64\n",
      " 10  Longitude                        3663 non-null   float64\n",
      " 11  Neighborhood                     3661 non-null   object \n",
      " 12  CouncilDistrictCode              3663 non-null   int64  \n",
      " 13  YearBuilt                        3663 non-null   int64  \n",
      " 14  NumberofFloors                   3663 non-null   int64  \n",
      " 15  NumberofBuildings                3663 non-null   int64  \n",
      " 16  PropertyGFATotal                 3663 non-null   int64  \n",
      " 17  PropertyGFABuilding(s)           3663 non-null   int64  \n",
      " 18  PropertyGFAParking               3663 non-null   int64  \n",
      " 19  ENERGYSTARScore                  3663 non-null   int64  \n",
      " 20  SiteEUIWN(kBtu/sf)               3663 non-null   float64\n",
      " 21  SiteEUI(kBtu/sf)                 3663 non-null   float64\n",
      " 22  SiteEnergyUse(kBtu)              3663 non-null   int64  \n",
      " 23  SiteEnergyUseWN(kBtu)            3663 non-null   int64  \n",
      " 24  SourceEUIWN(kBtu/sf)             3663 non-null   float64\n",
      " 25  SourceEUI(kBtu/sf)               3663 non-null   float64\n",
      " 26  EPAPropertyType                  3565 non-null   object \n",
      " 27  LargestPropertyUseType           3607 non-null   object \n",
      " 28  LargestPropertyUseTypeGFA        3663 non-null   int64  \n",
      " 29  SecondLargestPropertyUseType     1994 non-null   object \n",
      " 30  SecondLargestPropertyUseTypeGFA  3663 non-null   int64  \n",
      " 31  ThirdLargestPropertyUseType      722 non-null    object \n",
      " 32  ThirdLargestPropertyUseTypeGFA   3663 non-null   int64  \n",
      " 33  Electricity(kWh)                 3663 non-null   int64  \n",
      " 34  SteamUse(kBtu)                   3663 non-null   int64  \n",
      " 35  NaturalGas(therms)               3663 non-null   int64  \n",
      " 36  ComplianceStatus                 3663 non-null   object \n",
      " 37  ComplianceIssue                  3663 non-null   object \n",
      " 38  Electricity(kBtu)                3663 non-null   int64  \n",
      " 39  NaturalGas(kBtu)                 3663 non-null   int64  \n",
      " 40  TotalGHGEmissions                3663 non-null   float64\n",
      " 41  GHGEmissionsIntensity            3663 non-null   float64\n",
      "dtypes: float64(8), int64(21), object(13)\n",
      "memory usage: 1.2+ MB\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(   OSEBuildingID  DataYear                    BuildingName    BuildingType   \n",
       " 0              1      2021            MAYFLOWER PARK HOTEL  NonResidential  \\\n",
       " 1              2      2021                 PARAMOUNT HOTEL  NonResidential   \n",
       " 2              3      2021  WESTIN HOTEL (Parent Building)  NonResidential   \n",
       " 3              5      2021                       HOTEL MAX  NonResidential   \n",
       " 4              8      2021           WARWICK SEATTLE HOTEL  NonResidential   \n",
       " \n",
       "   TaxParcelIdentificationNumber         Address     City State  ZipCode   \n",
       " 0                     659000030   405 OLIVE WAY  SEATTLE    WA    98101  \\\n",
       " 1                     659000220     724 PINE ST  SEATTLE    WA    98101   \n",
       " 2                     659000475    1900 5TH AVE  SEATTLE    WA    98101   \n",
       " 3                     659000640  620 STEWART ST  SEATTLE    WA    98101   \n",
       " 4                     659000970   401 LENORA ST  SEATTLE    WA    98121   \n",
       " \n",
       "    Latitude  Longitude Neighborhood  CouncilDistrictCode  YearBuilt   \n",
       " 0  47.61220 -122.33799     DOWNTOWN                    1       1927  \\\n",
       " 1  47.61307 -122.33361     DOWNTOWN                    1       1996   \n",
       " 2  47.61367 -122.33822     DOWNTOWN                    1       1969   \n",
       " 3  47.61412 -122.33664     DOWNTOWN                    1       1926   \n",
       " 4  47.61375 -122.34047     DOWNTOWN                    1       1980   \n",
       " \n",
       "    NumberofFloors  NumberofBuildings  PropertyGFATotal   \n",
       " 0              12                  1             88434  \\\n",
       " 1              11                  1            103566   \n",
       " 2              41                  3            956110   \n",
       " 3              10                  1             61320   \n",
       " 4              18                  1            175580   \n",
       " \n",
       "    PropertyGFABuilding(s)  PropertyGFAParking  ENERGYSTARScore   \n",
       " 0                   88434                   0               78  \\\n",
       " 1                   88502               15064               96   \n",
       " 2                  759392              196718               96   \n",
       " 3                   61320                   0               76   \n",
       " 4                  113580               62000               90   \n",
       " \n",
       "    SiteEUIWN(kBtu/sf)  SiteEUI(kBtu/sf)  SiteEnergyUse(kBtu)   \n",
       " 0           73.800003         73.599998              6510477  \\\n",
       " 1           52.200001         52.200001              4617104   \n",
       " 2           46.500000         46.500000             43953212   \n",
       " 3           79.500000         79.500000              4873753   \n",
       " 4           92.400002         92.000000             11358936   \n",
       " \n",
       "    SiteEnergyUseWN(kBtu)  SourceEUIWN(kBtu/sf)  SourceEUI(kBtu/sf)   \n",
       " 0                6522024            144.399994          144.199997  \\\n",
       " 1                4617104             99.099998           99.099998   \n",
       " 2               43953212            105.300003          105.300003   \n",
       " 3                4873753            136.000000          136.000000   \n",
       " 4               11409090            161.500000          161.100006   \n",
       " \n",
       "   EPAPropertyType LargestPropertyUseType  LargestPropertyUseTypeGFA   \n",
       " 0           Hotel                  Hotel                      88434  \\\n",
       " 1           Hotel                  Hotel                      88502   \n",
       " 2           Hotel                  Hotel                     945349   \n",
       " 3           Hotel                  Hotel                      61320   \n",
       " 4           Hotel                  Hotel                     123445   \n",
       " \n",
       "   SecondLargestPropertyUseType  SecondLargestPropertyUseTypeGFA   \n",
       " 0                          NaN                                0  \\\n",
       " 1                      Parking                            15064   \n",
       " 2                      Parking                           117783   \n",
       " 3                          NaN                                0   \n",
       " 4                      Parking                            68009   \n",
       " \n",
       "   ThirdLargestPropertyUseType  ThirdLargestPropertyUseTypeGFA   \n",
       " 0                         NaN                               0  \\\n",
       " 1                         NaN                               0   \n",
       " 2               Swimming Pool                               0   \n",
       " 3                         NaN                               0   \n",
       " 4               Swimming Pool                               0   \n",
       " \n",
       "    Electricity(kWh)  SteamUse(kBtu)  NaturalGas(therms) ComplianceStatus   \n",
       " 0            944955         1798672               14876        Compliant  \\\n",
       " 1            657478               0               23738        Compliant   \n",
       " 2           8673722        10583473               37750        Compliant   \n",
       " 3            509497         1167770               19676        Compliant   \n",
       " 4           1333597               0               68087        Compliant   \n",
       " \n",
       "   ComplianceIssue  Electricity(kBtu)  NaturalGas(kBtu)  TotalGHGEmissions   \n",
       " 0        No Issue            3224187           1487620              241.6  \\\n",
       " 1        No Issue            2243315           2373790              135.4   \n",
       " 2        No Issue           29594739           3775000             1201.4   \n",
       " 3        No Issue            1738403           1967580              208.6   \n",
       " 4        No Issue            4550233           6808700              380.4   \n",
       " \n",
       "    GHGEmissionsIntensity  \n",
       " 0                    2.7  \n",
       " 1                    1.5  \n",
       " 2                    1.6  \n",
       " 3                    3.4  \n",
       " 4                    3.3  ,\n",
       " None)"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "pd.set_option('display.max_columns', None)\n",
    "pd.set_option('display.max_rows', None)\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "\n",
    "file_path = \"2021_Building_Energy_Benchmarking_20240322.csv\"\n",
    "data = pd.read_csv(file_path)\n",
    "\n",
    "data.head(), data.info()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>OSEBuildingID</th>\n",
       "      <th>DataYear</th>\n",
       "      <th>BuildingName</th>\n",
       "      <th>BuildingType</th>\n",
       "      <th>TaxParcelIdentificationNumber</th>\n",
       "      <th>Address</th>\n",
       "      <th>City</th>\n",
       "      <th>State</th>\n",
       "      <th>ZipCode</th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Neighborhood</th>\n",
       "      <th>CouncilDistrictCode</th>\n",
       "      <th>YearBuilt</th>\n",
       "      <th>NumberofFloors</th>\n",
       "      <th>NumberofBuildings</th>\n",
       "      <th>PropertyGFATotal</th>\n",
       "      <th>PropertyGFABuilding(s)</th>\n",
       "      <th>PropertyGFAParking</th>\n",
       "      <th>ENERGYSTARScore</th>\n",
       "      <th>SiteEUIWN(kBtu/sf)</th>\n",
       "      <th>SiteEUI(kBtu/sf)</th>\n",
       "      <th>SiteEnergyUse(kBtu)</th>\n",
       "      <th>SiteEnergyUseWN(kBtu)</th>\n",
       "      <th>SourceEUIWN(kBtu/sf)</th>\n",
       "      <th>SourceEUI(kBtu/sf)</th>\n",
       "      <th>EPAPropertyType</th>\n",
       "      <th>LargestPropertyUseType</th>\n",
       "      <th>LargestPropertyUseTypeGFA</th>\n",
       "      <th>SecondLargestPropertyUseType</th>\n",
       "      <th>SecondLargestPropertyUseTypeGFA</th>\n",
       "      <th>ThirdLargestPropertyUseType</th>\n",
       "      <th>ThirdLargestPropertyUseTypeGFA</th>\n",
       "      <th>Electricity(kWh)</th>\n",
       "      <th>SteamUse(kBtu)</th>\n",
       "      <th>NaturalGas(therms)</th>\n",
       "      <th>ComplianceStatus</th>\n",
       "      <th>ComplianceIssue</th>\n",
       "      <th>Electricity(kBtu)</th>\n",
       "      <th>NaturalGas(kBtu)</th>\n",
       "      <th>TotalGHGEmissions</th>\n",
       "      <th>GHGEmissionsIntensity</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>MAYFLOWER PARK HOTEL</td>\n",
       "      <td>NonResidential</td>\n",
       "      <td>659000030</td>\n",
       "      <td>405 OLIVE WAY</td>\n",
       "      <td>SEATTLE</td>\n",
       "      <td>WA</td>\n",
       "      <td>98101</td>\n",
       "      <td>47.61220</td>\n",
       "      <td>-122.33799</td>\n",
       "      <td>DOWNTOWN</td>\n",
       "      <td>1</td>\n",
       "      <td>1927</td>\n",
       "      <td>12</td>\n",
       "      <td>1</td>\n",
       "      <td>88434</td>\n",
       "      <td>88434</td>\n",
       "      <td>0</td>\n",
       "      <td>78</td>\n",
       "      <td>73.800003</td>\n",
       "      <td>73.599998</td>\n",
       "      <td>6510477</td>\n",
       "      <td>6522024</td>\n",
       "      <td>144.399994</td>\n",
       "      <td>144.199997</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>88434</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0</td>\n",
       "      <td>944955</td>\n",
       "      <td>1798672</td>\n",
       "      <td>14876</td>\n",
       "      <td>Compliant</td>\n",
       "      <td>No Issue</td>\n",
       "      <td>3224187</td>\n",
       "      <td>1487620</td>\n",
       "      <td>241.6</td>\n",
       "      <td>2.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>PARAMOUNT HOTEL</td>\n",
       "      <td>NonResidential</td>\n",
       "      <td>659000220</td>\n",
       "      <td>724 PINE ST</td>\n",
       "      <td>SEATTLE</td>\n",
       "      <td>WA</td>\n",
       "      <td>98101</td>\n",
       "      <td>47.61307</td>\n",
       "      <td>-122.33361</td>\n",
       "      <td>DOWNTOWN</td>\n",
       "      <td>1</td>\n",
       "      <td>1996</td>\n",
       "      <td>11</td>\n",
       "      <td>1</td>\n",
       "      <td>103566</td>\n",
       "      <td>88502</td>\n",
       "      <td>15064</td>\n",
       "      <td>96</td>\n",
       "      <td>52.200001</td>\n",
       "      <td>52.200001</td>\n",
       "      <td>4617104</td>\n",
       "      <td>4617104</td>\n",
       "      <td>99.099998</td>\n",
       "      <td>99.099998</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>88502</td>\n",
       "      <td>Parking</td>\n",
       "      <td>15064</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0</td>\n",
       "      <td>657478</td>\n",
       "      <td>0</td>\n",
       "      <td>23738</td>\n",
       "      <td>Compliant</td>\n",
       "      <td>No Issue</td>\n",
       "      <td>2243315</td>\n",
       "      <td>2373790</td>\n",
       "      <td>135.4</td>\n",
       "      <td>1.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2021</td>\n",
       "      <td>WESTIN HOTEL (Parent Building)</td>\n",
       "      <td>NonResidential</td>\n",
       "      <td>659000475</td>\n",
       "      <td>1900 5TH AVE</td>\n",
       "      <td>SEATTLE</td>\n",
       "      <td>WA</td>\n",
       "      <td>98101</td>\n",
       "      <td>47.61367</td>\n",
       "      <td>-122.33822</td>\n",
       "      <td>DOWNTOWN</td>\n",
       "      <td>1</td>\n",
       "      <td>1969</td>\n",
       "      <td>41</td>\n",
       "      <td>3</td>\n",
       "      <td>956110</td>\n",
       "      <td>759392</td>\n",
       "      <td>196718</td>\n",
       "      <td>96</td>\n",
       "      <td>46.500000</td>\n",
       "      <td>46.500000</td>\n",
       "      <td>43953212</td>\n",
       "      <td>43953212</td>\n",
       "      <td>105.300003</td>\n",
       "      <td>105.300003</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>945349</td>\n",
       "      <td>Parking</td>\n",
       "      <td>117783</td>\n",
       "      <td>Swimming Pool</td>\n",
       "      <td>0</td>\n",
       "      <td>8673722</td>\n",
       "      <td>10583473</td>\n",
       "      <td>37750</td>\n",
       "      <td>Compliant</td>\n",
       "      <td>No Issue</td>\n",
       "      <td>29594739</td>\n",
       "      <td>3775000</td>\n",
       "      <td>1201.4</td>\n",
       "      <td>1.6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>5</td>\n",
       "      <td>2021</td>\n",
       "      <td>HOTEL MAX</td>\n",
       "      <td>NonResidential</td>\n",
       "      <td>659000640</td>\n",
       "      <td>620 STEWART ST</td>\n",
       "      <td>SEATTLE</td>\n",
       "      <td>WA</td>\n",
       "      <td>98101</td>\n",
       "      <td>47.61412</td>\n",
       "      <td>-122.33664</td>\n",
       "      <td>DOWNTOWN</td>\n",
       "      <td>1</td>\n",
       "      <td>1926</td>\n",
       "      <td>10</td>\n",
       "      <td>1</td>\n",
       "      <td>61320</td>\n",
       "      <td>61320</td>\n",
       "      <td>0</td>\n",
       "      <td>76</td>\n",
       "      <td>79.500000</td>\n",
       "      <td>79.500000</td>\n",
       "      <td>4873753</td>\n",
       "      <td>4873753</td>\n",
       "      <td>136.000000</td>\n",
       "      <td>136.000000</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>61320</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0</td>\n",
       "      <td>509497</td>\n",
       "      <td>1167770</td>\n",
       "      <td>19676</td>\n",
       "      <td>Compliant</td>\n",
       "      <td>No Issue</td>\n",
       "      <td>1738403</td>\n",
       "      <td>1967580</td>\n",
       "      <td>208.6</td>\n",
       "      <td>3.4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>8</td>\n",
       "      <td>2021</td>\n",
       "      <td>WARWICK SEATTLE HOTEL</td>\n",
       "      <td>NonResidential</td>\n",
       "      <td>659000970</td>\n",
       "      <td>401 LENORA ST</td>\n",
       "      <td>SEATTLE</td>\n",
       "      <td>WA</td>\n",
       "      <td>98121</td>\n",
       "      <td>47.61375</td>\n",
       "      <td>-122.34047</td>\n",
       "      <td>DOWNTOWN</td>\n",
       "      <td>1</td>\n",
       "      <td>1980</td>\n",
       "      <td>18</td>\n",
       "      <td>1</td>\n",
       "      <td>175580</td>\n",
       "      <td>113580</td>\n",
       "      <td>62000</td>\n",
       "      <td>90</td>\n",
       "      <td>92.400002</td>\n",
       "      <td>92.000000</td>\n",
       "      <td>11358936</td>\n",
       "      <td>11409090</td>\n",
       "      <td>161.500000</td>\n",
       "      <td>161.100006</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>Hotel</td>\n",
       "      <td>123445</td>\n",
       "      <td>Parking</td>\n",
       "      <td>68009</td>\n",
       "      <td>Swimming Pool</td>\n",
       "      <td>0</td>\n",
       "      <td>1333597</td>\n",
       "      <td>0</td>\n",
       "      <td>68087</td>\n",
       "      <td>Compliant</td>\n",
       "      <td>No Issue</td>\n",
       "      <td>4550233</td>\n",
       "      <td>6808700</td>\n",
       "      <td>380.4</td>\n",
       "      <td>3.3</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   OSEBuildingID  DataYear                    BuildingName    BuildingType   \n",
       "0              1      2021            MAYFLOWER PARK HOTEL  NonResidential  \\\n",
       "1              2      2021                 PARAMOUNT HOTEL  NonResidential   \n",
       "2              3      2021  WESTIN HOTEL (Parent Building)  NonResidential   \n",
       "3              5      2021                       HOTEL MAX  NonResidential   \n",
       "4              8      2021           WARWICK SEATTLE HOTEL  NonResidential   \n",
       "\n",
       "  TaxParcelIdentificationNumber         Address     City State  ZipCode   \n",
       "0                     659000030   405 OLIVE WAY  SEATTLE    WA    98101  \\\n",
       "1                     659000220     724 PINE ST  SEATTLE    WA    98101   \n",
       "2                     659000475    1900 5TH AVE  SEATTLE    WA    98101   \n",
       "3                     659000640  620 STEWART ST  SEATTLE    WA    98101   \n",
       "4                     659000970   401 LENORA ST  SEATTLE    WA    98121   \n",
       "\n",
       "   Latitude  Longitude Neighborhood  CouncilDistrictCode  YearBuilt   \n",
       "0  47.61220 -122.33799     DOWNTOWN                    1       1927  \\\n",
       "1  47.61307 -122.33361     DOWNTOWN                    1       1996   \n",
       "2  47.61367 -122.33822     DOWNTOWN                    1       1969   \n",
       "3  47.61412 -122.33664     DOWNTOWN                    1       1926   \n",
       "4  47.61375 -122.34047     DOWNTOWN                    1       1980   \n",
       "\n",
       "   NumberofFloors  NumberofBuildings  PropertyGFATotal   \n",
       "0              12                  1             88434  \\\n",
       "1              11                  1            103566   \n",
       "2              41                  3            956110   \n",
       "3              10                  1             61320   \n",
       "4              18                  1            175580   \n",
       "\n",
       "   PropertyGFABuilding(s)  PropertyGFAParking  ENERGYSTARScore   \n",
       "0                   88434                   0               78  \\\n",
       "1                   88502               15064               96   \n",
       "2                  759392              196718               96   \n",
       "3                   61320                   0               76   \n",
       "4                  113580               62000               90   \n",
       "\n",
       "   SiteEUIWN(kBtu/sf)  SiteEUI(kBtu/sf)  SiteEnergyUse(kBtu)   \n",
       "0           73.800003         73.599998              6510477  \\\n",
       "1           52.200001         52.200001              4617104   \n",
       "2           46.500000         46.500000             43953212   \n",
       "3           79.500000         79.500000              4873753   \n",
       "4           92.400002         92.000000             11358936   \n",
       "\n",
       "   SiteEnergyUseWN(kBtu)  SourceEUIWN(kBtu/sf)  SourceEUI(kBtu/sf)   \n",
       "0                6522024            144.399994          144.199997  \\\n",
       "1                4617104             99.099998           99.099998   \n",
       "2               43953212            105.300003          105.300003   \n",
       "3                4873753            136.000000          136.000000   \n",
       "4               11409090            161.500000          161.100006   \n",
       "\n",
       "  EPAPropertyType LargestPropertyUseType  LargestPropertyUseTypeGFA   \n",
       "0           Hotel                  Hotel                      88434  \\\n",
       "1           Hotel                  Hotel                      88502   \n",
       "2           Hotel                  Hotel                     945349   \n",
       "3           Hotel                  Hotel                      61320   \n",
       "4           Hotel                  Hotel                     123445   \n",
       "\n",
       "  SecondLargestPropertyUseType  SecondLargestPropertyUseTypeGFA   \n",
       "0                          NaN                                0  \\\n",
       "1                      Parking                            15064   \n",
       "2                      Parking                           117783   \n",
       "3                          NaN                                0   \n",
       "4                      Parking                            68009   \n",
       "\n",
       "  ThirdLargestPropertyUseType  ThirdLargestPropertyUseTypeGFA   \n",
       "0                         NaN                               0  \\\n",
       "1                         NaN                               0   \n",
       "2               Swimming Pool                               0   \n",
       "3                         NaN                               0   \n",
       "4               Swimming Pool                               0   \n",
       "\n",
       "   Electricity(kWh)  SteamUse(kBtu)  NaturalGas(therms) ComplianceStatus   \n",
       "0            944955         1798672               14876        Compliant  \\\n",
       "1            657478               0               23738        Compliant   \n",
       "2           8673722        10583473               37750        Compliant   \n",
       "3            509497         1167770               19676        Compliant   \n",
       "4           1333597               0               68087        Compliant   \n",
       "\n",
       "  ComplianceIssue  Electricity(kBtu)  NaturalGas(kBtu)  TotalGHGEmissions   \n",
       "0        No Issue            3224187           1487620              241.6  \\\n",
       "1        No Issue            2243315           2373790              135.4   \n",
       "2        No Issue           29594739           3775000             1201.4   \n",
       "3        No Issue            1738403           1967580              208.6   \n",
       "4        No Issue            4550233           6808700              380.4   \n",
       "\n",
       "   GHGEmissionsIntensity  \n",
       "0                    2.7  \n",
       "1                    1.5  \n",
       "2                    1.6  \n",
       "3                    3.4  \n",
       "4                    3.3  "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>OSEBuildingID</th>\n",
       "      <th>DataYear</th>\n",
       "      <th>ZipCode</th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>CouncilDistrictCode</th>\n",
       "      <th>YearBuilt</th>\n",
       "      <th>NumberofFloors</th>\n",
       "      <th>NumberofBuildings</th>\n",
       "      <th>PropertyGFATotal</th>\n",
       "      <th>PropertyGFABuilding(s)</th>\n",
       "      <th>PropertyGFAParking</th>\n",
       "      <th>ENERGYSTARScore</th>\n",
       "      <th>SiteEUIWN(kBtu/sf)</th>\n",
       "      <th>SiteEUI(kBtu/sf)</th>\n",
       "      <th>SiteEnergyUse(kBtu)</th>\n",
       "      <th>SiteEnergyUseWN(kBtu)</th>\n",
       "      <th>SourceEUIWN(kBtu/sf)</th>\n",
       "      <th>SourceEUI(kBtu/sf)</th>\n",
       "      <th>LargestPropertyUseTypeGFA</th>\n",
       "      <th>SecondLargestPropertyUseTypeGFA</th>\n",
       "      <th>ThirdLargestPropertyUseTypeGFA</th>\n",
       "      <th>Electricity(kWh)</th>\n",
       "      <th>SteamUse(kBtu)</th>\n",
       "      <th>NaturalGas(therms)</th>\n",
       "      <th>Electricity(kBtu)</th>\n",
       "      <th>NaturalGas(kBtu)</th>\n",
       "      <th>TotalGHGEmissions</th>\n",
       "      <th>GHGEmissionsIntensity</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.0</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3663.00000</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3.663000e+03</td>\n",
       "      <td>3663.000000</td>\n",
       "      <td>3663.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>23922.210483</td>\n",
       "      <td>2021.0</td>\n",
       "      <td>96643.509419</td>\n",
       "      <td>47.598941</td>\n",
       "      <td>-122.268003</td>\n",
       "      <td>2.689326</td>\n",
       "      <td>1973.129675</td>\n",
       "      <td>5.056784</td>\n",
       "      <td>1.166530</td>\n",
       "      <td>1.070735e+05</td>\n",
       "      <td>9.098446e+04</td>\n",
       "      <td>16089.067977</td>\n",
       "      <td>52.903631</td>\n",
       "      <td>50.093885</td>\n",
       "      <td>50.728774</td>\n",
       "      <td>5.001209e+06</td>\n",
       "      <td>4.936628e+06</td>\n",
       "      <td>107.170543</td>\n",
       "      <td>108.467049</td>\n",
       "      <td>8.663969e+04</td>\n",
       "      <td>17906.29484</td>\n",
       "      <td>2743.545182</td>\n",
       "      <td>9.340146e+05</td>\n",
       "      <td>2.925918e+05</td>\n",
       "      <td>1.521752e+04</td>\n",
       "      <td>3.529597e+06</td>\n",
       "      <td>1.550744e+06</td>\n",
       "      <td>118.291974</td>\n",
       "      <td>1.217882</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>14331.756851</td>\n",
       "      <td>0.0</td>\n",
       "      <td>11933.848036</td>\n",
       "      <td>1.113665</td>\n",
       "      <td>2.858293</td>\n",
       "      <td>2.229124</td>\n",
       "      <td>34.712470</td>\n",
       "      <td>5.701001</td>\n",
       "      <td>2.180697</td>\n",
       "      <td>2.332236e+05</td>\n",
       "      <td>2.092554e+05</td>\n",
       "      <td>47304.388890</td>\n",
       "      <td>39.506607</td>\n",
       "      <td>53.522627</td>\n",
       "      <td>53.746031</td>\n",
       "      <td>1.556745e+07</td>\n",
       "      <td>1.546751e+07</td>\n",
       "      <td>111.018616</td>\n",
       "      <td>111.240183</td>\n",
       "      <td>2.832689e+05</td>\n",
       "      <td>47019.72011</td>\n",
       "      <td>14642.040441</td>\n",
       "      <td>2.570627e+06</td>\n",
       "      <td>4.454561e+06</td>\n",
       "      <td>7.531615e+04</td>\n",
       "      <td>1.756966e+07</td>\n",
       "      <td>7.594444e+06</td>\n",
       "      <td>568.240613</td>\n",
       "      <td>2.145308</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>1.000000</td>\n",
       "      <td>2021.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-122.414250</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1896.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.000000e+04</td>\n",
       "      <td>1.111100e+04</td>\n",
       "      <td>-3.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>20296.500000</td>\n",
       "      <td>2021.0</td>\n",
       "      <td>98104.000000</td>\n",
       "      <td>47.601425</td>\n",
       "      <td>-122.350505</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1953.000000</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>2.964100e+04</td>\n",
       "      <td>2.767550e+04</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>26.400000</td>\n",
       "      <td>26.750000</td>\n",
       "      <td>8.901380e+05</td>\n",
       "      <td>8.712810e+05</td>\n",
       "      <td>65.099998</td>\n",
       "      <td>65.900002</td>\n",
       "      <td>2.586550e+04</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.764520e+05</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>6.147670e+05</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>5.200000</td>\n",
       "      <td>0.100000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>23615.000000</td>\n",
       "      <td>2021.0</td>\n",
       "      <td>98112.000000</td>\n",
       "      <td>47.619330</td>\n",
       "      <td>-122.332320</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1980.000000</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>4.734300e+04</td>\n",
       "      <td>4.407500e+04</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>67.000000</td>\n",
       "      <td>35.500000</td>\n",
       "      <td>35.900002</td>\n",
       "      <td>1.806093e+06</td>\n",
       "      <td>1.780499e+06</td>\n",
       "      <td>82.900002</td>\n",
       "      <td>83.800003</td>\n",
       "      <td>4.160000e+04</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>3.304030e+05</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>3.501000e+03</td>\n",
       "      <td>1.148932e+06</td>\n",
       "      <td>3.574000e+05</td>\n",
       "      <td>29.200000</td>\n",
       "      <td>0.600000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>26976.000000</td>\n",
       "      <td>2021.0</td>\n",
       "      <td>98122.000000</td>\n",
       "      <td>47.657165</td>\n",
       "      <td>-122.318585</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>2002.000000</td>\n",
       "      <td>6.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.015550e+05</td>\n",
       "      <td>8.893200e+04</td>\n",
       "      <td>9581.500000</td>\n",
       "      <td>90.000000</td>\n",
       "      <td>55.049999</td>\n",
       "      <td>55.599998</td>\n",
       "      <td>4.248009e+06</td>\n",
       "      <td>4.201858e+06</td>\n",
       "      <td>113.850002</td>\n",
       "      <td>114.599998</td>\n",
       "      <td>8.230900e+04</td>\n",
       "      <td>14353.50000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>7.760675e+05</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.363000e+04</td>\n",
       "      <td>2.718398e+06</td>\n",
       "      <td>1.386725e+06</td>\n",
       "      <td>92.950000</td>\n",
       "      <td>1.400000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>50858.000000</td>\n",
       "      <td>2021.0</td>\n",
       "      <td>98199.000000</td>\n",
       "      <td>47.733870</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>2021.000000</td>\n",
       "      <td>76.000000</td>\n",
       "      <td>111.000000</td>\n",
       "      <td>9.320156e+06</td>\n",
       "      <td>9.320156e+06</td>\n",
       "      <td>686750.000000</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>908.900024</td>\n",
       "      <td>909.900024</td>\n",
       "      <td>5.180035e+08</td>\n",
       "      <td>5.132408e+08</td>\n",
       "      <td>2545.000000</td>\n",
       "      <td>2547.600098</td>\n",
       "      <td>1.521647e+07</td>\n",
       "      <td>750000.00000</td>\n",
       "      <td>480625.000000</td>\n",
       "      <td>7.925876e+07</td>\n",
       "      <td>1.680307e+08</td>\n",
       "      <td>3.737314e+06</td>\n",
       "      <td>9.191679e+08</td>\n",
       "      <td>3.737314e+08</td>\n",
       "      <td>20445.700000</td>\n",
       "      <td>47.400000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       OSEBuildingID  DataYear       ZipCode     Latitude    Longitude   \n",
       "count    3663.000000    3663.0   3663.000000  3663.000000  3663.000000  \\\n",
       "mean    23922.210483    2021.0  96643.509419    47.598941  -122.268003   \n",
       "std     14331.756851       0.0  11933.848036     1.113665     2.858293   \n",
       "min         1.000000    2021.0      0.000000     0.000000  -122.414250   \n",
       "25%     20296.500000    2021.0  98104.000000    47.601425  -122.350505   \n",
       "50%     23615.000000    2021.0  98112.000000    47.619330  -122.332320   \n",
       "75%     26976.000000    2021.0  98122.000000    47.657165  -122.318585   \n",
       "max     50858.000000    2021.0  98199.000000    47.733870     0.000000   \n",
       "\n",
       "       CouncilDistrictCode    YearBuilt  NumberofFloors  NumberofBuildings   \n",
       "count          3663.000000  3663.000000     3663.000000        3663.000000  \\\n",
       "mean              2.689326  1973.129675        5.056784           1.166530   \n",
       "std               2.229124    34.712470        5.701001           2.180697   \n",
       "min               0.000000  1896.000000        0.000000           0.000000   \n",
       "25%               1.000000  1953.000000        3.000000           1.000000   \n",
       "50%               1.000000  1980.000000        4.000000           1.000000   \n",
       "75%               4.000000  2002.000000        6.000000           1.000000   \n",
       "max               7.000000  2021.000000       76.000000         111.000000   \n",
       "\n",
       "       PropertyGFATotal  PropertyGFABuilding(s)  PropertyGFAParking   \n",
       "count      3.663000e+03            3.663000e+03         3663.000000  \\\n",
       "mean       1.070735e+05            9.098446e+04        16089.067977   \n",
       "std        2.332236e+05            2.092554e+05        47304.388890   \n",
       "min        2.000000e+04            1.111100e+04           -3.000000   \n",
       "25%        2.964100e+04            2.767550e+04            0.000000   \n",
       "50%        4.734300e+04            4.407500e+04            0.000000   \n",
       "75%        1.015550e+05            8.893200e+04         9581.500000   \n",
       "max        9.320156e+06            9.320156e+06       686750.000000   \n",
       "\n",
       "       ENERGYSTARScore  SiteEUIWN(kBtu/sf)  SiteEUI(kBtu/sf)   \n",
       "count      3663.000000         3663.000000       3663.000000  \\\n",
       "mean         52.903631           50.093885         50.728774   \n",
       "std          39.506607           53.522627         53.746031   \n",
       "min           0.000000            0.000000          0.000000   \n",
       "25%           0.000000           26.400000         26.750000   \n",
       "50%          67.000000           35.500000         35.900002   \n",
       "75%          90.000000           55.049999         55.599998   \n",
       "max         100.000000          908.900024        909.900024   \n",
       "\n",
       "       SiteEnergyUse(kBtu)  SiteEnergyUseWN(kBtu)  SourceEUIWN(kBtu/sf)   \n",
       "count         3.663000e+03           3.663000e+03           3663.000000  \\\n",
       "mean          5.001209e+06           4.936628e+06            107.170543   \n",
       "std           1.556745e+07           1.546751e+07            111.018616   \n",
       "min           0.000000e+00           0.000000e+00              0.000000   \n",
       "25%           8.901380e+05           8.712810e+05             65.099998   \n",
       "50%           1.806093e+06           1.780499e+06             82.900002   \n",
       "75%           4.248009e+06           4.201858e+06            113.850002   \n",
       "max           5.180035e+08           5.132408e+08           2545.000000   \n",
       "\n",
       "       SourceEUI(kBtu/sf)  LargestPropertyUseTypeGFA   \n",
       "count         3663.000000               3.663000e+03  \\\n",
       "mean           108.467049               8.663969e+04   \n",
       "std            111.240183               2.832689e+05   \n",
       "min              0.000000               0.000000e+00   \n",
       "25%             65.900002               2.586550e+04   \n",
       "50%             83.800003               4.160000e+04   \n",
       "75%            114.599998               8.230900e+04   \n",
       "max           2547.600098               1.521647e+07   \n",
       "\n",
       "       SecondLargestPropertyUseTypeGFA  ThirdLargestPropertyUseTypeGFA   \n",
       "count                       3663.00000                     3663.000000  \\\n",
       "mean                       17906.29484                     2743.545182   \n",
       "std                        47019.72011                    14642.040441   \n",
       "min                            0.00000                        0.000000   \n",
       "25%                            0.00000                        0.000000   \n",
       "50%                            0.00000                        0.000000   \n",
       "75%                        14353.50000                        0.000000   \n",
       "max                       750000.00000                   480625.000000   \n",
       "\n",
       "       Electricity(kWh)  SteamUse(kBtu)  NaturalGas(therms)   \n",
       "count      3.663000e+03    3.663000e+03        3.663000e+03  \\\n",
       "mean       9.340146e+05    2.925918e+05        1.521752e+04   \n",
       "std        2.570627e+06    4.454561e+06        7.531615e+04   \n",
       "min        0.000000e+00    0.000000e+00        0.000000e+00   \n",
       "25%        1.764520e+05    0.000000e+00        0.000000e+00   \n",
       "50%        3.304030e+05    0.000000e+00        3.501000e+03   \n",
       "75%        7.760675e+05    0.000000e+00        1.363000e+04   \n",
       "max        7.925876e+07    1.680307e+08        3.737314e+06   \n",
       "\n",
       "       Electricity(kBtu)  NaturalGas(kBtu)  TotalGHGEmissions   \n",
       "count       3.663000e+03      3.663000e+03        3663.000000  \\\n",
       "mean        3.529597e+06      1.550744e+06         118.291974   \n",
       "std         1.756966e+07      7.594444e+06         568.240613   \n",
       "min         0.000000e+00      0.000000e+00           0.000000   \n",
       "25%         6.147670e+05      0.000000e+00           5.200000   \n",
       "50%         1.148932e+06      3.574000e+05          29.200000   \n",
       "75%         2.718398e+06      1.386725e+06          92.950000   \n",
       "max         9.191679e+08      3.737314e+08       20445.700000   \n",
       "\n",
       "       GHGEmissionsIntensity  \n",
       "count            3663.000000  \n",
       "mean                1.217882  \n",
       "std                 2.145308  \n",
       "min                 0.000000  \n",
       "25%                 0.100000  \n",
       "50%                 0.600000  \n",
       "75%                 1.400000  \n",
       "max                47.400000  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.describe()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Outliers Count</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>ENERGYSTARScore</th>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SiteEUIWN(kBtu/sf)</th>\n",
       "      <td>322</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SiteEUI(kBtu/sf)</th>\n",
       "      <td>324</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SiteEnergyUse(kBtu)</th>\n",
       "      <td>406</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SiteEnergyUseWN(kBtu)</th>\n",
       "      <td>406</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SourceEUIWN(kBtu/sf)</th>\n",
       "      <td>335</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SourceEUI(kBtu/sf)</th>\n",
       "      <td>336</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                       Outliers Count\n",
       "ENERGYSTARScore                     0\n",
       "SiteEUIWN(kBtu/sf)                322\n",
       "SiteEUI(kBtu/sf)                  324\n",
       "SiteEnergyUse(kBtu)               406\n",
       "SiteEnergyUseWN(kBtu)             406\n",
       "SourceEUIWN(kBtu/sf)              335\n",
       "SourceEUI(kBtu/sf)                336"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "columns_to_check = [\n",
    "    'ENERGYSTARScore', 'SiteEUIWN(kBtu/sf)', 'SiteEUI(kBtu/sf)',\n",
    "    'SiteEnergyUse(kBtu)', 'SiteEnergyUseWN(kBtu)', 'SourceEUIWN(kBtu/sf)',\n",
    "    'SourceEUI(kBtu/sf)']\n",
    "\n",
    "# IQR method to detect outliers\n",
    "outliers_summary = {}\n",
    "\n",
    "for column in columns_to_check:\n",
    "    Q1 = data[column].quantile(0.25)\n",
    "    Q3 = data[column].quantile(0.75)\n",
    "    IQR = Q3 - Q1\n",
    "    lower_bound = Q1 - 1.5 * IQR\n",
    "    upper_bound = Q3 + 1.5 * IQR\n",
    "    \n",
    "    outliers_count = data.loc[(data[column] < lower_bound) | (data[column] > upper_bound), column].count()\n",
    "    outliers_summary[column] = outliers_count\n",
    "\n",
    "outliers_df = pd.DataFrame.from_dict(outliers_summary, orient='index', columns=['Outliers Count'])\n",
    "outliers_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABjUAAAJOCAYAAAD/KYUYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABjYUlEQVR4nO3debyWc+I//tfppE5JqMgoZW9EJaloGDTGNgyyZMo2jJCliZkxYWSJrJF9sg9lyzK2jJlmGGRpGORjjF0lUmGytOh0//6Yn/PtTMU56nR35/l8PHpwv9/Xfd2vc58ep+s6r/t6X2WFQqEQAAAAAACA5Vy9YgcAAAAAAACoCaUGAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQFQxwqFQrEjsJzydwMAAACgdpQaADXw2muvZeDAgfnBD36QzTbbLNtss01++ctf5tVXX6223UEHHZSDDjqo6vGdd96Z8847b4lf/5lnnkm7du2+9s/f//73JX6dYuvZs2d++9vfLna+Xbt2ueyyy5ZZnuuvvz6/+tWvkiR333132rVrl8mTJy92+549e1b7nnz/+99P9+7dc9RRRy30d2Xs2LE56aSTapVn7ty52WWXXfLCCy/U+msBAKD0/fa3v/3ac4If/OAH1bb7y1/+ssj9/O95y2WXXfaN5xtz5sz52m0333zz7Lzzzrn44oszb968hV7zgw8+yIUXXpjdd989nTt3TufOnbP33ntnxIgRmTVrVtV2gwYNSrt27fLEE08sMvvjjz+edu3a5cILL6waq8n5Wk2+xp49e1Z7rRNPPDHt2rXL9ddfX6vvR+fOnbPHHnvkhhtuWOg5NT23BGDx6hc7AMDy7vXXX0/v3r2z+eab59RTT03z5s3zwQcf5JZbbsn++++fP/zhD9l8882TJIMHD6723KuuuirdunVballOO+20bLrppouc22CDDZba65C8+eab+f3vf5/77ruvVs/bbrvt0r9//yTJvHnz8uGHH+b666/PIYcckoceeijNmzdPktx44421ztSgQYP86le/ykknnZQ//vGPqaioqPU+AAAobWussUYuv/zyRc6ttNJK1R4PHjw4W265ZVZbbbUa7fv2229f7FyDBg2+dtuPP/44DzzwQK6++urMmzcvv/71r6vmnnnmmRx//PFZddVV06dPn7Rr1y7z58/PM888k6uuuiqPPPJIRo4cmYYNG2bQoEF54oknctppp+WBBx5I48aNq/bz2Wef5bTTTku7du1y/PHHJ6n5+dp+++2Xbbfdtmpfd955Z0aPHl3t61jwa/z000/zl7/8JRtvvHFuv/32/PznP09ZWdlC78v/fj8KhUKmT5+e2267Leeee24aNmyYPn361CorAF9PqQHwDW644Yasvvrqueaaa1K//v/7sbnjjjtml112yZVXXpkRI0YkSTbccMM6zbLhhhs6yF1GLrjgguy+++5p2bJlrZ7XrFmzhb5HHTp0yI477piHH344ffv2XaJcO+64Yy655JLceuut+fnPf75E+wIAoPQ0aNCgRucEjRo1yn/+85+cddZZueiii2q079qcayxq2x122CGTJ0/O3XffXVVqfPTRRxk4cGDWXXfd3HDDDdVKih/84Af50Y9+lJ/97Ge56aab0q9fvzRt2jRnnHFGjj766Fx88cU55ZRTqra/6KKLMm3atFx11VVVBURNz9fWWmutrLXWWlXzjz/++Nd+zQ888ECS5JRTTskhhxySp59+OltvvfVC2y3u+7H99ttnxx13zN13311VatTm3BKAxbP8FMA3mD59egqFQubPn19tvHHjxjn55JOz6667Vo0teBl3z54989577+Wee+6ptmzRlClTcsIJJ6Rbt27p1KlTDjnkkLzyyitLLe/dd9+d9u3b58UXX0zv3r3ToUOH7LDDDrnuuuuqbTdnzpycf/752W677bLZZptljz32yEMPPVRtm549e+acc87JIYccko4dO1adULz55ps54ogjssUWW6RHjx65+OKLM2jQoKqv/fjjj88Pf/jDhd6zU045JTvvvPNS+1qffPLJ7L///uncuXO6du2ao48+Om+++Wa1bf7yl7+kV69e6dChQ37wgx9kyJAh+eKLL752v6+99loeffTR7L777ovdZubMmdlzzz3Ts2fPTJky5Wv3t+qqq1Z7fNBBB+XZZ5/Ns88+m3bt2uWZZ55Z7PJWi1qS66tL2efOnfu1rwsAwHdXs2bN0q9fvzzwwAMZO3bsMnvdJk2aVLuiYdSoUZkxY0aGDBlSrdD4ylfnRAvO9ezZM3vssUduueWWvPjii0mS5557LrfeemuOP/74fP/736/atjbna7Vx1113Zeutt85WW22Vtm3b5rbbbqvV81daaaU0atSo2ntR26z33ntv9t5773Tq1Cnbb799LrroomrnABMmTMjhhx+e7t27Z4sttshRRx2V119/vWr+q2WMb7vttuywww7ZYost8uSTTyZJ/vGPf+TAAw9Mp06d0q1bt5x00kn56KOPavU1AhSLUgPgG2y//faZMmVKDjjggIwcOTJvvvlm1Q2ed9lll+y9996LfN7ll1+eNdZYI9ttt11uv/32rLnmmvnoo49ywAEH5P/+7//yu9/9LhdddFHmz5+fvn37LvTL+EWZP39+5s2bt9CfysrKhbb75S9/md122y0jRozIFltskfPPP7/q00iFQiHHHHNMbrvttvz85z/PVVddlc6dO2fgwIG59957q+1r5MiR6dChQ6688srsu++++eijj3LggQfm/fffz9ChQ3Pqqafm4YcfrvokU5Lsu+++mTp1ap555pmqsdmzZ+fhhx9e7PtVW5MmTUr//v2z2Wab5aqrrsrZZ5+dt99+O/369as6Sbj//vtzzDHHZP31188VV1yRY489Nvfdd1/69+//tTfpvv/++7PGGmss9lNbn3/+eY444ojMnDkzf/jDH7L22mtXzRUKharvy9y5czNlypScffbZadGiRdVJyuDBg9O+ffu0b98+t99++2KXFFucXXbZJVOnTs2zzz5bq+cBALBiWNQ5wbx58xY6xj366KPTrl27DB48OJ988sm33u///hL+f7edO3dupk6dmmuuuSZPPvlk9txzz6rtxo4dm3bt2mWjjTZa7OuedNJJOfDAA6uNnXrqqWnWrFnOOuuszJ07N6effno233zzHH744dW2+7bna1/n9ddfz4QJE7LXXnslSfbaa6+MHTs206dPX+T2//teTJ48OUOHDs3bb79dtY/aZh05cmROOumkbLrpprn88svTr1+/3HzzzRkyZEiS5Omnn87PfvazJMk555yTIUOG5P33388BBxyw0Lnl5ZdfnpNOOimnnXZaOnfunPHjx+fQQw9NRUVFLrnkkpx88sl59tlnc/DBB2f27Nm1fr8AlrUVavmpuXPnplevXvnd736X7t271+g5t956a6699tp8/PHH6dy5c04//fSss846dZwUKCV9+vTJtGnTct111+XMM89Mkqy++urZZpttcvDBB6djx46LfF779u3ToEGDassR3XTTTfnkk09y6623plWrVkmSH/7wh9ltt90yfPjwXHrppV+b5dBDD13k+EYbbVStVCgUCunfv3/222+/JEmXLl3y5z//OY8++mi23XbbjBs3Lo8//nguvvji7LbbbkmSbbfdNrNmzaq6ed9Xl0OvvfbaVTfLTpLhw4fn888/z7333lu1NFOnTp2qXYGxzTbbZK211sq9995bdYn2n//853zxxRfVDuqXxEsvvZTZs2fnyCOPrMqx1lprZezYsfniiy+y8sor58ILL8y2225b7SaC6667bg499NA89thj2X777Re576effjodOnRY5Jq5c+bMydFHH52pU6fm5ptvTuvWravN33vvvQsVQ2VlZbngggvSrFmzJP9dRqxJkyZJaneJ/1fatm2bVVddNU899VS22WabWj8fAIDS9d577y32QzG/+c1vqv3Sf6WVVsq5556b/fbbL0OGDKl2XLwoi9tv3759c9ppp33jtmuvvXaOO+649OvXr2ps4sSJVTcwX9Cibia+4JJMq622Wk4//fQce+yxOeywwzJ58uTce++9KS8vr/acb3u+9nXuuuuurLbaalU3Dt97771z2WWXZfTo0TnqqKOqbbu478e6666bwYMHVxUPtck6f/78XHHFFdlxxx2rSowkmTVrVh588MF8+eWXueiii9K2bduMGDGi6j3ZZptt8uMf/ziXXnpphg8fXu11d9lll6rHF110UdZbb738/ve/r3pup06d8pOf/CR33XXXEi+ZC1DXVphSY86cOTnxxBOrXWb3TR5//PFccMEFueiii7Luuutm2LBhOeaYY2p9U1hgxTdgwIAceuihefzxx/PUU0/lmWeeyf33358HHnggJ598cg4++OAa7eepp57KJptskpYtW1YdxNerVy8//OEPa/Sz54wzzljkAfOibhjduXPnqv//qlz5atmlp556KmVlZdluu+2qnUz07Nkz9913X15//fVssskmSVL13688/fTT6dy5c7V7TbRq1ara69WrVy977713brrpppx++ulp1KhR7rnnnvTo0aPaOrbfxldFQ6dOndKwYcPsu+++2WWXXfLDH/4w3bt3rzoRePPNN/PBBx/kyCOPrPY1du3aNU2aNMmTTz652FJj0qRJ1b6eBf3mN7/Jyy+/nHPOOWeRJfgOO+yQY445Jsl/y6WPPvooY8aMya9+9avMmjUr+++//5J8+VXWXnvthZaqAgBgxbfGGmvkqquuWuTc9773vYXG2rdvnyOOOCJXXXVVdtttt6pf1C/K6NGjFznevHnzxW47a9as3HjjjXnmmWdy6qmn5kc/+lG17RZ3lceizmv+/e9/V3v84x//OLvttlseeuihnHbaaWnbtu0i8y2t87Uk+fLLL3Pfffdlxx13zOzZszN79uysvPLK6dKlS+64447069cv9er9v4VPFvx+zJw5M1deeWUmTpyYc889d5HnFDXJ+vbbb2fGjBn58Y9/XO25hx9+eA4//PB88cUXmTBhQo499thqJU/Tpk2zww475LHHHqv2vAXP6WbNmpUXX3wxhx9+eNVV5kmyzjrrZIMNNsiTTz6p1ACWeytEqfHGG2/kxBNP/NqlRBblscceyzbbbJMddtghSXLsscfmpz/9aT766KOqT9MCfGXVVVfN7rvvXnWfhVdeeSW//vWvc8EFF2SPPfbI6quv/o37+OSTT/Luu+8u9hNQs2bNSqNGjRb7/PXWWy8dOnSoUd7/LTrq1atX9XPyk08+SaFQyBZbbLHI53744YdVB77/u+7tRx99tMj8LVq0qHY59j777JOrr746jzzySLbaaqs89dRT3/jJsMaNGy/2PhFfjX/1/rRu3Tq33HJLRowYkdGjR+cPf/hDmjZtmj59+uSXv/xl1eX1Z5xxRs4444xFfo2L89lnny32+zB16tRsuummueKKK7LLLrtk5ZVXrja/2mqrLfQ92n777fPhhx/mggsuyD777LPQp8u+jUaNGuWzzz5b4v0AAFBaGjRoUONzgq/0798/Y8eOzWmnnZYuXbosdrva7HfBbbfccssceuihGTBgQG688cZsueWWVXOtWrXKe++9V+259evXr1ag3HHHHbnjjjsW+TrbbrttHnrooWy33XZfm2dpnK8lyaOPPpoZM2Zk9OjRiyx5Hn/88WpZ/vf7scUWW2SfffbJEUcckTvvvDPrrbderbN+dS6zqDIpST799NMUCoW0aNFiobkWLVrk008/rTa24DndzJkzM3/+/FxzzTW55pprFnp+w4YNF/maAMuTFaLUePbZZ9O9e/cMHDhwoWU8/vGPf+Scc87JG2+8kbZt2+bYY4+tWiJltdVWy4MPPpg333wzbdu2zb333ptWrVotdENX4Ltr6tSp2WeffTJgwICqpZy+0r59+wwcODDHHHNMJk2aVKOD5FVWWSXdunXLb37zm0XON2jQYKnkrkmOxo0b5w9/+MMi5xf3Cajkv0s8LWot2RkzZlR7vM4666Rbt24ZM2ZMPvnkkzRp0iQ77rjj1+Zq0aLFYsuGDz74oGqbr3Ts2DGXX3555s6dm+eeey633357rr766nz/+9/PhhtumOS/V1Z069Ztof193c/61VZbbaETga9cfvnladSoUXr16pWLL744p5566td+TV/ZbLPNMm7cuHz88ceLPPn46gqU//0k2+eff77I/c2cObPavTwAAGBxGjRokKFDh6Z37945++yzl/r+69Wrl6FDh+YnP/lJfvvb3+bBBx+s+uV4z549M2LEiEyaNKnalc4LFgGPPvrot3rdpX2+lvx36al11llnofepUCjk2GOPzW233fa1BUujRo1y7rnnpnfv3hk0aFBuvfXWlJWV1Spr06ZNk2ShG3d//PHHeeWVV9K5c+eUlZUt8rxs2rRpWW211Rabb+WVV05ZWVkOPfTQ/OQnP1lkfoDl3Qpxo/A+ffrk5JNPXugH77Rp03LkkUemV69euf/++/OLX/wiv/3tb/OPf/wjSXLQQQdl/fXXz2677ZaOHTvmjjvuyJVXXrlUPkELrBhatGiR+vXrZ9SoUZkzZ85C82+99VYaNmy42BJgwcuSk6Rbt255++23q664+OrPH//4x4wePXqZ/fzp1q1bvvjiixQKhWo5XnvttVxxxRWLXN/2K127ds0LL7yQadOmVY19+OGHeeGFFxbadt999824cePywAMPZLfddvvGT/1069YtL7744iKXVXr44YdTXl6erl27JkluvPHG7LDDDpk7d24aNGiQrbfeOmeddVaSZMqUKVl//fXTvHnzTJ48udrX2LJly1x00UV55ZVXFpujVatWef/99xc516JFi7Rr1y6HHnpoRo4cmRdffPFrv6avTJgwIauuumrVydT//t346h4bX5U3yX+X0FrUDR0LhUKmTp1adV8WAAD4Jptttll+8Ytf5I9//OPXHgt/W61atUr//v0zadKkalcA9O3bN6uttlp++9vfLvJK48rKyrz11lvf6jWX9Hztf02bNi2PP/54fvKTn6R79+7V/my11VbZZZdd8thjj2Xq1Klfu5+OHTtm//33zz//+c+q++3VJuv666+f1VdfPX/729+qbfPHP/4x/fr1y5dffpnNNtssY8aMSWVlZdX8p59+mkcfffRrr8Zp0qRJ2rdvn7feeqvaedJGG22Uyy67LM8880yN3iuAYlohrtRYnJEjR6ZHjx458MADk/z3k8f/+te/ctNNN2XLLbfMhx9+mDlz5uTCCy9M27Ztc+WVV+bXv/51Ro8e7XI7IElSXl6e008/Pcccc0z22Wef9O3bNxtssEFmzZqVJ598MiNHjsyAAQMW+6n/pk2b5pVXXsmzzz6bjh075tBDD80f//jHHHrooTnssMOy+uqr56GHHsodd9yRQYMGfWOeN954Y7E/n9ZYY40a/5J7u+22S9euXdO/f//0798/G2ywQV566aVceuml2Xbbbb92Cb6DDz44I0eOzOGHH15174grr7wyX3755UI31t55551z1lln5aWXXsrvfve7b8x14IEH5s4778xBBx2UI444IhtttFFmz56dcePG5eabb84vfvGLqqsTttpqq1x44YU55phjcuCBB6a8vDy33XZbGjRokB122CHl5eUZOHBgTjvttJSXl2eHHXaoWuP2qyWkFucHP/hBRo0alUKhsMibhSf/XbJwzJgxOfXUU3P33XdnpZVWSvLfT1MtWPDMmjUr9957b5566qmccMIJVcVV06ZN889//jNPPfVU2rdvn+7du6eioiLnnntuBgwYkM8//zyXXnrpIj9l9dprr+XTTz/Ntttu+43vKQAAK5a5c+cu8gNFX2nXrt1i54455piMHTt2sfcj/br9rrfeet+4ssWhhx6a0aNH55prrsnee++dVq1apWXLlrn88sszYMCA/PSnP03v3r2z6aabpl69enn55Zdz11135Z133slPf/rTr933oizp+dr/uvfeezNv3rxFXsGQJHvttVfuvPPO3HHHHTnuuOO+dl+//OUvM2bMmFx00UX58Y9/nCZNmtQq63HHHZczzzwzzZs3T8+ePfP222/n0ksvTd++fbPqqqvmxBNPzOGHH55+/fqlT58++fLLLzNixIjMnTu36jxtcU444YT069cvJ554Yn7605+msrIy119/fV588cX079+/Ru8VQDGt0KXGW2+9lb/97W/Vbsz05ZdfVq1nOHjw4Oy0007ZY489kiQXXXRRtt9++4wdOza77bZbUTIDy5/tt98+d9xxR6677rpcffXV+eijj9KgQYO0b98+F198cXbaaafFPvewww7LOeeck8MPPzw33HBDttxyy9x222256KKLcvrpp2fOnDlZd911c/bZZ2fffff9xixnnnnmYucOPvjgnHLKKTX6murVq5cRI0Zk+PDh+f3vf58ZM2akZcuW+fnPf/6NB8BNmzbNH/7wh5x99tn5zW9+k5VXXjl9+vRJo0aNFrr/RsOGDbPVVlvlrbfeqrqB9zfte/To0bnyyitz/fXX58MPP0xFRUXWX3/9DBkyJHvttVfVtt///vdz9dVX54orrsgJJ5yQysrKbLbZZrn++uuz/vrrJ0n222+/rLzyyrn22mtz++23p3Hjxtliiy1y4YUXLvIm31/ZaaedcsUVV+Sll15Kp06dFrlNo0aNctppp+XII4/MiBEjqt63xx57rNqN+Ro3bpz11lsvgwcPTp8+farG+/btm5dffjlHHHFEhg4dmj322COXXXZZLrroohxzzDFp1apVjj322KpPdi3o73//e9ZYY43F3hMFAIAV17Rp09K7d+/Fzi/q+PErCy5DtShft98rrrjiG5eTbdCgQU4++eQceeSROe+883LppZcm+e89N+6///7ceuutefjhh3PNNddk7ty5+d73vpetttoqF198cdq3b/+1+16cJTlf+1933313Ntpoo2y88caLnO/SpUtat26dO++88xt/+b/66qtnwIABOfPMM3PFFVfkpJNOqlXWvn37pnHjxrnuuuty++23Z6211soRRxyRI444Ikmy9dZb54Ybbsill16aE044IQ0aNMiWW26Z8847LxtttNHXZttmm21y3XXX5fLLL8/xxx+flVZaKZtuumluuOGGhZZ1B1gelRVqe3ft5Vy7du3yhz/8Id27d0///v2z6qqr5qijjqq2Tf369dOqVat06tQpl112WX74wx9Wze27777Zeeedq/6RAKC6F198MZ988km1dWTnzZuX7bffPj/5yU+qXXEye/bsbLfddunfv38OOeSQYsT91o466qisvvrqGTp0aLGjVFMoFLLzzjunT58+OfTQQ4sdBwAAAGCZWiHuqbE46623Xt599920bdu26s/YsWNz//33J0nWXHPNvPnmm1Xbz507N5MnT07r1q2LFRlguTdlypQceeSRVeutPvrooznuuOPy6aefZv/990+SvPfee7n88stz2GGHpaysLPvss0+RU9fewIED88gjj2TKlCnFjlLNI488ksrKyhxwwAHFjgIAAACwzK3QpUafPn3y8ssv5+KLL84777yT+++/P8OGDataj32//fbL1Vdfnb/97W9566238rvf/S4rr7xyevbsWeTkAMuvXXfdNYMHD84jjzySI488Mr/+9a/z5Zdf5pZbbskGG2yQ5L/LW91888354IMPcvHFF1fdBLuUtGvXLkceeWQuvPDCYkepMnfu3AwbNiznn39+Kioqih0HAAAAYJlboZefSpJx48blwgsvzGuvvVa1XvxXNw6vrKysWpvwk08+SefOnTN48OCvXWcdAAAAAAAojhWu1AAAAAAAAFZMK/TyUwAAAAAAwIpDqQEAAAAAAJSE+sUO8G3Nnz8/8+bNS7169VJWVlbsOAAAsMIpFAqZP39+6tevn3r1vtufh3L+AQAAdaum5x8lW2rMmzcvEyZMKHYMAABY4XXo0CENGjQodoyicv4BAADLxjedf5RsqfFVU9OhQ4eUl5cXOQ0AAKx4KisrM2HChO/8VRqJ8w8AAKhrNT3/KNlS46tLvsvLy51UAABAHbLckvMPAABYVr7p/MNHrgAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAYJkYN25cevfunXHjxhU7CgArgLlz52b33XfPM888UzU2adKkHHroodl8882z22675Yknnqj2nHHjxmX33XdPp06dcvDBB2fSpEnLOjYAALCElBoA1LnZs2dn2LBhmTp1aoYNG5bZs2cXOxIAJWzOnDk54YQT8vrrr1eNFQqFHHPMMWnRokXuuuuu7Lnnnjn22GMzZcqUJMmUKVNyzDHHpFevXhk9enSaNWuW/v37p1AoFOvLAAAAvgWlBgB1buTIkZkxY0aSZMaMGRk1alSREwFQqt54443sv//+mThxYrXxp59+OpMmTcqZZ56ZDTbYIEceeWQ233zz3HXXXUmSO++8M5tttlkOO+ywbLTRRhk6dGjee++9PPvss8X4MgAAgG+pfrEDALBimzx5ckaNGlX1SdhCoZBRo0Zlp512SuvWrYucDoBS8+yzz6Z79+4ZOHBgNt9886rxF198Me3bt0/jxo2rxrp06ZIXXnihan7LLbesmmvUqFE23XTTvPDCC+nevfuyig/UQKFQcGUvS8VX5yBlZWVFTkKpq6io8PcIliNKDQDqTKFQyPDhwxc7fv755zswBKBW+vTps8jxadOmZc0116w21rx583zwwQc1mq+pysrKWm0P1E6hUMiAAQPyf//3f8WOAlBls802yyWXXOL8FepYTY+1lRoA1JmJEydm/PjxC41XVlZm/PjxmThxYtq2bVuEZACsaGbNmpUGDRpUG2vQoEHmzp1bo/mamjBhwpIFBb5WoVDI559/XuwYANV89tlneeGFF5QasJxQagBQZ9q0aZOuXbvm+eefr9a2l5eXp0uXLmnTpk0R0wGwImnYsGE++eSTamNz585NRUVF1fz/Fhhz585N06ZNa/U6HTp0SHl5+RJlBb7etddea/kpltjs2bOz7777JklGjx5d9e8BfBuWn4Jlo7KyskYfIlJqAFBnysrKMmDAgBxyyCGLHHdQCMDS0rJly7zxxhvVxqZPn1615FTLli0zffr0heY32WSTWr1OeXm5UgOWgSZNmhQ7AiVuwZ/VK6+8cho1alTENAAsTfWKHQCAFVvr1q3Tp0+fqgKjrKwsffr0SatWrYqcDIAVSadOnfJ///d/1T7d/dxzz6VTp05V888991zV3KxZs/LKK69UzQMAAKVBqQFAnevbt2+aN2+eJGnRosVib/IKAN9Wt27d8r3vfS+DBg3K66+/nhEjRuSll16qWnpkn332yfPPP58RI0bk9ddfz6BBg9K6det07969yMkBAIDaUGoAUOcqKipywgknpGXLlhk4cKD1bAFY6srLy3PllVdm2rRp6dWrV+67775cccUVWXvttZP898rByy67LHfddVf23XfffPLJJ7niiisshQgAACXGPTUAWCZ69OiRHj16FDsGACuQf//739Uet23bNrfccstit99uu+2y3Xbb1XUsAACgDrlSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJRS01/vznP6ddu3bV/hx//PHFjAQAAAAAACyn6hfzxd94443ssMMOOeuss6rGGjZsWMREAAAAAADA8qqopcabb76ZjTfeOGussUYxYwAAAAAAACWgqMtPvfnmm1l33XWLGQEAAAAAACgRRSs1CoVC3n777TzxxBPZeeeds+OOO+bCCy/M3LlzixUJAAAAAABYjhVt+akpU6Zk1qxZadCgQS655JJMnjw5Q4YMyezZs3PqqafWeD+VlZV1mBIAAL67HGsDAADLm6KVGq1atcozzzyTVVddNWVlZdlkk00yf/78/PrXv86gQYNSXl5eo/1MmDChjpMCAAAAAADLg6LeKHy11Var9niDDTbInDlz8p///CfNmjWr0T46dOhQ4wIEAACoucrKSh8iAgAAlitFKzUef/zx/OpXv8qjjz6aRo0aJUn+9a9/ZbXVVqtxoZEk5eXlSg0AAAAAAPgOKNqNwjt37pyGDRvm1FNPzVtvvZXHHnss559/fn7xi18UKxIAAAAAALAcK9qVGk2aNMl1112Xc845J/vss09WXnnlHHDAAUoNAAAAAABgkYp6T42NNtooN9xwQzEjAAAAAAAAJaJoy08BAAAAAADUhlIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSUL/YAQAAACh9hUIhs2fPLnYMgCSp9vPIzyZgeVJRUZGysrJixyhpSg0AAACW2OzZs7PrrrsWOwbAQvbee+9iRwCoMmbMmDRq1KjYMUqa5acAAAAAAICS4EoNAAAAlqrPNv9ZCvWcbgJFVij897+WeQGKrGz+vDR54dZix1hhOMoEAABgqSrUq5+Ur1TsGAAAy4VCsQOsYCw/BQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAMAK4/3338+RRx6ZLbbYIj179syNN95YNffKK69kv/32S6dOnbLPPvvk5ZdfLl5QAADgW1FqAAAAK4xf/vKXady4ce6+++6cfPLJueSSS/LnP/85X3zxRfr165ctt9wyd999dzp37pwjjzwyX3zxRbEjAwAAtaDUAAAAVgj/+c9/8sILL+Too4/Ouuuumx133DHbbrttnnrqqTz00ENp2LBhfvOb32SDDTbIKaeckpVXXjkPP/xwsWMDAAC1oNQAAABWCBUVFWnUqFHuvvvufPnll3nrrbfy/PPPZ5NNNsmLL76YLl26pKysLElSVlaWLbbYIi+88EJxQwMAALVSv9gBAAAAloaGDRvmtNNOy1lnnZU//OEPqaysTK9evbLffvtl7Nix2XDDDatt37x587z++uu1eo3KysqlGXmF4r0BAPhmlZWVjpsWo6bvi1IDAABYYbz55pvZYYcd8vOf/zyvv/56zjrrrGy99daZNWtWGjRoUG3bBg0aZO7cubXa/4QJE5Zm3BXKnDlzih0BAGC599JLL6Vhw4bFjlHSlBoAAMAK4amnnsro0aPz2GOPpaKiIh06dMjUqVNz1VVXZZ111lmowJg7d24qKipq9RodOnRIeXn50oy9wpg1a1axIwAALPc6duyYRo0aFTvGcqmysrJGHyJSagAAACuEl19+OW3btq1WVLRv3z5XX311ttxyy0yfPr3a9tOnT8+aa65Zq9coLy9XaiyG9wUA4Js5nlxybhQOAACsENZcc828++671a7IeOutt9K6det06tQp//znP1MoFJIkhUIhzz//fDp16lSsuAAAwLeg1AAAAFYIPXv2zEorrZRTTz01b7/9dv7617/m6quvzkEHHZRddtklM2fOzNlnn5033ngjZ599dmbNmpVdd9212LEBAIBaUGoAAAArhFVWWSU33nhjpk2bln333TdDhw7N0Ucfnd69e6dJkyb5/e9/n+eeey69evXKiy++mBEjRqRx48bFjg0AANSCe2oAAAArjA033DA33HDDIuc6duyYe+65ZxknAgAAliZXagAAAAAAACVBqQEAAAAAAJSE5abU6NevX377298WOwYAAAAAALCcWi5KjQcffDCPPfZYsWMAAAAAAADLsaKXGp988knOP//8dOjQodhRAAAAAACA5Vj9Ygc477zzsueee+bDDz8sdhQAAAAAAGA5VtRS46mnnso//vGP3H///Tn99NO/1T4qKyuXbigAACCJY20AAGD5U7RSY86cORk8eHBOO+20VFRUfOv9TJgwYSmmAgAAAAAAlldFKzUuv/zybLbZZtl2222XaD8dOnRIeXn5UkoFAAB8pbKy0oeIAACA5UrRSo0HH3ww06dPT+fOnZMkc+fOTZL86U9/yj//+c8a76e8vFypAQAAAAAA3wFFKzVuvvnmzJs3r+rxhRdemCT51a9+VaxIAAAAAADAcqxopUarVq2qPV555ZWTJG3bti1GHAAAAAAAYDlXr9gBAAAAAAAAaqJoV2r8r3PPPbfYEQAAAAAAgOWYKzUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASkL9YgcAAABgBVP5ZbETAAAsPxwbLVVKDQAAAJZYoVCo+v9VXrytiEkAAJZfCx4z8e1YfgoAAAAAACgJrtQAAABgiZWVlVX9/6edDkjKVypiGgCA5Ujll1VXsi54zMS3o9QAAABg6SpfSakBAECdsPwUAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaACwT48aNS+/evTNu3LhiRwEAAACgRCk1AKhzs2fPzrBhwzJ16tQMGzYss2fPLnYkAAAAAEqQUgOAOjdy5MjMmDEjSTJjxoyMGjWqyIkAAAAAKEVKDQDq1OTJkzNq1KgUCoUkSaFQyKhRozJ58uQiJwMAAACg1Cg1AKgzhUIhw4cPX+z4V0UHAAAAANSEUgOAOjNx4sSMHz8+lZWV1cYrKyszfvz4TJw4sUjJAAAAAChFSg0A6kybNm3StWvXlJeXVxsvLy9Pt27d0qZNmyIlAwAAAKAUKTUAqDNlZWUZMGDAYsfLysqKkAoAAACAUqXUAKBOtW7dOn369KkqMMrKytKnT5+0atWqyMkAAAAAKDVKDQDqXN++fdO8efMkSYsWLdKnT58iJwIAAACgFCk1AKhzFRUVOeGEE9KyZcsMHDgwFRUVxY4EAAAAQAmqX+wAAHw39OjRIz169Ch2DAAAAABKmCs1AAAAAACAkqDUAAAAAAAASkKtlp+aPXt27r///jz++OP5v//7v3z00UcpKyvLGmuskfbt2+eHP/xhdtlllzRq1Kiu8gIAACzW3LlzM3To0DzwwANZaaWVsu+++2bgwIEpKyvLK6+8ksGDB+e1117LhhtumDPOOCObbbZZsSMDAAC1UKNSY+7cuRkxYkT+8Ic/ZN11102PHj2y0047ZbXVVsv8+fPz8ccf59///nduv/32nHvuuenTp0+OOuqoNGzYsK7zAwAAVBkyZEieeeaZXHfddfn8888zcODArL322vnpT3+afv36ZY899si5556bW2+9NUceeWT+/Oc/p3HjxsWODQAA1FCNSo0DDjggPXv2zEMPPZQWLVp87bbvvfde7rjjjvTu3Tv33nvv0sgIAADwjT755JPcddddueGGG9KxY8ckyWGHHZYXX3wx9evXT8OGDfOb3/wmZWVlOeWUU/L3v/89Dz/8cHr16lXk5AAAQE3VqNS4/vrrs9pqq9Voh61atcrAgQPz85//fElyAQAA1Mpzzz2XJk2apFu3blVj/fr1S5L87ne/S5cuXVJWVpYkKSsryxZbbJEXXnhBqQEAACWkRqVGTQuNJX0OAADAtzVp0qS0atUq9957b66++up8+eWX6dWrV44++uhMmzYtG264YbXtmzdvntdff71Wr1FZWbk0I69QvDcAAN+ssrLScdNi1PR9qdWNwpOkZ8+eVZ9uWpSxY8fWdpcAAABL7Isvvsi7776b2267LUOHDs20adNy2mmnpVGjRpk1a1YaNGhQbfsGDRpk7ty5tXqNCRMmLM3IK5Q5c+YUOwIAwHLvpZdeci/qJVTrUuO4446r9njevHmZNGlS7r777gwYMGCpBQMAAKiN+vXr57PPPstFF12UVq1aJUmmTJmSW2+9NW3btl2owJg7d24qKipq9RodOnRIeXn5Usu8Ipk1a1axIwAALPc6duyYRo0aFTvGcqmysrJGHyKqdamx9957L3K8U6dOuf7667PffvvVdpcAAABLbI011kjDhg2rCo0kWW+99fL++++nW7dumT59erXtp0+fnjXXXLNWr1FeXq7UWAzvCwDAN3M8ueTqLa0dbbjhhi7FBgAAiqZTp06ZM2dO3n777aqxt956K61atUqnTp3yz3/+M4VCIUlSKBTy/PPPp1OnTsWKCwAAfAu1vlJj/PjxC419/vnnufnmm7PRRhstlVAAAAC1tf7662f77bfPoEGDcvrpp2fatGkZMWJEjj766Oyyyy656KKLcvbZZ+eAAw7IbbfdllmzZmXXXXctdmwAAKAWal1qHHTQQQuNrbTSSunQoUOGDBlSq329++67OfPMM/P8889n1VVXzYEHHphf/OIXtY0EAACQJLnwwgtz1lln5Wc/+1kaNWqUvn375qCDDkpZWVl+//vfZ/DgwbnjjjvSrl27jBgxIo0bNy52ZAAAoBZqXWq8+uqrS+WF58+fn379+qVDhw6555578u677+aEE05Iy5Yts8ceeyyV1wAAAL5bVllllZx//vmLnOvYsWPuueeeZZwIAABYmmp9T40f/ehH+eSTTxYanzp1arbeeusa72f69OnZZJNNcvrpp2fdddfNdtttl6233jrPPfdcbSMBAAAAAADfATW6UuPhhx/OY489liR57733cuaZZ6Zhw4bVtnnvvfdqddf2NddcM5dcckmS/3eTvvHjx2fw4ME13gcAAAAAAPDdUaMrNbp161btcaFQWGibjTbaKFdeeeW3CtGzZ8/06dMnnTt3zs477/yt9gEAAAAAAKzYanSlRrNmzTJ06NAkSatWrXLYYYct1RvqXXrppZk+fXpOP/30DB06NKeeemqNn1tZWbnUcgAAAP+PY20AAGB5U+sbhR977LFV///RRx9lzJgxKRQK6dmzZ9Zee+1vFaJDhw5Jkjlz5uRXv/pVfvOb36RBgwY1eu6ECRO+1WsCAAAAAAClpUalxqxZs3L++efnoYceSpLsueeeOeigg3LAAQdk1qxZKRQKueCCC3Lttdema9euNXrh6dOn54UXXsiOO+5YNbbhhhvmyy+/zGeffZZmzZrVaD8dOnSo1b08AACAmqmsrPQhIgAAYLlSo1LjnHPOyUsvvZTBgwenoqIit9xyS/bff/9ss802Oeecc1JWVpYzzzwzl156aW6++eYavfDkyZNz7LHH5rHHHkvLli2TJC+//HKaNWtW40IjScrLy5UaAAAAAADwHVCjUmPs2LG5+uqr07FjxyTJ5ptvnh49euTAAw/MSiutlCQ57LDDsvfee9f4hTt06JBNN900J598cgYNGpT33nsvF1xwQY466qhv8WUAAAAAAAAruno12eijjz7KWmutVfW4WbNmadSoUVZfffWqsSZNmmT27Nk1fuHy8vJceeWVadSoUXr37p1TTjklBx10UA4++OBaxAcAAAAAAL4ranyj8EUt8VRWVrZEL96yZctcfvnlS7QPAAAAAADgu6HGpcY///nPrLrqqlWPC4VCXnrppXzwwQdJkv/85z9LPx0AAAAAAMD/r8alxrHHHrvQ2IknnrhUwwAAAAAAACxOjUqNV199ta5zAAAAAAAAfK0a3Sh8QU8++eQix6dMmZL+/fsvcSAAAAAAAIBFqXWp0b9///zpT3+qevzll1/myiuvzG677ZYZM2Ys1XAAAAAAAABfqfE9Nb5y0UUX5de//nVmzpyZtdZaK2eddVa++OKLDB48OHvvvXddZAQAAAAAAKh9qbHjjjvmuuuuy9FHH51PP/00P//5z3P00UenSZMmdZEPAAAAAAAgSQ1LjfHjxy80NnDgwJx99tmZN29eXn311RQKhSRJ165dl25CAAAAAACA1LDUOOiggxY7d9NNN+Wmm25KkpSVleVf//rX0kkGAAAAAACwgBqVGq+++mpd5wAAAAAAAPha9Wqy0ejRo6uWl6qJysrK3Hnnnd86FAAAAAAAwP+qUakxadKk7L777rnmmmvy9ttvL3a7d999N1deeWV22223TJw4camFBAAAAAAAqNHyUwMHDsyee+6Za6+9NnvvvXdWX331rL/++ll99dUzf/78fPLJJ3nttdcyc+bM/OQnP8mVV16ZDTbYoK6zAwAAAAAA3yE1KjWSZP31188555yTQYMG5dlnn80rr7ySjz76KGVlZdlggw1y0EEHpXv37mncuHFd5gUAAAAAAL6jalxqfGWVVVbJj370o/zoRz+qizwAAAAAAACLVKN7agAAAAAAABSbUgMAAAAAACgJSg0Alolx48ald+/eGTduXLGjAFAEn376aUaOHJkhQ4bko48+yt/+9rdMnDix2LEAAIASo9QAoM7Nnj07w4YNy9SpUzNs2LDMnj272JEAWIZee+217LTTTrnrrrty22235fPPP88jjzySPffcM88++2yx4wEAACXkW5Ua9913X3r16pUtt9wykyZNytlnn50RI0Ys7WwArCBGjhyZGTNmJElmzJiRUaNGFTkRAMvSkCFD8rOf/Sx33313VlpppSTJ0KFD06dPn5x//vlFTgcAAJSSWpcao0aNyvnnn59evXrlyy+/TJJsttlmue6663L55Zcv9YAAlLbJkydn1KhRKRQKSZJCoZBRo0Zl8uTJRU4GwLIyYcKE7LXXXguNH3DAAXnjjTeWfSAAAKBk1brUuPnmmzNkyJAceOCBqVfvv0/fc889c/755+fOO+9c6gEBKF2FQiHDhw9f7PhXRQcAK7ZmzZrl7bffXmj8+eefT/PmzYuQCAAAKFX1a/uEKVOmZIMNNlhofJ111sknn3yyNDIBsIKYOHFixo8fv9B4ZWVlxo8fn4kTJ6Zt27ZFSAbAsnTEEUfk1FNPzVFHHZVCoZCnn34699xzT2666aYMHDiw2PEAAIASUusrNTp16pR777232lihUMj111+fjh07Lq1cAKwA2rRpk65du6a8vLzaeHl5ebp165Y2bdoUKRkAy9IBBxyQM844Iw899FAqKipy/vnn54knnshZZ52VAw88sNjxAACAElLrKzVOPfXU9OvXL48++mjmzp2bM844I++8805mzZqVa6+9ti4yAlCiysrKMmDAgBxyyCGLHC8rKytSMgCWpWuvvTa77757Ro4cWewoAABAiav1lRobb7xx/vSnP+VnP/tZDj744Ky//vo5/PDD88gjj2STTTapi4wAlLDWrVunT58+VQVGWVlZ+vTpk1atWhU5GQDLytVXX50vv/yy2DEAAIAVQK1LjSRp2LBhevXqlZNOOim/+MUvstpqq+XDDz9c2tkAWEH07du36kawLVq0SJ8+fYqcCIBlaffdd89VV12Vd955J3Pnzi12HAAAoITVevmp5557Lr/85S9zwQUXZP3110+vXr0yZ86czJo1KxdccEF23XXXusgJQAmrqKjICSeckOHDh2fAgAGpqKgodiQAlqG///3vmTJlSu65555Fzv/rX/9axokAAIBSVetSY+jQodltt93SqVOnXHfddWnYsGH++te/5sEHH8yll16q1ABgkXr06JEePXoUOwYARXDuuecWOwIAALCCqPXyU6+99loOOeSQNGrUKH/961+z0047pUGDBunWrVumTJlSFxkBAIAS1q1bt3Tr1i1rrrlmPv300/znP/9Js2bNqsYBAABqqtZXarRo0SJvvPFGvvjii7zyyiv57W9/myQZN25cvve97y31gAAAQGmbOXNmBg0alLFjx2bVVVdNZWVlPv/883Tt2jVXXHFFVllllWJHBAAASkStr9Q49NBDc8wxx2SfffZJhw4d0q1bt1x99dU544wzcswxx9RFRgAAoIQNGTIkH3zwQR566KE888wz+cc//pH7778/X3zxRYYOHVrseAAAQAmp9ZUaBx98cLbccstMmTIl22yzTZJkq622yvbbb5/vf//7Sz0gAABQ2v7617/mhhtuyPrrr181tuGGG+a0007LEUccUcRkAABAqal1qZEk7du3T/v27aseb7755ksrDwAAsIJp2LBh6tVb+CLxsrKyVFZWFiERAABQqmpUavTs2TNlZWULP7l+/TRt2jSbbLJJDjzwwGy88cZLPSAAAFDaevbsmTPOOCMXXnhh2rRpkyR55513MmTIkGy33XZFTgcAAJSSGpUaxx133CLH58+fn08//TQvvvhievfunWuvvTZdunRZqgEBAIDS9utf/zrHHHNMdtppp6y66qpJkv/85z/54Q9/mN/97ndFTgcAAJSSGpUae++99zduc/nll+eSSy7JzTffvMShAACAFUfTpk1z880359///nfefPPNNGzYMOutt161e2wAAADUxLe6p8ai7LTTTrnxxhuX1u4AAIAVxNy5c3PJJZekVatW6du3b5KkV69e6dGjRwYMGJCVVlqpyAkBAIBSsfDd+r6lioqKFAqFpbU7AABgBTFkyJA89thj+f73v1811r9//zz66KM577zzipgMAAAoNUut1BgzZky1kxQAAIAkeeSRR3LhhRdWu//ejjvumKFDh+ahhx4qYjIAAKDU1Gj5qXvvvXeR41/dKPyFF17IX/7yl4wYMWJpZgMAAFYAhUIhc+bMWeT4l19+WYREAABAqapRqXHppZcucnyllVbKKqusknbt2uXWW2/NZptttlTDAQAApW/nnXfO7373uwwePDjt27dPkrz66qsZMmRIfvzjHxc5HQAAUEpqVGr89a9/rescAADACmrQoEE55ZRTcsghh2T+/PlJkvLy8uy55545+eSTi5wOAAAoJTUqNQAAAL6N6dOnZ/XVV8+wYcMyc+bMvPPOOxk/fnwaNmyYXr16pXHjxsWOCAAAlJCldqNwAACAr3z++ec56qijsu222+add95JkowdOzYHHHBARo4cmZEjR2aPPfbIBx98UNygAABASVFqAAAAS91ll12W9957L7fcckvWX3/9fPHFFxkyZEg6duyYP/3pTxkzZky22WabXHjhhcWOCgAAlBClBgDLxLhx49K7d++MGzeu2FEAWAYeeeSRnHLKKenSpUvKysryxBNP5PPPP89BBx2UlVZaKUnSq1evPPHEE0VOCgAAlJJalxoHHnhgbr311nz00Ud1kQeAFdDs2bMzbNiwTJ06NcOGDcvs2bOLHQmAOjZt2rS0adOm6vG4ceNSXl6ebbbZpmqsRYsWmTVrVjHiAQAAJarWpUbPnj1zzz335Ic//GEOP/zw3HXXXfn000/rIhsAK4iRI0dmxowZSZIZM2Zk1KhRRU4EQF1r2bJlJk2alCQpFAp57LHH0qlTp6y66qpV2/zzn//M9773vWJFBAAASlCtS43DDjssd9xxR/70pz9l6623zu23355tttkm/fv3z4MPPuiTVgBUM3ny5IwaNSqFQiHJf3+xNWrUqEyePLnIyQCoS3vuuWfOPvvsjB07Nuecc07ef//99OnTp2r+1VdfzbBhw7LLLrsUMSUAAFBq6n/bJ7Zq1Sq/+MUvstNOO+Wuu+7KjTfemL///e9ZaaWVsscee+SXv/xlmjVrtjSzAlBiCoVChg8fvtjx888/P2VlZUVIBkBdO/roo/PZZ5/l5JNPTllZWY4//vjsvvvuSZLzzjsvN9xwQ7bffvscffTRRU4KAACUkm9Varz77rt5+OGH8/DDD+e1115Lt27dcuqpp2bnnXfORx99lDPPPDP9+vXL6NGjl3ZeAErIxIkTM378+IXGKysrM378+EycODFt27YtQjIA6lr9+vUzaNCgDBo0aKG5vfbaK3vssUfat29fhGQAAEApq3Wp8dOf/jSvv/56OnTokL322iu77bZb1lhjjar5pk2bZv/998/vfve7pRoUgNLTpk2bdO3aNc8//3wqKyurxsvLy9OlS5dqN5AF4LujXbt2xY4AAACUqFqXGrvttlt+8pOfZJ111lnsNttuu20effTRJckFwAqgrKwsAwYMyCGHHLLIcUtPAQAAAFAbtS41unTpkg8++CAffPBBtfGysrKstNJKWWONNbL22msvtYAAlLbWrVunT58+ueWWW1IoFFJWVpY+ffqkVatWxY4GAAAAQImpdalxyimnZPLkyZk/f35WXXXVFAqFzJw5M2VlZSkrK0uhUEjHjh1z2WWXZc0116yLzACUmL59+2bMmDGZPn16WrRokT59+hQ7EgAAAAAlqF5tn7D33nunQ4cOGTNmTJ555pk8++yz+fOf/5wtt9wyv/71r/Pkk0+mZcuWGTJkSF3kBaAEVVRUZNddd029evWyyy67pKKiotiRAAAAAChBtS41brrpppxxxhlZb731qsbWWWednHLKKfn973+fZs2aZcCAAXnqqaeWalAAStfs2bMzZsyYzJ8/P2PGjMns2bOLHQkAAACAElTrUiNJPv7440WOVVZWVj1281cAvjJy5MjMmDEjSTJjxoyMGjWqyIkAAAAAKEW1LjX23XffnHTSSbnnnnvy+uuv57XXXss999yTQYMGZe+9987HH3+cCy64IN26dauLvACUmMmTJ2fUqFEpFApJkkKhkFGjRmXy5MlFTgYAAABAqan1jcJPPPHErLzyyrn44ovz4YcfJknWXHPNHHjggTn88MMzbty41K9fP6eddtpSDwtAaSkUChk+fPhix88//3xX9gEAAABQY7UuNR588MH87Gc/y9FHH52PP/449evXzyqrrFI1v+2222bbbbddqiEBKE0TJ07M+PHjFxqvrKzM+PHjM3HixLRt27YIyQAAAAAoRbVefuqMM87IRx99lCRZffXVqxUaALCgNm3apGvXrikvL682Xl5enm7duqVNmzZFSgYAAABAKap1qdG9e/c88MADmTt3bl3kAWAFUlZWlgEDBix23NJTAAAAANRGrZefmjFjRq688spcffXVadasWRo2bFhtfuzYsUstHAClr3Xr1unTp09uueWWFAqFlJWVpU+fPmnVqlWxowEAAABQYmpdauy///7Zf//96yILACuovn37ZsyYMZk+fXpatGiRPn36FDsSAAAAACWo1qXG3nvvXfX///nPf7LKKqukrKzMEiIALFZFRUVOOOGEDB8+PAMGDEhFRUWxIwEAAABQgmpdahQKhVx99dW58cYb8+mnn+ZPf/pThg8fnsaNG+fUU09NgwYN6iInACWuR48e6dGjR7FjAAAAAFDCan2j8CuuuCL33Xdfzj333KoCY++9986TTz6Z888/f6kHBAAAAAAASL5FqXHPPffkzDPPzA477FC15NQPfvCDnHfeeRkzZsxSDwgAAAAAAJB8i1JjxowZWXPNNRcab9q0ab744oulEgoAAAAAAOB/1brU2GqrrXLddddVG/vss88ybNiwdO/efakFAwAAAAAAWFCtS43TTz89r7zySn7wgx9kzpw56d+/f7bbbru89957OfXUU+siIwAAAAAAQOrX9glrrbVWRo8enaeeeipvvfVW5s2bl/XWWy/bbLNN6tWrdUcCAAAAAABQI7UuNb6y9dZbZ+utt16aWQAAAAAAABar1qXGK6+8kiFDhmTChAmZN2/eQvP/+te/lkowAAAAAACABdW61Dj55JOzyiqrZPjw4WnSpEldZAIAAAAAAFhIrUuNt956K/fff3/atm1bF3kAAAAAAAAWqdZ39t5kk03y5ptv1kUWAACApaJfv3757W9/W/X4lVdeyX777ZdOnTpln332ycsvv1zEdAAAwLdV6ys19txzz5x66qnp1atX2rZtm5VWWqna/F577bW0sgEAANTagw8+mMceeyx77713kuSLL75Iv379sscee+Tcc8/NrbfemiOPPDJ//vOf07hx4yKnBQAAaqPWpca1116bioqKPPTQQwvNlZWVKTUAAICi+eSTT3L++eenQ4cOVWMPPfRQGjZsmN/85jcpKyvLKaeckr///e95+OGH06tXryKmBQAAaqvWpcZf//rXusgBAACwxM4777zsueee+fDDD6vGXnzxxXTp0iVlZWVJ/vthrC222CIvvPCCUgMAAEpMjUqN8ePHp3Pnzqlff/Gbz5o1KzfddFOOOuqopRYOAACgpp566qn84x//yP3335/TTz+9anzatGnZcMMNq23bvHnzvP7667V+jcrKyiWNucLy3gAAfLPKykrHTYtR0/elRqXGwQcfnCeeeCLNmzevGttjjz0yYsSIfO9730uSfP755xk+fLhSAwAAWObmzJmTwYMH57TTTktFRUW1uVmzZqVBgwbVxho0aJC5c+fW+nUmTJiwRDlXZHPmzCl2BACA5d5LL72Uhg0bFjtGSatRqVEoFBYamzx5cubNm7fUAwEAANTW5Zdfns022yzbbrvtQnMNGzZcqMCYO3fuQuVHTXTo0CHl5eXfOueKbNasWcWOAACw3OvYsWMaNWpU7BjLpcrKyhp9iKjW99RYmqZOnZqzzz47Tz/9dBo2bJjddtstJ5xwgqYKAAColQcffDDTp09P586dk6SqxPjTn/6U3XffPdOnT6+2/fTp07PmmmvW+nXKy8uVGovhfQEA+GaOJ5dc0UqNQqGQ448/Pk2bNs3IkSPzn//8JyeffHLq1auXk046qVixAACAEnTzzTdXu5L8wgsvTJL86le/yvjx43PNNdekUCikrKwshUIhzz//vKVzAQCgBBWt1Hjrrbfywgsv5Mknn0yLFi2SJMcff3zOO+88pQYAAFArrVq1qvZ45ZVXTpK0bds2zZs3z0UXXZSzzz47BxxwQG677bbMmjUru+66azGiAgAAS6DGpcaYMWPSpEmTqsfz58/Pn//85zRr1ixJ8umnn9bqhddYY41ce+21VYXGVz777LNa7QcAAODrNGnSJL///e8zePDg3HHHHWnXrl1GjBiRxo0bFzsaAABQSzUqNdZee+1cf/311caaN2+eW265pdrY9773vRq/cNOmTavdxG/+/Pm55ZZbstVWW9V4H8l/bx4CAAAsfaV8rH3uuedWe9yxY8fcc889RUoDAAAsLTUqNf7617/WdY5ccMEFeeWVVzJ69OhaPa8md0MHAAAAAABKX9HuqbGgCy64IDfddFMuvvjibLzxxrV6bocOHdwtHgAA6kBlZaUPEQEAAMuVopcaZ511Vm699dZccMEF2XnnnWv9/PLycqUGAAAAAAB8BxS11Lj88stz2223ZdiwYdlll12KGQUAAAAAAFjO1SvWC7/55pu58sorc8QRR6RLly6ZNm1a1R8AVjzjxo1L7969M27cuGJHAQAAAKBEFe1KjbFjx6aysjJXXXVVrrrqqmpz//73v4uUCoC6MHv27AwbNizTp0/PsGHDssUWW6SioqLYsQAAAAAoMUUrNfr165d+/foV6+UBWIZGjhyZGTNmJElmzJiRUaNG5bDDDityKgAAAABKTdGWnwLgu2Hy5MkZNWpUCoVCkqRQKGTUqFGZPHlykZMBAAAAUGqUGgDUmUKhkOHDhy92/KuiAwAAAABqQqkBQJ2ZOHFixo8fn8rKymrjlZWVGT9+fCZOnFikZAAAAACUIqUGAHWmTZs26dq1a8rLy6uNl5eXp1u3bmnTpk2RkgEAAABQipQaANSZsrKyDBgwYLHjZWVlRUgFAAAAQKlSagBQp1q3bp0+ffpUFRhlZWXp06dPWrVqVeRkAAAAAJQapQYAda5v375p3rx5kqRFixbp06dPkRMBAAAAUIqUGgDUuYqKipxwwglp2bJlBg4cmIqKimJHAgAAAKAE1S92AAC+G3r06JEePXoUOwYAAAAAJcyVGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaAAAAAABASVBqAAAAAAAAJUGpAQAAAAAAlASlBgAAAAAAUBKUGgAAAAAAQElQagAAAAAAACVBqQEAAAAAAJQEpQYAAAAAAFASlBoAAAAAAEBJUGoAAAAAAAAlQakBAAAAAACUBKUGAAAAAABQEpQaAAAAAABASVBqALBMjBs3Lr179864ceOKHQUAAACAEqXUAKDOzZ49O8OGDcvUqVMzbNiwzJ49u9iRAAAAAChBSg0A6tzIkSMzY8aMJMmMGTMyatSoIicCAAAAoBQpNQCoU5MnT86oUaNSKBSSJIVCIaNGjcrkyZOLnAwAAACAUqPUAKDOFAqFDB8+PPPnz682XllZmeHDh1cVHQAAAABQE0oNAOrMxIkTM378+IXKi0KhkPHjx2fixIlFSgYAAABAKVJqAFBn1llnnTRt2nSRc02bNs0666yzjBMBAAAAUMqUGgDUmUmTJmXmzJmLnJs5c2YmTZq0jBMBAAAAUMqUGgDUmTZt2qRr164pKyurNl5WVpZu3bqlTZs2RUoGAAAAQClSagBQZ8rKyjJgwICFSo169eotchwAAAAAvo5SA4A61bp16/Tu3bvaWO/evdOqVasiJQIAAACgVCk1AAAAAACAkqDUAKBOTZ48OXfccUe1sTvuuCOTJ08uUiIAAAAASpVSA4A6UygUMnz48MWOFwqFIqQCAAAAoFQpNQCoMxMnTsz48eNTWVlZbbyysjLjx4/PxIkTi5QMAAAAgFKk1ACgzrRp0yZdu3ZNeXl5tfHy8vJ069Ytbdq0KVIyAAAAAEqRUgOAOlNWVpYBAwYsdrysrKwIqQAAAAAoVUoNAOpU69at06dPn6oCo6ysLH369EmrVq2KnAwAAACAUqPUAKDO9e3bN82bN0+StGjRIn369ClyIgAAAABKkVIDgDpXUVGRE044IS1btszAgQNTUVFR7EgAAAAAlKD6xQ4AwHdDjx490qNHj2LHAAAAAKCEuVIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAWCbGjRuX3r17Z9y4ccWOAgAAAECJUmoAUOdmz56dYcOGZerUqRk2bFhmz55d7EgAAAAAlCClBgB1buTIkZkxY0aSZMaMGRk1alSREwEAAABQipQaANSpyZMnZ9SoUSkUCkmSQqGQUaNGZfLkyUVOBgAAAECpUWoAUGcKhUKGDx+e+fPnVxuvrKzM8OHDq4oOAAAAAKgJpQYAdWbixIkZP378QuVFoVDI+PHjM3HixCIlAwAAAKAUKTUAqDPrrLNOmjZtusi5pk2bZp111lnGiQAAAAAoZUoNAOrMpEmTMnPmzEXOzZw5M5MmTVrGiQAAAAAoZUoNAOpMmzZt0rVr15SVlVUbLysrS7du3dKmTZsiJQNgRTR16tQcf/zx6datW7bddtsMHTo0c+bMSfLfov3QQw/N5ptvnt122y1PPPFEkdMCAADfhlIDgDpTVlaWAQMGLFRq1KtXb5HjAPBtFQqFHH/88Zk1a1ZGjhyZiy++OH/7299yySWXpFAo5JhjjkmLFi1y1113Zc8998yxxx6bKVOmFDs2AABQS/WLHQCAFVvr1q3Tt2/f3HzzzVVjffv2TatWrYqYCoAVzVtvvZUXXnghTz75ZFq0aJEkOf7443Peeeflhz/8YSZNmpTbbrstjRs3zgYbbJCnnnoqd911V4477rgiJwcAAGpjubhSY+7cudl9993zzDPPFDsKAHWgb9++Vb9gWmONNdKnT58iJwJgRbPGGmvk2muvrfr35iufffZZXnzxxbRv3z6NGzeuGu/SpUteeOGFZZwSAABYUkW/UmPOnDk58cQT8/rrrxc7CgB1pKKiIieccEKGDx+eAQMGpKKiotiRAFjBNG3aNNtuu23V4/nz5+eWW27JVlttlWnTpmXNNdestn3z5s3zwQcf1Pp1Kisrlzjrisp7AwDwzSorKx03LUZN35eilhpvvPFGTjzxxBQKhWLGAGAZ6NGjR3r06FHsGAB8R1xwwQV55ZVXMnr06Nx4441p0KBBtfkGDRpk7ty5td7vhAkTllbEFc5XN2UHAGDxXnrppTRs2LDYMUpaUUuNZ599Nt27d8/AgQOz+eabFzMKAACwgrjgggty00035eKLL87GG2+chg0b5pNPPqm2zdy5c7/VlYMdOnRIeXn5Ukq6Ypk1a1axIwAALPc6duyYRo0aFTvGcqmysrJGHyIqaqlhTXUAAGBpOuuss3LrrbfmggsuyM4775wkadmyZd54441q202fPn2hJalqory8XKmxGN4XAIBv5nhyyRX9nhpLyvpjAKXh+uuvz6233pqf/exnOeyww4odB4AaKLVj7csvvzy33XZbhg0bll122aVqvFOnThkxYkRmz55ddXXGc889ly5duhQrKgAA8C2VfKlhTVuA5d9nn32WUaNGpVAoZNSoUdl4443TpEmTYscCYAXy5ptv5sorr0y/fv3SpUuXTJs2rWquW7du+d73vpdBgwalf//++dvf/paXXnopQ4cOLWJiAADg2yj5UsOatgDLvwEDBqRQKCRJCoVCRo8enUsuuaS4oQD4RjVd03Z5MHbs2FRWVuaqq67KVVddVW3u3//+d6688sqccsop6dWrV9q2bZsrrrgia6+9dpHSAgAA31bJlxrWIANYvv3jH//Iyy+/XG1swoQJ+ec//5ktt9yySKkAWNH069cv/fr1W+x827Ztc8sttyzDRAAAQF2oV+wAAKy45s+fnzPPPHORc2eeeWbmz5+/jBMBAAAAUMqUGgDUmWeeeSYzZ85c5NzMmTPzzDPPLONEAAAAAJSy5Wb5qX//+9/FjgDAUta9e/c0bdp0kcXGqquumu7duxchFQAAAAClypUaANSZevXq5bTTTlvk3ODBg1Ovnn+GAAAAAKg5v00CoE5tueWW6dChQ7Wxjh07ZosttihSIgAAAABKlVIDgDp31llnVV2VUa9evcXePBwAAAAAvo5SA4A6t9pqq6Vv376pV69e+vbtm9VWW63YkQAAAAAoQcvNjcIBWLEdfvjhOfzww4sdAwAAAIAS5koNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgOAZWLcuHHp3bt3xo0bV+woAAAAAJQopQYAdW727Nk566yzMnXq1Jx11lmZPXt2sSMBAAAAUIKUGgDUueuuuy6zZs1KksyaNSvXX399kRMBAAAAUIqUGgDUqcmTJ+fOO++sNnbHHXdk8uTJRUoEAAAAQKlSagBQZwqFQs4888xFzp155pkpFArLOBEAAAAApUypAUCdeeedd/Laa68tcu61117LO++8s2wDAQAAAFDSlBoA1Jn3339/ieYBAAAAYEFKDQDqzFZbbZUmTZoscq5JkybZaqutlnEiAAAAAEqZUgOAOlOvXr2cfvrpi5w788wzU6+ef4YAAAAAqDm/TQKgTm255ZZp0aJFtbE11lgjW2yxRZESAQAAAFCqlBoA1KnJkydnxowZ1camT5+eyZMnFykRAAAAAKVKqQFAnSkUChk+fHgKhUKNxgEAAADg6yg1AKgzEydOzPjx4xc5N378+EycOHEZJwIAAACglCk1AKgzrVu3Tnl5+SLnysvL07p162WcCAAAAIBSptQAoM48++yzqaysXORcZWVlnn322WWcCAAAAIBSptQAoM5069bta6/U6Nat2zJOBAAAAEApU2oAUGcmTZr0tVdqTJo0aRknAgAAAKCUKTUAqDPz589fonkAAAAAWJBSA4A68/777y/RPAAAAAAsSKkBQJ1Ze+21l2geAAAAABak1ACgzrRp0yZlZWWLnCsrK0ubNm2WcSIAAAAASplSA4A68/TTT6dQKCxyrlAo5Omnn17GiQAAAAAoZUoNAOrM4gqNms4DAAAAwIKUGgDUGffUAAAAAGBpUmoAAAAAAAAlQakBQJ15//33l2geAAAAABak1ACgzqy11lpLNA8AAAAAC1JqAFBnpk6dukTzAAAAALAgpQYAdWarrbZKkyZNFjnXpEmTbLXVVss4EQAAAAClTKkBQJ2pV69edt9990XO7b777qlXzz9DAAAAANSc3yYBUGcqKytzxx13LHLujjvuSGVl5TJOBAAAAEApU2oAUGfuu+++zJ8/f5Fz8+fPz3333beMEwEAAABQypQaANSZ5s2bL9E8AAAAACxIqQFAnSkrK1uieQAAAABYkFIDgDqzuKWnajoPAAAAAAtSagBQZ1577bUlmgcAAACABSk1AKgz8+bNW6J5AAAAAFiQUgOAOlNeXr5E8wAAAACwIKUGAHVmk002WaJ5AAAAAFiQUgOAOtO1a9clmgcAAACABSk1AKgz55133hLNAwAAAMCClBoA1JnevXsv0TwAAAAALEipAUCd+de//rVE8wAAAACwIKUGAHVm0003XaJ5AAAAAFiQUgOAOnPvvfcu0TwAAAAALEipAUCdmTdv3hLNAwAAAMCClBoA1Jnp06cv0TwAAAAALEipAUCd2WOPPZZoHgAAAAAWpNQAoM6Ul5cv0TwAAAAALEipAUCdmTp16hLNAwAAAMCClBoA1JmPPvpoieYBAAAAYEFKDQDqTLNmzZZoHgAAAAAWpNQAoM40b958ieYBAAAAYEFKDQDqzAMPPLBE8wAAAACwIKUGAHWmRYsWSzQPAAAAAAtSagBQZz777LMlmgcAAACABSk1AKgzlZWVSzQPAAAAAAtSagBQZz744IMlmgcAAACABSk1AKgza6yxxhLNAwAAAMCClBoA1Jl///vfSzQPAAAAAAtSagBQZ8rKypZoHgAAAAAWpNQAoM7MnDlzieYBAAAAYEFKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAk1C92AGD5VygUMnv27GLHYAU1a9asYkegBFVUVKSsrKzYMQAAAIBlTKkBfK1CoZDjjjsuL7/8crGjsILaddddix2BErTZZpvlsssuU2wAAADAd4zlpwAAAAAAgJLgSg3ga5WVleWyyy6z/BTfyqRJk9KvX7/Fzo8YMSLrrLPOMkzEisLyUwAAAPDdpNQAvlFZWVkaNWpU7BiUoI033niJ5gEAAABgQUqNFZQbOwPLizFjxizyvhljxoxxk3Cg6Fz1AwAAUFqKWmrMmTMnZ5xxRh555JFUVFTksMMOy2GHHVbMSCuM2bNnu/kusFzzMwpYHowZM8bViN8hzj8AAKD0FbXUOP/88/Pyyy/npptuypQpU3LSSSdl7bXXzi677FLMWAAAwArI+QcAAJS+opUaX3zxRe68885cc8012XTTTbPpppvm9ddfz8iRI51ULAWFQqHYEQAAlnuOmb47nH8AAMCKoWilxquvvpp58+alc+fOVWNdunTJ1Vdfnfnz56devXrFigYAAKxgnH8sW2Xz50VlyLdWKCTz5xU7BcD/U69+4j5sLIEy/64tVUUrNaZNm5bVV189DRo0qBpr0aJF5syZk08++STNmjWr0X4qKyvrKmJJmz9/frEjAAAs9+bPn+948musSO+N84+6t+B70+SFW4uYBABg+VVZWemYcjFq+r4UrdSYNWtWtROKJFWP586dW+P9TJgwYanmWlEUCoWcc845xY7BCqBQKGTEiBF59913ix0FoErbtm3Tr1+/lPm0FEvo1Vdf9ffoO8L5R92bM2dOsSMAACz3XnrppTRs2LDYMUpa0UqNhg0bLnTy8NXjioqKGu+nQ4cOKS8vX6rZgOq6d++e2bNnFzsGJe6rdev98pCloaKiwt8lWAYqKytXmF/iO/+oe4VCIQ888ECxY7ACKBQKSjJgudKwYUPnHyw1zmcXr6bnH0UrNVq2bJmPP/448+bNS/36/40xbdq0VFRUpGnTpjXeT3l5uZMKWAaaNGlS7AgAAN+a849lwzEjS8sqq6xS7AgAwHKqaHfD22STTVK/fv288MILVWPPPfdcOnTo4CZ9AADAUuX8AwAAVgxFO3pv1KhR9tprr5x++ul56aWX8pe//CXXX399Dj744GJFAgAAVlDOPwAAYMVQtOWnkmTQoEE5/fTTc8ghh6RJkyY57rjjstNOOxUzEgAAsIJy/gEAAKWvrPDVnVtLTGVlZV544YVsvvnm1rQFAIA64Jj7//FeAABA3arpMbfFYwEAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSoNQAAAAAAABKglIDAAAAAAAoCUoNAAAAAACgJCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAEqCUgMAAAAAACgJSg0AAAAAAKAkKDUAAAAAAICSUL/YAb6tQqGQJKmsrCxyEgAAWDF9daz91bH3d5nzDwAAqFs1Pf8o2VJj/vz5SZIJEyYUOQkAAKzYvjr2/i5z/gEAAMvGN51/lBVK9GNX8+fPz7x581KvXr2UlZUVOw4AAKxwCoVC5s+fn/r166deve/2yrXOPwAAoG7V9PyjZEsNAAAAAADgu+W7/XErAAAAAACgZCg1AAAAAACAkqDUAAAAAAAASoJSAwAAAAAAKAlKDQAAAAAAoCQoNQAAAAAAgJKg1AAAAAAAAErC/wdFih06bf0TwAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 1600x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# bloxplot\n",
    "# Set the style of seaborn\n",
    "sns.set_style(\"whitegrid\")\n",
    "\n",
    "# Creating boxplots for 'SiteEnergyUse(kBtu)' and 'ENERGYSTARScore'\n",
    "fig, ax = plt.subplots(1, 2, figsize=(16, 6))\n",
    "\n",
    "# Site Energy Use (kBtu)\n",
    "sns.boxplot(data=data, y='SiteEnergyUse(kBtu)', ax=ax[0])\n",
    "ax[0].set_title('Site Energy Use (kBtu)')\n",
    "ax[0].set_ylabel('Energy Use (kBtu)')\n",
    "\n",
    "# ENERGYSTARScore\n",
    "sns.boxplot(data=data, y='ENERGYSTARScore', ax=ax[1])\n",
    "ax[1].set_title('ENERGYSTARScore')\n",
    "ax[1].set_ylabel('Score')\n",
    "\n",
    "plt.tight_layout()\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th>ComplianceStatus</th>\n",
       "      <th>Compliant</th>\n",
       "      <th>Contact the Help Desk</th>\n",
       "      <th>Not Compliant</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EPAPropertyType</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Adult Education</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Automobile Dealership</th>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Bank Branch</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>College/University</th>\n",
       "      <td>37</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Convenience Store without Gas Station</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Convention Center</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Courthouse</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Data Center</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Distribution Center</th>\n",
       "      <td>41</td>\n",
       "      <td>0</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Enclosed Mall</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Energy/Power Station</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Financial Office</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fire Station</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fitness Center/Health Club/Gym</th>\n",
       "      <td>7</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Food Service</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Hospital (General Medical &amp; Surgical)</th>\n",
       "      <td>10</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Hotel</th>\n",
       "      <td>76</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>K-12 School</th>\n",
       "      <td>147</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Laboratory</th>\n",
       "      <td>14</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Library</th>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Mailing Center/Post Office</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Manufacturing/Industrial Plant</th>\n",
       "      <td>13</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Medical Office</th>\n",
       "      <td>36</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Mixed Use Property</th>\n",
       "      <td>118</td>\n",
       "      <td>0</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Movie Theater</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily HR (10+)</th>\n",
       "      <td>117</td>\n",
       "      <td>2</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily Housing</th>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily LR (1-4)</th>\n",
       "      <td>1024</td>\n",
       "      <td>0</td>\n",
       "      <td>25</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily MR (5-9)</th>\n",
       "      <td>699</td>\n",
       "      <td>3</td>\n",
       "      <td>29</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Museum</th>\n",
       "      <td>9</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Non-Refrigerated Warehouse</th>\n",
       "      <td>143</td>\n",
       "      <td>0</td>\n",
       "      <td>9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Office</th>\n",
       "      <td>469</td>\n",
       "      <td>1</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other</th>\n",
       "      <td>50</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Education</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Entertainment/Public Assembly</th>\n",
       "      <td>14</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Lodging/Residential</th>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Mall</th>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Public Services</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Recreation</th>\n",
       "      <td>23</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Restaurant/Bar</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Services</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Utility</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other/Specialty Hospital</th>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Outpatient Rehabilitation/Physical Therapy</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Parking</th>\n",
       "      <td>18</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Performing Arts</th>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Personal Services (Health/Beauty, Dry Cleaning, etc)</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Police Station</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Pre-school/Daycare</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Prison/Incarceration</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Refrigerated Warehouse</th>\n",
       "      <td>14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Repair Services (Vehicle, Shoe, Locksmith, etc)</th>\n",
       "      <td>11</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Residence Hall/Dormitory</th>\n",
       "      <td>16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Residential Care Facility</th>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Restaurant</th>\n",
       "      <td>7</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Retail Store</th>\n",
       "      <td>60</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Self-Storage Facility</th>\n",
       "      <td>46</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Social/Meeting Hall</th>\n",
       "      <td>12</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Stadium (Open)</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Strip Mall</th>\n",
       "      <td>9</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Supermarket/Grocery Store</th>\n",
       "      <td>36</td>\n",
       "      <td>0</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Transportation Terminal/Station</th>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Urgent Care/Clinic/Other Outpatient</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Wholesale Club/Supercenter</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Worship Facility</th>\n",
       "      <td>64</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "ComplianceStatus                                    Compliant   \n",
       "EPAPropertyType                                                 \n",
       "Adult Education                                             3  \\\n",
       "Automobile Dealership                                       6   \n",
       "Bank Branch                                                 1   \n",
       "College/University                                         37   \n",
       "Convenience Store without Gas Station                       0   \n",
       "Convention Center                                           1   \n",
       "Courthouse                                                  3   \n",
       "Data Center                                                 3   \n",
       "Distribution Center                                        41   \n",
       "Enclosed Mall                                               2   \n",
       "Energy/Power Station                                        1   \n",
       "Financial Office                                            1   \n",
       "Fire Station                                                3   \n",
       "Fitness Center/Health Club/Gym                              7   \n",
       "Food Service                                                1   \n",
       "Hospital (General Medical & Surgical)                      10   \n",
       "Hotel                                                      76   \n",
       "K-12 School                                               147   \n",
       "Laboratory                                                 14   \n",
       "Library                                                     4   \n",
       "Mailing Center/Post Office                                  1   \n",
       "Manufacturing/Industrial Plant                             13   \n",
       "Medical Office                                             36   \n",
       "Mixed Use Property                                        118   \n",
       "Movie Theater                                               1   \n",
       "Multifamily HR (10+)                                      117   \n",
       "Multifamily Housing                                         5   \n",
       "Multifamily LR (1-4)                                     1024   \n",
       "Multifamily MR (5-9)                                      699   \n",
       "Museum                                                      9   \n",
       "Non-Refrigerated Warehouse                                143   \n",
       "Office                                                    469   \n",
       "Other                                                      50   \n",
       "Other - Education                                           1   \n",
       "Other - Entertainment/Public Assembly                      14   \n",
       "Other - Lodging/Residential                                 5   \n",
       "Other - Mall                                                5   \n",
       "Other - Public Services                                     2   \n",
       "Other - Recreation                                         23   \n",
       "Other - Restaurant/Bar                                      1   \n",
       "Other - Services                                            3   \n",
       "Other - Utility                                             3   \n",
       "Other/Specialty Hospital                                    4   \n",
       "Outpatient Rehabilitation/Physical Therapy                  1   \n",
       "Parking                                                    18   \n",
       "Performing Arts                                             6   \n",
       "Personal Services (Health/Beauty, Dry Cleaning,...          1   \n",
       "Police Station                                              3   \n",
       "Pre-school/Daycare                                          2   \n",
       "Prison/Incarceration                                        1   \n",
       "Refrigerated Warehouse                                     14   \n",
       "Repair Services (Vehicle, Shoe, Locksmith, etc)            11   \n",
       "Residence Hall/Dormitory                                   16   \n",
       "Residential Care Facility                                   4   \n",
       "Restaurant                                                  7   \n",
       "Retail Store                                               60   \n",
       "Self-Storage Facility                                      46   \n",
       "Social/Meeting Hall                                        12   \n",
       "Stadium (Open)                                              1   \n",
       "Strip Mall                                                  9   \n",
       "Supermarket/Grocery Store                                  36   \n",
       "Transportation Terminal/Station                             4   \n",
       "Urgent Care/Clinic/Other Outpatient                         2   \n",
       "Wholesale Club/Supercenter                                  1   \n",
       "Worship Facility                                           64   \n",
       "\n",
       "ComplianceStatus                                    Contact the Help Desk   \n",
       "EPAPropertyType                                                             \n",
       "Adult Education                                                         0  \\\n",
       "Automobile Dealership                                                   0   \n",
       "Bank Branch                                                             0   \n",
       "College/University                                                      0   \n",
       "Convenience Store without Gas Station                                   0   \n",
       "Convention Center                                                       0   \n",
       "Courthouse                                                              0   \n",
       "Data Center                                                             0   \n",
       "Distribution Center                                                     0   \n",
       "Enclosed Mall                                                           0   \n",
       "Energy/Power Station                                                    0   \n",
       "Financial Office                                                        0   \n",
       "Fire Station                                                            0   \n",
       "Fitness Center/Health Club/Gym                                          0   \n",
       "Food Service                                                            0   \n",
       "Hospital (General Medical & Surgical)                                   0   \n",
       "Hotel                                                                   0   \n",
       "K-12 School                                                             0   \n",
       "Laboratory                                                              0   \n",
       "Library                                                                 0   \n",
       "Mailing Center/Post Office                                              0   \n",
       "Manufacturing/Industrial Plant                                          0   \n",
       "Medical Office                                                          0   \n",
       "Mixed Use Property                                                      0   \n",
       "Movie Theater                                                           0   \n",
       "Multifamily HR (10+)                                                    2   \n",
       "Multifamily Housing                                                     0   \n",
       "Multifamily LR (1-4)                                                    0   \n",
       "Multifamily MR (5-9)                                                    3   \n",
       "Museum                                                                  0   \n",
       "Non-Refrigerated Warehouse                                              0   \n",
       "Office                                                                  1   \n",
       "Other                                                                   0   \n",
       "Other - Education                                                       0   \n",
       "Other - Entertainment/Public Assembly                                   0   \n",
       "Other - Lodging/Residential                                             0   \n",
       "Other - Mall                                                            0   \n",
       "Other - Public Services                                                 0   \n",
       "Other - Recreation                                                      0   \n",
       "Other - Restaurant/Bar                                                  0   \n",
       "Other - Services                                                        0   \n",
       "Other - Utility                                                         0   \n",
       "Other/Specialty Hospital                                                0   \n",
       "Outpatient Rehabilitation/Physical Therapy                              0   \n",
       "Parking                                                                 0   \n",
       "Performing Arts                                                         0   \n",
       "Personal Services (Health/Beauty, Dry Cleaning,...                      0   \n",
       "Police Station                                                          0   \n",
       "Pre-school/Daycare                                                      0   \n",
       "Prison/Incarceration                                                    0   \n",
       "Refrigerated Warehouse                                                  0   \n",
       "Repair Services (Vehicle, Shoe, Locksmith, etc)                         0   \n",
       "Residence Hall/Dormitory                                                0   \n",
       "Residential Care Facility                                               0   \n",
       "Restaurant                                                              0   \n",
       "Retail Store                                                            0   \n",
       "Self-Storage Facility                                                   0   \n",
       "Social/Meeting Hall                                                     0   \n",
       "Stadium (Open)                                                          0   \n",
       "Strip Mall                                                              0   \n",
       "Supermarket/Grocery Store                                               0   \n",
       "Transportation Terminal/Station                                         0   \n",
       "Urgent Care/Clinic/Other Outpatient                                     0   \n",
       "Wholesale Club/Supercenter                                              0   \n",
       "Worship Facility                                                        0   \n",
       "\n",
       "ComplianceStatus                                    Not Compliant  \n",
       "EPAPropertyType                                                    \n",
       "Adult Education                                                 0  \n",
       "Automobile Dealership                                           0  \n",
       "Bank Branch                                                     1  \n",
       "College/University                                              2  \n",
       "Convenience Store without Gas Station                           1  \n",
       "Convention Center                                               0  \n",
       "Courthouse                                                      0  \n",
       "Data Center                                                     0  \n",
       "Distribution Center                                             6  \n",
       "Enclosed Mall                                                   0  \n",
       "Energy/Power Station                                            0  \n",
       "Financial Office                                                0  \n",
       "Fire Station                                                    0  \n",
       "Fitness Center/Health Club/Gym                                  1  \n",
       "Food Service                                                    0  \n",
       "Hospital (General Medical & Surgical)                           0  \n",
       "Hotel                                                           2  \n",
       "K-12 School                                                     1  \n",
       "Laboratory                                                      1  \n",
       "Library                                                         0  \n",
       "Mailing Center/Post Office                                      0  \n",
       "Manufacturing/Industrial Plant                                  1  \n",
       "Medical Office                                                 10  \n",
       "Mixed Use Property                                              6  \n",
       "Movie Theater                                                   0  \n",
       "Multifamily HR (10+)                                            6  \n",
       "Multifamily Housing                                             3  \n",
       "Multifamily LR (1-4)                                           25  \n",
       "Multifamily MR (5-9)                                           29  \n",
       "Museum                                                          0  \n",
       "Non-Refrigerated Warehouse                                      9  \n",
       "Office                                                          8  \n",
       "Other                                                           1  \n",
       "Other - Education                                               0  \n",
       "Other - Entertainment/Public Assembly                           1  \n",
       "Other - Lodging/Residential                                     0  \n",
       "Other - Mall                                                    0  \n",
       "Other - Public Services                                         0  \n",
       "Other - Recreation                                              0  \n",
       "Other - Restaurant/Bar                                          0  \n",
       "Other - Services                                                0  \n",
       "Other - Utility                                                 0  \n",
       "Other/Specialty Hospital                                        0  \n",
       "Outpatient Rehabilitation/Physical Therapy                      0  \n",
       "Parking                                                         2  \n",
       "Performing Arts                                                 0  \n",
       "Personal Services (Health/Beauty, Dry Cleaning,...              0  \n",
       "Police Station                                                  0  \n",
       "Pre-school/Daycare                                              1  \n",
       "Prison/Incarceration                                            0  \n",
       "Refrigerated Warehouse                                          0  \n",
       "Repair Services (Vehicle, Shoe, Locksmith, etc)                 0  \n",
       "Residence Hall/Dormitory                                        0  \n",
       "Residential Care Facility                                       0  \n",
       "Restaurant                                                      0  \n",
       "Retail Store                                                   10  \n",
       "Self-Storage Facility                                           1  \n",
       "Social/Meeting Hall                                             0  \n",
       "Stadium (Open)                                                  0  \n",
       "Strip Mall                                                      0  \n",
       "Supermarket/Grocery Store                                       3  \n",
       "Transportation Terminal/Station                                 0  \n",
       "Urgent Care/Clinic/Other Outpatient                             0  \n",
       "Wholesale Club/Supercenter                                      0  \n",
       "Worship Facility                                                2  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# compliance status distribution 'EPAPropertyType' and 'ComplianceStatus'\n",
    "compliance_status_grouped = data.groupby(['EPAPropertyType', 'ComplianceStatus']).size().unstack(fill_value=0)\n",
    "compliance_status_grouped"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th>ComplianceStatus</th>\n",
       "      <th>Compliant</th>\n",
       "      <th>Contact the Help Desk</th>\n",
       "      <th>Not Compliant</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EPAPropertyType</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Convenience Store without Gas Station</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>100.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Bank Branch</th>\n",
       "      <td>50.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>50.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily Housing</th>\n",
       "      <td>62.500000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>37.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Pre-school/Daycare</th>\n",
       "      <td>66.666667</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>33.333333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Medical Office</th>\n",
       "      <td>78.260870</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>21.739130</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Retail Store</th>\n",
       "      <td>85.714286</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>14.285714</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Distribution Center</th>\n",
       "      <td>87.234043</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>12.765957</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fitness Center/Health Club/Gym</th>\n",
       "      <td>87.500000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>12.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Parking</th>\n",
       "      <td>90.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>10.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Supermarket/Grocery Store</th>\n",
       "      <td>92.307692</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>7.692308</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Manufacturing/Industrial Plant</th>\n",
       "      <td>92.857143</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>7.142857</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Laboratory</th>\n",
       "      <td>93.333333</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>6.666667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Entertainment/Public Assembly</th>\n",
       "      <td>93.333333</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>6.666667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Non-Refrigerated Warehouse</th>\n",
       "      <td>94.078947</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>5.921053</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>College/University</th>\n",
       "      <td>94.871795</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>5.128205</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Mixed Use Property</th>\n",
       "      <td>95.161290</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>4.838710</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily HR (10+)</th>\n",
       "      <td>93.600000</td>\n",
       "      <td>1.600000</td>\n",
       "      <td>4.800000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily MR (5-9)</th>\n",
       "      <td>95.622435</td>\n",
       "      <td>0.410397</td>\n",
       "      <td>3.967168</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Worship Facility</th>\n",
       "      <td>96.969697</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>3.030303</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Hotel</th>\n",
       "      <td>97.435897</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.564103</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Multifamily LR (1-4)</th>\n",
       "      <td>97.616778</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.383222</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Self-Storage Facility</th>\n",
       "      <td>97.872340</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.127660</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other</th>\n",
       "      <td>98.039216</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.960784</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Office</th>\n",
       "      <td>98.117155</td>\n",
       "      <td>0.209205</td>\n",
       "      <td>1.673640</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>K-12 School</th>\n",
       "      <td>99.324324</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.675676</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Mailing Center/Post Office</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Data Center</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Police Station</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Enclosed Mall</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Prison/Incarceration</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Refrigerated Warehouse</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Repair Services (Vehicle, Shoe, Locksmith, etc)</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Residence Hall/Dormitory</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Residential Care Facility</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Restaurant</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Courthouse</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Performing Arts</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Social/Meeting Hall</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Stadium (Open)</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Strip Mall</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Convention Center</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Transportation Terminal/Station</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Urgent Care/Clinic/Other Outpatient</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Wholesale Club/Supercenter</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Personal Services (Health/Beauty, Dry Cleaning, etc)</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Energy/Power Station</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Movie Theater</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Financial Office</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Library</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Hospital (General Medical &amp; Surgical)</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Food Service</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Museum</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Fire Station</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Automobile Dealership</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Education</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Lodging/Residential</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Outpatient Rehabilitation/Physical Therapy</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Mall</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Public Services</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Recreation</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Restaurant/Bar</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Services</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other - Utility</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other/Specialty Hospital</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Adult Education</th>\n",
       "      <td>100.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "ComplianceStatus                                     Compliant   \n",
       "EPAPropertyType                                                  \n",
       "Convenience Store without Gas Station                 0.000000  \\\n",
       "Bank Branch                                          50.000000   \n",
       "Multifamily Housing                                  62.500000   \n",
       "Pre-school/Daycare                                   66.666667   \n",
       "Medical Office                                       78.260870   \n",
       "Retail Store                                         85.714286   \n",
       "Distribution Center                                  87.234043   \n",
       "Fitness Center/Health Club/Gym                       87.500000   \n",
       "Parking                                              90.000000   \n",
       "Supermarket/Grocery Store                            92.307692   \n",
       "Manufacturing/Industrial Plant                       92.857143   \n",
       "Laboratory                                           93.333333   \n",
       "Other - Entertainment/Public Assembly                93.333333   \n",
       "Non-Refrigerated Warehouse                           94.078947   \n",
       "College/University                                   94.871795   \n",
       "Mixed Use Property                                   95.161290   \n",
       "Multifamily HR (10+)                                 93.600000   \n",
       "Multifamily MR (5-9)                                 95.622435   \n",
       "Worship Facility                                     96.969697   \n",
       "Hotel                                                97.435897   \n",
       "Multifamily LR (1-4)                                 97.616778   \n",
       "Self-Storage Facility                                97.872340   \n",
       "Other                                                98.039216   \n",
       "Office                                               98.117155   \n",
       "K-12 School                                          99.324324   \n",
       "Mailing Center/Post Office                          100.000000   \n",
       "Data Center                                         100.000000   \n",
       "Police Station                                      100.000000   \n",
       "Enclosed Mall                                       100.000000   \n",
       "Prison/Incarceration                                100.000000   \n",
       "Refrigerated Warehouse                              100.000000   \n",
       "Repair Services (Vehicle, Shoe, Locksmith, etc)     100.000000   \n",
       "Residence Hall/Dormitory                            100.000000   \n",
       "Residential Care Facility                           100.000000   \n",
       "Restaurant                                          100.000000   \n",
       "Courthouse                                          100.000000   \n",
       "Performing Arts                                     100.000000   \n",
       "Social/Meeting Hall                                 100.000000   \n",
       "Stadium (Open)                                      100.000000   \n",
       "Strip Mall                                          100.000000   \n",
       "Convention Center                                   100.000000   \n",
       "Transportation Terminal/Station                     100.000000   \n",
       "Urgent Care/Clinic/Other Outpatient                 100.000000   \n",
       "Wholesale Club/Supercenter                          100.000000   \n",
       "Personal Services (Health/Beauty, Dry Cleaning,...  100.000000   \n",
       "Energy/Power Station                                100.000000   \n",
       "Movie Theater                                       100.000000   \n",
       "Financial Office                                    100.000000   \n",
       "Library                                             100.000000   \n",
       "Hospital (General Medical & Surgical)               100.000000   \n",
       "Food Service                                        100.000000   \n",
       "Museum                                              100.000000   \n",
       "Fire Station                                        100.000000   \n",
       "Automobile Dealership                               100.000000   \n",
       "Other - Education                                   100.000000   \n",
       "Other - Lodging/Residential                         100.000000   \n",
       "Outpatient Rehabilitation/Physical Therapy          100.000000   \n",
       "Other - Mall                                        100.000000   \n",
       "Other - Public Services                             100.000000   \n",
       "Other - Recreation                                  100.000000   \n",
       "Other - Restaurant/Bar                              100.000000   \n",
       "Other - Services                                    100.000000   \n",
       "Other - Utility                                     100.000000   \n",
       "Other/Specialty Hospital                            100.000000   \n",
       "Adult Education                                     100.000000   \n",
       "\n",
       "ComplianceStatus                                    Contact the Help Desk   \n",
       "EPAPropertyType                                                             \n",
       "Convenience Store without Gas Station                            0.000000  \\\n",
       "Bank Branch                                                      0.000000   \n",
       "Multifamily Housing                                              0.000000   \n",
       "Pre-school/Daycare                                               0.000000   \n",
       "Medical Office                                                   0.000000   \n",
       "Retail Store                                                     0.000000   \n",
       "Distribution Center                                              0.000000   \n",
       "Fitness Center/Health Club/Gym                                   0.000000   \n",
       "Parking                                                          0.000000   \n",
       "Supermarket/Grocery Store                                        0.000000   \n",
       "Manufacturing/Industrial Plant                                   0.000000   \n",
       "Laboratory                                                       0.000000   \n",
       "Other - Entertainment/Public Assembly                            0.000000   \n",
       "Non-Refrigerated Warehouse                                       0.000000   \n",
       "College/University                                               0.000000   \n",
       "Mixed Use Property                                               0.000000   \n",
       "Multifamily HR (10+)                                             1.600000   \n",
       "Multifamily MR (5-9)                                             0.410397   \n",
       "Worship Facility                                                 0.000000   \n",
       "Hotel                                                            0.000000   \n",
       "Multifamily LR (1-4)                                             0.000000   \n",
       "Self-Storage Facility                                            0.000000   \n",
       "Other                                                            0.000000   \n",
       "Office                                                           0.209205   \n",
       "K-12 School                                                      0.000000   \n",
       "Mailing Center/Post Office                                       0.000000   \n",
       "Data Center                                                      0.000000   \n",
       "Police Station                                                   0.000000   \n",
       "Enclosed Mall                                                    0.000000   \n",
       "Prison/Incarceration                                             0.000000   \n",
       "Refrigerated Warehouse                                           0.000000   \n",
       "Repair Services (Vehicle, Shoe, Locksmith, etc)                  0.000000   \n",
       "Residence Hall/Dormitory                                         0.000000   \n",
       "Residential Care Facility                                        0.000000   \n",
       "Restaurant                                                       0.000000   \n",
       "Courthouse                                                       0.000000   \n",
       "Performing Arts                                                  0.000000   \n",
       "Social/Meeting Hall                                              0.000000   \n",
       "Stadium (Open)                                                   0.000000   \n",
       "Strip Mall                                                       0.000000   \n",
       "Convention Center                                                0.000000   \n",
       "Transportation Terminal/Station                                  0.000000   \n",
       "Urgent Care/Clinic/Other Outpatient                              0.000000   \n",
       "Wholesale Club/Supercenter                                       0.000000   \n",
       "Personal Services (Health/Beauty, Dry Cleaning,...               0.000000   \n",
       "Energy/Power Station                                             0.000000   \n",
       "Movie Theater                                                    0.000000   \n",
       "Financial Office                                                 0.000000   \n",
       "Library                                                          0.000000   \n",
       "Hospital (General Medical & Surgical)                            0.000000   \n",
       "Food Service                                                     0.000000   \n",
       "Museum                                                           0.000000   \n",
       "Fire Station                                                     0.000000   \n",
       "Automobile Dealership                                            0.000000   \n",
       "Other - Education                                                0.000000   \n",
       "Other - Lodging/Residential                                      0.000000   \n",
       "Outpatient Rehabilitation/Physical Therapy                       0.000000   \n",
       "Other - Mall                                                     0.000000   \n",
       "Other - Public Services                                          0.000000   \n",
       "Other - Recreation                                               0.000000   \n",
       "Other - Restaurant/Bar                                           0.000000   \n",
       "Other - Services                                                 0.000000   \n",
       "Other - Utility                                                  0.000000   \n",
       "Other/Specialty Hospital                                         0.000000   \n",
       "Adult Education                                                  0.000000   \n",
       "\n",
       "ComplianceStatus                                    Not Compliant  \n",
       "EPAPropertyType                                                    \n",
       "Convenience Store without Gas Station                  100.000000  \n",
       "Bank Branch                                             50.000000  \n",
       "Multifamily Housing                                     37.500000  \n",
       "Pre-school/Daycare                                      33.333333  \n",
       "Medical Office                                          21.739130  \n",
       "Retail Store                                            14.285714  \n",
       "Distribution Center                                     12.765957  \n",
       "Fitness Center/Health Club/Gym                          12.500000  \n",
       "Parking                                                 10.000000  \n",
       "Supermarket/Grocery Store                                7.692308  \n",
       "Manufacturing/Industrial Plant                           7.142857  \n",
       "Laboratory                                               6.666667  \n",
       "Other - Entertainment/Public Assembly                    6.666667  \n",
       "Non-Refrigerated Warehouse                               5.921053  \n",
       "College/University                                       5.128205  \n",
       "Mixed Use Property                                       4.838710  \n",
       "Multifamily HR (10+)                                     4.800000  \n",
       "Multifamily MR (5-9)                                     3.967168  \n",
       "Worship Facility                                         3.030303  \n",
       "Hotel                                                    2.564103  \n",
       "Multifamily LR (1-4)                                     2.383222  \n",
       "Self-Storage Facility                                    2.127660  \n",
       "Other                                                    1.960784  \n",
       "Office                                                   1.673640  \n",
       "K-12 School                                              0.675676  \n",
       "Mailing Center/Post Office                               0.000000  \n",
       "Data Center                                              0.000000  \n",
       "Police Station                                           0.000000  \n",
       "Enclosed Mall                                            0.000000  \n",
       "Prison/Incarceration                                     0.000000  \n",
       "Refrigerated Warehouse                                   0.000000  \n",
       "Repair Services (Vehicle, Shoe, Locksmith, etc)          0.000000  \n",
       "Residence Hall/Dormitory                                 0.000000  \n",
       "Residential Care Facility                                0.000000  \n",
       "Restaurant                                               0.000000  \n",
       "Courthouse                                               0.000000  \n",
       "Performing Arts                                          0.000000  \n",
       "Social/Meeting Hall                                      0.000000  \n",
       "Stadium (Open)                                           0.000000  \n",
       "Strip Mall                                               0.000000  \n",
       "Convention Center                                        0.000000  \n",
       "Transportation Terminal/Station                          0.000000  \n",
       "Urgent Care/Clinic/Other Outpatient                      0.000000  \n",
       "Wholesale Club/Supercenter                               0.000000  \n",
       "Personal Services (Health/Beauty, Dry Cleaning,...       0.000000  \n",
       "Energy/Power Station                                     0.000000  \n",
       "Movie Theater                                            0.000000  \n",
       "Financial Office                                         0.000000  \n",
       "Library                                                  0.000000  \n",
       "Hospital (General Medical & Surgical)                    0.000000  \n",
       "Food Service                                             0.000000  \n",
       "Museum                                                   0.000000  \n",
       "Fire Station                                             0.000000  \n",
       "Automobile Dealership                                    0.000000  \n",
       "Other - Education                                        0.000000  \n",
       "Other - Lodging/Residential                              0.000000  \n",
       "Outpatient Rehabilitation/Physical Therapy               0.000000  \n",
       "Other - Mall                                             0.000000  \n",
       "Other - Public Services                                  0.000000  \n",
       "Other - Recreation                                       0.000000  \n",
       "Other - Restaurant/Bar                                   0.000000  \n",
       "Other - Services                                         0.000000  \n",
       "Other - Utility                                          0.000000  \n",
       "Other/Specialty Hospital                                 0.000000  \n",
       "Adult Education                                          0.000000  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "compliance_status_percentage = compliance_status_grouped.apply(lambda x: (x / x.sum()) * 100, axis=1)\n",
    "\n",
    "compliance_status_percentage_sorted = compliance_status_percentage.sort_values(by='Not Compliant', ascending=False)\n",
    "compliance_status_percentage_sorted"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "ComplianceIssue\n",
       "No Issue                                3466\n",
       "Portfolio Manager Account Not Shared      58\n",
       "Account Requires Verification             49\n",
       "Default Data                              43\n",
       "Missing 2021 EUI or Electricity Data      35\n",
       "Unknown - Contact the Help Desk           12\n",
       "Name: count, dtype: int64"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Calculate the frequency of each category in 'ComplianceIssue' column\n",
    "compliance_issue_freq = data['ComplianceIssue'].value_counts()\n",
    "\n",
    "# Display the frequency distribution\n",
    "compliance_issue_freq\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'Outliers in SiteEnergyUse(kBtu)': 35, 'Inconsistencies in GFA Reporting': 0}\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>OSEBuildingID</th>\n",
       "      <th>BuildingName</th>\n",
       "      <th>SiteEnergyUse(kBtu)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>26</td>\n",
       "      <td>KING COUNTY COURTHOUSE</td>\n",
       "      <td>63318804</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>43</td>\n",
       "      <td>BOEING PLANT 2 ISAL  (BLDG 2-122)-North Boeing...</td>\n",
       "      <td>518003488</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48</th>\n",
       "      <td>63</td>\n",
       "      <td>SHERATON GRAND SEATTLE - (Parent Building)</td>\n",
       "      <td>56395556</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95</th>\n",
       "      <td>147</td>\n",
       "      <td>SWEDISH MEDICAL CENTER BALLARD - CAMPUS</td>\n",
       "      <td>63914268</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>116</th>\n",
       "      <td>198</td>\n",
       "      <td>SEATTLE CHILDREN'S HOSPITAL  (MAIN CAMPUS)</td>\n",
       "      <td>230595392</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>145</th>\n",
       "      <td>245</td>\n",
       "      <td>CENTURYLINK FIELD-STADIUM</td>\n",
       "      <td>59420744</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>159</th>\n",
       "      <td>268</td>\n",
       "      <td>SWEDISH MEDICAL CENTER CHERRY HILL - CAMPUS</td>\n",
       "      <td>141621536</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>162</th>\n",
       "      <td>276</td>\n",
       "      <td>HARBORVIEW MEDICAL CENTER - Campus</td>\n",
       "      <td>240119280</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>179</th>\n",
       "      <td>295</td>\n",
       "      <td>CENTURYLINK: 3RD AVE 3-D</td>\n",
       "      <td>54317328</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>209</th>\n",
       "      <td>338</td>\n",
       "      <td>SEATTLE CHILDREN'S HOSPITAL  (JACK R MCDONALD ...</td>\n",
       "      <td>58116100</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>296</th>\n",
       "      <td>435</td>\n",
       "      <td>SEATTLE CONVENTION CENTER</td>\n",
       "      <td>54539952</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>329</th>\n",
       "      <td>477</td>\n",
       "      <td>FHCRC - Weintraub/Hutchinson/Thomas Bldgs comb...</td>\n",
       "      <td>144551776</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>457</th>\n",
       "      <td>640</td>\n",
       "      <td>STARBUCKS (SODO) CENTER - Campus</td>\n",
       "      <td>61798728</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>475</th>\n",
       "      <td>675</td>\n",
       "      <td>HARBORVIEW MEDICAL CENTER - Ninth &amp; Jefferson ...</td>\n",
       "      <td>61187836</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>538</th>\n",
       "      <td>753</td>\n",
       "      <td>THE WESTIN BUILDING EXCHANGE (Office)</td>\n",
       "      <td>270897440</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>559</th>\n",
       "      <td>775</td>\n",
       "      <td>COLUMBIA CENTER</td>\n",
       "      <td>74236416</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>597</th>\n",
       "      <td>828</td>\n",
       "      <td>SWEDISH MEDICAL CENTER FIRST HILL - Campus</td>\n",
       "      <td>190562112</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1266</th>\n",
       "      <td>21476</td>\n",
       "      <td>HARBOR STEPS NE SITE</td>\n",
       "      <td>56869144</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1464</th>\n",
       "      <td>22062</td>\n",
       "      <td>SEATTLE UNIVERSITY (MAIN CAMPUS)</td>\n",
       "      <td>95224416</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1644</th>\n",
       "      <td>23071</td>\n",
       "      <td>UW NORTHWEST HOSPITAL &amp; MEDICAL CENTER - CAMPUS</td>\n",
       "      <td>95649288</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1662</th>\n",
       "      <td>23113</td>\n",
       "      <td>KAISER PERMANENTE CAPITOL HILL - Campus</td>\n",
       "      <td>110285840</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1707</th>\n",
       "      <td>23271</td>\n",
       "      <td>PORT OF SEATTLE- TERMINAL 18 (CAMPUS)</td>\n",
       "      <td>66171828</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2285</th>\n",
       "      <td>25251</td>\n",
       "      <td>SEATTLE PACIFIC UNIVERSITY - MASTER CAMPUS</td>\n",
       "      <td>57673908</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2778</th>\n",
       "      <td>27100</td>\n",
       "      <td>FOUR SEASONS HOTEL &amp; 99 UNION STREET PRIVATE R...</td>\n",
       "      <td>52529956</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3138</th>\n",
       "      <td>49725</td>\n",
       "      <td>VIRGINIA MASON - Jones Pavilion</td>\n",
       "      <td>62716216</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3143</th>\n",
       "      <td>49732</td>\n",
       "      <td>SEATTLE CHILDREN'S HOSPITAL BUILDING HOPE</td>\n",
       "      <td>59513360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3152</th>\n",
       "      <td>49770</td>\n",
       "      <td>University of Washington Medical Center (UWMC)</td>\n",
       "      <td>325567936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3162</th>\n",
       "      <td>49783</td>\n",
       "      <td>UW Medicine Phase 3</td>\n",
       "      <td>68276376</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3225</th>\n",
       "      <td>49940</td>\n",
       "      <td>VIRGINIA MASON - Central Pavilion</td>\n",
       "      <td>74310120</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3241</th>\n",
       "      <td>49973</td>\n",
       "      <td>UW - FOEGE BUILDING (genome sciences)</td>\n",
       "      <td>111853328</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3244</th>\n",
       "      <td>49976</td>\n",
       "      <td>UW - HEALTH SCIENCES RR-WING</td>\n",
       "      <td>66197036</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3417</th>\n",
       "      <td>50294</td>\n",
       "      <td>DAY 1 TOWER AND SPHERES (Amazon)</td>\n",
       "      <td>52713236</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3476</th>\n",
       "      <td>50376</td>\n",
       "      <td>KOMO PLAZA EAST &amp; WEST BUILDINGS</td>\n",
       "      <td>134579488</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3495</th>\n",
       "      <td>50416</td>\n",
       "      <td>PORT OF SEATTLE- TERMINAL 91 (CAMPUS)</td>\n",
       "      <td>55399684</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3566</th>\n",
       "      <td>50523</td>\n",
       "      <td>HYATT REGENCY SEATTLE</td>\n",
       "      <td>75943616</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      OSEBuildingID                                       BuildingName   \n",
       "20               26                             KING COUNTY COURTHOUSE  \\\n",
       "34               43  BOEING PLANT 2 ISAL  (BLDG 2-122)-North Boeing...   \n",
       "48               63         SHERATON GRAND SEATTLE - (Parent Building)   \n",
       "95              147            SWEDISH MEDICAL CENTER BALLARD - CAMPUS   \n",
       "116             198         SEATTLE CHILDREN'S HOSPITAL  (MAIN CAMPUS)   \n",
       "145             245                          CENTURYLINK FIELD-STADIUM   \n",
       "159             268        SWEDISH MEDICAL CENTER CHERRY HILL - CAMPUS   \n",
       "162             276                 HARBORVIEW MEDICAL CENTER - Campus   \n",
       "179             295                           CENTURYLINK: 3RD AVE 3-D   \n",
       "209             338  SEATTLE CHILDREN'S HOSPITAL  (JACK R MCDONALD ...   \n",
       "296             435                          SEATTLE CONVENTION CENTER   \n",
       "329             477  FHCRC - Weintraub/Hutchinson/Thomas Bldgs comb...   \n",
       "457             640                   STARBUCKS (SODO) CENTER - Campus   \n",
       "475             675  HARBORVIEW MEDICAL CENTER - Ninth & Jefferson ...   \n",
       "538             753              THE WESTIN BUILDING EXCHANGE (Office)   \n",
       "559             775                                    COLUMBIA CENTER   \n",
       "597             828         SWEDISH MEDICAL CENTER FIRST HILL - Campus   \n",
       "1266          21476                               HARBOR STEPS NE SITE   \n",
       "1464          22062                   SEATTLE UNIVERSITY (MAIN CAMPUS)   \n",
       "1644          23071    UW NORTHWEST HOSPITAL & MEDICAL CENTER - CAMPUS   \n",
       "1662          23113            KAISER PERMANENTE CAPITOL HILL - Campus   \n",
       "1707          23271              PORT OF SEATTLE- TERMINAL 18 (CAMPUS)   \n",
       "2285          25251         SEATTLE PACIFIC UNIVERSITY - MASTER CAMPUS   \n",
       "2778          27100  FOUR SEASONS HOTEL & 99 UNION STREET PRIVATE R...   \n",
       "3138          49725                    VIRGINIA MASON - Jones Pavilion   \n",
       "3143          49732          SEATTLE CHILDREN'S HOSPITAL BUILDING HOPE   \n",
       "3152          49770     University of Washington Medical Center (UWMC)   \n",
       "3162          49783                                UW Medicine Phase 3   \n",
       "3225          49940                  VIRGINIA MASON - Central Pavilion   \n",
       "3241          49973              UW - FOEGE BUILDING (genome sciences)   \n",
       "3244          49976                       UW - HEALTH SCIENCES RR-WING   \n",
       "3417          50294                   DAY 1 TOWER AND SPHERES (Amazon)   \n",
       "3476          50376                   KOMO PLAZA EAST & WEST BUILDINGS   \n",
       "3495          50416              PORT OF SEATTLE- TERMINAL 91 (CAMPUS)   \n",
       "3566          50523                              HYATT REGENCY SEATTLE   \n",
       "\n",
       "      SiteEnergyUse(kBtu)  \n",
       "20               63318804  \n",
       "34              518003488  \n",
       "48               56395556  \n",
       "95               63914268  \n",
       "116             230595392  \n",
       "145              59420744  \n",
       "159             141621536  \n",
       "162             240119280  \n",
       "179              54317328  \n",
       "209              58116100  \n",
       "296              54539952  \n",
       "329             144551776  \n",
       "457              61798728  \n",
       "475              61187836  \n",
       "538             270897440  \n",
       "559              74236416  \n",
       "597             190562112  \n",
       "1266             56869144  \n",
       "1464             95224416  \n",
       "1644             95649288  \n",
       "1662            110285840  \n",
       "1707             66171828  \n",
       "2285             57673908  \n",
       "2778             52529956  \n",
       "3138             62716216  \n",
       "3143             59513360  \n",
       "3152            325567936  \n",
       "3162             68276376  \n",
       "3225             74310120  \n",
       "3241            111853328  \n",
       "3244             66197036  \n",
       "3417             52713236  \n",
       "3476            134579488  \n",
       "3495             55399684  \n",
       "3566             75943616  "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Anomaly Detection for 'SiteEnergyUse(kBtu)' using Z-score to identify outliers\n",
    "from scipy.stats import zscore\n",
    "\n",
    "# Compute Z-scores for 'SiteEnergyUse(kBtu)'\n",
    "data['SiteEnergyUse_zscore'] = zscore(data['SiteEnergyUse(kBtu)'].dropna())\n",
    "\n",
    "# Identify outliers as those with a Z-score greater than 3 or less than -3 (typical threshold for outliers)\n",
    "outliers_energy_use = data.loc[data['SiteEnergyUse_zscore'].abs() > 3]\n",
    "\n",
    "# Consistency Check for Gross Floor Area\n",
    "# Checking if the sum of 'PropertyGFABuilding(s)' and 'PropertyGFAParking' equals 'PropertyGFATotal'\n",
    "data['GFA_Check'] = (data['PropertyGFABuilding(s)'] + data['PropertyGFAParking']) == data['PropertyGFATotal']\n",
    "\n",
    "# Summary of findings\n",
    "outliers_summary = {\n",
    "    'Outliers in SiteEnergyUse(kBtu)': outliers_energy_use.shape[0],\n",
    "    'Inconsistencies in GFA Reporting': (~data['GFA_Check']).sum()\n",
    "}\n",
    "\n",
    "print(outliers_summary)\n",
    "\n",
    "# Display the rows corresponding to the 30 outliers identified in 'SiteEnergyUse(kBtu)'\n",
    "outliers_energy_use[['OSEBuildingID', 'BuildingName', 'SiteEnergyUse(kBtu)']]\n",
    "\n",
    "#SWEDISH MEDICAL CENTER CHERRY HILL - CAMPUS (OSEBuildingID: 268) is listed as \"Not Compliant.\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.1"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}