hw2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "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>bad</th>\n",
       "      <th>loan</th>\n",
       "      <th>mortdue</th>\n",
       "      <th>value</th>\n",
       "      <th>reason</th>\n",
       "      <th>job</th>\n",
       "      <th>yoj</th>\n",
       "      <th>derog</th>\n",
       "      <th>delinq</th>\n",
       "      <th>clage</th>\n",
       "      <th>ninq</th>\n",
       "      <th>clno</th>\n",
       "      <th>debtinc</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>153</th>\n",
       "      <td>0</td>\n",
       "      <td>18200</td>\n",
       "      <td>94727.0</td>\n",
       "      <td>136877.0</td>\n",
       "      <td>DebtCon</td>\n",
       "      <td>Mgr</td>\n",
       "      <td>15.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>168.96</td>\n",
       "      <td>2.0</td>\n",
       "      <td>26.0</td>\n",
       "      <td>36.056</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>154</th>\n",
       "      <td>0</td>\n",
       "      <td>21700</td>\n",
       "      <td>79240.0</td>\n",
       "      <td>96784.0</td>\n",
       "      <td>DebtCon</td>\n",
       "      <td>Mgr</td>\n",
       "      <td>5.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.51</td>\n",
       "      <td>6.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>38.079</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>155</th>\n",
       "      <td>0</td>\n",
       "      <td>34100</td>\n",
       "      <td>241931.0</td>\n",
       "      <td>36486.0</td>\n",
       "      <td>DebtCon</td>\n",
       "      <td>Mgr</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>196.01</td>\n",
       "      <td>3.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>42.459</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>156</th>\n",
       "      <td>0</td>\n",
       "      <td>8400</td>\n",
       "      <td>62989.0</td>\n",
       "      <td>76718.0</td>\n",
       "      <td>HomeImp</td>\n",
       "      <td>Mgr</td>\n",
       "      <td>3.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>131.47</td>\n",
       "      <td>0.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>29.200</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>157</th>\n",
       "      <td>0</td>\n",
       "      <td>17400</td>\n",
       "      <td>25859.0</td>\n",
       "      <td>43684.0</td>\n",
       "      <td>DebtCon</td>\n",
       "      <td>Mgr</td>\n",
       "      <td>16.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>95.36</td>\n",
       "      <td>1.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>27.108</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     bad   loan   mortdue     value   reason  job   yoj  derog  delinq  \\\n",
       "153    0  18200   94727.0  136877.0  DebtCon  Mgr  15.0    0.0     0.0   \n",
       "154    0  21700   79240.0   96784.0  DebtCon  Mgr   5.0    0.0     0.0   \n",
       "155    0  34100  241931.0   36486.0  DebtCon  Mgr   1.0    0.0     2.0   \n",
       "156    0   8400   62989.0   76718.0  HomeImp  Mgr   3.0    0.0     2.0   \n",
       "157    0  17400   25859.0   43684.0  DebtCon  Mgr  16.0    1.0     0.0   \n",
       "\n",
       "      clage  ninq  clno  debtinc  \n",
       "153  168.96   2.0  26.0   36.056  \n",
       "154   64.51   6.0  24.0   38.079  \n",
       "155  196.01   3.0  50.0   42.459  \n",
       "156  131.47   0.0  22.0   29.200  \n",
       "157   95.36   1.0  17.0   27.108  "
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "\n",
    "bank = pd.read_csv(\"hmelq.csv\")\n",
    "bank = bank.dropna() \n",
    "bank.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "          count           mean           std       min         25%        50%  \\\n",
      "bad      2018.0       0.088206      0.283665      0.00      0.0000      0.000   \n",
      "loan     2018.0   19010.257681  10768.607345   1700.00  12000.0000  16900.000   \n",
      "mortdue  2018.0   76337.035183  44476.368189   5076.00  49577.2500  67595.000   \n",
      "value    2018.0  106899.244797  54026.707096  23031.00  71287.2500  94514.000   \n",
      "yoj      2018.0       9.003964      7.573012      0.00      3.0000      7.000   \n",
      "derog    2018.0       0.150149      0.580421      0.00      0.0000      0.000   \n",
      "delinq   2018.0       0.264123      0.789816      0.00      0.0000      0.000   \n",
      "clage    2018.0     180.645020     84.420116      0.49    118.2050    176.670   \n",
      "ninq     2018.0       1.037661      1.578414      0.00      0.0000      0.000   \n",
      "clno     2018.0      22.265114      9.457238      0.00     16.0000     21.000   \n",
      "debtinc  2018.0      34.105829      8.438411      4.03     29.2955     34.949   \n",
      "\n",
      "                 75%         max  \n",
      "bad           0.0000       1.000  \n",
      "loan      23600.0000   89900.000  \n",
      "mortdue   92922.7500  399412.000  \n",
      "value    121220.7500  512650.000  \n",
      "yoj          13.0000      41.000  \n",
      "derog         0.0000      10.000  \n",
      "delinq        0.0000      10.000  \n",
      "clage       229.6925    1168.230  \n",
      "ninq          2.0000      13.000  \n",
      "clno         27.0000      64.000  \n",
      "debtinc      38.9520     144.189  \n"
     ]
    }
   ],
   "source": [
    "print(bank.describe().T)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def describeColumn(columnName):\n",
    "    print(\"Mean:\", bank[columnName].mean())\n",
    "    print(\"Median: \", bank[columnName].median())\n",
    "    print(\"Mode: \", bank[columnName].mode())\n",
    "    print(\"Standard Deviation: \",bank[columnName].std())\n",
    "    print(\"Sum: \", bank[columnName].sum())\n",
    "    print(\"Count: \", bank[columnName].count())\n",
    "    print(\"Min: \", bank[columnName].min())\n",
    "    print(\"Max: \", bank[columnName].max())\n",
    "    print(\"Abs: \", bank[columnName].abs())\n",
    "    print(\"Prod: \", bank[columnName].prod())\n",
    "    print(\"Cum sum : \", bank[columnName].cumsum())\n",
    "    print(\"Cum prod: \", bank[columnName].cumprod())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Mean: 19010.25768087215\n",
      "Median:  16900.0\n",
      "Mode:  0    10600\n",
      "1    11000\n",
      "2    12900\n",
      "3    13900\n",
      "4    15200\n",
      "dtype: int64\n",
      "Standard Deviation:  10768.60734485882\n",
      "Sum:  38362700\n",
      "Count:  2018\n",
      "Min:  1700\n",
      "Max:  89900\n",
      "Abs:  153     18200\n",
      "154     21700\n",
      "155     34100\n",
      "156      8400\n",
      "157     17400\n",
      "158     23900\n",
      "159     19700\n",
      "160     17300\n",
      "161     11700\n",
      "162     25300\n",
      "163     11300\n",
      "164     12000\n",
      "165     17400\n",
      "166     15800\n",
      "167      9800\n",
      "168     25900\n",
      "169      8500\n",
      "170     21800\n",
      "171     28400\n",
      "172     14300\n",
      "173     19800\n",
      "174     15200\n",
      "175     12300\n",
      "176     18200\n",
      "177     32400\n",
      "178     19500\n",
      "179     32300\n",
      "180     21800\n",
      "181     17000\n",
      "182     22600\n",
      "        ...  \n",
      "3411    32000\n",
      "3412    27700\n",
      "3413    26800\n",
      "3414    22800\n",
      "3415     9300\n",
      "3416    20300\n",
      "3417    15200\n",
      "3418    23300\n",
      "3419     8200\n",
      "3420    21900\n",
      "3421    20900\n",
      "3422    42400\n",
      "3424    42800\n",
      "3425     5500\n",
      "3426    13900\n",
      "3428    36600\n",
      "3519    28900\n",
      "3520    18200\n",
      "3521     5100\n",
      "3522    27200\n",
      "3523     9000\n",
      "3524     6900\n",
      "3525    16300\n",
      "3542    20500\n",
      "3543    38000\n",
      "3544     7200\n",
      "3546    53400\n",
      "3547    27300\n",
      "3548    37300\n",
      "3549    13300\n",
      "Name: loan, Length: 2018, dtype: int64\n",
      "Prod:  0\n",
      "Cum sum :  153        18200\n",
      "154        39900\n",
      "155        74000\n",
      "156        82400\n",
      "157        99800\n",
      "158       123700\n",
      "159       143400\n",
      "160       160700\n",
      "161       172400\n",
      "162       197700\n",
      "163       209000\n",
      "164       221000\n",
      "165       238400\n",
      "166       254200\n",
      "167       264000\n",
      "168       289900\n",
      "169       298400\n",
      "170       320200\n",
      "171       348600\n",
      "172       362900\n",
      "173       382700\n",
      "174       397900\n",
      "175       410200\n",
      "176       428400\n",
      "177       460800\n",
      "178       480300\n",
      "179       512600\n",
      "180       534400\n",
      "181       551400\n",
      "182       574000\n",
      "          ...   \n",
      "3411    37716500\n",
      "3412    37744200\n",
      "3413    37771000\n",
      "3414    37793800\n",
      "3415    37803100\n",
      "3416    37823400\n",
      "3417    37838600\n",
      "3418    37861900\n",
      "3419    37870100\n",
      "3420    37892000\n",
      "3421    37912900\n",
      "3422    37955300\n",
      "3424    37998100\n",
      "3425    38003600\n",
      "3426    38017500\n",
      "3428    38054100\n",
      "3519    38083000\n",
      "3520    38101200\n",
      "3521    38106300\n",
      "3522    38133500\n",
      "3523    38142500\n",
      "3524    38149400\n",
      "3525    38165700\n",
      "3542    38186200\n",
      "3543    38224200\n",
      "3544    38231400\n",
      "3546    38284800\n",
      "3547    38312100\n",
      "3548    38349400\n",
      "3549    38362700\n",
      "Name: loan, Length: 2018, dtype: int64\n",
      "Cum prod:  153                   18200\n",
      "154               394940000\n",
      "155          13467454000000\n",
      "156      113126613600000000\n",
      "157    -5398539246922022912\n",
      "158    -8559949911743594496\n",
      "159     9121060503909302272\n",
      "160      897911119424782336\n",
      "161    -9084024744491089920\n",
      "162     2158378722728607744\n",
      "163     3083901389240270848\n",
      "164     2648059021889634304\n",
      "165    -3739715246823047168\n",
      "166    -2579631712451428352\n",
      "167    -8351401041912135680\n",
      "168     5234022793888137216\n",
      "169    -4352957738272161792\n",
      "170    -4427179171193552896\n",
      "171     1119144507401568256\n",
      "172    -8007400137464741888\n",
      "173     3242591731706757120\n",
      "174    -2305843009213693952\n",
      "175    -9223372036854775808\n",
      "176                       0\n",
      "177                       0\n",
      "178                       0\n",
      "179                       0\n",
      "180                       0\n",
      "181                       0\n",
      "182                       0\n",
      "               ...         \n",
      "3411                      0\n",
      "3412                      0\n",
      "3413                      0\n",
      "3414                      0\n",
      "3415                      0\n",
      "3416                      0\n",
      "3417                      0\n",
      "3418                      0\n",
      "3419                      0\n",
      "3420                      0\n",
      "3421                      0\n",
      "3422                      0\n",
      "3424                      0\n",
      "3425                      0\n",
      "3426                      0\n",
      "3428                      0\n",
      "3519                      0\n",
      "3520                      0\n",
      "3521                      0\n",
      "3522                      0\n",
      "3523                      0\n",
      "3524                      0\n",
      "3525                      0\n",
      "3542                      0\n",
      "3543                      0\n",
      "3544                      0\n",
      "3546                      0\n",
      "3547                      0\n",
      "3548                      0\n",
      "3549                      0\n",
      "Name: loan, Length: 2018, dtype: int64\n",
      "Mean: 76337.03518334986\n",
      "Median:  67595.0\n",
      "Mode:  0      44097.0\n",
      "1      45279.0\n",
      "2      46406.0\n",
      "3      49367.0\n",
      "4      58606.0\n",
      "5      61213.0\n",
      "6      62933.0\n",
      "7      69697.0\n",
      "8      94727.0\n",
      "9      95351.0\n",
      "10    206126.0\n",
      "dtype: float64\n",
      "Standard Deviation:  44476.368188641594\n",
      "Sum:  154048137.0\n",
      "Count:  2018\n",
      "Min:  5076.0\n",
      "Max:  399412.0\n",
      "Abs:  153      94727.0\n",
      "154      79240.0\n",
      "155     241931.0\n",
      "156      62989.0\n",
      "157      25859.0\n",
      "158     140328.0\n",
      "159      70709.0\n",
      "160      79412.0\n",
      "161      78157.0\n",
      "162      51261.0\n",
      "163      46120.0\n",
      "164      76674.0\n",
      "165      43110.0\n",
      "166      41661.0\n",
      "167      78022.0\n",
      "168     108319.0\n",
      "169      50569.0\n",
      "170     107271.0\n",
      "171      36484.0\n",
      "172      45615.0\n",
      "173      65021.0\n",
      "174     137128.0\n",
      "175      81590.0\n",
      "176      68743.0\n",
      "177     128523.0\n",
      "178     100446.0\n",
      "179     235343.0\n",
      "180      80092.0\n",
      "181      82662.0\n",
      "182      96485.0\n",
      "          ...   \n",
      "3411     33148.0\n",
      "3412     97834.0\n",
      "3413    138216.0\n",
      "3414    129697.0\n",
      "3415     84900.0\n",
      "3416     63803.0\n",
      "3417    105328.0\n",
      "3418     71734.0\n",
      "3419     87979.0\n",
      "3420     52735.0\n",
      "3421     66575.0\n",
      "3422    220757.0\n",
      "3424      9133.0\n",
      "3425     85845.0\n",
      "3426    144263.0\n",
      "3428    397299.0\n",
      "3519     63748.0\n",
      "3520     67809.0\n",
      "3521    129319.0\n",
      "3522     61118.0\n",
      "3523     52257.0\n",
      "3524     53328.0\n",
      "3525     79087.0\n",
      "3542    113963.0\n",
      "3543     94083.0\n",
      "3544     15637.0\n",
      "3546    228236.0\n",
      "3547    106221.0\n",
      "3548     68789.0\n",
      "3549     25698.0\n",
      "Name: mortdue, Length: 2018, dtype: float64\n",
      "Prod:  inf\n",
      "Cum sum :  153         94727.0\n",
      "154        173967.0\n",
      "155        415898.0\n",
      "156        478887.0\n",
      "157        504746.0\n",
      "158        645074.0\n",
      "159        715783.0\n",
      "160        795195.0\n",
      "161        873352.0\n",
      "162        924613.0\n",
      "163        970733.0\n",
      "164       1047407.0\n",
      "165       1090517.0\n",
      "166       1132178.0\n",
      "167       1210200.0\n",
      "168       1318519.0\n",
      "169       1369088.0\n",
      "170       1476359.0\n",
      "171       1512843.0\n",
      "172       1558458.0\n",
      "173       1623479.0\n",
      "174       1760607.0\n",
      "175       1842197.0\n",
      "176       1910940.0\n",
      "177       2039463.0\n",
      "178       2139909.0\n",
      "179       2375252.0\n",
      "180       2455344.0\n",
      "181       2538006.0\n",
      "182       2634491.0\n",
      "           ...     \n",
      "3411    151132746.0\n",
      "3412    151230580.0\n",
      "3413    151368796.0\n",
      "3414    151498493.0\n",
      "3415    151583393.0\n",
      "3416    151647196.0\n",
      "3417    151752524.0\n",
      "3418    151824258.0\n",
      "3419    151912237.0\n",
      "3420    151964972.0\n",
      "3421    152031547.0\n",
      "3422    152252304.0\n",
      "3424    152261437.0\n",
      "3425    152347282.0\n",
      "3426    152491545.0\n",
      "3428    152888844.0\n",
      "3519    152952592.0\n",
      "3520    153020401.0\n",
      "3521    153149720.0\n",
      "3522    153210838.0\n",
      "3523    153263095.0\n",
      "3524    153316423.0\n",
      "3525    153395510.0\n",
      "3542    153509473.0\n",
      "3543    153603556.0\n",
      "3544    153619193.0\n",
      "3546    153847429.0\n",
      "3547    153953650.0\n",
      "3548    154022439.0\n",
      "3549    154048137.0\n",
      "Name: mortdue, Length: 2018, dtype: float64\n",
      "Cum prod:  153      9.472700e+04\n",
      "154      7.506167e+09\n",
      "155      1.815975e+15\n",
      "156      1.143864e+20\n",
      "157      2.957919e+24\n",
      "158      4.150788e+29\n",
      "159      2.934981e+34\n",
      "160      2.330727e+39\n",
      "161      1.821626e+44\n",
      "162      9.337838e+48\n",
      "163      4.306611e+53\n",
      "164      3.302051e+58\n",
      "165      1.423514e+63\n",
      "166      5.930502e+67\n",
      "167      4.627096e+72\n",
      "168      5.012024e+77\n",
      "169      2.534531e+82\n",
      "170      2.718816e+87\n",
      "171      9.919330e+91\n",
      "172      4.524702e+96\n",
      "173     2.942007e+101\n",
      "174     4.034315e+106\n",
      "175     3.291598e+111\n",
      "176     2.262743e+116\n",
      "177     2.908145e+121\n",
      "178     2.921115e+126\n",
      "179     6.874640e+131\n",
      "180     5.506037e+136\n",
      "181     4.551400e+141\n",
      "182     4.391419e+146\n",
      "            ...      \n",
      "3411              inf\n",
      "3412              inf\n",
      "3413              inf\n",
      "3414              inf\n",
      "3415              inf\n",
      "3416              inf\n",
      "3417              inf\n",
      "3418              inf\n",
      "3419              inf\n",
      "3420              inf\n",
      "3421              inf\n",
      "3422              inf\n",
      "3424              inf\n",
      "3425              inf\n",
      "3426              inf\n",
      "3428              inf\n",
      "3519              inf\n",
      "3520              inf\n",
      "3521              inf\n",
      "3522              inf\n",
      "3523              inf\n",
      "3524              inf\n",
      "3525              inf\n",
      "3542              inf\n",
      "3543              inf\n",
      "3544              inf\n",
      "3546              inf\n",
      "3547              inf\n",
      "3548              inf\n",
      "3549              inf\n",
      "Name: mortdue, Length: 2018, dtype: float64\n",
      "Mean: 106899.24479682854\n",
      "Median:  94514.0\n",
      "Mode:  0     57280.0\n",
      "1     65481.0\n",
      "2     85424.0\n",
      "3     98088.0\n",
      "4    100332.0\n",
      "5    102287.0\n",
      "6    106233.0\n",
      "7    114494.0\n",
      "8    114906.0\n",
      "9    116160.0\n",
      "dtype: float64\n",
      "Standard Deviation:  54026.70709615272\n",
      "Sum:  215722676.0\n",
      "Count:  2018\n",
      "Min:  23031.0\n",
      "Max:  512650.0\n",
      "Abs:  153     136877.0\n",
      "154      96784.0\n",
      "155      36486.0\n",
      "156      76718.0\n",
      "157      43684.0\n",
      "158     183533.0\n",
      "159      97394.0\n",
      "160     114906.0\n",
      "161     103016.0\n",
      "162      87025.0\n",
      "163      52911.0\n",
      "164      97936.0\n",
      "165      73614.0\n",
      "166      80457.0\n",
      "167     103466.0\n",
      "168     135884.0\n",
      "169      60940.0\n",
      "170     161219.0\n",
      "171      73478.0\n",
      "172      80818.0\n",
      "173      95188.0\n",
      "174     184660.0\n",
      "175      90050.0\n",
      "176      92112.0\n",
      "177     179390.0\n",
      "178     121516.0\n",
      "179      45270.0\n",
      "180      98561.0\n",
      "181     111357.0\n",
      "182     120470.0\n",
      "          ...   \n",
      "3411     69965.0\n",
      "3412    132396.0\n",
      "3413    170052.0\n",
      "3414    167475.0\n",
      "3415    105798.0\n",
      "3416     89079.0\n",
      "3417    113931.0\n",
      "3418    100072.0\n",
      "3419     94854.0\n",
      "3420     75372.0\n",
      "3421     79416.0\n",
      "3422    297444.0\n",
      "3424     44083.0\n",
      "3425    122205.0\n",
      "3426    167164.0\n",
      "3428    512650.0\n",
      "3519     36392.0\n",
      "3520     89339.0\n",
      "3521    141875.0\n",
      "3522     35184.0\n",
      "3523     58275.0\n",
      "3524     58088.0\n",
      "3525    102921.0\n",
      "3542    131135.0\n",
      "3543    132791.0\n",
      "3544     39660.0\n",
      "3546    305514.0\n",
      "3547    137727.0\n",
      "3548     86172.0\n",
      "3549     58760.0\n",
      "Name: value, Length: 2018, dtype: float64\n",
      "Prod:  inf\n",
      "Cum sum :  153        136877.0\n",
      "154        233661.0\n",
      "155        270147.0\n",
      "156        346865.0\n",
      "157        390549.0\n",
      "158        574082.0\n",
      "159        671476.0\n",
      "160        786382.0\n",
      "161        889398.0\n",
      "162        976423.0\n",
      "163       1029334.0\n",
      "164       1127270.0\n",
      "165       1200884.0\n",
      "166       1281341.0\n",
      "167       1384807.0\n",
      "168       1520691.0\n",
      "169       1581631.0\n",
      "170       1742850.0\n",
      "171       1816328.0\n",
      "172       1897146.0\n",
      "173       1992334.0\n",
      "174       2176994.0\n",
      "175       2267044.0\n",
      "176       2359156.0\n",
      "177       2538546.0\n",
      "178       2660062.0\n",
      "179       2705332.0\n",
      "180       2803893.0\n",
      "181       2915250.0\n",
      "182       3035720.0\n",
      "           ...     \n",
      "3411    212036852.0\n",
      "3412    212169248.0\n",
      "3413    212339300.0\n",
      "3414    212506775.0\n",
      "3415    212612573.0\n",
      "3416    212701652.0\n",
      "3417    212815583.0\n",
      "3418    212915655.0\n",
      "3419    213010509.0\n",
      "3420    213085881.0\n",
      "3421    213165297.0\n",
      "3422    213462741.0\n",
      "3424    213506824.0\n",
      "3425    213629029.0\n",
      "3426    213796193.0\n",
      "3428    214308843.0\n",
      "3519    214345235.0\n",
      "3520    214434574.0\n",
      "3521    214576449.0\n",
      "3522    214611633.0\n",
      "3523    214669908.0\n",
      "3524    214727996.0\n",
      "3525    214830917.0\n",
      "3542    214962052.0\n",
      "3543    215094843.0\n",
      "3544    215134503.0\n",
      "3546    215440017.0\n",
      "3547    215577744.0\n",
      "3548    215663916.0\n",
      "3549    215722676.0\n",
      "Name: value, Length: 2018, dtype: float64\n",
      "Cum prod:  153      1.368770e+05\n",
      "154      1.324750e+10\n",
      "155      4.833484e+14\n",
      "156      3.708152e+19\n",
      "157      1.619869e+24\n",
      "158      2.972995e+29\n",
      "159      2.895518e+34\n",
      "160      3.327124e+39\n",
      "161      3.427471e+44\n",
      "162      2.982756e+49\n",
      "163      1.578206e+54\n",
      "164      1.545632e+59\n",
      "165      1.137802e+64\n",
      "166      9.154410e+68\n",
      "167      9.471701e+73\n",
      "168      1.287053e+79\n",
      "169      7.843299e+83\n",
      "170      1.264489e+89\n",
      "171      9.291211e+93\n",
      "172      7.508971e+98\n",
      "173     7.147639e+103\n",
      "174     1.319883e+109\n",
      "175     1.188555e+114\n",
      "176     1.094802e+119\n",
      "177     1.963964e+124\n",
      "178     2.386531e+129\n",
      "179     1.080383e+134\n",
      "180     1.064836e+139\n",
      "181     1.185769e+144\n",
      "182     1.428496e+149\n",
      "            ...      \n",
      "3411              inf\n",
      "3412              inf\n",
      "3413              inf\n",
      "3414              inf\n",
      "3415              inf\n",
      "3416              inf\n",
      "3417              inf\n",
      "3418              inf\n",
      "3419              inf\n",
      "3420              inf\n",
      "3421              inf\n",
      "3422              inf\n",
      "3424              inf\n",
      "3425              inf\n",
      "3426              inf\n",
      "3428              inf\n",
      "3519              inf\n",
      "3520              inf\n",
      "3521              inf\n",
      "3522              inf\n",
      "3523              inf\n",
      "3524              inf\n",
      "3525              inf\n",
      "3542              inf\n",
      "3543              inf\n",
      "3544              inf\n",
      "3546              inf\n",
      "3547              inf\n",
      "3548              inf\n",
      "3549              inf\n",
      "Name: value, Length: 2018, dtype: float64\n",
      "Mean: 9.00396432111001\n",
      "Median:  7.0\n",
      "Mode:  0    0.0\n",
      "dtype: float64\n",
      "Standard Deviation:  7.573011610414845\n",
      "Sum:  18170.0\n",
      "Count:  2018\n",
      "Min:  0.0\n",
      "Max:  41.0\n",
      "Abs:  153     15.0\n",
      "154      5.0\n",
      "155      1.0\n",
      "156      3.0\n",
      "157     16.0\n",
      "158      8.0\n",
      "159      0.0\n",
      "160     17.0\n",
      "161      1.0\n",
      "162      4.0\n",
      "163     19.0\n",
      "164      5.0\n",
      "165      1.0\n",
      "166      3.0\n",
      "167      1.0\n",
      "168      8.0\n",
      "169      6.0\n",
      "170     15.0\n",
      "171     13.0\n",
      "172      4.0\n",
      "173      1.0\n",
      "174     11.0\n",
      "175      6.0\n",
      "176      0.0\n",
      "177     13.0\n",
      "178      0.0\n",
      "179      3.0\n",
      "180      5.0\n",
      "181     17.0\n",
      "182      9.0\n",
      "        ... \n",
      "3411    13.0\n",
      "3412    18.0\n",
      "3413    10.0\n",
      "3414     5.0\n",
      "3415     2.0\n",
      "3416    21.0\n",
      "3417     3.0\n",
      "3418     6.0\n",
      "3419    10.0\n",
      "3420     5.0\n",
      "3421    13.0\n",
      "3422    30.0\n",
      "3424     5.0\n",
      "3425     5.0\n",
      "3426     8.0\n",
      "3428     4.0\n",
      "3519     9.0\n",
      "3520     0.0\n",
      "3521     7.0\n",
      "3522    10.0\n",
      "3523     5.0\n",
      "3524     2.0\n",
      "3525    23.0\n",
      "3542     3.0\n",
      "3543    10.0\n",
      "3544     6.0\n",
      "3546     6.0\n",
      "3547     3.0\n",
      "3548     2.0\n",
      "3549    15.0\n",
      "Name: yoj, Length: 2018, dtype: float64\n",
      "Prod:  0.0\n",
      "Cum sum :  153        15.0\n",
      "154        20.0\n",
      "155        21.0\n",
      "156        24.0\n",
      "157        40.0\n",
      "158        48.0\n",
      "159        48.0\n",
      "160        65.0\n",
      "161        66.0\n",
      "162        70.0\n",
      "163        89.0\n",
      "164        94.0\n",
      "165        95.0\n",
      "166        98.0\n",
      "167        99.0\n",
      "168       107.0\n",
      "169       113.0\n",
      "170       128.0\n",
      "171       141.0\n",
      "172       145.0\n",
      "173       146.0\n",
      "174       157.0\n",
      "175       163.0\n",
      "176       163.0\n",
      "177       176.0\n",
      "178       176.0\n",
      "179       179.0\n",
      "180       184.0\n",
      "181       201.0\n",
      "182       210.0\n",
      "         ...   \n",
      "3411    17924.0\n",
      "3412    17942.0\n",
      "3413    17952.0\n",
      "3414    17957.0\n",
      "3415    17959.0\n",
      "3416    17980.0\n",
      "3417    17983.0\n",
      "3418    17989.0\n",
      "3419    17999.0\n",
      "3420    18004.0\n",
      "3421    18017.0\n",
      "3422    18047.0\n",
      "3424    18052.0\n",
      "3425    18057.0\n",
      "3426    18065.0\n",
      "3428    18069.0\n",
      "3519    18078.0\n",
      "3520    18078.0\n",
      "3521    18085.0\n",
      "3522    18095.0\n",
      "3523    18100.0\n",
      "3524    18102.0\n",
      "3525    18125.0\n",
      "3542    18128.0\n",
      "3543    18138.0\n",
      "3544    18144.0\n",
      "3546    18150.0\n",
      "3547    18153.0\n",
      "3548    18155.0\n",
      "3549    18170.0\n",
      "Name: yoj, Length: 2018, dtype: float64\n",
      "Cum prod:  153        15.0\n",
      "154        75.0\n",
      "155        75.0\n",
      "156       225.0\n",
      "157      3600.0\n",
      "158     28800.0\n",
      "159         0.0\n",
      "160         0.0\n",
      "161         0.0\n",
      "162         0.0\n",
      "163         0.0\n",
      "164         0.0\n",
      "165         0.0\n",
      "166         0.0\n",
      "167         0.0\n",
      "168         0.0\n",
      "169         0.0\n",
      "170         0.0\n",
      "171         0.0\n",
      "172         0.0\n",
      "173         0.0\n",
      "174         0.0\n",
      "175         0.0\n",
      "176         0.0\n",
      "177         0.0\n",
      "178         0.0\n",
      "179         0.0\n",
      "180         0.0\n",
      "181         0.0\n",
      "182         0.0\n",
      "         ...   \n",
      "3411        0.0\n",
      "3412        0.0\n",
      "3413        0.0\n",
      "3414        0.0\n",
      "3415        0.0\n",
      "3416        0.0\n",
      "3417        0.0\n",
      "3418        0.0\n",
      "3419        0.0\n",
      "3420        0.0\n",
      "3421        0.0\n",
      "3422        0.0\n",
      "3424        0.0\n",
      "3425        0.0\n",
      "3426        0.0\n",
      "3428        0.0\n",
      "3519        0.0\n",
      "3520        0.0\n",
      "3521        0.0\n",
      "3522        0.0\n",
      "3523        0.0\n",
      "3524        0.0\n",
      "3525        0.0\n",
      "3542        0.0\n",
      "3543        0.0\n",
      "3544        0.0\n",
      "3546        0.0\n",
      "3547        0.0\n",
      "3548        0.0\n",
      "3549        0.0\n",
      "Name: yoj, Length: 2018, dtype: float64\n",
      "Mean: 180.6450198216057\n",
      "Median:  176.67000000000002\n",
      "Mode:  0    219.09\n",
      "dtype: float64\n",
      "Standard Deviation:  84.42011595239947\n",
      "Sum:  364541.65\n",
      "Count:  2018\n",
      "Min:  0.49\n",
      "Max:  1168.23\n",
      "Abs:  153     168.96\n",
      "154      64.51\n",
      "155     196.01\n",
      "156     131.47\n",
      "157      95.36\n",
      "158     322.79\n",
      "159      94.86\n",
      "160     247.36\n",
      "161      69.81\n",
      "162     249.84\n",
      "163     229.44\n",
      "164     147.91\n",
      "165     114.29\n",
      "166     137.32\n",
      "167      83.84\n",
      "168     188.38\n",
      "169     189.95\n",
      "170     187.30\n",
      "171     100.98\n",
      "172     115.14\n",
      "173      73.08\n",
      "174     254.36\n",
      "175     136.75\n",
      "176      95.34\n",
      "177     150.08\n",
      "178     306.11\n",
      "179     210.92\n",
      "180      73.32\n",
      "181     252.31\n",
      "182     112.84\n",
      "         ...  \n",
      "3411    112.51\n",
      "3412    303.95\n",
      "3413    199.74\n",
      "3414    202.60\n",
      "3415     41.81\n",
      "3416    323.81\n",
      "3417    138.13\n",
      "3418    122.80\n",
      "3419    108.48\n",
      "3420    217.02\n",
      "3421    107.27\n",
      "3422    237.59\n",
      "3424    333.91\n",
      "3425    172.21\n",
      "3426    111.59\n",
      "3428    197.52\n",
      "3519    124.25\n",
      "3520    102.46\n",
      "3521    157.86\n",
      "3522    134.65\n",
      "3523     80.25\n",
      "3524    101.53\n",
      "3525    256.25\n",
      "3542    243.68\n",
      "3543    212.99\n",
      "3544    354.24\n",
      "3546     11.15\n",
      "3547    300.57\n",
      "3548    235.64\n",
      "3549    147.91\n",
      "Name: clage, Length: 2018, dtype: float64\n",
      "Prod:  inf\n",
      "Cum sum :  153        168.96\n",
      "154        233.47\n",
      "155        429.48\n",
      "156        560.95\n",
      "157        656.31\n",
      "158        979.10\n",
      "159       1073.96\n",
      "160       1321.32\n",
      "161       1391.13\n",
      "162       1640.97\n",
      "163       1870.41\n",
      "164       2018.32\n",
      "165       2132.61\n",
      "166       2269.93\n",
      "167       2353.77\n",
      "168       2542.15\n",
      "169       2732.10\n",
      "170       2919.40\n",
      "171       3020.38\n",
      "172       3135.52\n",
      "173       3208.60\n",
      "174       3462.96\n",
      "175       3599.71\n",
      "176       3695.05\n",
      "177       3845.13\n",
      "178       4151.24\n",
      "179       4362.16\n",
      "180       4435.48\n",
      "181       4687.79\n",
      "182       4800.63\n",
      "          ...    \n",
      "3411    359259.79\n",
      "3412    359563.74\n",
      "3413    359763.48\n",
      "3414    359966.08\n",
      "3415    360007.89\n",
      "3416    360331.70\n",
      "3417    360469.83\n",
      "3418    360592.63\n",
      "3419    360701.11\n",
      "3420    360918.13\n",
      "3421    361025.40\n",
      "3422    361262.99\n",
      "3424    361596.90\n",
      "3425    361769.11\n",
      "3426    361880.70\n",
      "3428    362078.22\n",
      "3519    362202.47\n",
      "3520    362304.93\n",
      "3521    362462.79\n",
      "3522    362597.44\n",
      "3523    362677.69\n",
      "3524    362779.22\n",
      "3525    363035.47\n",
      "3542    363279.15\n",
      "3543    363492.14\n",
      "3544    363846.38\n",
      "3546    363857.53\n",
      "3547    364158.10\n",
      "3548    364393.74\n",
      "3549    364541.65\n",
      "Name: clage, Length: 2018, dtype: float64\n",
      "Cum prod:  153     1.689600e+02\n",
      "154     1.089961e+04\n",
      "155     2.136432e+06\n",
      "156     2.808768e+08\n",
      "157     2.678441e+10\n",
      "158     8.645740e+12\n",
      "159     8.201349e+14\n",
      "160     2.028686e+17\n",
      "161     1.416225e+19\n",
      "162     3.538298e+21\n",
      "163     8.118270e+23\n",
      "164     1.200773e+26\n",
      "165     1.372364e+28\n",
      "166     1.884530e+30\n",
      "167     1.579990e+32\n",
      "168     2.976385e+34\n",
      "169     5.653643e+36\n",
      "170     1.058927e+39\n",
      "171     1.069305e+41\n",
      "172     1.231198e+43\n",
      "173     8.997592e+44\n",
      "174     2.288628e+47\n",
      "175     3.129698e+49\n",
      "176     2.983854e+51\n",
      "177     4.478169e+53\n",
      "178     1.370812e+56\n",
      "179     2.891317e+58\n",
      "180     2.119914e+60\n",
      "181     5.348754e+62\n",
      "182     6.035534e+64\n",
      "            ...     \n",
      "3411             inf\n",
      "3412             inf\n",
      "3413             inf\n",
      "3414             inf\n",
      "3415             inf\n",
      "3416             inf\n",
      "3417             inf\n",
      "3418             inf\n",
      "3419             inf\n",
      "3420             inf\n",
      "3421             inf\n",
      "3422             inf\n",
      "3424             inf\n",
      "3425             inf\n",
      "3426             inf\n",
      "3428             inf\n",
      "3519             inf\n",
      "3520             inf\n",
      "3521             inf\n",
      "3522             inf\n",
      "3523             inf\n",
      "3524             inf\n",
      "3525             inf\n",
      "3542             inf\n",
      "3543             inf\n",
      "3544             inf\n",
      "3546             inf\n",
      "3547             inf\n",
      "3548             inf\n",
      "3549             inf\n",
      "Name: clage, Length: 2018, dtype: float64\n",
      "Mean: 1.037661050545094\n",
      "Median:  0.0\n",
      "Mode:  0    0.0\n",
      "dtype: float64\n",
      "Standard Deviation:  1.5784144093595993\n",
      "Sum:  2094.0\n",
      "Count:  2018\n",
      "Min:  0.0\n",
      "Max:  13.0\n",
      "Abs:  153     2.0\n",
      "154     6.0\n",
      "155     3.0\n",
      "156     0.0\n",
      "157     1.0\n",
      "158     0.0\n",
      "159     0.0\n",
      "160     4.0\n",
      "161     2.0\n",
      "162     2.0\n",
      "163     1.0\n",
      "164     2.0\n",
      "165     1.0\n",
      "166     0.0\n",
      "167     2.0\n",
      "168     5.0\n",
      "169     0.0\n",
      "170     0.0\n",
      "171     6.0\n",
      "172     0.0\n",
      "173     0.0\n",
      "174     0.0\n",
      "175     1.0\n",
      "176     0.0\n",
      "177     1.0\n",
      "178     1.0\n",
      "179     3.0\n",
      "180     7.0\n",
      "181     4.0\n",
      "182     2.0\n",
      "       ... \n",
      "3411    3.0\n",
      "3412    2.0\n",
      "3413    1.0\n",
      "3414    1.0\n",
      "3415    1.0\n",
      "3416    4.0\n",
      "3417    2.0\n",
      "3418    0.0\n",
      "3419    3.0\n",
      "3420    6.0\n",
      "3421    7.0\n",
      "3422    1.0\n",
      "3424    5.0\n",
      "3425    1.0\n",
      "3426    5.0\n",
      "3428    0.0\n",
      "3519    6.0\n",
      "3520    2.0\n",
      "3521    0.0\n",
      "3522    7.0\n",
      "3523    0.0\n",
      "3524    0.0\n",
      "3525    0.0\n",
      "3542    1.0\n",
      "3543    0.0\n",
      "3544    3.0\n",
      "3546    0.0\n",
      "3547    3.0\n",
      "3548    2.0\n",
      "3549    1.0\n",
      "Name: ninq, Length: 2018, dtype: float64\n",
      "Prod:  0.0\n",
      "Cum sum :  153        2.0\n",
      "154        8.0\n",
      "155       11.0\n",
      "156       11.0\n",
      "157       12.0\n",
      "158       12.0\n",
      "159       12.0\n",
      "160       16.0\n",
      "161       18.0\n",
      "162       20.0\n",
      "163       21.0\n",
      "164       23.0\n",
      "165       24.0\n",
      "166       24.0\n",
      "167       26.0\n",
      "168       31.0\n",
      "169       31.0\n",
      "170       31.0\n",
      "171       37.0\n",
      "172       37.0\n",
      "173       37.0\n",
      "174       37.0\n",
      "175       38.0\n",
      "176       38.0\n",
      "177       39.0\n",
      "178       40.0\n",
      "179       43.0\n",
      "180       50.0\n",
      "181       54.0\n",
      "182       56.0\n",
      "         ...  \n",
      "3411    2030.0\n",
      "3412    2032.0\n",
      "3413    2033.0\n",
      "3414    2034.0\n",
      "3415    2035.0\n",
      "3416    2039.0\n",
      "3417    2041.0\n",
      "3418    2041.0\n",
      "3419    2044.0\n",
      "3420    2050.0\n",
      "3421    2057.0\n",
      "3422    2058.0\n",
      "3424    2063.0\n",
      "3425    2064.0\n",
      "3426    2069.0\n",
      "3428    2069.0\n",
      "3519    2075.0\n",
      "3520    2077.0\n",
      "3521    2077.0\n",
      "3522    2084.0\n",
      "3523    2084.0\n",
      "3524    2084.0\n",
      "3525    2084.0\n",
      "3542    2085.0\n",
      "3543    2085.0\n",
      "3544    2088.0\n",
      "3546    2088.0\n",
      "3547    2091.0\n",
      "3548    2093.0\n",
      "3549    2094.0\n",
      "Name: ninq, Length: 2018, dtype: float64\n",
      "Cum prod:  153      2.0\n",
      "154     12.0\n",
      "155     36.0\n",
      "156      0.0\n",
      "157      0.0\n",
      "158      0.0\n",
      "159      0.0\n",
      "160      0.0\n",
      "161      0.0\n",
      "162      0.0\n",
      "163      0.0\n",
      "164      0.0\n",
      "165      0.0\n",
      "166      0.0\n",
      "167      0.0\n",
      "168      0.0\n",
      "169      0.0\n",
      "170      0.0\n",
      "171      0.0\n",
      "172      0.0\n",
      "173      0.0\n",
      "174      0.0\n",
      "175      0.0\n",
      "176      0.0\n",
      "177      0.0\n",
      "178      0.0\n",
      "179      0.0\n",
      "180      0.0\n",
      "181      0.0\n",
      "182      0.0\n",
      "        ... \n",
      "3411     0.0\n",
      "3412     0.0\n",
      "3413     0.0\n",
      "3414     0.0\n",
      "3415     0.0\n",
      "3416     0.0\n",
      "3417     0.0\n",
      "3418     0.0\n",
      "3419     0.0\n",
      "3420     0.0\n",
      "3421     0.0\n",
      "3422     0.0\n",
      "3424     0.0\n",
      "3425     0.0\n",
      "3426     0.0\n",
      "3428     0.0\n",
      "3519     0.0\n",
      "3520     0.0\n",
      "3521     0.0\n",
      "3522     0.0\n",
      "3523     0.0\n",
      "3524     0.0\n",
      "3525     0.0\n",
      "3542     0.0\n",
      "3543     0.0\n",
      "3544     0.0\n",
      "3546     0.0\n",
      "3547     0.0\n",
      "3548     0.0\n",
      "3549     0.0\n",
      "Name: ninq, Length: 2018, dtype: float64\n",
      "Mean: 22.265113974231912\n",
      "Median:  21.0\n",
      "Mode:  0    19.0\n",
      "dtype: float64\n",
      "Standard Deviation:  9.45723783975282\n",
      "Sum:  44931.0\n",
      "Count:  2018\n",
      "Min:  0.0\n",
      "Max:  64.0\n",
      "Abs:  153     26.0\n",
      "154     24.0\n",
      "155     50.0\n",
      "156     22.0\n",
      "157     17.0\n",
      "158     17.0\n",
      "159     17.0\n",
      "160     30.0\n",
      "161     24.0\n",
      "162     25.0\n",
      "163     17.0\n",
      "164     27.0\n",
      "165      8.0\n",
      "166      8.0\n",
      "167     24.0\n",
      "168     18.0\n",
      "169     22.0\n",
      "170     26.0\n",
      "171     16.0\n",
      "172      9.0\n",
      "173     16.0\n",
      "174     21.0\n",
      "175     26.0\n",
      "176     16.0\n",
      "177     28.0\n",
      "178     20.0\n",
      "179     49.0\n",
      "180     24.0\n",
      "181     30.0\n",
      "182     24.0\n",
      "        ... \n",
      "3411    27.0\n",
      "3412    26.0\n",
      "3413    29.0\n",
      "3414    18.0\n",
      "3415     3.0\n",
      "3416    28.0\n",
      "3417    32.0\n",
      "3418    21.0\n",
      "3419    19.0\n",
      "3420    23.0\n",
      "3421    28.0\n",
      "3422    25.0\n",
      "3424    15.0\n",
      "3425    26.0\n",
      "3426    38.0\n",
      "3428    21.0\n",
      "3519    23.0\n",
      "3520    19.0\n",
      "3521    36.0\n",
      "3522    23.0\n",
      "3523    21.0\n",
      "3524    21.0\n",
      "3525    45.0\n",
      "3542    26.0\n",
      "3543    64.0\n",
      "3544    42.0\n",
      "3546     2.0\n",
      "3547    45.0\n",
      "3548    25.0\n",
      "3549    10.0\n",
      "Name: clno, Length: 2018, dtype: float64\n",
      "Prod:  nan\n",
      "Cum sum :  153        26.0\n",
      "154        50.0\n",
      "155       100.0\n",
      "156       122.0\n",
      "157       139.0\n",
      "158       156.0\n",
      "159       173.0\n",
      "160       203.0\n",
      "161       227.0\n",
      "162       252.0\n",
      "163       269.0\n",
      "164       296.0\n",
      "165       304.0\n",
      "166       312.0\n",
      "167       336.0\n",
      "168       354.0\n",
      "169       376.0\n",
      "170       402.0\n",
      "171       418.0\n",
      "172       427.0\n",
      "173       443.0\n",
      "174       464.0\n",
      "175       490.0\n",
      "176       506.0\n",
      "177       534.0\n",
      "178       554.0\n",
      "179       603.0\n",
      "180       627.0\n",
      "181       657.0\n",
      "182       681.0\n",
      "         ...   \n",
      "3411    44177.0\n",
      "3412    44203.0\n",
      "3413    44232.0\n",
      "3414    44250.0\n",
      "3415    44253.0\n",
      "3416    44281.0\n",
      "3417    44313.0\n",
      "3418    44334.0\n",
      "3419    44353.0\n",
      "3420    44376.0\n",
      "3421    44404.0\n",
      "3422    44429.0\n",
      "3424    44444.0\n",
      "3425    44470.0\n",
      "3426    44508.0\n",
      "3428    44529.0\n",
      "3519    44552.0\n",
      "3520    44571.0\n",
      "3521    44607.0\n",
      "3522    44630.0\n",
      "3523    44651.0\n",
      "3524    44672.0\n",
      "3525    44717.0\n",
      "3542    44743.0\n",
      "3543    44807.0\n",
      "3544    44849.0\n",
      "3546    44851.0\n",
      "3547    44896.0\n",
      "3548    44921.0\n",
      "3549    44931.0\n",
      "Name: clno, Length: 2018, dtype: float64\n",
      "Cum prod:  153     2.600000e+01\n",
      "154     6.240000e+02\n",
      "155     3.120000e+04\n",
      "156     6.864000e+05\n",
      "157     1.166880e+07\n",
      "158     1.983696e+08\n",
      "159     3.372283e+09\n",
      "160     1.011685e+11\n",
      "161     2.428044e+12\n",
      "162     6.070110e+13\n",
      "163     1.031919e+15\n",
      "164     2.786180e+16\n",
      "165     2.228944e+17\n",
      "166     1.783155e+18\n",
      "167     4.279573e+19\n",
      "168     7.703232e+20\n",
      "169     1.694711e+22\n",
      "170     4.406248e+23\n",
      "171     7.049997e+24\n",
      "172     6.344998e+25\n",
      "173     1.015200e+27\n",
      "174     2.131919e+28\n",
      "175     5.542990e+29\n",
      "176     8.868784e+30\n",
      "177     2.483260e+32\n",
      "178     4.966519e+33\n",
      "179     2.433594e+35\n",
      "180     5.840626e+36\n",
      "181     1.752188e+38\n",
      "182     4.205251e+39\n",
      "            ...     \n",
      "3411             NaN\n",
      "3412             NaN\n",
      "3413             NaN\n",
      "3414             NaN\n",
      "3415             NaN\n",
      "3416             NaN\n",
      "3417             NaN\n",
      "3418             NaN\n",
      "3419             NaN\n",
      "3420             NaN\n",
      "3421             NaN\n",
      "3422             NaN\n",
      "3424             NaN\n",
      "3425             NaN\n",
      "3426             NaN\n",
      "3428             NaN\n",
      "3519             NaN\n",
      "3520             NaN\n",
      "3521             NaN\n",
      "3522             NaN\n",
      "3523             NaN\n",
      "3524             NaN\n",
      "3525             NaN\n",
      "3542             NaN\n",
      "3543             NaN\n",
      "3544             NaN\n",
      "3546             NaN\n",
      "3547             NaN\n",
      "3548             NaN\n",
      "3549             NaN\n",
      "Name: clno, Length: 2018, dtype: float64\n",
      "Mean: 34.10582854311199\n",
      "Median:  34.949\n",
      "Mode:  0    28.172\n",
      "1    33.691\n",
      "2    40.206\n",
      "3    42.930\n",
      "dtype: float64\n",
      "Standard Deviation:  8.438411321479823\n",
      "Sum:  68825.56199999999\n",
      "Count:  2018\n",
      "Min:  4.03\n",
      "Max:  144.189\n",
      "Abs:  153      36.056\n",
      "154      38.079\n",
      "155      42.459\n",
      "156      29.200\n",
      "157      27.108\n",
      "158      45.413\n",
      "159      25.557\n",
      "160      36.536\n",
      "161      38.448\n",
      "162      31.659\n",
      "163      36.335\n",
      "164      41.226\n",
      "165      26.379\n",
      "166      27.434\n",
      "167      37.435\n",
      "168      40.505\n",
      "169      29.202\n",
      "170      31.557\n",
      "171      34.764\n",
      "172      28.738\n",
      "173      27.443\n",
      "174      28.815\n",
      "175      41.733\n",
      "176      26.693\n",
      "177      43.144\n",
      "178      43.022\n",
      "179      41.265\n",
      "180      37.138\n",
      "181      37.316\n",
      "182      31.866\n",
      "         ...   \n",
      "3411     36.022\n",
      "3412      6.887\n",
      "3413     76.421\n",
      "3414     22.044\n",
      "3415     26.957\n",
      "3416     45.876\n",
      "3417     61.891\n",
      "3418     34.530\n",
      "3419     41.271\n",
      "3420     34.016\n",
      "3421     28.382\n",
      "3422     44.471\n",
      "3424     13.418\n",
      "3425     32.677\n",
      "3426     91.613\n",
      "3428     43.458\n",
      "3519     44.579\n",
      "3520     38.345\n",
      "3521     25.651\n",
      "3522     45.255\n",
      "3523     34.145\n",
      "3524     35.764\n",
      "3525     36.987\n",
      "3542    143.950\n",
      "3543     62.802\n",
      "3544     33.624\n",
      "3546     34.558\n",
      "3547     32.560\n",
      "3548    114.051\n",
      "3549     30.962\n",
      "Name: debtinc, Length: 2018, dtype: float64\n",
      "Prod:  inf\n",
      "Cum sum :  153        36.056\n",
      "154        74.135\n",
      "155       116.594\n",
      "156       145.794\n",
      "157       172.902\n",
      "158       218.315\n",
      "159       243.872\n",
      "160       280.408\n",
      "161       318.856\n",
      "162       350.515\n",
      "163       386.850\n",
      "164       428.076\n",
      "165       454.455\n",
      "166       481.889\n",
      "167       519.324\n",
      "168       559.829\n",
      "169       589.031\n",
      "170       620.588\n",
      "171       655.352\n",
      "172       684.090\n",
      "173       711.533\n",
      "174       740.348\n",
      "175       782.081\n",
      "176       808.774\n",
      "177       851.918\n",
      "178       894.940\n",
      "179       936.205\n",
      "180       973.343\n",
      "181      1010.659\n",
      "182      1042.525\n",
      "          ...    \n",
      "3411    67508.417\n",
      "3412    67515.304\n",
      "3413    67591.725\n",
      "3414    67613.769\n",
      "3415    67640.726\n",
      "3416    67686.602\n",
      "3417    67748.493\n",
      "3418    67783.023\n",
      "3419    67824.294\n",
      "3420    67858.310\n",
      "3421    67886.692\n",
      "3422    67931.163\n",
      "3424    67944.581\n",
      "3425    67977.258\n",
      "3426    68068.871\n",
      "3428    68112.329\n",
      "3519    68156.908\n",
      "3520    68195.253\n",
      "3521    68220.904\n",
      "3522    68266.159\n",
      "3523    68300.304\n",
      "3524    68336.068\n",
      "3525    68373.055\n",
      "3542    68517.005\n",
      "3543    68579.807\n",
      "3544    68613.431\n",
      "3546    68647.989\n",
      "3547    68680.549\n",
      "3548    68794.600\n",
      "3549    68825.562\n",
      "Name: debtinc, Length: 2018, dtype: float64\n",
      "Cum prod:  153     3.605600e+01\n",
      "154     1.372976e+03\n",
      "155     5.829521e+04\n",
      "156     1.702220e+06\n",
      "157     4.614378e+07\n",
      "158     2.095527e+09\n",
      "159     5.355540e+10\n",
      "160     1.956700e+12\n",
      "161     7.523120e+13\n",
      "162     2.381745e+15\n",
      "163     8.654069e+16\n",
      "164     3.567726e+18\n",
      "165     9.411306e+19\n",
      "166     2.581898e+21\n",
      "167     9.665333e+22\n",
      "168     3.914943e+24\n",
      "169     1.143242e+26\n",
      "170     3.607728e+27\n",
      "171     1.254191e+29\n",
      "172     3.604293e+30\n",
      "173     9.891261e+31\n",
      "174     2.850167e+33\n",
      "175     1.189460e+35\n",
      "176     3.175026e+36\n",
      "177     1.369833e+38\n",
      "178     5.893296e+39\n",
      "179     2.431869e+41\n",
      "180     9.031474e+42\n",
      "181     3.370185e+44\n",
      "182     1.073943e+46\n",
      "            ...     \n",
      "3411             inf\n",
      "3412             inf\n",
      "3413             inf\n",
      "3414             inf\n",
      "3415             inf\n",
      "3416             inf\n",
      "3417             inf\n",
      "3418             inf\n",
      "3419             inf\n",
      "3420             inf\n",
      "3421             inf\n",
      "3422             inf\n",
      "3424             inf\n",
      "3425             inf\n",
      "3426             inf\n",
      "3428             inf\n",
      "3519             inf\n",
      "3520             inf\n",
      "3521             inf\n",
      "3522             inf\n",
      "3523             inf\n",
      "3524             inf\n",
      "3525             inf\n",
      "3542             inf\n",
      "3543             inf\n",
      "3544             inf\n",
      "3546             inf\n",
      "3547             inf\n",
      "3548             inf\n",
      "3549             inf\n",
      "Name: debtinc, Length: 2018, dtype: float64\n"
     ]
    }
   ],
   "source": [
    "colums=['loan','mortdue','value','yoj','clage','ninq','clno','debtinc']\n",
    "for c in colums:\n",
    "    describeColumn(c)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEGCAYAAACkQqisAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAWSklEQVR4nO3df5Bd5X3f8fdXuwJJUIRYKQxeQURmNU7xjzZmwaT2dFxLwGLZCA92g5uJth5m5ImxUJJOWjtlasZWMo7bqYsUm5EGFAvXY8qQdJDHgJCwkzqpTViwK4zllFtYoZWIkVeCYIkf+vHtH/csrFarZe8R957dve/XzM7e5znPufu9ml197nPOuc+JzESSpDJmVV2AJGn6MkQkSaUZIpKk0gwRSVJphogkqbTOqgtotYULF+aSJUuqLkOSpo3HHnvsF5m5aLxtbRciS5YsYWBgoOoyJGnaiIjdp9rm4SxJUmmGiCSpNENEklSaISJJKs0QUSm1Wo0VK1ZQq9WqLkVShQwRlbJu3ToOHTrEunXrqi5FUoUMETWsVqsxODgIwODgoLMRqY0ZImrY2NmHsxGpfRkiatjILORUbUntwxBRwyJiwrak9mGIqGFj74bp3TGl9mWISJJKM0TUsEsvvXTCtqT2YYioYQsXLpywLal9GCJq2Pe///0J25LahyGihi1fvpzOzvqtaDo7O7nyyisrrkhSVQwRNay/v59Zs+q/Oh0dHaxatariiiRVxRBRw7q6uujr6yMi6Ovro6urq+qSJFWk7W6Pq7dGf38/g4ODzkKkNmeIqJSuri7Wr19fdRmSKubhLElSaYaIJKk0Q0SSVJohIkkqzRCRJJVmiEiSSjNEJEmlNTVEIuL3I+LJiPhJRHwrIuZExMUR8UhEPBUR/yMizijGnlm0a8X2JaOe53NF/99HxNWj+vuKvlpEfLaZr0WSdLKmhUhEdAM3A72Z+U6gA7gB+FPgK5m5FDgI3FjsciNwMDN7gK8U44iIS4r93gH0AV+LiI6I6AC+ClwDXAJ8ohgrSWqRZh/O6gTmRkQnMA94DvggcG+xfQtwXfF4ZdGm2L4s6jfvXgncnZmvZuYzQA24vPiqZebTmfkacHcxVpLUIk0LkczcC/wX4Fnq4fEi8BjwQmYeLYYNAd3F425gT7Hv0WJ81+j+Mfucqv8kEbE6IgYiYmD//v2n/+IkSUBzD2ctoD4zuBh4G3AW9UNPY+XILqfY1mj/yZ2ZmzKzNzN7Fy1a9GalS5ImqZmHs5YDz2Tm/sw8Avwl8C+Ac4vDWwCLgX3F4yHgQoBi+3zgwOj+Mfucql+S1CLNDJFngSsiYl5xbmMZ8FPge8DHijH9wH3F461Fm2L7dzMzi/4biqu3LgaWAn8HPAosLa72OoP6yfetTXw9kqQxmrYUfGY+EhH3Ao8DR4EfAZuA7wB3R8S6ou/OYpc7gW9ERI36DOSG4nmejIh7qAfQUeCmzDwGEBGfAbZRv/Jrc2Y+2azXI0k6WdTf7LeP3t7eHBgYqLoMSZo2IuKxzOwdb5ufWJcklWaISJJKM0QkSaUZIpKk0gwRSVJphogkqTRDRJJUmiEiSSrNEJEklWaISJJKM0RUyvDwMDfffDPDw8NVlyKpQoaIStm4cSM7d+5k06ZNVZciqUKGiBo2PDzM9u3bAXjooYecjUhtzBBRwzZu3MjI6s+Z6WxEamOGiBq2Y8eOE9ojsxJJ7ccQUcPG3oOm3e5JI+kNhogaVr/b8anbktqHIaKGzZkzZ8K2pPZhiKhhhw8fnrAtqX0YIpKk0gwRNeyCCy6YsC2pfRgiatgLL7wwYVtS+zBE1LArr7zyhPZVV11VUSWSqmaIqGH9/f0ntFetWlVRJdL4XCC0dQwRNezgwYMTtqWqbdmyhSeeeIK77rqr6lJmPENEDVu3bt2EbalKw8PDPPjgg2QmDz74oLORJjNE1LDBwcEJ21KVtmzZwvHjxwE4duyYs5EmM0TUMJc90VS2Y8cOjh49CsDRo0ddILTJDBE1zAUYNZUtX76czs5OADo7O0+6mlBvLUNEDRv5Az1VW6pSf38/s2bV/2vr6Ojw6sEmM0TUsI6OjgnbUpW6urro6+sjIujr66Orq6vqkmY0Q0QNu/rqq09o9/X1VVSJNL7+/n7e9a53OQtpAUNEDevv73999tHZ2ekfqqacrq4u1q9f7yykBQwRNayrq4sVK1YQEaxYscI/VKmNGSIq5dprr2XevHl85CMfqboUSRUyRFTK1q1bOXz4MN/+9rerLkVShZoaIhFxbkTcGxE/i4hdEfGbEXFeRGyPiKeK7wuKsRER6yOiFhE7I+I9o56nvxj/VET0j+q/NCKeKPZZH37qrSVGLyvxwAMPuKyE1MaaPRO5DXgwM38d+GfALuCzwMOZuRR4uGgDXAMsLb5WA7cDRMR5wOeB9wKXA58fCZ5izOpR+3mZUAts2bKFI0eOAHDkyBGXlZDaWNNCJCLOAf4lcCdAZr6WmS8AK4EtxbAtwHXF45XAXVn3Q+DciLgAuBrYnpkHMvMgsB3oK7adk5k/yPpHpu8a9Vxqou3bt7/+KfXM5KGHHqq4IklVaeZM5NeA/cCfR8SPIuKOiDgLOD8znwMovv9KMb4b2DNq/6Gib6L+oXH6TxIRqyNiICIG9u/ff/qvrM2NvRpr4cKFFVUiqWrNDJFO4D3A7Zn5G8Ah3jh0NZ7xzmdkif6TOzM3ZWZvZvYuWrRo4qr1pvbt23dCe+/evRVVIqlqzQyRIWAoMx8p2vdSD5WfF4eiKL4/P2r8haP2Xwzse5P+xeP0q8lG1iU6VVtS+2jaX39m/gOwJyLeXnQtA34KbAVGrrDqB+4rHm8FVhVXaV0BvFgc7toGXBURC4oT6lcB24ptL0XEFcVVWatGPZea6LLLLpuwLal9NHv51TXANyPiDOBp4JPUg+ueiLgReBb4eDH2fuBDQA04XIwlMw9ExBeBR4txX8jMA8Xj3wW+DswFHii+1GR79uw5oT00NHSKkZJmuqaGSGb+GOgdZ9OyccYmcNMpnmczsHmc/gHgnadZpho0NjTGhoqk9uHBbDVsyZIlE7YltQ9DRA275ZZbJmxLah+GiBrW09Pz+uxjyZIl9PT0VFuQpMoYIirllltu4ayzznIWIrU5Q0SSVJoholLWrVvHoUOHWLduXdWlSCep1WqsWLGCWq1WdSkzniGihtVqNQYHBwEYHBz0D1VTjm9yWscQUcPG/mH6h6qpxDc5rWWIqGEjf6CnaktV8k1OaxkikmYU3+S0liEiaUZxRYXWMkQkzSiuqNBahoikGaWnp4cLLrgAgLe97W2uqNBkhoikGad+iyG1giEiaUap1Wqv38J53759XuLbZIaIpBnFS3xbyxCRNKN4iW9rGSKSZhQv8W0tQ0QNe/e73z1hW6qSl/i21qTusR4RZwLXA0tG75OZX2hOWZrKrrvuOnbu3Pl6+6Mf/WiF1UgnGrlp2uDgoDdNa4HJzkTuA1YCR4FDo77Uhr785S9P2Jaq5k3TWmdSMxFgcWb2NbUSTRuvvPLKCe2XX365okqk8fX09PCd73yn6jLawmRnIv87It7V1EokSdPOZGci7wf+bUQ8A7wKBJCZ6RlVSWpjkw2Ra5pahSRpWprU4azM3J2Zu4GXgRz1JUlTzvDwMDfffDPDw8NVlzLjTSpEIuLaiHgKeAb4a2AQeKCJdWkK6+jomLAtVW3jxo3s3LmTTZs2VV3KjDfZE+tfBK4A/m9mXgwsA/62aVVpSlu2bNkJ7eXLl1dUiXSy4eFhduzYAcD27dudjTTZZEPkSGYOA7MiYlZmfg/4502sS1PYpz71qdeX2o4IVq9eXXFF0hs2btzI8ePHATh+/LizkSabbIi8EBFnA98HvhkRt1H/4KHaUFdXF3PmzAFgzpw5dHV1VVyR9IaHH374hPbIrETNMdkQWUn9pPrvAQ8C/w/4SLOK0tRWq9Ve/4Dhyy+/7P0aNKWMvSGVN6hqrslenXUIWAj0AcPA3cXhLbWhW2+9dcK2VKWx5+zGtvXWmuwCjP8a+M/AX1H/oOGGiPjDzLy3ibVpihoaGpqwrfa1YcOGymemR44cOaG9Z88e1q5dW0ktPT09rFmzppKf3SqT/bDhfwQuy8znASJiEbADMEQkTSmzZ8+mo6ODY8eOsWDBAmbPnl11STPaZENk1kiAFIbxXiSSxpgq77o//elPs3v3bu644w4v/GiyyYbIgxGxDfhW0f4t4P7mlCRJp2f27Nn09PQYIC0wqRDJzD+MiOuB91E/J7IpM//nZPaNiA5gANibmR+OiIuBu4HzgMeB38nM14obX90FXEp9pvNbmTlYPMfngBuBY8DNmbmt6O8DbgM6gDsy80uTe9nT21Q47jxWVcecoT2OO0tT1aQPSWXmX2TmH2Tm7082QAprgV2j2n8KfCUzlwIHqYcDxfeDmdkDfKUYR0RcAtwAvIP61WFfi4iOIpy+Sn1xyEuATxRj1WTz58+fsC2pfUw4E4mIlxh/ocWRpeDPeZP9FwMrgD8G/iDqF2x/EPg3xZAtwK3A7dQ/i3Jr0X8v8GfF+JXULyl+FXgmImrA5cW4WmY+Xfysu4uxP52oppmg6nfdw8PDXH/99a+3N2/e7GEDqU1NGCKZ+U9O8/n/G/DvgZHn6QJeyMyRT7sPAd3F425gT/Fzj0bEi8X4buCHo55z9D57xvS/d7wiImI1sBrgoosuOo2XI6h/Yn3+/Pm8+OKLfOADHzBApDbWtCusIuLDwPOZ+djo7nGG5ptsa7T/5M7MTZnZm5m9ixYtmqBqTVZ3dzdnnXVW5bMiSdWa7NVZZbwPuDYiPgTMAc6hPjM5NyI6i9nIYmBfMX4IuBAYiohOYD5wYFT/iNH7nKpfTebVL5KgiTORzPxcZi7OzCXUT4x/NzN/G/ge8LFiWD9wX/F4a9Gm2P7dzMyi/4aIOLO4smsp8HfAo8DSiLg4Is4ofsbWZr0eSdLJmjkTOZX/ANwdEeuAHwF3Fv13At8oTpwfoB4KZOaTEXEP9RPmR4GbMvMYQER8BthG/RLfzZn5ZEtfiSS1uZaESGb+FfV1tyiuprp8nDGvAB8/xf5/TP0Kr7H99+OHHiWpMi5dIkkqzRCRJJVmiEiSSjNEJEmlGSKSpNIMEUlSaYaIJKk0Q0SSVJohIkkqzRCRJJVmiEiSSjNEJEmlGSKSpNIMEUlSaYaIJKk0Q0SSVJohIkkqzRCRJJVmiEiSSjNEJEmlGSKSpNIMEUlSaYaIJKk0Q0SSVJohIkkqzRCRJJVmiEiSSjNEJEmlGSKSpNIMEUlSaZ1VFzCdbNiwgVqtVnUZU8LIv8PatWsrrmRq6OnpYc2aNVWXIbWcIdKAWq3Gj3+yi2Pzzqu6lMrNei0BeOzpn1dcSfU6Dh+ougTANzmj+SbnRM18k2OINOjYvPN4+dc/VHUZmkLm/uz+qksA6v9xPvXkj7jo7GNVl1K5M47Uj9S/unug4kqq9+wvO5r6/IaININcdPYx/ug9/1h1GZpC/uTxc5r6/J5YlySVZohIkkprWohExIUR8b2I2BURT0bE2qL/vIjYHhFPFd8XFP0REesjohYROyPiPaOeq78Y/1RE9I/qvzQinij2WR8R0azXI0k6WTNnIkeBf5eZ/xS4ArgpIi4BPgs8nJlLgYeLNsA1wNLiazVwO9RDB/g88F7gcuDzI8FTjFk9ar++Jr4eSdIYTQuRzHwuMx8vHr8E7AK6gZXAlmLYFuC64vFK4K6s+yFwbkRcAFwNbM/MA5l5ENgO9BXbzsnMH2RmAneNei5JUgu05JxIRCwBfgN4BDg/M5+DetAAv1IM6wb2jNptqOibqH9onP7xfv7qiBiIiIH9+/ef7suRJBWaHiIRcTbwF8DvZeZE1x6Odz4jS/Sf3Jm5KTN7M7N30aJFb1ayJGmSmhoiETGbeoB8MzP/suj+eXEoiuL780X/EHDhqN0XA/vepH/xOP2SpBZp5tVZAdwJ7MrM/zpq01Zg5AqrfuC+Uf2riqu0rgBeLA53bQOuiogFxQn1q4BtxbaXIuKK4metGvVckqQWaOYn1t8H/A7wRET8uOj7I+BLwD0RcSPwLPDxYtv9wIeAGnAY+CRAZh6IiC8CjxbjvpCZI4sV/S7wdWAu8EDxJUlqkaaFSGb+DeOftwBYNs74BG46xXNtBjaP0z8AvPM0ypQknQbXzmrA3r176Tj84pRZcE9TQ8fhYfbuPVp1GVIlDBFphti7dy+HXupo+oJ7ml52v9TBWXv3Nu35DZEGdHd38w+vdroUvE4w92f30919ftVlSJUwRKQZoru7m1ePPudS8DrBnzx+Dmd2j/s57LeEq/hKkkozRCRJpRkikqTSDBFJUmmGiCSpNENEklSaISJJKs3PiTSo4/ABlz0BZr1S/yzC8Tl+Orrj8AHADxuqPRkiDejp6am6hCmjVnsJgJ5f8z9PON/fDbUtQ6QBa9asqbqEKWPt2rUA3HbbbRVXIqlKnhORJJVmiEiSSjNEJEmleU5EmkGe/aX3EwH4+eH6++Pz5x2vuJLqPfvLDpY28fkNEWmG8AqxN7xWqwFw5q/6b7KU5v5uGCLSDOHVg2/w6sHW8ZyIJKk0Q0SSVJohIkkqzRCRJJVmiEiSSjNEJEmlGSKSpNIMEUlSaYaIJKk0Q0SSVJohIkkqzRCRJJVmiEiSSnMV32low4YN1Iqlrqsy8vNHVkutUk9PjyvYShUxRFTK3Llzqy5BU9BUeIMDU+dNTju8wTFEpqGZ/kspnS7f5LTOtA+RiOgDbgM6gDsy80sVlyS1Ld/gtJ9pfWI9IjqArwLXAJcAn4iIS6qtSpLax7QOEeByoJaZT2fma8DdwMqKa5KktjHdQ6Qb2DOqPVT0nSAiVkfEQEQM7N+/v2XFSdJMN91DJMbpy5M6MjdlZm9m9i5atKgFZUlSe5juITIEXDiqvRjYV1EtktR2pnuIPAosjYiLI+IM4AZga8U1SVLbmNaX+Gbm0Yj4DLCN+iW+mzPzyYrLkqS2Ma1DBCAz7wfur7oOSWpHkXnSeegZLSL2A7urrmOGWAj8ouoipFPw9/Ot86uZOe5VSW0XInrrRMRAZvZWXYc0Hn8/W2O6n1iXJFXIEJEklWaI6HRsqroAaQL+fraA50QkSaU5E5EklWaISJJKM0RUSkT0RcTfR0QtIj5bdT3SiIjYHBHPR8RPqq6lHRgiapg3A9MU93Wgr+oi2oUhojK8GZimrMz8X8CBqutoF4aIypjUzcAkzXyGiMqY1M3AJM18hojK8GZgkgBDROV4MzBJgCGiEjLzKDByM7BdwD3eDExTRUR8C/gB8PaIGIqIG6uuaSZz2RNJUmnORCRJpRkikqTSDBFJUmmGiCSpNENEklSaISK1QEQsKbuq7OnsKzWbISJJKs0QkVqnMyK2RMTOiLg3IuZFxH+KiEcj4icRsSkiAiAiLo2I/xMRPwBuqrhu6ZQMEal13g5sysx3A/8IfBr4s8y8LDPfCcwFPlyM/XPg5sz8zWpKlSbHEJFaZ09m/m3x+L8D7wf+VUQ8EhFPAB8E3hER84FzM/Ovi7HfqKBWaVI6qy5AaiNj1xhK4GtAb2buiYhbgTnUl9p3PSJNC85EpNa5KCJGDk99Avib4vEvIuJs4GMAmfkC8GJEvL/Y/tutLVOaPGciUuvsAvojYiPwFHA7sAB4AhikvsT+iE8CmyPiMPXVkqUpyVV8JUmleThLklSaISJJKs0QkSSVZohIkkozRCRJpRkikqTSDBFJUmn/H6qZg18YoJVKAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEGCAYAAACpXNjrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAe2klEQVR4nO3dfZBX1Z3n8feHRo3EGLXtWKbRYELvJOgaoh3DVrSWKGDrRnEmplZ2KvZmrYXN+JSH2lXH7Bh1kjJblXHFJJYkMjaZJIQySdmmENI8uMYpH2gSAyK69PjYwGrboNFhhAG++8c9rb9ufjQN3vu7/fB5Vd363fu95957fkr3t889596jiMDMzCxP48qugJmZjT5OLmZmljsnFzMzy52Ti5mZ5c7JxczMcje+7AoMF8cff3xMmjSp7GqYmY0oa9eufS0iGgbGnVySSZMm0dnZWXY1zMxGFEkvVov7tpiZmeXOycXMzHLn5GJmZrlzcjEzs9w5uZjZmNHb28s111xDb29v2VUZ9QpPLpLqJP1B0m/S9imSHpe0SdIvJB2e4kek7a60f1LFOW5I8WclnV8Rb0mxLknXV8SrXsPMxra7776bdevWsWDBgrKrMurVouVyLbCxYvu7wO0R0QRsB65I8SuA7RExGbg9lUPSFOAy4FSgBfhhSlh1wA+AC4ApwJxUdrBrmNkY1dvby4oVKwDo6Ohw66VghSYXSROB/wD8OG0LOBe4LxVpAy5J67PTNmn/ean8bGBxROyMiOeBLuCstHRFxHMRsQtYDMw+wDXMbIy6++672bt3LwB79+5166VgRbdc/jfwP4C9abseeD0idqftbqAxrTcCLwOk/W+k8u/EBxyzv/hg1+hH0lxJnZI6e3p6DvU7mtkIsHLlyn7bfa0YK0ZhyUXS54FXI2JtZbhK0TjAvrzi+wYjFkREc0Q0NzTs8/YCMxtFspsa+9+2fBXZcvkscLGkF8huWZ1L1pI5RlLfa2cmAlvSejdwEkDa/0FgW2V8wDH7i782yDXMbIw677zzBt22fBWWXCLihoiYGBGTyDrkV0XEXwKrgUtTsVbg/rTenrZJ+1dFNgdzO3BZGk12CtAEPAGsAZrSyLDD0zXa0zH7u4aZjVFz585l3LjsV964ceOYO3duyTUa3cp4zuU64OuSusj6R+5J8XuA+hT/OnA9QERsAJYATwPLgCsjYk/qU7kKWE42Gm1JKjvYNcxsjKqvr2fmzJkAzJw5k/r6+pJrNLop+0Pfmpubw29FNhvdent7ufnmm7npppucXHIiaW1ENA+M+5X7ZjZm1NfXM3/+/LKrMSb49S9mZpY7JxczM8udk4uZmeXOycXMzHLn5GK582vNzczJxXLX1tbG+vXrWbRoUdlVMbOSOLlYrnp7e1m2bBkRwbJly9x6MRujnFwsV21tbezZsweA3bt3u/ViNkY5uViuVqxY8U5y2bNnDx0dHSXXyMzK4ORiuTr77LP7bZ9zzjkl1cTMyuTkYrnatWtXv+2dO3eWVBMzK5OTi+XqkUceGXTbzMYGJxfL1cC3bPut22Zjk5OL5Wrg7H4zZswoqSZmViYnF8vVvHnzPNufmRWXXCS9T9ITkv4oaYOkm1P8XknPS3oyLVNTXJLmS+qStE7SGRXnapW0KS2tFfEzJa1Px8yXpBQ/TlJHKt8h6diivqf1V19f/05rxbP9mY1dRbZcdgLnRsQngalAi6Rpad9/j4ipaXkyxS4AmtIyF7gLskQB3AR8BjgLuKkiWdyVyvYd15Li1wMrI6IJWJm2rUbmzZvH6aef7laL2RhWWHKJzFtp87C0DNa7OxtYlI57DDhG0onA+UBHRGyLiO1AB1miOhE4OiIejazXeBFwScW52tJ6W0XcaqBvtj+3WszGrkL7XCTVSXoSeJUsQTyedn073fq6XdIRKdYIvFxxeHeKDRbvrhIHOCEitgKkzw/tp35zJXVK6uzp6Tnk72lmZv0VmlwiYk9ETAUmAmdJOg24Afg48GngOOC6VFzVTnEI8YOp34KIaI6I5oaGhoM51MzMBlGT0WIR8TrwENASEVvTra+dwN+T9aNA1vI4qeKwicCWA8QnVokDvJJum5E+X831C5mZ2aCKHC3WIOmYtH4kMAN4puKXvsj6Qp5Kh7QDl6dRY9OAN9ItreXALEnHpo78WcDytO9NSdPSuS4H7q84V9+ostaKuJmZ1cD4As99ItAmqY4siS2JiN9IWiWpgey21pPAf0vllwIXAl3ADuDLABGxTdKtwJpU7paI2JbWvwLcCxwJPJgWgNuAJZKuAF4CvljYtzQzs33Ir+fINDc3R2dnZ9nVMDMbUSStjYjmgXE/oW9mZrlzcjEzs9w5uZiZWe6cXMzMLHdOLmZmljsnFzMzy52Ti5mZ5c7JxczMcufkYmZmuXNyMTOz3Dm5mJlZ7pxczMwsd04uZmaWOycXMzPLnZOLmZnlzsnFzMxyV+Q0x++T9ISkP0raIOnmFD9F0uOSNkn6haTDU/yItN2V9k+qONcNKf6spPMr4i0p1iXp+op41WuYmVltFNly2QmcGxGfBKYCLZKmAd8Fbo+IJmA7cEUqfwWwPSImA7enckiaAlwGnAq0AD+UVJemT/4BcAEwBZiTyjLINczMrAYKSy6ReSttHpaWAM4F7kvxNuCStD47bZP2nydJKb44InZGxPNAF3BWWroi4rmI2AUsBmanY/Z3DTMzq4FC+1xSC+NJ4FWgA/gn4PWI2J2KdAONab0ReBkg7X8DqK+MDzhmf/H6Qa4xsH5zJXVK6uzp6XkvX9XMzCoUmlwiYk9ETAUmkrU0PlGtWPrUfvblFa9WvwUR0RwRzQ0NDdWKmJnZIajJaLGIeB14CJgGHCNpfNo1EdiS1ruBkwDS/g8C2yrjA47ZX/y1Qa5hZmY1UORosQZJx6T1I4EZwEZgNXBpKtYK3J/W29M2af+qiIgUvyyNJjsFaAKeANYATWlk2OFknf7t6Zj9XcPMzGpg/IGLHLITgbY0qmscsCQifiPpaWCxpL8F/gDck8rfA/xEUhdZi+UygIjYIGkJ8DSwG7gyIvYASLoKWA7UAQsjYkM613X7uYaZmdWAsj/0rbm5OTo7O8uuhpnZiCJpbUQ0D4z7CX0zM8udk4uZmeXOycXMzHLn5GJmZrlzcjEzs9w5uZiZWe6cXMzMLHdOLmZmljsnFzMzy52Ti5mZ5c7JxczMcufkYmZmuXNyMTOz3Dm5mJlZ7pxczMwsd04uZmaWuyKnOT5J0mpJGyVtkHRtin9L0mZJT6blwopjbpDUJelZSedXxFtSrEvS9RXxUyQ9LmmTpF+k6Y5JUyL/IpV/XNKkor6nmZntq8iWy27gGxHxCWAacKWkKWnf7RExNS1LAdK+y4BTgRbgh5Lq0jTJPwAuAKYAcyrO8910riZgO3BFil8BbI+IycDtqZyZmdVIYcklIrZGxO/T+pvARqBxkENmA4sjYmdEPA90AWelpSsinouIXcBiYLYkAecC96Xj24BLKs7VltbvA85L5c3MrAZq0ueSbkt9Cng8ha6StE7SQknHplgj8HLFYd0ptr94PfB6ROweEO93rrT/jVR+YL3mSuqU1NnT0/OevqO9a9WqVUyfPp3Vq1eXXRUzK0nhyUXSUcAvga9GxJ+Au4CPAVOBrcD3+opWOTwOIT7YufoHIhZERHNENDc0NAz6PWzobr31VgBuueWWkmtiZmUpNLlIOowssfw0In4FEBGvRMSeiNgL/IjsthdkLY+TKg6fCGwZJP4acIyk8QPi/c6V9n8Q2Jbvt7NqVq1aRUSWxyPCrRezMarI0WIC7gE2RsTfVcRPrCj258BTab0duCyN9DoFaAKeANYATWlk2OFknf7tkf0GWw1cmo5vBe6vOFdrWr8UWBV9v/GsUH2tlj5uvdhw0tvbyzXXXENvb2/ZVRn1imy5fBb4EnDugGHH/0vSeknrgM8BXwOIiA3AEuBpYBlwZWrh7AauApaTDQpYksoCXAd8XVIXWZ/KPSl+D1Cf4l8H3hm+bMUamMOd0204aWtrY/369SxatKjsqox68g9/prm5OTo7O8uuxog3ffr0fWIPPfRQzethNlBvby9z5sxh165dHHHEEfzsZz+jvn6fcT52kCStjYjmgfEhtVwk/RtJKyU9lbZPl/TNvCtpZlaUtrY29u7dC8CePXvceinYUG+L/Qi4AfhXgIhYR9b3YWY2IqxYsYLdu7MnF3bv3k1HR0fJNRrdhppcJkTEEwNiu6uWNDMbhmbMmEHfs9SSmDlzZsk1Gt2Gmlxek/Qx0rMiki4le0bFzGxEuPjii/sNk7/oootKrtHoNtTkciVwN/BxSZuBrwJfKaxWZmY5a29v79dyeeCBB0qu0eg2pOSS3us1A2gAPh4RZ0fEC4XWzMwsRytWrOjXcnGfS7HGH7gISPqbAdsARISfkLN+Zs2axW9/+9t3ts8///xBSpvVzowZM3jggQeICPe51MBQb4v9c8Wyh+z195MKqpONYPPmzeu3PXfu3JJqYtaf+1xqa6i3xb5XsXwbmM7gr8+3Maq+vp5Zs2YBWavFD6nZcOE+l9o61Ne/TAA+mmdFbPSYN28ep59+ulstNqy4z6W2hvqE/vo0/8o6SRuAZ4E7iq2ajVT19fXMnz/frRYbVmbMmMH48Vk38/jx493nUrChtlw+D1yUllnAhyPi+4XVykY0v3nWhqPW1lbGjct+5dXV1XH55ZeXXKPRbdDkIuk4SccBb1Ys/wIcneJm+7jxxhtZt24d3/ymXz9nw0d9fT0tLS1IoqWlxS3rgh2o5bIW6EyfPcD/BTal9bXFVs1Got7eXp555hkANm7c6NaLDSsXX3wxEyZM8EixGhg0uUTEKRHxUbK5VC6KiOMjop7sNtmvalFBG1luvPHGfttuvdhw0t7ezo4dOzxSrAaG2ufy6YhY2rcREQ8C/76YKtlI1tdq6bNx48aSamLWX29vL8uWLSMiWLZsmVvVBTuYF1d+U9IkSR+RdCMw6P8ZSSdJWi1po6QNkq5N8eMkdUjalD6PTXFJmi+pK41KO6PiXK2p/CZJrRXxM9NItq50rAa7hpmNXZ7PpbaGmlzmkL1X7NdpaeDA87nsBr4REZ8ApgFXSppCNuXwyohoAlby7hTEFwBNaZkL3AVZogBuAj4DnAXcVJEs7kpl+45rSfH9XcPMxijP51JbQ00u50XEtRHxqYg4IyK+CswY7ICI2BoRv0/rbwIbyZ7qnw20pWJtwCVpfTawKDKPAcdIOhE4H+iIiG0RsR3oAFrSvqMj4tHInoxaNOBc1a5hZmPUOeecM+i25WuoyeWGIcaqkjQJ+BTwOHBCRGyFLAEBH0rFGoGXKw7rTrHB4t1V4gxyjYH1miupU1JnT0/PUL+OmY1AfU/nW20c6DmXCyTdCTSmPo2+5V6GOBOlpKOAXwJfjYg/DVa0SiwOIT5kEbEgIpojormhoeFgDjWzEeaRRx7pt/273/2upJqMDQdquWwhe87lbbLnWvqWdrLbVYOSdBhZYvlpRPQNXX4l3dIifb6a4t3ASRWHT0zXHyw+sUp8sGuY2Rg1Y8YM6urqgOwJfb/+pVgHes7lj8A/AI9ERFvF8qvU/7FfaeTWPcDGiPi7il3tQN+Ir1bg/or45WnU2DTgjXRLazkwS9KxqSN/FrA87XtT0rR0rcsHnKvaNcxsjGptbX0nuYwfP96vfynYAftcImIPUC/p8IM892eBLwHnSnoyLRcCtwEzJW0CZqZtgKXAc0AX8CPgr9L1twG3AmvSckuKQTbV8o/TMf8EPJji+7uGmY1R9fX1TJ8+HYDp06f79S8FG9JMlMCLwD9KaiebMAyAAS2SfiLiEar3iwCcV6V8AFfu51wLgYVV4p3AaVXivdWuYWZjW998Lla8oY4W2wL8JpX/QMViZjYi9Pb2snr1agAeeughP6FfsCG1XCLiZgBJH8g2461Ca2VmlrNqT+h/7WtfK7lWo9dQJws7TdIfgKeADZLWSjq12KqZmeXHT+jX1lBviy0Avh4RH4mIjwDfIOt0NzMbETwTZW0NNbm8PyJW921ExEPA+wupkZlZATwTZW0NNbk8J+l/prciT5L0TeD5IitmZpYnz0RZW0NNLv+F7E3IvySbJOx44D8XVCcbwfr+MtzftlmZPBNl7Qz1J/9jZK9gGQccRvYMycNFVcpGrr7ROPvbNiuTZ6KsnaEml5+SPcT4F2RTHH8ecOo3sxHDM1HW1lCTS09EPBARz0fEi31LoTUzM8uRZ6KsraEml5sk/VjSHEl/0bcUWjMbkSZMmDDotllZ/JxLbQ01uXwZmEo2jfBFafl8UZWykWvHjh2DbpuVxc+51NZQX1z5yYj4t4XWxEaFSZMm8cILL/TbNhsOWltbefDB7MXp48aN83MuBRtqy+UxSVMKrYmNCldddVW/7auvvrqkmpj1V19fT2NjNhP6hz/8YT/nUrChJpezgSclPStpnaT1ktYVWTEbmR5++OFBt83K0tvby5Yt2WS1W7Zs8Wixgg01ubQATWSzQPb1t3gosu1jxYoV/bbdaWrDReVosb1793q0WMGUzdFVwImlhWRJ6NWIOC3FvgX8V6AnFfvriFia9t0AXAHsAa6JiOUp3gLcAdQBP46I21L8FGAxcBzwe+BLEbFL0hHAIuBMoBf4jxHxwoHq29zcHJ2dnTl88/LceeeddHV1lVqH7u7ufn8R1tfXM3HixFLqMnnyZN+Ws3dceOGF/QaYTJgwgaVLl5ZYo9FB0tqIaB4YL/LdHPeStXgGuj0ipqalL7FMAS4DTk3H/FBSnaQ64AfABcAUYE5F389307magO1kiYn0uT0iJgO3p3JWIyeccMI765L6bZuVyaPFamuoo8UOWkQ8LGnSEIvPBhZHxE7geUldwFlpX1dEPAcgaTEwW9JG4FzgP6UybcC3gLvSub6V4vcB35ekKKqJNowMl7/Sv/CFL9Db28vFF1/syZhs2GhtbWXZsmWA34pcC4Ull0FcJelyoBP4RkRsBxqBxyrKdKcYwMsD4p8B6oHXI2J3lfKNfcdExG5Jb6Tyrw2siKS5wFyAk08++b1/MwOy1svbb7/tH17rZzjctpUEwFFHHcUtt9xSal1G+23bWr+y9i6yl2BOBbYC30txVSkbhxAf7Fz7BiMWRERzRDQ3NDQMVm87CIcddhiTJ0/2UE8bdsaNG8e4ceN8u7YGatpyiYhX+tYl/Qj4TdrsJnvrcp+JwJa0Xi3+GnCMpPGp9VJZvu9c3ZLGAx8EtuX8VczsIA2Hv9KvvfZaAO64446SazL61bTlIunEis0/B55K6+3AZZKOSKPAmoAngDVAk6RTJB1O1unfnvpPVgOXpuNbgfsrztWa1i8FVo2F/hYzs+GksJaLpJ8D04HjJXUDNwHTJU0lu031AjAPICI2SFoCPA3sBq6MiD3pPFcBy8mGIi+MiA3pEtcBiyX9LfAH4J4Uvwf4SRoUsI0sIZmZWQ0VOVpsTpXwPVVifeW/DXy7SnwpsM9g9DSC7Kwq8beBLx5UZc3MLFeeg9bMzHLn5GJmZrlzcjEzs9w5uZiZWe6cXMzMLHdOLmZmljsnFzMzy52Ti5mZ5c7JxczMcufkYmZmuXNyMTOz3Dm5mJlZ7pxczMwsd04uZmaWOycXMzPLnZOLmZnlrrDkImmhpFclPVURO05Sh6RN6fPYFJek+ZK6JK2TdEbFMa2p/CZJrRXxMyWtT8fMl6TBrmFmZrVTZMvlXqBlQOx6YGVENAEr0zbABUBTWuYCd0GWKMimR/4M2ayTN1Uki7tS2b7jWg5wDTMzq5HCkktEPEw2h32l2UBbWm8DLqmIL4rMY8Axkk4Ezgc6ImJbRGwHOoCWtO/oiHg0IgJYNOBc1a5hZmY1Uus+lxMiYitA+vxQijcCL1eU606xweLdVeKDXWMfkuZK6pTU2dPTc8hfyszM+hsuHfqqEotDiB+UiFgQEc0R0dzQ0HCwh5uZ2X7UOrm8km5pkT5fTfFu4KSKchOBLQeIT6wSH+waZmZWI+NrfL12oBW4LX3eXxG/StJiss77NyJiq6TlwHcqOvFnATdExDZJb0qaBjwOXA7ceYBrFOrOO++kq6urFpca9vr+O1x77bUl12R4mDx5MldffXXZ1TCrqcKSi6SfA9OB4yV1k436ug1YIukK4CXgi6n4UuBCoAvYAXwZICWRW4E1qdwtEdE3SOArZCPSjgQeTAuDXKNQXV1dPPnURvZMOK4WlxvWxu3K7lCufe6VkmtSvrodA8e0mI0NhSWXiJizn13nVSkbwJX7Oc9CYGGVeCdwWpV4b7Vr1MKeCcfxLx+/sIxL2zB15DNLy66CWSmGS4e+mZmNIrXuczGzGnN/4LvcH9hfkf2BTi5mo1xXVxebNvyBk4/aU3ZVSnf4v2Y3a3a+2FlyTcr30lt1hZ7fycVsDDj5qD389Rl/KrsaNox85/dHF3p+97mYmVnunFzMzCx3Ti5mZpY7JxczM8udk4uZmeXOo8VysnnzZup2vOEnsq2fuh29bN68u+xqmNWcWy5mZpY7t1xy0tjYyP/bOd7vFrN+jnxmKY2NJ5RdDbOac3IxG+U2b97MP79ZV/hDczayvPhmHe/fvLmw8/u2mJmZ5c4tF7NRrrGxkZ27t/r1L9bPd35/NEc0NhZ2frdczMwsd6UkF0kvSFov6UlJnSl2nKQOSZvS57EpLknzJXVJWifpjIrztKbymyS1VsTPTOfvSseq9t/SzGzsKrPl8rmImBoRzWn7emBlRDQBK9M2wAVAU1rmAndBlozIpk7+DHAWcFNfQkpl5lYc11L81zEzsz7Dqc9lNjA9rbcBDwHXpfiiNBXyY5KOkXRiKtsREdsAJHUALZIeAo6OiEdTfBFwCfBg0V+gbsc2P0QJjHs7u7e/930enVS3Yxvgocg29pSVXAL4raQA7o6IBcAJEbEVICK2SvpQKtsIvFxxbHeKDRbvrhLfh6S5ZC0cTj755Pf0hSZPnvyejh9NurreBGDyR/1LFU7wvw0bk8pKLp+NiC0pgXRIemaQstX6S+IQ4vsGs6S2AKC5ublqmaEqaqrQkahvCtk77rij5JqYWVlK6XOJiC3p81Xg12R9Jq+k212kz1dT8W7gpIrDJwJbDhCfWCVuZmY1UvPkIun9kj7Qtw7MAp4C2oG+EV+twP1pvR24PI0amwa8kW6fLQdmSTo2deTPApanfW9KmpZGiV1ecS4zM6uBMm6LnQD8Oo0OHg/8LCKWSVoDLJF0BfAS8MVUfilwIdAF7AC+DBAR2yTdCqxJ5W7p69wHvgLcCxxJ1pFfeGe+mZm9q+bJJSKeAz5ZJd4LnFclHsCV+znXQmBhlXgncNp7rqzZKPHSW363GMArO7KbNSdM2FtyTcr30lt1NBV4/uE0FNnMCuDRau/a1dUFwBEf8X+TJor9t+HkYjbKeSTjuzySsXb8bjEzM8udk4uZmeXOycXMzHLn5GJmZrlzcjEzs9w5uZiZWe6cXMzMLHdOLmZmljsnFzMzy52Ti5mZ5c7JxczMcufkYmZmuXNyMTOz3PmtyGZWE3feeSdd6ZX3Zem7ft/bkcs0efLkUf3G6lGbXCS1AHcAdcCPI+K2kqtUuOHwwwvD5wd4tP/w2sE78sgjy67CmDEqk4ukOuAHwEygG1gjqT0ini63ZmODf4CtGif6sWVUJhfgLKArTamMpMXAbGBUJxf/8JrZcDFaO/QbgZcrtrtTrB9JcyV1Surs6empWeXMzEa70ZpcVCUW+wQiFkREc0Q0NzQ01KBaZmZjw2hNLt3ASRXbE4EtJdXFzGzMGa3JZQ3QJOkUSYcDlwHtJdfJzGzMGJUd+hGxW9JVwHKyocgLI2JDydUyMxszRmVyAYiIpcDSsuthZjYWjdbbYmZmViInFzMzy50i9hmhOyZJ6gFeLLseo8jxwGtlV8KsCv/bzNdHImKfZzmcXKwQkjojornsepgN5H+bteHbYmZmljsnFzMzy52TixVlQdkVMNsP/9usAfe5mJlZ7txyMTOz3Dm5mJlZ7pxcLFeSWiQ9K6lL0vVl18esj6SFkl6V9FTZdRkLnFwsNxXTS18ATAHmSJpSbq3M3nEv0FJ2JcYKJxfL0zvTS0fELqBvemmz0kXEw8C2susxVji5WJ6GNL20mY1+Ti6WpyFNL21mo5+Ti+XJ00ubGeDkYvny9NJmBji5WI4iYjfQN730RmCJp5e24ULSz4FHgT+T1C3pirLrNJr59S9mZpY7t1zMzCx3Ti5mZpY7JxczM8udk4uZmeXOycXMzHLn5GJWMkmTDvVNve/lWLMiObmYmVnunFzMhofxktokrZN0n6QJkv5G0hpJT0laIEkAks6U9EdJjwJXllxvs6qcXMyGhz8DFkTE6cCfgL8Cvh8Rn46I04Ajgc+nsn8PXBMR/66cqpodmJOL2fDwckT8Y1r/B+Bs4HOSHpe0HjgXOFXSB4FjIuL/pLI/KaGuZgc0vuwKmBmw79QEAfwQaI6IlyV9C3gf2bQGfmeTDXtuuZgNDydL6rvNNQd4JK2/Juko4FKAiHgdeEPS2Wn/X9a2mmZD45aL2fCwEWiVdDewCbgLOBZYD7xANp1Bny8DCyXtIHsDtdmw47cim5lZ7nxbzMzMcufkYmZmuXNyMTOz3Dm5mJlZ7pxczMwsd04uZmaWOycXMzPL3f8HAfC0N+0N4gYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEGCAYAAACpXNjrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAXIklEQVR4nO3df5BdZZ3n8fc3HUF+iEATUk4HDTudchct15EeZVe3dkYCBHcEq/yxulNLr0VNKBcDU27VGrFqYbSkmCpqXZIVarsWlmZndpByd8o4hGCLij8qKo3jJCI66UGEBMaEDiBOGCTJd/+4T+Cmuek04bn3dPq+X1Vdfc/3Pueeb6guPve555znRmYiSVJNi5puQJK08BgukqTqDBdJUnWGiySpOsNFklTd4qYbmC9OO+20XL58edNtSNJR5f77738iM5fMrBsuxfLly5mcnGy6DUk6qkTELzrV/VhMklSd4SJJqs5wkSRVZ7hIkqozXCT1jenpaa644gqmp6ebbmXBM1wk9Y2xsTG2bNnC2NhY060seIaLpL4wPT3NxMQEABMTE85eusxwkdQXxsbG2L9/PwD79+939tJlhoukvnDPPffMuq26DBdJfWHmFyP6RYnd1dVwiYiHI2JrRPwoIiZL7dSImIiIbeX3KaUeEbEuIqYiYktEvK3tdUbL+G0RMdpWP7u8/lTZN2Y7hqT+de655x60vXLlyoY66Q+9mLn8fma+NTNHyvZa4J7MXAHcU7YBLgRWlJ/VwE3QCgrgauAdwNuBq9vC4qYy9sB+qw5zDEl96rLLLmPRotb/8hYtWsTq1asb7mhha+JjsYuB8fJ4HHhfW/22bPkecHJEvA64AJjIzN2Z+SQwAawqz52UmZuzNb+9bcZrdTqGpD41ODj4wmzlvPPOY3BwsOGOFrZuh0sCX42I+yPiwNuEpZn5OED5fXqpDwGPtu27vdRmq2/vUJ/tGAeJiNURMRkRk7t27TrCf6Kko8Vll13GW97yFmctPdDtJfffmZmPRcTpwERE/HSWsdGhlkdQn7PMHAPGAEZGRjy7Jy1wg4ODrFu3ruk2+kJXZy6Z+Vj5vRP4S1rnTH5ZPtKi/N5Zhm8HzmjbfRnw2GHqyzrUmeUYkqQe6Fq4RMQJEfGaA4+B84EfAxuAA1d8jQJfLo83AJeUq8bOAZ4uH2ndDZwfEaeUE/nnA3eX556JiHPKVWKXzHitTseQJPVANz8WWwr8Zbk6eDHwfzJzU0TcB9wREZcCjwAfLOM3Au8BpoA9wEcBMnN3RHwWuK+M+0xm7i6PPwbcChwH3FV+AK47xDEkST0Q3kjUMjIykn7NsSS9PBFxf9utJi/wDn1JUnWGiySpOsNFklSd4aLq/LY/SYaLqlu/fj1btmxh/fr1TbciqSGGi6qanp7m3nvvBeDee+919iL1KcNFVa1fv/6F78nITGcvUp8yXFTVgVnLobYl9QfDRVX5bX+SwHBRZcuWLZt1W1J/MFxU1TXXXDPrtqT+YLioquHh4RdmK8uWLWN4eLjhjiQ1wXBRdddccw0nnHCCsxapj3X7myjVh4aHh7nzzjubbkNSg5y5SJKqM1wkSdUZLpKk6gwXSVJ1hoskqTrDRZJUneEiSarOcJEkVWe4SJKqM1wkSdUZLpKk6gwXSVJ1hoskqTrDRZJUneEiSarOcJEkVWe4SJKqM1wkSdV1PVwiYiAi/joi/qpsnxkR34+IbRHxxYg4ptSPLdtT5fnlba/xqVL/WURc0FZfVWpTEbG2rd7xGJKk3ujFzOVK4MG27T8FPp+ZK4AngUtL/VLgycwcBj5fxhERZwEfBt4ErAJuLIE1AHwBuBA4C/hIGTvbMSRJPdDVcImIZcC/Af5n2Q7g3cCXypBx4H3l8cVlm/L8uWX8xcDtmflcZv4cmALeXn6mMvOhzPwNcDtw8WGOIUnqgW7PXP4b8J+B/WV7EHgqM/eW7e3AUHk8BDwKUJ5/uox/oT5jn0PVZzvGQSJidURMRsTkrl27jvTfKEmaoWvhEhF/AOzMzPvbyx2G5mGeq1V/aTFzLDNHMnNkyZIlnYZIko7A4i6+9juBiyLiPcCrgZNozWROjojFZWaxDHisjN8OnAFsj4jFwGuB3W31A9r36VR/YpZjSJJ6oGszl8z8VGYuy8zltE7Ifz0z/xD4BvCBMmwU+HJ5vKFsU57/emZmqX+4XE12JrAC+AFwH7CiXBl2TDnGhrLPoY4hSeqBJu5z+STwiYiYonV+5OZSvxkYLPVPAGsBMvMB4A7gJ8Am4PLM3FdmJR8H7qZ1NdodZexsx5Ak9UC03uhrZGQkJycnm25Dko4qEXF/Zo7MrHuHviSpOsNFklSd4SJJqs5wkSRVZ7hIkqozXCRJ1RkukqTqDBdJUnWGiySpOsNFklSd4SJJqs5wkSRVZ7hIkqozXCRJ1RkukqTqDBdJUnWGiySpOsNFklSd4SJJqs5wkSRVZ7hIkqozXCRJ1RkukqTqDBdJUnWGiySpOsNFklSd4SJJqs5wkSRVZ7hI6hvT09NcccUVTE9PN93Kgme4SOob4+PjbN26ldtuu63pVhY8w0XV+e5Q89H09DSbNm0iM7nrrrv8++yyroVLRLw6In4QEX8TEQ9ExJ+U+pkR8f2I2BYRX4yIY0r92LI9VZ5f3vZanyr1n0XEBW31VaU2FRFr2+odj6HeuP7669myZQvXX399061ILxgfH+f5558H4Pnnn3f20mXdnLk8B7w7M/858FZgVUScA/wp8PnMXAE8CVxaxl8KPJmZw8Dnyzgi4izgw8CbgFXAjRExEBEDwBeAC4GzgI+UscxyDHXZ9PQ0mzdvBmDz5s2+O9S8MTExQWYCkJl89atfbbijha1r4ZItvy6bryo/Cbwb+FKpjwPvK48vLtuU58+NiCj12zPzucz8OTAFvL38TGXmQ5n5G+B24OKyz6GOoS6bOVtx9qL5YnBw8KDt0047raFO+kNXz7mUGcaPgJ3ABPB3wFOZubcM2Q4MlcdDwKMA5fmngcH2+ox9DlUfnOUYM/tbHRGTETG5a9euV/JPVXFg1nKobakpjz322EHbO3bsaKiT/tDVcMnMfZn5VmAZrZnGP+s0rPyOQzxXq96pv7HMHMnMkSVLlnQaImmBWLRo0azbqqsn/3Uz8yngm8A5wMkRsbg8tQw48HZiO3AGQHn+tcDu9vqMfQ5Vf2KWY0jqU+9617tm3VZdhw2XiFgaETdHxF1l+6yIOOwJ8ohYEhEnl8fHASuBB4FvAB8ow0aBL5fHG8o25fmvZ+vs2wbgw+VqsjOBFcAPgPuAFeXKsGNonfTfUPY51DEk9aljjjn4otFjjz22oU76w1xmLrcCdwO/Vbb/FvjjOez3OuAbEbGFVhBMZOZfAZ8EPhERU7TOj9xcxt8MDJb6J4C1AJn5AHAH8BNgE3B5+bhtL/Dx0tuDwB1lLLMcQ1Kf+s53vnPQ9re//e2GOukPiw8/hNMy846I+BS0TrZHxL7D7ZSZW4Df6VB/iNb5l5n1fwQ+eIjX+hzwuQ71jcDGuR5DUv9auXIld955J/v27WNgYIDzzjuv6ZYWtLnMXP4hIgYpJ8XLvSpPd7UrSapsdHSUgYEBABYvXswll1zScEcL21zC5RO0znv8dkR8F7gNWNPVriSpssHBQVatWkVEsGrVqpfc96K6DhsumflD4F8D/xK4DHhT+chLeonTTz/9oO2lS5c21In0UhdddBHHH388733ve5tuZcGby9VilwD/DjgbeButZVacT6qja6+99qDtz33uJafKpMZs2LCBPXv28JWvfKXpVha8uXws9rttP/8KuAa4qIs96Sj21FNPHbT99NOentP80L4q8qZNm1z3rsvm8rHYmrafP6J1BZirDKujT3/60wdtX3XVVQ11Ih1sfHyc/fv3A7Bv3z5XRe6yI7lDfw+tGxmll3juuedm3Zaa8rWvfY29e1tLDu7du5eJiYmGO1rYDnufS0R8hRfX5lpEa3n7O7rZlCTVtnLlSjZu3MjevXtZvHix97l02VxuomxfM30v8IvM3N6lfiSpK0ZHR9m0aRMAAwMD3ufSZXM553Jv2893DRZJRyPvc+mtQ85cIuIZOi9VH7S+C+ykrnUlSV0wOjrKww8/7KylBw4ZLpn5ml42IkndNjg4yLp165puoy/M5ZwLABFxOvDqA9uZ+UhXOpIkHfXmcof+RRGxDfg5cC/wMHBXl/uSpOqmp6e54oorvIGyB+Zyn8tnaX2D5N9m5pnAucB3u9qVJHXB+Pg4W7du9QbKHphLuDyfmdPAoohYlJnfAN7a5b4kqSqXf+mtuYTLUxFxIvBt4M8j4gZa97tI0lHD5V96ay7h8i3gZOBKWl8z/HeA61VLOqq4/EtvzSVcgtb31H8TOBH4YvmYTJKOGitXrmTx4tYFsi7/0n1zuUP/TzLzTcDlwG8B90bE17remSRVNDo6yqJFrf/lufxL972cVZF3An8PTAOnH2asJM0rLv/SW3NZFfljwL8FlgBfAv4oM3/S7cYkqTaXf+mdudyh/wbgjzPzR91uRpK6yeVfeuew4ZKZa3vRiCRp4TiSb6KUDunYY4+ddVtSfzBcVNWJJ5540PZrXuPi2lI/MlxU1cwlNZ544omGOpHUJMNFklSd4SJJqs5wkSRVN+dvopSkI7V+/XqmpqaaboMdO3YAMDQ01Ggfw8PDrFmzptEeuq1rM5eIOCMivhERD0bEAxFxZamfGhETEbGt/D6l1CMi1kXEVERsiYi3tb3WaBm/LSJG2+pnR8TWss+6iIjZjiGpvz377LM8++yzTbfRF7o5c9kL/KfM/GFEvAa4PyImgP8A3JOZ10XEWmAt8EngQmBF+XkHcBPwjog4FbgaGAGyvM6GzHyyjFkNfA/YCKyi9RXMaw9xDEkNmC/v0q+88koAbrjhhoY7Wfi6NnPJzMcz84fl8TPAg8AQcDEwXoaNA+8rjy8GbsuW7wEnR8TrgAuAiczcXQJlAlhVnjspMzdnZgK3zXitTseQJPVAT07oR8Ry4HeA7wNLM/NxaAUQL66wPAQ82rbb9lKbrb69Q51ZjjGzr9URMRkRk7t27TrSf54kaYauh0v5iuT/S2vxy1/NNrRDLY+gPmeZOZaZI5k5smTJkpezqyRpFl29WiwiXkUrWP48M/9fKf8yIl6XmY+Xj7Z2lvp24Iy23ZcBj5X6782of7PUl3UYP9sxFrT5ckXOTAc+5+61frgiR5qvunm1WAA3Aw9m5n9te2oDcOCKr1Hgy231S8pVY+cAT5ePtO4Gzo+IU8pVX+cDd5fnnomIc8qxLpnxWp2OIUnqgW7OXN4J/Htga0Qc+C6Yq4DrgDsi4lLgEeCD5bmNwHuAKWAP8FGAzNwdEZ8F7ivjPpOZu8vjjwG3AsfRukrsrlI/1DEWtPnwLv1DH/oQO3e+OFFcunSpV+ZIfahr4ZKZ36HzeRGAczuMT+DyQ7zWLcAtHeqTwJs71Kc7HUPdd9NNN/H+97//he0bb7yxwW4kNcXlX1TV4OAgAwMDQGvW4veUS/3JcFF1b3zjGznhhBOctUh9zHBRda961asYHh521iL1McNFklSd4SJJqs5wkSRVZ7hIkqozXCRJ1RkukqTqDBdJUnWGiySpOsNFklSd4SJJqs5wkSRVZ7hIkqozXCRJ1RkukqTqDBdJUnWGiySpOsNFklSd4SJJqs5wkSRVZ7hIkqozXCRJ1RkukqTqIjOb7mFeGBkZycnJySPef/369UxNTVXs6Oh14L/D8PBww500b3h4mDVr1jTag3+bL/Jv82A1/j4j4v7MHJlZX/yKXlUvmJqa4kc/fpB9x5/adCuNW/Sb1huW+x/6ZcOdNGtgz+6mWwBaf5vbHvhrXn/ivqZbadwxz7c+rHnuF0f+RnKheOTXA119fcOlon3Hn8qz//Q9TbeheeK4n25suoUXvP7EfVz1tl813YbmkWt/eFJXX99zLpKk6gwXSVJ1hoskqTrDRZJUXdfCJSJuiYidEfHjttqpETEREdvK71NKPSJiXURMRcSWiHhb2z6jZfy2iBhtq58dEVvLPusiImY7hiSpd7o5c7kVWDWjtha4JzNXAPeUbYALgRXlZzVwE7SCArgaeAfwduDqtrC4qYw9sN+qwxxDktQjXQuXzPwWMPNC/4uB8fJ4HHhfW/22bPkecHJEvA64AJjIzN2Z+SQwAawqz52UmZuzdRfobTNeq9MxJEk90utzLksz83GA8vv0Uh8CHm0bt73UZqtv71Cf7RgvERGrI2IyIiZ37dp1xP8oSdLB5stNlNGhlkdQf1kycwwYg9byLy93/3Y7duxgYM/T8+rGOTVrYM80O3bsbboNqRG9nrn8snykRfm9s9S3A2e0jVsGPHaY+rIO9dmOIUnqkV7PXDYAo8B15feX2+ofj4jbaZ28fzozH4+Iu4Fr207inw98KjN3R8QzEXEO8H3gEmD9YY7RVUNDQ/z9c4td/kUvOO6nGxkaWtp0G1IjuhYuEfEXwO8Bp0XEdlpXfV0H3BERlwKPAB8swzcC7wGmgD3ARwFKiHwWuK+M+0xmHrhI4GO0rkg7Drir/DDLMSRJPdK1cMnMjxziqXM7jE3g8kO8zi3ALR3qk8CbO9SnOx1D6lc7duzgH54Z6PpChTq6/OKZAU7YsaNrr+8d+pKk6ubL1WKSumRoaIjn9j7ukvs6yLU/PIljh4YOP/AIOXORJFVnuEiSqvNjsYoG9uz2Jkpg0T+2Pn7Z/+r+PoHc+ppjL0VWfzJcKhkeHm66hXljauoZAIb/Sb//j3WpfxfqW4ZLJWvWrGm6hXnjyiuvBOCGG25ouBNJTfGciySpOmcuUh945NfeRAnwyz2t99NLj9/fcCfNe+TXA6zo4usbLtIC53mfF/1magqAY9/gf5MVdPdvw3CRFjjPB77I84G94zkXSVJ1hoskqTrDRZJUneEiSarOcJEkVWe4SJKqM1wkSdUZLpKk6gwXSVJ1hoskqTrDRZJUneEiSarOhSsXkPXr1zNVVn1t0oEeDiwS2JTh4WEXbZQaYriouuOOO67pFjTP+MbnYP3wxsdwWUAW+h+r9Er5xqd3DBdJXecbn/7jCX1JUnWGiySpOsNFklSd4SJJqm7BhktErIqIn0XEVESsbbofSeonCzJcImIA+AJwIXAW8JGIOKvZriSpfyzIcAHeDkxl5kOZ+RvgduDihnuSpL6xUMNlCHi0bXt7qR0kIlZHxGRETO7atatnzUnSQrdQb6KMDrV8SSFzDBgDiIhdEfGLbjfWR04Dnmi6CakD/zbrekOn4kINl+3AGW3by4DHZtshM5d0taM+ExGTmTnSdB/STP5t9sZC/VjsPmBFRJwZEccAHwY2NNyTJPWNBTlzycy9EfFx4G5gALglMx9ouC1J6hsLMlwAMnMjsLHpPvrYWNMNSIfg32YPROZLznNLkvSKLNRzLpKkBhkukqTqDBdV5Zpumq8i4paI2BkRP266l35guKga13TTPHcrsKrpJvqF4aKaXNNN81ZmfgvY3XQf/cJwUU1zWtNN0sJnuKimOa3pJmnhM1xU08te003SwmS4qCbXdJMEGC6qKDP3AgfWdHsQuMM13TRfRMRfAJuBN0bE9oi4tOmeFjKXf5EkVefMRZJUneEiSarOcJEkVWe4SJKqM1wkSdUZLlLDImL5ka7U+0r2lbrJcJEkVWe4SPPD4ogYj4gtEfGliDg+Iv5LRNwXET+OiLGICICIODsi/iYiNgOXN9y31JHhIs0PbwTGMvMtwK+A/wj898z83cx8M3Ac8Adl7P8CrsjMf9FMq9LhGS7S/PBoZn63PP4z4F3A70fE9yNiK/Bu4E0R8Vrg5My8t4z93w30Kh3W4qYbkAS89KsJErgRGMnMRyPiGuDVtL7WwDWbNO85c5Hmh9dHxIGPuT4CfKc8fiIiTgQ+AJCZTwFPR8S7yvN/2Ns2pblx5iLNDw8CoxHxP4BtwE3AKcBW4GFaX2dwwEeBWyJiD60VqKV5x1WRJUnV+bGYJKk6w0WSVJ3hIkmqznCRJFVnuEiSqjNcJEnVGS6SpOr+P9XVbBW225AiAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import seaborn as sns\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "sns.boxplot(x = \"bad\", y = \"loan\", data = bank)\n",
    "plt.show()\n",
    "sns.boxplot(x = \"bad\", y = \"mortdue\", data = bank)\n",
    "plt.show()\n",
    "sns.boxplot(x = \"bad\", y = \"value\", data = bank)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "loan\n",
      "Statistics=0.836, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "mortdue\n",
      "Statistics=0.862, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "value\n",
      "Statistics=0.844, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "yoj\n",
      "Statistics=0.912, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "derog\n",
      "Statistics=0.277, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "delinq\n",
      "Statistics=0.380, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "clage\n",
      "Statistics=0.905, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "ninq\n",
      "Statistics=0.669, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "clno\n",
      "Statistics=0.964, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "debtinc\n",
      "Statistics=0.824, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# SHAPIRO\n",
    "from scipy.stats import shapiro\n",
    "def shapiroTest(columnName):\n",
    "    stat, p  = shapiro(bank[columnName]) \n",
    "    print(\"Statistics=%.3f, p=%.3f\" % (stat,p))\n",
    "    alpha = 0.05\n",
    "    if p > alpha:\n",
    "        print(\"Örneklem normal dağılımından gelmektedir\")\n",
    "    else:\n",
    "        print(\"Örneklem normal dağılımından gelmemektedir\")\n",
    "        \n",
    "for column in bank.columns:\n",
    "    if(column != 'reason' and column != 'job' and column != 'bad'):\n",
    "        print(column)\n",
    "        shapiroTest(column)\n",
    "        print()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "loan\n",
      "Statistics=2.395, p=0.122\n",
      "Örneklem normal dağılımından gelmektedir\n",
      "\n",
      "mortdue\n",
      "Statistics=2.307, p=0.129\n",
      "Örneklem normal dağılımından gelmektedir\n",
      "\n",
      "value\n",
      "Statistics=3.991, p=0.046\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "yoj\n",
      "Statistics=4.210, p=0.040\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "derog\n",
      "Statistics=118.586, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "delinq\n",
      "Statistics=172.844, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "clage\n",
      "Statistics=0.008, p=0.927\n",
      "Örneklem normal dağılımından gelmektedir\n",
      "\n",
      "ninq\n",
      "Statistics=19.275, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "clno\n",
      "Statistics=9.247, p=0.002\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n",
      "debtinc\n",
      "Statistics=115.704, p=0.000\n",
      "Örneklem normal dağılımından gelmemektedir\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# LEVENE\n",
    "from scipy.stats import levene\n",
    "def leveneTest(column):\n",
    "    stat, p = levene(bank[column][bank['bad'] == 0], bank[column][bank['bad'] == 1])\n",
    "    print(\"Statistics=%.3f, p=%.3f\" % (stat,p))\n",
    "    alpha = 0.05\n",
    "    if p > alpha:\n",
    "        print(\"Örneklem normal dağılımından gelmektedir\")\n",
    "    else:\n",
    "        print(\"Örneklem normal dağılımından gelmemektedir\")\n",
    "        \n",
    "for column in bank.columns:\n",
    "    if(column != 'reason' and column != 'job' and column != 'bad'):\n",
    "        print(column)\n",
    "        leveneTest(column)\n",
    "        print()"
   ]
  }
 ],
 "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.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}