{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "from datascience import *\n", "%matplotlib inline\n", "path_data = 'https://raw.githubusercontent.com/ChemeketaCS/datasci-textbook/main/assets/data/'\n", "import matplotlib.pyplot as plots\n", "plots.style.use('fivethirtyeight')\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Confidence Intervals\n", "We have developed a method for estimating a parameter by using random sampling and the bootstrap. Our method produces an interval of estimates, to account for chance variability in the random sample. By providing an interval of estimates instead of just one estimate, we give ourselves some wiggle room.\n", "\n", "In the previous example we saw that our process of estimation produced a good interval about 95% of the time, a \"good\" interval being one that contains the parameter. We say that we are *95% confident* that the process results in a good interval. Our interval of estimates is called a *95% confidence interval* for the parameter, and 95% is called the *confidence level* of the interval.\n", "\n", "The method is called the *boostrap percentile method* because the interval is formed by picking off two percentiles of the bootstrapped estimates.\n", "\n", "The situation in the previous example was a bit unusual. Because we happened to know the value of the parameter, we were able to check whether an interval was good or a dud, and this in turn helped us to see that our process of estimation captured the parameter about 95 out of every 100 times we used it.\n", "\n", "But usually, data scientists don't know the value of the parameter. That is the reason they want to estimate it in the first place. In such situations, they provide an interval of estimates for the unknown parameter by using methods like the one we have developed. Because of statistical theory and demonstrations like the one we have seen, data scientists can be confident that their process of generating the interval results in a good interval a known percent of the time." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Estimating a Population Median\n", "\n", "We will now use the bootstrap method to estimate an unknown population median. You have encountered the dataset before. It comes from a sample of newborns in a large hospital system. we will treat it as if it were a simple random sample though the sampling was done in multiple stages. [Stat Labs](https://www.stat.berkeley.edu/~statlabs/) by Deborah Nolan and Terry Speed has details about a larger dataset from which this set is drawn. \n", "\n", "The table `births` contains the following variables for mother-baby pairs: the baby's birth weight in ounces, the number of gestational days (the number of days the mother was pregnant), the mother's age in completed years, the mother's height in inches, pregnancy weight in pounds, and whether or not the mother smoked during pregnancy." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "births = Table.read_table(path_data + 'baby.csv')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
Birth Weight | Gestational Days | Maternal Age | Maternal Height | Maternal Pregnancy Weight | Maternal Smoker | \n", "
---|---|---|---|---|---|
120 | 284 | 27 | 62 | 100 | False | \n", "
113 | 282 | 33 | 64 | 135 | False | \n", "
128 | 279 | 28 | 64 | 115 | True | \n", "
... (1171 rows omitted)
" ], "text/plain": [ "Birth Weight | Gestational Days | Ratio BW:GD | \n", "
---|---|---|
120 | 284 | 0.422535 | \n", "
113 | 282 | 0.400709 | \n", "
128 | 279 | 0.458781 | \n", "
108 | 282 | 0.382979 | \n", "
136 | 286 | 0.475524 | \n", "
138 | 244 | 0.565574 | \n", "
132 | 245 | 0.538776 | \n", "
120 | 289 | 0.415225 | \n", "
143 | 299 | 0.478261 | \n", "
140 | 351 | 0.39886 | \n", "
... (1164 rows omitted)
" ], "text/plain": [ "Birth Weight | Gestational Days | Ratio BW:GD\n", "120 | 284 | 0.422535\n", "113 | 282 | 0.400709\n", "128 | 279 | 0.458781\n", "108 | 282 | 0.382979\n", "136 | 286 | 0.475524\n", "138 | 244 | 0.565574\n", "132 | 245 | 0.538776\n", "120 | 289 | 0.415225\n", "143 | 299 | 0.478261\n", "140 | 351 | 0.39886\n", "... (1164 rows omitted)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ratios" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here is a histogram of the ratios." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAEcCAYAAAC4WdEMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de1RU9d4/8Pc4gOCgwoFhjDAVmVA4KoqG4kkFj/qkj6KkgRoVHfPWRSpJyUJFExHsEQ0Rj3rCW2lKJcaxUPGSKJlKGaahBinKICgIExcZ5veHy/1r5OIMzgWY92st1mr2/s7+fj6MrHd7z76ISktL1SAiIjIT7UxdABERkTEx+IiIyKww+IiIyKww+IiIyKww+IiIyKww+IiIyKww+IiIyKyYLPj69OkDOzu7ej8vvPACAECtViM6Ohq9evVCly5dMG7cOPz6668a26iurkZ4eDhcXV3h7OyM4OBgFBQUmKIdIiJqJUwWfBkZGbh06ZLwc/ToUYhEIkycOBEAEB8fj4SEBMTExODw4cOQSqWYNGkSysvLhW1EREQgNTUVmzdvRlpaGsrLyxEUFASVSmWqtkwuNzfX1CWYhLn2DZhv7+baN2C+veurb5MFn6OjI2QymfCTnp6Ojh07YuLEiVCr1UhMTERYWBgCAgLg4eGBxMREVFRUYM+ePQCAsrIybNu2DVFRUfDz84OXlxeSkpKQk5ODI0eOmKotIiJq4VrEd3xqtRrbtm1DUFAQOnTogPz8fCgUCvj7+wtjbGxs4Ovri6ysLABAdnY27t27pzHGxcUF7u7uwhgiIqKHWZi6AOD+Yc/8/HyEhIQAABQKBQBAKpVqjJNKpbh58yYAoKioCGKxGA4ODvXGFBUVNTlfWz9M0Nb7a4y59g2Yb+/m2jdgvr031LdcLtdpGy0i+JKTkzFgwAD07dtXY7lIJNJ4rVar6y17mDZjdP0ltSa5ubltur/GmGvfgPn2bq59A+bbu776Nvmhzlu3biEtLQ0vv/yysEwmkwFAvT234uJiYS/QyckJKpUKJSUljY4hIiJ6mMmDb+fOnWjfvj0CAwOFZd26dYNMJkNGRoawrKqqCidPnoSPjw8AwMvLC5aWlhpjCgoKcOnSJWEMERHRw0x6qFOtVmPr1q0IDAxEx44dheUikQhz5szB6tWrIZfL4ebmhri4OEgkEkyePBkA0LlzZ4SEhCAyMhJSqRT29vZYtGgRPD09MWLECBN1RERELZ1Jg+/48eO4cuUKNm7cWG/dvHnzUFlZifDwcJSWlsLb2xspKSkaAblixQqIxWKEhoaiqqoKw4YNw4YNGyAWi43ZBhERtSIiPoG9bTHnL7072Tmi+M5dU5cCR/tOkEntjTafOX/m5tg3YL6966vvFnFWJ5E+FN+5i6i1O0xdBiLfmm7U4CMi3Zj85BYiIiJjYvAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZYfAREZFZMWnwFRYWYvbs2ejZsydkMhl8fHzw/fffC+vVajWio6PRq1cvdOnSBePGjcOvv/6qsY3q6mqEh4fD1dUVzs7OCA4ORkFBgbFbISKiVsJkwVdaWooxY8ZArVZj9+7dyMrKwqpVqyCVSoUx8fHxSEhIQExMDA4fPgypVIpJkyahvLxcGBMREYHU1FRs3rwZaWlpKC8vR1BQEFQqlSnaIiKiFs7CVBOvXbsWXbp0QVJSkrCse/fuwn+r1WokJiYiLCwMAQEBAIDExETI5XLs2bMHoaGhKCsrw7Zt25CQkAA/Pz8AQFJSEvr06YMjR45g5MiRRu2JiIhaPpPt8X3zzTfw9vZGaGgo3Nzc8I9//AMbN26EWq0GAOTn50OhUMDf3194j42NDXx9fZGVlQUAyM7Oxr179zTGuLi4wN3dXRhDRET0Vybb48vLy8PmzZsxd+5chIWF4fz581iwYAEAYObMmVAoFACgcejzweubN28CAIqKiiAWi+Hg4FBvTFFRUaNz5+bm6rOVFqet99eYCmUFlEqlqctAhbLC6J+BuX7m5to3YL69N9S3XC7XaRsmC766ujr0798fixcvBgD069cPV69exaZNmzBz5kxhnEgk0nifWq2ut+xhjxqj6y+pNcnNzW3T/TUmNzcXthJbSCQSU5cCW4kt5PJuRpvPnD9zc+wbMN/e9dW3yQ51ymQyuLu7ayx7+umncf36dWE9gHp7bsXFxcJeoJOTE1QqFUpKShodQ0RE9FcmC77Bgwfj8uXLGssuX76Mrl27AgC6desGmUyGjIwMYX1VVRVOnjwJHx8fAICXlxcsLS01xhQUFODSpUvCGCIior8y2aHOuXPnYvTo0YiLi0NgYCB+/vlnbNy4ER9++CGA+4c458yZg9WrV0Mul8PNzQ1xcXGQSCSYPHkyAKBz584ICQlBZGQkpFIp7O3tsWjRInh6emLEiBGmao2IiFowkwXfgAEDsGPHDkRFRSE2NhYuLi54//33MWPGDGHMvHnzUFlZifDwcJSWlsLb2xspKSno2LGjMGbFihUQi8UIDQ1FVVUVhg0bhg0bNkAsFpuiLSIiauFMFnwAMGbMGIwZM6bR9SKRCBEREYiIiGh0jLW1NWJjYxEbG2uIEomIqI3hvTqJiMisMPiIiMisMPiIiMisMPiIiMisMPiIiMisaB18J06cQHFxcaPrS0pKcOLECb0URUREZChaB9/48eM17pDysKNHj2L8+PF6KYqIiMhQtA6+B48LakxNTQ3ateORUyIiatmavID97t27KCsrE17fvn0b165dqzeutLQUe/fuxRNPPKH/ComIiPSoyeBbv349Vq1aBeDRd1FRq9XCfTaJiIhaqiaDb8SIEbC2toZarUZUVBQCAwPRp08fjTEikQgdOnRA//79MXDgQIMWS0RE9LiaDL7Bgwdj8ODBAIDq6mqMHz8enp6eRimMiIjIELS+SfXChQsNWQcREZFRNBp8n332GQAgODgYIpFIeP0oU6dO1U9lREREBtBo8M2dOxcikQjPP/88rKysMHfu3EduTCQSMfiIiKhFazT4fvrpJwCAlZWVxmsiIqLWrNHge+qpp5p8TURE1BrxVitERGRWtD6rEwCOHDmC5ORk5OXl4c6dO/VuYyYSiZCdna3XAomIiPRJ6+BLTEzEokWL4OjoiIEDB6J3796GrIuIiMggtA6+hIQEDB06FHv37hVOeCEiImpttP6Or6SkBIGBgXoLvejoaNjZ2Wn8PP3008J6tVqN6Oho9OrVC126dMG4cePw66+/amyjuroa4eHhcHV1hbOzM4KDg1FQUKCX+oiIqG3SOvi8vLzwxx9/6HVyuVyOS5cuCT+ZmZnCuvj4eCQkJCAmJgaHDx+GVCrFpEmTUF5eLoyJiIhAamoqNm/ejLS0NJSXlyMoKAgqlUqvdRIRUduhdfB99NFH2LlzJ44dO6a3yS0sLCCTyYQfR0dHAPf39hITExEWFoaAgAB4eHggMTERFRUV2LNnDwCgrKwM27ZtQ1RUFPz8/ODl5YWkpCTk5OTgyJEjequRiIjaFq2/44uOjkanTp0wceJE9OzZE127doVYLNYYIxKJsHv3bq0nz8vLQ+/evWFpaYmBAwciMjIS3bt3R35+PhQKBfz9/YWxNjY28PX1RVZWFkJDQ5GdnY179+5pjHFxcYG7uzuysrIwcuRIresgIiLzoXXwXbx4ESKRCC4uLqiursbly5frjRGJRFpPPHDgQKxfvx5yuRzFxcWIjY3F6NGjcerUKSgUCgCAVCrVeI9UKsXNmzcBAEVFRRCLxXBwcKg3pqioqMm5c3Nzta6zNWrr/TWmQlkBpVJp6jJQoaww+mdgrp+5ufYNmG/vDfUtl8t12obWwXf+/HmdNvwoo0aN0ng9cOBAeHl5YefOnRg0aBCA+kGqVqsfGa7ajNH1l9Sa5Obmtun+GpObmwtbiS0kEompS4GtxBZyeTejzWfOn7k59g2Yb+/66rvF3LnF1tYWvXr1wtWrVyGTyQCg3p5bcXGxsBfo5OQElUqFkpKSRscQERE9TOvgu3btmlY/zVVVVYXc3FzIZDJ069YNMpkMGRkZGutPnjwJHx8fAPfPMrW0tNQYU1BQgEuXLgljiIiIHqb1oc6+fftq9R3e7du3tdreBx98gP/5n/+Bi4uL8B3fn3/+ialTp0IkEmHOnDlYvXo15HI53NzcEBcXB4lEgsmTJwMAOnfujJCQEERGRkIqlcLe3h6LFi2Cp6cnRowYoW1bRERkZrQOvk8++aRe8KlUKuTn5+Pzzz+Hk5MTZsyYofXEN27cwIwZM1BSUiLcBi09PV14CsS8efNQWVmJ8PBwlJaWwtvbGykpKejYsaOwjRUrVkAsFiM0NBRVVVUYNmwYNmzYUO9sUyIioge0Dr7p06c3ui4sLAz+/v6oqKjQeuItW7Y0uV4kEiEiIgIRERGNjrG2tkZsbCxiY2O1npeIiMybXk5usbW1xfTp07F+/Xp9bI6IiMhg9HZWp6WlpXCNHRERUUull+A7f/48NmzYAHd3d31sjoiIyGAe+6zOsrIy3L17F7a2tkhISNBrcURERPqmdfANHTq0XvCJRCLY2dnB1dUVzz//POzs7PReIBERkT7p9AR2IiKi1q7F3LKMiIjIGBh8RERkVhh8RERkVhh8RERkVrQ+uYWoMYpbd1B8565Ja6hQ1sDS0qQlEFEroVXwVVVVIT4+HoMGDYK/v7+ha6JWpvjOXUSt3WHSGpRKJRa9+aJJayCi1kGrQ53W1tb4v//7P1y/ft3Q9RARERmU1t/x9enTB1evXjVkLURERAandfBFRkZi69at+Pbbbw1ZDxERkUFpfXLL2rVrYWdnh6lTp8LZ2Rndu3eHjY2NxhiRSITdu3frvUgiIiJ90Tr4Ll68CJFIBBcXFwDAH3/8UW9MQzexJiIiakm0Dr7z588bsg4iIiKj4AXsRERkVnQKPpVKhd27d+ONN95AUFAQfvnlFwBAaWkpvvzySxQWFhqkSCIiIn3ROvjKysowevRozJo1C19//TXS09NRUlICAOjYsSMWLVqEjRs3GqxQIiIifdA6+JYuXYqLFy/iiy++QHZ2NtRqtbBOLBZj/PjxSE9PN0iRRERE+qJ18H3zzTeYOXMm/vnPfzZ49mbPnj1x7dq1ZheyevVq2NnZITw8XFimVqsRHR2NXr16oUuXLhg3bhx+/fVXjfdVV1cjPDwcrq6ucHZ2RnBwMAoKCppdBxERtW1aB19paSl69OjR6Hq1Wo2amppmFXH69GkkJyfD09NTY3l8fDwSEhIQExODw4cPQyqVYtKkSSgvLxfGREREIDU1FZs3b0ZaWhrKy8sRFBQElUrVrFqIiKht0zr4nnrqKVy4cKHR9SdOnICbm5vOBZSVleG1117DunXrYGdnJyxXq9VITExEWFgYAgIC4OHhgcTERFRUVGDPnj3Ce7dt24aoqCj4+fnBy8sLSUlJyMnJwZEjR3SuhYiI2j6tg2/KlCnYunUrTpw4ISx7cMgzKSkJ+/fvx7Rp03Qu4EGwDR8+XGN5fn4+FAqFxtMgbGxs4Ovri6ysLABAdnY27t27pzHGxcUF7u7uwhgiIqK/0voC9rfffhs//vgjJkyYADc3N4hEIixcuBC3b9+GQqHAuHHjMGvWLJ0mT05OxtWrV5GUlFRvnUKhAABIpVKN5VKpFDdv3gQAFBUVQSwWw8HBod6YoqKiRufNzc3Vqc7Wxtj9VShroFQqjTpnQ2pVtS2ijgplhdE/g7b+b7ox5to3YL69N9S3XC7XaRtaB5+lpSV2796NL774Al999RVEIhFqa2vRr18/BAYG4oUXXtDplmW5ubmIiorCf//7X1hZWTU67uFtqtXqR87zqDG6/pJak9zcXKP3l/NbPiQSiVHnfJhSqYSF2MLkdQCArcQWcnk3o81nis+8JTDXvgHz7V1ffev8BPYpU6ZgypQpjz3xDz/8gJKSEgwZMkRYplKpkJmZiS1btuDUqVMA7u/VPbg/KAAUFxcLe4FOTk5QqVQoKSmBo6OjxhhfX9/HrpGIiNoenYMPAH755Rfh0oWuXbvC09NT5xtUjxs3Dv3799dY9vrrr6Nnz55455134ObmBplMhoyMDAwYMADA/SfBnzx5ElFRUQAALy8vWFpaIiMjQwjjgoICXLp0CT4+Ps1pjYiI2jidgm/v3r1YvHgxbty4IVzALhKJ4OzsjMWLF+u0J2hnZ6dxFicAdOjQAfb29vDw8AAAzJkzB6tXr4ZcLoebmxvi4uIgkUgwefJkAEDnzp0REhKCyMhISKVS2NvbY9GiRfD09MSIESN0aY2IiMyE1sG3Y8cOvPHGG5DL5Vi6dCnc3NygVqtx5coVbN26FbNmzUJNTQ2mT5+ut+LmzZuHyspKhIeHo7S0FN7e3khJSUHHjh2FMStWrIBYLEZoaCiqqqowbNgwbNiwAWKxWG91EBFR2yEqLS1VP3oY4O3tDXt7e+zfvx/W1tYa66qqqjB27FiUlZXhzJkzBimUtGOqk1ui1u4w6pwPUyqVWPTmi1iz5UuT1gEAi+dNh1qrvyr9qFBWwFZiW2+5o30nyKT2xivEyMz1BA/AfHs3+sktBQUFmDlzZr3QAwBra2sEBQVhyZIlj10QUWtXVv6nUQNYqVQ2eDZr5FvT23TwETWX1hew9+rVS7h+riE3btyAu7u7XooiIiIyFK2DLyoqCsnJyfjyy/r/J7t3715s3boVy5Yt02txRERE+qb1oc5169bBwcEB//rXv7Bw4UL06NEDIpEIV69exa1bt9CzZ0+sXbsWa9euFd4jEomwe/dugxRORETUHFoH38WLFyESiYSLyW/cuAEAaN++PVxcXFBdXY1Lly5pvEfXa/uIiIgMTevgO3/+vCHrICIiMgqtv+MjIiJqCxh8RERkVhh8RERkVhh8RERkVhh8RERkVhh8RERkVrQOvn79+iEtLa3R9QcOHEC/fv30UhQREZGhaB18f/zxB5RKZaPrlUql8HBaIiKilkqnQ51N3Ynl8uXLGs/JIyIiaomavHPLzp078dlnnwmv4+LikJycXG9caWkpLly4gDFjxui/QiIiIj1qMviUSiUUCoXwuqysDHV1dRpjRCIROnTogJdffhkLFy40TJVERER60mTwvfbaa3jttdcAAH379sXKlSsxduxYoxRGRERkCFrfpPrnn382ZB1ERERGoXXwPVBeXo7r16/jzp07UKvV9dYPHTpUL4UREREZgtbBd+fOHSxYsABffvklVCpVvfVqtRoikQi3b9/Wa4FERET6pHXwvf3229i/fz9ee+01DB06FHZ2doasi4iIyCC0Dr6DBw9i1qxZ+Oijj/Qy8b///W/85z//ES5679WrF+bPny9cEqFWq7Fy5UokJyejtLQU3t7eiIuLQ+/evYVtVFdX44MPPsDevXtRVVWFYcOGYfXq1XjyySf1UiMREbU9Wl/AbmVlhZ49e+ptYmdnZyxduhRHjx5FRkYGhg0bhunTp+OXX34BAMTHxyMhIQExMTE4fPgwpFIpJk2ahPLycmEbERERSE1NxebNm5GWloby8nIEBQU1eCiWiIgI0CH4AgICkJ6erreJx40bh1GjRsHV1RVubm748MMPYWtri9OnT0OtViMxMRFhYWEICAiAh4cHEhMTUVFRgT179gC4f03htm3bEBUVBT8/P3h5eSEpKQk5OTk4cuSI3uokIqK2Revge/PNN1FYWIjZs2fj9OnTKCwsxK1bt+r9NIdKpcLevXuhVCrxzDPPID8/HwqFAv7+/sIYGxsb+Pr6IisrCwCQnZ2Ne/fuaYxxcXGBu7u7MIaIiOhhWn/H5+3tDZFIhOzsbOzevbvRcbqc1ZmTk4PRo0ejqqoKEokE27dvh6enpxBcUqlUY7xUKsXNmzcBAEVFRRCLxXBwcKg3pqioqMl5c3Nzta6xNTJ2fxXKmiZvYG4stapas62jofkqlBX8t96GmWvvDfUtl8t12obWwffee+81eZPq5pDL5Th+/DjKysqwb98+zJkzB/v37xfWPzzfg0smmqLNGF1/Sa1Jbm6u0fvL+S0fEonEqHM+TKlUwkJsYfI6ABi9DqVS2eB8thJbyOXdjFaHsZni33pLYa6966tvrYMvIiLisSd7mJWVFVxdXQEA/fv3x9mzZ7F+/XrMnz8fwP29OhcXF2F8cXGxsBfo5OQElUqFkpISODo6aozx9fXVe61ERNQ2NOsJ7CqVCrdv30Ztba1ei6mrq0NNTQ26desGmUyGjIwMYV1VVRVOnjwJHx8fAICXlxcsLS01xhQUFODSpUvCGCIioofpFHxnz57FxIkT4ezsDDc3N5w4cQIAUFJSghdeeAFHjx7VeltLlixBZmYm8vPzkZOTg6VLl+L777/HlClTIBKJMGfOHKxZswb79u3DhQsXMHfuXEgkEkyePBkA0LlzZ4SEhCAyMhJHjhzBTz/9hFmzZsHT0xMjRozQpS0iIjIjWh/q/OGHHzBhwgTIZDIEBwdj69atwjoHBwdUVFRg27ZtGD58uFbbUygUmDlzJoqKitCpUyd4enpiz549GDlyJABg3rx5qKysRHh4uHABe0pKisbDblesWAGxWIzQ0FDhAvYNGzZALBZr2xYREZkZrYNv2bJl6NmzJw4dOgSlUqkRfADw7LPPYteuXVpPnJiY2OR6kUiEiIiIJr9btLa2RmxsLGJjY7Wel4iIzJvWhzrPnj2LF198EdbW1g2eNfnkk09qPLSWiIioJdI6+Nq1a4d27RofrlAoYGNjo5eiiIiIDEXr4PPy8sKBAwcaXFdTU4MvvvgCzzzzjN4KIyIiMgStg++dd97BsWPH8MYbb+D8+fMAgMLCQhw8eBATJkzA77//jnfffddghRIREemD1ie3+Pn5ISkpCeHh4di5cycAYM6cOVCr1ejcuTM2bdqEQYMGGaxQIiIifdA6+ABg8uTJGDt2LDIyMnDlyhXU1dWhR48eGDlyJGxtbQ1VIxERkd7oFHwA0KFDB4wbN84QtRARERmc1t/xpaWlITw8vNH14eHhjZ78QkRE1FJoHXzr1q3Dn3/+2ej6qqoqxMfH66UoIiIiQ9E6+C5cuAAvL69G1/fr1w8XL17US1FERESGonXw1dbWorKystH1lZWVqK6u1ktRREREhqJ18Hl4eGDfvn2oq6urt66urg779u1Dr1699FocERGRvmkdfLNnz8aZM2cwdepUZGdno7q6GtXV1cjOzsa0adNw5swZzJo1y5C1EhERPTatL2d4/vnn8fvvvyM6Ohrp6ekA7j9BQa1WQyQSYcGCBQgKCjJYoURERPqg03V88+fPx+TJk5Gamoq8vDyo1Wr06NED48ePR/fu3Q1UIhERkf5oFXyVlZV44YUXEBQUhBdffBFvvvmmoesiIiIyCK2+47OxscFPP/0ElUpl6HqIiIgMSuuTW/7xj38gMzPTkLUQEREZnNbBFxMTg7Nnz+LDDz9EXl5eg5c1EBERtXRan9wyaNAgqNVqJCQkICEhAe3atYOlpaXGGJFIhBs3bui9SCIiIn3ROvgmTZoEkUhkyFqIiIgMTuvgS0xM1OvEH3/8MVJTU3H58mVYWVlh4MCBWLx4MTw8PIQxarUaK1euRHJyMkpLS+Ht7Y24uDj07t1bGFNdXY0PPvgAe/fuRVVVFYYNG4bVq1fjySef1Gu9RETUNmj9HZ++ff/99/jXv/6Fb7/9Fvv27YOFhQUmTpyIO3fuCGPi4+ORkJCAmJgYHD58GFKpFJMmTUJ5ebkwJiIiAqmpqdi8eTPS0tJQXl6OoKAgnoFKREQN0in4/vjjD7z11lvw8vJC165d8f333wMASkpK8O677yI7O1vrbaWkpODFF1+Eh4cHPD09kZSUhOLiYpw6dQrA/b29xMREhIWFISAgAB4eHkhMTERFRQX27NkDACgrK8O2bdsQFRUFPz8/eHl5ISkpCTk5OThy5IgurRERkZnQOvguXbqE4cOH4+uvv0bPnj2hVCqFvSoHBwecPn0amzZtanYhFRUVqKurg52dHQAgPz8fCoUC/v7+whgbGxv4+voiKysLAJCdnY179+5pjHFxcYG7u7swhoiI6K+0/o5v8eLF6NixIw4ePAixWAw3NzeN9aNHj8ZXX33V7EIWLlyIPn364JlnngEAKBQKAIBUKtUYJ5VKcfPmTQBAUVERxGIxHBwc6o0pKipqdK7c3Nxm19kaGLu/CmUNlEqlUedsSK2q1mzraGi+CmUF/623Yebae0N9y+VynbahdfBlZmZi/vz5cHJywu3bt+ut79q1qxBIunr//fdx6tQpHDhwAGKxWGPdw2eSPrgpdlMeNUbXX1Jrkpuba/T+cn7Lh0QiMeqcD1MqlbAQW5i8DgBGr0OpVDY4n63EFnJ5N6PVYWym+LfeUphr7/rqW6cH0Tb1x3znzp16oaWNiIgI7N27F/v27dO40bVMJgOAentuxcXFwl6gk5MTVCoVSkpKGh1DRET0Vzo9iPb48eMNrlOr1UhNTYWXl5dOky9YsAB79uzBvn378PTTT2us69atG2QyGTIyMoRlVVVVOHnyJHx8fAAAXl5esLS01BhTUFCAS5cuCWOIiIj+SutDnXPmzMGMGTOwatUqBAYGArj/5PXffvsN0dHROHfuHHbt2qX1xPPnz8euXbuwfft22NnZCd/pSSQS2NraQiQSYc6cOVi9ejXkcjnc3NwQFxcHiUSCyZMnAwA6d+6MkJAQREZGQiqVwt7eHosWLYKnpydGjBihw6+BiIjMhU4Por127Ro++ugjrFy5UlgGAGKxGMuXL8eoUaO0nvjBGaABAQEayxcsWICIiAgAwLx581BZWYnw8HDhAvaUlBR07NhRGL9ixQqIxWKEhoYKF7Bv2LChWYddiYio7dPpQbRhYWGYPHky9u3bh6tXr6Kurg49evTAhAkT0K2bbl+il5aWPnKMSCRCRESEEIQNsba2RmxsLGJjY3Wan4iIzNMjg6+6utjI7WgAABXcSURBVBppaWnIy8vD3/72N4wZMwZz5841Rm1ERER612TwKRQKjB07Fr///jvUajWA+9/B7dq1C0OHDjVKgURERPrU5Fmdy5cvR15eHubOnYtdu3YhOjoa7du3x3vvvWes+oiIiPSqyT2+w4cPY+rUqVi+fLmwzMnJCTNmzEBBQQGfgEBERK1Ok3t8CoWi3vVwgwcPhlqtxvXr1w1aGBERkSE0GXwqlQrW1tYayx68rqqqMlxVREREBvLIszrz8vJw5swZ4fXdu3cB3L9nmq2tbb3x3t7eeiyPiIhIvx4ZfNHR0YiOjq63/OETXB7cGLqhG1gTERG1FE0GX0JCgrHqICIiMoomg2/atGnGqoOIiMgotH46AxERUVvA4CMiIrPC4CMiIrPC4CMiIrOi02OJiKj1EImAnN/yTV0GHO07QSa1N3UZRAIGH1EbVVb+J9Zs+dLUZSDyrekMPmpReKiTiIjMCoOPiIjMCg91tmKKW3dQfOeuxrIKZY3Rv9eprqkx6nxERI+DwdeKFd+5i6i1OzSWKZVKSCQSo9YR9uoko85HRPQ4eKiTiIjMCoOPiIjMikmD78SJEwgODkbv3r1hZ2eHHTs0D9up1WpER0ejV69e6NKlC8aNG4dff/1VY0x1dTXCw8Ph6uoKZ2dnBAcHo6CgwJhtEBFRK2LS4FMqlfDw8MDKlSthY2NTb318fDwSEhIQExODw4cPQyqVYtKkSSgvLxfGREREIDU1FZs3b0ZaWhrKy8sRFBQElUplzFaIiKiVMGnwjR49GpGRkQgICEC7dpqlqNVqJCYmIiwsDAEBAfDw8EBiYiIqKiqwZ88eAEBZWRm2bduGqKgo+Pn5wcvLC0lJScjJycGRI0dM0BEREbV0LfY7vvz8fCgUCvj7+wvLbGxs4Ovri6ysLABAdnY27t27pzHGxcUF7u7uwhgiIqK/arGXMygUCgCAVCrVWC6VSnHz5k0AQFFREcRiMRwcHOqNKSoqanTbubm5eq7WNCqUNVAqlfWWN7TMkGpVtUafk3Voami+lvL7qFBWGOxvrq38LTeHufbeUN9yuVynbbTY4HtAJBJpvFar1fWWPexRY3T9JbVUOb/l17tmzxTX8VmILYw+58OUSmWLqAMw/u+jsc+8pfw+bCW2kMu76X27ubm5beZvWVfm2ru++m6xhzplMhkA1NtzKy4uFvYCnZycoFKpUFJS0ugYIiKiv2qxwdetWzfIZDJkZGQIy6qqqnDy5En4+PgAALy8vGBpaakxpqCgAJcuXRLGEBER/ZVJD3VWVFTg6tWrAIC6ujpcv34dP//8M+zt7dG1a1fMmTMHq1evhlwuh5ubG+Li4iCRSDB58mQAQOfOnRESEoLIyEhIpVLY29tj0aJF8PT0xIgRI0zYGRERtVQmDb5z585h/Pjxwuvo6GhER0dj6tSpSExMxLx581BZWYnw8HCUlpbC29sbKSkp6Nixo/CeFStWQCwWIzQ0FFVVVRg2bBg2bNgAsVhsipaIiKiFM2nwPfvssygtLW10vUgkQkREBCIiIhodY21tjdjYWMTGxhqiRCIiamNa7Hd8REREhsDgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis2LSB9ESUdsnEgE5v+XrfbsVyhqdt+to3wkyqb3ea6HWhcFHRAZVVv4n1mz5Uu/bVSqVkEgkOr0n8q3pDD5i8OlKcesOiu/cNXUZAIDqmhpTl0BE1Oow+HRUfOcuotbuMHUZAICwVyeZugQiolaHJ7cQEZFZaTPBt2nTJvTt2xcymQzDhw9HZmamqUsiIqIWqE0EX0pKChYuXIh3330Xx44dwzPPPIMpU6bg2rVrpi6NiIhamDYRfAkJCZg2bRpefvlluLu7IzY2FjKZDFu2bDF1aURE1MKISktL1aYu4nHU1NTgiSeewObNmzFx4kRh+fz583HhwgWkpaWZsDoiImppWv0eX0lJCVQqFaRSqcZyqVSKoqIiE1VFREQtVasPvgdEIpHGa7VaXW8ZERFRqw8+BwcHiMXient3xcXF9fYCiYiIWn3wWVlZwcvLCxkZGRrLMzIy4OPjY6KqiIiopWoTd255/fXXMWvWLHh7e8PHxwdbtmxBYWEhQkNDTV0aERG1MK1+jw8AAgMDER0djdjYWDz77LM4deoUdu/ejaeeesrUpemdLhfqX7x4Ef/7v/8LuVwOmUyGfv36ISoqCjWt8B6fzb1BwZUrV+Di4oInn3zSwBUaji695+fnw87Ort7PwYMHjVixfuj6mavVaqxfvx6DBg2Ck5MT3N3dsWTJEuMUq2e69B4dHd3gZ25nZ4dbt24ZserHp+tnfujQIYwaNQouLi5wdXXF1KlTcfny5UfO0+ovZzAnKSkpmDlzJlavXo3Bgwdj06ZN2LlzJ06dOoWuXbvWG3/16lWcPHkSffv2RefOnfHLL79g3rx5mDp1KqKiokzQQfPo2vcDNTU1GDVqFGQyGU6cOIGCggIjVq0fuvaen5+Pfv36Ye/evfj73/8uLLe3t4eVlZUxS38szfnM33//fXz77beIioqCp6cnysrKoFAoMHr0aCNX/3h07b2iogJKpVJj2auvvgqRSIT9+/cbq+zHpmvfeXl58PHxwaxZs/DKK6+goqICixcvRl5eHs6dO9fkXAy+VmTkyJHw9PTE2rVrhWUDBgxAQEAAFi9erNU23n//fZw+fRrp6emGKlPvmtt3REQEysrKMHToULz33nutMvh07f1B8GVkZKB///7GLFWvdO07NzcXQ4YMwYkTJ+Du7m7MUvXucf/Or1+/jr59+yIpKQlTpkwxZKl6pWvfX3/9NUJDQ3Hr1i2IxWIAwLFjxzBhwgRcuXIFDg4Ojc7VJg51moOamhpkZ2fD399fY7m/vz+ysrK02sbVq1dx6NAhDB061BAlGkRz+/7222/x7bffIiYmxtAlGszjfOYhISFwc3PDmDFj8PXXXxuyTL1rTt9paWno3r07Dh48iH79+qFPnz6YPXt2qzvUp4+/823btqFz586YMGGCIUo0iOb07eXlBUtLS2zduhUqlQrl5eX47LPPMGDAgCZDD2DwtRqPc6H+6NGjIZPJMGDAAAwePBiRkZGGLFWvmtN3YWEh5s2bh6SkJHTs2NEYZRpEc3q3tbXFsmXL8J///AdffPEFhg0bhtDQUOzatcsYJetFc/rOy8vDtWvXkJKSgvXr1yMpKQm5ubkIDg5GXV2dMcrWi8e9IUddXR127NiB4OBgtG/f3lBl6l1z+u7WrRu+/PJLREdHw8nJCU899RQuXLig1b91Bl8r05wL9bds2YKjR49i06ZNSE9Px5o1awxZokHo0vfMmTPx6quvYtCgQcYozeB06d3BwQFvvvkmBg0ahP79+2PRokUIDQ1FfHy8MUrVK136rqurQ3V1NZKSkjB06FD4+voiKSkJZ86cwdmzZ41Rrl4194Yc6enpuH79Ol566SVDlWZQuvStUCjw5ptvIjg4GIcPH8b+/ftha2uLV1555ZH/s9MmLmcwB49zob6LiwsAoFevXlCpVHjrrbfw1ltvwcKi5X/8zen72LFjOHHihHCYU61Wo66uDg4ODli9ejVeeeUVQ5etF/q6OYO3tzd27GgZD0/WRnP6lslksLCwgJubm7CsZ8+esLCwwPXr1zFw4ECD1qwvj/uZf/rpp/Dx8UHv3r0NVaJBNKfvf//73+jQoYPGiXobN26Ep6cnsrKyMGTIkEbn4x5fK6GvC/Xr6upQW1sLlUql7xINojl9Z2Zm4vjx48LP+++/DxsbGxw/flzjRuYtnb4+8/Pnz0Mmk+m7PINpTt+DBw9GbW0tfv/9d2FZXl4eamtrmzzzt6V5nM/85s2b+O6771rl3l5z+q6srBROanngwWvu8bUhj7pQf+nSpThz5gz27dsHAPj8889hbW0NDw8PWFlZ4dy5c4iKikJAQECrOv6va98eHh4a7z937hzatWtXb3lroGvvO3fuhKWlJfr27Yt27drhwIED2LRpU6u7nk3XvkeMGIF+/frh9ddfR3R0NID7Z/UOHDiw1Z3dqmvvD2zfvh0SiQSTJk0yRdmPTde+R48ejfXr12PlypWYMmUKysvLsWzZMri4uMDLy6vJuRh8rUhgYCBu376N2NhYKBQK9O7dW+NC/cLCQo3/47WwsMDHH3+Mq1evQq1Wo2vXrpgxYwbmzp1rqhaaRde+25Lm9B4XF4dr165BLBajZ8+e+OSTTxAUFGSK8ptN177btWuHXbt2YcGCBRg3bhysra3h5+eHjz76CO3ata4DW835zNVqNbZt24YpU6agQ4cOpij7sena9/Dhw7Fp0ybEx8dj3bp1sLa2xsCBA7Fnzx5IJJIm5+J1fEREZFZa1/8KERERPSYGHxERmRUGHxERmRUGHxERmRUGHxERmRUGHxERmRUGH9FDduzYofEwT6lUCi8vL0RFRaGqqqpZ24yOjsbRo0frLZ8zZw769OnzuCUDaPghtK6urhg7diwOHTokjAsMDES/fv3qvX/NmjWws7Nr8Cbmfn5+GDFixCNr+OGHH/Dqq6/Cw8MDUqkUXbt2hZ+fH5YvX47CwkKNsX+t09HRET179sRzzz2HVatWtbqnKlDrwuAjakRycjLS09Oxe/dujBw5Eh9//DE+/PDDZm0rJiYGx48fr7f8vffew/bt2x+3VA3vvPMO0tPTkZ6ejnXr1sHKygpBQUE4ffo0AMDX1xf5+fn1nk+YmZmJDh061HvqdUVFBX7++Wf4+vo2Oe+6deswZswYFBcXY9GiRfjqq6+wefNmjBw5EsnJyXjjjTfqvWfatGlIT0/HN998g08++QS+vr7YuHEjBg8erPVjeIh0xTu3EDWiT58+cHV1BXB/j+fKlSvYvn07YmJi9HY3kB49euhlO3/VvXt3jSdT+Pn5oXv37ti/fz8GDRokBFhmZqbwoNK6ujqcOnUKL730EjZt2gSlUinc/SIrKwsqlarJ5zgeO3YMkZGRmD17tnDLsAdGjx6Nt99+G1999VW99zk7O2vU+txzz2HWrFl47rnnEBISgnPnzj3yLhxEuuIeH5GW+vXrh8rKSpSUlAjLDh8+jClTpsDd3R1PPPEEhgwZgnXr1mncBNzOzg7A/VuJPTi09yAcGjrUWVhYiFmzZsHV1RVOTk7w9fV9rOfpWVtbw9LSEvfu3QNw/2kN1tbWGnt258+fx927d/H6669DLBYLe4fA/YAUiURN7vHFx8fDwcEBS5cubXC9RCLB9OnTtarXyckJUVFRKCoqwt69e7V6D5EuGHxEWvrjjz/QqVMn/O1vfxOW5eXlYdiwYfjkk0+wa9cuTJ06FTExMVi2bJkwJj09HcD/P6yXnp7e6B30lUolxo0bh4MHDyIyMhI7duyAp6cnZs2ahU8//VSrOh88gaO2thY3b97EokWLUFVVJTyRu3379vD29tYIvszMTLi5uaFr164YMGAATpw4obGud+/esLe3B3D/+0o7Ozvk5+cDAGpra3HixAn4+fnByspKqxofxd/fHxYWFjzcSQbBQ51EjVCpVKitrUVFRQVSU1Oxb98+REdHazwK5dVXXxX+W61Ww9fXFzU1NVi3bh0iIyPRrl074VDew4f1GrJjxw5cuXIFqampePbZZwEAo0aNQlFREZYvX46QkJB6j2J5WFhYGMLCwoTX7du3x5o1azB48GBhma+vL2JjY1FcXAxHR0dkZmYKzy8bMmSIEIpVVVU4e/asRlC3a9cOYrFYeEDo7du3UVVVJTz38a9qa2s1Xmv7DEgbGxs4ODhAoVBoNZ5IF9zjI2rEoEGD4OjoiO7du+PNN99EaGgoZs6cqTGmsLAQYWFh+Pvf/w6pVApHR0csX74cZWVlzTozMTMzE87OzkLoPfDCCy+guLgYFy9efOQ25s+fj4yMDGRkZCAlJQUvvfQSwsLCNA4b/uMf/xDmA4BTp04JwTd48GCcOXMGNTU1+PHHH1FdXa1xmHPBggUoKSkR7pqvVjd8n3uFQgFHR0eNn4eDsCnaPnWcSFfc4yNqxPbt2/Hkk0+iuLgYCQkJ2LRpE7y9vTF16lQA9w8pTp06FYWFhVi4cCHkcjlsbGzwzTffIC4urlmXPty5c6fBh8Y+WHbnzp1HbuOpp57SeAadv78/8vPzERERgcDAQIhEIgwaNAiWlpbIzMyEu7s7bt26JYSbj48PampqcObMGSEYm/p+z8HBAdbW1rh+/Xq95Q8eLPrpp58iOTn5kbU/8OC71Nb0AF1qPRh8RI3w8PAQzuocNmwYhg4disjISEyYMAESiQS///47zp07h6SkJI3n3f33v/9t9pz29va4fPlyveUPDvn99ftFXfTq1Qvfffcdbt26BScnJ3To0AFeXl5C8D3xxBPo3r07AKBz587w8PBAZmam8N1fUwFkYWEBX19fZGRkoKamRviez8LCQgjgAwcO6FTvoUOHoFKpNA7PEukLD3USaaF9+/aIiorCrVu3sHnzZgDAn3/+CQCwtLQUxt27dw9ffPFFvfdbWVmhsrLykfMMHToUBQUFOHXqlMbyPXv2QCqVwt3dvVn15+TkwMrKCp06ddKY65dffsGBAweEw5wPDBkyBMeOHcPp06cfef0eAMybNw8lJSVYvHhxs+r7q1u3bmHx4sXo0qULnn/++cfeHtHDuMdHpKWxY8diwIABWLduHV577TW4u7uja9euWLZsGcRiMSwsLLB+/foG3+vu7o7vvvsO//znP2FnZ4cuXbrgiSeeqDdu2rRp2LBhA0JCQvDhhx/C2dkZu3fvRkZGBtasWfPIE1uA+2eaPrgcobS0FAcOHMChQ4cwY8YMWFtbC+N8fX2xZs0afPfdd1i1apXGNoYMGYJNmzYJJ+z8VUxMDFatWoVz584J3/MNHz4cS5YswZIlS5CTk4Pg4GB069YN1dXVuHz5MlJSUiCRSOp9Z3fjxg2cPn0adXV1uHPnDn788UckJydDrVbj888/h42NzSP7JdIVg49IBx988AECAwOxZcsWvP7669ixYwfee+89zJ49G/b29pg+fTq6du2Kt956S+N9sbGxWLBgAYKDg1FdXY0FCxYgIiKi3vYlEgm++eYbREZGYsmSJaioqICbm1u9w6lN+fjjj/Hxxx8DADp16oTu3bsjLi4Or7zyisa4wYMHQywWQ6VSNbjH9+CklYeDr66uDiqVqt5JLfPmzYOPjw82bNiAZcuWobi4GNbW1nBzc8OkSZPw6quv1gvunTt3YufOnbCwsECnTp3w9NNPY+bMmQgNDYWjo6NW/RLpSlRaWtrwKVlERERtEL/jIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis8LgIyIis/L/ANUsJThMRyXZAAAAAElFTkSuQmCC\n", "text/plain": [ "Birth Weight | Gestational Days | Ratio BW:GD | \n", "
---|---|---|
116 | 148 | 0.783784 | \n", "