Skip to content

Instantly share code, notes, and snippets.

@ryonsd
Created December 17, 2019 12:00
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save ryonsd/1bd0402800419ef88067b0a322ef5b6b to your computer and use it in GitHub Desktop.
Save ryonsd/1bd0402800419ef88067b0a322ef5b6b to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"from scipy import stats\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"t分布は,母集団が正規分布に従うときのt値の標本分布のこと"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# t値"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$$t=\\frac{\\hat{\\mu}-\\mu}{\\hat{\\sigma}/\\sqrt{N}}$$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"ここで,$\\hat{\\mu}$は標本平均,$\\mu$は母平均,$\\hat{\\sigma}$は不偏分散の平方根,$N$はサンプルサイズ"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"簡単に言うと,t値は標本平均の母平均からのズレを標準化(平均0,分散1に)したものである.(分散は1にはならないので厳密には標準化ではない)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# t分布"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"サンプルサイズが$N$のとき,$n=N-1$を自由度と呼ぶ."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"平均:0 \n",
"分散:$\\frac{n}{n-2} (n>2)$"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[]"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD4CAYAAADiry33AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3zV5d34/9d1TjbZi5AB2YEEwgoJEEbY4AB3Rf0qal3Faut9t9qf1tXaW4u19121rVZtragoWBUUGcqeSViBEEJCCGRPssdZ1++Pc8CAIIeQ5JNxPR+PPDzn+ozz/iC8z5X35/pcl5BSoiiKovRfOq0DUBRFUbqXSvSKoij9nEr0iqIo/ZxK9IqiKP2cSvSKoij9nIPWAVzI399fhoeHax2GoihKn7Jv375qKWXAxbb1ukQfHh5OZmam1mEoiqL0KUKIU5fapko3iqIo/ZxK9IqiKP2cSvSKoij9nF01eiHEfOD/AD3wjpTy5UvsdzOwCpggpcy0tf0GuB8wA49JKdd3ReCKogwcRqOR4uJi2tratA5Fcy4uLoSGhuLo6Gj3MZdN9EIIPfAmMAcoBjKEEKullEcv2M8DeBzY26EtHrgdSACCgW+FELFSSrPdESqKMuAVFxfj4eFBeHg4Qgitw9GMlJKamhqKi4uJiIiw+zh7SjfJQL6UskBKaQBWAIsust/vgFeAjl+5i4AVUsp2KeVJIN92PkVRFLu1tbXh5+c3oJM8gBACPz+/K/7Nxp5EHwIUdXhfbGvr+OHjgDAp5ddXeqzt+AeFEJlCiMyqqiq7AlcUZWAZ6En+rM78OVz1OHohhA54DVjS2XNIKd8G3gZISkpS8yYrvUphfSHfnPwGs63i6OHkwa2xt+Lm6KZxZIpiH3sSfQkQ1uF9qK3tLA9gJLDF9k0TBKwWQiy041hF6bXMFjMfHP2ANw6+Qbu5HYG1JyWRfHzsY16Y/AIpQ1I0jlJRLs+e0k0GECOEiBBCOGG9ubr67EYpZb2U0l9KGS6lDAf2AAtto25WA7cLIZyFEBFADJDe5VehKF2suLGYu9fdzZ/2/YnU4FQ237aZrHuyyLoni38v+DcOOgd+uuGn/H7P7zGajVqHqyg/6rI9eimlSQjxKLAe6/DK96SU2UKIF4FMKeXqHzk2WwjxKXAUMAFL1YgbpbdrMbbw6HePUtlayStTX2FBxILz6qJjA8ey8vqVvHHgDf591Jr0n0p+SsOIlb6kra2NadOm0d7ejslk4pZbbuGFF17o1s+0q0YvpVwLrL2g7dlL7Jt2wfuXgJc6GZ+i9CgpJb/b8zsK6gt4a85bTAqeBKUH4ODHcLaP4uKF66RH+dWEX2GRFpbnLGds4Fjmhc/TNnilT3B2dmbTpk24u7tjNBqZMmUKCxYsYOLEid32mb1uUjNF0dKqvFV8VfAVPxvzMyYFJcO2V2HL/4DOEZxsN19b6+DAh3DDX3li/BNkVWfx3K7niPOJI9wrXNP4le510003ER8fz7Zt2ygsLOS9995j9uzZV3QOIQTu7u6A9UEwo9HY7SOKVKJXFJucmhxe3vsyk4Mn89Cw6+Cf10DRHki4Ea77M7j6WHcsOwSfPQDLb8Ix5RH+NOUVbl27mCe2PsGH13yIq4OrthfSz72wJpujpQ1des74YE+euz7hsvsdPnyYyZMns23bNj7//HM+/PDD8xL91KlTaWxs/MFxr7766nn7mc1mxo8fT35+PkuXLiUlpXtv6qtEryhYSza/3/N7vJy9+J9JL6L76CdQnQc3vg2Jt0HHHteQ0fDQVtj4LOz9G0HO7rw89WUe/vZhlh9dzgOJD2h3IUq3aWlpob6+nl/+8peAtTfu7e193j7bt2+361x6vZ6DBw9SV1fHjTfeyJEjRxg5cmSXx3yWSvSKAmwp2kJWdRbPT3oe311vWOvyt30A8QsvfoCjK1yzDAwtsO1VUiOmkxaWxj+P/JPb4m7Dy9mrZy9gALGn590djh49yvjx49Hr9QBkZWX9IDnb26M/y9vbmxkzZrBu3TqV6BWlO1mkhdcPvs4wz2EsxAN2/QWS7rt0ku9owSvW8s5/HuTRxf/m1qKtvJ/9Po+Ne6z7A1d61OHDhxkzZsy591lZWSxadP5sMPb06KuqqnB0dMTb25vW1lY2btzIk08+2eXxdqSmKVYGvHUn15F3Jo+lw/8fjl8uhYARMO8P9h3s7A63vAct1cRt/TPzw+ezPGc51a3V3Ru00uMuTPSdLbeUlZUxY8YMEhMTmTBhAnPmzOG6667rylB/QPXolQHNaDHy5sE3ifWJZV72Rmirh7u/tJZm7DVkNMx+Adb/hqUxy9hg3sA7h99RY+v7mddee+289wUFBZ06T2JiIgcOHOiKkOymevTKgPZl/pecbjzNY5E3ozuyCiY+DIPjr/xEKQ9DwAiG7XmHG6IW8Wnup5Q2lXZ9wIrSCSrRKwOWlJL3s99npN9Iph37ztqLn9zJ2rpOB2lPQnUuD7sMRUrrfDiK0huoRK8MWBnlGRQ2FLI4ZAbiyGeQ/AAM8u/8CUcsgoARBO1+ixlhaXyR/wXt5vauC1hROknV6JUB65PcT/By9mLuiT0/2pv/8mAJ/9pViNlinUHby9WR314XT+xgj/N3PNurX7mE25xuYWN7HRsKN3B91PXdfSmK8qNUj14ZkKpbq9l0ehOLgqfhcuQ/F+3Nm8wWfvfVUR5fcZCWdjN+g5zwG+TE0dIGbnxzJ+uOlP/wxLZeffKBlQzzGMrK4yt76IoU5dJUj14ZkD7P+xyTNHFLZfFFe/Nnmg0s/Wg/u07UsGRyOE9fOwJHvbVfVF7fxsPL9/Hw8n38fGY0T8yJ/X6uEluvXrdyCbfGp/Fq0TqOnzlOrE9sT1+iopyjevTKgGO2mFl1fBUpgeOIOLYeRi8+rzcvpeSJTw+SWXiGV28dzfMLE84leYAgLxc+eWgityWF8vqmfJbvPX3+B4xYCF5hLCw6hpPOiZW5qlevaEslemXA2Vm6k9LmUm51CARTG4y/57ztqw+Vsjm3iicXDOeW8aEXPYezg55Xbk4kNdqPV745Rll96/cbdXoYexc+J7cxN3gyawrW0GJs6c5LUpQfpRK9MuCsOr4KPxc/ZubvtD7sNGT0uW21zQZeWHOUMWHeLJkc/qPnEULwPzcmYrJY+O0XR5Cyw3LHY+8CoeM2g45mYzPrCtd109UofdF9991HYGDgD56sXbduHXFxcURHR/Pyyy932efZleiFEPOFELlCiHwhxA8e9xNCPCyEOCyEOCiE2CGEiLe1hwshWm3tB4UQf++yyBWlE+rb69lesp1rA5NwrMiGcef35l9ck01jm5FXbk5Er7v8HOFD/dz4rzlxfJtTyVdZZd9v8AqF6NmMObqBoR5hrD259tInUQacJUuWsG7d+V/+ZrOZpUuX8s0333D06FE+/vhjjh492iWfd9lEL4TQA28CC4B4YPHZRN7BR1LKUVLKMcAfgY7PCp+QUo6x/TzcJVErSidtOr0Jk8XEgjPV4OgGo249t21LbiVfHCzlkbRo4oI8fuQs57s3NZxRIV48vzqbuhbD9xvG3Y1oLGW+ZywZ5Rlq/pt+4KabbuKZZ55h2rRpDB06lG+//bZT55k2bRq+vr7ntaWnpxMdHU1kZCROTk7cfvvtfPnll10Rtl2jbpKBfCllAYAQYgWwCOs6sABIKTuuAjAIkChKL7SucB2hg0JIOLbRuqCIiydgvQH76oZcwv3cWDoj6orO6aDX8fLNo7j2Lzv4165CfjHbNsImdj4MCmRBxSnelhY2ntrI4uGLu/qSBp5vnoLyw117zqBRsODypZKuWnjkYkpKSggLCzv3PjQ0lL17917BRVyaPYk+BCjq8L4Y+MFyKEKIpcATgBMws8OmCCHEAaABeEZK+YN5PIUQDwIPAgwdOtTu4BXlStS21bK3bC/3+ichDLvPK9vsKajlSEkDf7hxFM4O+is+d0KwFzOHB/LB7lM8PD0KF0c96B1hzB1E73qd6MRU1p1cpxJ9H9aVC4/0tC4bRy+lfBN4UwhxB/AMcA9QBgyVUtYIIcYDXwghEi74DQAp5dvA2wBJSUnqtwGlW3x76lvM0sz8sgLwj4Ow5HPb3tlegN8gJ24aF9Lp8z8wNZLF/9jD5wdKWJxs67CMuxt2/i/zHXx5o3I/5c3lBA0KutpLGdjs6Hl3h+5YeKSjkJAQioq+71MXFxcTEtL5v48d2ZPoS4CwDu9DbW2XsgL4G4CUsh1ot73eJ4Q4AcQCmZ2KVlGuwjcnvyHCPYzYwzthxjPnlgfMr2zku2OV/GJ2jLUn3kkTI30ZGeLJP7YX8JOkMHQ6AX5REDyO+WUneMMZ1heu556Eey5/MqXX6aqFRy5lwoQJ5OXlcfLkSUJCQlixYgUfffRRp8/XkT2jbjKAGCFEhBDCCbgdWN1xByFETIe31wJ5tvYA281chBCRQAzQuUmcFeUqVLZUsq9iHwucgxAACTec2/bujpM4O+j4fxOHXdVnCCF4YGokBVXNbM6t/H5Dwo0MKz1MvFc06wvXX9VnKNrpqoVHABYvXsykSZPIzc0lNDSUd999FwcHB9544w3mzZvHiBEjuO2220hI6JplEy/bo5dSmoQQjwLrAT3wnpQyWwjxIpAppVwNPCqEmA0YgTNYyzYA04AXhRBGwAI8LKWs7ZLIFeUKbCjcgEQyr6IQBo8Ef2vfpKqxnc/2l3DL+FD83J2v+nOuGTWEV745xtvbCpg1YrC1MX4RbPwt8/U+vFadQVFjEWEeYT9+IqXX6aqFRwA+/vjiU1hfc801XHPNNZ0+76XYNY5eSrlWShkrpYySUr5ka3vWluSRUj4upUywDaGcIaXMtrV/1qF9nJRyTZdfgaLYYX3hemI9I4gs2gfx3/fml+85hcFk4f4pEV3yOY56HfemRrD3ZC2Hi+utjT7DIHgc88ryzsWiKD1JPRmr9Hs1rTUcqjrEbMcAa4OtbGOxSFbtK2Z6bABRAe5d9nk/SQ7DxVHHp5kdBqsl3EBwaRYJ3tFsKdrSZZ+lKPZQiV7p97YVb0MiSas4CYEJ58o2e0/WUlLXelUjbS7G08WROfFBrMkqxWCyWBvjrTftpgtPsqqyqGmt6dLPVJQfoxK90u9tKdrCYFd/hhftO+8m7OcHihnkpGdufNcPd7xxbDB1LUa2nL0p6xMOwWOZUXECiWRb8bYu/0xFuRSV6JV+rd3czu6y3aSdHW1jq8+3Gc18c7ic+SOH4OrU+SGVlzI1JgC/QU58cbDDSOSEG4krOUyQqz9bi7d2+WcqyqWoRK/0a3vL9tJqaiWtptxatgmwTk/wbU4Fje2mLi/bnOWo13H96GC+zamkvtVobYxfhACmOwawq3SXWk9W6TEq0Sv92tairbjqXZlQfBBGXHeu/fP9JQz2dGZipF+3ffaNY0MwmCysPWyb1dInHIJGkXamilZTK+ll6d322YrSkUr0Sr8lpWRL8RZS3YfhbLFA7DwAapra2Xq8ihvGhNg1FXFnJYZ6ERkwiM/3dyjfxM5nQvFBXPWuavSN0mNUolf6rZzaHCpbKpne2gbug2HIWAC+yirDZJHc2E1lm7OEENw0NoT0wlqKam0rTMXOx9liYbL7MLYUbzl/sRJlwAgPD2fUqFGMGTOGpKSkc+2aLjyiKH3RlqItCATTTmdBzFzrwt3AlwdLGB7kwfAgz26PYdEY65fJmqxSa0PwOBgUQFprO5UtleTU5nR7DErvtHnzZg4ePEhmpnXqL00XHlGUvmpL0RZGe0bg21pnnRseqGxoY//pOq4ZNaRHYgjzdSMx1Iv12RXWBp0OYuYxtegQAsHWIjX6pi/pqoVHLkbrhUcUpc+pbq0mpzaHx9yHg94JItMA2HDUmnDnJfTcVMHzEoJYtj6X8vo2grxcIHYefgeXM8ojkR2lO3hkzCM9Fkt/8Er6KxyrPdal5xzuO5wnk5+87H5dtfCIEIK5c+cihOChhx7iwQcf1HzhEUXpc3aX7gZgclk+REwDZ+sUB+uzywn3cyN2cNdNeXA58xIGs2x9LhuOlnP3pHCImgF6JyZbnHi7+gj17fV4OXv1WDxK53TlwiM7duwgJCSEyspK5syZw/Dhw7s83o5Uolf6pZ2lO/F18mRE9RFI/hkA9a1Gdp+o4b4pEQjRfaNtLhQd6EFkwCA2ZFdYE72zB4RPIbXyJH8fZGF32W7mh8/vsXj6Ont63t2hKxceObugSGBgIDfeeCPp6emkpqZ228Ijqkav9DsWaWF36W4mOgda/4LbhlVuya3EZJHMSxjc4zHNSwhiT0EN9S22h6di5zOy8gQeDm7sKtnV4/EoV+5iC48kJiaet8/27ds5ePDgD346Jvnm5uZzXwbNzc1s2LCBkSNHnrfwiMFgYMWKFSxcuLBLYleJXul3jtUeo7atltSGM9anYb2ty/qtzy4nwMOZsWE+PR7TvIQgTBbJd8dsN2Vj5+EATHQJYmfpTjXMsg/oqoVHKioqmDJlCqNHjyY5OZlrr72W+fPna7vwiKL0NbtKrT3kycXZkGIt27QZzWzJreKGsSHWJf56WGKIF4M9nVmfXc5N40KtT8n6xzG5qZGNVHGi7gTRPtE9Hpdiv65aeCQyMpJDhw5ddJumC48IIeYLIXKFEPlCiKcusv1hIcRhIcRBIcQOIUR8h22/sR2XK4SY15XBK8rF7CzZSZzbEPxNBoieZW3Lr6bFYO7R0TYd6XSCufFBbD1eRavBbG2MnkVqiXWc9M7SnZrEpQwMl030tjVf3wQWAPHA4o6J3OYjKeUoKeUY4I/Aa7Zj47GuMZsAzAf+enYNWUXpDs3GZg5WHiTV4gyOgyBsImAt23g4OzCpG+e2uZx5CUG0GS1sz6uyNkTNYoihlUjXwed+C1GU7mBPjz4ZyJdSFkgpDcAK4Lylz6WUDR3eDgLOFhwXASuklO1SypNAvu18itIt0svSMUkTqdWnIWIqODhhsUg2HaskbXggTg7a3ZZKifTF08WBb3Nsdfphk0HvzGThRmZ5Jq2mVs1i6wvUfQyrzvw52PO3PgTosCYaxba28wghlgohTmDt0T92hcc+KITIFEJkVlVV2Ru7ovzAztKduOpdGFtVCFHWsk1WST3VTQZmDQ/UNDZHvY5psQFszq3CYpHg5AbDJpNaXYLBYmBfxT5N4+vNXFxcqKmpGfDJXkpJTU0NLi4uV3Rcl92MlVK+CbwphLgDeAa45wqOfRt4GyApKWlg/59UrsrOkp0kuwbhyPFz9fnNxyoRAqbFBmgcHcyIC+SrrDKOljUwMsQLomcxfuNvcYqMZGfJTqaETNE6xF4pNDSU4uJiVEfQ+qUXGhp6RcfYk+hLgLAO70NtbZeyAvhbJ49VlE4rbiymuKmYu/RDwHsY+EYCsDm3krFh3vgOctI4QpgeF4AQsOlYpTXRR83CdcMzjHMNZk/ZHq3D67UcHR2JiIjQOow+y57STQYQI4SIEEI4Yb25urrjDkKImA5vrwXybK9XA7cLIZyFEBFADKBWW1C6xd4y67wgE0tzrL15IahqbCeruJ4ZcdqWbc7yd3cmMdSbzWfXkg0cAR7BpLQbyK/Lp7q1WtsAlX7psoleSmkCHgXWAznAp1LKbCHEi0KIs49tPSqEyBZCHASewFa2kVJmA58CR4F1wFIppbkbrkNR2FO2hwAnLyJbGs7V588uzj1D4/p8RzPjAjlYVEdNUzsIAVEzmVRu7Rud/bJSlK5k1xAEKeVaKWWslDJKSvmSre1ZKeVq2+vHpZQJUsoxUsoZtgR/9tiXbMfFSSm/6Z7LUAY6i7SQXp5Oit4LoXOwTmQGbMmtItDDmYTg7p973l4zhgcgJWw7O8wyeibDG2vxcHBTiV7pFmoKBKVfyDuTR21bLSkN1RCaDC6eGM0Wth2vYkZcYI9OYnY5I4O98Hd3ZtMxW6KPnIEeQYqjH3vK9gz4kSVK11OJXukXzt7InFiWB1EzAdh36gyN7aZeVbYB61OyaXEBbDtehclsATdfCBlHSlM9Zc1lFDUWXf4kinIFVKJX+oW9ZXsJd/YjyGy2zveOdVilo14wJcZf4+h+aObwQOpbjRwoqrM2RKZZv6RAjb5RupxK9EqfZzQbyazIJAVncPaCIdYZBjfnVpIc4Yu7c++bu29KjD8OOsGmY7bRN5FpDDO2M9jJSyV6pcupRK/0eYerD9NqamVidal12gO9A6V1rRyvaCIttneVbc7ydHFk3DAfth231enDUhAOrkwU7mSUZ2CRFm0DVPoVleiVPm9v2V4Eggm1xefWhj2bQKfHaf807KVMjw0gu7SBqsZ2cHCGYZNJqaukrr2O3NpcrcNT+hGV6JU+b0/ZHuJdAvGyWCDSWp/fllfFEC8XYgJ7bm3YKzXdNiXDudksI9OYWHUSUHV6pWupRK/0aS3GFrKqskgxSvAMBb8oTGYL2/OqmRYT0KuGVV4ofogn/u5ObD3+faIPMFuIcvZX4+mVLqUSvdKn7a/cj0maSKk8CVFpIAQHi+pobDP16rINWIdZTosJYHtetXU2y8Ejwc2fZIsD+yv3YzQbtQ5R6SdUolf6tPTydByEnjGNtd+XbY5XoROQGtX7hlVeaFpsALXNBo6U1oNOB5HTSakuotXUypGaI1qHp/QTKtErfVp6WTqJzv64SXlu2oOtx6sYO9QHLzdHjaO7vKkx/ggBW3O/L98k1VUgEKp8o3QZleiVPqvB0EBObQ7JbW3Wsod7ILXNBrJK6s/d6Ozt/NydGRXidV6d3stiYbizH+nlaqJXpWuoRK/0WfvK92GRFpIrTp4bVrk9rwope8ciI/aaFhPAgaI66luN4D0UfKOYYLBwqPIQbaY2rcNT+gGV6JU+K708HWedI6Nbms4l+q3Hq/Bxc2RUiJemsV2J6XEBmC2SXfm2uegjp5NSdQqDxcChqkPaBqf0CyrRK31Wenk6Yxy8cdI5wtBJWCySbcermRITgF7Xe4dVXmhsmDceLg7nlW/GNdWhFzpVp1e6hEr0Sp9U21bL8TPHSWlpgtAJ4OxOTnkD1U3tfaY+f5aDXkdqlD/bjldZpygOn4q7hARHXzLKM7QOT+kH7Er0Qoj5QohcIUS+EOKpi2x/QghxVAiRJYT4TggxrMM2sxDioO1n9YXHKkpnZJZnAjChshAipwOc6xFP64WzVV7OtNgASuvbyK9ssk5bPCSR5PZ2jlQfodnYrHV4Sh932UQvhNADbwILgHhgsRAi/oLdDgBJUspEYBXwxw7bWm0rT42RUi5EUbpAenk6bjonEtrbz5vfZniQB4GeLprG1hnTYq1fTh3LNxMqCzFJE/sr9msXmNIv2NOjTwbypZQFUkoDsAJY1HEHKeVmKWWL7e0eILRrw1SU86WXpzNe74mjkzuEjKep3cS+U2d6/dOwlxLq40ZUwCC25dluyEZMZ2xrCw5Cr8o3ylWzJ9GHAB2XvCm2tV3K/UDHtWFdhBCZQog9QogbLnaAEOJB2z6ZVVVVdoSkDGRVLVWcrD9JcuMZGDYZ9I7sPlGD0SyZHtM3Ez3A9NhA9hbU0GY0w9BJuOocSXT0Zm+5uiGrXJ0uvRkrhLgLSAKWdWgeJqVMAu4A/lcIEXXhcVLKt6WUSVLKpICAvvsPVekZZx8kSq45f1piV0c948N9tAvsKk2L9afdZGHvyVpwcoOwFFKamzhWe4z69nqtw1P6MHsSfQkQ1uF9qK3tPEKI2cDTwEIpZfvZdillie2/BcAWYOxVxKsopJen46F3Ic5ghAjrjdhteVVMjvLD2UGvcXSdlxLhh5OD7vvFSCKmM6G6CIu0sK9in7bBKX2aPYk+A4gRQkQIIZyA24HzRs8IIcYCb2FN8pUd2n2EEM621/5AKnC0q4JXBqb0snSShCt6N38IjKewuplTNS196mnYi3F10pMS4dvhhux0Rre14SwcVJ1euSqXTfRSShPwKLAeyAE+lVJmCyFeFEKcHUWzDHAHVl4wjHIEkCmEOARsBl6WUqpEr3RaaVMpxU3FJNdVWYdV6nRssy3c0dcTPVgXI8mvbKK0rhWCx+Hk5MEYvYea90a5KnatmiylXAusvaDt2Q6vZ1/iuF3AqKsJUFE6Olefr6uCibayzfEqhvq6Ee7npmVoXWJabAB8ncO241XcnjwUwqeQ3HCU101nONN2Bh+XvnsPQtGOejJW6VMyyjPw0bsQbTRCZBoGk4VdJ2qYFuvfq1eTsldMoDtDvFzOG0+ffKYMQJVvlE5TiV7pM6SUpJenk2RxQucTAT7DyCyspcVgZnpsoNbhdQkhBNNjA9iRX43JbIHINBLaDbjqHFX5Ruk0leiVPqOosYjy5nJSzpSdG1a55XgVTnodk6P8NI2tK6XFBdDYZmL/6ToIiMPRPYhxuKoevdJpKtErfcbZHu2EpvrvE31uJRMifBjkbNftpj5hcrQ/DjrBltxKEAIi00ipr6agvoDq1mqtw1P6IJXolT4jvTwdf70LEUYzREyjtK6V4xVNpPWTss1Zni6OjBvmw5YOywsmN9QCqk6vdI5K9EqfIKUkvSydZKNADBkNbr7nEmFfnd/mx6TFBXC0rIHKhjaInM5wgwEPnZOan17pFJXolT7hZP1JatpqrCNQItMAa9km2MuFmEB3TWPrDmd/S9lyvAo8g9H7xzHe4qh69EqnqESv9AlnJ/ZKbmk+N6xyZ3410+MC+8WwyguNGOJBoIczWzuUb1LOVHC68TRlTWWaxqb0PSrRK31Celk6wTpXQqUDDJ1I5qlamg1m0vph2Qa+H2a5Pa/q3DDL5OZGADXMUrliKtErvZ7ZYia9PJ2UdhNi6ERwdGXr8Soc9YLU6L63mpS90uICaWgzcbCoDsJTiTFZ8NU5qzq9csVUold6vdwzuTQYGs6rz2/NrSJpmC/u/WhY5YWmxPij1wnrTWcXL0TIeJKN1jKWlFLr8JQ+RCV6pddLL7OWKlLa2iAyjbL6Vo6VN/bbss1ZXq6OjBvqzeZc24SwUTNIPlNOZUslpxpOaRuc0qeoRK/0envK9xCpcyXAyRuGjGHTMWvimzm8f42fv5gZwwPJLm2gvL4NomYysbUVQJVvlCuiEr3Sq9SYo4sAACAASURBVBnNRvZX7Ce5uQmiZoBOx6acSsJ8XYnuh8MqLzRr+GAAa68+ZDyhejeG6JzV8oLKFVGJXunVDlcfptXUysSGWoiaSavBzI78amYNH9wvh1VeKHawOyHernyXUwl6R0TENJJb28koz8AiLVqHp/QRKtErvdre8r0IIKmtHSJnsLugmnaTZUCUbcA6zHLWiEB25ldbFw2PmkFKQzV17XUcP3Nc6/CUPsKuRC+EmC+EyBVC5AshnrrI9ieEEEeFEFlCiO+EEMM6bLtHCJFn+7mnK4NX+r/0snSG44yXXyx4hfBdTiVuTnpSIn21Dq3HzBweSKvRzO6CGoiaSXKrdUlmVadX7HXZRC+E0ANvAguAeGCxECL+gt0OAElSykRgFfBH27G+wHNACpAMPCeEUEvkKHZpNbVyqOrQubKNlJJNxyqZGuPfpxcBv1ITI/1wddSzKacSfCMZ7DmUcNS8N4r97OnRJwP5UsoCKaUBWAEs6riDlHKzlLLF9nYPEGp7PQ/YKKWslVKeATYC87smdKW/O1BxAKPFaJ32IGomOWWNlNW3nbtBOVC4OOqZEuPPpmOV1vHzUTNJaWpgX8U+jBaj1uEpfYA9iT4EKOrwvtjWdin3A99cybFCiAeFEJlCiMyqqio7QlIGgj1le3BExzijhGGT2XSsAoC04f17/PzFzBoeSEldK7kVjdZhls2NtJhaOFx1WOvQlD6gS2/GCiHuApKAZVdynJTybSllkpQyKSBg4P0jVi5ud9luxpgFbmEp4DSI745VMjrUi0APF61D63EzbDefv8uphIipTGg3okOwu2y3xpEpfYE9ib4ECOvwPtTWdh4hxGzgaWChlLL9So5VlAvVtNZwrPYYE+utNyCrm9o5WFTHzAFWtjlrsKcLo0K8+C6nAly88ApOIsGiZ3epSvTK5dmT6DOAGCFEhBDCCbgdWN1xByHEWOAtrEm+ssOm9cBcIYSP7SbsXFubovyoszM0Tmq1PhH6XU4FUsKsEQNjWOXFzB4xmANFdVQ1tlvLNw01HKk+TKOhUevQlF7usoleSmkCHsWaoHOAT6WU2UKIF4UQC227LQPcgZVCiINCiNW2Y2uB32H9ssgAXrS1KcqP2l26Gw/0xDv5QlAi67MrCPF2JSHYU+vQNDM3YTBSwsajFRAzh0mtbZilRS1GolyWXVP/SSnXAmsvaHu2w+vZP3Lse8B7nQ1QGXiklOwu3U1Kaxv66Nk02Z6GvStl2IB4GvZShgd5MNTXjQ1Hy7ljQhKj9Z64omN36W5mDp2pdXhKL6aejFV6nVMNpyhvKWdScyPEzGFrbhUGk4V5CQOzPn+WEIJ5CYPZlV9Do8GMU/Rsxrcb2FO2R+vQlF5OJXql1zk7kmRimwEi01ifXY7vICeSwgfO07CXMjchCIPZwubcKoiezcSmRgobCilvLtc6NKUXU4le6XV2l+4mxCIIC07C4OjJ5mOVzB4RiF43cMs2Z40b6oO/uxPrs8shaiaT2qwD3NToG+XHqESv9Comi4mMsr1MbGpERM9hd0ENje0m5iUEaR1ar6DXCebED2bLsUranbyIGTwGP6nG0ys/TiV6pVc5Un2EJlMLE9vaIGYu67PLcXPS9+u1Ya/U3Pggmg1mduXXIGLmMrG5ib2lu9W0xcolqUSv9Cq7SnchgIl6bywB8Ww8WkFaXAAujgNnErPLmRzth7uzg7V8Ez2Hya1t1LbXcaz2mNahKb2USvRKr7KzZDujDCa8o2dzoNj6cJAq25zP2UFPWlwAG49WYAocyWRhXWlrZ8lOjSNTeiuV6JVeo66tjsPV2aQ2N0P0HNYcKsPJQTdgFhm5EteMGkJNs4G9hXX4R81mhNHMjpLtWoel9FIq0Su9xq7SXUgkU9qMWCKms/ZwGWmxAXi4OGodWq8zIy4QNyc9X2WVQsxcpjQ3cajyEA2GBq1DU3ohleiVXmNnyQ68LJAQPJGMMhOVje1cNzpY67B6JVcnPbNHDOabI+UYI2aQ2m7CjEUtRqJclEr0Sq9gkRZ2Fm9ncksz+uHX8FVWGS6OOmapss0lXZc4hLoWIzuL2kkcMhF3CTuLd2gdltILqUSv9Aq5tbnUGOpIbW3DFD2Pb46UMWv4YAY52zUd04A0PS4AD2cHvsoqw3H4NUxsaWFH8VbrKlSK0oFK9EqvsLPUOmIk1T2CvbVuVDcZuC5xiMZR9W7ODnrmJAxmfXY57VFzmdLSSkVbDSfqTmgdmtLLqESv9ArbT29mRLsB/7jr+CqrlEFO+nOrKimXdn1iMI1tJraXO5PqEQnAjhJVvlHOpxK9orlGQyOHqg+T2tqKMWYB3xwpZ3b8YPWQlB1So/3xcnXkq6xSguKuJ9pgZMfpzVqHpfQyKtErmttbthczklThzo4m6w3G6xLVaBt7ODnomJ8QxMajFbRFzSO1tZX91YdoMbZoHZrSi9iV6IUQ84UQuUKIfCHEUxfZPk0IsV8IYRJC3HLBNrNt1alzK08pSkdbT2/Cw2JhdOQ8/nOgFC9XR6bFqrlt7LVobDDNBjPrq/2ZijtGaVaTnCnnuWyiF0LogTeBBUA8sFgIEX/BbqeBJcBHFzlFq5RyjO1n4UW2KwOY2WJm2+lNTGlpxRA5nw3Z5SwcHYyzgyrb2GtihB8h3q58dqCUcZHz8LBY2HLqO63DUnoRe3r0yUC+lLJASmkAVgCLOu4gpSyUUmYBavo85Yocrj5MramZ6Ub4uj6KdpOFm8eHah1Wn6LTCW4aF8KOvCqahs1jSksr205vUrNZKufYk+hDgKIO74ttbfZyEUJkCiH2CCFuuNgOQogHbftkVlVVXcGplb5u6+lN6KVkStgMVh6qJCpgEKNDvbQOq8+5aVwoFgmfVocz3aSn1tTM4erDWoel9BI9cTN2mJQyCbgD+F8hRNSFO0gp35ZSJkkpkwICAnogJKW32HJyHePa2mkfei0ZhWe4eXzogF4AvLMi/AcxfpgPqw6UM2XoLPRSslWVbxQbexJ9CRDW4X2orc0uUsoS238LgC3A2CuIT+nHihuLyW8pY7rBwse1MQgBN469kl8WlY5uGhdCXmUTNUMWMLatnS0F32gdktJL2JPoM4AYIUSEEMIJuB2wa/SMEMJHCOFse+0PpAJHOxus0r9sPb0JgOlBk1h5sIrUKH+GeLlqHFXfdV1iME4OOj4oH0aaEfJayylpsrtPpvRjl030UkoT8CiwHsgBPpVSZgshXhRCLAQQQkwQQhQDtwJvCSGybYePADKFEIeAzcDLUkqV6BUAtuSvJtxgpD3wWorPtHLzeNWbvxpero7MiR/MF1mVTB0yBYAtp77VOCqlN7Brxigp5Vpg7QVtz3Z4nYG1pHPhcbuAUVcZo9IPNRmayKw7zl3tJt4pjcDDuV6tJNUFbh0fytdZZVR5zCX8zG625q3mzoR7tA5L0Zh6MlbRxM7i7ZiQpPqO5svsM9w0LgQ3JzVT5dWaFhNAqI8rfykMYYbBQkZ9Hk2GJq3DUjSmEr2iie+OfYK32cwZhzkYzBbunDhM65D6BZ1OcEfKUHaebGSC9xhMSLaruW8GPJXolR7Xbm5na9VBZrYa+XNhOMnhvsQO9tA6rH7jtqQwHPWC05ZZBJhMbMz5WOuQFI2pRK/0uN2nt9KCmQmDRnK81sydE4dqHVK/4u/uzPyRQ/jTiVBmtlvYUZutJjkb4FSiV3rcxuwP8DBbyG2fje8gJ+aPVDdhu9qdKUOpbYME17G0YmFn4UatQ1I0pBK90qOMZiOba7JIM0jeOjWUW5NC1QRm3SAlwpfoQHd21s/C12xmY/aHWoekaEgleqVH7T31LY1YiNSPwmDRcUeyKtt0ByEEd6YMZVX5EKYb9WytP0a7uV3rsBSNqESv9KiNh99nkMXChtLpzB4xmGF+g7QOqd+6ZXwoHs6O+JpH04JkV94arUNSNKISvdJjTBYTm87kMNmgJ711KA9MjdA6pH7Nw8WR25PD+LJ4Op6qfDOgqUSv9JjMvK+oExYcm+MZHepNcoSv1iH1e/emRlBECMkGF7Y05mM0G7UOSdGASvRKj1l/5H1cLRa2Vc3mgWmRajriHhDs7cp1iUNoq0ugUcAu1asfkFSiV3qEwdjK+oZ8xre5oPOKZL6a16bH/HRqJDvrFuBttvC1SvQDkkr0So/Yvu9vNOqgrmY890+JwEGv/ur1lJEhXoyMimREixeb2stoaq7UOiSlh6l/bUqP+CrvM3zNFo6ZruO2CWGXP0DpUg9Mi6S0ZirtQvDdnte0DkfpYSrRK92uvvo4W831hDUEsmTKCNyd1SyVPS0tNgDhey1BRslXp9VTsgONSvRKt9u4ZxlGIahsmc+9U8K1DmdAEkLwxJzh+DYMY69op7Jwm9YhKT3IrkQvhJgvhMgVQuQLIZ66yPZpQoj9QgiTEOKWC7bdI4TIs/2oFRAGGouFNWW7CTLouD7lRjxdHLWOaMBKiwvA4vwTpBB8nfEXrcNRetBlE70QQg+8CSwA4oHFQoj4C3Y7DSwBPrrgWF/gOSAFSAaeE0L4XH3YSl9RenQV+50EHs3DWZKqHpDSkhCCpbPnE9rmyNd12WBQM1oOFPb06JOBfCllgZTSAKwAFnXcQUpZKKXMAiwXHDsP2CilrJVSngE2AvO7IG6lj1iV8VcAUkY8gofqzWtuemwAPpaJ5Do5kLP3r1qHo/QQexJ9CFDU4X2xrc0edh0rhHhQCJEphMisqqqy89RKb2euOs46cznBrR4sTZuidTgKtsnOpv4CvYSVRz4AKbUOSekBveJmrJTybSllkpQyKSAgQOtwlC7y9brnKHJ0ZELI7WqkTS9yTUIMEcYw1jsbqT7ytdbhKD3AnkRfAnQc+Bxqa7PH1Ryr9GGGxho2NWTiZtbx1LwHtA5H6UAIwR0Tf0GDXs8XW1/WOhylB9iT6DOAGCFEhBDCCbgdWG3n+dcDc4UQPrabsHNtbUo/t/2Ll9k6yJlU/zTcnVy1Dke5wC0j5+BvcWO7YwVFxzK1DkfpZpdN9FJKE/Ao1gSdA3wqpcwWQrwohFgIIISYIIQoBm4F3hJCZNuOrQV+h/XLIgN40dam9GN1jU0cqf4SkxA8Ou1xrcNRLkIIwc0j7mC/iwu7v/m91uEo3cyuGr2Ucq2UMlZKGSWlfMnW9qyUcrXtdYaUMlRKOUhK6SelTOhw7HtSymjbzz+75zKU3uTbVX9nnYeOUW7hRHpHah2Ocgl3jPl/6KUgX5/FnqwcrcNRulGvuBmr9B+HTtdiLv8XxY6O3DnuIa3DUX6Er4svM4dM4it3V3LWvEKb0ax1SEo3UYle6TJGs4U1K95ij1crXnpX5oTP1Tok5TJuT7yfRr0ON6fNvL0uQ+twlG6iEr3SZd7ZdoIpbR/ynZsbtwxfjJPeSeuQlMuYEDSBaPcwPvF2QZ/+JkdLG7QOSekGKtErXaKwupljm5az06cJvU7HnfF3aR2SYgchBPeMfpA8JydGuG/hD6u2Y7aoh6j6G5Xolatmtkh+89lB7nH4jC88PLgu8noC3NSDb33FNRHXEODswyeezkyuXMG7Owq0DknpYirRK1ft7W0F+J5aR6ZnPW0C7k5Qk5T2JU56J+5IuJtdbq5MdfuOd9dncKSkXuuwlC6kEr1yVQ4V1fF/G7J5ctAXfOTtTWpwKjE+MVqHpVyhW2NvxVXvwsdeLjzuvIbHVxygxWDSOiyli6hEr3Rac7uJX3xykAfdtpDpXEONDpaMXKJ1WEoneDl7cVPszXzj7s4M3QZkTT6/+0qNre8vVKJXOkVKyXOrs6mrKednulX8MyCI4T7DSQlK0To0pZPuGnEXFgT/9vHmrcDP+Tj9NN8cLtM6LKULqESvdMryPadYta+Y98I3sc7ZTCFGHh79MEIIrUNTOinUI5Tro67nU49BeDbu5u7BBfz3ykMcr2jUOjTlKqlEr1yxPQU1vLDmKHdFtpJQvoq/BwYT7xfPzKEztQ5NuUoPj34YixD8IzCE3zosx91R8OC/M6lvMWodmnIVVKJXrkjxmRZ+9uF+hvq68rzzh3zh7UOJpY2fj/256s33A6EeodwcczOfueqpqM9j5YQcSupaefTj/ZjMFy4gp/QVKtErdmtoM/LAv/dhNFv4eOJpTIWbecs/gLGBY0kNTtU6PKWLPDDqAXRCz1thcQzd/yqvzvVne141v/86B6lWpOqTVKJX7NJmNPPA+5nkVTTy1o1DGbzzOVaGxVNpala9+X5m8KDB/GT4T1hNC4V6yaLiZdw3OZx/7Srkr1tOaB2e0gkq0SuXZTJbeOzjA+w9WcufbhvN5NyXqTe18LYLpAxJYULQBK1DVLrY/SPvx9nBmT9Hj4e8DTwz9DA3jAlm2fpcVqSf1jo85QqpRK/8KItF8vTnR9hwtILnro9nkfN+OPoFr4+YRoOpmV8l/UrrEJVu4Ofqx4OJD7Kp6SQ7wkajW/8UyxYMIS0ugP/v88Nq2GUfY1eiF0LMF0LkCiHyhRBPXWS7sxDiE9v2vUKIcFt7uBCiVQhx0Pbz964NX+lOZovkyc+y+CSziMdmRnNvoht89QQ5Q+JZ2ZTH7XG3E+cbp3WYSje5O/5uwj3DednTBYOhGce1v+Svd4xl7FAfHv34AGsOlWodomKnyyZ6IYQeeBNYAMQDi4UQ8Rfsdj9wRkoZDfwZeKXDthNSyjG2n4e7KG6lmxnNFn75yUFW7ivm8Vkx/HJWFPznp1jaG3kpaAjezt4sHbtU6zCVbuSkd+Kp5Kc41VLGv8ctgty1uO1/m/fvS2b8MB8eX3GAVfuKtQ5TsYM9PfpkIF9KWSClNAArgEUX7LMIeN/2ehUwS6i7c31Wm9HMox/tZ/WhUp6cP5xfzolFbFsGJ7exZtI9HKrL45fjf4mnk6fWoSrdLDUklVlDZ/F23SHK4+bCxmdxrzzA+/cmMznKn/9eeYgPdhdqHaZyGfYk+hCgqMP7YlvbRfexLSZeD/jZtkUIIQ4IIbYKIaZe7AOEEA8KITKFEJlVVVVXdAFK16puaueOf+xhfba1Jv9IWhQUbIGtr1Az6mZeq95NYkAiC6MWah2q0kN+PeHXSAkv+fshPYNh1b24mup5554kZo8I5LdfZvPS10fVPPa9WHffjC0DhkopxwJPAB8JIX7QDZRSvi2lTJJSJgUEqHnMtZJX0cgNb+7kaFkDf7tzHPemRsCZU/DZT5H+MTzv5UKToZnnJz2PTqj7+ANFsHswj417jC2lO/l88r3QWA6f/RQXneTvd43nnknD+Mf2kzy8fJ+a8bKXsudfawkQ1uF9qK3tovsIIRwAL6BGStkupawBkFLuA04AsVcbtNL11h0p56a/7qLdZOGTByexYNQQaD0DH94KZgP/mbyELSXbeXzc42oa4gHozhF3khKUwst5Kyia/Vs48R2s/S8cdIIXFo3k+evj+S6ngpv/tpvC6matw1UuYE+izwBihBARQggn4HZg9QX7rAbOrjZxC7BJSimFEAG2m7kIISKBGEAtX9OLGEwWXliTzcPL9xEZMIgvlqYyOswbTO2w4i6oLaBo0f/xSs77pASlcJdaInBA0gkdv5/yexyEA7+py8CU+gvY9y/Y8WcAlqRG8N6SCZTVt3Ld6zv4KkuNyOlNLpvobTX3R4H1QA7wqZQyWwjxohDibKH2XcBPCJGPtURzdgjmNCBLCHEQ603ah6WUtV19EUrnnKhq4ta/7+KfOwu5NzWclQ9PJsTbFSwW+PJROLUD48LX+c3Jz3DQOfD7Kb9XJZsBLGhQEE9PfJpDVYd4JyAIRt4M370AWSsBSIsL5OvHphI72J1HPzrA058fprldlXJ6A9Hb5q5ISkqSmZmZWofRr5ktkvd2nOTVDbk4O+j44y2JzB85xLrRYoE1P4cDy5EznuZ3Tm2sPL6SZdOXMT98vraBK73Cb7b/hq8LvuYv0/9E2uY/w+k9cMu7kHAjYB2a++r6XN7eXkCojyuv3JzI5Ch/jaPu/4QQ+6SUSRfbprpnA8yx8gZu/fsuXlqbw9SYAL59YnqHJG+GL5fCgeUw7desCAxj5fGV3D/yfpXklXOenfQsI/xG8OTOZ8i/9n8gLBlW3Q+HVwHgqNfxm2tG8MmDk9ALwR3/2MvTnx9WUx1rSPXoB4j6FiOvbczlgz2n8HR15PnrE1g0Jvj7ycjMRmu5JmsFpP2GvXEzeGjjQ0wNmcr/zfw/VbJRzlPeXM7irxfjonfh4znv4P3ZA3B6N9zwNxh9+7n9Wg1mlq3P5V+7TuLt5sSv5sVxW1IYep16zKar/ViPXiX6fq7NaObj9NO8vimfuhYDd6YM47/mxuLt5vT9Tq11sHIJFGyGGc9wPHER9667l0C3QD5Y8AHuTu6axa/0XoeqDnHfuvsY6T+Sv03/E24r74WT22DG0zDtV9Dhmcns0npeWH2U9MJaRoZ48ut5w5ka469mPe1CKtEPQEazhc/2FfOX7/IorW9jUqQfv70unvjgCx5jqC2Aj34CtSfh+v/lRMQk7lt/Hw46B/41/1+EeYRd/AMUBVhXuI4ntz3JhMETeH36n3Bd+2vrb4WjboOFr4Ojy7l9pZSsPlTKH9flUlLXSnKEL7+aF8eEcF8Nr6D/UIl+AGluN7Eio4h3txdQWt/GmDBvfjUvjtToi9wMO74BPn8IkPCT5RT4BHPfuvvQCR3vzXuPcK/wng5f6YPWnFjD0zueZuKQibw+8y8473oDNv0OQpKsN2l9ws/bv91k5pOMIl7flE9VYzvJ4b48ND2SGXGB6FRJp9NUoh8ATtU089He06zIKKK+1UhyhC+PTI8iLS7gh78eG9tg47OQ/hYEJsBPPiBXmHjk20cwSzP/nP9PIr0itbkQpU/6PO9znt31LJOGTOK1tNdwz99kvbEPcO1rkHjrD45pNVjLiu/uOElJXSsxge7cPWkYN4wNwcPFsYevoO9Tib6fajOa2XSskhUZRWw7XoVeJ5gbP5gHpkUybqjPxQ8qPQBfLIXKbEh5BGY/z5byPfx626/xcPLgrdlvEe0T3aPXofQPn+d9zgu7XyDCK4I3Z71JsNEI/3kAivZaSznzX4ZBfj84zmi2sPZwGf/YXsCRkgbcnPQsGhPCrUmhjA3zVnV8O6lE34+YzBbSC2tZfbCUrw+X0dhmIsjThcXJQ/nJhDCCvFwufmBrnfXX6Yx3YVAALHoTGTOH5TnLWZaxjBF+I3h95usEugX27AUp/cqu0l3895b/xlHvyF9m/oXRvgmwbRlsfxWcPWD28zD2btBdfBTXoaI6lu85xZqsUtqMFsL93LhhbAjXJQYTHagGBfwYlej7uFaDmd0F1WzIrmDD0Qpqmw24OemZnxDEjeNCmBzlf+nhaiYDHPgAtrwMLdUw4QGY+TRnsPDi7hf59vS3zB46m5emvISbo1vPXpjSLxXUFbD0u6WUt5Tz6JhHWZKwBH31cfj6v+DUTmvtfvbzEHHRyWwBaGwz8s2Rcv6zv5g9BdaH6WMC3VkwMojZ8YMZGeyl6vkXUIm+j7FYJLkVjew6UcO241XsLqjBYLLg7uzAzOGBzB8ZRFpcAG5ODpc+idkIWZ/A1leg7jSETYQFr0DwGLYXb+fZXc9S117Hz8f+nCUJS9Q4eaVL1bXV8eKeF9l4aiPjAsfxh6l/IGRQMBxaYZ02obEMIqbBjGdgaMqPnqusvpX1R8r55kg5GYW1WCT4uzsxLTaAaTEBTIryY7DnJX6THUBUou/ljGYL2aUN7Dt1hn2natlTUEttswGACP9BzIgLZMbwAJIjfHF20P/4yVpqrZNNZbwDDSUQPNb6jyl6FqXNZfx5359ZV7iOaO9oXp76sloKUOk2UkrWFKzhD3v/gEVauH/k/dyTcA8uUkLmP2HHa9BcBaHJMPERGLEQ9D/SeQFqmtrZnlfN5txKth6vos72tG1kwCBSInwZP8yXpGE+DPNzG3C1fZXoexGDycKJqiZyyhrIKq4nq7iO7NIG2k0WAEK8XUmO8GVylB+TovwI9bGjnGKxQOF2a28p+3MwtULEdJj4M4idR4OxkX8d+RfvZ7+PTuhYMnIJPx31U5z1zt18tYoCpU2lLMtYxrenvyVoUBCPj3ucBeEL0JvaYP+/Ye9bcOYkeIbA6MXWJ2v9Lz8VttkiySlrYNeJanadqGHfqTM0tlknUfMd5ERiqBeJod6MDPZkxBBPQn1c+3XyV4leA+0mM0W1LZyoaiavopG8yiaOVzSRX9mI0Wz9M3d11DMyxJNRId6MG+bN+GE+DPFyte8DzCbraIbctZD9BTQUg7OndWKplIdgcALlzeV8cPQDVh1fRYuphesir+PxcY8TNCioG69cUS4uozyDZRnLyKnN+f/bu9PYOM7zgOP/Z2b2vsjlJZEiJVIHZdZ1rdgQAudwFCeFqggWCgSwY6Sw0RRB0QRNgBRF03wp+qlo0aYFUqAw0gZtcxhuk6BCYcNOasEJEqe1I9mMJVmybJkiKVHLS9yLe87TDzOkKIWiKXmpNVfvDxjM7Ozs7jNc6XmPeWdf+hP9PD7yOEd2HSFsBeDsc14r9O1joC703QcjR2D4U9C5vlFgrquczeR45Z15Xhu/wujEAm9mcixNfJUIO+zdkmB3T4Ld3XF2dccZ6oqzNRluif5+k+g3gOsqM4Uyk/OLTF5ZZHJ+kQtzxeVlfK7IypnV+toi7OqOs3drgpGtXg1jqDOGY99E3/j8mHeL+fkX4dz/wOIc2EEYOgC/9QgMH6JqOfx08qccfesoL46/iKIcHDzIE7/xBHvTexv/hzCMm+CqywsXXuBbr3+L0ZlR2kJtHBo8xMO7HmYkPYLkL8Ov/gNGn4apUe9FnXtg58e9Vur2ByDStu7PK1ZqnL6U4/SlLKcuZTk7lePs5RzZ0tWfTw4HLHZ0xBhIR72lI0pfW4S+9gh9bZFNM6bfJPqbUK27sJI+dQAAClBJREFUzBcqzOQrzOTLzOTLTOfKZPzl8kKJS9lFLi+UqdTda17bHg0wkI7Sn44y1BljsCvGYKdXc4iH1u57/DWLV+DySW/c+8TLMPlLWPCn7o33wNDHYPgQ7HqIgmXx0sWX+MnETzg2fowr5St0hDs4PHSYx+56jN54b0P+NobRKKrK8cxxvnv6uxwbP0bVrbIztZMDAwf46LaPck/nPdjZi3DmWTj7LIy95HVJigXdI7Dtfm/0Tu+9XkHgrL8bUlWZzpU5l8lzfrbA+ekC52cKy5W0pW7UJfGQw9ZUmC2pMN2JMD3JEN2JEJ2JEJ1xb+mIBUlFAk1tGdxRiV5VKddcCuUa+XKNXGlpqZIt1cguVllYsVwpVpgveuu5QuWakn6laNCmJxmmKxGiNxVmSyrClmSIbe1Rr+Rvj5C82ZK/XoPcRa+mPnvOW2behMxpWLhw9bi2Ae8f9cAHYfBBZuKdvD77OicyJ3g18yqjM6PU3BqJQIIP932YwzsP80DvAzjWTRYuhtEEC+UFnh97nmfefoYTmRPUtU4ymOQDPR9gX/c+9nXvYzg5SHTqpNeinXgZJl+B0oL3BpYDncPQNez17XfshvQQtG+HaMc1P672bpYKgQm/lX7xyiKXFkpMLZS4lC0xnS2RyZWprTIRum0J7dEA7dEg7dEgbdEAqUhgeZ2KBEhGAiTDARJhh3jYIR5ySIQCxEL2zbXuV/GeE72IHAT+AbCBb6rqX133fAj4N+A+YBZ4RFXf8Z/7KvA5oA78sao+t9Zn3Wqin86V+eTXX6RQri33ga8lEXaWv4SlL6Y9GiAdC5GOB+mIBemMh+hKhOiMB9fffFOFcs7rVinOQnEeChnIZ7wRBtmL3tCy7KS37a4oWJwwpHdC916KXcNcbutjMpriQi3HeG6ct668xZn5M8yVvHHFjuUwkh7hvp77+Mi2j3Bv970ErM3RzDSM1WQrWX5+8ef8bPJnnMicYCw7BoAgDCQH2NO+hx3JHfTHt9GvQm9ulq75cQKZUzB9xhtKzIr//4EYpLZBcqt3sTfeA/Fu76bBWBdE0xBJe2snvK5CwXWV+eLVVv90rsxcwasozhYqfuWxwnzBr0wuVihV3Xd933DA4v7tab79B2sPN72R95To/TlfzwKfBCbw5pD9jKqeWnHMHwH3qOofisijwO+q6iMiMgJ8D9gP9AI/Bvaoav1Gn3eriX6xXOOvnxklHrRIhCxiDiRCFomgEA8KyaAQDwjxoBJzFAcX3KpXq3arUK94Y89rZaiXvRuNaiVvf3XR264Wve1KESp5qBS8dTmPW16gVspSq+So4lIVoSJCBaFsCSURSk6YxViaYjRNIZwgH06QDYTJOQHmLWGuXmKuNMf04jS5Su6a84s4EXYkdzCcHmZP+x7uSt/F3Z13E3bM+GGjdc0uzjI6Pcob82/w5vybnJ0/y0RugvqKFCIIHZEOOiOdpINttItDm6skamWS5SLxUo5IaYFoYZ5IcY6wWyesSkiVoCoBf+3YAZxQCieUQEJJ707eYByCMQhGvUIjELm6OGGvy8gOgRNcsQ6CFfDWtgOWQ9m1yFchX4VcBfJVJV9RshWXXFkpVJR81aUjEeVzD97akOe1Ev162vb7gXOq+rb/Zk8BR4BTK445AvyFv/2fwDfEG8d0BHhKVcvAeX9O2f3AS7dyImupFMY4PvPY8uOVxddqRZmKt19XHKPI1W0Bl6vHuAgqQl287boILlC3oB4FNwqQ9Je1uMAMlGegDI44xINx2kJtpMNphlJD7N+yn55YDz3RHnrjvWxPbqcj3NHSQ8MMYzUdkQ4ODBzgwMCB5X1Vt8pUfoqx3BhThSkyxQyXi5eZXZxlrjTHWGmObDlLrrqisuQAKRtSXev41Bo2c9g6h10BqwKWKrYqFooo3hpvij4BRGHpf6f4zy0t+M+vJDfYHs7E4cHGX6NcT6LvA8ZXPJ4Arm9bLB+jqjURWQA6/P2/uO61fdd/gIh8Hvg8wMDAwHpjv4YdSvg/xiV+80sQ8b8Gwd+2QPz9IiA2Ipb/2EIsGxEbsWwQG9sOIJaDWA6WFcCybATBtmwssbCwcCwH27KxxSZgBXAsB8dyCFpBgnYQx3KIOBFCdoiwEyYaiBJ1vCURTBBxWntsr2E0WsAK0J/spz+59lwJdbdOvpqnWC1SqBYo1ooUa0XKtTKleolyvUy1XqXiVqjUK9Tcmrdojbpbp6715bWrLnWto6q4bh11K6hbx3WruG4N3DqqddSto+r62y7ggrqoul637tI2Cqp4PSq6/LgvtnVD/mbvi6t1qvok8CR4XTe38h7xWDd/++n/bmhchmFsXrZlkwqlSIVSzQ6l6dZzmXcSWFl0bvP3rXqMiDhACu+i7HpeaxiGYWyg9ST6l4HdIjIoIkHgUeDodcccBR73tz8NvKBem+Qo8KiIhERkENgN/F9jQjcMwzDW4127bvw+9y8Cz+ENr/wXVT0pIn8JvKKqR4F/Bv7dv9g6h1cY4B/3NN6F2xrwhbVG3BiGYRiN13I3TBmGYdyJ1hpeaX6E3DAMo8WZRG8YhtHiTKI3DMNocSbRG4ZhtLj33cVYEZkGxpodxy3oBGaaHcRtZs75zmDOeXPYrqqr/sbD+y7Rb1Yi8sqNrni3KnPOdwZzzpuf6boxDMNocSbRG4ZhtDiT6BvnyWYH0ATmnO8M5pw3OdNHbxiG0eJMjd4wDKPFmURvGIbR4kyi3wAi8hURURHpbHYsG01E/kZE3hCRURH5oYi0NTumjSAiB0XkjIicE5E/a3Y8G01E+kXkmIicEpGTIvKlZsd0u4iILSInRKRlZjIyib7BRKQf+G3gQrNjuU1+BNytqvfgTSL/1SbH03AiYgP/CPwOMAJ8xp/4vpXVgK+o6gjwQeALd8A5L/kScLrZQTSSSfSN93XgT1l9TvKWo6rPq2rNf/gLvFnEWs1+4Jyqvq2qFeApvInvW5aqXlLV4/52Di/x/dp8z61GRLYBnwK+2exYGskk+gYSkSPApKq+1uxYmuT3gWebHcQG6APGVzxedZL7ViUiO4B9wP82N5Lb4u/xKmpuswNppPfF5OCbiYj8GNiyylNfA/4cr9umpax1zqr6X/4xX8Nr7n/ndsZmbCwRiQPfB76sqtlmx7ORROQwkFHVX4rIx5odTyOZRH+TVPUTq+0Xkd8EBoHXRAS8LozjIrJfVaduY4gNd6NzXiIiTwCHgYe0NW/MuCMnuReRAF6S/46q/qDZ8dwGHwIeFpFDQBhIisi3VfWzTY7rPTM3TG0QEXkHuF9VN9sv4N0UETkI/B3woKpONzuejSAiDt6F5ofwEvzLwGOqerKpgW0g8Wor/wrMqeqXmx3P7ebX6P9EVQ83O5ZGMH30xnv1DSAB/EhEXhWRf2p2QI3mX2z+IvAc3kXJp1s5yfs+BPwe8HH/e33Vr+kam5Cp0RuGYbQ4U6M3DMNocSbRG4ZhtDiT6A3DMFqcSfSGYRgtziR6wzCMFmcSvWEYRoszid4wDKPF/T8wtMqxS1cbAAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# 自由度を変えてt分布をプロット\n",
"x = np.linspace(-5, 5, 100)\n",
"\n",
"plt.plot(x, stats.t.pdf(x, 3), label=r\"$n=3$\")\n",
"plt.plot(x, stats.t.pdf(x, 10), label=r\"$n=10$\")\n",
"plt.plot(x, stats.t.pdf(x, 50), label=r\"$n=50$\")\n",
"plt.legend()\n",
"plt.plot()"
]
}
],
"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.6.9"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment