{ "cells": [ { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "btw4g5", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "Problem 1 Instructions\n", "----\n", "Compute the requested confidence interval for the population mean given the data below. Do not use the function you create in problem set 2" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "t4wq5q", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 1.1 \n", "Compute a 95% double-sided confidence interval (2.5% to 97.5%)" ] }, { "cell_type": "code", "execution_count": 39, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "bth4as5ga4gf", "locked": false, "points": 2, "schema_version": 1, "solution": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-0.41 +/- 1.74\n" ] } ], "source": [ "data = [-4.4,-10.1,1.6,4.0,3.4,-2.0,-5.8,-1.9,-3.1,1.8,-2.9,4.9,2.7,-9.5,2.8,4.3,-1.1,3.9,-9.1,5.6,1.1,-1.0,0.7,9.7,-1.6,-2.3,1.4,2.5,-9.0,1.1]\n", "\n", "### BEGIN SOLUTION\n", "import scipy.stats as stats\n", "x = np.mean(data)\n", "s = np.std(data, ddof=1)\n", "Z = stats.norm.ppf(0.975)\n", "print('{:.2f} +/- {:.2f}'.format(x, s * Z / np.sqrt(len(data))))\n", "### END SOLUTION" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "btrg5r4qgq5", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 1.2\n", "Compute an 80% double-sided confidence interval using the data above" ] }, { "cell_type": "code", "execution_count": 40, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "bvtrabtn", "locked": false, "points": 2, "schema_version": 1, "solution": true }, "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-0.41 +/- 1.14\n" ] } ], "source": [ "Z = stats.norm.ppf(0.90)\n", "print('{:.2f} +/- {:.2f}'.format(x, s * Z / np.sqrt(len(data))))" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "btranbtyabtrwath", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 1.3\n", "Compute a 95% double-sided confidence interval using the data below." ] }, { "cell_type": "code", "execution_count": 41, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "bgtrata", "locked": false, "points": 2, "schema_version": 1, "solution": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "12.88 +/- 3.57\n" ] } ], "source": [ "data = [11.5,10.1,20.5,12.4,7.9,9.1,18.1,13.4]\n", "\n", "### BEGIN SOLUTION\n", "x = np.mean(data)\n", "s = np.std(data, ddof=1)\n", "T = stats.t.ppf(0.975, df=len(data))\n", "print('{:.2f} +/- {:.2f}'.format(x, s * T / np.sqrt(len(data))))\n", "### END SOLUTION" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "btaeratt", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 1.4\n", "Using the same data above, compute a 95% double-sided confidence assuming the population standard deviation is 5.0" ] }, { "cell_type": "code", "execution_count": 42, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "btafraf", "locked": false, "points": 2, "schema_version": 1, "solution": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "12.88 +/- 3.46\n" ] } ], "source": [ "Z = stats.norm.ppf(0.975)\n", "print('{:.2f} +/- {:.2f}'.format(x, 5 * Z / np.sqrt(len(data))))" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "ntrwtrw", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "Problem 2 Instructions\n", "-----\n", "\n", "Plot the following quantities" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "bteragrae", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 2.1\n", "Plot a standard normal distribution. Typically, plotting from $\\mu-5\\sigma$ to $\\mu + 5\\sigma$ is a good choice when plotting normal distributions. Label your axes." ] }, { "cell_type": "code", "execution_count": 43, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "trsrfge", "locked": false, "points": 4, "schema_version": 1, "solution": true } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEKCAYAAADjDHn2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4XOWZ5/3vrd2LvGnxLi+yvOPYIAw2YHYwDbGTNAmG\nTkO6M0OYQJaLzqRJJ5N0kzfvdJJOh0xCT8LbIW8yCXFISGh3MBjMDrbxAsYrtiWv8qLNtixbttZ7\n/qgSFLLskmQdnVLp97muulxnq7qFKP3qnOc5z2PujoiIyPmkhF2AiIgkPoWFiIjEpbAQEZG4FBYi\nIhKXwkJEROJSWIiISFwKCxERiUthISIicSksREQkrrSwC+guubm5Pn78+LDLEBHpVTZs2FDl7nnx\n9kuasBg/fjzr168PuwwRkV7FzPZ1ZD9dhhIRkbgUFiIiEpfCQkRE4lJYiIhIXIGGhZktNLMdZlZi\nZg+dZ7/bzczNrDhm3deix+0ws5uDrFNERM4vsN5QZpYKPArcCJQB68xsmbtva7NfNvBF4K2YddOB\nJcAMYBSw0swmu3tzUPWKiMi5BXlmMRcocffd7t4ALAUWt7Pft4HvAWdi1i0Glrp7vbvvAUqiryci\nIiEI8j6L0cCBmOUy4LLYHcxsDjDW3f9sZl9pc+yaNseODqpQkSCcqm/i9V1VlFTUkpWeypyCIcwZ\nO5SUFAu7NJFOCzIs2vtEvD/ht5mlAD8EPtPZY2Ne417gXoCCgoIuFSnS3VpanMff3MOPXyqh5nTj\nh7YV5Q/knxbNYP6k3JCqE+maIMOiDBgbszwGOBSznA3MBF4xM4ARwDIzW9SBYwFw98eAxwCKi4vP\nChORnlbX0MR9v36b13ZWcu2UPO5dUMjssUM43djMKzsqeGTlLv7q52/x1Zunct/VE4n+vy+S8IIM\ni3VAkZlNAA4SabC+q3Wju9cA73+9MrNXgK+4+3ozOw08YWb/SqSBuwhYG2CtIhfsdEMz9zy+lg37\njvGdj8/krrkF74dBv4xUPnHxGP7iopH89z9s4rvPvUdTcwtfuL4o5KpFOiawsHD3JjN7AFgBpAKP\nu/tWM3sYWO/uy85z7FYzexLYBjQB96snlCQyd+erT21i/b5j/K8lc/joR0a1u19Weio/umM26SnG\nD17Yyeih/fjExWN6uFqRzjP35Lh6U1xc7BpIUMLyy1V7+dayrXx14RQ+f82kuPs3Nrfw1z9/i3cP\n1PCfX7iCSfnZPVClyNnMbIO7F8fbT3dwi1yg3ZUn+Z/PbueaKXn8t6sLO3RMemoK/2vJHPpnpPLg\nk+/S3JIcX9okeSksRC6Au/M//mMLGakpfPcvZ3WqwTp/UBb/tHgGm8pq+MWbewKsUuTCKSxELsDK\n7RW8WVLNgzdOZvigrE4ff+tFI7l2Sh6PrNzF0VMNAVQo0j0UFiJd1NDUwnee2cak/IH81eXjuvQa\nZsbXb51GXUMTP3mppJsrFOk+CguRLlq6bj97q+v4+q3TSE/t+kdpUn42n7xkLL9es4+yY3XdWKFI\n91FYiHRBY3MLP3t1N8XjhnLN5LjTF8f15RuLMIMfrdzVDdWJdD+FhUgXPP3OQQ4eP839107qlruw\nRw7ux51zC3h640GO1JyJf4BID1NYiHRSc4vzv18tZcaoQVwz5cLPKlp99soJNLc4v1ilnlGSeBQW\nIp300nsV7K48xX+7prBbx3YaO6w/t1w0kifW7Kf2TGP8A0R6kMJCpJN+tXovIwZlsXDGiG5/7c8t\nmEhtfRNPri/r9tcWuRAKC5FO2FN1itd3VXHXZQWkXUAPqHOZNWYIFxcM4Tdv7SNZhuKR5KCwEOmE\nX6/ZR3qqsWTu2Pg7d9GdcwvYXXmKtXuOBvYeIp2lsBDpoDONzfx+/QEWzhxJfnbn79buqNtmjSI7\nK40n1u4P7D1EOkthIdJBz28r58SZJu68NLizCojOfTFnNM9uPsIxDQEiCUJhIdJBT20oY/SQflw+\nMSfw91oyt4CG5hae3ngw8PcS6QiFhUgHlJ84w+u7Kvn4nNGkpAQ/Feq0kYOYPnIQT288azZhkVAE\nGhZmttDMdphZiZk91M72+8xss5ltNLM3zGx6dP14MzsdXb/RzH4aZJ0i8Tz9zkFaHD5x8egee8+P\nzxnNuweOs7vyZI+9p8i5BBYWZpYKPArcAkwH7mwNgxhPuPtF7j4b+B7wrzHbSt19dvRxX1B1inTE\nH98+yMUFQ5iYN7DH3nPR7FGYobMLSQhBnlnMBUrcfbe7NwBLgcWxO7j7iZjFAYA6lkvC2VVey47y\nWhbP7rmzCoDhg7K4ojCXp985qHsuJHRBhsVo4EDMcll03YeY2f1mVkrkzOKLMZsmmNk7ZvaqmV3V\n3huY2b1mtt7M1ldWVnZn7SLve2bzYczglpndf8d2PB+bM5r9R+vYeOB4j7+3SKwgw6K9VsCzvh65\n+6PuXgj8PfCN6OrDQIG7zwEeBJ4ws0HtHPuYuxe7e3FeXvcN6CYSa/nmw1w6fhj5XZgJ70LdOH04\n6anGc1uO9Ph7i8QKMizKgNgO6WOA8118XQp8DMDd6929Ovp8A1AKTA6oTpFz2lley87yk9w2a2Qo\n7z+4XzrzC3N5dssRXYqSUAUZFuuAIjObYGYZwBJgWewOZlYUs3grsCu6Pi/aQI6ZTQSKgN0B1irS\nrmc2RS5BLQzhElSrW2aOYP/ROrYeOhF/Z5GABBYW7t4EPACsALYDT7r7VjN72MwWRXd7wMy2mtlG\nIpeb7omuXwBsMrN3gT8A97m7BsqRHrd882Hmjh8W6PAe8dw0YwSpKboUJeFKC/LF3X05sLzNum/G\nPP/SOY57CngqyNpE4tlZXsuuipM8vHhGqHUMG5DBZROG8eyWw3zl5imh1iJ9l+7gFjmHF7aVA3Bz\nAPNWdNYtM0dQWnmKXeW1YZcifZTCQuQcVm4vZ9aYwQwPoRdUWzfPGIEZLN+sS1ESDoWFSDuqTtaz\n8cBxrp86POxSAMgflMXFBUN58b3ysEuRPkphIdKOl96rwB2un5Yfdinvu25qPpvKaqg4cSbsUqQP\nUliItOPF7eWMHJzFjFFn3Qsamtbgeum9ipArkb5IYSHSxpnGZl7fVcV1U/MxC3448o6aMjyb0UP6\n8aLCQkKgsBBpY83uauoamrlhWmK0V7QyM66fls8bu6o409gcdjnSxygsRNp4cXsF/dJTmVcY/Ix4\nnXXd1HxONzazend12KVIH6OwEInh7ry4vZwri3LJSk8Nu5yzXD4xh37pqby0XZeipGcpLERibD9c\ny6GaM9yQQL2gYmWlp3JlUS4vbi/XwILSoxQWIjFe3hH5xn7tlMQMC4hcijpUc4ZdFZpuVXqOwkIk\nxms7K5k2clAoc1d01ILJkblbXtupCb+k5ygsRKJO1jexYd8xFkzODbuU8xo9pB+FeQN4VWEhPUhh\nIRK1urSaphbn6qLEn3VxweQ81u45qi600mMUFiJRr+2spF96KpeMHxp2KXEtmJxHfVMLb+3RNC/S\nMxQWIlGv7apkXmEOmWmJ12W2rcsn5JCRlqJ2C+kxgYaFmS00sx1mVmJmD7Wz/T4z22xmG83sDTOb\nHrPta9HjdpjZzUHWKbKv+hT7qutYUJTY7RWt+mWkMnf8MIWF9JjAwiI6h/ajwC3AdODO2DCIesLd\nL3L32cD3gH+NHjudyJzdM4CFwL+1zsktEoTWP7qtPY16gwWTc9lVcZJDx0+HXYr0AUGeWcwFStx9\nt7s3AEuBxbE7uHvsDPQDgNa7jBYDS9293t33ACXR1xMJxKs7qxgztB8TcgeEXUqHtQbb67t0diHB\nCzIsRgMHYpbLous+xMzuN7NSImcWX+zksfea2XozW19ZqQ+MdE1DUwurS6tYMDkvoUaZjWfK8GyG\nD8rktZ1VYZcifUCQYdHep+6s8Qnc/VF3LwT+HvhGJ499zN2L3b04L6/3XD6QxPLO/mOcamhmQS/o\nMhvLzLiqKI83SqpobtHQHxKsIMOiDBgbszwGOHSe/ZcCH+visSJd9mZJFSkG8ycl3iiz8VxVlEvN\n6Ua2HKwJuxRJckGGxTqgyMwmmFkGkQbrZbE7mFlRzOKtwK7o82XAEjPLNLMJQBGwNsBapQ9bVVrN\nrDFDGJSVHnYpnTa/MNJ7681SXYqSYAUWFu7eBDwArAC2A0+6+1Yze9jMFkV3e8DMtprZRuBB4J7o\nsVuBJ4FtwHPA/e6uW1Wl252qb2LjgeMJOXdFR+RlZzJleDarSjS/hQQrLcgXd/flwPI2674Z8/xL\n5zn2O8B3gqtOBNbtPUpTizO/l4YFRC6fPfHWfs40NifkHBySHHQHt/Rpq0urSU81iscNC7uULptf\nmEt9Uwvv7D8edimSxBQW0qetKq1mTsFQ+mX03m/kl00cRorBKrVbSIAUFtJn1dQ1suVQTa++BAUw\nKCudWWOG8GaJwkKCo7CQPmvNnmrcP+hR1JtdMSmHd8tqqD3TGHYpkqQUFtJnrS6tJis9hdljh4Rd\nygW7ojCX5hZnrYYsl4AoLKTPWlVaxaXjh5GR1vs/BhePG0pGWgqrStWFVoLR+z8lIl1QWVvPzvKT\nSXEJCiArPZXicUPVbiGBUVhIn7R6d+QbeG9v3I51xaRc3jtSS9XJ+rBLkSSksJA+aXVpFdlZacwY\nNSjsUrpNa/Ct1qUoCYDCQvqkVaXVXDYhh7TU5PkIXDR6MNmZabrfQgKRPJ8UkQ4qO1bHvuq6pLoE\nBZCWmsJlE3N4U+NESQAUFtLntF6m6Y1DksdzxaQc9h+t48DRurBLkSSjsJA+Z3VpNTkDMpicnx12\nKd2utXeX2i2kuykspE9xd1bvrubywhxSUnrPFKodNXn4QHIHZqjdQrqdwkL6lL3VdRyuOZN07RWt\nzIx5hbmsKq3GXVOtSvcJNCzMbKGZ7TCzEjN7qJ3tD5rZNjPbZGYvmtm4mG3NZrYx+ljW9liRrmj9\nxp0sN+O1Z35hDhW19ZRWngq7FEkigYWFmaUCjwK3ANOBO81sepvd3gGK3X0W8AfgezHbTrv77Ohj\nESLdYFVpNSMHZzE+p3/YpQTmg/stdClKuk+QZxZzgRJ33+3uDcBSYHHsDu7+sru3dttYA4wJsB7p\n41panDWl1cwrzMEs+dorWhUM68/oIf00TpR0qyDDYjRwIGa5LLruXD4LPBuznGVm681sjZl9LIgC\npW/ZWVFL9amGpL4EBa3tFjms3l1NS4vaLaR7BBkW7X11a/f/XDP7NFAMfD9mdYG7FwN3AY+YWWE7\nx90bDZT1lZWV3VGzJLFV0ZvV5iVp43as+YU5HK9rZPuRE2GXIkkiyLAoA8bGLI8BDrXdycxuAL4O\nLHL390dAc/dD0X93A68Ac9oe6+6PuXuxuxfn5eV1b/WSdFaVVjMuJ3KJJtnN0zhR0s2CDIt1QJGZ\nTTCzDGAJ8KFeTWY2B/gZkaCoiFk/1Mwyo89zgSuAbQHWKkmuqbmFt3ZXJ22X2bZGDu7HxNwBareQ\nbhNYWLh7E/AAsALYDjzp7lvN7GEza+3d9H1gIPD7Nl1kpwHrzexd4GXgn91dYSFdtvXQCWrrm5iX\n5O0VseYV5vDW7moam1vCLkWSQFqQL+7uy4HlbdZ9M+b5Dec4bhVwUZC1Sd/S+g173sS+cWYBkXtJ\nfvPWfjYfrOHigqFhlyO9nO7glj5hVWkVk4cPJC87M+xSeszlE4cBareQ7qGwkKTX0NTCur1Hk77L\nbFs5AzOZOiJb40RJt1BYSNLbeOA4Zxpb+kSX2bbmF+ayfu8xzjQ2h12K9HIKC0l6q0qrMIPLJ/TF\nsMihvqmFd/YfD7sU6eUUFpL0VpVWM3PUYAb3Tw+7lB43d+IwUkzjRMmFU1hIUjvd0Mw7+4/1yUtQ\nAIOy0rlozBDdbyEXTGEhSW39vqM0NnufuRmvPfMLc9h44Din6pvCLkV6MYWFJLVVpdWkpRiXjh8W\ndimhmV+YQ1OLs27v0bBLkV5MYSFJbVVJFXMKhjAgM9D7TxNa8bhhpKea7reQC6KwkKRVc7qRzQdr\n+tQQH+3pl5HKnIKhareQC6KwkKS1ds9RWpw+3V7Ran5hDlsO1VBT1xh2KdJLKSwkaa0qrSIrPYU5\nBUPCLiV08wtzcYc1e3R2IV2jsJCktaqkmkvHDyMzLTXsUkI3e+wQstJT1G4hXdahVj8zG0NkPoqr\ngFHAaWAL8AzwrLtrDGRJKJW19ewor2XxnFFhl5IQMtJSuHT8MI0TJV0W98zCzH4BPA40AN8F7gQ+\nD6wEFgJvmNmCIIsU6aw1uyPfoPva4IHnM78wl53lJ6msrY+/s0gbHTmz+IG7b2ln/Rbgj9FZ8Aq6\ntyyRC7OqtJrszDRmjhoUdikJo7Whf/XuahZ9RGdc0jlxzyxag8LM8ttuM7Mp7t7g7iXtHWtmC81s\nh5mVmNlD7Wx/0My2mdkmM3vRzMbFbLvHzHZFH/d07seSvm51aRWXTcwhLVXNcq1mjBpEdlaaxomS\nLunMJ+l1M/tU64KZ/R3wp3PtbGapwKPALcB04E4zm95mt3eAYnefBfwB+F702GHAt4DLgLnAt8xM\nU31Jhxw8fpq91XXqMttGWmoKl03I0f0W0iWdCYtrgL82s9+b2WvAZCJ/yM9lLlDi7rvdvQFYCiyO\n3cHdX3b3uujiGmBM9PnNwAvuftTdjwEvEGkfEYlrVUnkm/P8SQqLtuYX5rCvuo6yY3XxdxaJ0eGw\ncPfDwHPAPGA88Ct3P3meQ0YDB2KWy6LrzuWzwLNdPFbkfatLq8kZkMHk/OywS0k4rQGqLrTSWR0O\nCzN7gchloZnAXwA/NLN/Od8h7azzc7z2p4Fi4PudOdbM7jWz9Wa2vrKy8nzlSx/h7qwqrWZeYQ4p\nKe39b9S3Tc7PJmdAhsJCOq0zl6Eedfe73f14tNF7PlBznv3LgLExy2OAQ213MrMbgK8Di9y9vjPH\nuvtj7l7s7sV5eXmd+FEkWe2pOsWRE2fUZfYcUlKMywsj7Rbu7X53E2lXR+6zMAB3fzp2vbs3ufu3\nY/dpYx1QZGYTot1rlwDL2rz2HOBnRIKiImbTCuAmMxsabdi+KbpO5LzeLG29v0LtFecyvzCHIyfO\nsKfqVNilSC/SkTOLl83sC2b2oXspzCzDzK4zs18CZ3Vtdfcm4AEif+S3A0+6+1Yze9jMFkV3+z4w\nEPi9mW00s2XRY48C3yYSOOuAh6PrRM5rdWkVowZnMS6nf9ilJKzWsy71ipLO6MhNeQuBvwV+a2YT\ngONAFpAKPA/80N03tneguy8HlrdZ982Y5zec603d/XEid46LdEhzS6S94oZpw2n/ZFcAxuf0Z+Tg\nLFaXVvPpy8fFP0CEDoSFu58B/g34NzNLB3KB0+5+POjiRDpjy8Eajtc1clWR2ivOx8yYV5jDKzsq\naWlxdQSQDulIm0WWmX3ZzH4C/A1QqaCQRPRG9P6KKyYpLOKZX5jL0VMN7CivDbsU6SU60mbxSyLd\nWjcT6TL7g0ArEumi13ZWMmPUIHIHZoZdSsKbF+0AoHYL6aiOhMV0d/+0u/8MuJ3IMOUiCeVUfRNv\n7z/GVUXqQt0Ro4f0Y3xO//fvdheJpyNh8f48jNEeTiIJ56091TQ2u9orOuHKolzW7K6moUnT0Uh8\nHQmLj5jZieijFpjV+tzMTgRdoEhHvLYzMoXqJeM03mRHLSjK41RDMxv2HQu7FOkFOjJEeaq7D4o+\nst09Lea5JguQhPD6rkrmTsghK11TqHbUvMIc0lKM13ZpqByJT4P9S6936PhpSitPsUCXoDolOyud\ni8cN5bWdCguJT2Ehvd4buyKNtGrc7ryrJ+ex9dAJTbUqcSkspNd7bVcl+dmZTB4+MOxSep0F0YB9\nXZeiJA6FhfRqLS3OmyVVXFmUqyE+umDGqEHkDMjQpSiJS2EhvdrWQyc4piE+uiwlxbiqKJfXd1XR\n0qIhy+XcFBbSq7X25NEQH123YHIe1aca2HZYPeHl3BQW0qu9sqOCmaMHkZ+dFXYpvVZrx4BXdSlK\nzkNhIb3W8boGNuw7xrVT8sMupVfLy85k+shBareQ81JYSK/12q4qWhyuUVhcsAWT89iw7xgn6zWi\nj7Qv0LAws4VmtsPMSszsoXa2LzCzt82sycxub7OtOTp73vsz6InEeuW9Cob2T2f22CFhl9LrLZic\nS1OLa2BBOafAwsLMUoFHgVuA6cCdZja9zW77gc8AT7TzEqfdfXb0said7dKHtbQ4r+ys5OrJeaRq\n8p4LVjxuGAMz03h5R0XYpUiC6si0ql01Fyhx990AZrYUWAxsa93B3fdGt2nYS+mUTQdrOHqqgWun\n6hJUd8hIS2HB5Fxe3F6Bu+ueFTlLkJehRgMHYpbLous6KsvM1pvZGjP7WPeWJr3dS+9VYPbBHchy\n4a6bOpyK2nq2HlIXWjlbkGHR3leTztz1U+DuxcBdwCNmVnjWG5jdGw2U9ZWV6snRl7yyo4I5Y4cw\ndEBG2KUkjWun5GEGK7eXh12KJKAgw6IMGBuzPAY41NGD3f1Q9N/dwCvAnHb2eczdi929OC9P3zD7\nisraejaV1XCdLkF1q5yBmcwZO4SX3lO7hZwtyLBYBxSZ2QQzywCWAB3q1WRmQ80sM/o8F7iCmLYO\n6dtabx5Tl9nud/204Wwqq6HixJmwS5EEE1hYRKdgfQBYAWwHnnT3rWb2sJktAjCzS82sDPgk8DMz\n2xo9fBqw3szeBV4G/tndFRYCwMpt5QwflMmMUZp7q7u1nq3p7ELaCrI3FO6+HFjeZt03Y56vI3J5\nqu1xq4CLgqxNeqczjc28urOSv7xktHrsBGDqiGxGD+nHi+9VsGRuQdjlSALRHdzSq7xZUsXpxmZu\nmj4i7FKSkplx3dR83thVxZnG5rDLkQSisJBe5fmt5WRnpnH5xJywS0la103L53RjM6t3V4ddiiQQ\nhYX0Gs0tzsrt5VwzNZ+MNP2vG5R5E3Pon5HK81vVhVY+oE+c9Bpv7z9G9akGbpo+POxSklpWeirX\nTs3nhW1HaNaESBKlsJBe4/mtR0hPNa6ZontqgnbLzBFUnWxg/d6jYZciCUJhIb2Cu/P8tnLmF+aS\nnZUedjlJ75opkUt9z209EnYpkiAUFtIr7Ko4yb7qOm6aoUtQPWFgZhoLivJYseUI7roUJQoL6SVW\nbIl8w71xmsKipyycOYJDNWfYVFYTdimSABQW0is8s/kwFxcMIX+Q5truKTdMyyctxXh2iy5FicJC\neoGSipO8d6SW22aNCruUPmVI/wzmFebw3JbDuhQlCgtJfMs3H8YM/uKikWGX0ucsnDmCvdV17Civ\nDbsUCZnCQhLenzcd4tJxwxgxWJegetpN00dgBss3HQ67FAmZwkIS2s7yWnaWn+TWWTqrCENedibz\nJuaw7N1DuhTVxyksJKE9sylyCeqWizRwYFgWzx7F3uo6Nh9Ur6i+TGEhCcvd+fOmQ1w2YRj52boE\nFZaFM0aSkZrCf2zs8ESXkoQUFpKwth46QWnlKfWCCtng/ulcPSWP/3z3kMaK6sMCDQszW2hmO8ys\nxMweamf7AjN728yazOz2NtvuMbNd0cc9QdYpiemPbx8kIzWF29ReEbrFs0dRUVvPWxq2vM8KLCzM\nLBV4FLgFmA7caWbT2+y2H/gM8ESbY4cB3wIuA+YC3zKzoUHVKomnqbmFZe8e5Lqp+QzpnxF2OX3e\n9VOHMyAjlWXv6lJUXxXkmcVcoMTdd7t7A7AUWBy7g7vvdfdNQEubY28GXnD3o+5+DHgBWBhgrZJg\nXt9VRdXJBj5x8eiwSxGgX0YqN88YwTObD2sGvT4qyLAYDRyIWS6Lruu2Y83sXjNbb2brKysru1yo\nJJ6n3i5jaP90rpmSH3YpEvWXl4yh9kwTKzQSbZ8UZFhYO+s62jrWoWPd/TF3L3b34rw8zXGQLE6c\naeSFbeV89COjNCNeApk3MYfRQ/rxhw1lYZciIQjyk1gGjI1ZHgN09ILnhRwrvdzyTYepb2rh43N0\nCSqRpKQYt18yhjdKqig7Vhd2OdLDggyLdUCRmU0wswxgCbCsg8euAG4ys6HRhu2bouukD/jtugMU\n5Q9k9tghYZcibdx+yRjc4akNB8MuRXpYYGHh7k3AA0T+yG8HnnT3rWb2sJktAjCzS82sDPgk8DMz\n2xo99ijwbSKBsw54OLpOktzWQzW8e+A4d84twKy9q5ESprHD+nPFpBx+v+EALbrnok9JC/LF3X05\nsLzNum/GPF9H5BJTe8c+DjweZH2SeJauPUBGWop6QSWwT14yli//biNrdlczf1Ju2OVID1HroSSM\nuoYmnn7nILdeNFL3ViSwhTNHMKR/Or9+a1/YpUgPUlhIwvjzpsPU1jdx12UFYZci55GVnsqnisey\nYms5R2rOhF2O9BCFhSSMJ97az6T8gRSP0836ie7Tl42jxZ0ndHbRZygsJCG8s/8YGw8c59OXqWG7\nNyjI6c81k/N4Yu0BGpraDsAgyUhhIQnh8Tf3kp2VxieLx8bfWRLC3fPGU3Wynud0R3efoLCQ0B06\nfprlmw+z5NKxDMgMtIOedKOrJ+dRMKw/v1y1N+xSpAcoLCR0v1q9D3fnnvnjwy5FOiElxfibK8az\nYd8xNuzTbVDJTmEhoapraOK3a/dzy8yRjBnaP+xypJPuuHQsQ/qn89NXd4ddigRMYSGh+t26A9Sc\nbuRvrxwfdinSBf0z0rh73nhe2FZOSUVt2OVIgBQWEpozjc389NVSLpswjEvGDQu7HOmie+aNIys9\nhcde09lFMlNYSGh+v6GM8hP1fPH6orBLkQuQMzCTTxWP5U/vHNRNeklMYSGhaGhq4aevlHJxwRDm\nF+aEXY5coP961UTc4dGXS8IuRQKisJBQ/OmdMg4eP80Xri/STXhJYOyw/txx6ViWrtvPgaOa6yIZ\nKSykx9U3NfPjl0q4aPRgrpmsGQ6TxQPXTcLM+PFLu8IuRQKgsJAe939W76Ps2Gm+unCKziqSyMjB\n/fj0ZeN46u2D7K48GXY50s0CDQszW2hmO8ysxMweamd7ppn9Lrr9LTMbH10/3sxOm9nG6OOnQdYp\nPaemrpHBsQ4AAAANX0lEQVQfv1TCgsl5XFWks4pk8/lrC8lITeEHL+wMuxTpZoGFhZmlAo8CtwDT\ngTvNbHqb3T4LHHP3ScAPge/GbCt199nRx31B1Sk9699eKeHEmUYeWjg17FIkALkDM/nc1RN5ZtNh\n1u7RXd3JJMgzi7lAibvvdvcGYCmwuM0+i4FfRp//AbjedF0iaR04WscvVu3lE3PGMH3UoLDLkYB8\nbkEhowZn8Y/LttKsqVeTRpBhMRo4ELNcFl3X7j7RObtrgNZ+lBPM7B0ze9XMrgqwTukB7s63lm0l\nLcX4u5smh12OBKhfRir/cOs0th0+we/WHYh/gPQKQYZFe2cIbb9mnGufw0CBu88BHgSeMLOzvoqa\n2b1mtt7M1ldWVl5wwRKcFVuP8NJ7FTx442RGDekXdjkSsFsvGsncCcP4l+d3cOxUQ9jlSDcIMizK\ngNjJCcYAh861j5mlAYOBo+5e7+7VAO6+ASgFzvo66u6PuXuxuxfn5amxNFGdrG/iH5dtY9rIQXxG\nI8v2CWbGPy2awYnTjXz7z9vCLke6QZBhsQ4oMrMJZpYBLAGWtdlnGXBP9PntwEvu7maWF20gx8wm\nAkWABp7ppf5lxQ7Ka8/wnY/PJC1VvbX7imkjB/H5awr54zsHeem98rDLkQsU2Cc32gbxALAC2A48\n6e5bzexhM1sU3e3nQI6ZlRC53NTavXYBsMnM3iXS8H2fu6trRS/0+q5K/v9Ve7n78nFcXKC5tfua\nB64rYsrwbL72x83UnG4Muxy5AOaeHL0ViouLff369WGXITGOnWpg4Y9eIzsrnf984Er6ZaSGXZKE\nYFPZcT726Jss+sgofnjHbN2ImWDMbIO7F8fbT9cEJBDuzj/8aTNHTzXwyB2zFRR92KwxQ/jyDZN5\neuMhnlyv3lG9lcJCAvHzN/bw7JYjPHjjFGaOHhx2ORKy+6+dxJWTcvnmf2xl++ETYZcjXaCwkG73\nxq4q/t/l27l5xnA+t2Bi2OVIAkhNMX54x2wG9Uvn/t+8TU2d2i96G4WFdKv91XU88Nu3mZQ/kB98\najYpKbo+LRF52Zk8etfFHDhWx+d+vZ76puawS5JOUFhIt6msrefux9/CHR7762IGZqaFXZIkmLkT\nhvH92z/Cmt1HeeipzSRLB5u+QJ9m6RY1pxu5+/G1lJ+o59f/ZS7jcweEXZIkqI/NGc2Bo3X84IWd\n5AzI4Ou3TlMPqV5AYSEXrOZ0I3/zi7WUVNTy7/dcyiXjhoVdkiS4B66bRNXJev79jT048A0FRsJT\nWMgFqag9wz2Pr6OkopYf3zmHqzXznXSAmfGPi2ZgZvz8jT00tzjfvG262rgSmMJCumxP1Sk+84u1\nVJyo5+f3XMoCBYV0gpnxrY9OJzUlEhiHjp/mkSWz6Z+hP0uJSA3c0iUrt5Wz6CdvcOJ0I7/5r5cp\nKKRLzIxv3DqNb310Oiu3l/Opn63m4PHTYZcl7VBYSKc0NLXw3efe47/8aj3jcwbwn1+4UmM+yQUx\nM/7mign8f3cXs7eqjlseeY3lmw+HXZa0obCQDttUdpyP/vgN/vcrpSy5dCy/v28eY4b2D7ssSRLX\nTxvOM1+8kgl5A/n8b97m7558l+qT9WGXJVG6OChxVZ+s55GVu3hi7X5yB2bw+GeKuW7q8LDLkiQ0\nLmcAf7hvHj9auYufvlrKyu3lfHXhFO4oHqvh7UOmUWflnI7XNfCr1ft47LXdnG5s5q65BXzl5ikM\n7pcedmnSB+wqr+UbT2/hrT1HmZg7gC/dUMRts0aRqh5T3aqjo84qLOQsO47UsnTdfn637gB1Dc3c\nOH04f79wKpPyB4ZdmvQx7s6KrUd4ZOUu3jtSS8Gw/tw9bxyfvGQsg/vrS0t3UFhIpxw6fpoXtpXz\n1NtlbCqrIS3FWPSRUdx79USmjjhr+nORHtXS4jy39Qi/eHMP6/YeIyMthWun5HHbrFFcPy1f3W0v\nQEKEhZktBH4EpAL/7u7/3GZ7JvAr4BKgGrjD3fdGt30N+CzQDHzR3Vec770UFp1z7FQDGw8cZ+3e\no7z8XgXvHakFIlNh3n7JGBbPHkXuwMyQqxQ527ZDJ3hy/QGWbz5MRW09mWkpzJ0wjPmFucwvzGHa\nyEFkpKl9o6NCD4voHNo7gRuBMiJzct/p7tti9vk8MMvd7zOzJcDH3f0OM5sO/BaYC4wCVgKT3f2c\nw1QqLNp3uqGZPVWnKK08ye7KyL+bD9awp+oUEBk6unjcUK6bms91U/MpGp4dcsUiHdPc4qzdc5Tn\ntx1hVUk1O8ojX3jSU40pI7K5aPRgpgzPZlzuAMYN68+Yof0VIu3oaFgEee42Fyhx993RgpYCi4Ft\nMfssBv4x+vwPwE8sMkDMYmCpu9cDe6JzdM8FVgdYb0JpaXEamlsij6YWGqP/NjS1UN/Uwsn6JmrP\nNFF7pvH95ydON1J5sp6KE/VU1J6h/ET9h+Y9NoNRg/sxY9QgPlU8ljkFQ5g1ZrBO4aVXSk0x5hXm\nMK8wB4gMPbN2z1G2HDzBloM1LN98hN+u/WBmvhSDEYOyyMvOJHdg5JEzMINhAzIYmJnGgMy09/8d\nkJnKgIw00tNSSE81MlJTSH//YX1yHKsg/0qMBmLnUCwDLjvXPu7eZGY1QE50/Zo2x44OosjjdQ3c\n/tPVuDsO4OCReogu4g6OR/6NORFrPSZ2e+wxtB7T3uvFvN8H650Wh8bmFppaOn/Gl5GaQl52JnnZ\nmUzIHcDlE3MYPiiLcTn9KcwbyITcAWSla3pTSU752VncNmsUt80aBUQ+T5Un69lfXce+6jr2VZ+i\n7NhpKk/Wc7jmDJsP1lB9qoHmLnzW0lON9NQUUs0wi9xYmGKQErNsfLD8wfrocz7Y5yx23sXIujZh\nNW3kIH5855xO/xydEWRYtPcztv2tnGufjhyLmd0L3AtQUFDQ2fqAyLeTKa2XXowP/RJbfx8f+sUa\nGNZm2wfrPvgdRpfbbI99PWK2tT5PMUhPTSEjLfItJjPtg+cZMeuzs9Kij/T3n2emKQhEWpkZ+dlZ\n5GdnUTy+/ZGQW1qc2vomTkUfJ+ubOFXfzMn6JuoamiJn9M1OY/TsvqnFP3Sm3+z+oS9/LdEvfO+v\n88i61m2t61u8nT9ocNb8Hu3GWDsrxw7t18n/Op0XZFiUAWNjlscAh86xT5mZpQGDgaMdPBZ3fwx4\nDCJtFl0pMjsrnUf/6uKuHCoivVxKijG4X7ruHeqAIFt71gFFZjbBzDKAJcCyNvssA+6JPr8deMkj\n0boMWGJmmWY2ASgC1gZYq4iInEdgZxbRNogHgBVEus4+7u5bzexhYL27LwN+DvyfaAP2USKBQnS/\nJ4k0hjcB95+vJ5SIiARLN+WJiPRhHe06q07HIiISl8JCRETiUliIiEhcCgsREYlLYSEiInElTW8o\nM6sE9oVdRxfkAlVhF9HD9DP3DfqZe4dx7p4Xb6ekCYveyszWd6TbWjLRz9w36GdOLroMJSIicSks\nREQkLoVF+B4Lu4AQ6GfuG/QzJxG1WYiISFw6sxARkbgUFgnEzL5iZm5muWHXEjQz+76ZvWdmm8zs\nT2Y2JOyagmBmC81sh5mVmNlDYdcTNDMba2Yvm9l2M9tqZl8Ku6aeYmapZvaOmf057FqCoLBIEGY2\nFrgR2B92LT3kBWCmu88CdgJfC7mebmdmqcCjwC3AdOBOM5seblWBawL+zt2nAZcD9/eBn7nVl4Dt\nYRcRFIVF4vgh8FXOMZNisnH35929Kbq4hshsiMlmLlDi7rvdvQFYCiwOuaZAufthd387+ryWyB/P\n0eFWFTwzGwPcCvx72LUERWGRAMxsEXDQ3d8Nu5aQ/C3wbNhFBGA0cCBmuYw+8IezlZmNB+YAb4Vb\nSY94hMiXvZawCwlKkHNwSwwzWwmMaGfT14F/AG7q2YqCd76f2d3/I7rP14lcuvhNT9bWQ6yddX3i\nzNHMBgJPAV929xNh1xMkM7sNqHD3DWZ2Tdj1BEVh0UPc/Yb21pvZRcAE4F0zg8jlmLfNbK67H+nB\nErvduX7mVmZ2D3AbcL0nZx/uMmBszPIY4FBItfQYM0snEhS/cfc/hl1PD7gCWGRmfwFkAYPM7Nfu\n/umQ6+pWus8iwZjZXqDY3XvbYGSdYmYLgX8Frnb3yrDrCYKZpRFpvL8eOAisA+5y962hFhYgi3zj\n+SVw1N2/HHY9PS16ZvEVd78t7Fq6m9osJCw/AbKBF8xso5n9NOyCulu0Af8BYAWRht4nkzkooq4A\n/hq4Lvp73Rj9xi29nM4sREQkLp1ZiIhIXAoLERGJS2EhIiJxKSxERCQuhYWIiMSlsBARkbgUFiIi\nEpfCQiQgZnZpdL6OLDMbEJ3fYWbYdYl0hW7KEwmQmf0/RMYL6geUufv/DLkkkS5RWIgEyMwyiIwJ\ndQaY7+7NIZck0iW6DCUSrGHAQCLjYGWFXItIl+nMQiRAZraMyAx5E4CR7v5AyCWJdInmsxAJiJnd\nDTS5+xPR+bhXmdl17v5S2LWJdJbOLEREJC61WYiISFwKCxERiUthISIicSksREQkLoWFiIjEpbAQ\nEZG4FBYiIhKXwkJEROL6v+qqQTGUeDjBAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "#make some ponits\n", "x = np.linspace(-5, 5, 500)\n", "#plot them using the scipystats norm pdf\n", "plt.plot(x, stats.norm.pdf(x))\n", "plt.xlabel('x')\n", "plt.ylabel('P(x)')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "getrarefa", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 2.2\n", "Plot the uncertainty in the population mean, $\\mu - \\bar{x}$ vs $P(\\mu - \\bar{x})$, given that we know the population standard deviation is 0.25 and we have 7 samples" ] }, { "cell_type": "code", "execution_count": 44, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "bterarfare", "locked": false, "points": 8, "schema_version": 1, "solution": true } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEMCAYAAADNtWEcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4XFd5P/DvOzNarN0arZYly5Jsy7udON7tJM5mk5CQ\nhEBCgPBr2lAolBZoIQW6UAoEfiwNpZQ0KaQQAiEkbUhMguN4jRTbsp3IsmVtXiXLkkb7Yi0z8/YP\njRzZkaWRPHPvzNzv53n8aEZzNeedK3m+c8+59xxRVRARkfXYzC6AiIjMwQAgIrIoBgARkUUxAIiI\nLIoBQERkUQwAIiKLYgAQEVkUA4CIyKIYAEREFuUwu4DxpKWlaX5+vtllEBGFjYMHD7pUNd2fbUM6\nAPLz81FWVmZ2GUREYUNETvu7LbuAiIgsigFARGRRDAAiIotiABARWRQDgIjIohgAREQWxQAgIrKo\nkL4OgCiUuD1elNS14ui5LogA8zITsa4oDdEOfo6i8MQAIPLDqxWN+MYrlahvv3DJ99MSYvClzfPw\nwWtnQkRMqo5oahgAROPwehVff/kYfl5yCgtnJOGrt8/HmsI02G2C/Sdb8W9v1OJvni/HvpNt+Obd\ni3k0QGGFAUB0BaqKv3vxCH594Cz+ZN1sPPq+YkTZ332D31SciRvmZuDxN2rww9drcGHQg8cfWA67\njUcCFB4YAERX8B+7TuDXB87iMzcW4Qu3zh2zi8dmE/zVzXMRH+3Av2ytxIyUWHzl9gUmVEs0eQwA\nojG8daIV33ntON6/dMYV3/xH+7ONBTjb3of/3HMS1+Wn4taFWQZVSjR17LAkukzPgBtf/O07yEuN\nw2P3LvZ7cPcrt8/H4pxkfPmFI2jrHQxylURXjwFAdJnHt9egoeMCvv+hpYiL9v8gOcZhx/+/bym6\nLgzhG68cC2KFRIFheACIiF1EDovIy0a3TTSRk65e/OzNk7jv2pm4dlbqpH9+XlYiPnl9AV441IDD\nZ9qDUCFR4JhxBPA5AJUmtEs0oW9urUS03YYv3jZvys/x6RuKkJYQjcdePQ5VDWB1RIFlaACIyEwA\ntwN40sh2ifxRdqoN24414dM3FiEjMXbKzxMf48BnN83BWyfasKu6JYAVEgWW0UcAPwTwtwC8BrdL\nNKEf76hFanw0/mTd7Kt+rgdW5iEvNQ7ffa2KRwEUsgwLABG5A0Czqh6cYLtHRKRMRMpaWvjpiYxR\n0dCJHVUteHj9bEyLtl/180U7bPjMjUU4eq4Lb9a2BqBCosAz8ghgHYA7ReQUgF8D2CQiv7x8I1V9\nQlVXqOqK9HS/FrYnumo/2VmHxFgHPrZmVsCe867lM5CeGIOf7q4L2HMSBZJhAaCqj6rqTFXNB3A/\ngDdU9aNGtU90JWfb+rC1ohEfWz0LSbFRAXveGIcdn1ibjz01LlQ2dgXseYkChdcBkOX9ct9p2EQC\n+ul/xEdXzUJctB1P7jkZ8OcmulqmBICq7lTVO8xom2i0/iEPnjtwFrfMz0R28rSAP39yXBTuuSYH\nL5efQ2ffUMCfn+hq8AiALO2V8ka09w3h40H49D/igZV5GHB78cLh+qC1QTQVDACytF/uO43C9His\nKXQGrY2FM5KxNDcFv9p3hqeEUkhhAJBl1Tb34PCZDjywMi/oq3l9ZGUuapp7cPA0p4eg0MEAIMv6\n3aF62G2Cu5blBL2tO5bMQEKMA785cDbobRH5iwFAluTxKl481IDr56YjPTEm6O3FxziweVEWXq04\nj/4hT9DbI/IHA4AsqaTOhfNd/bj3mpmGtXn38hx0D7ixvbLZsDaJxsMAIEt64VADkmIduGl+hmFt\nri5wIjMpBi8ebjCsTaLxMADIcvqHPNh2rAlbFmUjNurq5/3x18h4w86qZq4YRiGBAUCWs7u6BT0D\nbty+JNvwtj+wLAdur+IPFY2Gt010OQYAWc7WI41IiYsK6rn/VzI/OxGz0+LxasV5w9smuhwDgCxl\npPtn88IsRNmN//MXEWxelIWSula0sxuITMYAIEvZVd2C3kEP3rfY+O6fEVsWZcHjVWyrbDKtBiKA\nAUAWs/VII6ab1P0zYnFOMnJSprEbiEzHACDL6B/y4PVjTbjNpO6fESKCLYuysLfGhe5+zhBK5mEA\nkGWU1rWid9CDzYuyzC4FWxZnYdDjxRvHeVEYmYcBQJbxemUT4qLtWF1gXvfPiOW505GZFIOtR3g6\nKJmHAUCWoKp443gzNsxJM/Tiryux2QS3LsjCnhoX5wYi0zAAyBKOnutCY2c/bpqfaXYpF22an4G+\nQQ/2nWwzuxSyKAYAWcL2ymaIAJuKjZv7ZyJrCpyYFmXHdp4OSiZhAJAlbD/ehGW5KUhLCP7Uz/6K\njbJjXVEatlc2c6UwMgUDgCJeU1c/yus7cXMIdf+MuHl+Bho6LqC6qcfsUsiCGAAU8UZOtTRy6md/\n3ejrknqd3UBkAgYARbztlU3ISZmGeZmJZpfyHplJsVick8zrAcgUDACKaP1DHuytdeHm+RlBX/h9\nqjYVZ+DQmXa09gyYXQpZDAOAIlrpiVb0D3mxKQT7/0dsKs6AKrCnxmV2KWQxDACKaLurWxDjsGHV\n7FSzS7miRTnJmB4Xhd3VLWaXQhbDAKCItru6BStnp4bE1b9XYrcJ1s9Jx+4aF7xeng5KxmEAUMRq\n6LiAupZeXD833exSJrRxThpcPQOoPN9ldilkIQwAilgjXSobwyEAfDXuruY4ABmHAUARa3d1C7KS\nYjEnI8HsUiaUmRSL4qxEjgOQoRgAFJHcHi/21rqwcW5ayJ7+ebmNc9NRdroNvQNus0shi2AAUER6\np74D3f3usOj+GbFxTjqGPIq3TrSaXQpZBAOAItKuahdsAqwvSjO7FL+tyJ+O2Cgbu4HIMAwAiki7\nq1uwZGYKUuKizS7Fb7FRw6uV7eYFYWQQBgBFnI6+QZTXd4RV98+IjXPScdLVi7NtfWaXQhbAAKCI\n89aJVnh1+Nz6cLNx7nDNe2t5FEDBxwCgiFNS14q4aDuW5qaYXcqkFaYnICMxBm8yAMgADACKOCV1\nrbguPxVR9vD78xYRrCtKQ2ldK6eFoKAz7H+IiMSKyH4ReUdEjorIPxnVNllHc1c/apt7sLbQaXYp\nU7a20InW3kFUNXWbXQpFOCM/Ig0A2KSqSwEsA7BZRFYb2D5ZQKnvHPq1heHX/z9ire/U1ZI6Xg9A\nwWVYAOiwkYVPo3z/eIxLAVVS24qkWAcWzEgyu5Qpy0mZhnxnHEo4DkBBZmgnqYjYReRtAM0Atqnq\nvjG2eUREykSkrKWFF8TQ5JSccGF1gRN2W3hM/3Ala4vSsO9kG9wer9mlUAQzNABU1aOqywDMBLBS\nRBaNsc0TqrpCVVekp4ffedxknrNtfTjbdiGs+/9HrCtMQ8+AG+/Ud5pdCkUwU06TUNUOADsBbDaj\nfYpMpb4+87VhNP3DlazxhRi7gSiYjDwLKF1EUny3pwG4GcBxo9qnyFdS50JaQnRYTP88kdT4aCzI\nTuJAMAWVkUcA2QB2iEg5gAMYHgN42cD2KYKpKkrqWrGmMHymf57I2kInDp5pR/+Qx+xSKEIZeRZQ\nuaouV9UlqrpIVb9uVNsU+epaetHcPRAR/f8j1hWlYdDtRdmpdrNLoQgVfpdKEo2htG64rzySAmDl\n7FQ4bII36zgOQMHBAKCIUFLXipyUachLjTO7lICJj3FgWW4KB4IpaBgAFPa8XkXpiVasKXRGTP//\niLVFaTjS0InOC0Nml0IRiAFAYa/yfBc6+oYiqvtnxLpCJ7wK7OMykRQEDAAKeyPn/6+JwABYlpeC\naVF2ng5KQcEAoLBXWteKgrR4ZCdPM7uUgItx2LEifzpKOBBMQcAAoLDm9nix72RbRH76H7G2MA3V\nTT1o6R4wuxSKMAwACmtHGjrRM+AO6+mfJzIytlHKcQAKMAYAhbWRvvHVBakmVxI8C2ckITHWcfFa\nB6JAYQBQWCuta0VxViKcCTFmlxI0DrsNq2Y7ORBMAccAoLA14PbgwKm2iO7+GbG20InTrX2ob+8z\nuxSKIAwACluHz3RgwO2NyPP/L7e2yDcOwKMACiAGAIWtkrpW2ARYGcH9/yPmZiTCGR/NAKCAYgBQ\n2Cqtc2FxTjKSYqPMLiXobDbB6sLhcQBVLqVNgcEAoLDUN+jG4TMdWGOB/v8RawudON/Vj5OuXrNL\noQjBAKCwdOBUO9xetUT//4iRwW6eDUSBwgCgsFRS50KUXbAif7rZpRgm3xmH7ORYjgNQwDAAKCyV\n1rViee50xEU7zC7FMCKCNYVOlJ5ohdfLcQC6epMKABGJFxF7sIoh8kdn3xAqGjojev6fK1lbmIa2\n3kFUNXWbXQpFgHEDQERsIvIREXlFRJoBHAfQKCJHReS7IjLHmDKJ3rXvZCu8GlnLP/prJPQ4DkCB\nMNERwA4AhQAeBZClqrmqmgFgA4C3AHxbRD4a5BqJLlFS14rYKBuW5aWYXYrhclKmId8Zx3mBKCAm\n6kC9WVXfsxadqrYB+B2A34lI5J+ETSGltK4VK2alIsZhzd7INYVpePmdc3B7vHDYOYxHUzfuX8/I\nm7+IHBGRZ0TkSyKyRURmishXRm9DZISW7gFUNXVfnBrBitYWOtE94EbFuS6zS6Ew5+/Hh+sB/CeA\nCwDuB1AB4PZgFUV0JW/55sS3wgRwV7K6YGQcgN1AdHX8CgBVbVPVnar6uKo+BOA6ADXBLY3ovUrq\nXEiMdWDRjCSzSzFNemIM5mUm8noAump+BcDlZ/uoag2AJUGpiGgcJXWtWDXbafm+7zWFThw41YYB\nt8fsUiiM+fu/6AkROSMipSLyUxF5GkCFiMQFszii0erb+3C6tc+Sp39ebm2hE/1DXrx9psPsUiiM\n+dsFdKOq5gH4MIBXANQCiAPwjogcD2J9RBeNdHlYeQB4xKoCJ2zC6wHo6kzqOnpVPQPgDICXRr4n\nIgmBLopoLCV1rXDGR2NeZqLZpZgueVoUFuUko7SuFX99i9nVULi66o5UVe0JRCFE41FVlNS5sKbQ\nCRExu5yQsKbQicNn29E36Da7FApT1h5Jo7BxwtWLpq4BS5/+ebm1hWkY8ijKTrWbXQqFqclOBvf+\nYBVCNJ6Rvm4OAL/ruvzpcNiE4wA0ZZM9AviXoFRBNIHSOhdyUqZhlpMnno2Ii3ZgeV4K5wWiKZts\nALDzlQzn9SpK61rZ/z+GNYVpONLQic4LnJGFJm+yAcBVKMhwlee70N43xO6fMawtdMKrwP6TbWaX\nQmGIg8AU8i6e/88B4PdYnpeCGIeN8wLRlDAAKOSV1LWiID0eWcmxZpcScmIcdlyXn8p5gWhKJhsA\nTUGpgugKhjxe7DvRyu6fcawpdOL4+W64egbMLoXCzKQCQFWnfM2hiOSKyA4RqfQtKfm5qT4XWUd5\nfSd6Bz3s/hnHSDiOTJVN5C8ju4DcAL6gqvMBrAbwFyKywMD2KQyNnOI4Mgc+vdfinGQkxDh4PQBN\nmmEBoKqNqnrId7sbQCWAHKPap/C0t9aFBdlJSI2PNruUkOWw27BqNscBaPKuGAAisn/U7R9KAE/A\nFpF8AMsB7BvjsUdEpExEylpaWgLVJIWh3gE3Dp5ux4a57P6ZyJpCJ066enGu44LZpVAYGe8IYPRi\n7z0AXhKReAAQkVtF5M2pNOibPfR3AP5KVd+zqKmqPqGqK1R1RXp6+lSaoAix/2QbhjyKDUX8O5jI\nyBgJjwJoMsYLgO6RG6r6VQDPAtgpInsBfAHAlyfbmIhEYfjN/xlVfWGyP0/WsrumBTEOG1bkTze7\nlJBXnJWI6XFRHAegSbniegCqunHktojcBODPAPQCyAbwsKpWTaYhXxfSUwAqVfX7UyuXrGRPjQur\nCpyIjbKbXUrIs9kEawqdKKlzQVU5ZQb5xd9B4K8A+Jqq3gDggwB+IyKbJtnWOgAfA7BJRN72/Xvf\nJJ+DLKKx8wJqm3uwoYj9//5aX5SOxs5+1LVwiQ7yz7grgomI6LCLb/aqekREtgB4HsAbI9tM1JCq\n7gUnkyM/7akZPv2TA8D+2+jbV7uqXSjK4KppNLGJjgB2iMhnRSTvsu+3Avi2b3H4h4JTGlnZnhoX\n0hNjuPzjJMycHoeC9HjsrubZc+SfiQJgMwAPgGdFpFFEjonICQA1AO4E8ANV/XmQaySL8XoVb9a6\nsKEojX3Zk7RxTjr2nWxF/5DH7FIoDIwbAKrar6r/rqrrAOQBuAnANao6S1X/TFXfNqRKspRjjV1o\n6x1k988UXD83Hf1DXk4PTX4ZNwBE5CERcYlIG4AnAfSoaocxpZFV7a4Z7sJYxwHgSVtVkIpou43d\nQOSXibqAvgbgFgDFAM4A+GbQKyLL21vjQnFWIjISOf3zZMVFO3Dd7OkXQ5RoPBMFQJeqHlbVZlX9\nGoCVRhRF1tU36EbZqXZsmMNP/1N1/dx0VDf1oLGT00LQ+CYKgGzf3DwbRCQdl04PQRRw+060YdDj\nxfo5nP5hqjbOHd53e6q5ShiNb6IA+AcASwB8A0AVgEUislVEviUiDwS9OrKcnVXNiI0ant2SpmZe\nZiIyk2Kwi91ANIFxLwRT1SdG3xeRmRgOhMUA3ofh+YGIAkJV8UZVM9YVpnH6h6sgItgwJx3bjjXB\n41XYbTyVlsY22RXB6lV1q6o+pqofC1ZRZE11Lb0423YBNxRnmF1K2Ns4Nx2dF4ZQXs+T9ujKuCg8\nhYydVc0AgBvnsf//ag1fRAfs4umgNA4GAIWMHVXNmJuZgJnT48wuJexNj4/GstwU7DjebHYpFMIY\nABQSegbc2H+yDTfOY/dPoNxUnIF36jvR3N1vdikUohgAFBL21rgw5FHcwAAImE3FmQCAncfZDURj\nYwBQSNhZ1YzEGAdX/wqg+dmJmJEci9crm8wuhUIUA4BMp6rYUdWMDXPTEGXnn2SgiAg2zc/A3loX\nZwelMfF/G5musrEbTV0D7P4JgpuKM9E36ME+zg5KY2AAkOler2yCCHADT/8MuDWFTsRG2bCd3UA0\nBgYAme6Px85jeW4KZ/8MgtgoO9YXpWN7ZTP8WLmVLIYBQKY613EBFQ1duHVhltmlRKyb5megoeMC\nqpu4WDxdigFAptp2bLhr4tYFmSZXErlGrq3g2UB0OQYAmeqPx86jMD0eBekJZpcSsbKSY7FkZjL+\neIwBQJdiAJBpOvuG8NaJNnb/GOC2hVl452wHznVwkRh6FwOATPNG1fB0xez+Cb4ti4ZD9rWj502u\nhEIJA4BMs+1YEzISY7B0ZorZpUS8gvQEzM1MwKsVDAB6FwOATNE/5MHOqhbcsiATNi5YYojNi7Jx\n4FQbXD0DZpdCIYIBQKZ4s9aFvkEP+/8NtHlhFrz67plXRAwAMsUrRxqRGOvA6gKu/WuU+dmJmOWM\nwx/YDUQ+DAAy3IDbg21Hm3DbwizEOLj2r1FEBJsXZqGk1oXOviGzy6EQwAAgw+2pdqF7wI07lmSb\nXYrlbF6UBbdXsf04u4GIAUAmeLn8HJKnRWFdUZrZpVjO0pkpyE6OxdYjjWaXQiGAAUCG6h/yYNux\nJmxemMW5/01gswluX5yNXdUt6OgbNLscMhn/B5KhdlW3oHfQgzuWsvvHLHcty8GQR3lNADEAyFgv\nlzdielwU1hQ4zS7FshblJKEgLR7/+/Y5s0shkzEAyDA9A25sO3YeWxZnw8HuH9OICN6/dAbeOtmK\n8539ZpdDJuL/QjLMqxXn0T/kxb3X5JhdiuXduWwGVIcH5Mm6GABkmBcO1WOWMw7X5E03uxTLK0xP\nwOKcZLz0DgPAygwLABH5LxFpFpEKo9qk0HGu4wJKT7Ti7uU5EOHcP6HgrmUzUF7fidpmrhRmVUYe\nAfwcwGYD26MQ8j9vN0AVuHs5u39CxZ3LZsBuE/z24FmzSyGTGBYAqrobQJtR7VHoUFW8eKgBK2ZN\nxyxnvNnlkE9GYixunJeBFw41wO3xml0OmYBjABR0Rxo6UdPcg7s5+Bty7lsxEy3dA9hV3WJ2KWSC\nkAsAEXlERMpEpKylhX+UkeDZ/WcRG2XDHUtmmF0KXWZTcQbSEqLxXBm7gawo5AJAVZ9Q1RWquiI9\nPd3scugq9Qy48dLbDbhjyQwkT4syuxy6TJTdhruX52B7ZTMXirGgkAsAiiy/f+ccegc9eGBlntml\n0BXctyIXbq/ifw43mF0KGczI00CfBVAKYJ6I1IvIw0a1TeZ5dv8ZzMtMxDV5XPc3VM31/X6e2XcG\nXq+aXQ4ZyMizgB5Q1WxVjVLVmar6lFFtkzkqGjpRXt+Jj6zK47n/Ie7ja/Jx0tWLvbUus0shA7EL\niILmV/vPIMZhwwd47n/I27I4C874aPx36SmzSyEDMQAoKDr6BvHioQbctYyDv+EgxmHH/Stzsf14\nM8629ZldDhmEAUBB8ez+s7gw5MGfrJ9tdinkp4+smgUB8My+M2aXQgZhAFDADXm8eLrkFNYVOVGc\nlWR2OeSnnJRpuHl+Jn5z4AwuDHrMLocMwACggPtDxXmc7+rHw/z0H3YeXj8b7X1DeJ7zA1kCA4AC\nSlXx1N6TKEiLxw1zM8wuhyZp5exULMtNwX/uOcn5gSyAAUABtf9kG94524H/ty4fNhtP/Qw3IoI/\nv74QZ9r68AeuGRzxGAAUUD96oxZpCTG4b0Wu2aXQFN2yIBMFafH46e46qPLCsEjGAKCAOXi6HXtr\nXfjkxgLERtnNLoemyG4TPLKxABUNXbwwLMIxAChgfvRGDVLjo/Hgas77E+7uviYHWUmx+OHrNTwK\niGAMAAqI8voO7KxqwcPrZyMu2mF2OXSVYhx2fPamIhw83Y6dXCsgYjEAKCB+sK0aSbEOfHzNLLNL\noQC579pc5KZOw/f+WMWjgAjFAKCrVlLrwo6qFvzFjUVIjOW0D5Ei2mHD526ai4qGLrx2tMnscigI\nGAB0VbxexTf/UImclGl4aG2+2eVQgN29PAeF6fH43h+reF1ABGIA0FX5ffk5VDR04Yu3zeWZPxHI\nbhN8aXMxapp7OEdQBGIA0JT1D3nw3deqsCA7CXct5ZTPkeqWBZlYX5SG72+rRnvvoNnlUAAxAGjK\nfryjFvXtF/DVO+bzqt8IJiL4+/cvQM+AG9/bVmV2ORRADACaktrmbvzHrjrcszwHawvTzC6Hgmxu\nZiI+tnoWfrXvDCoaOs0uhwKEAUCTpqr4yosVmBZlx9/dPt/scsggf33zXKTGx+BLvyvHEAeEIwID\ngCbttwfrse9kG768ZT7SEmLMLocMkhwXhW98YCGOnuvCT3fVmV0OBQADgCblTGsf/umlo1g5OxX3\nX8cJ36xm86Js3L4kG49vr0V1U7fZ5dBVYgCQ39weL/76ubdhE8H3P7SUA78W9fU7FyIh1oHPP/c2\nBtxcOSycMQDIbz/ZWYeDp9vxzx9YhJnT48wuh0ziTIjBY/cuQUVDF775SqXZ5dBVYACQX/bUtOAH\nr1fj/Utn4K5lM8wuh0x2y4JM/On62Xi69DReKW80uxyaIgYATehMax8+86vDmJORiG/fsxgi7Poh\n4G83F2NZbgq+9Lty1DZzPCAcMQBoXL0DbjzyizIAwBMfvxbxMZzqmYZFO2z48YPXIDbKjk/87ABa\nugfMLokmiQFAVzTg9uDPf3kQ1U3d+NEDyzHLGW92SRRiclKm4amHVsDVM4A/ffoALgxyUDicMABo\nTG6PF5979m3sqXHhsXuXYOPcdLNLohC1NDcFj9+/HOUNnfjUMwfRP8QQCBcMAHoPt8eLv32+HK8e\nPY+v3bGAC7zThG5dmIVv3r0YO6ta8MlfMATCBQOALtE/5MGnnjmEFw434Au3zMXD62ebXRKFiQdW\n5uFb9yzGruoWPPKLg+gdcJtdEk2AAUAXdfQN4hM/249tx5rwT3cuxGdvmmN2SRRmHliZh+/cuwR7\na1pw33+UorHzgtkl0TgYAAQAOHquE+//t704eLod/3r/Mq7uRVP2oety8dRD1+F0ay8+8OM38fbZ\nDrNLoitgAFicquLZ/Wdwz7+XYMit+M0n1+CuZVzcha7OjcUZeP5Ta+Gw2fDBn5TgJzvr4PVyYflQ\nwwCwsIaOC/j4f+3Hoy8cwbWzpuP3n12Pa/Kmm10WRYj52UnY+pcbcOvCTDz26nE8+OQ+nGjpMbss\nGkVUQzeVV6xYoWVlZWaXEXH6hzx4au9J/PuOWiiAR7cU48FVszi5GwWFquK5srP4xiuVGBjy4pPX\nF+BTNxQiLpoXFQaDiBxU1RV+bcsAsI7+IQ9eONSAH71Rg8bOfty6IBNfu2MBclM5sRsFX3N3P761\n9ThePNwAZ3w0PnVDIR5cNQvTou1mlxZRGAB0icbOC3i+rB5Pl56Gq2cAS3NT8OiWYqwucJpdGlnQ\nwdPt+MG2auytdSE1Phofvi4XD67K4wyzAcIAIHT0DWJHVTNeONSAvbUuqAIb5qThU9cXYk2hkxO6\nken2n2zDk3tO4PXKJgDAqtlO3LE0G1sWZSM1Ptrk6sJXyAaAiGwG8K8A7ACeVNVvj7c9A8B/vQNu\nlNd34tCZduysasbB0+3w6vBcLfdek4N7rpmJ/DTO5UOhp6HjAp47cBa/Lz+HEy29sMnw9BLri9Kw\ntjANS3OTOV4wCSEZACJiB1AN4BYA9QAOAHhAVY9d6WcYAO816PbiTFsfTrT0oK6lFydaenD0XBeO\nn+/CyFl2C2ckYVNxBm4szsCymSkc3KWwoKqobOzGqxWN2FPrQnl9JzxehQhQmJ6ARTOSsGBGEmY5\n4zHLGYe81DgGwxhCNQDWAPhHVb3Nd/9RAFDVb13pZyIpAFQVgx4vBt1eDHkUg+7h24MeDwbcXvQN\netDdP4Tufje6+t3o6Xeju38IrT2DaOruR3PXAJq7+9HaO4jRv7L0xBgUZyVied50LM9LwbKZKZjO\nw2eKAF39Qyg71Yby+k5UNHShoqET57v6L9kmPTEGGYkxSEuIgTMhGukJMUiNj0ZCrAMJMQ7ERzsQ\nH+O7HWNHtMOGaLsNUXYbohw2RNkFUTZbRH1ImkwAGBmfOQDOjrpfD2BVMBr68E9L4eoZgAKAAiPv\nl6oKBaAPBDGlAAAG3UlEQVQKKHT466g308sfx8XbI9vpxfvv2V714vdx8WeGvzfkGX7Tnyy7TeCM\nj0ZmUiyyk2OxNDcFGYkxmOWMQ0F6AgrS45EUGzXp5yUKB0mxUdhUnIlNxZkXv9fZN4TTbb041dqH\nM629ONPWB1fPIFw9A6hp6oarZxCDHu+k23LYBFF2Gxw2gQggIrCN+goMf7XJ8OM23xiazQYI3t12\nzBgZ45tjbTd6XC7absPWz22Y9OuYLCMDYKzX/J53RRF5BMAjAJCXlzelhgrSE5CWEAPIcKOjfzFy\n+fdk+Bd46WPD3xv5fYhvw8sfv/z53t323Z8XAFG+Tx3Ro75GXfwqiHHYEBftQGKsA4mxUUjyfY2N\nsnGwlmiU5LgoLIlLwZKZKWM+rqroHfSgp9+NngE3en3/egbc6B10Y8itGPJ6MTRyJO7x+j6gDd8f\n8ngv+UDn1eEPfV7fB0Cv1/c932PQ0ffHruc93xuz8Evv2g06IjEyAOoBjJ5XeCaAc5dvpKpPAHgC\nGO4CmkpD37pn8VR+jIjCnIggwdflQxMzciqIAwDmiMhsEYkGcD+Alwxsn4iIRjEsJlXVLSKfAfAa\nhk8D/S9VPWpU+0REdClDj5NUdSuArUa2SUREY+NsoEREFsUAICKyKAYAEZFFMQCIiCyKAUBEZFEh\nPR20iLQAOG12HWNIA+AyuwiTWX0f8PVb+/UDobsPZqlquj8bhnQAhCoRKfN3sqVIZfV9wNdv7dcP\nRMY+YBcQEZFFMQCIiCyKATA1T5hdQAiw+j7g66ew3wccAyAisigeARARWRQDwA8ikioi20Skxvd1\n+jjbJolIg4j8m5E1Bps/+0BElolIqYgcFZFyEfmwGbUGkohsFpEqEakVkS+P8XiMiPzG9/g+Eck3\nvsrg8eP1f15Ejvl+39tFZJYZdQbTRPtg1HYfFBEVkbA5M4gB4J8vA9iuqnMAbPfdv5J/BrDLkKqM\n5c8+6APwcVVdCGAzgB+KyNhLN4UBEbED+DGALQAWAHhARBZcttnDANpVtQjADwA8ZmyVwePn6z8M\nYIWqLgHwPIDvGFtlcPm5DyAiiQD+EsA+Yyu8OgwA/9wF4Gnf7acBfGCsjUTkWgCZAP5oUF1GmnAf\nqGq1qtb4bp8D0AzArwtSQtRKALWqekJVBwH8GsP7YbTR++V5ADdJ5KzjOeHrV9Udqtrnu/sWhlf6\niyT+/A0Awx/8vgOgf4zHQhYDwD+ZqtoIAL6vGZdvICI2AN8D8DcG12aUCffBaCKyEkA0gDoDaguW\nHABnR92v931vzG1U1Q2gE4DTkOqCz5/XP9rDAP4Q1IqMN+E+EJHlAHJV9WUjCwsELpzpIyKvA8ga\n46Gv+PkUnwawVVXPhusHwADsg5HnyQbwCwAPqao3ELWZZKxf5OWnzfmzTbjy+7WJyEcBrABwfVAr\nMt64+8D3we8HAD5hVEGBxADwUdWbr/SYiDSJSLaqNvre3JrH2GwNgA0i8mkACQCiRaRHVccbLwgp\nAdgHEJEkAK8A+KqqvhWkUo1SDyB31P2ZAM5dYZt6EXEASAbQZkx5QefP64eI3IzhDwnXq+qAQbUZ\nZaJ9kAhgEYCdvg9+WQBeEpE7VbXMsCqniF1A/nkJwEO+2w8B+N/LN1DVB1U1T1XzAXwRwH+H05u/\nHybcByISDeBFDL/23xpYW7AcADBHRGb7Xtv9GN4Po43eLx8E8IZGzsU1E75+X/fHTwHcqapjfigI\nc+PuA1XtVNU0Vc33/d9/C8P7IuTf/AEGgL++DeAWEakBcIvvPkRkhYg8aWplxvFnH3wIwEYAnxCR\nt33/lplT7tXz9el/BsBrACoBPKeqR0Xk6yJyp2+zpwA4RaQWwOcx/hliYcXP1/9dDB/x/tb3+748\nIMOan/sgbPFKYCIii+IRABGRRTEAiIgsigFARGRRDAAiIotiABARWRQDgIjIohgAREQWxQAg8hGR\nHSJyi+/2N0Tk8Sk+z04Rmee77RSRikDWSRQonAuI6F3/AODrIpIBYDmAqV7pWQSgxnd7CYAjAaiN\nKOAYAEQ+qrrbN5f/5wHcoKqe0Y+PN1uqqv6vb5tZABpGzYK6BEB5EMsmmjIGAJGPiCwGkA3Apard\nlz8+3mypoyzDpW/41wL4TWAqJAosjgEQ4eIaBs9geLWnXhG5bYpPtRRArO855/iej11AFJIYAGR5\nIhIH4AUAX1DVSgwv7/ePU3y6ZQBsIvIOgL/H8AySD43/I0Tm4GygRAHkmxZ6+VhdSEShhkcARAEi\nIokAvHzzp3DBIwAiIoviEQARkUUxAIiILIoBQERkUQwAIiKLYgAQEVkUA4CIyKIYAEREFsUAICKy\nqP8Du5QMg/KbqxcAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#get standard error from this info\n", "se = 0.25 / np.sqrt(7)\n", "#make 500 points from -5 sigma to 5 sigma \n", "x = np.linspace(-5 * se, 5 * se, 500)\n", "plt.plot(x, stats.norm.pdf(x, scale=se))\n", "plt.xlabel('$x - \\mu$')\n", "#prefix with r to allow curly braces\n", "plt.ylabel(r'$P(\\bar{x} - \\mu)$')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "btasrfreasf", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 2.3\n", "Plot the probability of the population mean with a sample mean of 3.5 and sample standard deviation of 1.2. Plot it for 3, 5, 10, and 25 samples. Be sure to create a legend" ] }, { "cell_type": "code", "execution_count": 45, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "bnaetrgfretas", "locked": false, "points": 10, "schema_version": 1, "solution": true } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEMCAYAAADqG+D0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8ldX9wPHPuSt7bzIg7CE7IoooOBEVZxW1uKHuOlpb\nW22t1lGtxaKto+rPRV24UHECFpQNYc+ws5MbErJzx/n9cW9CyCDr3psb+L5fr/tK8jznOc+5YvK9\n5znnfI/SWiOEEEJ0lKG7GyCEEKJnkgAihBCiUySACCGE6BQJIEIIITpFAogQQohOkQAihBCiUySA\nCCGE6BQJIEIIITpFAogQQohOMXV3A7wpNjZW9+nTp7ubIYQQPcbatWuLtdZx7Sl7XAeQPn36sGbN\nmu5uhhBC9BhKqf3tLSuPsIQQQnSKBBAhhBCdIgFECCFEpxzXYyBCCAFgs9nIzs6mpqamu5viNwID\nA0lJScFsNne6DgkgQojjXnZ2NmFhYfTp0welVHc3p9tprbFarWRnZ5Oent7peuQRlhDiuFdTU0NM\nTIwEDzelFDExMV3ukUkAEUKcECR4HM0T/z0kgAjhQ9rppHTePEo/+wzZTlr0dBJAhPAh6+uvk/fw\nI+T9/iEOvTu3u5sjfEgpxQMPPNDw89///nceffTRLtc7ZcoURo4cybBhw7jttttwOBxdrrO9JIAI\n4SPOykqsL71M6FlnEXzqeIpfegmnzAo6YQQEBPDJJ59QXFzs0Xo//PBDNmzYwObNmykqKuKjjz7y\naP3HIgFECB85/M23OKuqiJl5KzG33IqjpITKn37q7mYJHzGZTMyaNYvZs2d7tN7w8HAA7HY7dXV1\nPh3rkWm8QvhIxY8/YkpKImjUKLDbMUZGcvi77wg755zubtoJ5S9fbGFr7mGP1jm0Vzh/vnhYm+Xu\nvPNORowYwYMPPthqmcWLF3Pfffc1Ox4cHMyyZctavOb8889n1apVXHDBBVx55ZXtb3gXSQARwge0\n3U7l8uWEX3CB6xOi2UzIaadRtXIVWmuZIXSCCA8P5/rrr2fOnDkEBQW1WGby5MmsX7++Q/V+++23\n1NTUcN1117Fo0SLOPfdcTzS3TRJAhPCB2qwsnBUVBI8b13AseNzJHF6wAFt2NpbU1G5s3YmlPT0F\nb7r33nsZM2YMN910U4vnO9MDAdfK8mnTpvH5559LABHieFKzeTMAQcNPajgWNGYMANWZmRJATiDR\n0dFcddVVvP7669x8883NznekB1JRUUF5eTlJSUnY7XYWLFjAxIkTPd3kVskguhA+UL1pM4awMMxp\naQ3HAvr2RVks1Gzf0Y0tE93hgQce8MhsrMrKSqZNm8aIESMYOXIk8fHx3HbbbR5oYftID0QIH6jZ\nvJnAk4ahDEc+symTiYABA6jdvq0bWyZ8paKiouH7hIQEqqqqulxnQkICq1ev7nI9nSU9ECG8TNts\n1OzcSeDQoc3OBQwZTM227bIqXfRIEkCE8LK6g9lgsxHQf0Czc4GDh+A4dAh7YWE3tEyIrpEAIoSX\n1e3dA0BAv77NzgUMdAWV2p27fNomITxBAogQXla72xVALH1bCCDuvRjq9u/3aZuE8AQJIEJ4Wd2e\nPZji4zGGhjY7Z4yNxRASQt2+fb5vmBBd5LMAopR6QylVqJTa3Mr5SUqpMqXUevfrT43OTVFK7VBK\nZSmlfu+rNgvhCbV79mBp4fEVuDK0Wnr3lgAieiRf9kDeBKa0UWap1nqU+/UYgFLKCPwLuAAYClyj\nlGo+nUUIP1W3bx+WPn1aPW/p00cCyAnAW+ncJ02axKBBgxg1ahSjRo2i0IcTMnwWQLTWS4CSTlw6\nDsjSWu/RWtcB7wOXeLRxQniJ4/BhnIcPY0lpfaW5pU8fbDk5OOvqfNgy4WveSucOMHfuXNavX8/6\n9euJj4/3eP2t8bcxkFOVUhuUUl8rpeoT1iQDBxuVyXYfE8Lv2XJyADCnpLRaxtKnN2iN7eDBVsuI\nns9b6dy7kz+tRF8H9NZaVyilpgKfAQOAltKUtrrqSik1C5gFkNYobYQQ3aEuOxsAc0rrn3nMya5z\nttxcAvr180m7Tmhf/x7yN3m2zsThcMHTbRbzVjr3m266CaPRyBVXXMHDDz/ss+zOfhNAtNaHG32/\nQCn1b6VULK4eR+P+fwqQe4x6XgVeBcjIyJDlvaJb2bJdPRBLcjsCSE6r/1uL44Q30rnPnTuX5ORk\nysvLueKKK3jnnXe4/vrrPdXkY/KbAKKUSgQKtNZaKTUO1+M1K1AKDFBKpQM5wHTg2u5rqRDtZ8vO\nxhAaiiEiotUyprg4MJkaHncJL2tHT8GbPJ3OPdn9ASQsLIxrr72WVatWHX8BRCn1HjAJiFVKZQN/\nBswAWuuXgSuB25VSdqAamK5dCYLsSqm7gG8BI/CG1nqLr9otRFfYsrMxp6Qc85GCMhoxJyZiy5Ue\nyInAk+nc7XY7paWlxMbGYrPZ+PLLLznHhztc+iyAaK2vaeP8i8CLrZxbACzwRruE8Ka6nOxjTuGt\nZ05OlgByAnnggQd48cUW/9x1SG1tLeeffz42mw2Hw8E555zDzJkzPdDC9vGbR1hCHG+01thycgmd\ncHqbZc29elF5jN3mRM/njXTuISEhrF27tsv1dJa/TeMV4rjhKClBV1cfcwpvPXOvXtgLC9GyFkT0\nIBJAhPASW14+AOakxDbLmpOTXWtB8vO93SwhPEYCiBBeYi8sAMCU0I4A0qsXgIyDiB5FAogQXmIv\nqA8gbaeWMCfXB5A8r7ZJCE+SACKEl9jyC8BoxBQT02ZZkzt/UX2vRYieQAKIEF5iLyjAFB+PMhrb\nLGsICMAYEYGtQAKI6DkkgAjhJbaCfMwdyIxqSkjAXljkxRaJ7uStdO5//OMfSU1NJbTJhmW1tbVc\nffXV9O/fn1NOOYV9XtgyQAKIEF5iLyjElJDQ7vKmhISGcRNx/PFWOveLL76YVatWNTv++uuvExUV\nRVZWFvfddx+/+93vPHpfkAAihNfYCwowJXYggMTHSQA5jnkrnfv48eNJSkpqdvzzzz/nhhtuAODK\nK69k4cKFuLJDeY6sRBfCCxwVFTgrKzF3oAdiTkjAbrWi7XaUSX41veVvq/7G9pLtHq1zcPRgfjeu\n7U/43krn3pKcnBxSU12JzE0mExEREVitVmJjY9tdR1vk/1IhvODIFN6214DUM8XHg9OJ3WrtUOAR\nPYc30rm3pqXehqf3CZEAIoQX1K8oN7djDUg9U7wraNgLCiSAeFF7egre5Ol07q1JSUnh4MGDpKSk\nYLfbKSsrIzo6utPtbokEECG8wF5QCNCxQfSGtSCFXmmT8A+eTOd+LNOmTeOtt97i1FNPZd68eZx1\n1lke74HIILoQXnAkjUlHxkBcAUTWghz/HnjgAY/NxnrwwQdJSUmhqqqKlJSUhqnBt9xyC1arlf79\n+/OPf/yDp5/2/EZa0gMRwgts+fkYIyMxBAS0+xpjdDQYjQ29F3F88UY6d4BnnnmGZ555ptnxwMBA\nPvroI4/cozXSAxHCC+xFxa6tajtAGY2Y4uLkEZboMSSACOEFjuJiTHEdny5pSoiXfFiix5AAIoQX\n2K1WjDEdDyDm+Hhs8ghL9BASQITwMK019uJiTJ1YsGWKi8deJPmwRM8gAUQID3NWVqFrajDFtp3G\nvSljbAzOw4dxyta2ogfwWQBRSr2hlCpUSm1u5fx1SqmN7tcypdTIRuf2KaU2KaXWK6XW+KrNQnSG\nw+qanmlsxz4gTdX3WhweTrgnhDf4sgfyJjDlGOf3AmdqrUcAjwOvNjk/WWs9Smud4aX2CeERdvcf\nf1Nsx2ZhNb7GbrV6tE2i+3kjnXtVVRUXXnghgwcPZtiwYfz+979vOPfmm28SFxfHqFGjGDVqFK+9\n9lqX7tUSnwUQrfUSoOQY55dprQ+5f1wBpPikYUJ4mL3Y9ce/M4+w6q+xF0kP5HjjrXTuv/nNb9i+\nfTuZmZn8/PPPfP311w3nrr76atavX8/69eu59dZbPXpf8N8xkFuArxv9rIHvlFJrlVKzuqlNQrSL\nvdg1CN6pQXT3NfV1iOOHN9K5BwcHM3nyZAAsFgtjxowhOzvbY/W3xe9WoiulJuMKIKc3OjxBa52r\nlIoHvldKbXf3aFq6fhYwCyAtLc3r7RWiKYfVCgYDxqioDl9rrB8DkUdYXpP/5JPUbvNsOveAIYNJ\n/MMf2iznzXTupaWlfPHFF/z6179uOPbxxx+zZMkSBg4cyOzZsxvSu3uKXwUQpdQI4DXgAq11w2+Q\n1jrX/bVQKfUpMA5oMYBorV/FPX6SkZHh2d1ThGgHe7EVY1RUu/ZCb8pgsWAID5dHWMcpb6Vzt9vt\nXHPNNdxzzz307dsXcO1UeM011xAQEMDLL7/MDTfcwKJFi7r8HhrzmwCilEoDPgFmaK13NjoeAhi0\n1uXu788DHuumZgrRps6uAalnio1tGIgXnteenoI3eSOd+6xZsxgwYAD33ntvw7GYRrMAZ86c6ZUt\nbX0WQJRS7wGTgFilVDbwZ8AMoLV+GfgTEAP8251y2O6ecZUAfOo+ZgL+q7X+xlftFqKj7NZiTJ2Y\nwlvPFBsrs7COY55O5/7www9TVlbWbJZVXl5ew1a38+fPZ8iQIV1reAt8FkC01te0cf5WoNk0Aa31\nHmBk8yuE8E+OYiuW3r07fb0pNobqLVs82CLhbx544AFefPHFLteTnZ3NE088weDBgxkzZgwAd911\nF7feeitz5sxh/vz5mEwmoqOjefPNN7t8v6b85hGWEMeDI2lMOr4GpJ4xNhaHjIEcd7yRzj0lJaXF\nrWsBnnrqKZ566qku3+NY/HUarxA9krOyEl1b28VHWHE4q6pwemi/CCG8RQKIEB5UnwixM6nc69UH\nHxkHEf5OAogQHlS/fqMzebDq1QcfmcrrWa096jlReeK/hwQQITzoSBqTrk3jBddsLuEZgYGBWK1W\nCSJuWmusViuBgYFdqkcG0YXwoCOJFDsfQOo3opKMvJ6TkpJCdnY2RbLXSoPAwEBSUrqWclACiBAe\nZLcWu9KYREZ2ug5TdBQoJY+wPMhsNpOent7dzTjuyCMsITzIUWzFGB3dqTQm9ZTZjDEqSgbRhd+T\nACKEB3U1jUk9U0yMpDMRfk8CiBAeZLdau7QGpJ4pLlZSugu/JwFECA9yeKgHYoyNxVEsj7CEf5MA\nIoSH1KcxMXZiJ8KmTDGujLwy7VT4MwkgQniIs6ICXVeHKcYDYyCxsejaWpyN8icJ4W8kgAjhIfXT\nbruSxqRew97oMpAu/JgEECE8xOFeOe6JQfSGxYQylVf4MQkgQniIvSEPlgd6IPX5sKQHIvyYBBAh\nPMSjj7DqM/LKTCzhxySACOEhdmsxGI1dSmNSzxgVBQaDJFQUfk0CiBAe4rBaMUVHowxd/7VSRiPG\n6GhZCyL8mgQQITzEXlSM0QOLCOtJOhPh7ySACOEhnkpjUs8UEyMJFYVf82kAUUq9oZQqVEptbuW8\nUkrNUUplKaU2KqXGNDp3g1Jql/t1g+9aLUT7tJVIscpWxePLH+fU/57KhPcm8LdVf6PaXt1qeVNc\nrOwJIvyar/cDeRN4EXi7lfMXAAPcr1OAl4BTlFLRwJ+BDEADa5VS87XWh7zeYiHaQWvtzoPVcg+k\nylbFrd/dyhbrFi7qexF2p5252+ay69AuXjz7RQJNzXeGMzZKZ6KU8vZbEKLDfBpAtNZLlFJ9jlHk\nEuBt7UoAtEIpFamUSgImAd9rrUsAlFLfA1OA97zbYiHax1lejrbZWl0D8sTKJ9hcvJnZk2Zzdu+z\nATg9+XT+8NMfmL12Ng+d8lCza0wxMei6OpwVFRjDwrzafiE6w9/GQJKBg41+znYfa+24EH7hWFvZ\nrs5fzfzd85k5YmZD8AC4uN/FXDfkOv67/b9kFmY2u04WEwp/528BpKV+uj7G8eYVKDVLKbVGKbVG\n9j8WvnIkgBz9CMupnTy7+lmSQpKYOXxms+vuGX0PsUGxzF47u1nmXaN7QF7GQYS/8rcAkg2kNvo5\nBcg9xvFmtNavaq0ztNYZcXFxXmuoEI3V56xq2gNZfHAx20q2cffou1sc5wg2B3P7yNvJLMxkdf7q\no87V1yUzsYS/8rcAMh+43j0bazxQprXOA74FzlNKRSmlooDz3MeE8Av1aUyargOZu20uSSFJXJB+\nQavXXtL/EqIConh327tHHW8IILKYUPgpX0/jfQ9YDgxSSmUrpW5RSt2mlLrNXWQBsAfIAv4D3AHg\nHjx/HFjtfj1WP6AuhD+wW62uNCYREQ3HdpTsYHX+aqYPno7J0Pp8lQBjAFcOvJIfD/7IwfIjQ33G\nyEhJZyL8mq9nYV3TxnkN3NnKuTeAN7zRLiG6yl5chCkm5qg0Jp9lfYbZYOby/pe3ef3Vg67m9c2v\n81nWZ9w9+m6gcToTCSDCP/nbIywheiRHsfWorWztTjtf7/2aM1POJDKw7eSKCSEJjE8az1d7vjpq\nMN2VzkQeYQn/JAFECA9wpTE5Mv6xKn8V1horU/tObXcdF/W9iJyKnKOm9JpiY2UQXfgtCSBCeIC9\nuPioPFgL9iwgzBzGGSlntLuOs9POJsgUxIK9CxqOmWJjsBfLdHThnySACNFFWmtXD8Q9bdzmsLHw\nwELO7n02AcaAdtcTbA7mtF6nsfjg4obHWMaYWBzF1mZrRITwBxJAhOgiZ1kZ2GwNiwjXFa6jwlbB\nWalndbiuyamTKawqZKt1K+B6hFWfzkQIfyMBRIguql+FXr8G5H/Z/8NisHBK0ikdruuMlDMwKAOL\nDi4Cjqxsr19nIoQ/kQAiRBfVz5IyxboeYS3JXsLJSScTbA7ucF1RgVGMjh/N4oOLgUbpTGQtiPBD\nHQ4gSqkQpZTRG40RoidqyIMVF8u+sn3sP7yfM1PO7HR9k1Mns+vQLnIrchuCkszEEv6ozQCilDIo\npa5VSn2llCoEtgN5SqktSqlnlVIDvN9MIfxX/SwpU0wMS7KXAHRo9lVTpyefDsDy3OVHHmHJWhDh\nh9rTA1kM9AMeAhK11qla63hgIrACeFop9UsvtlEIv+YoLkaZzRjCw1mSs4T+kf1JDu38bgN9I/oS\nHxTP8rzlR9KZyFRe4Yfak8rkHK21relBdy6qj4GPlVJmj7dMiB7CXmzFGBdLnbOOzIJMrh58dZfq\nU0oxvtd4lmQvQRsUxpjohmy/QviTNnsgLQWPzpQR4njl2gs9jg2FG6hz1nFKYsdnXzU1Pmk8pbWl\nbCvZhikmVh5hCb/U7kF0pdQmpdRcpdTvlFIXKKVSlFJ/9GbjhOgJ6lehr8xfiVEZGZswtst1ntrr\nVMA9DhITI7sSCr/UkVlYZ+JKsV4NTAc2Axd6o1FC9CSuHkgsq/JWMSxmGKGW0C7XGRsUy4CoAazI\nXeHOhyUBRPifdgcQrXWJ1vpHrfUcrfUNwMnALu81TQj/px0OHCUlOKMj2Fy8mXFJ4zxW96lJp7oS\nK0ZHSjoT4Zc68gjrqOm6WutdwAiPt0iIHsRx6BA4neRaqrBre6dWn7cmIyGDOmcdRUE2VzqT8nKP\n1S2EJ3RkQ6lXlVL9gBxgIxAIbFZKhWitK73SOiH8XP3YxE5VgNlgZlTcKI/VPSZhDAB7jCUMxj3b\nKzzcY/UL0VXtWUioALTWk7XWacDVwFe4tp0NBjKVUju82koh/FR9jqp19r2Mih9FoCnQY3VHBEQw\nIGoAW5w5gKQzEf6nXQsJlVJ3K6XSALTWB7TW84G/Af/GtZjwX15soxB+q35we5M+SEZChsfrHxs/\nlnW23e57yVRe4V/aE0CmAA7gPaVUnlJqq1JqD64B9OnAc1rrOd5spBD+qn6/8pJgJ6PiPff4ql5G\nYgYFgbWAZOQV/qfNMRCtdQ2unsa/3SvOY4FqrXWptxsnhL+zFxVjDzBjC1CMjBvp8frHJoylPAi0\nQclUXuF32jMGcoNSqlgpVQK8BlR0NngopaYopXYopbKUUr9v4fxspdR692unUqq00TlHo3PzO3N/\nITzNXlxMeZiRQVGDCDGHeLz+2KBYekemUxVqlnQmwu+0ZxbWI8C5uGZf3Q086f7aIe4U8P9y15UN\nrFZKzddab60vo7W+r1H5u4HRjaqo1lp7/hmBEF1gKy6mKMjmlcdX9cYmjMUatJf4IkmoKPxLe8ZA\nDmutM7XWhVrrR4DOrpQaB2RprfdoreuA94FLjlH+GuC9Tt5LCJ+oKsihJNjJ6PjRbRfupLEJYykJ\ndlJZkOu1ewjRGe0JIElKqVlKqYlKqTigs5l3k4GDjX7Odh9rRinVG0gHFjU6HKiUWqOUWqGUurST\nbRDCoxzFVspC8HoAKQuBuuJCr91DiM5ozyOsP+NacX4dMBwIVUotADYAG7XW7e0lqBaOtZabYTow\nT2vtaHQsTWudq5TqCyxSSm3SWu9udhOlZgGzANLS0trZNCE6TtfVYaqoxh4VTmJIotfu0yu0F7bI\nEAzby9Fa416aJUS3a08691e11ndprc/UWkcDfYEXgVJgagfulQ2kNvo5BWitTz6dJo+vtNa57q97\ngB85enykaXsztNYZcXFxHWieEB1jcw9qR/bq4/V7hSelYbI7cRw+7PV7CdFeHd4TXWudrbVeoLX+\nm9Z6RgcuXQ0MUEqlK6UsuIJEs9lUSqlBQBSwvNGxKKVUgPv7WGACsLXptUL4Uv7BbQAkpQ7x+r0S\nUgYBkOu+pxD+oMMBpLO01nbgLuBbYBvwodZ6i1LqMaXUtEZFrwHe10enHh0CrFFKbcC1xe7TjWdv\nCdEdsnavAaBv3zFev1efPq5ZXjuyVnj9XkK0V0eSKXaZ1noBsKDJsT81+fnRFq5bhmv8RQi/kbtv\nC0lAeno7AojTCVk/wPYvoDjLdSxuIAy9FPpOgjbGNdL6jGQ/kL1vUxdbLYTn+DSACHE8Kc3Zg1YQ\nEJ9w7IJ5G2H+3ZC3HgIiIPEk0Bo2fwJr34TUU+Ci5yFhaKtVBCQlAWA9mOXBdyBE10gAEaITDtcd\nhiIrtvBglPkYM9s3fADz74KgaLj0ZTjpCjBZXOdsNbDxfVj4OLx2Nlz6EgxreYa6ISwMh8WEs6iI\nstoyIgIivPCuhOgYn42BCHE82Vi0kagKMMYfY6bfmjfg01mQNh5uXwajrjkSPADMgTD2Rrj9Z0gc\nDvNugk3zWqxKKYWKiyG6HNYXrvfsmxGikySACNEJmYWZRJdDaK9W1hptnQ9f3gcDzodrP4KQmNYr\nC0uEGZ9C2mnwyUzY9X2LxYKTkomugHWF6zzwDoToOgkgQnRCZmEmcZUGAhJ7NT+ZtwE+/RWknAxX\nve3qabTFEgLXfgAJw2DeLUcG2hsXSUgiocrCugIJIMI/SAARooNsThvb8jYSUunA1PQRVl0lfHQT\nBEXB9P+2L3jUCwh1XWM0wYczwF571GlTfDwRh+1sLt5Ejb3GA+9EiK6RACJEB+0o2UHgYdcfcHNC\nkxlY3z0CJXvgslcgNL7jlUemuQbbC7fC4iePOmVKiMdocxBQZWeLdUtnmy+Ex0gAEaKD6sc/wNUr\naJC1ENa8DqfeCekTO3+DgefBmOth2Rw4sLLhsNl9r+gKVxuE6G4SQITooMzCTPrbo4BGAcReCwt+\nAzH94axHun6T856A8GT44tfgsLnu5e7tDHMmyTiI8AsSQIToAK01mYWZDHW6su/W/1Fn2RzXo6up\nz3Zs3KM1geEw5Wko2garXzvqXifRi/WF63FqZ9fvI0QXSAARogOyK7Ipri4m3RaBMpsxRkbCof2w\n5DkYegn0O8tzNxt8oau+xU9BRREmd3bpfrYoym3lZJXKqnTRvSSACNEB9WMP8ZUmTPHxrr05Fv7F\nlcvq/CfbuLqDlIIpfwNbJSz8C4aAAIyRkSRWuRYjZhbIOIjoXhJAhOiAdQXrCDOHEVRa7Rr/yNsI\nmz+G8XdARIrnbxg3EMb9CtbPhaIdmOLjCSitJD4oXhYUim4nAUSIDlhfuJ6R8SOxFxS6xiQWPuZa\n8zHhHu/ddOL9YA6GRX/FlJCAvaCQ0QmjZSaW6HYSQIRop7LaMnaX7WZM/BjshYWYAm2Q9T2cfj8E\nejG5YUgsnHoXbJuPKcSAvbCQ0fGjyavMI68iz3v3FaINEkCEaKf6JIajQgbhrKrCXLIGwnrBuJne\nv/mpd0JQNOaKzdiLixkTPRKQvFiie0kAEaKd1hWuw2QwMdDuSoxoqtsHZz4I5iDv3zwwHCbej6lm\nNzidpDujCTGHyGMs0a0kgAjRTusL1zM0eijGokMAmONjYPQvfdeAjFswR4cCoAsKGRU3SnogoltJ\nABGiHWodtWwq3sTo+NHYNi0FwDzpZjAeYzMpT7MEY5pwDQC2zT8zOn40WYeyKKst810bhGhEAogQ\n7bDVuhWb08bo+NHY1y4ApTFN8sHYRxPmc24HwLbyE8YkjEGj2VC0weftEAJkS1sh2qV+rGGUNlGb\nfQBTZBwqOLzd1x+usbFu/yF2FpRjrazDoBTRwRYGJ4UxJi2KkID2/SoaY5IwBJmx7dvJSQ4DJmUi\nszCTM1LO6NT7EqIrfBpAlFJTgH8CRuA1rfXTTc7fCDwL5LgPvai1fs197gbgYffxv2qt3/JJo4XA\nteq7d3hvYla+xv6aAMxp/dq8RmvN0l3F/HflARZuL8Dm0ABYjAacWmN3un4OMBmYPCieWyamc3Kf\n6DbrNaekYaupImjZiwyNGSqJFUW38VkAUUoZgX8B5wLZwGql1Hyt9dYmRT/QWt/V5Npo4M9ABqCB\nte5rD/mg6eIE59RO1hetZ1LcaFj6Njb7AIJSUo95zfqDpTz51TZW7SshNtTC9af24azB8ZyUHEF4\noOvXrrTKxubcMhZuK2T+hly+2ZLPKenRPHbJSQxKDGu1bnOvFGy7rbDlU0afdSfv7f+GOkcdFqOl\n1WuE8AZf9kDGAVla6z0ASqn3gUuApgGkJecD32utS9zXfg9MAd7zUluFaLC3bC+ltaWMLj6INgVh\nK60hvFdSi2VrbA5m/7CT/yzZQ0xoAI9fehJXZ6RiMTUfbowKsTBxQBwTB8TxuymDeX/1AeYs3MWF\nc5Zy+6SsNhKNAAAgAElEQVR+/PrsAZiMza8z9Uqien0mmAIZnb+Lt5x1bLVuZVT8KI+/dyGOxZeD\n6MnAwUY/Z7uPNXWFUmqjUmqeUqr+Y157rxXC41bnrwbg5N0/Yx8wHex2TEnNA0heWTVXv7KcV/63\nh6tPTmXRA2cyY3zvFoNHU0EWIzdNSGfhA5OYNqoXLyzK4trXVlJwuPnWteakXjjKDuM86ZeM3rEI\ngLUFa7v4LoXoOF8GENXCMd3k5y+APlrrEcAPQP04R3uudRVUapZSao1Sak1RUVGnGytEvdX5q0k0\nBJJid2JPuRAAc5MAsu7AIS5+4WeyCit4ZcZYnrp8BGGBHZ/iGx1i4R9XjWL21SPZlF3Gpf/6mR35\n5UeVqb+3rc9lRGvoY5QFhaJ7+DKAZAONHxynALmNC2itrVrrWveP/wHGtvfaRnW8qrXO0FpnxLn3\nTxCis7TWrMlfxcnlh1CjrsFWbgfA3KtXQ5mfdhVz3X9WEhJg5LM7J3D+sMQu3/ey0Sl8fPtpOJya\nK19exoo91oZzZvfjM1u5hpHTGVNWTGbBOtlgSvicLwPIamCAUipdKWUBpgPzGxdQSjX+WDcN2Ob+\n/lvgPKVUlFIqCjjPfUwIr9pdupuS2lJOrq6GCfdiy3UlL6zvBXy3JZ+b31xN75hg5t12GgMSWh/8\n7qihvcL59M4JJIQHcuP/rWJZVrHr3u7gZcvLhdPvY3RNNYdt5ewp3eOxewvRHj4LIFprO3AXrj/8\n24APtdZblFKPKaWmuYvdo5TaopTaANwD3Oi+tgR4HFcQWg08Vj+gLoQ3rT64BICM1DMhph+2vDwM\nISEYwsJYvL2QO+auY2ivcN6fNZ64sACP3z85MogPZo2nd3QIN7+1muW7ra6dCY1GbHl5ENOPse41\nIOtyfvL4/YU4Fp+uRNdaL9BaD9Ra99NaP+E+9iet9Xz39w9prYdprUdqrSdrrbc3uvYNrXV/9+v/\nfNluceJaveMTEu12UiY+BIAtNxdzryTW7D/Ebe+uZUhSOO/cMo7IYO9NoY0JDWDuzFNIjQrm5jdX\nk5lbjikhHluOa7lUyhkPEWt3kLn9U6+1QYiWSCoTIVqhaytYU7GPceYYVNJwAGwHDlATm8jNb64m\nOSqIN286uVOD5R0VGxrAf2eOJz48gFvfWoMjPgnbwWwAVNIIRpsjySzLgrpKr7dFiHoSQIRoxe4V\n/+SQQZEx6HLANaBee/Ag3x4yERpg4p1bTiEm1POPrVoTFxbAWzeNQwFLKwOoOXBkZvuYfheSYzKQ\nv+JfPmuPEBJAhGiJvZZVm+cCcPLQqwCoyCuAmhpygmN455ZxJEf6YB+QJvrEhvD6jSezPyAKbS2m\n/NBhADKGXAHAqo1vgr32GDUI4TkSQIRoyYb3Walq6RUQTXJoMlprnn97MQCXXTSO/vGem23VUaNS\nI5k6JQOAp95YiNOpGRg1kGhzKMtUDWx4v9vaJk4sEkCEaMphx/bTP1gVHMKpaZNRSvHioiz2bdwJ\nQMb44d3cQBh9ykkA7N2wk+cX7sKgDJyacibLQ0Jx/vQPcNi7uYXiRCABRIimtn7Gpqo8KpRmQvIE\nvtuSz3Pf72RyuA2UwpzS/Vl0LKmudbXnRzuYs3AXX23M47Tk0yhRmp2VObD1s25uoTgRSAARojGt\nYek/WBaTjEEZiDYM5b4P1jMyJYIzQmsxJSVisHR/1ltjZCSG8HDOjbQxtncUD3y0nig1DIBlsamw\n9B+u9yKEF0kAEaKxnd9C4RaWR8UzJOok7ntvJyEBJl6ZkYEzOxtLG2ncfcmSkoIjO5uXfzmW6GAL\nv/tgH+nh/VkWkwKFW1zvRQgvkgAiRD2tYenfKYtKY3N1PkVFvckvq+HlGWNJjAik7uBBzGn+E0DM\naWnUHTxAXFgAr16fQUlVHaXWdNZV5VAdkQZL/y69EOFVEkCEqJe1ELJXs+KkC3FqJ3sPpvDk5cMZ\nkxaFs7ISh9WKJTWtu1vZwJKaii0nF+1wcFJyBH//xUhy8lKxOW2sGTkNslfDPklvIrxHAogQ4Pqk\nvvgJiExjbnkt2hHIjWMmcuXYFADqsl2rvi2pKd3ZyqOY01LBbnflxAIuGtGLmSefjXaaeOOwA0Li\nYelz3dxKcTyTACIEwM5vIHcdu4fewTrrSqLUSfxh6rCG03V79wJgSU/vrhY2E9CnD3CkbQAPnjec\nSMNgVhWtYO/AG2DPYjiwoptaKI53EkCEcPc+7BF9uDrTgDKVc8/4S4/aTrZ2925QCov7j7Y/sPTv\nD7jb5mYwKGaOvRiDxcr0LYk4guNg4WMyFiK8QgKIENu+gPxNzHFcRm3ANgzKyHl9Jx1VpG7PXsy9\nemEI8n36ktaYoqIwRkVR1yiAAJyffhYAVUG7eFlfDvt/ht0Lu6OJ4jgnAUSc2JwO9I9PUWBO5d/W\nMST22k1GwlgiAiKOKla7Zw+Wfn27qZGts/TrS+3uozeSSgxJZGjMUNLT9jGndALFpkS09EKEF0gA\nESe2De+jCrfy18pL+NU50eRX72NS6qSjimink7q9ewlI978AEtC3H3W7d6ObBIfJqZPZV7GNO85L\n5amqS1F5G2Dr593USnG8kgAiTlx1VVR/+xfWO/tiGfkLEhJdn+SbBhB7Xh66psYveyAB/friKCvD\nUXL0Bp2TUyej0aQk74WRV7HLmUzFN49KjizhURJAxAmr4LvZBNUUMC/6Vzx5xXAWH1xM/8j+pIYd\nvViwdo8rsAT09b8AYunXfCAdYGDUQHqF9GLxwcU8cfko5kXeRGj5XvJ/fLU7mimOUxJAxAmpKD+b\n0DUvsMQwjl/fchOH60pYW7CWc3qf06xsbZbrj7OlXz9fN7NNAe5eUd2eo8dBlFKc3ftsluUuo05X\ncvMtd7FWDSNw6VOUWou6o6niOCQBRJxwqursrH7zQQJ0LUlX/o24sAC+2/8dGs0FfS5oVr52xw6M\ncbGYoqK6obXHZkpMxBAc3BDkGpuaPhWb08bC/QtJiAgieNrfCNfl/PTGb6m1O7qhteJ4IwFEnFDs\nDid/+78POb96AbkDrmXA0DEAfL33awZFDaJvZPPHVDU7dhA4aLCvm9ouSiks/ftTu2tXs3PDYoaR\nGpbKgr0LABgyeiIH+lzB+RXzefrdL3E6ZVaW6BqfBhCl1BSl1A6lVJZS6vctnL9fKbVVKbVRKbVQ\nKdW70TmHUmq9+zXfl+0WxwetNY98upFLcp6jLiCKtCueACCnIocNRRuYkj6l+TU2G3VZWQQOHuTr\n5rZb4ODB1Gzb1mwmllKKC9IvYFX+KoqriwHo84un0KZATt89mye+2todzRXHEZ8FEKWUEfgXcAEw\nFLhGKTW0SbFMIENrPQKYBzzT6Fy11nqU+zXNJ40Wx5U5C7OwZ85ljCGLoKlPQlAkAN/uc6U9vyC9\nhcdXe/aibTYC/LQHAhA4dAjOw4ex5+Y2Ozc1fSpO7Wx4j4TGYz7r95xtzCRn+Ye8tnRPs2uEaC9f\n9kDGAVla6z1a6zrgfeCSxgW01ou11lXuH1cA/pO5TvRoc1fu540f1vHnwA/QaafCyOmAq1cyP2s+\nI+NGkhzafKfB2h3bAfy7BzJkCAA127Y1O9cvsh8DowY2PMYCUOPvQCcO5+mgd/jnV2v5fH2Oz9oq\nji++DCDJwMFGP2e7j7XmFuDrRj8HKqXWKKVWKKUu9UYDxfHp47XZPPzZZl6MmUeIswI19e+gFAAb\nijawu2w3lw+4vMVra7bvQFksfpVEsamAgQPBYKBma/MAAnBR34vYWLSR3aXugXajCXXxHCKch3gu\n+jPu/3AD32zO82GLxfHClwFEtXCsxVE8pdQvgQzg2UaH07TWGcC1wPNKqRbnVCqlZrkDzZqiIpmu\neKL7amMev523gTt6ZTGx8jvU6fdC4kkN5z/N+pQgUxDn9zm/xetrtmwhYOBAlMnkqyZ3mCEoCEt6\neos9EIBp/aZhMpj4eNfHRw4mj0GdchvnVX3J1fHZ3P1eJgu3FfioxeJ44csAkg00XqGVAjR7aKuU\nOgf4IzBNa11bf1xrnev+ugf4ERjd0k201q9qrTO01hlxcXGea73ocb7fWsCv38/kjFQTD9T+G+KH\nwpm/azhfaavk671fc0H6BYSYQ5pdrx0OajZtImjkSF82u1MChwyhZsuWFs/FBMVwVupZzN89n1pH\n7ZETk/8IEak8rl5idKKF299dx/92yocu0X6+DCCrgQFKqXSllAWYDhw1m0opNRp4BVfwKGx0PEop\nFeD+PhaYAMgUEtGqrzbmcfu7axnWK5z/xH+MobIILv03mAIaynyz9xuq7dVc1v+yFuuozdqNs6qK\noJEjfNXsTgsaORJ7YWHD5lJNXTnwSspqy/hh/w9HDgaEwqUvYTy0l3dS5tM/PpRZb6/hxx2FLdYh\nRFM+CyBaaztwF/AtsA34UGu9RSn1mFKqflbVs0Ao8FGT6bpDgDVKqQ3AYuBprbUEENGiD9cc5O73\n1jE6LZL3Tz2AefMHMPEB6HWk06q15t1t7zIwaiAj41ruYVRvWA9A0IgeEEBGjQKgOjOzxfOnJJ1C\nSmgKH+748OgT6RNhwj0EbHiLDyaV0j8+lJlvr+GLDc1ndAnRlE/XgWitF2itB2qt+2mtn3Af+5PW\ner77+3O01glNp+tqrZdprYdrrUe6v77uy3aLnuPNn/fy4LyNTOgfy9uXRBP07W8h7dSjHl0BLM9b\nTlZpFjOGzkCplobnoHrjRowREZh7927xvD8JHDwIFRhI1fr1LZ43KAPXDL6GdYXr2Fi08eiTk/8I\nicMJ+/Ze3p+exujUKO55P5P/rjzgg5aLnkxWoovjgtOpeXLBNh79YivnDU3gtetOIujzW8FogSte\nB+PRg+DvbH2HmMAYpqZPbbXOmg0bCBwxotUA40+U2UzQ8OFUZ7YcQACuGHgFYeYw3try1tEnTAFw\nxRtgryFs/i28dcMoJg2M4w+fbuKFhbuaLVAUop4EENHjVdc5uH3uWl5dsocZ43vz72tHE/D1A5C/\nCS59CSKOni2+u3Q3P+X8xPTB07EYLS3Wabdaqd2VRXBGhi/egkcEjR5NzbZtOKurWzwfYg7hF4N+\nwQ8HfuBg+cGjT8YNhEv+BdmrCVr0MK/MyOCy0ck89/1O7v9wAzU2yZ0lmpMAInq0gsM1XP3qcr7b\nWsCfLhrKY5cMw7TyRdjwHkz6Awxqnp7klQ2vEGQK4qpBV7Vab9WqVQCEjD/Fa233tOCMsWC3U7Vu\nXatlrh18LQZl4M3NbzY/OexSOO1uWP0alo3v8o+rRvKb8wbyaWYO1722kuKK2ubXiBOaBBDRYy3L\nKubCOUvJKqzgPzMyuPn0dNTOb+D7P8Owy+DMB5tds+vQLr7Z9w3XDbmO6MDoVuuuXLESQ0gIgcOG\nefMteFRwRgaYzVQuW9ZqmYSQBC7vfzmf7PqE7PLs5gXOfhT6ToYv70PtXsRdZw3g39eNYUtuGZe8\n+DPrDhzy3hsQPY4EENHjOJ2aFxbu4pevryQy2MLnd07gnKEJsH85fHQTJI2ES/7dsNq8sZc2vESw\nOZgbh914zHtUrVhB8Mkn+/UCwqYMwcEEjxpF5bLlxyw3a8QsjAYjL214qflJowmuehviBsOH10Pe\nRqYOT+KjX52GUnDVy8t5bekeGRcRgAQQ0cPklVVzw/+t4rnvd3LxyF58fucEBiSEQd5G+O/VrvGO\n6+aBJbjZtZuKNvH9/u+ZMXQGEQERrd7DlpND3f79BJ/Scx5f1QuZcBq127Zht1pbLZMQksA1g6/h\nyz1fHklv0lhgOFz3EQRGwLtXQNFOhqdE8NU9EzlnSAJ//WobM99eg1UeaZ3wJICIHkFrzcdrszlv\n9hLW7DvEE5edxPNXjyIkwAQFW+HdyyEgDGZ8BqHNMxA4tZMnVz5JXFBcm72P8oWLAAibPMkL78S7\nQiZMAKDy55+PWe7mk24mxBTC06uebrk3Ed4LZnwKaHjrIijeRUSQmZd+OYZHLx7Kkp3FnDd7CV9t\nlBxaJzIJIMLv5ZZWM+udtTzw0QYGJ4bx9a8nct0pvV3Ta3PWwptTwWCC6z+HyNQW6/gs6zM2Wzdz\nf8b9LaYtaax80SIs/fph6dPHC+/GuwKHDcMUH0/5998fs1xUYBR3j7mbFXkrjqR6bypuENzwJTgd\n8OZFULgdpRQ3Tkjni7tPp1dkEHf+dx13zF1LUbn0Rk5EEkCE36q1O/jX4izOfu5/LNlZxB+mDub9\nWafSJ9YdAPYugbemQUA43PQ1xPZvsZ7i6mKeX/s8Y+LHcGH6hce8p6O0lKrVqwk7+2xPvx2fUAYD\nYeeeS8WSpTgrK49Z9qqBVzEkegjPrH6G8rrylgvFD4YbvwTthDfOg/2uAfpBiWF8esdpPDhlED9s\nLeSs537kjZ/2YnM4Pf2WhB+TACL8jtaa77bkM+X5pTz77Q7OGBjLD/efyawz+mE0uAfG1/wfvHMZ\nhCfDzd9AdMvp1rXW/GXZX6i0VfLI+EfaXBRYvnAhOByEndMzAwhA2PnnoWtrqVi69JjljAYjj4x/\nBGuNladWPtV6wfghcOv3EBIPb18Km11ZfU1GA3dM6s+CX09kdFoUj325lan/XMoSSch4wpAAIvzK\nst3FXP7SMma9sxYFvHXzOF6ZkUFqtHtQ3F4HC34LX94L6WfCLd+5nte34rOsz/gx+0fuGXMP/aNa\n7qE0Vvrpp1jS0wkcPtxD78j3gseOxRgXS9nnbe/8PDxuOL8a8Su+2PMFX+/9uvWCUX1c/617jYZ5\nN8O3fwSHDYD+8aG8ddPJvHZ9BrV2J9e/sYrpry5n1d4SD70j4a8kgIhup7Xmxx2FXPPqCq79z0ry\ny2p4+vLhfHffGZw5sNGAeHEWvH4urHoVxt8B137YsC1tS7ZZt/Hkyic5OfFkZgyd0WY76vbvp3rN\nWiIuu6xHpC9pjTIaibz0Mir+9z9sBW3v8TFrxCxGxI3g8eWPs69sX+sFg6Phhvlw8kxY/iK8dTGU\nuXYzVEpxztAEvr//DB69eCi7iyq56pXlzHh9Jav2lsi03+OUOp7/YTMyMvSaNWu6uxmiFXV2J19u\nzOXVJXvYnl9OQngAMyf25ZfjexNoNh4p6HTCurfg2z+48jZNewGGXHzMuq3VVqZ/5dq29r0L3yM2\nKLbN9hQ+9xzW19+g/+JFmBMSuvTeulvdgQPsPu984n59D7G3395m+ZyKHK796lrCLGHMnTr3mNOc\nAdg0D+bf45q8cP4TMPqXR627qa5z8O6K/bz8v91YK+sYkRLBLaenM3V4EmajfG71Z0qpte7N+9ou\nKwFE+FpWYQUfrD7Ax+tyKKmsY2BCKDMn9uWSUclYTE3+uBRsgS/vh4MroM9EuOyVZrmtmqqoq2Dm\ndzPZVbqLty94m6ExQ9tsk6OigqzJZxEyYQIpz8/uytvzGwduvoXaXbvo98P3GAIC2iy/rmAdt353\nKyPjRvLSOS8RaAo89gXW3TD/btj/M/Q7C6b+HWKO3ii0us7Bx+uyeePnvewpqiQxPJCrMlK4cmwq\naTHN1+qI7icBxE0CiP+wVtTy7ZYCPsvMYdW+EkwGxdlD4rlmXBpnDIjDYGjyyKiiEJb8Hda87ppl\ndd7jMPJaMBz702uVrYrbf7idjUUbmT15NpNSJ7Wvfa+/TuGzf6fPvHkEndRz0pccS+WKFRy48SYS\n/vQI0dde265rvtrzFQ8tfYjxSeOZc9actoOI0+n6N/rhUbDXwrhZcOZvISiqSTHNjzsLeXPZfpbu\nKkJrOCU9ml9kpHLesATCA82dfJfC0ySAuEkA6V5F5bV8v7WABZvyWL7HisOp6Rsbwi8yUrlibDLx\nYS38caoqgWUvwMqXXX+QxsyAs//sev7eBmu1lbsX3c0W6xaeOeOZVvc5b8pRVsbu86cQOHQoaW8c\nP1vNaK3Zf90vseXm0u/rBRiCgtp13WdZn/Gnn/9ERmIGsyfNbvtxFkB5ASz+K6x7x7WS/ZTbXK8W\n/t3yyqr5ZF0OH605yD5rFWajYkL/WKYMS+TcoQnEhLbdWxLeIwHETQKIb9kcTtbtP8T/dhaxZFcR\nm3MOA5AeG8LU4YlMHZ7E0KTwlgeoi3bAipdgw/tgr4HhV8Kkh5o9EmnNjpId3Lv4Xoqri3n6jKc5\nO63903Dzn3ySQ+/OJf2TjwkcPLjd1/UEVatXs3/G9cT86lfE33dvu6/7cs+XPPLzI6SEpvDCWS/Q\nJ6JP+y7M3wQ/Pg3bvwRzCIy9EU6+pcV/R601mQdL+WZzPl9vzuNgSTUGBaNSIzl9QBwTB8QyKjVS\nxkx8TAKImwQQ76qotZN54BBr97te6/YforLOgcmgGNM7ijMHxjFpUFzrQaPmMGz9HDZ+APuWgjEA\nRlzlmmGV0Pa4BbhSlHy440OeXf0s4QHhzJk8h+Fx7Z+CW7V2LftnXE/kVb8g6dFH231dT5L7u99T\ntmAB6R992KEAubZgLfcuvpdaRy2/yfgNvxj4i/bPTivcBj/Ndg22aweknwFjboBBF4CleSYArTXb\n8sr5Zks+S3YWsTG7FKeGEIuRU/rGMLZ3FKPTIhmZEulKXyO8RgKImwQQzzlcY2Nb7mG25R1mW145\nm3LK2J5/GKd2Tb4ZnBjO2N6RTBwQx2n9Yghr7Zl2pRWyvocdX8POb1y9jei+MOpaGHsThLQ9W6re\n9pLtPLHiCdYXref05NP564S/EhMU0+7r7SUl7L3ySpTRRPqnn2AMDW33tT2JvaSEvZdehgoKJH3e\nPIxhYe2+tqCygId/fpgVeSsYnzSeB09+kAFRA9p/8/J8yHwH1r4NZQfAHAwDzoWhl0L/c1yPu1pQ\nVmVj+Z5ilu4qZvkeK3uKXKvqDe7/10alRTIkKZyhSWEMSgwnVIKKx0gAcZMA0jFaawrLa9lbXMm+\n4kr2Fleyu6iS7fmHyT50ZJe76BALQ5PCGdM7igz3J8NWA0ZtBWSvhgMrYPci1/do16rmodNgxHRI\nyWgx9Xprsg5l8Z9N/+Gbfd8QGRDJ/WPv5+J+F2NQ7X/U4aio4MCNN1G7axe9332HoB68cLA9qtau\nZf8NNxI8ahSpr76CIbj9M6Cc2skHOz7ghcwXqLRVckm/S7hx2I30jezb/gY4Ha7ZWls+g23zobII\nlNH1b993MvSd5FqkaG550L60qo7MA6WsO3CIdQcOsTG7jPIae8P51OgghiSG0y8+lPSYEPrEhtAn\nNpi40IAevaanO0gAcZMAcrQ6u5OCwzXkldWQV1ZNXlkN+WU15JZWc/BQNfutlVTVHdm61GI0kBYT\nzODEMPenvXCG9gonPqyVX8q6Kteji4JNkL/Zlegwb4PrEQbKtU/HwCkw8DxIGt3mjKrGauw1LDqw\niE+zPmVF3gqCTEFMHzydW066pX2DvI3Y8vI4eMed1O7cScoLLxB21uQOXd9TlX31Fbm/fZCgESNI\nnvNPzPHxHbq+tKaUlze+zLyd86h11HJ68ulc1v8yzkg5o+3ZWo05He4PFAth92LIzQS0a01Jwkmu\noJKcAQnDIHYAmJsP/mutyS2rYXueu1ecX872vMPst1Zhdx75mxYaYKJ3TDCpUcEkRgTSKzKQpIgg\nekUGkhgRREJYACYZYzmK3wYQpdQU4J+AEXhNa/10k/MBwNvAWMAKXK213uc+9xBwC+AA7tFat5JC\n9IjjOYBoramsc1BRY6e8xoa1so6SyjrX14o6SiprG46VVNZRXFHX4pakYYEmkiICSY4Mok9sCOmx\nIfSJcX3tFRl0JPdUPVs1lGVDyV44tBdK9ri+L9nt+l67k+kFhEPiCEgbD2mnQurJrv0lOvD+DpQf\nYHX+av6X/T9W5q2k2l5Ncmgyl/a/lOmDphMZ2Poq9BbrdDgo/eQTCp95FpxOkp9/ntCJp3eojp7u\n8DffkvvQQxiCgoh/8LdETJuG6kAgByipKeGDHR/w4Y4PKa4uJsQcwsTkiYxPGs/4XuNJDj32Op1m\nqkpcSRpz1kLOGsjJhIbkjsqVRiVuEMQOdH0fmeZ6RaQ22/fF7nCSU1rd0IveZ61ib3ElOaXV5JVW\nU1l39N7uBgXRIQHEhlqIDrEQExpATIjF9QoNICbUQlSwhfAgE+GBZsKDzIRYjMd1r8YvA4hSygjs\nBM4FsoHVwDVa662NytwBjNBa36aUmg5cprW+Wik1FHgPGAf0An4ABmqtHU3v05g/BBCtNbV2J9V1\nDmrsDqrrHFTbHNTYHFTXOV1fbUeO1R+vstmprLVTXlP/sh31fUWtHecx/unCA03EhAYQHeL6xYgN\ntZAYHkRSRGDDJ7HE8EBCqYaaUqguhZqyRt+XQkWBa3pmRf6RrzVlR9/IEgpR6a5khvFDIXE4JJ4E\nkb3b/ViqtKaUfYf3caD8APvK9rG1ZCubizdTVuu6V1JIEmeknME5vc9hXOK4Dj2qArDl5nL42+84\n9N572A4cIDgjg6S/Pt4j07V7Qs3OneT/6c9Ur1+PpW9foq69lrBzz+nw6nuH08HqgtV8vfdrlmYv\npajalUQxNiiWIdFDGBIzhAGRA0gJSyElNIWIgIj2/eF1OqB4FxRtc83OK9ru+lq8C5y2o8sGx7hy\noYXENXrFur4Gx7pS3QSEQ2A42hJKOUHkldWRW1ZNXqmrJ15cUYu1wvXhy1rh+uDV+PFYUwYF4UFm\nwgLdQSXQTHiQiRCLiSCLkWCLkSCLieD6781Ggi0mgiwGgsxHjgeajVhMBixGg+uryYDJoLo9OPlr\nADkVeFRrfb7754cAtNZPNSrzrbvMcqWUCcgH4oDfNy7buNyx7tnZAPLkgm1U1zmwOZzUOZzYHBqb\n3dnoZyc2uxOnw47N4cThcGB32HHYHdgdThwOO3an67jD4cCAxoBGoTHgbPheKe0+52z4asaBCQeB\nRgfhZgizQLhZE2bRhJogxKwJMblewUYnwUYnQUbtOm50EGyw/X97dxcjV1nHcfz7O+fMdna3u20p\nTaPBnqIAAAnASURBVFsolJeWXmArRVJjSIgREAgVjXJREg0xGrwAAzHGqBcavfLKeGdCeAkqQhAk\naZSIRCRqCK+VF0sxtk2xS7eWttDtzm67c878vThnd2Z3Z/blNLNndvl/kpPzPuc/2+n5P89znpmH\nMmOE8WhaW6iOZNNow3wUxobTXlAz5eBwGfSthb71sHwt9K0j6V1D3H8B8YqLiFduoFpeQWwxca0+\nnU3OUqlW0imuMFIdYaQ6QiWucOrsKY6PHufE6AlOnDnB8dHjVKr1nx0PFHDZisvYev5Wtq7ZyvY1\n27l85eUt/1NZrUZtZJRapUJtpELy4YdU3z9CdXCQsQP7GX3zLcYOHQKgvG0bq7/5DfpuuGHepe6l\nxmo1hv74DCcffpgz76RluK5Nl9N95ZUsu+IKSuvXE61bT3TeKoLeXoLly1G53PrfwYyDpw7y8uDL\n7D2xl30n93Hwo4MkDZ+v5aXlrO1Zy6ryKs4rn8eq8ipWl1fT19VHT6mH7qibnqiHnlIPPVEP5ahM\nKSgRBVE6R0QjJymdPko0NEj40X/R0EBauKkcS5+rDH8A8WjTGFNKBx5b1p/Ou3rTJrKonD5/ibqh\nVCYOy5yxEiO1EpVaiRErMRoHVBIxUhXDsRiuikoVTldhaAyGq2I0EcNVqMQBMSFVQhJCaoja+B3A\nlK2nd4KEYGI/CojCkCgKKYUhURhRKkXpchTRVQonEk4pFFEQEIaiFIgorG9b0V3iuzdtyfXZmE8C\nWciuCxcChxvWB4CpY4ZOHGNmsaRTwOps+0tTzp1nPXnudvzsy0TTCiBpotV4vtWkzRM0aZPqx7fQ\nuN8aXnYuxwOcyaYPW+yfdv5EfN1ANzJN7Gh+bgIMgB2evqvJe2/UYzD1Ua0QgYQUECACAgJFBAoJ\nlK6jQWAQ7FkS0mprer3JF7Q4xkZb3yyiNWsob9vGytu/Qt+NN9K1cWPLYz9uFASs+MJO+nfeytiB\nA5z+y/OM7tnD8Isvtv4V3zBEXV0oDFEYQhRl8xCFEQRiB2kzAYCxdlLBomoxSe0wib1HzRJiS6hZ\nffwQE4yRTh81ubw1+c+h7FM3PkcNvfBs+mdS4zswYAgYqh/TtDA9fVs5m+beX/DcNbsvZHekprei\ns2Xgpn1tjmphE0jrv8Hsx8zl3PQFpLuAuwAuvvji+cQ34cy6XhTXJjfBaGp4qq82HJfeoNVweMNx\nqGHW+NpTlqVsvxrWqW9vmAthGq/2jh8/9R3N9D7q1xdgU0qY9UsHBA2TCNJEEARZEghAIghCQoWU\ngogoKz1GQSkrSYZoavPT1BLtrOsNi0GYlo57etJ5bw/hypWULriA0rp1BL0zjzzo0l/RXbZpE8s2\n1X/qPhkaojp4lPjoIMmpUyTDw2kNb7iCVauQxFicYEkMSZIux/HkG/DUm/G0m3O6XjMjrsUktZik\nlpBYQmwxSRKTWELNDLMaZkbNahjZ3Ayjlu4fvxXY+NLEFmoT152yzxqPsuY3k0lxW7Zs07enLzax\nz6adM/191xdn2D/pOs2PadWCZN0L89MwC5lABoDG8UY3AEdaHDOQNWGtAE7O8VwAzOx+4H5Im7Dy\nBLrziaX54N25uQj7+wn7+2HLFUWH4jrcQjYEvwpslnSppC5gFzC1rrwbuDNbvh143tIUuxvYJWmZ\npEuBzcArCxS3c865JhasBpI907gHeJa0G+9DZrZX0k+B18xsN/Ag8GtJ+0lrHruyc/dKegJ4B4iB\nu2frgeWcc669/IuEzjnnJsynF9bHuy+jc8653DyBOOecy8UTiHPOuVw8gTjnnMvFE4hzzrlclnQv\nLEkfAO/N87TzgeNtCKddPN728njby+NtrzzxbjSzNXM5cEknkDwkvTbXLmydwONtL4+3vTze9mp3\nvN6E5ZxzLhdPIM4553LxBDLd/UUHME8eb3t5vO3l8bZXW+P1ZyDOOedy8RqIc865XDyBNJB0s6R/\nS9ov6ftFxzMTSQ9JOibpX0XHMheSLpL0V0n7JO2VdG/RMc1EUlnSK5LezOL9SdExzUZSKOmfkv5Q\ndCyzkXRI0tuS3pDU8b94KmmlpCclvZt9hj9TdEytSNqS/V3HpyFJ97XlWt6ElZIUko6eeiPpAFav\nAneY2TuFBtaCpOuAYeBXZvaJouOZjaT1wHoz2yOpD3gd+FIH/30F9JrZsKQS8A/gXjN7aZZTCyPp\nO8A1QL+Z7Sw6nplIOgRcY2aL4jsVkh4B/m5mD2TjGfWYWbORdztKdl97H/i0mc33O3Gz8hpI3Q5g\nv5kdNLMx4HHgiwXH1JKZ/Y10zJRFwcwGzWxPtnwa2Ecbx7U/V5YazlZL2dSxpS1JG4BbgQeKjmWp\nkdQPXEc6XhFmNrYYkkfmeuBAO5IHeAJpdCFwuGF9gA6+wS1mki4BtgMvFxvJzLImoTeAY8BzZtbJ\n8f4C+B5QKzqQOTLgz5Jel3RX0cHM4jLgA+DhrInwAUm9RQc1R7uAx9r14p5A6tRkW8eWOBcrScuB\np4D7zGyo6HhmYmaJmV0FbAB2SOrIpkJJO4FjZvZ60bHMw7VmdjVwC3B31iTbqSLgauCXZrYdqAAd\n/YwUIGtquw34Xbuu4QmkbgC4qGF9A3CkoFiWpOxZwlPAo2b2+6LjmausueIF4OaCQ2nlWuC27LnC\n48DnJP2m2JBmZmZHsvkx4GnSJuRONQAMNNRAnyRNKJ3uFmCPmf2vXRfwBFL3KrBZ0qVZ5t4F7C44\npiUjeyj9ILDPzH5edDyzkbRG0spsuRu4AXi32KiaM7MfmNkGM7uE9HP7vJl9teCwWpLUm3WkIGsK\n+jzQsb0JzewocFjSlmzT9UBHdv6Y4g7a2HwFadXMAWYWS7oHeBYIgYfMbG/BYbUk6THgs8D5kgaA\nH5vZg8VGNaNrga8Bb2fPFQB+aGbPFBjTTNYDj2S9WALgCTPr+O6xi8Ra4Om0TEEE/NbM/lRsSLP6\nNvBoVrg8CHy94HhmJKmHtEfpt9p6He/G65xzLg9vwnLOOZeLJxDnnHO5eAJxzjmXiycQ55xzuXgC\ncc45l4snEOecc7l4AnHOOZeLJxDnFpCkF8a/0Sxp9WIZz8W5ZjyBOLewNgH/yZa3AW8XGItz58QT\niHMLRNJG4H0zG//J9W3AWwWG5Nw58QTi3MK5iskJ41N4AnGLmCcQ5xbOJ4EygKTNpCNeehOWW7Q8\ngTi3cK4CAklvAj8iHdb3zmJDci4//zVe5xaIpP3A9mxMeOcWPa+BOLcAsgGUap483FLiNRDnnHO5\neA3EOedcLp5AnHPO5eIJxDnnXC6eQJxzzuXiCcQ551wunkCcc87l4gnEOedcLp5AnHPO5fJ/S691\njyrXDocAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#make points consistent for all graphs\n", "#use widest possible one\n", "x = np.linspace(-5 * 1.2 / np.sqrt(3) + 3.5, 5 * 1.2 / np.sqrt(3) + 3.5, 1000)\n", "#iterate over sample count and plot\n", "for N in [3, 5, 10, 25]:\n", " se = 1.2 / np.sqrt(N) \n", " plt.plot(x, stats.t.pdf(x, loc=3.5, scale=se,df=N), label='N = {}'.format(N))\n", "plt.legend(loc='best')\n", "plt.xlabel('$\\mu$')\n", "plt.ylabel(r'$P(\\mu)$')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "nbgrader": { "grade": false, "grade_id": "2344g4", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 2.4\n", "Plot the probability of the population mean with a sample mean of 5, a sample standard deviation of 3.2, and 8 samples. Using the `fill_between` command, also color the area corresponding to a 95% double-sided confidence interval" ] }, { "cell_type": "code", "execution_count": 60, "metadata": { "collapsed": false, "nbgrader": { "grade": true, "grade_id": "hb4qw5", "locked": false, "points": 8, "schema_version": 1, "solution": true } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEMCAYAAAA8vjqRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl83PV95/HXZ0bH6MKXZMvyjTE2wpYv4XAkJAQCJqGQ\n7SYFkrSkm12WR0KbNsk2dJMlfZBtH0nTdpNNyUETcnQDhJCjTjBxOAzG2IBPfMuSfEo+dNiWD8k6\nP/vHjEAI2dbIM/PTaN7Px2MenvnN76f5jHW85/f7XubuiIiIxCMUdAEiIpJ+FB4iIhI3hYeIiMRN\n4SEiInFTeIiISNwUHiIiEjeFh4iIxE3hISIicVN4iIhI3LKCLiBZiouLffr06UGXISKSVjZs2NDk\n7iUX2m/Ehsf06dNZv3590GWIiKQVM9s/mP102UpEROKm8BARkbgpPEREJG4KDxERiZvCQ0RE4qbw\nEBGRuKU0PMxsqZlVmVmNmT0wwPP3mdlWM9tsZqvNrDy2fbqZtcW2bzaz76WybhERebuUjfMwszDw\nMPABoA5YZ2bL3H1Hn90ec/fvxfa/HfgXYGnsuVp3X5CqekUSrbWji9XVTVQ3nCaSHWbh1NEsmDya\nUMiCLk0kbqkcJLgEqHH3PQBm9gRwB/BmeLj7yT77FwBaYF3SXk+P86M1+/j289WcaOt823Mzi/P5\n6ofnce1lxQFVJzI0qQyPScDBPo/rgHf138nMPgN8DsgB3t/nqRlmtgk4CXzZ3V9OYq0iCdHa0cWn\nf7aRF6saWVyWx3++biyzxuXQ3u1sqG/j8a0tfPwHr/GFD1zGp99/OWY6C5H0kMrwGOi34h1nFu7+\nMPCwmX0M+DJwD3AYmOruzWa2GPiNmV3Z70wFM7sXuBdg6tSpia5fJC5nO7v55I/WsX7fMT69ZCy3\nzip8MxzysuHGmYVcNy2fb61t5hvP1tDR1c1f31IecNUig5PKBvM6YEqfx5OBQ+fZ/wngwwDu3u7u\nzbH7G4Ba4PL+B7j7I+5e6e6VJSUXnNdLJGncnb95agvr9h7j89cV88HLiwY8q4hkhfgf7y7mhhkF\nfGvlXn7+2t4AqhWJXyrDYx0wy8xmmFkOcBewrO8OZjarz8MPAdWx7SWxBnfM7FJgFrAnJVWLDMG/\nv7qfZW8c4k8XjOa90wvOu2/IjM9eM46543N58Lc72XXoRIqqFBm6lIWHu3cB9wMrgJ3Ak+6+3cwe\nivWsArjfzLab2Wai7R73xLZfD2wxszeAp4D73P1YqmoXiceextP8w/KdLC7L46NXXjKoY7JCxv94\ndzG5YeOvHt9IV3dPkqsUuTjmPjI7NFVWVrqmZJdUc3f+9Ievs/nAMR6+bSLj8uNrVnxp3xm+sbqJ\nL9w4g/s/oPYPST0z2+DulRfaTyPMRRLohV0NrK5p4uMVo+MODoDrp+WzuCzCd1bto/FkWxIqFEkM\nhYdIgnR09fDV3+1gyqhsbr28cEhfw8z4r4vHcrbL+edntiW4QpHEUXiIJMjP1x1gX3Mrn1o8hqyL\nGDU+ZVQ2N15awFNvNLC3oSWBFYokjsJDJAE6u3v43kt7KC/JZfHEyEV/vY9VjMaAf3pmxwX3FQmC\nwkMkAf5j8yHqT7TxJ/NGJWSUeElBFrfMKmJF1TEONp288AEiKabwELlI3T3Od16sYebYnIScdfT6\n8Jwiehy+v7IqYV9TJFEUHiIXaeWuBvY0nuEjVybmrKNXaVE2107J59dbGmlpbU/Y1xVJBIWHyEX6\n6av7Kc7P4topeQn/2n9cfglnOp0fr9qd8K8tcjEUHiIXYW/TGVbtbmTprELCSViX4/LiXOYU5/LU\npsP09GjUuQwfCg+Ri/CzV/cTDsHNM88/f9XFWDqrkIMtnby4oz5pryESL4WHyBCd7ezmFxvquG5K\nPmOHMJp8sN49LZ+CbOOxV/cl7TVE4qXwEBmiZ3ccpaWtk1tmFSX1dSJZIW64tJCX9pykoeVMUl9L\nZLAUHiJD9MuNdZQUZDFvQm7SX+uWywrp7IGfr61N+muJDIbCQ2QIjp48y6rdjdwwo4BQCpaOnTEm\nhxljsnl6WwMjdSZsSS8KD5Eh+M2menoc3j8jeQ3l/d0wo4BdTe3srGtO2WuKnIvCQ2QIfrWpnjkl\nuUwelZ2y17x+egEG/OL1fSl7TZFzUXiIxKmm4RRVR07xvgssL5toxflZzC+N8PtdzRrzIYFTeIjE\n6ektRzDgmiSMKL+Q980o4PCpLtbuPpzy1xbpS+EhEqentx7iygmRIa0UeLGunpxH2OB3m+tS/toi\nfaU0PMxsqZlVmVmNmT0wwPP3mdlWM9tsZqvNrLzPc38bO67KzG5JZd0ivaqPnmL30dO8Z1pqL1n1\nKswNs2BihJU1J+ju7g6kBhFIYXiYWRh4GLgVKAfu7hsOMY+5+zx3XwD8I/AvsWPLgbuAK4GlwHdi\nX08kpZ7eejh6yWpy4qZej9e1U/I5crqLjXsaAqtBJJVnHkuAGnff4+4dwBPAHX13cPe+q94UAL0d\n2u8AnnD3dnffC9TEvp5ISi3fepgrJ0SSOh3JhVw9JZ+QwW83HwysBpFUhsckoO9Pe11s29uY2WfM\nrJbomcdfxnOsSDIFfcmq16hImLkTIrxQfUIDBiUwqQyPgYbhvuMn390fdveZwBeBL8dzrJnda2br\nzWx9Y2PjRRUr0t8fdhwF4OoAL1n1um5KPnUnO9myT5euJBipDI86YEqfx5OBQ+fZ/wngw/Ec6+6P\nuHulu1eWlJRcZLkib/fcjqPMGpcbSC+r/q6ZkocBv1WvKwlIKsNjHTDLzGaYWQ7RBvBlfXcws1l9\nHn4IqI7dXwbcZWa5ZjYDmAW8noKaRQBoOt3O5oMnWDI59WM7BjI2P4vZxbm8VHNcl64kECkLD3fv\nAu4HVgA7gSfdfbuZPWRmt8d2u9/MtpvZZuBzwD2xY7cDTwI7gN8Dn3F39VOUlHlhVwMOLJk0PMID\nYMnkPKqb2znQ2BJ0KZKBUnr+7e7LgeX9tj3Y5/5nz3Ps3wN/n7zqRM7t+Z1HKS7I4tIxqZvL6kKW\nTMrjp5tPsGJLHffeNDrociTDaIS5yAWc7ezm5eomlkzKw1Iw/fpgTRudTUlBmJW7m4IuRTKQwkPk\nAl7d00xrR/ewae/oZWYsmZTHxvpWTreeDbocyTAKD5ELeH5nA5Eso2JC8F10+7tqUj7t3c7KHfVB\nlyIZRuEhch7uznM7j7KgNEJOePhcsupVURohN2w8v/No0KVIhlF4iJzHzsOnONxylndNCXZU+bnk\nhI2FEyO8su+kJkqUlFJ4iJzHi7ujI7gXT8wNuJJzq5yUR+OZbrYdUMO5pI7CQ+Q8Vu1u5NIxOYFO\nhHghi8qiDfnP79ACUZI6Cg+RczjT3sWG/cdZWDb8Gsr7Gl+QxeRLsnil9ljQpUgGUXiInMPa2mY6\nu53FE4dXF92BLCrLY+vRs5w80xZ0KZIhFB4i57CqupFIlnFFyfBt7+i1aGIeHd3Oqp3nm2tUJHEU\nHiLn8NLuRuZNiJA9DLvo9jd3Qi7ZIXhpt6Zol9RQeIgM4EBzK/ubW1lcNvwvWQFEskJcOT7C6/tP\naZZdSQmFh8gAXqqOLia2cOLwbizva1FZhP0tnew5ooZzST6Fh8gAVu1uZEJhFmVFw7eLbn+LYg37\nL6jLrqSAwkOkn87uHtbUNrFoYmRYzaJ7IdNGZzM2L8zqmuagS5EMoPAQ6Wfj/uOcae9m8aT8oEuJ\ni5mxaGKEDfWtdHR2BV2OjHAKD5F+XqlpImQwd/zw76Lb38KJeZzu6GF9rSZKlORSeIj0s6a2mVnj\ncinMSb9fj4rSaAP/S1VHAq5ERrr0++0QSaIz7V1sPnhiWK7dMRhj8sJMHZXNuv1a11ySK6XhYWZL\nzazKzGrM7IEBnv+cme0wsy1m9ryZTevzXLeZbY7dlqWybskc6/Ydo6vHmV+anuEBML80wvajZ7W6\noCRVysLDzMLAw8CtQDlwt5mV99ttE1Dp7hXAU8A/9nmuzd0XxG63p6RoyThra5vJCsGckpygSxmy\nBaUR2rudtdW6dCXJk8ozjyVAjbvvcfcO4Angjr47uPtKd2+NPXwVmJzC+kRYU9vEnOJcIlnpe0V3\n7oQIIYOXNVWJJFEqf0MmAQf7PK6LbTuXTwHP9HkcMbP1ZvaqmX04GQVKZmtp7WRb/Unmp8EsuudT\nkBNi1rgc1h04qalKJGlSGR4DjbYa8CfbzD4BVALf6LN5qrtXAh8DvmlmMwc47t5YwKxvbGxMRM2S\nQV7d24wD8yakXxfd/uaXRqhqaqe55UzQpcgIlcrwqAOm9Hk8GXjH/NFmdhPwJeB2d2/v3e7uh2L/\n7gFeBBb2P9bdH3H3SnevLCkpSWz1MuKtrW0mN2zMHpf+4VExIUKPw+oqTVUiyZHK8FgHzDKzGWaW\nA9wFvK3XlJktBL5PNDga+mwfY2a5sfvFwHXAjpRVLhlhTU0TV47PTYsp2C/kipLoFO2ra7SuuSRH\nysLD3buA+4EVwE7gSXffbmYPmVlv76lvAIXAL/p1yb0CWG9mbwArga+5u8JDEqbxVDu7G06nfXtH\nr9ysEFeMj7ChTlO0S3KkdMpQd18OLO+37cE+9286x3FrgHnJrU4y2at7opMJzkvDKUnOZUFphJ9u\nPkFdUwtTSkYHXY6MMOnbH1EkgdbUNlOQE2Lm2PQd39Ff70DHVbvU7iGJp/AQITq+Y+74XMKh9G/v\n6HXZ2Bzys401tZqiXRJP4SEZr/5EG/ubW0dMe0evcMiYOz7Cxvoz9PT0BF2OjDAKD8l4a2OfzCtG\nUHtHr/kTIxw+1UXt4eNBlyIjjMJDMt6a2iZGRcJMHZ0ddCkJ19vu8eIuzXMliaXwkIzm7qypbWbe\nhFxCabTk7GBNG5XN6EiI1/YeC7oUGWEUHpLR9jW3cqTlLPNLR1Z7Ry8zo2JChE31Z+ju7g66HBlB\nFB6S0dbURkdgV4yA+azOpaI0QnNbNzvq1OtKEkfhIRltTW0zxflhyopSOl42pd4a76F1zSVxFB6S\nsXp6nLW1zVSURrAR2N7Rq7Qwi5KCMK+q3UMSSOEhGWt3wymOnekYse0dvcyM+RMivHG4lc7OrqDL\nkRFC4SEZa03NyJvP6lwqSiOcbO/hjf1a50YSQ+EhGWtNbTMTi7IYXzhy2zt6VUyItXtUqd1DEkPh\nIRmpq7uHV/c0v/lHdaQrLshiUlEWr+07EXQpMkIoPCQjbT90ktPtXSwYYfNZnU9FaYQtR1pp7+gM\nuhQZARQekpF6Z5qdmwHtHb3ml0Zo63TW1erSlVw8hYdkpDW1TUwbnc2YvHDQpaTMvNglutW7Gy6w\np8iFKTwk43R09bBu37GMae/oNSoSZvrobF7br3YPuXgKD8k4mw+e4Gxnz5sjrzPJ/NII24+e5czZ\njqBLkTSX0vAws6VmVmVmNWb2wADPf87MdpjZFjN73sym9XnuHjOrjt3uSWXdMrKsqW3CyKz2jl4V\npRE6up3XqjVFu1yclIWHmYWBh4FbgXLgbjMr77fbJqDS3SuAp4B/jB07FvgK8C5gCfAVMxuTqtpl\nZFlT28xl43IozM2c9o5ec8dHCBm8rHYPuUipPPNYAtS4+x537wCeAO7ou4O7r3T31tjDV4HJsfu3\nAM+6+zF3Pw48CyxNUd0ygrR1dLPpwPE3G48zTUFOiMvG5vD6gZagS5E0l8rwmAQc7PO4LrbtXD4F\nPDPEY0UGtH7/MTq7PSPbO3pVlEbY1XCWk61ngy5F0lgqw2OgaUt9wB3NPgFUAt+I51gzu9fM1pvZ\n+sZGzeEj77SmtplwCMpLMq+9o9f80gjdDqur1O4hQ5fK8KgDpvR5PBk41H8nM7sJ+BJwu7u3x3Os\nuz/i7pXuXllSUpKwwmXkWFPTxOziXPKyM7ej4RUluWSFYHW12j1k6FL5G7QOmGVmM8wsB7gLWNZ3\nBzNbCHyfaHD0/cleAdxsZmNiDeU3x7aJDFpLWydb61sy+pIVQCQrxJziXNbtPxl0KZLGUhYe7t4F\n3E/0j/5O4El3325mD5nZ7bHdvgEUAr8ws81mtix27DHgq0QDaB3wUGybyKC9vvcYPc6IX79jMCpK\nI9Qea6fpZOuFdxYZQNxzUZtZAXDW3bvjPdbdlwPL+217sM/9m85z7KPAo/G+pkivNbVN5ISN2eNy\ngi4lcBWlER7b0sLLuw7zn5bMDLocSUMXPPMws5CZfczMnjazBmAXcNjMtpvZN8xsVvLLFLl4a2qb\nKR+fS3Z45C45O1izx+WSGzZeqVHHEhmawVy2WgnMBP4WKHX3Ke4+HngP0bEYX4v1jhIZtppOt1N1\n5FTGt3f0yg4b5eNzWX9A7R4yNIO5bHWTu79jAYBYm8MvgV+aWXbCKxNJoLWxKdjV3vGW+aURfrzp\nBIeOnaJsbFHQ5UiaueCZx0DBMZR9RIK0praZgpwQM8foc06vN5em3aXxHhK/Qfe2MrOtZvYzM/ui\nmd1qZpPN7EvJLE4kUdbWNnHl+FzCIbV39Jo5NoeCbLV7yNDE01X3vcC/AW1Ex2hsAz6UjKJEEqn+\nRBv7mlszasnZwQiHjLkTImw4eAr3ASd7EDmnQYdHbFLCF939/7r7PcBVQHXyShNJjDU1TQBUZOAU\n7BdSURrh0Kku9h7VRIkSn3guW72tS667VwMVCa9IJMHW1jYzOhJm6mi1d/Q3/812j8MBVyLpJp5B\ngo+Y2UygHtgCRIBtZlbg7meSUp3IRXJ31tQ2MW9CLiFTe0d/U0dnMyo3xJo9zXzyfUFXI+lkMIME\nDcDdb3D3qcCdwNNADZAPbDKzqqRWKTJEe5vOcORku9o7ziFkRkVphI11p+np6Qm6HEkjgxokaGZ/\nYWZTAdz9gLsvA74OfIfoQMGHk1ijyJCtiY3vmDdB7R3nUjEhQlNrN7vqNV2cDN5gwmMp0A08bmaH\nY2uM7yHaWH4X8M/u/n+TWaTIUK2pbaKkIIuJhXFP45YxKmKj7ldpfQ+JwwV/o9z9LNEzjO/ERpIX\nA23ufiLZxYlcjO4eZ01NM4sn5mJq7zinsqIsivPDvLrnGPcFXYykjcG0edxjZk1mdgz4AXBawSHp\nYFt9CyfaOllUlh90KcOaxdo9NtWfoasr7smyJUMN5rLV/wI+AMwBDgD/kNSKRBJkdWx8x/xStXdc\nyPwJEVrae9hyQKPNZXAGEx4n3X2Tuze4+/8CliS7KJFEWLW7kZljcxgdCQddyrA3r7fdY9fRgCuR\ndDGY8JhoZvea2XvMrATQSCsZ9s60d7HxwHEWTNQU7IMxviCLiUVZrN17POhSJE0MpgvKV4iOJP84\nMA8oNLPlwBvAFnd/PIn1iQzJa3ub6ex2Fqu9Y9AWlEZ4ce8ZWs92kB/RaotyfoOZkv0Rd7/f3d/r\n7mOBS4F/BU4AH0x2gSJDsWp3dMnZOcX6IzhYi8vyaOtyXqnSVCVyYfHMqguAu9e5+3J3/7q7/2k8\nx5rZUjOrMrMaM3tggOevN7ONZtZlZh/p91y3mW2O3ZbFW7dkltU1TcybECFHS84OWsWECGGDF6vU\n7iEXFnd4DJWZhYmORL8VKAfuNrPyfrsdAD4JPDbAl2hz9wWx2+1JLVbS2qETbdQ0nGZhmaYkiUd+\nTogrSnJ5dV+LpmiXC0pZeBDtpVXj7nvcvQN4Arij7w7uvs/dtwCaZEeGbHV1tIvuQq1XHrdFZXnU\nHuugvllrm8v5pTI8JgEH+zyui20brIiZrTezV83sw4ktTUaSl2uaGJsXZuooTUkSr0Wx3mkv7FC7\nh5xfKsNjoIvP8ZwbT3X3SuBjwDdj08O//QWiXYrXm9n6xkYNdspEPT3O6upGFk6MaEqSIbh0bA6j\nckOs2q3fHzm/VIZHHTClz+PJwKHBHuzuh2L/7gFeBBYOsM8j7l7p7pUlJSUXV62kpe2HTnK8tZOF\n6qI7JCEzFkyMsO7gaU1VIueVyvBYB8wysxlmlkN0Rt5B9ZoyszFmlhu7XwxcB+xIWqWStlZVRz8x\nz9cU7EO2qCyPlvYeNu5tCLoUGcZSFh7u3gXcD6wAdgJPuvt2M3vIzG4HMLOrzKwO+CjwfTPbHjv8\nCmC9mb0BrAS+5u4KD3mHF6sauGxcDmPyNCXJUC2KLZyldg85n5S2KLr7cmB5v20P9rm/jujlrP7H\nrSE6ul3knFpaO9mw/zh/MndU0KWktTF5YS4dk80aTVUi55HKy1YiSbWqupEeh8pJGt9xsRaV5bH9\n6FmOnWoNuhQZphQeMmKs3NXAJblhZo3VlCQXa9HEPLodVm6vD7oUGaYUHjIi9PQ4L1Y1sKgsQjik\nLroX64qSXPKyjZVVajSXgSk8ZETYUt/CsdZOrpqkLrqJkB02Fk3MY+2+U3R3q8uuvJPCQ0aElbsa\nMGDhRHXRTZQlk/Jobutmwx6dfcg7KTxkRFhZ1cCcklwuyVUX3URZPCkPA1ZsG/RYXskgCg9Je42n\n2tlS16JeVgk2OhJmdnEuq2qOa5ZdeQeFh6S9l2LzMF2l8Ei4JZPzqG5up675VNClyDCj8JC099yO\no4zLz2LG6OygSxlxegN5xZa6gCuR4UbhIWntbGc3q6obedfkPM2imwTTR2dTUhDmhSrNsitvp/CQ\ntPZKTROtHd1cO1VddJPBzFgyKY8NdWc409YedDkyjCg8JK39YftRCnJCXFmiLrrJctWkfNq7nRc0\n2lz6UHhI2urucZ7beZTFZXlkh3XJKlkqSiNEsow/7DgSdCkyjCg8JG1tPHCc5jMdXDNFl6ySKSds\nXDUpj5f3tNCpBaIkRuEhaevZHUfJChmLY+tuS/JcOzWfE2d7WL1LAwYlSuEhacndWbHtCPNLI+Tn\n6Mc42RaX5ZEdguVbFB4Spd86SUvVDafZf6yVa9TLKiXys0MsKsvjxZrj9PT0BF2ODAMKD0lLf9ge\nbbxdUqZLVqly7ZR8Glu7ea1GDeei8JA09fTWw1xRksvY/JSupJzRlkzOI2zwu80abS4pDg8zW2pm\nVWZWY2YPDPD89Wa20cy6zOwj/Z67x8yqY7d7Ule1DDe1jafZefgU75leEHQpGaUoN0xFaYSV1bp0\nJSkMDzMLAw8DtwLlwN1mVt5vtwPAJ4HH+h07FvgK8C5gCfAVMxuT7JpleHp6y2EMuHaKJkJMtWun\n5HPoVBdv7NN0JZkulWceS4Aad9/j7h3AE8AdfXdw933uvgXo/7HmFuBZdz/m7seBZ4GlqShahp+n\ntxziyvERinXJKuWunpKPAcs2HQy6FAlYKsNjEtD3J64uti3Zx8oIUn30FFVHT/PuaeplFYQxeWHm\nlUZYsatZl64yXCrDY6D5Iwa7wsygjjWze81svZmtb2zUafVI9LvYJavr1EU3MO+bHr109XrN0aBL\nkQClMjzqgCl9Hk8GBjviaFDHuvsj7l7p7pUlJSVDLlSGJ3fn6a2HmTchwpg8LTcblGun5JMVgl9v\nOBB0KRKgVIbHOmCWmc0wsxzgLmDZII9dAdxsZmNiDeU3x7ZJBtl+6CQ1DafVyypghblhKsvyeHb3\ncc11lcFSFh7u3gXcT/SP/k7gSXffbmYPmdntAGZ2lZnVAR8Fvm9m22PHHgO+SjSA1gEPxbZJBvnV\nxnqyQ8a7p6qXVdDeO6OAY23dvLRD07RnqpR2V3H35cDyftse7HN/HdFLUgMd+yjwaFILlGGrq7uH\n/9hcz1WT8yjK1SWroF01KY+8LOPXm+q4qWJq0OVIADTCXNLCy9VNNJ/p4MZLC4MuRYBIVoirp+Tz\nYs0JzpztCLocCYDCQ9LCrzbVc0kkzCJNvz5s3HhpAWc6nf9Yvy/oUiQACg8Z9k6e7eQP249w/bR8\nrRg4jFSURhhfEOaXm9TukYkUHjLsLd9ymPauHm6YoUtWw0nIjBsvLWRjfSt7G1qCLkdSTOEhw94T\n6w4ydVQ2l4/LDroU6eemmYU48Nia2qBLkRRTeMiwtuPQSTYfPMHSWUWY6ZLVcDOhMIv5pRF+u62R\n7m5NV5JJFB4yrD3++gFywsYNMzQdyXB108wCjpzuYuUOrfORSRQeMmy1dnTx6031XDc1X2M7hrFr\np+RTlBPi39fuC7oUSSGFhwxbv9tymNPtXSydVRR0KXIeuVkhPjCzkJf3nuJg86mgy5EUUXjIsPXY\naweYOjqH8pKcoEuRC/jg5YW4w49XVQddiqSIwkOGpU0HjrP54AluvaxQDeVpoLQom8VlefxqSwPt\nnZosMRMoPGRYevSVfRTkhLhxpmbQTRcfml3E8bZufr1uT9ClSAooPGTYOXSijeVbD3PzZYXkZ+tH\nNF0sLotQWpjFT1/VOh+ZQL+ZMuz8dO1+3J0/mq2G8nQSMuP2OUXsaDjLK1WHgy5HkkzhIcNKa0cX\nj79+gGun5jO+IKUrBkgC3HxZIUU5Ib67Ug3nI53CQ4aVn687SEtbJ3fMuSToUmQIIlkhbptdxOp9\np9hRr/XaRjKFhwwbZzu7+d5LtcydEOGKktygy5Ehum12Eblh4+HnqoIuRZJI4SHDxi821HH0ZDt3\nzRsVdClyEUZFwtw0s5AVVceoO3Y66HIkSRQeMix0dPXw3RdruKIkl/kTdNaR7v64/BLc4V+e2R50\nKZIkKQ0PM1tqZlVmVmNmDwzwfK6Z/Tz2/GtmNj22fbqZtZnZ5tjte6msW5Lv15vqOHTiLHfOG6VB\ngSPAhMIsbr6skGXbm9jXqClLRqKUhYeZhYGHgVuBcuBuMyvvt9ungOPufhnwf4Cv93mu1t0XxG73\npaRoSYn2rm6+/UINs8blsljLzI4Yd84dhQH/9My2oEuRJEjlmccSoMbd97h7B/AEcEe/fe4AfhK7\n/xRwo+lj6Ij372v3U3e8jXsWjNZZxwhSXJDFBy8vYvnOY+w+fCLociTBUhkek4CDfR7XxbYNuI+7\ndwEtwLjYczPMbJOZvWRm70l2sZIaLa2dfPuFGhaX5bFAZx0jzkfnjiI7ZHztaZ19jDSpDI+BPlL6\nIPc5DEx194XA54DHzOwdAwHM7F4zW29m6xsbGy+6YEm+77xUw8m2Tj65cEzQpUgSjI6E+c9XXsIL\nNS2s1qhkbAyCAAANtklEQVTzESWV4VEHTOnzeDJw6Fz7mFkWMAo45u7t7t4M4O4bgFrg8v4v4O6P\nuHulu1eWlJQk4S1IIh081sqPXtnH+y8tYMYYrU8+Uv1x+SWU5If5u2Xb6e7p/3lR0lUqw2MdMMvM\nZphZDnAXsKzfPsuAe2L3PwK84O5uZiWxBnfM7FJgFqCpO9OYu/OVZdsJG3xi/uigy5EkimSF+PNF\nY6hpbucnL2vakpEiZeERa8O4H1gB7ASedPftZvaQmd0e2+2HwDgzqyF6eaq3O+/1wBYze4NoQ/p9\n7q65D9LYiu1HeWFXAx+vGEWJ5rAa8d4zLZ8rx+fyrRdqaT59NuhyJAHMfWSeRlZWVvr69euDLkMG\ncLq9i5v++SXys5z/s3QC4ZB6WGWCPcc7+Ovlh7l5zji+e8/VQZcj52BmG9y98kL7aYS5pNw/raji\n6MmzfGbJGAVHBrl0TA4fnTuKZ3Y28/stdUGXIxdJ4SEp9XJ1Iz9es48PzS5idrGmIck0d84dxbTR\n2XzpN9toae0Iuhy5CAoPSZnjZzr4wi/eYOroHD65UI3kmSg7bPzVNeM43tbNF5/cwEi9bJ4JFB6S\nEu7O//z1VppPd/D5a8cRydKPXqaaNS6Xu+eN4ve7jvH/1qjTZLrSb7CkxA9X7+WZbUf4xPzRzByb\nE3Q5ErA/mTuKBaURvrq8iu31mrokHSk8JOlWVzfxD8t3cu3UfP64XOuSC4RDxuevK6YgJ8R9P13H\nCbV/pB2FhyTVgeZW7n9sI1NH5/BX14wjpIkPJWZMXpgvvqeYw6c6+C+Pvkp7V3fQJUkcFB6SNI2n\n2vmzR1+ju6eHL11fQn62ftzk7eaOj/DZa8axse4Un3tcDejpRL/NkhQtbZ382aOvc6TlLA++r4SJ\nRRpFLgO7YUYhn5g/iqe3N/J3v9miAEkT+o2WhGtp6+TPf/Q61UdP8eAN47miROM55PzunDuKE2d7\n+MlrdWRlZfHl28q1tsswp/CQhGo4dZZ7Hl1HTcMpvv7h2VxR2E5PT0/QZckwZ2b898oxGPDDV/bR\n7fDgbeWENAPBsKXwkITZ23SGT/7odRpOtvPDe65i/oQc6uvrgy5L0oSZcW/lGEqKx/HoK/s4dKKN\nb961gPwc/ZkajtTmIQnx3I6j3P6vqznZ1snP/tu7uP5yraci8TMzvvTBOXzlj8p5budR/uT7a6k/\n0RZ0WTIAhYdclI6uHr7++13815+uZ/q4An77F+9m0VStCihDZ2b8+XUz+Lc/q2RfUyu3fnMVy7dq\nFcLhRuEhQ7al7gR/9O3VfPfFWu66agq/uO8aJo/JD7osGSFuvGICT//lu5lRUsinf7aRzz/5Bs2n\n24MuS2J0MVHi1nS6nW89V81jrx+guDCHRz9ZyfvnTAi6LBmBpo0r4Kn7ruFbz1XzvZdqeW7nUf5m\n6WzurJxCVliffYOk8JBBO9HawU/W7OffXt5DW2c3H1sylS/cMptReVp/XJInOxziC7fM5o4FZXz5\nN9v40q+38cOX9/LZm2ZxW0WZ1oQJiMJDLqjqyCkef/0AT64/SGtHNx8on8AXl87hsvGFQZcmGWTW\nhCKeuPdqVmw/wjefq+azT2zmn/+wmz+7ZhofXTyFUfn6EJNKCg8Z0KETbTy74yi/3FjHlroWskLG\n7fPLuPe9lzKn9JKgy5MMZWYsnTuRm8tLWbH9CD96ZR//++md/OOKKm6YXcJtFWXceMV4de9NgZT+\nD5vZUuBbQBj4gbt/rd/zucBPgcVAM3Cnu++LPfe3wKeAbuAv3X1FCksf8Y6d6WDzweO8vvc4L1Y1\nsOvIKQDKJ17Cg7eVc8eCMsYVaqS4DA+hkHHrvIncOm8iOw6d5Mn1B1m+9TArth8lNyvEVdPHct1l\nxVw7cxxzJhaRmxUOuuQRx1I1j4yZhYHdwAeAOmAdcLe77+izz6eBCne/z8zuAv6Tu99pZuXA48AS\noAx4Drjc3c85DWdlZaWvX78+eW8oTbV1dLO36Qy1jadjtzNsrTvBvuZWIDpVduW0Mbx/znjeP2c8\nsyYMfQr1lpYW6uvrNcJc4lJeXk4oFH9jeHeP8/reY/xhxxHW1DRTdTT6ASg7bMwuLWJu2SjmlBYx\nrbiAaWPzmTwmnxwtSvYOZrbB3SsvtF8qzzyWADXuvgfAzJ4A7gB29NnnDuDvYvefAv7VohPc3AE8\n4e7twF4zq4l9vbUpqj1wPT1OR3dP9NbVQ2fs346uHtq7ejjd3sWps12cOtv55v2TbZ00nm6n4WQ7\nDafOcvRkOy1tnW9+TTMoG5XHlWWXcOdVU1k4dTQVk0fplF/SUjhkXDNzHNfMHAdEp8p5fe8xttWf\nZFt9C89sO8IT6w6+uX/IoPSSCCVFuRQXRm/jCnMYW5BDYW4WBblZb/5bkBumICeL7KwQ2WEjJxwi\n+82bZeQ8XKn8KzEJONjncR3wrnPt4+5dZtYCjIttf7XfsZOSUeSJ1g4+8r21uDsO4ODReog9xB0c\nj/7b58St95i+z/c9ht5jBvp6fV7vre1Oj0Nndw9dPfGfIeaEQ5QU5VJSlMuM4gKuvnQcEy6JMG1c\nPjNLCplRXEAkO7mn85ohVYIyvijCbRVl3FZRBkR/FhtPt3OguZX9za3sbz5D3fE2Gk+3c7jlLFvr\nW2g+00H3EH7XssNGdjhE2AyzaNtMyCDU57Hx1uO3tsfu89Y+72DnfRjd1i+8rph4Cd++e2Hc7yMe\nqQyPgd5z/+/SufYZzLGY2b3AvQBTp06Ntz4g+ulldu+lGuNt39Te78/bvtEGhvV77q1tb31PY4/7\nPd/369Hnud77IYt2VczJin7Kyc16635On+1FkazYLfvN+0Ff583KysLdh3QJQjKTuyftU7yZMb4o\nwviiCJXTxw64T0+Pc6q9izOx2+n2Ls60d3O6vYvWjq7oGX+30xk7++/q8bddCeh2f9uHwZ7YB8A3\nt3l0W+9zvdt7fIA/aLzzw9eAsTbAxilj8uL834lfKsOjDpjS5/Fk4NA59qkzsyxgFHBskMfi7o8A\nj0C0zWMoRRZFsnn444uGcqj0U1BQwNy5c4MuQ2TQQiFjVF62xi4NQio/Eq4DZpnZDDPLAe4ClvXb\nZxlwT+z+R4AXPBq9y4C7zCzXzGYAs4DXU1S3iIj0k7Izj1gbxv3ACqJddR919+1m9hCw3t2XAT8E\n/j3WIH6MaMAQ2+9Joo3rXcBnztfTSkREkitlXXVTTV11RUTiN9iuumrJFBGRuCk8REQkbgoPERGJ\nm8JDRETipvAQEZG4jdjeVmbWCOwPuIxioCngGpJB7yu96H2lj+Hwnqa5e8mFdhqx4TEcmNn6wXR5\nSzd6X+lF7yt9pNN70mUrERGJm8JDRETipvBIrkeCLiBJ9L7Si95X+kib96Q2DxERiZvOPEREJG4K\njyQxs6VmVmVmNWb2QND1JIKZTTGzlWa208y2m9lng64pUcwsbGabzOx3QdeSKGY22syeMrNdse/Z\nNUHXlAhm9texn79tZva4mUWCrmkozOxRM2sws219to01s2fNrDr275ggazwfhUcSmFkYeBi4FSgH\n7jaz8mCrSogu4PPufgVwNfCZEfK+AD4L7Ay6iAT7FvB7d58DzGcEvD8zmwT8JVDp7nOJLu9wV7BV\nDdmPgaX9tj0APO/us4DnY4+HJYVHciwBatx9j7t3AE8AdwRc00Vz98PuvjF2/xTRP0ZJWUs+lcxs\nMvAh4AdB15IoZnYJcD3RNXJw9w53PxFsVQmTBeTFVhvNZ4BVRdOBu68ium5RX3cAP4nd/wnw4ZQW\nFQeFR3JMAg72eVzHCPgj25eZTQcWAq8FW0lCfBP4G6An6EIS6FKgEfhR7HLcD8ysIOiiLpa71wP/\nBBwADgMt7v6HYKtKqAnufhiiH9aA8QHXc04Kj+SwAbaNmG5tZlYI/BL4K3c/GXQ9F8PMbgMa3H1D\n0LUkWBawCPiuuy8EzjCML4EMVqwN4A5gBlAGFJjZJ4KtKjMpPJKjDpjS5/Fk0vTUuj8zyyYaHD9z\n918FXU8CXAfcbmb7iF5efL+Z/b9gS0qIOqDO3XvPDJ8iGibp7iZgr7s3unsn8Cvg2oBrSqSjZjYR\nIPZvQ8D1nJPCIznWAbPMbIaZ5RBt0FsWcE0XzcyM6DX0ne7+L0HXkwju/rfuPtndpxP9Pr3g7mn/\nSdbdjwAHzWx2bNONwI4AS0qUA8DVZpYf+3m8kRHQEaCPZcA9sfv3AP8RYC3nlRV0ASORu3eZ2f3A\nCqK9QR519+0Bl5UI1wF/Cmw1s82xbf/T3ZcHWJOc218AP4t9gNkD/HnA9Vw0d3/NzJ4CNhLt/beJ\nNBqV3ZeZPQ68Dyg2szrgK8DXgCfN7FNEg/KjwVV4fhphLiIicdNlKxERiZvCQ0RE4qbwEBGRuCk8\nREQkbgoPERGJm8JDRETipvAQEZG4KTxEUsjMXuwd9W1m4/qu5SCSThQeIql1GVAdu18BbA2wFpEh\nU3iIpIiZTQPq3b136vcKYEuAJYkMmcJDJHUW8PawWIzCQ9KUwkMkdeYDEQAzm0V0XQpdtpK0pPAQ\nSZ0FQMjM3gAeJDqV+D3nP0RkeNKsuiIpYmY1wMLY+u8iaU1nHiIpYGZFQI+CQ0YKnXmIiEjcdOYh\nIiJxU3iIiEjcFB4iIhI3hYeIiMRN4SEiInFTeIiISNwUHiIiEjeFh4iIxO3/A8TL3nq/1agTAAAA\nAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#get standard error from this info\n", "se = 3.2 / np.sqrt(8)\n", "\n", "#make 500 points from -5 sigma to 5 sigma \n", "x = np.linspace(-5 * se + 5, 5 * se + 5, 500)\n", "plt.plot(x, stats.norm.pdf(x, scale=se, loc=5))\n", "y = stats.t.ppf(0.975, df=8)\n", "x2 = np.linspace(5 - y * se, 5 + y * se, 500)\n", "plt.fill_between(x2,stats.norm.pdf(x2, scale=se, loc=5), color='lightgray')\n", "plt.xlabel('$\\mu$')\n", "#prefix with r to allow curly braces\n", "plt.ylabel(r'$P(\\mu)$')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "nbtrstr", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "Problem 3 Instructions\n", "----\n", "\n", "State what distribution the following data follows and compute or plot the requested property." ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "ntyrss", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 3.1\n", "\n", "The proportion of students graduating at colleges follows a Beta distribution. Each year, a college rankings agency compiles these rates for 100 colleges and also reports the average. If the sample mean and sample standard deviation are 0.74 and 0.21, respectively, what is the population mean for all colleges with 95% confidence" ] }, { "cell_type": "code", "execution_count": 46, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "btrqegrqegr", "locked": false, "points": 3, "schema_version": 1, "solution": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.74 +/- 0.04\n" ] } ], "source": [ "Z = stats.norm.ppf(0.975)\n", "print('{} +/- {:.2f}'.format(0.74, 0.21 * Z / np.sqrt(100)))" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true, "nbgrader": { "grade": false, "grade_id": "45thtrye", "locked": true, "schema_version": 1, "solution": false } }, "source": [ "#### 3.2\n", "The number of fatalities due to work related injuries follows a Poisson distribution at each job site. OSHA computes the average for each state taken from 45 random job sites. All numbers are in fatalities per thousand workers. Alaska has a sample mean of 7.4 a sample standard deviation of 6.0 and Texas has a sample mean of 5.3 with a sample standard deviation of 2.3. Plot the probability for the population mean number of accidents per thousand workers for these two states." ] }, { "cell_type": "code", "execution_count": 47, "metadata": { "collapsed": false, "deletable": true, "editable": true, "nbgrader": { "grade": true, "grade_id": "asebtrsgtr", "locked": false, "points": 8, "schema_version": 1, "solution": true } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VfX9+PHXOztkQSABJIwwZMgSIi5UFLeCbdWvq1pH\ntWqt1q/6rdaKs62tba3+tLXUOqugOLFgrdSBiihDQDaIIGFlMDIgkOS+f3+cc5MQMm6Se+7JeD8f\nj/s49557xvvewH2fzzyiqhhjjDEAUX4HYIwxpvWwpGCMMaaKJQVjjDFVLCkYY4ypYknBGGNMFUsK\nxhhjqlhSMMYYU8WSgjHGmCqWFIwxxlSJ8erAIvIMcC6Qp6rD63j/MuAX7ssS4AZVXdrYcbt166b9\n+vULZ6jGGNPuLVq0qEBVMxrbzrOkADwHPAG8UM/73wInqeouETkLmAoc3dhB+/Xrx8KFC8MWpDHG\ndAQisimU7TxLCqo6V0T6NfD+vBov5wNZXsVijDEmNK2lTeEa4F2/gzDGmI7Oy+qjkIjIyThJYXwD\n21wHXAfQp0+fCEVmjDEdj69JQURGAk8DZ6lqYX3bqepUnDYHcnJybK5vY0yV8vJycnNzKSsr8zuU\nViEhIYGsrCxiY2Obtb9vSUFE+gBvAJer6lq/4jDGtG25ubmkpKTQr18/RMTvcHylqhQWFpKbm0t2\ndnazjuFll9RpwASgm4jkAvcCsQCq+hQwBegK/MX9Q1aoao5X8Rhj2qeysjJLCC4RoWvXruTn5zf7\nGF72Prqkkfd/DPzYq/MbYzoOSwjVWvpdtJbeR6ajqdgPX/4d1linM2NaE0sKxh/v3Q2zb4dpF8Pq\nWX5HY0yLvfnmm4gIq1evBmDjxo0MH37IZA4hSU5ODmdoTWJJwUTerk2w4GnIuRq6DYY594FapzLT\ntk2bNo3x48czffp0v0NpEUsKJvKWvOQsT7gNjvsZFKyFrYv9jcmYFigpKeGzzz7jH//4R51JYePG\njZxwwgmMGTOGMWPGMG+eM6HDtm3bOPHEExk9ejTDhw/nk08+OWi/goICjj32WGbNmkVJSQkTJ05k\nzJgxjBgxgrffftuTz+L74DXTAa2ZDX2Pg7QsGHou/OtWWPk29Brrd2Smjbv/nRWs3FoU1mMOOyyV\neycd0eA2b731FmeeeSaHH3446enpLF68mPT09Kr3MzMzef/990lISGDdunVccsklLFy4kJdffpkz\nzjiDu+++m8rKSvbu3Vu1z44dO5g8eTIPPfQQp512GhUVFbz55pukpqZSUFDAMcccw+TJk8PeyG5J\nwURWSR5s/xomTnFeJ3aBrKPg27n+xmVMC0ybNo2f//znAFx88cVMmzaNn/70p1Xvl5eXc9NNN7Fk\nyRKio6NZu9YZmnXUUUdx9dVXU15ezve+9z1Gjx5dtf3EiRN58sknOemkkwBnDMIvf/lL5s6dS1RU\nFFu2bGHHjh306NEjrJ/FkoKJrO/mO8vsk6rXZZ8Acx+Bsj2QkOZPXKZdaOyK3guFhYV88MEHLF++\nHBGhsrISEeHGG2+s2ubRRx+le/fuLF26lEAgQEJCAgAnnngic+fOZdasWVx++eXccccdXHHFFcTE\nxDB27Fjee++9qqTw0ksvkZ+fz6JFi4iNjaVfv36ejOK2NgUTWVsXQ1Qs9BhRva730aAB2LrEv7iM\naabXXnuNK664gk2bNrFx40Y2b95MdnY2ubm5Vdvs2bOHnj17EhUVxYsvvkhlZSUAmzZtIjMzk2uv\nvZZrrrmGxYudtjUR4ZlnnmH16tU8/PDDVcfIzMwkNjaWDz/8kE2bQpoJu8mspGAia8ti6D4MYuKr\n1/Uc5Sy3L4P+J9W9nzGt1LRp07jzzjsPWnf++efzm9/8pur1jTfeyPnnn8+MGTM4+eSTSUpKAuCj\njz7ikUceITY2luTkZF54ofr2M9HR0UyfPp1JkyaRmprKZZddxqRJk8jJyWH06NEMGTLEk88j2sa6\nAubk5KjdZKeNUoWH+8Lw78Okxw5+749Dod94OP/v/sRm2qxVq1YxdOhQv8NoVer6TkRkUShTCVn1\nkYmcXd/C/j3Qc/Sh7/Uc6ZQUjDG+sqRgIqdgnbPMrOOqrsdIZ7zCgb2HvmeMiRhLCiZyCtwZ0rsd\nfuh7mUOcxuadGyIbkzHmIJYUTOQUrIPEdOiUfuh7XQc5y8J1kY3JGHMQSwomcgrXQ7dBdb/XdYCz\nLFgfuXiMMYewpGAip2BddYmgtrgkSO1lJQVjfGbjFExklO2B0jzoNrD+bboOrG6MNqaNKCwsZOLE\niQBs376d6OhoMjIyAPjyyy+Ji4vzM7wms6RgIiNYLVRXI3NQt0GwbIYznsHupGXaiK5du7JkiTMa\n/7777iM5OZnbb7/d56iaz6qPTGQUukmhawMlhS79nHEMZbsjEpIxXnv++ecZN24co0eP5sYbbyQQ\nCLBhwwYGDRrEzp07qays5LjjjuODDz4AYNKkSYwdO5YjjjiCp59+GoCKigouv/xyRowYwfDhw3n8\n8cc9jdlKCiYydn/nLDv3qX+b4Hu7v3NmTzWmqd6905mFN5x6jICzHm7ybsuXL+fNN99k3rx5xMTE\ncN111zF9+nQuvfRSbrvtNm688UZGjRrFkUceySmnnAI4SSQ9PZ29e/eSk5PD+eefz9q1aykoKODr\nr53PtXu3txdNlhRMZOzeBEmZEJtY/zY1k0JwPiRj2qg5c+awYMECcnKcmSX27dtH7969Abj++uuZ\nMWMGzz77LF999VXVPo8++igzZ84EIDc3l2+++YaBAweyZs0abrnlFs4++2xOP/10T+O2pGAiY89m\n6Ny74W0693WWuzd7H49pn5pxRe8VVeXqq6/mwQcfPOS9kpIStm3bRmVlJSUlJSQlJTFnzhzmzp3L\n/PnzSUxMZPz48ZSVldG1a1eWLVvGu+++y+OPP87rr7/O1KlTPYvb2hRMZOze3HDVEThVRnHJ1VVN\nxrRhp556Kq+++ioFBQWA00vpu++cf9t33HEHV155JVOmTOEnP/kJ4EyNnZ6eTmJiIitWrGDBggUA\n5Ofno6pceOGF3H///VXTa3vFSgrGe4EA7MmFIec0vJ2IkzgsKZh2YMSIEdx7772ceuqpBAIBYmNj\neeqpp1i3bh1Lly7liSeeIDo6mtdff50XX3yRCy+8kKlTpzJq1CiGDBnC0UcfDcDmzZu55pprUFVE\nhN/97neexm1TZxvvFW+HPw6Gs/8A465teNuXL4I9W+CGTyMTm2nzbOrsQ7XKqbNF5BkRyROR5fW8\nLyLyuIisF5FlIjLGq1iMz4JtBI1VHwW32WMlBWP84mWbwnPAmQ28fxYwyH1cB/zVw1iMn4I/8mmN\nNDQHtynbA/tsrIIxfvAsKajqXGBnA5ucB7ygjvlAZxHp6VU8xkdVJYUQkkKwNLHHeiCZ0LW1anAv\ntfS78LP3US+g5v/8XHedaW/2bIaEzhCf0vi2qe4/gaJt3sZk2o2EhAQKCwstMeAkhMLCQhISEpp9\nDD97H9U1uU2df1URuQ6niok+fUKolzatS9G26h/7xqS6hcXird7FY9qVrKwscnNzyc/P9zuUViEh\nIYGsrKxm7+9nUsgFatYnZAF1/hKo6lRgKji9j7wPzYRV8VZI6RHatsndAbGSgglZbGws2dnZfofR\nbvhZfTQTuMLthXQMsEdV7ZegPSreXl0CaEx0LCRlQLH9UzDGD56VFERkGjAB6CYiucC9QCyAqj4F\nzAbOBtYDe4GrvIrF+KiyAkp2QEoT+hCk9rSkYIxPPEsKqnpJI+8r8FOvzm9aidJ80EDTkkJKT2cA\nmzEm4mzuI+OtYINxU5OCNTQb4wtLCsZbxdudZahtCgCph8HeQqjY701Mxph6WVIw3ipqZkkBqhOK\nMSZiLCkYbxVvA4l2ehSFqiopWGOzMZFmScF4q3i7M/YgKjr0fYJVTUXWrmBMpFlSMN4q2tq09gSw\nkoIxPrKkYLxVvL1p7Qng3IEtJsGSgjE+sKRgvFW8telJQcSZFsOmujAm4iwpGO8c2OvcGyHUeY9q\nSjnMeh8Z4wNLCsY7Je6PenOSQnIGlOaFNx5jTKMsKRjvlLhTGSd3b/q+yd2dOZOMMRFlScF4J3il\n35QxCkFJmU7Vk41qNiaiLCkY75S4SSE5s+n7BvcpsSokYyLJkoLxTqlbfdSckkIwKVi7gjERZUnB\neKckzxlzEB3b9H2tpGCMLywpGO+U5jltA82RZEnBGD9YUjDeKclvXnsCWEnBGJ9YUjDeKc1rXnsC\nQEw8JKRZm4IxEWZJwXinJSUFsLEKxvjAkoLxRvk+OFDc/JICuEkhP3wxGWMaZUnBeKMlYxSCkjKs\npGBMhFlSMN6oGqPQwuqjUispGBNJlhSMN6pKCi2pPsqA/UXObKvGmIiwpGC8UTXvUQtLCjWPZYzx\nnCUF442SFkxxEVQ1gM2qkIyJFE+TgoicKSJrRGS9iNxZx/t9RORDEflKRJaJyNlexmMiqDQP4tMg\nNqH5x6gawGaNzcZEimdJQUSigSeBs4BhwCUiMqzWZr8CXlXVI4GLgb94FY+JsJK8lrUngE2KZ4wP\nvCwpjAPWq+oGVT0ATAfOq7WNAqnu8zRgq4fxmEgqzW9ZewJUVz3ZVBfGRIyXSaEXsLnG61x3XU33\nAT8UkVxgNvAzD+MxkRSOkkJ0LHTqaknBmAjyMilIHeu01utLgOdUNQs4G3hRRA6JSUSuE5GFIrIw\nP98aHduEkhbMkFpTUqa1KRgTQV4mhVygd43XWRxaPXQN8CqAqn4OJADdah9IVaeqao6q5mRktPDq\n03ivvAz272nZaOag5EwbwGZMBHmZFBYAg0QkW0TicBqSZ9ba5jtgIoCIDMVJCvYL0Na15I5rtSVb\nScGYSPIsKahqBXAT8B6wCqeX0QoReUBEJrub3QZcKyJLgWnAlapau4rJtDWlYZj3KCgp08YpGBNB\nMV4eXFVn4zQg11w3pcbzlcDxXsZgfBD8EQ+OSG6J5AwoL4UDpRCX1PLjGWMaZCOaTfhVTXERhuoj\nuy2nMRFlScGEXzimzQ6qGsBmVUjGRIIlBRN+pfkQlwKxiS0/lg1gMyaiLCmY8AvHwLUgm+rCmIiy\npGDCLxxTXARVlRSs+siYSLCkYMIvnCWF6FhI7GIlBWMixJKCCb/SME1xEZSUYW0KxkSIJQUTXpXl\nsG9XeHoeBSXZVBfGREpISUFEXheRc+qarM6Yg4Rzioug5AxLCsZESKg/8n8FLgXWicjDIjLEw5hM\nWxbOMQpBNtWFMRETUlJQ1TmqehkwBtgIvC8i80TkKhGJ9TJA08ZUlRTCmBSSM5xZV8vLwndMY0yd\nQq4OEpGuwJXAj4GvgMdwksT7nkRm2qaqkkIYq4+SbFSzMZES0oR4IvIGMAR4EZikqtvct14RkYVe\nBWfaoKp5j8JZUqgxgK1z74a3Nca0SKizpD7tznhaRUTiVXW/quZ4EJdpq0ryIbYTxCeH75hVk+JZ\nScEYr4VaffRQHes+D2cgpp0ozQtvzyOoroqyAWzGeK7BkoKI9AB6AYkiciTV911OBTp5HJtpi0rz\nw9vzCGz6bGMiqLHqozNwGpezgD/VWF8M/NKjmExbVpIPXfqG95ixCRCfag3NxkRAg0lBVZ8HnheR\n81X19QjFZNqy0jzI8qCZyaa6MCYiGqs++qGq/hPoJyL/W/t9Vf1THbuZjipQCXsLw199BM4xraRg\njOcaqz4K3hQ3jF1JTLu1txA0EN7uqEFJGZC/JvzHNcYcpLHqo7+5y/sjE45p07wYuBaUnAkbPwn/\ncY0xB2ms+ujxht5X1ZvDG45p07yY4iIoKdOZfbWy3LnHgjHGE41VHy2KSBSmffBihtSgqrEK+ZB6\nWPiPb4wBQut9ZExovKw+qjlWwZKCMZ5prProz6r6cxF5B9Da76vqZM8iM21PaR5Ex0FC5/AfO9km\nxTMmEhqrPnrRXf6hOQcXkTNxZlONxpk/6eE6tvkf4D6cpLNUVS9tzrlMK1CS71QdiTS+bVMFq6Rs\nrIIxnmqs+miRu/xYROJwZkpVYI2qHmhoXxGJBp4ETgNygQUiMlNVV9bYZhBwF3C8qu4SEQ9aKE3E\neDHvUVCSzX9kTCSEejvOc4BvgMeBJ4D1InJWI7uNA9ar6gY3gUwHzqu1zbXAk6q6C0BV7X98W+bF\nvEdB8cnO7Ks2U6oxngp1ltQ/Aier6gRVPQk4GXi0kX16AZtrvM5119V0OHC4iHwmIvPd6ibTVgWr\nj7ySlGElBWM8Fur9FPJUdX2N1xuAxv531lWxXLuxOgYYBEzAmXTvExEZrqq7DzqQyHXAdQB9+vQJ\nMWQTUapOScHLpJCcaW0Kxnissd5HP3CfrhCR2cCrOD/sFwILGjl2LlDzNllZwNY6tpmvquXAtyKy\nBidJHHRsVZ0KTAXIyck5pBeUaQX27YJAuXfVR+B0S9210bvjG2MarT6a5D4SgB3ASThX9flAl0b2\nXQAMEpFst5H6YmBmrW3ewqmKQkS64VQnbWhC/Ka18HI0c1CyVR8Z47XGeh9d1dwDq2qFiNwEvIfT\nJfUZVV0hIg8AC1V1pvve6SKyEqgE7lDVwuae0/jIy4FrQUmZzqR7gUqIivbuPMZ0YCG1KYhIAnAN\ncAROqQEAVb26of3c+zrPrrVuSo3nCvyv+zBtWfAK3tOSQqYzC6tX03MbY0LuffQi0APnTmwf47QP\nFHsVlGmDSgucpde9j8Aam43xUKhJYaCq3gOUuvMhnQOM8C4s0+aU5IFEQad0785RNdWFJQVjvBJq\nUih3l7tFZDiQBvTzJCLTNpXmQadu3tb1V02KZwPYjPFKqOMUpopIF+AenB5Eye5zYxwlHo5mDkq2\nqS6M8VpISUFVn3affgz09y4c02Z5Oe9RUHwqRMdbm4IxHgp17qOuIvL/RGSxiCwSkT+LSFevgzNt\nSCRKCiLOOWz6bGM8E2qbwnScaS3OBy4ACoBXvArKtDGRmOIiKCnDSgrGeCjUpJCuqg+q6rfu4yHA\ngzupmDbpQAlU7ItMUkjOtDYFYzwUalL4UEQuFpEo9/E/wCwvAzNtSNVo5ggMKEvKsN5HxniosQnx\ninEmwBOcUcf/dN+KAkqAez2NzrQNJREYzRwUbFMIBCAq1GsaY0yoGpv7KCVSgZg2rGS7s0zp7v25\nkjJBK51ZWZOsr4Mx4RbqOAVEZDJwovvyI1X9lzchmTaneIezTOnp/blqjlWwpGBM2IXaJfVh4BZg\npfu4xV1njFNSiIqBRA+nuAiy+Y+M8VSoJYWzgdGqGgAQkeeBr4A7vQrMtCHFOyC5e2Tq+IPtFjZW\nwRhPNOV/cc0uqGnhDsS0YSXbnaQQCcEeTlZSMMYToZYUfgt8JSIf4vREOhG4y7OoTNtSvAO69I3M\nuRI6O1VVNlbBGE80mhRERIBPgWOAo3CSwi9UdbvHsZm2omQ79B4XmXNFRdlYBWM81GhSUFUVkbdU\ndSyH3mPZdHQVB5w7oaX0iNw5k+xezcZ4JdQ2hfkicpSnkZi2qcTtjhqpNgVwur4WW0HVGC+E2qZw\nMnC9iGwESnGqkFRVR3oVmGkjSiI4RiEopQdsXRy58xnTgYSaFM7yNArTdhVHcDRzUOphTpfUynKI\njo3cecMgv3g/874p4Ju8EgpKDxAXHUVWl0SG90ojp28XYqJt6g7jr8bmPkoArgcGAl8D/1DVikgE\nZtqI4BQXyRFsUwi2XxRvh869I3feZlJVPl6bzz8+/ZZP1hUAECXQpVMcByoCFO93/kt1TYrjgpws\nfjy+Pxkp8X6GbDqwxkoKz+Pcn/kTnNLCMJyRzcY4incAEplps4NSDnPP3fqTwob8Eu6duYJP1hXQ\nPTWeWyYO4tSh3RnSM4VYt1RQULKfBd/uZObSrfx97gZemLeJW08bxNXHZ1vJwURcY0lhmKqOABCR\nfwBfeh+SaVNKtjsJITrkabRarqqksDVy52yG1xflcs/by4mJEu6dNIzLju5LXMyhP/LdkuM5a0RP\nzhrRkw35Jfxm9mp+M3s1/16+nb/+cCzdUxN8iN50VI1dhpQHn1i1kalT8Y7ItieA06YArbYHUiCg\n3P/OCm6bsZThvdJ479YTuer47DoTQm39M5L5+xVjeezi0azeXsw5j3/K4u92RSBqYxyN/SsdJSJF\n7qMYGBl8LiJFjR1cRM4UkTUisl5E6p0nSUQuEBEVkZymfgDjs5LtkW1PAGfivahYKGp9JYUDFQF+\nNv0rnv1sI1cfn820a4+hZ1pik44hIpw3uhdv/fR4kuKj+eHTX/DZ+gKPIjbmYA0mBVWNVtVU95Gi\nqjE1nqc2tK+IRANPUt0WcYmIDKtjuxTgZuCL5n8M4xs/SgpRUa1yrEJlQLn1lSXMWraNX549hCmT\nhhEdJc0+3uHdU5hx/bH07tKJq55dwNy1NorbeM/LVqxxwHpV3aCqB4DpwHl1bPcg8HugzMNYjBcC\nlc7I4kiXFMBpV2hFbQqBgHLn68uY9fU27j57KNedOCAsx81MSeCVnxzDgMxkrv/nIpZu3h2W4xpT\nHy+TQi9gc43Xue66KiJyJNDbbtjTRpXsAA1U1/FHUmrrKin8ec5aZizK5eaJg7j2xP5hPXbnTnE8\nf9VRpCfFcdVzC/i2oDSsxzemJi+TQl3lZq16UyQKeBS4rdEDiVwnIgtFZGF+vhWhW41gnX5qr4a3\n80JKTyjaFvnz1mHWsm08/sF6Lhybxa2nDvLkHJmpCbx4zdEAXPvCQkr2W78P4w0vk0IuULMTeRZQ\ns7yfAgwHPnKnzzgGmFlXY7OqTlXVHFXNyciIYH9407CiLc4yzaekcKAY9hdH/tw1rNi6h9tnLGVM\nn8489P3hOJMKeyO7WxJPXHokG/JLuGPGUlS18Z2MaSIvk8ICYJCIZItIHHAxNWZZVdU9qtpNVfup\naj9gPjBZVRd6GJMJpz1uUvCrpAC+ViEVlZVz/T8XkZYYy1OXjyU+Jtrzcx43oBt3nTWUd5dv56mP\nN3h+PtPxeJYU3HENNwHvAauAV1V1hYg8ICKTvTqviaCiLRCTAIldIn/u1GBS8KcKSVW5563lbNm1\njycuPZLMlMgNMPvxCdmcM7Inf/jPGhvDYMLO0zH0qjpbVQ9X1QGq+mt33RRVPeS+DKo6wUoJbUzR\nVqeU4GGVSb2CJQWf2hXeWLyFt5ds5ZaJh5PTLz2i5xYRfvuDEfRITeDn05dY+4IJK5tYxTRf0RZ/\neh6Br1NdbCwoZcrbyxnXL52bThkY8fMDpCbE8uhFo8ndtZcH3lnhSwymfbKkYJovWFLwQ3wKxKdG\nfFRzIKDc8dpSoqOERy8e3aLBaS01LjudGyYM4NWFufx7eevoiWXaPksKpnkClU59vh89j4LSsqob\nuyPkxfmbWLBxF/ecO4xenZs2fYUXfn7q4Yzolcbdby5nV+kBv8Mx7YAlBdM8JXkQqPCv+gjcpPBd\nxE63eedefvfv1Zx4eAYXjM2K2HkbEhsdxe8vGMmefeU8OGul3+GYdsCSgmmeqoFrPv44pvWGPbkR\nOZWqctcbXyPAb38wwtPxCE01tGcqN0wYwBuLt/DRmjy/wzFtnCUF0zzBgWt+lxT27YL9JZ6fasbC\nXD5dX8CdZw9tFdVGtd10ykAGZCRx95vLrTeSaRFLCqZ5inwcuBbUuY+z9Li0UFiyn4dmrWRcdjqX\njevj6bmaKz4mmt9fMJKte/bxh/fW+B2OacMsKZjm2ZPrDFzrFNk++gdJy6qOxUMPv7uafeWV/Ob7\nw4nysbdRY8b2TeeHR/flhc83snzLHr/DMW2UJQXTPLu/c+r0/axbT3On1vKwsXnRpp3MWJTLNeP7\nMzAzxbPzhMvtZwwmPSmOX721nEDA5kYyTWdJwTTP7u+gS19/Y0jpARLtWUmhojLAPW+toGdaAj/z\naZBaU6UlxnLXWUNZsnk3ry7c3PgOxtRiScE0z+5N1XX6fomKdto0dnvz4/fP+ZtYua2Ie84dRlJ8\njCfn8MIPxvRiXL90fvfv1TZ2wTSZJQXTdGVFTq+fzj6XFAA6e9MtNb94P3/8z1pOGNSNs4b7cGe5\nFhARHvjeERSVVfD791b7HY5pYywpmKbb7dbh+119BO4AtvCXFH777irKKiq5f/IRrWpMQqiG9Ejl\n6uP7MX3BZr6ymVRNE1hSME0XTAp+Vx+B09hctBUqw9c3/8tvd/LG4i1cd2J/+mckh+24kXbLqYeT\nmRLPPW8vp9IanU2ILCmYptu9yVl27udrGIBTfaSV1eMmWqi8MsA9by2nV+dEbjrZm1trRkpyfAz3\nnDuM5VuKeOmLTX6HY9oISwqm6XZtgtgkf8coBKX3d5a7vg3L4V74fBNrdhQzZdIwEuO8v5Oa184Z\n0ZPxA7vxh/fWUFCy3+9wTBtgScE0XbA7amuoa++S7Sx3tvzWlDuKynj0/bVMGJzB6cO6t/h4rYGI\ncP95R7CvvJKH37VGZ9M4Swqm6VpDd9Sg1F4QHQ87W15S+M3sVRyoDHDfpLbZuFyfARnJXHtCf15b\nlMvCjTv9Dse0cpYUTNOoOiWF1tAdFSAqyim1tLCk8Pk3hby9ZCvXn9ifft2SwhRc63HTKQM5LC2B\nX721nIrKgN/hmFbMkoJpmn27YH9R6+iOGpTeH3ZtbPbu5ZUBpry9nKwuidx4ctsYudxUneJimDJp\nGKu3F/PifGt0NvWzpGCapvAbZ5k+wN84auqS7VQfafO6XT732UbW5ZVw36QjSIht+43L9TnjiB6c\neHgGf/rPWvKKy/wOx7RSlhRM0xSud5ZdW9EVdXp/KC917gbXRNv3lPHnOWuZOCSTU9tJ43J9RIT7\nJx/B/ooAD8+2RmdTN0sKpmkK1zuT0LWq6qPm90D69exVVASUeycdEeagWqfsbkn85KT+vPHVFr7Y\nUOh3OKYVsqRgmqZwPXTpB9GxfkdSrZljFeatL+CdpVu5ccJA+nTt5EFgrdONEwbSq3MiU95eQbk1\nOptaLCmYpin8pnVVHYF7X4foJpUUDlQEuOft5fRJ78RPTurvYXCtT2JcNPdOGsaaHcU8P2+j3+GY\nVsbTpCAZIISjAAAX30lEQVQiZ4rIGhFZLyJ31vH+/4rIShFZJiL/FZFWVCdhDhEIwM5WmBRi4pxx\nEwXrQt7lmc++5Zv8Uu6bPKxdNy7X57Rh3TllSCZ/nrOOHUXW6GyqeZYURCQaeBI4CxgGXCIiw2pt\n9hWQo6ojgdeA33sVjwmD4m1Qvhe6tqKeR0EZQ6BgbUib5u7ay+P/Xef+MLbvxuX6iAj3ThrGgcoA\nv561yu9wTCviZUlhHLBeVTeo6gFgOnBezQ1U9UNV3eu+nA9keRiPaalC90q8VSaFwU5JoZHZUlWV\nX721HIB7J9W+RulY+nZN4oaTBjBz6VbmfVPgdzimlfAyKfQCak50n+uuq881wLsexmNaKn+Ns8wY\n4m8cdckYAoHyRhubZy7dykdr8rn99MFkdek4jcv1uWHCAHqnJ/Krt5ZTVl7pdzimFfAyKdQ1eUyd\no4tE5IdADvBIPe9fJyILRWRhfn5+GEM0TbJjBSSmQ3IrrHLJGOws8+vvf7+r9AAPvLOSUb0786Pj\n+kUmrlYuITaaX39vBBvyS3n8v6G3yZj2y8ukkAv0rvE6C9haeyMRORW4G5isqnXO7auqU1U1R1Vz\nMjIyPAnWhCBvFWQOax2zo9bW7XBn2UBS+PXsVezZV87DPxhBdFQr/Aw+OfHwDC4cm8Xf5m5g+ZY9\nfodjfOZlUlgADBKRbBGJAy4GZtbcQESOBP6GkxCaPhzVRI6qmxSG+h1J3eKTIa1PdRVXLZ+uK+C1\nRbn85KT+DO2ZGuHgWr9fnTOM9KQ47nhtmY1d6OA8SwqqWgHcBLwHrAJeVdUVIvKAiEx2N3sESAZm\niMgSEZlZz+GM3/bkwoFi6N6KG2czBjuJq5bisnLufGMZ2d2S+Nkpbftual5J6xTLQ98bzqptRTz1\n0Td+h2N8FOPlwVV1NjC71ropNZ6f6uX5TRjlrXSWma04KfQYDhs+hIr9EBNftfqhf61i6+59zLj+\nuA45JiFUZxzRg3NG9uTxD9Zx+hE9GNwjxe+QjA9sRLMJzY4VzrI19jwK6jkKAhUHlRbmrNzBKws3\nc/1JAxjbt4uPwbUN908+gpSEWG59ZQn7K6w3UkdkScGEZttS58Y6iZ39jqR+PUY6y21LAdhZeoA7\n3/iaIT1SuOVUqzYKRbfkeH53/khWbiviT/8JbTCgaV8sKZjQbF0Mhx3pdxQN65IN8amwfRmqyl1v\nLGPPvgM8etFo4mOs2ihUpw3rzqVH92HqJxuYt94GtXU0lhRM40oLnVtw9hrjdyQNi4qCHiNg21Je\n+HwT763YwR1nDLbeRs3wq3OGkt0tif99dSm79x7wOxwTQZYUTOO2fuUsW3tJAaDnKALbvubhWcuZ\nOCSTH4/vWDOghkunuBgev/hICkv3c/uMZQQCzburnWl7LCmYxgWTQs/R/sYRgr0ZI4mqLOPoTlv5\nw4WjiLJBas02vFcavzx7KHNW7eCvH1s31Y7CkoJp3JZF0HUQJLTuaphAQHlwaRoAD4wpoUtSnM8R\ntX1XHtePyaMO44//WcOn66x9oSOwpGAaFgjAd59Dn6P9jqRRj/13HdPWKiXx3elT8rXf4bQLIsJv\nfzCCARnJ3Dz9K7bs3ud3SMZjlhRMw/JWQNlu6Dve70gaNGvZNh777zouGJtF0qDx8N18Z2oO02JJ\n8TE8dflYDlQEuPb5hZTsb3h6ctO2WVIwDdv4mbPsd7y/cTRg+ZY93DZjCWP7duHX3x+O9D4Girc6\nPaZMWAzISOaJS49kzY5ibnp5MRU2P1K7ZUnBNGzTp86tLjv38TuSOn1bUMqVz35J16R4nvrhWGc8\nQt/j3Dfn+htcOzNhcCYPnjecj9bkc+/MFaiVxNolSwqmfoFKp6TQSquOdhSVcfk/viCg8PzV48hI\ncec76n4EpPSE9XP8DbAduvToPlx/0gBe+uI7HrP7L7RLnk6IZ9q43IWwbycMan3zFu4qPcAV//iS\nXaUHmHbdMQzMTK5+UwQGToRV7zi354y2f+bh9H9nDKawZD9/nrOOuJgobpww0O+QTBhZScHUb81s\niIqBga0rKeQX7+fiqfP5trCUqVfkMDKrjvmYBp4GZXtgy8LIB9jORUUJD58/kvNGH8bv/72Gpz/Z\n4HdIJozsEsrUb8270Pd4SEjzO5Iq2/bs47K/f8G2PWU8e+VRHD+wW90b9p/gJLTVs6DPMZEMsUOI\njhL+eOEoyisDPDRrFaX7K7l54kCkNd6VzzSJlRRM3QrWQcEaGHyW35FUWZ9XzP/87XPyivfzwjXj\n6k8I4MzmOmAiLH/DGWthwi4mOorHLj6S88dk8eictUx5ewWVNh1Gm2dJwdRtycsg0XDE9/2OBHBu\np/n9v8xj34FKXvrx0RzVL73xnUZcAEW5sHm+9wF2ULHRUfzhwpFcf9IAXpy/iRv+ucjGMbRxlhTM\noQKVsHSa05aQ0sPXUFSVZz79lh89+yWHpSXy1k+PZ1TvEO/pMPhsiEmEZa94G2QHJyLcedYQ7p00\njP+uzuO8Jz5lfV6J32GZZrKkYA61fg4Ub4MjL/M1jN17D3DtC4t44F8rOXlwJq/dcCxZXTqFfoD4\nZBj+A1j2Kuzb5V2gBoCrjs/mn9ccze695Xzvyc94e8kWG8vQBllSMIf67DFIOQwO9689Yc7KHZz5\n50/4eG0e95w7jL9fMZaUhNimH+iYG6B8Lyx6LuwxmkMdO6Ar/7p5PId3T+aW6Uv46cuLKSzZ73dY\npgksKZiDbfocNn0Gx98MMZGfZTSvuIyfvryYH7+wkLTEWF67/jiuGZ/d/F4tPUZA9okw/yk4UBre\nYE2deqYl8upPjuX/zhzMnJV5nP7oXGYs3Gz3ZGgjLCmYaoEAzLkPOnWDMT+K6KlL9lfwp/fXMuGR\nj3h/xQ5uO+1w3vnZ+NDbDxoy4ZdQsh3m/b+WH8uEJCbaGdT2zs/G06drJ+54bRnnPfkZX2wo9Ds0\n0wgbp2CqffWi01PnvL9AXBPq7ltgz95y/vnFJp759FsKSw9wzoie3H7GYLK7JYXvJH2PhWHfc6rF\nRv4PpNvd2CJlcI8U3rjhOGYu3crD767moqnzOaZ/OjedPIjjB3a1cQ2tkLS1hqCcnBxduNBGqYZd\n3mp4eiL0HAVXznKmivCIqrJiaxEzFm5mxqJc9h6o5IRB3bjt9MGMDkfJoC57cuGvx0OXfnDNfyAm\n3pvzmHrtO1DJS19s4u+fbGBH0X6G9UzlknG9Oe/IXqQ2p73INImILFLVnEa3s6Rg2LMFnjvHqXO/\n7iNI6xX2U6gq3+SXMmfVDt5cvIU1O4qJjRYmjTqMa0/oz9CeEbir2+pZMP1SGDoZLnjW5kTyyf6K\nSl5ftIUX529i1bYiEmKjOH1YD844ogcnDc4gOd7+Ll5oFUlBRM4EHgOigadV9eFa78cDLwBjgULg\nIlXd2NAxLSmE2ZbF8OqPnC6bl78BvceF5bCqytY9ZXz13S6+2LCTD9fkkbvLuWvXmD6d+cGYLM4d\n2ZPOnSLcmP35X+C9u5wxGD/4O3QKYRCc8YSqsnxLEdMXfMe7y7ezs/QAcdFRHN0/nWP6d2Vcdjoj\ns9Kc6dBNi/meFEQkGlgLnAbkAguAS1R1ZY1tbgRGqur1InIx8H1Vvaih41pSCJPdm2He47DwGUjK\nhEtehsOObPJhVJWisgq+LSjlm7wSvskvYe2OEpbm7ia/2OmKmBgbzfEDuzFhcAYTBmc0bayBFxY+\nA+/+AuJTYcKdMOpiiE/xN6YOrjKgLNy4k/+s3MGn6wpYs6MYgLjoKAZ1T2ZIj1SG9kxhYGYyvdM7\n0atzIgmxliyaojUkhWOB+1T1DPf1XQCq+tsa27znbvO5iMQA24EMbSAoSwpNVLEf9u6EvYVQuB52\nrIANH0LuQjQqmsqRl7J3whQOxKRyoCJAeWWAAxUB9lcEKNlfQXFZBUX7yikuK3eel5WTX7yf7UVl\n5BU5y70HKqtOFxMl9O3aiZFZnRnduzNH9unMkB6pxMW0so5u25bCu3fCd/MgNsmZQK/vsZAxBLpk\nQ1JXJ2lE2Q+PH3aVHmDBxp0s2rSLVduLWbWtqOoiIygzJZ7DOifSLTmOLp3iSE+Oo2tSHJ07xZEc\nH0NiXDRJcTF0iosmMS7aWcZGExMdRUyUEBMlREdJh2nsDjUpeFl51wvYXON1LlD77u9V26hqhYjs\nAboCBeEOZtlHr5M2916gOt+I+1xq5aCq9Sha47WzLkgP3lbrXl/7ee3j13XMut6rPkZ9+x76XgwB\nEjn4P1KlCqvI5v3KC3itYjxb5mfA/C8JVWJsNN1S4uiRmsDQw1I5eUgm3VPj6dc1iQGZyfRJ70Rs\ndCtLAHXpOQqumu3cM2LJS7DhI1gz69DtYpMgNtFJDhLtLqOqX3v+g+Lx8VvpD2IX4HT3AUBnqEhV\nyisD7sN5XrFLqSxUKjVAZUDrvS13JVDsPmoTABGE6q9DDnrzoCeN/kW8/Eq3D7iQYy6717sT4G1S\nqOurqf0nC2UbROQ64DqAPn2ad1vIuKQ0Cjv1d09Q47RS/bN78HtyUHTB9QJoQ/vUOKZS899UPfvU\ns16kev/6tq0+jxzyT1ZFUIR90ansi01jX0wau+KzKEjoSyA6gbiYKC6IjiIuJor4GGcZGx1FnLsu\nLiaKlPgYUhJiSU10likJMW3jBz9UItD7KOcBUJIHhd/Aro1OG8v+Iigrgop9znxQWumM5dDK6tde\n8rwTSNvqZBLjPhLreV+BimCyCCiV7qMiEKh6HkwcAZylKlRq8LlzERhwv/faX3/N6z6t47vTGht5\n9c3GpHT36Mg1zuHhsXOB3jVeZwFb69km160+SgN21j6Qqk4FpoJTfdScYIYcdSoc1bpuFmNameRM\n59H3WL8jMc0gQKz7MM3n5WXfAmCQiGSLSBxwMTCz1jYzgeDQ2QuADxpqTzDGGOMtz0oKbhvBTcB7\nOF1Sn1HVFSLyALBQVWcC/wBeFJH1OCWEi72KxxhjTOM8HSWiqrOB2bXWTanxvAy40MsYjDHGhK4d\ntRoaY4xpKUsKxhhjqlhSMMYYU8WSgjHGmCqWFIwxxlRpc1Nni0g+sKmZu3fDgyk0Wjn7zB2DfeaO\noSWfua+qZjS2UZtLCi0hIgtDmRCqPbHP3DHYZ+4YIvGZrfrIGGNMFUsKxhhjqnS0pDDV7wB8YJ+5\nY7DP3DF4/pk7VJuCMcaYhnW0koIxxpgGdIikICK9ReRDEVklIitE5Ba/Y4oUEYkWka9E5F9+xxIJ\nItJZRF4TkdXu37vd3xxBRG51/10vF5FpIpLgd0zhJiLPiEieiCyvsS5dRN4XkXXusoufMYZbPZ/5\nEfff9jIReVNEOof7vB0iKQAVwG2qOhQ4BvipiAzzOaZIuQVY5XcQEfQY8G9VHQKMop1/dhHpBdwM\n5KjqcJxp6tvjFPTPAWfWWncn8F9VHQT8133dnjzHoZ/5fWC4qo4E1gJ3hfukHSIpqOo2VV3sPi/G\n+aHo5W9U3hORLOAc4Gm/Y4kEEUkFTsS5TweqekBVd/sbVUTEAInu3Qs7cegdDts8VZ3LoXdlPA94\n3n3+PPC9iAblsbo+s6r+R1Ur3Jfzce5oGVYdIinUJCL9gCOBL/yNJCL+DPwfEPA7kAjpD+QDz7pV\nZk+LSJLfQXlJVbcAfwC+A7YBe1T1P/5GFTHdVXUbOBd+QKbP8UTa1cC74T5oh0oKIpIMvA78XFWL\n/I7HSyJyLpCnqov8jiWCYoAxwF9V9UiglPZXpXAQtx79PCAbOAxIEpEf+huV8ZqI3I1TLf5SuI/d\nYZKCiMTiJISXVPUNv+OJgOOBySKyEZgOnCIi//Q3JM/lArmqGiwFvoaTJNqzU4FvVTVfVcuBN4Dj\nfI4pUnaISE8Ad5nnczwRISI/As4FLvPinvYdIimIiODUM69S1T/5HU8kqOpdqpqlqv1wGh4/UNV2\nfQWpqtuBzSIy2F01EVjpY0iR8B1wjIh0cv+dT6SdN67XMBP4kfv8R8DbPsYSESJyJvALYLKq7vXi\nHB0iKeBcNV+Oc7W8xH2c7XdQxhM/A14SkWXAaOA3PsfjKbdU9BqwGPga5/90uxvpKyLTgM+BwSKS\nKyLXAA8Dp4nIOuA093W7Uc9nfgJIAd53f8eeCvt5bUSzMcaYoI5SUjDGGBMCSwrGGGOqWFIwxhhT\nxZKCMcaYKpYUjDHGVLGkYA4iIpU1uu0ucacFqW/bfiJyaQjH7Bec6VFEckTkcff5BBE5rsZ214vI\nFS3/FOEhInfX+B5qfi83i8hzInKB3zHWR0TuE5Hba63rLCKF7ngGRORYEVF3jixEJE1EdopIyL8L\nrf17ME0X43cAptXZp6qjQ9y2H3Ap8HKoB1fVhcBC9+UEoASY574X9j7XTSUi0apa6cbza+DX7vqS\nmt+LiDznT4TNp6q7RWQ7MBRnUN9xwFfu8lWcGYS/UNWQ5spyJ+BrNhGJqTG5m2klrKRgGuVe6X8i\nIovdR/Dq/mHgBPfq+dYGtqt5rAki8i+3BHI9cKu7/wk1r25FZICI/FtEFrnHHOKuv9C9b8BSEZlb\nz/HnunPNrxSRp4JXviJyuoh87sY2w50LCxHZKCJTRORT4MImfDUnisg8EdkQvFoWxyNujF+LyEU1\nP3eNOJ8QkSvd5w+7sS4TkT+46yaJyBfiTOw3R0S6u+vvE2ee/Y/c895c45h3i8gaEZkDDKZun1E9\nDcZxwKO1Xs9zjzVaROZL9bz9Xdz1H4nIb0TkY5xp2Wt+9w+6JYcoERkrIh+7f7/3pHo6ioP2b+zv\naXygqvawR9UDqASWuI833XWdgAT3+SBgoft8AvCvGvvWt10/YHntfYD7gNtr7F/1Gmd+/EHu86Nx\npukAZ9RuL/d55zrinwCU4cyYGo0z//wFQDdgLpDkbvcLYIr7fCPwf418LyW1Xj8HzMC5sBoGrHfX\nn++eMxrojjMNRc86vqsngCuBdGAN1QNJO7vLLjXW/Rj4Y43vaB4Q736mQiAWGOt+N52AVGB9ze+2\nxnmvBJ5xn38FJACfuq/fB05xny8DTnKfPwD82X3+EfCXWt/DBcDvgb8B4sYzD8hwt7moxjlr79/g\n39MekX9Y9ZGpra7qo1jgCREZjZM0Dq9n31C3a5B7BX8cMMOt/gbnRxCcK93nRORVnMnf6vKlqm5w\njzUNGI+TKIYBn7nHjMOZQiDolWaE+pY6VS0rg1fy7rmmqVMFtcO9Ij4KqG9W3iI3tqdFZBYQLE1k\nAa+4V9hxwLc19pmlqvuB/SKSh5N8TsBJ4nvdzz2znvN9BtwpItnARlUtc0s3yTiJ5UsRScP5gf7Y\n3ed5nAQYVPu7ugen2uk699yDgeE4UzGAkyC31bN/KH9PE0GWFEwobgV24NzJLArnR6wl2zUmCthd\nR3JCVa8XkaNxbh60RERGq2ph7c3qeC3A+6p6ST3nLG1GnPtrPJday9oqOLi6NgFAVStEZBzORHYX\nAzcBpwD/D/iTqs4UkQk4JYS6zltJ9f/jRuesUdV1blXQJKqT4iLgKpzZVkvcpNCQ2t/VAmCsiKSr\n6k6c72CFqtZ3K9Sq/UP8e5oIsjYFE4o0YJt7VXw5zpUfQDHO5FyNbVef2vsDoM69Lr4VkQuhqp5+\nlPt8gKp+oapTgAKgdx3HHSci2W5bwkXApzh3qTpeRAa6x+kkIs0qyTRiLnCROPfGzsC5E9yXwCZg\nmIjEuz+6E904koE0VZ0N/BxnEj9wvsst7vMf0bi5wPdFJFFEUnB+9OvzOU57wOc1Xv+c6gb/PcAu\nETnBff9y4OPaB6nh3zjtS7Pcc68BMsS9P7aIxIrIEXXtGOLf00SQlRRMKP4CvO7+SH9I9ZXeMqBC\nRJbi1C3Xt1193gFeE5HzcGY3reky4K8i8iucaqnpwFLgEREZhHM1+l93XW2f4/xIjcD5sXxTVQNu\nw+40EQlWRf0K5z634fQmcKwbl+K0VWwHcKtIlgHrcOrzwUmKb4tIgvuZbnXX34dTfbYFJ6FlN3RS\nVV0sIq/gtAVtAj5pYPPPgLOp7gX2OU4bzLwa2/wIeEpEOgEbcEoSDZ1/hpsQZrrHvgB43E2AMTh3\nAVxRx66h/D1NBNksqaZdcatablfVc/2OxZi2yKqPjDHGVLGSgjHGmCpWUjDGGFPFkoIxxpgqlhSM\nMcZUsaRgjDGmiiUFY4wxVSwpGGOMqfL/ARWLTbJkWpjRAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#get standard error from this info\n", "se1 = 6.0 / np.sqrt(45)\n", "se2 = 2.3 / np.sqrt(45)\n", "#some trial and error to get this number\n", "x = np.linspace(2, 12, 1000)\n", "plt.plot(x, stats.norm.pdf(x, scale=se1, loc=7.4), label='Alaska')\n", "plt.plot(x, stats.norm.pdf(x, scale=se2, loc=5.3), label='Texas')\n", "plt.legend(loc='best')\n", "plt.xlabel('Fatalities per Thousand Workers')\n", "#prefix with r to allow curly braces\n", "plt.ylabel('Probability')\n", "plt.show()" ] } ], "metadata": { "celltoolbar": "Create Assignment", "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.5.3" } }, "nbformat": 4, "nbformat_minor": 2 }