diff --git a/your-code/.ipynb_checkpoints/main-checkpoint.ipynb b/your-code/.ipynb_checkpoints/main-checkpoint.ipynb
new file mode 100644
index 0000000..8e8dcf6
--- /dev/null
+++ b/your-code/.ipynb_checkpoints/main-checkpoint.ipynb
@@ -0,0 +1,1729 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Understanding Descriptive Statistics\n",
+ "\n",
+ "Import the necessary libraries here:"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "from sklearn.linear_model import LinearRegression\n",
+ "import pandas as pd\n",
+ "import numpy as np\n",
+ "import matplotlib.pyplot as plt\n",
+ "from pylab import rcParams # library for plotting histograms"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Challenge 1\n",
+ "#### 1.- Define a function that simulates rolling a dice 10 times. Save the information in a dataframe.\n",
+ "**Hint**: you can use the *choices* function from module *random* to help you with the simulation."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "import random\n",
+ "\n",
+ "def rolling_dice(times):\n",
+ " \n",
+ " dice = [1,2,3,4,5,6]\n",
+ " result = []\n",
+ " counter = 0\n",
+ " while counter < times:\n",
+ " result.append(random.choices(dice))\n",
+ " counter += 1\n",
+ " result_df = pd.DataFrame(result)\n",
+ " result_df[0] = result_df.rename(columns={0: \"value\"}, inplace=True)\n",
+ " result_df = result_df[[\"value\"]]\n",
+ " return result_df"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 5 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 6 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 7 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 8 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 9 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 6 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 7 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 8 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ " 5 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 9 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# plotting the value\n",
+ "result.plot(kind='bar', alpha=0.7)\n",
+ "plt.xlabel('roll(index)')\n",
+ "plt.ylabel('value')\n",
+ "plt.title('')\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 3.- Calculate the frequency distribution and plot it. What is the relation between this plot and the plot above? Describe it with words."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 6,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "4 4\n",
+ "2 3\n",
+ "6 2\n",
+ "5 1\n",
+ "Name: value, dtype: int64"
+ ]
+ },
+ "execution_count": 6,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "# calculating FD\n",
+ "frequency_distribution = result['value'].value_counts()\n",
+ "frequency_distribution"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "metadata": {
+ "scrolled": true
+ },
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAHACAYAAABeV0mSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA1oUlEQVR4nO3de1xVdb7/8fdOZCMKKigIeS3LC+YNnMLESyQmHnS8nJOTeUltxtHSJE7ngc2ccmpCJ3XMLNBJxUveipx0MEdLISs9ieKtTK1UyMBrgVJubuv3hw/3zy14YQvszfL1fDz2H+u7vmutz9qr4t13fddeFsMwDAEAAJjEXa4uAAAAoDIRbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQboAqlJycLIvFUu4nLi7O1eWZyrXftZeXl5o0aaI+ffooISFBp0+fLrPNyy+/LIvFUqHj/PLLL3r55ZeVlpZWoe3KO1bLli31H//xHxXaz82sXLlSc+fOLXedxWLRyy+/XKnHA9yRh6sLAO4ES5YsUdu2bR3agoODXVSNuV35rouKinT69Gl99tlnmjlzpmbNmqU1a9bo0UcftfcdP368HnvssQrt/5dfftH06dMlSb17977l7Zw5ljNWrlypgwcP6rnnniuzbseOHWratGmV1wC4GuEGqAYdOnRQWFjYLfUtKiqSxWKRhwf/ejrj2u966NChmjp1qnr06KEhQ4bo6NGjCgwMlCQ1bdq0yv/Y//LLL/L29q6WY93MQw895NLjA9WF21KAC6WlpclisWj58uV6/vnndffdd8tqterbb7+VJH388ceKjIyUr6+vvL299fDDD+uTTz4ps5/U1FR17txZVqtVrVq10qxZs8rcBjl+/LgsFouSk5PLbF/e7YqjR4/qiSeeUEBAgKxWq9q1a6e33nqr3PpXrVqlF198UcHBwfL19dWjjz6qw4cPlznOpk2bFBkZqfr168vb21vt2rVTQkKCJGn58uWyWCzasWNHme3+8pe/qHbt2vrxxx9v+p2Wp3nz5po9e7YuXLigBQsW2NvLu1W0detW9e7dW/7+/qpTp46aN2+uoUOH6pdfftHx48fVuHFjSdL06dPtt8DGjBnjsL89e/Zo2LBhatiwoe69997rHuuKdevWqWPHjvLy8tI999yjefPmOay/csvt+PHjDu1Xvv8rt8h69+6t1NRUnThxwuEW3RXlXeeDBw9q0KBBatiwoby8vNS5c2ctXbq03OPc6nUGXI1wA1SDkpISFRcXO3yuFh8fr6ysLCUlJWnDhg0KCAjQihUrFBUVJV9fXy1dulRr166Vn5+f+vXr5xBwPvnkEw0aNEg+Pj5avXq1Xn/9da1du1ZLlixxut6vv/5a3bp108GDBzV79mz961//0oABAzR58mT7LZmrTZs2TSdOnNA777yjhQsX6ujRo4qJiVFJSYm9z6JFixQdHa3S0lL7eU6ePFk//PCDJOnxxx9XkyZNygSo4uJiLViwQIMHD76tW3nR0dGqVauWPv300+v2OX78uAYMGCBPT08tXrxYmzZt0owZM1S3bl0VFhYqKChImzZtkiSNGzdOO3bs0I4dO/TnP//ZYT9DhgxR69at9d577ykpKemGde3du1fPPfecpk6dqnXr1ql79+6aMmWKZs2aVeFzfPvtt/Xwww+rSZMm9trKC4tXHD58WN27d9dXX32lefPm6YMPPlD79u01ZswY/e1vfyvT/1auM+AWDABVZsmSJYakcj9FRUXGtm3bDElGz549HbYrKCgw/Pz8jJiYGIf2kpISo1OnTsZvfvMbe9uDDz5oBAcHG7/++qu9LT8/3/Dz8zOu/lf82LFjhiRjyZIlZeqUZLz00kv25X79+hlNmzY18vLyHPo988wzhpeXl3H+/HnDMAx7/dHR0Q791q5da0gyduzYYRiGYVy4cMHw9fU1evToYZSWll73+3rppZcMT09P49SpU/a2NWvWGJKM9PT0625nGP//u961a9d1+wQGBhrt2rVzON7V39H7779vSDL27t173X2cOXOmzPd17f7+93//97rrrtaiRQvDYrGUOV7fvn0NX19fo6CgwOHcjh075tDvyve/bds2e9uAAQOMFi1alFv7tXUPHz7csFqtRlZWlkO//v37G97e3sbPP//scJybXWfAXTByA1SDZcuWadeuXQ6fq+fUDB061KH/F198ofPnz2v06NEOoz2lpaV67LHHtGvXLhUUFKigoEC7du3SkCFD5OXlZd/ex8dHMTExTtV66dIlffLJJxo8eLC8vb0djh8dHa1Lly5p586dDtsMHDjQYbljx46SpBMnTtjPJz8/XxMnTrzh00l//OMfJUn/+Mc/7G3z58/XAw88oJ49ezp1PlczDOOG6zt37ixPT0/9/ve/19KlS/X99987dZxrr+eNhISEqFOnTg5tTzzxhPLz87Vnzx6njn+rtm7dqsjISDVr1syhfcyYMfrll1/KjPrc7DoD7oJwA1SDdu3aKSwszOFztaCgIIflU6dOSZKGDRum2rVrO3xmzpwpwzB0/vx5/fTTTyotLVWTJk3KHLO8tltx7tw5FRcX68033yxz7OjoaEnS2bNnHbbx9/d3WLZarZKkX3/9VZJ05swZSbrphNrAwEA9/vjjWrBggUpKSrR//35t375dzzzzjFPncrWCggKdO3fuhre27r33Xn388ccKCAjQpEmTdO+99+ree+/VG2+8UaFjXXs9b+RG1+7cuXMVOm5FnTt3rtxar3xH1x7/ZtcZcBc8jgG4gWtHMxo1aiRJevPNN6/7hEtgYKD9yarc3Nwy669tuzKyY7PZHNqv/QPWsGFD1apVSyNHjtSkSZPKPXarVq1ucDZlXZmEe2V+zY1MmTJFy5cv14cffqhNmzapQYMGGjFiRIWOV57U1FSVlJTc9PHtiIgIRUREqKSkRBkZGXrzzTf13HPPKTAwUMOHD7+lY1Xkt3NudO2uhInrXbtrQ2ZF+fv7Kycnp0z7lYnbV/45BGoawg3ghh5++GE1aNBAX3/99Q1HLTw9PfWb3/xGH3zwgV5//XX7H8ELFy5ow4YNDn0DAwPl5eWl/fv3O7R/+OGHDsve3t7q06ePMjMz1bFjR3l6et72+XTv3l3169dXUlKShg8ffsM//qGhoerevbtmzpypgwcP6ve//73q1q17W8fPyspSXFyc6tevrz/84Q+3tE2tWrX04IMPqm3btnr33Xe1Z88eDR8+vNJHK7766ivt27fP4dbUypUr5ePjo65du0q6/GN/krR//361adPG3m/9+vVl9me1Wm+5tsjISK1bt04//vijw4jWsmXL5O3tzaPjqLEIN4Abqlevnt58802NHj1a58+f17BhwxQQEKAzZ85o3759OnPmjBITEyVJr7zyih577DH17dtXzz//vEpKSjRz5kzVrVtX58+ft+/TYrHoySef1OLFi3XvvfeqU6dO+vLLL7Vy5coyx3/jjTfUo0cPRURE6I9//KNatmypCxcu6Ntvv9WGDRu0devWCp/P7NmzNX78eD366KN6+umnFRgYqG+//Vb79u3T/PnzHfpPmTJFjz/+uCwWiyZOnFihYx08eNA+R+j06dPavn27lixZolq1amndunX2UaTyJCUlaevWrRowYICaN2+uS5cuafHixZJk//E/Hx8ftWjRQh9++KEiIyPl5+enRo0a2QNIRQUHB2vgwIF6+eWXFRQUpBUrVmjLli2aOXOmvL29JUndunVTmzZtFBcXp+LiYjVs2FDr1q3TZ599VmZ/DzzwgD744AMlJiYqNDRUd91113V/Y+mll17Sv/71L/Xp00f/+7//Kz8/P7377rtKTU3V3/72N9WvX9+pcwJcztUzmgEzu9kTPFeeQnnvvffKXZ+enm4MGDDA8PPzM2rXrm3cfffdxoABA8r0X79+vdGxY0fD09PTaN68uTFjxoxyn87Jy8szxo8fbwQGBhp169Y1YmJijOPHj5f79M+xY8eMsWPHGnfffbdRu3Zto3Hjxkb37t2NV1999ab1X+/JrI0bNxq9evUy6tata3h7exvt27c3Zs6cWea8bTabYbVajccee6zc76U81z6Z5unpaQQEBBi9evUyXnvtNeP06dNltrn2O9qxY4cxePBgo0WLFobVajX8/f2NXr16GevXr3fY7uOPPza6dOliWK1WQ5IxevRoh/2dOXPmpscyjMtPSw0YMMB4//33jZCQEMPT09No2bKlMWfOnDLbHzlyxIiKijJ8fX2Nxo0bG88++6yRmppa5mmp8+fPG8OGDTMaNGhgWCwWh2OWd50PHDhgxMTEGPXr1zc8PT2NTp06lbluFb3OgKtZDOMmjw8AqJFefvllTZ8+/aZPCLmjDRs2aODAgUpNTbVPYgaAW8VtKQBu4+uvv9aJEyf0/PPPq3Pnzurfv7+rSwJQA/EoOAC3MXHiRA0cOFANGzbUqlWrKvzGbgCQJG5LAQAAU2HkBgAAmArhBgAAmArhBgAAmMod97RUaWmpfvzxR/n4+DBZEQCAGsIwDF24cEHBwcG6664bj83cceHmxx9/LPMGXAAAUDNkZ2ff9CW8d1y48fHxkXT5y/H19XVxNQAA4Fbk5+erWbNm9r/jN3LHhZsrt6J8fX0JNwAA1DC3MqWECcUAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBU3CbcJCQkyGKx6Lnnnrthv/T0dIWGhsrLy0v33HOPkpKSqqdAAABQI7hFuNm1a5cWLlyojh073rDfsWPHFB0drYiICGVmZmratGmaPHmyUlJSqqlSAADg7lwebi5evKgRI0boH//4hxo2bHjDvklJSWrevLnmzp2rdu3aafz48Ro7dqxmzZpVTdUCAAB35/JwM2nSJA0YMECPPvroTfvu2LFDUVFRDm39+vVTRkaGioqKqqpEAABQg3i48uCrV6/Wnj17tGvXrlvqn5ubq8DAQIe2wMBAFRcX6+zZswoKCiqzjc1mk81msy/n5+ffXtEAAMCtuSzcZGdna8qUKdq8ebO8vLxueTuLxeKwbBhGue1XJCQkaPr06c4XWsnGJd9akDObRWO6uboEAMAdwmW3pXbv3q3Tp08rNDRUHh4e8vDwUHp6uubNmycPDw+VlJSU2aZJkybKzc11aDt9+rQ8PDzk7+9f7nHi4+OVl5dn/2RnZ1fJ+QAAAPfgspGbyMhIHThwwKHtqaeeUtu2bfU///M/qlWrVpltwsPDtWHDBoe2zZs3KywsTLVr1y73OFarVVartfIKBwAAbs1l4cbHx0cdOnRwaKtbt678/f3t7fHx8Tp58qSWLVsmSZowYYLmz5+v2NhYPf3009qxY4cWLVqkVatWVXv9AADAPbn8aakbycnJUVZWln25VatW2rhxo9LS0tS5c2e98sormjdvnoYOHerCKgEAgDtx6dNS10pLS3NYTk5OLtOnV69e2rNnT/UUBAAAahy3HrkBAACoKMINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFZeGm8TERHXs2FG+vr7y9fVVeHi4Pvroo+v2T0tLk8ViKfP55ptvqrFqAADgzjxcefCmTZtqxowZat26tSRp6dKlGjRokDIzMxUSEnLd7Q4fPixfX1/7cuPGjau8VgAAUDO4NNzExMQ4LP/1r39VYmKidu7cecNwExAQoAYNGlRxdQAAoCZymzk3JSUlWr16tQoKChQeHn7Dvl26dFFQUJAiIyO1bdu2aqoQAADUBC4duZGkAwcOKDw8XJcuXVK9evW0bt06tW/fvty+QUFBWrhwoUJDQ2Wz2bR8+XJFRkYqLS1NPXv2LHcbm80mm81mX87Pz6+S8wAAAO7B5eGmTZs22rt3r37++WelpKRo9OjRSk9PLzfgtGnTRm3atLEvh4eHKzs7W7NmzbpuuElISND06dOrrH4AAOBeXH5bytPTU61bt1ZYWJgSEhLUqVMnvfHGG7e8/UMPPaSjR49ed318fLzy8vLsn+zs7MooGwAAuCmXj9xcyzAMh9tIN5OZmamgoKDrrrdarbJarZVRGgAAqAFcGm6mTZum/v37q1mzZrpw4YJWr16ttLQ0bdq0SdLlUZeTJ09q2bJlkqS5c+eqZcuWCgkJUWFhoVasWKGUlBSlpKS48jQAAIAbcWm4OXXqlEaOHKmcnBzVr19fHTt21KZNm9S3b19JUk5OjrKysuz9CwsLFRcXp5MnT6pOnToKCQlRamqqoqOjXXUKAADAzVgMwzBcXUR1ys/PV/369ZWXl+fwQ4DVZVzyrmo/pjtYNKabq0sAANRgFfn77fIJxQAAAJWJcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEzFpeEmMTFRHTt2lK+vr3x9fRUeHq6PPvrohtukp6crNDRUXl5euueee5SUlFRN1QIAgJrApeGmadOmmjFjhjIyMpSRkaFHHnlEgwYN0ldffVVu/2PHjik6OloRERHKzMzUtGnTNHnyZKWkpFRz5QAAwF15uPLgMTExDst//etflZiYqJ07dyokJKRM/6SkJDVv3lxz586VJLVr104ZGRmaNWuWhg4dWh0lAwAAN+c2c25KSkq0evVqFRQUKDw8vNw+O3bsUFRUlENbv379lJGRoaKiouooEwAAuDmXjtxI0oEDBxQeHq5Lly6pXr16Wrdundq3b19u39zcXAUGBjq0BQYGqri4WGfPnlVQUFCZbWw2m2w2m305Pz+/ck8AAAC4FZeHmzZt2mjv3r36+eeflZKSotGjRys9Pf26AcdisTgsG4ZRbvsVCQkJmj59euUWDdyiccm7XF2CSywa083VJQC4g7n8tpSnp6dat26tsLAwJSQkqFOnTnrjjTfK7dukSRPl5uY6tJ0+fVoeHh7y9/cvd5v4+Hjl5eXZP9nZ2ZV+DgAAwH24fOTmWoZhONxGulp4eLg2bNjg0LZ582aFhYWpdu3a5W5jtVpltVorvU4AAOCeXDpyM23aNG3fvl3Hjx/XgQMH9OKLLyotLU0jRoyQdHnUZdSoUfb+EyZM0IkTJxQbG6tDhw5p8eLFWrRokeLi4lx1CgAAwM24dOTm1KlTGjlypHJyclS/fn117NhRmzZtUt++fSVJOTk5ysrKsvdv1aqVNm7cqKlTp+qtt95ScHCw5s2bx2PgAADAzqXhZtGiRTdcn5ycXKatV69e2rNnTxVVBAAAajqXTygGAACoTIQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKi4NNwkJCerWrZt8fHwUEBCg3/72tzp8+PANt0lLS5PFYinz+eabb6qpagAA4M5cGm7S09M1adIk7dy5U1u2bFFxcbGioqJUUFBw020PHz6snJwc++e+++6rhooBAIC783DlwTdt2uSwvGTJEgUEBGj37t3q2bPnDbcNCAhQgwYNqrA6AABQE7nVnJu8vDxJkp+f3037dunSRUFBQYqMjNS2bduqujQAAFBDuHTk5mqGYSg2NlY9evRQhw4drtsvKChICxcuVGhoqGw2m5YvX67IyEilpaWVO9pjs9lks9nsy/n5+VVSPwAAcA9uE26eeeYZ7d+/X5999tkN+7Vp00Zt2rSxL4eHhys7O1uzZs0qN9wkJCRo+vTplV4vAABwT25xW+rZZ5/V+vXrtW3bNjVt2rTC2z/00EM6evRouevi4+OVl5dn/2RnZ99uuQAAwI05NXJz7NgxtWrV6rYPbhiGnn32Wa1bt05paWlO7zMzM1NBQUHlrrNarbJarbdTJgAAqEGcCjetW7dWz549NW7cOA0bNkxeXl5OHXzSpElauXKlPvzwQ/n4+Cg3N1eSVL9+fdWpU0fS5ZGXkydPatmyZZKkuXPnqmXLlgoJCVFhYaFWrFihlJQUpaSkOFUDAAAwF6duS+3bt09dunTR888/ryZNmugPf/iDvvzyywrvJzExUXl5eerdu7eCgoLsnzVr1tj75OTkKCsry75cWFiouLg4dezYUREREfrss8+UmpqqIUOGOHMqAADAZCyGYRjOblxcXKwNGzYoOTlZH330ke677z6NGzdOI0eOVOPGjSuzzkqTn5+v+vXrKy8vT76+vtV+/HHJu6r9mO5g0Zhuri7BJbjeAFA5KvL3+7YmFHt4eGjw4MFau3atZs6cqe+++05xcXFq2rSpRo0apZycnNvZPQAAQIXdVrjJyMjQxIkTFRQUpDlz5iguLk7fffedtm7dqpMnT2rQoEGVVScAAMAtcWpC8Zw5c7RkyRIdPnxY0dHRWrZsmaKjo3XXXZezUqtWrbRgwQK1bdu2UosFAAC4GafCTWJiosaOHaunnnpKTZo0KbdP8+bNtWjRotsqDgAAoKKcCjfX+8G8q3l6emr06NHO7B4AAMBpTs25WbJkid57770y7e+9956WLl1620UBAAA4y6lwM2PGDDVq1KhMe0BAgF577bXbLgoAAMBZToWbEydOlPuqhBYtWjj84B4AAEB1cyrcBAQEaP/+/WXa9+3bJ39//9suCgAAwFlOhZvhw4dr8uTJ2rZtm0pKSlRSUqKtW7dqypQpGj58eGXXCAAAcMucelrq1Vdf1YkTJxQZGSkPj8u7KC0t1ahRo5hzAwAAXMqpcOPp6ak1a9bolVde0b59+1SnTh098MADatGiRWXXBwAAUCFOhZsr7r//ft1///2VVQsAAMBtcyrclJSUKDk5WZ988olOnz6t0tJSh/Vbt26tlOIAAAAqyqlwM2XKFCUnJ2vAgAHq0KGDLBZLZdcFAADgFKfCzerVq7V27VpFR0dXdj0AAAC3xalHwT09PdW6devKrgUAAOC2ORVunn/+eb3xxhsyDKOy6wEAALgtTt2W+uyzz7Rt2zZ99NFHCgkJUe3atR3Wf/DBB5VSHAAAQEU5FW4aNGigwYMHV3YtAAAAt82pcLNkyZLKrgMAAKBSODXnRpKKi4v18ccfa8GCBbpw4YIk6ccff9TFixcrrTgAAICKcmrk5sSJE3rssceUlZUlm82mvn37ysfHR3/729906dIlJSUlVXadAAAAt8SpkZspU6YoLCxMP/30k+rUqWNvHzx4sD755JNKKw4AAKCinH5a6vPPP5enp6dDe4sWLXTy5MlKKQwAAMAZTo3clJaWqqSkpEz7Dz/8IB8fn9suCgAAwFlOhZu+fftq7ty59mWLxaKLFy/qpZde4pUMAADApZy6LfX3v/9dffr0Ufv27XXp0iU98cQTOnr0qBo1aqRVq1ZVdo0AAAC3zKlwExwcrL1792rVqlXas2ePSktLNW7cOI0YMcJhgjEAAEB1cyrcSFKdOnU0duxYjR07tjLrAQAAuC1OhZtly5bdcP2oUaOcKgYAAOB2ORVupkyZ4rBcVFSkX375RZ6envL29ibcAAAAl3HqaamffvrJ4XPx4kUdPnxYPXr0YEIxAABwKaffLXWt++67TzNmzCgzqnMjCQkJ6tatm3x8fBQQEKDf/va3Onz48E23S09PV2hoqLy8vHTPPffwugcAAGBXaeFGkmrVqqUff/zxlvunp6dr0qRJ2rlzp7Zs2aLi4mJFRUWpoKDgutscO3ZM0dHRioiIUGZmpqZNm6bJkycrJSWlMk4BAADUcE7NuVm/fr3DsmEYysnJ0fz58/Xwww/f8n42bdrksLxkyRIFBARo9+7d6tmzZ7nbJCUlqXnz5vYfEWzXrp0yMjI0a9YsDR06tGInAgAATMepcPPb3/7WYdlisahx48Z65JFHNHv2bKeLycvLkyT5+fldt8+OHTsUFRXl0NavXz8tWrRIRUVFql27ttPHBwAANZ9T4aa0tLSy65BhGIqNjVWPHj3UoUOH6/bLzc1VYGCgQ1tgYKCKi4t19uxZBQUFOayz2Wyy2Wz25fz8/MotHAAAuBWnf8Svsj3zzDPav3+/Pvvss5v2tVgsDsuGYZTbLl2etDx9+vTKKRIAbmBc8i5Xl+ASi8Z0c3UJgAOnwk1sbOwt950zZ85N+zz77LNav369Pv30UzVt2vSGfZs0aaLc3FyHttOnT8vDw0P+/v5l+sfHxzvUm5+fr2bNmt1i9QAAoKZxKtxkZmZqz549Ki4uVps2bSRJR44cUa1atdS1a1d7v/JGUq5mGIaeffZZrVu3TmlpaWrVqtVNjx0eHq4NGzY4tG3evFlhYWHlzrexWq2yWq23cloAAMAEnAo3MTEx8vHx0dKlS9WwYUNJl3/Y76mnnlJERISef/75W9rPpEmTtHLlSn344Yfy8fGxj8jUr1/f/gLO+Ph4nTx50v7KhwkTJmj+/PmKjY3V008/rR07dmjRokX8eCAAAJDk5O/czJ49WwkJCfZgI0kNGzbUq6++WqGnpRITE5WXl6fevXsrKCjI/lmzZo29T05OjrKysuzLrVq10saNG5WWlqbOnTvrlVde0bx583gMHAAASHJy5CY/P1+nTp1SSEiIQ/vp06d14cKFW97PlYnAN5KcnFymrVevXtqzZ88tHwcAANw5nBq5GTx4sJ566im9//77+uGHH/TDDz/o/fff17hx4zRkyJDKrhEAAOCWOTVyk5SUpLi4OD355JMqKiq6vCMPD40bN06vv/56pRYIAABQEU6FG29vb7399tt6/fXX9d1338kwDLVu3Vp169at7PoAAAAq5LZenJmTk6OcnBzdf//9qlu37i3NoQEAAKhKToWbc+fOKTIyUvfff7+io6OVk5MjSRo/fvwtPwYOAABQFZwKN1OnTlXt2rWVlZUlb29ve/vjjz9e5k3fAAAA1cmpOTebN2/Wv//97zKvSrjvvvt04sSJSikMAADAGU6N3BQUFDiM2Fxx9uxZXnUAAABcyqlw07NnT/vrEKTL75AqLS3V66+/rj59+lRacQAAABXl1G2p119/Xb1791ZGRoYKCwv1wgsv6KuvvtL58+f1+eefV3aNAAAAt8ypkZv27dtr//79+s1vfqO+ffuqoKBAQ4YMUWZmpu69997KrhEAAOCWVXjkpqioSFFRUVqwYIGmT59eFTUBAAA4rcIjN7Vr19bBgwdlsViqoh4AAIDb4tRtqVGjRmnRokWVXQsAAMBtc2pCcWFhod555x1t2bJFYWFhZd4pNWfOnEopDgAAoKIqFG6+//57tWzZUgcPHlTXrl0lSUeOHHHow+0qAADgShUKN/fdd59ycnK0bds2SZdftzBv3jwFBgZWSXEAAAAVVaE5N9e+9fujjz5SQUFBpRYEAABwO5yaUHzFtWEHAADA1SoUbiwWS5k5NcyxAQAA7qRCc24Mw9CYMWPsL8e8dOmSJkyYUOZpqQ8++KDyKgQAAKiACoWb0aNHOyw/+eSTlVoMAADA7apQuFmyZElV1QEAAFApbmtCMQAAgLsh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFNxabj59NNPFRMTo+DgYFksFv3zn/+8Yf+0tDT7m8mv/nzzzTfVUzAAAHB7FXq3VGUrKChQp06d9NRTT2no0KG3vN3hw4fl6+trX27cuHFVlAcAAGogl4ab/v37q3///hXeLiAgQA0aNKj8ggAAQI1XI+fcdOnSRUFBQYqMjNS2bdtcXQ4AAHAjLh25qaigoCAtXLhQoaGhstlsWr58uSIjI5WWlqaePXuWu43NZpPNZrMv5+fnV1e5AADABWpUuGnTpo3atGljXw4PD1d2drZmzZp13XCTkJCg6dOnV1eJAADAxWrkbamrPfTQQzp69Oh118fHxysvL8/+yc7OrsbqAABAdatRIzflyczMVFBQ0HXXW61WWa3WaqwIAAC4kkvDzcWLF/Xtt9/al48dO6a9e/fKz89PzZs3V3x8vE6ePKlly5ZJkubOnauWLVsqJCREhYWFWrFihVJSUpSSkuKqUwAAAG7GpeEmIyNDffr0sS/HxsZKkkaPHq3k5GTl5OQoKyvLvr6wsFBxcXE6efKk6tSpo5CQEKWmpio6OrraawcAAO7JpeGmd+/eMgzjuuuTk5Mdll944QW98MILVVwVAACoyWr8hGIAAICrEW4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpuDTcfPrpp4qJiVFwcLAsFov++c9/3nSb9PR0hYaGysvLS/fcc4+SkpKqvlAAAFBjuDTcFBQUqFOnTpo/f/4t9T927Jiio6MVERGhzMxMTZs2TZMnT1ZKSkoVVwoAAGoKD1cevH///urfv/8t909KSlLz5s01d+5cSVK7du2UkZGhWbNmaejQoVVUJQAAqElq1JybHTt2KCoqyqGtX79+ysjIUFFRkYuqAgAA7sSlIzcVlZubq8DAQIe2wMBAFRcX6+zZswoKCiqzjc1mk81msy/n5+dXeZ0AAMB1alS4kSSLxeKwbBhGue1XJCQkaPr06VVeFwDgzjIueZerS3CJRWO6ubqEm6pRt6WaNGmi3Nxch7bTp0/Lw8ND/v7+5W4THx+vvLw8+yc7O7s6SgUAAC5So0ZuwsPDtWHDBoe2zZs3KywsTLVr1y53G6vVKqvVWh3lAQAAN+DSkZuLFy9q79692rt3r6TLj3rv3btXWVlZki6PuowaNcref8KECTpx4oRiY2N16NAhLV68WIsWLVJcXJwrygcAAG7IpSM3GRkZ6tOnj305NjZWkjR69GglJycrJyfHHnQkqVWrVtq4caOmTp2qt956S8HBwZo3bx6PgQMAADuXhpvevXvbJwSXJzk5uUxbr169tGfPniqsCgAA1GQ1akIxAADAzRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqbg83Lz99ttq1aqVvLy8FBoaqu3bt1+3b1pamiwWS5nPN998U40VAwAAd+bScLNmzRo999xzevHFF5WZmamIiAj1799fWVlZN9zu8OHDysnJsX/uu+++aqoYAAC4O5eGmzlz5mjcuHEaP3682rVrp7lz56pZs2ZKTEy84XYBAQFq0qSJ/VOrVq1qqhgAALg7l4WbwsJC7d69W1FRUQ7tUVFR+uKLL264bZcuXRQUFKTIyEht27atKssEAAA1jIerDnz27FmVlJQoMDDQoT0wMFC5ubnlbhMUFKSFCxcqNDRUNptNy5cvV2RkpNLS0tSzZ89yt7HZbLLZbPbl/Pz8yjsJAADgdlwWbq6wWCwOy4ZhlGm7ok2bNmrTpo19OTw8XNnZ2Zo1a9Z1w01CQoKmT59eeQUDAAC35rLbUo0aNVKtWrXKjNKcPn26zGjOjTz00EM6evToddfHx8crLy/P/snOzna6ZgAA4P5cFm48PT0VGhqqLVu2OLRv2bJF3bt3v+X9ZGZmKigo6LrrrVarfH19HT4AAMC8XHpbKjY2ViNHjlRYWJjCw8O1cOFCZWVlacKECZIuj7qcPHlSy5YtkyTNnTtXLVu2VEhIiAoLC7VixQqlpKQoJSXFlacBAADciEvDzeOPP65z587pL3/5i3JyctShQwdt3LhRLVq0kCTl5OQ4/OZNYWGh4uLidPLkSdWpU0chISFKTU1VdHS0q04BAAC4GZdPKJ44caImTpxY7rrk5GSH5RdeeEEvvPBCNVQFAABqKpe/fgEAAKAyEW4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpuDzcvP3222rVqpW8vLwUGhqq7du337B/enq6QkND5eXlpXvuuUdJSUnVVCkAAKgJXBpu1qxZo+eee04vvviiMjMzFRERof79+ysrK6vc/seOHVN0dLQiIiKUmZmpadOmafLkyUpJSanmygEAgLtyabiZM2eOxo0bp/Hjx6tdu3aaO3eumjVrpsTExHL7JyUlqXnz5po7d67atWun8ePHa+zYsZo1a1Y1Vw4AANyVy8JNYWGhdu/eraioKIf2qKgoffHFF+Vus2PHjjL9+/Xrp4yMDBUVFVVZrQAAoObwcNWBz549q5KSEgUGBjq0BwYGKjc3t9xtcnNzy+1fXFyss2fPKigoqMw2NptNNpvNvpyXlydJys/Pv91TcErhrxddclxXc9X37Wpc7zsL1/vOwvV2zXENw7hpX5eFmyssFovDsmEYZdpu1r+89isSEhI0ffr0Mu3NmjWraKm4DSsmuroCVCeu952F631ncfX1vnDhgurXr3/DPi4LN40aNVKtWrXKjNKcPn26zOjMFU2aNCm3v4eHh/z9/cvdJj4+XrGxsfbl0tJSnT9/Xv7+/jcMUWaTn5+vZs2aKTs7W76+vq4uB1WM631n4XrfWe7U620Yhi5cuKDg4OCb9nVZuPH09FRoaKi2bNmiwYMH29u3bNmiQYMGlbtNeHi4NmzY4NC2efNmhYWFqXbt2uVuY7VaZbVaHdoaNGhwe8XXYL6+vnfUvwx3Oq73nYXrfWe5E6/3zUZsrnDp01KxsbF65513tHjxYh06dEhTp05VVlaWJkyYIOnyqMuoUaPs/SdMmKATJ04oNjZWhw4d0uLFi7Vo0SLFxcW56hQAAICbcemcm8cff1znzp3TX/7yF+Xk5KhDhw7auHGjWrRoIUnKyclx+M2bVq1aaePGjZo6dareeustBQcHa968eRo6dKirTgEAALgZl08onjhxoiZOLH92UnJycpm2Xr16ac+ePVVclflYrVa99NJLZW7RwZy43ncWrvedhet9cxbjVp6pAgAAqCFc/m4pAACAykS4AQAApkK4AQAApkK4AQAApkK4uYMwdxwAcCcg3NxBrFarDh065OoyAACoUi7/nRtUvqvfpXW1kpISzZgxw/4erjlz5lRnWahCv/76q3bv3i0/Pz+1b9/eYd2lS5e0du1ah1/7Rs2VmZmpBg0aqFWrVpKkFStWKDExUVlZWWrRooWeeeYZDR8+3MVVoqr89NNPWrp0qY4ePaqgoCCNHj2aF0GXg9+5MaG77rpLnTp1KvMOrfT0dIWFhalu3bqyWCzaunWrawpEpTpy5IiioqKUlZUli8WiiIgIrVq1SkFBQZKkU6dOKTg4WCUlJS6uFJWha9eumj17tvr06aN33nlHkydP1tNPP6127drp8OHDeuedd/TGG29o7Nixri4VlSA4OFgHDhyQv7+/jh07pu7du0uSHnjgAR06dEgXLlzQzp071bZtWxdX6l4INyaUkJCgf/zjH3rnnXf0yCOP2Ntr166tffv2lfk/e9RsgwcPVnFxsZYsWaKff/5ZsbGxOnjwoNLS0tS8eXPCjcnUrVtXhw4dUvPmzdW1a1dNmDBBv//97+3rV65cqb/+9a/66quvXFglKstdd92l3NxcBQQE6He/+51yc3OVmpoqb29v2Ww2DRs2TF5eXnrvvfdcXapbYc6NCcXHx2vNmjX64x//qLi4OBUVFbm6JFShL774Qq+99poaNWqk1q1ba/369erfv78iIiL0/fffu7o8VLI6derozJkzkqSTJ0/qwQcfdFj/4IMP6tixY64oDVXs//7v//TnP/9Z3t7eki7Po/zTn/6knTt3urgy90O4Malu3bpp9+7dOnPmjMLCwnTgwAFZLBZXl4Uq8Ouvv8rDw3H63FtvvaWBAweqV69eOnLkiIsqQ1Xo37+/EhMTJV1+197777/vsH7t2rVq3bq1K0pDFbny326bzabAwECHdYGBgfawi/+PCcUmVq9ePS1dulSrV69W3759uS1hUm3btlVGRobatWvn0P7mm2/KMAwNHDjQRZWhKsycOVMPP/ywevXqpbCwMM2ePVtpaWn2OTc7d+7UunXrXF0mKlFkZKQ8PDyUn5+vI0eOKCQkxL4uKytLjRo1cmF17olwcwcYPny4evTood27d6tFixauLgeVbPDgwVq1apVGjhxZZt38+fNVWlqqpKQkF1SGqhAcHKzMzEzNmDFDGzZskGEY+vLLL5Wdna2HH35Yn3/+ucLCwlxdJirJSy+95LB85ZbUFRs2bFBERER1llQjMKEYAACYCnNuAACAqRBuAACAqRBuAACAqRBuAJhGy5YtNXfuXFeXAcDFCDcAAMBUCDcAAMBUCDcA3MKCBQt09913q7S01KF94MCBGj16tL777jsNGjRIgYGBqlevnrp166aPP/74uvs7fvy4LBaL9u7da2/7+eefZbFYlJaWZm/7+uuvFR0drXr16ikwMFAjR47U2bNnK/v0AFQjwg0At/Cf//mfOnv2rLZt22Zv++mnn/Tvf/9bI0aM0MWLFxUdHa2PP/5YmZmZ6tevn2JiYpSVleX0MXNyctSrVy917txZGRkZ2rRpk06dOqX/+q//qoxTAuAi/EIxALfg5+enxx57TCtXrlRkZKQk6b333pOfn58iIyNVq1YtderUyd7/1Vdf1bp167R+/Xo988wzTh0zMTFRXbt21WuvvWZvW7x4sZo1a6YjR47o/vvvv72TAuASjNwAcBsjRoxQSkqKbDabJOndd9/V8OHDVatWLRUUFOiFF15Q+/bt1aBBA9WrV0/ffPPNbY3c7N69W9u2bVO9evXsn7Zt20qSvvvuu0o5JwDVj5EbAG4jJiZGpaWlSk1NVbdu3bR9+3bNmTNHkvTf//3f+ve//61Zs2apdevWqlOnjoYNG6bCwsJy93XXXZf/3+3qN8wUFRU59CktLVVMTIxmzpxZZvugoKDKOi0A1YxwA8Bt1KlTR0OGDNG7776rb7/9Vvfff79CQ0MlSdu3b9eYMWM0ePBgSdLFixd1/Pjx6+6rcePGki7Pq+nSpYskOUwulqSuXbsqJSVFLVu2lIcH/zkEzILbUgDcyogRI5SamqrFixfrySeftLe3bt1aH3zwgfbu3at9+/bpiSeeKPNk1dXq1Kmjhx56SDNmzNDXX3+tTz/9VH/6058c+kyaNEnnz5/X7373O3355Zf6/vvvtXnzZo0dO1YlJSVVdo4AqhbhBoBbeeSRR+Tn56fDhw/riSeesLf//e9/V8OGDdW9e3fFxMSoX79+6tq16w33tXjxYhUVFSksLExTpkzRq6++6rA+ODhYn3/+uUpKStSvXz916NBBU6ZMUf369e23tQDUPBbj6hvSAAAANRz/awIAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEzl/wFwxbjt1Il4RgAAAABJRU5ErkJggg==",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# plotting FD\n",
+ "frequency_distribution.plot(kind='bar', alpha=0.7)\n",
+ "plt.xlabel('value')\n",
+ "plt.ylabel('Frequency')\n",
+ "plt.title('Frequency Distribution')\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : Horizontal and vertical axes are interchanged. "
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Challenge 2\n",
+ "Now, using the dice results obtained in *challenge 1*, your are going to define some functions that will help you calculate the mean of your data in two different ways, the median and the four quartiles. \n",
+ "\n",
+ "#### 1.- Define a function that computes the mean by summing all the observations and dividing by the total number of observations. You are not allowed to use any methods or functions that directly calculate the mean value. "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 8,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "def calculate_mean(input_series):\n",
+ " count = 0\n",
+ " \n",
+ " for key,value in input_series.items():\n",
+ " count += key*value\n",
+ " return count/sum(input_series.values)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 2.- First, calculate the frequency distribution. Then, calculate the mean using the values of the frequency distribution you've just computed. You are not allowed to use any methods or functions that directly calculate the mean value. "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "4 4\n",
+ "2 3\n",
+ "6 2\n",
+ "5 1\n",
+ "Name: value, dtype: int64"
+ ]
+ },
+ "execution_count": 9,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "# calculating FD\n",
+ "frequency_distribution = result['value'].value_counts()\n",
+ "frequency_distribution"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "3.9"
+ ]
+ },
+ "execution_count": 10,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_mean(frequency_distribution)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 3.- Define a function to calculate the median. You are not allowed to use any methods or functions that directly calculate the median value. \n",
+ "**Hint**: you might need to define two computation cases depending on the number of observations used to calculate the median."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 11,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "def calculate_median(input_series):\n",
+ "\n",
+ " values_list = input_series.index.tolist()\n",
+ " values_list.sort()\n",
+ "\n",
+ " n = len(values_list)\n",
+ " if n % 2 == 0:\n",
+ " middle1 = values_list[n // 2 - 1]\n",
+ " middle2 = values_list[n // 2]\n",
+ " median = (middle1 + middle2) / 2\n",
+ " else:\n",
+ " median = values_list[n // 2]\n",
+ "\n",
+ " return median"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "4.5"
+ ]
+ },
+ "execution_count": 12,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_median(frequency_distribution)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 4.- Define a function to calculate the four quartiles. You can use the function you defined above to compute the median but you are not allowed to use any methods or functions that directly calculate the quartiles. "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 13,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "def calculate_quartiles(df, column):\n",
+ " \n",
+ " values_list = df[column].to_list()\n",
+ " values_list.sort()\n",
+ " \n",
+ " q75, q50, q25 = np.percentile(values_list, [75, 50, 25])\n",
+ " return q75, q50, q25"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "(4.75, 4.0, 2.5)"
+ ]
+ },
+ "execution_count": 14,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_quartiles(result, \"value\")"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Challenge 3\n",
+ "Read the csv `roll_the_dice_hundred.csv` from the `data` folder.\n",
+ "#### 1.- Sort the values and plot them. What do you see?"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 15,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " Unnamed: 0 \n",
+ " roll \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 0 \n",
+ " 0 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 1 \n",
+ " 1 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 2 \n",
+ " 2 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 3 \n",
+ " 3 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 4 \n",
+ " 4 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " Unnamed: 0 \n",
+ " roll \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 0 \n",
+ " 0 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 47 \n",
+ " 47 \n",
+ " 47 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 56 \n",
+ " 56 \n",
+ " 56 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 9 \n",
+ " 9 \n",
+ " 9 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 73 \n",
+ " 73 \n",
+ " 73 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 17 \n",
+ " 17 \n",
+ " 17 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 11 \n",
+ " 11 \n",
+ " 11 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 24 \n",
+ " 24 \n",
+ " 24 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 21 \n",
+ " 21 \n",
+ " 21 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 99 \n",
+ " 99 \n",
+ " 99 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
100 rows × 3 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# plotting the value\n",
+ "dice_hundred[\"value\"].plot(kind='bar', alpha=0.7)\n",
+ "plt.xlabel('roll(index)')\n",
+ "plt.ylabel('value')\n",
+ "plt.title('')\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : The number of times that each value appears (the frequency) differs by value (ex. 5 less vs. 4 and 6)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 2.- Using the functions you defined in *challenge 2*, calculate the mean value of the hundred dice rolls."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 18,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "3.74"
+ ]
+ },
+ "execution_count": 18,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "frequency_distribution = dice_hundred['value'].value_counts()\n",
+ "calculate_mean(frequency_distribution)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 3.- Now, calculate the frequency distribution.\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 19,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "6 23\n",
+ "4 22\n",
+ "2 17\n",
+ "3 14\n",
+ "1 12\n",
+ "5 12\n",
+ "Name: value, dtype: int64"
+ ]
+ },
+ "execution_count": 19,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "# already calculated in above cell\n",
+ "frequency_distribution"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 4.- Plot the histogram. What do you see (shape, values...) ? How can you connect the mean value to the histogram? "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 20,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAASIAAAEiCAYAAABdvt+2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAASuUlEQVR4nO3df2zU9R3H8de37e0rdXeYgv0Vuq6amrnVX6MOB05KtlYbY4L4K2u2gZlOY2HDi2EwQjymtgvJCFmITdgfqH808I8ocWh7ZlpUgkI3o6JTjFU7pGMg65XWHAf33R+kN2u7tle+X970fD6S/nHffve59312ffrttT0cz/M8AYChPOsBAIAQATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwFyB9QBflU6n9dlnnykcDstxHOtxAEyR53kaGBhQeXm58vLGv+Y570L02WefqaKiwnoMAD7p7e3VnDlzxj3nvAtROByWdGb4SCRiPI3/UqmUOjs71dDQoFAoZD1OzmF/g5XN/iYSCVVUVGS+psdz3oVo+NuxSCSSsyEqLCxUJBLhCyUA7G+wprK/k3mJhRerAZgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYO+9+sxqAP769+i++r+nme9rwA9+X5YoIgD1CBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5frMaOakm1qHkaX//OaqP/3Czr+vhf7giAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIC5rELU2tqqa6+9VuFwWMXFxVq8eLHef//9Eed4nqdYLKby8nLNmDFDdXV1OnDggK9DA8gtWYWoq6tLzc3N2rt3r+LxuE6dOqWGhgYNDg5mztmwYYM2btyozZs3a9++fSotLVV9fb0GBgZ8Hx5AbsjqHRpfeOGFEbe3bt2q4uJidXd364YbbpDnedq0aZPWrl2rJUuWSJKefPJJlZSUqL29Xffdd59/kwPIGWf1VrH9/f2SpKKiIklST0+P+vr61NDQkDnHdV0tXLhQe/bsGTNEyWRSyWQyczuRSEiSUqmUUqnU2Yx3Xhp+TLn42M4Hw/vq5nmBrT1duPn+78Hwvk5mL7LZrymHyPM8RaNRXX/99aqpqZEk9fX1SZJKSkpGnFtSUqJPPvlkzHVaW1u1fv36Ucc7OztVWFg41fHOe/F43HqEnPZIbdr3NXft2uX7mkHa8IPg1p7M83doaGjS6005RMuXL9dbb72lV199ddTnHGfkm5Z7njfq2LA1a9YoGo1mbicSCVVUVKihoUGRSGSq4523UqmU4vG46uvrFQqFrMfJOcP7u25/npJpf988/53Yjb6uF7SaWIfva7p5nh6pTU/q+Tv83c1kTClEK1as0M6dO7V7927NmTMnc7y0tFTSmSujsrKyzPEjR46Mukoa5rquXNcddTwUCuX0F2quPz5rybTj+7/iMd3+//L78X/ZZJ6/2exXVj818zxPy5cv19NPP62//vWvqqqqGvH5qqoqlZaWjrhsO3nypLq6ujR//vxs7grA10hWV0TNzc1qb2/Xs88+q3A4nHlNaObMmZoxY4Ycx9HKlSvV0tKi6upqVVdXq6WlRYWFhWpqagrkAQCY/rIKUVtbmySprq5uxPGtW7dq2bJlkqRVq1bpiy++0AMPPKDjx49r3rx56uzsVDgc9mVgALknqxB53sQ/DnQcR7FYTLFYbKozAfia4W/NAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMJfVP7AI/9TEOpQ87fi65sd/uNnX9YBzhSsiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJjLOkS7d+/WLbfcovLycjmOo2eeeWbE55ctWybHcUZ8XHfddX7NCyAHZR2iwcFBXXXVVdq8efP/Peemm27S4cOHMx+7du06qyEB5Las36GxsbFRjY2N457juq5KS0unPBSAr5dA3ir25ZdfVnFxsS666CItXLhQjz32mIqLi8c8N5lMKplMZm4nEglJUiqVUiqVCmI8U8OPyc3zAlv764z9/R833/89GN7XyexFNvvleJ435Wkdx9GOHTu0ePHizLHt27frm9/8piorK9XT06N169bp1KlT6u7uluu6o9aIxWJav379qOPt7e0qLCyc6mgAjA0NDampqUn9/f2KRCLjnut7iL7q8OHDqqys1LZt27RkyZJRnx/riqiiokJHjx6dcPjpKJVKKR6Pa93+PCXT/r55/juxG31dbzpif/+nJtbh+5punqdHatOqr69XKBQa99xEIqHZs2dPKkSB/yseZWVlqqys1MGDB8f8vOu6Y14phUKhCR/odJZMO77/Kx65vF/ZYn/l++P/ssl8fWazX4H/HtGxY8fU29ursrKyoO8KwDSV9RXRiRMn9OGHH2Zu9/T06M0331RRUZGKiooUi8V02223qaysTB9//LF+97vfafbs2br11lt9HRxA7sg6RPv379eiRYsyt6PRqCRp6dKlamtr09tvv62nnnpK//nPf1RWVqZFixZp+/btCofD/k0NIKdkHaK6ujqN9/p2R4f/L5AByG38rRkAc4QIgDlCBMAcIQJgLvBfaMT09+3Vfwlk3Y//cHMg62L64YoIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmBu2r9ndVDvpyzxnso4N4J8Dk8XXBEBMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwFzWIdq9e7duueUWlZeXy3EcPfPMMyM+73meYrGYysvLNWPGDNXV1enAgQN+zQsgB2UdosHBQV111VXavHnzmJ/fsGGDNm7cqM2bN2vfvn0qLS1VfX29BgYGznpYALkp63dobGxsVGNj45if8zxPmzZt0tq1a7VkyRJJ0pNPPqmSkhK1t7frvvvuO7tpAeQkX98qtqenR319fWpoaMgcc11XCxcu1J49e8YMUTKZVDKZzNxOJBKSpFQqpVQqNeF9uvmeD5OPbTL3P9U13Tz/5w5iXim4PWZ/zwjyOey34X2dzF5ks1+O53lT3gXHcbRjxw4tXrxYkrRnzx4tWLBAhw4dUnl5eea8X/3qV/rkk0/U0dExao1YLKb169ePOt7e3q7CwsKpjgbA2NDQkJqamtTf369IJDLuuYG8eb7jOCNue5436tiwNWvWKBqNZm4nEglVVFSooaFhwuElqSY2Om5+eSd2o+9rplIpxeNxrdufp2R67D2ZqiDmlYLbY/b3jCCfw35z8zw9UptWfX29QqHQuOcOf3czGb6GqLS0VJLU19ensrKyzPEjR46opKRkzP+N67pyXXfU8VAoNOEDlaTkaX+fbF+dISjJtOP77EHNG9Qes79nBPkcDspkvj6z2S9ff4+oqqpKpaWlisfjmWMnT55UV1eX5s+f7+ddAcghWV8RnThxQh9++GHmdk9Pj958800VFRXpW9/6llauXKmWlhZVV1erurpaLS0tKiwsVFNTk6+DA8gdWYdo//79WrRoUeb28Os7S5cu1RNPPKFVq1bpiy++0AMPPKDjx49r3rx56uzsVDgc9m9qADkl6xDV1dVpvB+0OY6jWCymWCx2NnMB+Brhb80AmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATDne4hisZgcxxnxUVpa6vfdAMghBUEs+r3vfU8vvvhi5nZ+fn4QdwMgRwQSooKCAq6CAExaICE6ePCgysvL5bqu5s2bp5aWFl1yySVjnptMJpVMJjO3E4mEJCmVSimVSk14X26+58/QY5jM/U91TTfP/7mDmFcKbo/Z3zOCfA77bXhfJ7MX2eyX43mer7vw/PPPa2hoSJdddpn+9a9/6dFHH9U//vEPHThwQLNmzRp1fiwW0/r160cdb29vV2FhoZ+jATiHhoaG1NTUpP7+fkUikXHP9T1EXzU4OKhLL71Uq1atUjQaHfX5sa6IKioqdPTo0QmHl6SaWIev837ZO7EbfV8zlUopHo9r3f48JdOOr2sHMa8U3B6zv2cE+Rz2m5vn6ZHatOrr6xUKhcY9N5FIaPbs2ZMKUSDfmn3ZhRdeqCuuuEIHDx4c8/Ou68p13VHHQ6HQhA9UkpKn/X2yfXWGoCTTju+zBzVvUHvM/p4R5HM4KJP5+sxmvwL/PaJkMqn33ntPZWVlQd8VgGnK9xA99NBD6urqUk9Pj15//XXdfvvtSiQSWrp0qd93BSBH+P6t2T//+U/99Kc/1dGjR3XxxRfruuuu0969e1VZWen3XQHIEb6HaNu2bX4vCSDH8bdmAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYCyxEjz/+uKqqqnTBBRdo7ty5euWVV4K6KwDTXCAh2r59u1auXKm1a9fq73//u370ox+psbFRn376aRB3B2CaCyREGzdu1C9/+Uvdc889uvzyy7Vp0yZVVFSora0tiLsDMM0V+L3gyZMn1d3drdWrV4843tDQoD179ow6P5lMKplMZm739/dLkj7//HOlUqkJ76/g1OBZTvz/HTt2zPc1U6mUhoaGVJDK0+m04+vaQcwrBbfH7O8ZQT6H/VaQ9jQ0lNaxY8cUCoXGPXdgYECS5HnexAt7Pjt06JAnyXvttddGHH/ssce8yy67bNT5Dz/8sCeJDz74yNGP3t7eCbvh+xXRMMcZ+V8jz/NGHZOkNWvWKBqNZm6n02l9/vnnmjVr1pjnT3eJREIVFRXq7e1VJBKxHifnsL/BymZ/Pc/TwMCAysvLJ1zX9xDNnj1b+fn56uvrG3H8yJEjKikpGXW+67pyXXfEsYsuusjvsc47kUiEL5QAsb/Bmuz+zpw5c1Lr+f5i9Te+8Q3NnTtX8Xh8xPF4PK758+f7fXcAckAg35pFo1H9/Oc/V21trX74wx9qy5Yt+vTTT3X//fcHcXcAprlAQnTXXXfp2LFj+v3vf6/Dhw+rpqZGu3btUmVlZRB3N624rquHH3541Lej8Af7G6yg9tfxvMn8bA0AgsPfmgEwR4gAmCNEAMwRIgDmCNE50NraqmuvvVbhcFjFxcVavHix3n//feuxclZra6scx9HKlSutR8kphw4d0s9+9jPNmjVLhYWFuvrqq9Xd3e3L2oToHOjq6lJzc7P27t2reDyuU6dOqaGhQYOD0+ePHaeLffv2acuWLbryyiutR8kpx48f14IFCxQKhfT888/r3Xff1R//+Eff/gqCH98b+Pe//63i4mJ1dXXphhtusB4nZ5w4cULf//739fjjj+vRRx/V1VdfrU2bNlmPlRNWr16t1157LbA3OOSKyMDwW50UFRUZT5JbmpubdfPNN+snP/mJ9Sg5Z+fOnaqtrdUdd9yh4uJiXXPNNfrzn//s2/qE6BzzPE/RaFTXX3+9ampqrMfJGdu2bdPf/vY3tba2Wo+Skz766CO1tbWpurpaHR0duv/++/XrX/9aTz31lC/rB/Y2IBjb8uXL9dZbb+nVV1+1HiVn9Pb26je/+Y06Ozt1wQUXWI+Tk9LptGpra9XS0iJJuuaaa3TgwAG1tbXpF7/4xVmvzxXRObRixQrt3LlTL730kubMmWM9Ts7o7u7WkSNHNHfuXBUUFKigoEBdXV3605/+pIKCAp0+fdp6xGmvrKxM3/3ud0ccu/zyy317H3quiM4Bz/O0YsUK7dixQy+//LKqqqqsR8opP/7xj/X222+POHb33XfrO9/5jn77298qPz/faLLcsWDBglG/cvLBBx/49ofshOgcaG5uVnt7u5599lmFw+HMm8bNnDlTM2bMMJ5u+guHw6Neb7vwwgs1a9YsXofzyYMPPqj58+erpaVFd955p9544w1t2bJFW7Zs8WV9fnx/Dvy/t7zdunWrli1bdm6H+Zqoq6vjx/c+e+6557RmzRodPHhQVVVVikajuvfee31ZmxABMMeL1QDMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAuf8CkE34rsQ/UD0AAAAASUVORK5CYII=",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "dice_hundred[\"value\"].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show() "
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : As larger numbers (4 and 6) appears more frequently compared to smaller values (2 and 3), mean is larger than the median (1 and 5 are found with same frequency, meaning they make no difference)."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 5.- Read the `roll_the_dice_thousand.csv` from the `data` folder. Plot the frequency distribution as you did before. Has anything changed? Why do you think it changed?"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 21,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " Unnamed: 0 \n",
+ " roll \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 0 \n",
+ " 0 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 1 \n",
+ " 1 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 2 \n",
+ " 2 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 3 \n",
+ " 3 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 4 \n",
+ " 4 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "dice_thousand[\"value\"].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show() "
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comments : The frequency distribution is more consistent for all values. It becomes more stable as the number of rolls increases."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Challenge 4\n",
+ "In the `data` folder of this repository you will find three different files with the prefix `ages_population`. These files contain information about a poll answered by a thousand people regarding their age. Each file corresponds to the poll answers in different neighbourhoods of Barcelona.\n",
+ "\n",
+ "#### 1.- Read the file `ages_population.csv`. Calculate the frequency distribution and plot it as we did during the lesson. Try to guess the range in which the mean and the standard deviation will be by looking at the plot. "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 23,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " observation \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 68.0 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 12.0 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 45.0 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 38.0 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 49.0 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 995 \n",
+ " 27.0 \n",
+ " \n",
+ " \n",
+ " 996 \n",
+ " 47.0 \n",
+ " \n",
+ " \n",
+ " 997 \n",
+ " 53.0 \n",
+ " \n",
+ " \n",
+ " 998 \n",
+ " 33.0 \n",
+ " \n",
+ " \n",
+ " 999 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
1000 rows × 1 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "population['observation'].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Assumption\n",
+ "- range of mean : 30-40\n",
+ "- std : 15 "
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 2.- Calculate the exact mean and standard deviation and compare them with your guesses. Do they fall inside the ranges you guessed?"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 25,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "36.56"
+ ]
+ },
+ "execution_count": 25,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "np.mean(population['observation'])"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 26,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "12.81008977329979"
+ ]
+ },
+ "execution_count": 26,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "np.std(population['observation'])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Diffrence vs. assumption\n",
+ "- mean: falls inside the range assumed\n",
+ "- std: lower than assumption"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 3.- Now read the file `ages_population2.csv` . Calculate the frequency distribution and plot it."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 27,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " observation \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 25.0 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 29.0 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 29.0 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 995 \n",
+ " 26.0 \n",
+ " \n",
+ " \n",
+ " 996 \n",
+ " 22.0 \n",
+ " \n",
+ " \n",
+ " 997 \n",
+ " 21.0 \n",
+ " \n",
+ " \n",
+ " 998 \n",
+ " 19.0 \n",
+ " \n",
+ " \n",
+ " 999 \n",
+ " 28.0 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
1000 rows × 1 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "population_2['observation'].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 4.- What do you see? Is there any difference with the frequency distribution in step 1?"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : Distribution is centered in the range of 26-29 and std is assumed to be smaller compared to that of step 2 (meaning most of the values are found between smaller range of values)."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 5.- Calculate the mean and standard deviation. Compare the results with the mean and standard deviation in step 2. What do you think?"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 29,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "27.155"
+ ]
+ },
+ "execution_count": 29,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "np.mean(population_2['observation'])"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 30,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "2.9683286543103704"
+ ]
+ },
+ "execution_count": 30,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "np.std(population_2['observation'])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : Same to above."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Challenge 5\n",
+ "Now is the turn of `ages_population3.csv`.\n",
+ "\n",
+ "#### 1.- Read the file `ages_population3.csv`. Calculate the frequency distribution and plot it."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 31,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " observation \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 21.0 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 21.0 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 24.0 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 54.0 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 995 \n",
+ " 16.0 \n",
+ " \n",
+ " \n",
+ " 996 \n",
+ " 55.0 \n",
+ " \n",
+ " \n",
+ " 997 \n",
+ " 30.0 \n",
+ " \n",
+ " \n",
+ " 998 \n",
+ " 35.0 \n",
+ " \n",
+ " \n",
+ " 999 \n",
+ " 43.0 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
1000 rows × 1 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "population_3['observation'].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 2.- Calculate the mean and standard deviation. Compare the results with the plot in step 1. What is happening?"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 33,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "41.989"
+ ]
+ },
+ "execution_count": 33,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "np.mean(population_3['observation'])"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 34,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "16.136631587788084"
+ ]
+ },
+ "execution_count": 34,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "np.std(population_3['observation'])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : Mean is larger than mode and it is assume to be positively skewed. Distribution is varied in wider range resulting in larger std."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 3.- Calculate the four quartiles. Use the results to explain your reasoning for question in step 2. How much of a difference is there between the median and the mean?"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 35,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "(53.0, 40.0, 30.0)"
+ ]
+ },
+ "execution_count": 35,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_quartiles(population_3, \"observation\")"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : mean is slightly larger than median and the difference between median(q2) and q3 is larger than that of q1 and median, meaning it is positive distribution."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "#### 4.- Calculate other percentiles that might be useful to give more arguments to your reasoning."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 36,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "array([22., 67.])"
+ ]
+ },
+ "execution_count": 36,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "values_list = population_3[\"observation\"].to_list()\n",
+ "values_list.sort()\n",
+ " \n",
+ "np.percentile(values_list, [10, 90])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : If we see the 10th and 90th percentile and compared the difference between the median, we find that the difference vs. 90th(67-40=23) is larger than that vs. 10th (40-22=18) , meaning it is positively skewed."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Bonus challenge\n",
+ "Compare the information about the three neighbourhoods. Prepare a report about the three of them. Remember to find out which are their similarities and their differences backing your arguments in basic statistics."
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 37,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "# your code here"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "\"\"\"\n",
+ "your comments here\n",
+ "\"\"\""
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 3 (ipykernel)",
+ "language": "python",
+ "name": "python3"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 3
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython3",
+ "version": "3.11.3"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
diff --git a/your-code/ages_population.csv b/your-code/ages_population.csv
new file mode 100644
index 0000000..64d8a0a
--- /dev/null
+++ b/your-code/ages_population.csv
@@ -0,0 +1,1001 @@
+observation
+68.0
+12.0
+45.0
+38.0
+49.0
+27.0
+39.0
+12.0
+42.0
+33.0
+30.0
+25.0
+30.0
+44.0
+53.0
+46.0
+50.0
+22.0
+6.0
+29.0
+29.0
+27.0
+35.0
+38.0
+28.0
+26.0
+60.0
+41.0
+38.0
+41.0
+44.0
+52.0
+46.0
+39.0
+44.0
+46.0
+32.0
+23.0
+15.0
+40.0
+42.0
+32.0
+45.0
+29.0
+22.0
+41.0
+39.0
+63.0
+39.0
+31.0
+34.0
+28.0
+45.0
+33.0
+32.0
+61.0
+64.0
+37.0
+56.0
+44.0
+33.0
+38.0
+40.0
+38.0
+56.0
+14.0
+52.0
+34.0
+14.0
+34.0
+31.0
+46.0
+50.0
+37.0
+13.0
+12.0
+25.0
+28.0
+51.0
+13.0
+36.0
+52.0
+13.0
+30.0
+36.0
+35.0
+26.0
+34.0
+51.0
+52.0
+35.0
+44.0
+23.0
+29.0
+25.0
+30.0
+27.0
+42.0
+18.0
+39.0
+42.0
+48.0
+30.0
+40.0
+34.0
+28.0
+48.0
+48.0
+42.0
+53.0
+34.0
+37.0
+32.0
+29.0
+18.0
+35.0
+58.0
+37.0
+32.0
+49.0
+20.0
+42.0
+29.0
+22.0
+39.0
+41.0
+54.0
+20.0
+19.0
+39.0
+39.0
+39.0
+58.0
+23.0
+45.0
+13.0
+44.0
+39.0
+30.0
+37.0
+62.0
+45.0
+33.0
+55.0
+33.0
+39.0
+34.0
+32.0
+25.0
+21.0
+39.0
+43.0
+18.0
+40.0
+50.0
+35.0
+46.0
+36.0
+30.0
+44.0
+34.0
+58.0
+14.0
+27.0
+62.0
+42.0
+27.0
+50.0
+29.0
+41.0
+30.0
+37.0
+43.0
+43.0
+49.0
+31.0
+17.0
+42.0
+48.0
+29.0
+38.0
+31.0
+20.0
+50.0
+26.0
+45.0
+56.0
+35.0
+48.0
+35.0
+55.0
+44.0
+23.0
+39.0
+39.0
+45.0
+21.0
+43.0
+38.0
+40.0
+35.0
+25.0
+7.0
+26.0
+40.0
+52.0
+46.0
+47.0
+38.0
+1.0
+16.0
+54.0
+45.0
+35.0
+54.0
+41.0
+12.0
+37.0
+50.0
+37.0
+24.0
+47.0
+39.0
+30.0
+41.0
+31.0
+23.0
+42.0
+24.0
+23.0
+31.0
+45.0
+53.0
+21.0
+59.0
+63.0
+41.0
+53.0
+59.0
+48.0
+57.0
+39.0
+32.0
+36.0
+21.0
+40.0
+51.0
+44.0
+59.0
+25.0
+52.0
+42.0
+8.0
+34.0
+43.0
+46.0
+36.0
+42.0
+46.0
+44.0
+29.0
+17.0
+56.0
+41.0
+30.0
+67.0
+41.0
+34.0
+42.0
+22.0
+19.0
+57.0
+16.0
+24.0
+26.0
+32.0
+35.0
+36.0
+26.0
+21.0
+20.0
+28.0
+19.0
+56.0
+41.0
+38.0
+42.0
+31.0
+30.0
+41.0
+31.0
+33.0
+42.0
+45.0
+62.0
+29.0
+41.0
+19.0
+42.0
+42.0
+54.0
+2.0
+40.0
+28.0
+42.0
+28.0
+26.0
+43.0
+32.0
+54.0
+34.0
+25.0
+32.0
+14.0
+66.0
+24.0
+43.0
+33.0
+29.0
+8.0
+43.0
+30.0
+34.0
+37.0
+40.0
+42.0
+48.0
+19.0
+35.0
+41.0
+26.0
+36.0
+23.0
+46.0
+35.0
+28.0
+39.0
+21.0
+4.0
+73.0
+43.0
+48.0
+20.0
+49.0
+28.0
+26.0
+34.0
+20.0
+39.0
+46.0
+37.0
+20.0
+29.0
+63.0
+36.0
+49.0
+36.0
+36.0
+34.0
+46.0
+44.0
+15.0
+38.0
+82.0
+48.0
+29.0
+49.0
+57.0
+16.0
+12.0
+36.0
+59.0
+49.0
+17.0
+25.0
+33.0
+37.0
+40.0
+43.0
+57.0
+43.0
+38.0
+35.0
+30.0
+14.0
+48.0
+24.0
+32.0
+47.0
+29.0
+50.0
+43.0
+55.0
+36.0
+49.0
+46.0
+45.0
+52.0
+36.0
+30.0
+33.0
+34.0
+18.0
+32.0
+40.0
+37.0
+36.0
+63.0
+44.0
+57.0
+35.0
+28.0
+57.0
+15.0
+40.0
+47.0
+17.0
+53.0
+39.0
+29.0
+47.0
+37.0
+30.0
+19.0
+66.0
+56.0
+8.0
+22.0
+43.0
+39.0
+21.0
+41.0
+54.0
+51.0
+37.0
+23.0
+56.0
+70.0
+39.0
+14.0
+60.0
+26.0
+30.0
+47.0
+52.0
+30.0
+54.0
+5.0
+22.0
+23.0
+34.0
+2.0
+34.0
+45.0
+31.0
+42.0
+47.0
+35.0
+36.0
+39.0
+41.0
+60.0
+42.0
+26.0
+45.0
+25.0
+32.0
+47.0
+36.0
+37.0
+40.0
+57.0
+40.0
+59.0
+31.0
+32.0
+63.0
+38.0
+41.0
+43.0
+17.0
+34.0
+28.0
+43.0
+51.0
+8.0
+30.0
+43.0
+24.0
+1.0
+16.0
+43.0
+27.0
+71.0
+50.0
+50.0
+9.0
+30.0
+15.0
+32.0
+50.0
+39.0
+24.0
+55.0
+38.0
+17.0
+36.0
+43.0
+40.0
+42.0
+37.0
+55.0
+31.0
+31.0
+31.0
+21.0
+43.0
+45.0
+36.0
+23.0
+32.0
+14.0
+16.0
+69.0
+11.0
+33.0
+36.0
+39.0
+19.0
+31.0
+26.0
+52.0
+41.0
+30.0
+31.0
+32.0
+35.0
+43.0
+52.0
+41.0
+17.0
+38.0
+28.0
+54.0
+42.0
+24.0
+39.0
+44.0
+31.0
+48.0
+6.0
+50.0
+27.0
+45.0
+28.0
+29.0
+43.0
+44.0
+14.0
+28.0
+15.0
+37.0
+61.0
+36.0
+35.0
+32.0
+34.0
+26.0
+37.0
+41.0
+59.0
+45.0
+46.0
+32.0
+30.0
+33.0
+27.0
+10.0
+42.0
+30.0
+54.0
+36.0
+44.0
+45.0
+54.0
+37.0
+40.0
+23.0
+66.0
+32.0
+39.0
+34.0
+46.0
+39.0
+38.0
+38.0
+50.0
+55.0
+30.0
+47.0
+41.0
+56.0
+39.0
+41.0
+27.0
+23.0
+30.0
+17.0
+45.0
+29.0
+42.0
+26.0
+22.0
+39.0
+23.0
+48.0
+51.0
+39.0
+52.0
+39.0
+41.0
+56.0
+38.0
+44.0
+33.0
+39.0
+44.0
+29.0
+49.0
+13.0
+53.0
+58.0
+27.0
+23.0
+27.0
+45.0
+12.0
+30.0
+40.0
+30.0
+49.0
+40.0
+32.0
+53.0
+37.0
+31.0
+46.0
+31.0
+22.0
+19.0
+12.0
+49.0
+47.0
+42.0
+38.0
+58.0
+43.0
+38.0
+42.0
+30.0
+45.0
+26.0
+38.0
+26.0
+46.0
+20.0
+29.0
+39.0
+37.0
+42.0
+43.0
+48.0
+55.0
+22.0
+26.0
+30.0
+55.0
+35.0
+51.0
+46.0
+38.0
+65.0
+55.0
+41.0
+30.0
+37.0
+27.0
+36.0
+42.0
+38.0
+20.0
+47.0
+31.0
+34.0
+35.0
+36.0
+5.0
+40.0
+41.0
+27.0
+47.0
+63.0
+35.0
+32.0
+27.0
+50.0
+39.0
+38.0
+36.0
+17.0
+27.0
+45.0
+30.0
+28.0
+43.0
+33.0
+45.0
+20.0
+29.0
+48.0
+20.0
+44.0
+42.0
+60.0
+25.0
+48.0
+41.0
+39.0
+11.0
+35.0
+43.0
+51.0
+24.0
+33.0
+22.0
+21.0
+45.0
+48.0
+56.0
+25.0
+33.0
+36.0
+53.0
+51.0
+42.0
+15.0
+32.0
+22.0
+59.0
+54.0
+32.0
+39.0
+47.0
+17.0
+28.0
+29.0
+56.0
+34.0
+28.0
+44.0
+67.0
+26.0
+25.0
+23.0
+56.0
+52.0
+25.0
+49.0
+26.0
+39.0
+63.0
+48.0
+48.0
+30.0
+37.0
+46.0
+14.0
+15.0
+44.0
+39.0
+40.0
+22.0
+30.0
+39.0
+43.0
+46.0
+35.0
+16.0
+44.0
+18.0
+26.0
+27.0
+37.0
+20.0
+35.0
+50.0
+55.0
+45.0
+41.0
+54.0
+56.0
+46.0
+56.0
+24.0
+33.0
+12.0
+31.0
+12.0
+25.0
+49.0
+52.0
+55.0
+27.0
+49.0
+44.0
+49.0
+40.0
+34.0
+16.0
+19.0
+36.0
+32.0
+67.0
+52.0
+38.0
+33.0
+68.0
+14.0
+43.0
+10.0
+11.0
+39.0
+13.0
+54.0
+41.0
+40.0
+18.0
+16.0
+45.0
+35.0
+39.0
+34.0
+46.0
+53.0
+23.0
+55.0
+37.0
+29.0
+48.0
+35.0
+45.0
+68.0
+29.0
+40.0
+33.0
+64.0
+45.0
+10.0
+47.0
+23.0
+39.0
+20.0
+41.0
+36.0
+41.0
+25.0
+46.0
+49.0
+28.0
+40.0
+58.0
+46.0
+37.0
+41.0
+38.0
+35.0
+31.0
+30.0
+35.0
+34.0
+15.0
+58.0
+41.0
+59.0
+32.0
+12.0
+27.0
+41.0
+37.0
+47.0
+49.0
+47.0
+53.0
+35.0
+38.0
+24.0
+50.0
+27.0
+27.0
+26.0
+24.0
+30.0
+43.0
+35.0
+22.0
+22.0
+32.0
+37.0
+55.0
+27.0
+36.0
+50.0
+21.0
+45.0
+27.0
+53.0
+41.0
+28.0
+45.0
+39.0
+43.0
+29.0
+18.0
+34.0
+19.0
+26.0
+32.0
+65.0
+38.0
+24.0
+27.0
+24.0
+19.0
+49.0
+25.0
+41.0
+22.0
+31.0
+24.0
+21.0
+8.0
+56.0
+42.0
+43.0
+62.0
+22.0
+67.0
+38.0
+43.0
+41.0
+34.0
+38.0
+21.0
+35.0
+28.0
+27.0
+21.0
+40.0
+44.0
+25.0
+40.0
+29.0
+39.0
+25.0
+33.0
+37.0
+35.0
+26.0
+34.0
+43.0
+39.0
+32.0
+33.0
+24.0
+24.0
+25.0
+36.0
+49.0
+31.0
+12.0
+9.0
+35.0
+21.0
+36.0
+37.0
+34.0
+24.0
+35.0
+38.0
+33.0
+42.0
+32.0
+29.0
+36.0
+40.0
+27.0
+47.0
+53.0
+33.0
+31.0
diff --git a/your-code/ages_population2.csv b/your-code/ages_population2.csv
new file mode 100644
index 0000000..00860cb
--- /dev/null
+++ b/your-code/ages_population2.csv
@@ -0,0 +1,1001 @@
+observation
+25.0
+31.0
+29.0
+31.0
+29.0
+29.0
+28.0
+30.0
+24.0
+26.0
+30.0
+29.0
+24.0
+26.0
+30.0
+25.0
+25.0
+30.0
+21.0
+29.0
+25.0
+28.0
+27.0
+27.0
+31.0
+27.0
+26.0
+24.0
+29.0
+23.0
+31.0
+26.0
+28.0
+22.0
+28.0
+31.0
+29.0
+22.0
+28.0
+25.0
+26.0
+32.0
+26.0
+28.0
+29.0
+31.0
+27.0
+25.0
+24.0
+30.0
+27.0
+22.0
+25.0
+28.0
+25.0
+31.0
+27.0
+30.0
+28.0
+27.0
+29.0
+26.0
+28.0
+26.0
+22.0
+30.0
+27.0
+31.0
+27.0
+24.0
+27.0
+28.0
+23.0
+20.0
+23.0
+25.0
+29.0
+33.0
+27.0
+25.0
+25.0
+26.0
+31.0
+28.0
+29.0
+23.0
+30.0
+26.0
+27.0
+24.0
+29.0
+26.0
+28.0
+24.0
+25.0
+21.0
+19.0
+26.0
+28.0
+28.0
+24.0
+28.0
+30.0
+27.0
+29.0
+26.0
+28.0
+31.0
+32.0
+25.0
+30.0
+30.0
+24.0
+22.0
+29.0
+25.0
+21.0
+30.0
+30.0
+29.0
+29.0
+22.0
+33.0
+31.0
+26.0
+25.0
+28.0
+31.0
+28.0
+28.0
+31.0
+23.0
+33.0
+25.0
+33.0
+23.0
+30.0
+27.0
+24.0
+29.0
+24.0
+28.0
+23.0
+24.0
+32.0
+29.0
+31.0
+34.0
+31.0
+27.0
+26.0
+28.0
+26.0
+28.0
+26.0
+33.0
+29.0
+22.0
+26.0
+25.0
+28.0
+29.0
+28.0
+28.0
+31.0
+25.0
+30.0
+26.0
+26.0
+30.0
+27.0
+23.0
+29.0
+30.0
+25.0
+32.0
+27.0
+29.0
+28.0
+27.0
+27.0
+23.0
+23.0
+27.0
+25.0
+27.0
+35.0
+30.0
+27.0
+27.0
+28.0
+30.0
+26.0
+32.0
+21.0
+29.0
+28.0
+31.0
+31.0
+28.0
+27.0
+31.0
+26.0
+24.0
+28.0
+28.0
+26.0
+26.0
+24.0
+28.0
+25.0
+29.0
+26.0
+24.0
+30.0
+25.0
+29.0
+30.0
+29.0
+26.0
+28.0
+25.0
+29.0
+29.0
+31.0
+26.0
+23.0
+25.0
+24.0
+27.0
+23.0
+28.0
+31.0
+26.0
+25.0
+31.0
+29.0
+24.0
+28.0
+26.0
+27.0
+25.0
+27.0
+31.0
+29.0
+29.0
+29.0
+32.0
+24.0
+30.0
+28.0
+27.0
+22.0
+29.0
+21.0
+25.0
+27.0
+24.0
+23.0
+33.0
+23.0
+30.0
+25.0
+35.0
+26.0
+23.0
+27.0
+30.0
+27.0
+28.0
+24.0
+27.0
+20.0
+29.0
+31.0
+23.0
+29.0
+29.0
+31.0
+25.0
+24.0
+23.0
+23.0
+20.0
+28.0
+28.0
+28.0
+27.0
+36.0
+20.0
+31.0
+28.0
+29.0
+30.0
+33.0
+28.0
+25.0
+29.0
+28.0
+26.0
+23.0
+26.0
+34.0
+29.0
+30.0
+28.0
+29.0
+27.0
+26.0
+30.0
+33.0
+25.0
+25.0
+24.0
+25.0
+29.0
+34.0
+24.0
+34.0
+27.0
+27.0
+28.0
+30.0
+27.0
+26.0
+21.0
+21.0
+19.0
+29.0
+28.0
+32.0
+29.0
+28.0
+27.0
+24.0
+28.0
+30.0
+26.0
+29.0
+24.0
+27.0
+32.0
+29.0
+31.0
+31.0
+24.0
+30.0
+26.0
+24.0
+27.0
+23.0
+25.0
+28.0
+28.0
+31.0
+24.0
+21.0
+28.0
+23.0
+28.0
+25.0
+26.0
+24.0
+28.0
+23.0
+24.0
+24.0
+31.0
+30.0
+26.0
+28.0
+28.0
+25.0
+28.0
+25.0
+28.0
+29.0
+27.0
+27.0
+33.0
+25.0
+29.0
+25.0
+25.0
+25.0
+32.0
+26.0
+28.0
+32.0
+29.0
+25.0
+30.0
+28.0
+26.0
+27.0
+20.0
+26.0
+30.0
+31.0
+29.0
+23.0
+27.0
+29.0
+28.0
+27.0
+23.0
+26.0
+29.0
+31.0
+27.0
+26.0
+33.0
+25.0
+26.0
+26.0
+30.0
+30.0
+26.0
+28.0
+26.0
+29.0
+28.0
+31.0
+29.0
+31.0
+23.0
+23.0
+28.0
+27.0
+26.0
+29.0
+26.0
+25.0
+30.0
+27.0
+25.0
+27.0
+24.0
+28.0
+29.0
+21.0
+24.0
+27.0
+26.0
+28.0
+30.0
+22.0
+25.0
+22.0
+25.0
+31.0
+29.0
+22.0
+20.0
+30.0
+28.0
+30.0
+29.0
+32.0
+27.0
+31.0
+26.0
+30.0
+23.0
+25.0
+28.0
+26.0
+29.0
+34.0
+26.0
+27.0
+33.0
+27.0
+25.0
+30.0
+30.0
+26.0
+25.0
+26.0
+30.0
+21.0
+25.0
+27.0
+29.0
+24.0
+22.0
+30.0
+23.0
+33.0
+28.0
+29.0
+22.0
+29.0
+28.0
+27.0
+26.0
+25.0
+29.0
+27.0
+32.0
+30.0
+24.0
+24.0
+26.0
+30.0
+26.0
+24.0
+28.0
+30.0
+20.0
+20.0
+28.0
+29.0
+26.0
+24.0
+25.0
+29.0
+26.0
+30.0
+25.0
+29.0
+27.0
+28.0
+28.0
+25.0
+36.0
+25.0
+30.0
+26.0
+20.0
+26.0
+22.0
+23.0
+27.0
+25.0
+24.0
+28.0
+27.0
+25.0
+28.0
+26.0
+24.0
+29.0
+22.0
+29.0
+30.0
+25.0
+26.0
+31.0
+22.0
+29.0
+32.0
+27.0
+28.0
+30.0
+29.0
+27.0
+27.0
+26.0
+27.0
+32.0
+30.0
+29.0
+28.0
+29.0
+28.0
+31.0
+33.0
+22.0
+30.0
+27.0
+30.0
+25.0
+29.0
+30.0
+28.0
+28.0
+25.0
+26.0
+26.0
+24.0
+25.0
+26.0
+23.0
+32.0
+31.0
+22.0
+28.0
+22.0
+29.0
+27.0
+29.0
+27.0
+24.0
+27.0
+32.0
+26.0
+24.0
+24.0
+26.0
+33.0
+23.0
+23.0
+27.0
+30.0
+31.0
+27.0
+24.0
+26.0
+31.0
+27.0
+26.0
+30.0
+27.0
+28.0
+31.0
+35.0
+30.0
+30.0
+28.0
+26.0
+22.0
+32.0
+28.0
+28.0
+26.0
+22.0
+27.0
+24.0
+27.0
+30.0
+32.0
+30.0
+28.0
+27.0
+28.0
+28.0
+20.0
+32.0
+26.0
+31.0
+27.0
+25.0
+26.0
+24.0
+29.0
+31.0
+29.0
+27.0
+26.0
+23.0
+27.0
+25.0
+28.0
+28.0
+24.0
+29.0
+30.0
+26.0
+29.0
+22.0
+28.0
+28.0
+28.0
+25.0
+28.0
+25.0
+29.0
+31.0
+33.0
+28.0
+26.0
+28.0
+31.0
+29.0
+25.0
+27.0
+28.0
+30.0
+26.0
+21.0
+29.0
+22.0
+31.0
+30.0
+33.0
+29.0
+21.0
+26.0
+30.0
+32.0
+25.0
+26.0
+31.0
+24.0
+24.0
+27.0
+29.0
+28.0
+28.0
+23.0
+30.0
+29.0
+25.0
+24.0
+30.0
+22.0
+29.0
+23.0
+26.0
+30.0
+27.0
+24.0
+33.0
+31.0
+25.0
+24.0
+30.0
+23.0
+27.0
+26.0
+26.0
+25.0
+30.0
+28.0
+33.0
+23.0
+28.0
+25.0
+33.0
+28.0
+31.0
+26.0
+22.0
+30.0
+28.0
+28.0
+24.0
+29.0
+28.0
+28.0
+30.0
+21.0
+32.0
+30.0
+27.0
+24.0
+26.0
+25.0
+25.0
+31.0
+29.0
+27.0
+23.0
+29.0
+29.0
+28.0
+26.0
+24.0
+29.0
+32.0
+25.0
+25.0
+26.0
+29.0
+27.0
+27.0
+28.0
+26.0
+29.0
+25.0
+29.0
+23.0
+27.0
+31.0
+27.0
+28.0
+28.0
+25.0
+21.0
+23.0
+29.0
+24.0
+29.0
+30.0
+28.0
+30.0
+29.0
+25.0
+25.0
+28.0
+26.0
+27.0
+27.0
+28.0
+25.0
+32.0
+26.0
+29.0
+28.0
+24.0
+28.0
+27.0
+24.0
+31.0
+27.0
+29.0
+26.0
+33.0
+26.0
+30.0
+32.0
+28.0
+25.0
+25.0
+27.0
+28.0
+30.0
+25.0
+33.0
+21.0
+31.0
+30.0
+26.0
+28.0
+29.0
+27.0
+24.0
+27.0
+27.0
+27.0
+26.0
+27.0
+32.0
+25.0
+30.0
+22.0
+25.0
+34.0
+26.0
+27.0
+33.0
+26.0
+27.0
+28.0
+24.0
+26.0
+25.0
+26.0
+25.0
+24.0
+29.0
+27.0
+25.0
+26.0
+22.0
+24.0
+22.0
+26.0
+32.0
+29.0
+27.0
+31.0
+26.0
+27.0
+34.0
+32.0
+26.0
+30.0
+24.0
+27.0
+24.0
+30.0
+31.0
+28.0
+27.0
+27.0
+29.0
+30.0
+28.0
+29.0
+24.0
+24.0
+29.0
+30.0
+31.0
+28.0
+27.0
+29.0
+28.0
+30.0
+26.0
+20.0
+28.0
+24.0
+28.0
+26.0
+20.0
+28.0
+32.0
+24.0
+27.0
+28.0
+24.0
+26.0
+29.0
+26.0
+28.0
+24.0
+29.0
+29.0
+30.0
+26.0
+24.0
+27.0
+24.0
+29.0
+26.0
+22.0
+26.0
+31.0
+27.0
+24.0
+30.0
+27.0
+24.0
+27.0
+28.0
+31.0
+28.0
+27.0
+22.0
+27.0
+32.0
+26.0
+30.0
+25.0
+32.0
+25.0
+26.0
+28.0
+26.0
+27.0
+26.0
+29.0
+29.0
+24.0
+22.0
+22.0
+26.0
+27.0
+32.0
+29.0
+27.0
+27.0
+25.0
+30.0
+30.0
+27.0
+28.0
+27.0
+26.0
+23.0
+22.0
+31.0
+28.0
+25.0
+28.0
+27.0
+29.0
+28.0
+24.0
+31.0
+28.0
+28.0
+24.0
+29.0
+26.0
+30.0
+27.0
+20.0
+25.0
+22.0
+28.0
+25.0
+29.0
+21.0
+27.0
+27.0
+25.0
+27.0
+25.0
+26.0
+25.0
+27.0
+26.0
+22.0
+21.0
+19.0
+28.0
diff --git a/your-code/ages_population3.csv b/your-code/ages_population3.csv
new file mode 100644
index 0000000..6339a1d
--- /dev/null
+++ b/your-code/ages_population3.csv
@@ -0,0 +1,1001 @@
+observation
+21.0
+21.0
+24.0
+31.0
+54.0
+52.0
+37.0
+69.0
+36.0
+30.0
+64.0
+30.0
+76.0
+73.0
+42.0
+52.0
+67.0
+28.0
+24.0
+50.0
+70.0
+57.0
+20.0
+27.0
+51.0
+67.0
+56.0
+30.0
+52.0
+44.0
+49.0
+31.0
+64.0
+65.0
+32.0
+37.0
+37.0
+18.0
+65.0
+42.0
+45.0
+40.0
+66.0
+72.0
+43.0
+64.0
+43.0
+54.0
+48.0
+30.0
+35.0
+41.0
+17.0
+44.0
+32.0
+37.0
+45.0
+41.0
+23.0
+45.0
+69.0
+36.0
+47.0
+31.0
+33.0
+51.0
+27.0
+20.0
+15.0
+26.0
+28.0
+67.0
+36.0
+54.0
+17.0
+51.0
+75.0
+41.0
+29.0
+55.0
+37.0
+63.0
+42.0
+49.0
+32.0
+39.0
+46.0
+50.0
+32.0
+43.0
+45.0
+26.0
+22.0
+15.0
+26.0
+24.0
+71.0
+55.0
+68.0
+35.0
+40.0
+55.0
+21.0
+60.0
+41.0
+22.0
+34.0
+30.0
+56.0
+61.0
+51.0
+33.0
+47.0
+48.0
+26.0
+30.0
+56.0
+42.0
+66.0
+50.0
+57.0
+12.0
+20.0
+69.0
+47.0
+38.0
+68.0
+41.0
+33.0
+29.0
+63.0
+38.0
+55.0
+39.0
+44.0
+74.0
+40.0
+35.0
+49.0
+37.0
+34.0
+25.0
+15.0
+70.0
+20.0
+35.0
+23.0
+54.0
+29.0
+27.0
+37.0
+27.0
+32.0
+34.0
+71.0
+38.0
+32.0
+65.0
+50.0
+23.0
+72.0
+48.0
+32.0
+52.0
+43.0
+34.0
+26.0
+42.0
+53.0
+54.0
+48.0
+23.0
+24.0
+37.0
+64.0
+70.0
+21.0
+50.0
+47.0
+67.0
+44.0
+63.0
+44.0
+12.0
+62.0
+48.0
+62.0
+70.0
+39.0
+25.0
+41.0
+59.0
+32.0
+43.0
+35.0
+63.0
+70.0
+36.0
+61.0
+38.0
+44.0
+48.0
+30.0
+31.0
+29.0
+31.0
+35.0
+50.0
+19.0
+45.0
+68.0
+18.0
+47.0
+32.0
+36.0
+36.0
+27.0
+37.0
+77.0
+64.0
+35.0
+48.0
+29.0
+37.0
+48.0
+39.0
+46.0
+72.0
+32.0
+49.0
+39.0
+20.0
+33.0
+41.0
+39.0
+52.0
+41.0
+35.0
+69.0
+25.0
+50.0
+14.0
+48.0
+30.0
+39.0
+69.0
+53.0
+43.0
+37.0
+37.0
+39.0
+32.0
+24.0
+19.0
+50.0
+17.0
+53.0
+36.0
+16.0
+29.0
+64.0
+36.0
+32.0
+1.0
+43.0
+26.0
+32.0
+40.0
+72.0
+35.0
+33.0
+24.0
+32.0
+38.0
+67.0
+29.0
+52.0
+55.0
+67.0
+22.0
+18.0
+41.0
+34.0
+65.0
+17.0
+50.0
+16.0
+45.0
+65.0
+46.0
+45.0
+31.0
+52.0
+71.0
+53.0
+52.0
+11.0
+23.0
+71.0
+52.0
+31.0
+27.0
+34.0
+51.0
+74.0
+19.0
+16.0
+25.0
+36.0
+13.0
+59.0
+50.0
+40.0
+45.0
+50.0
+43.0
+32.0
+36.0
+40.0
+42.0
+27.0
+33.0
+29.0
+75.0
+67.0
+42.0
+19.0
+47.0
+16.0
+46.0
+45.0
+35.0
+28.0
+36.0
+37.0
+59.0
+40.0
+45.0
+30.0
+35.0
+40.0
+44.0
+55.0
+28.0
+9.0
+28.0
+65.0
+31.0
+37.0
+28.0
+28.0
+40.0
+10.0
+48.0
+39.0
+69.0
+66.0
+54.0
+39.0
+65.0
+56.0
+19.0
+38.0
+21.0
+21.0
+24.0
+41.0
+50.0
+21.0
+69.0
+67.0
+27.0
+43.0
+57.0
+27.0
+39.0
+30.0
+51.0
+35.0
+20.0
+39.0
+35.0
+59.0
+60.0
+47.0
+24.0
+32.0
+46.0
+72.0
+69.0
+49.0
+33.0
+50.0
+35.0
+45.0
+52.0
+28.0
+70.0
+27.0
+67.0
+31.0
+28.0
+19.0
+71.0
+29.0
+41.0
+67.0
+32.0
+63.0
+19.0
+66.0
+71.0
+67.0
+18.0
+32.0
+2.0
+53.0
+71.0
+43.0
+70.0
+25.0
+66.0
+22.0
+33.0
+48.0
+38.0
+72.0
+24.0
+32.0
+39.0
+32.0
+22.0
+19.0
+59.0
+44.0
+36.0
+31.0
+69.0
+70.0
+46.0
+12.0
+31.0
+52.0
+25.0
+48.0
+17.0
+24.0
+40.0
+24.0
+41.0
+15.0
+46.0
+41.0
+44.0
+33.0
+64.0
+48.0
+68.0
+36.0
+41.0
+48.0
+32.0
+69.0
+66.0
+38.0
+25.0
+24.0
+43.0
+8.0
+60.0
+17.0
+23.0
+32.0
+68.0
+59.0
+39.0
+51.0
+64.0
+51.0
+25.0
+36.0
+56.0
+67.0
+56.0
+42.0
+48.0
+45.0
+14.0
+17.0
+15.0
+24.0
+60.0
+73.0
+65.0
+40.0
+35.0
+70.0
+34.0
+48.0
+45.0
+32.0
+38.0
+48.0
+49.0
+44.0
+43.0
+73.0
+69.0
+34.0
+24.0
+52.0
+46.0
+30.0
+38.0
+36.0
+34.0
+35.0
+68.0
+20.0
+45.0
+34.0
+70.0
+44.0
+56.0
+65.0
+69.0
+58.0
+72.0
+34.0
+21.0
+48.0
+20.0
+49.0
+34.0
+23.0
+67.0
+61.0
+32.0
+64.0
+49.0
+39.0
+63.0
+74.0
+67.0
+65.0
+50.0
+69.0
+34.0
+68.0
+38.0
+31.0
+37.0
+26.0
+23.0
+67.0
+29.0
+30.0
+26.0
+38.0
+68.0
+45.0
+54.0
+34.0
+32.0
+66.0
+43.0
+66.0
+27.0
+46.0
+28.0
+32.0
+32.0
+22.0
+40.0
+43.0
+63.0
+37.0
+17.0
+23.0
+27.0
+40.0
+36.0
+30.0
+29.0
+30.0
+30.0
+35.0
+42.0
+31.0
+37.0
+20.0
+30.0
+57.0
+17.0
+18.0
+68.0
+41.0
+39.0
+13.0
+69.0
+39.0
+28.0
+53.0
+42.0
+36.0
+47.0
+63.0
+42.0
+37.0
+18.0
+69.0
+35.0
+46.0
+35.0
+18.0
+49.0
+47.0
+38.0
+27.0
+34.0
+57.0
+29.0
+38.0
+39.0
+57.0
+33.0
+36.0
+40.0
+49.0
+48.0
+54.0
+39.0
+49.0
+37.0
+49.0
+66.0
+70.0
+29.0
+61.0
+4.0
+37.0
+36.0
+18.0
+38.0
+66.0
+44.0
+35.0
+31.0
+31.0
+38.0
+67.0
+52.0
+35.0
+51.0
+54.0
+43.0
+40.0
+51.0
+21.0
+40.0
+61.0
+31.0
+22.0
+41.0
+73.0
+68.0
+63.0
+30.0
+46.0
+40.0
+12.0
+15.0
+26.0
+39.0
+70.0
+36.0
+67.0
+45.0
+38.0
+68.0
+43.0
+28.0
+58.0
+52.0
+53.0
+44.0
+73.0
+49.0
+26.0
+34.0
+37.0
+38.0
+47.0
+52.0
+35.0
+29.0
+51.0
+36.0
+4.0
+28.0
+30.0
+38.0
+37.0
+22.0
+74.0
+35.0
+39.0
+73.0
+65.0
+30.0
+66.0
+70.0
+58.0
+35.0
+21.0
+61.0
+66.0
+41.0
+42.0
+64.0
+31.0
+67.0
+38.0
+46.0
+68.0
+65.0
+37.0
+36.0
+40.0
+31.0
+48.0
+46.0
+50.0
+47.0
+53.0
+28.0
+41.0
+68.0
+56.0
+26.0
+41.0
+63.0
+69.0
+67.0
+24.0
+25.0
+43.0
+67.0
+32.0
+32.0
+30.0
+65.0
+44.0
+52.0
+11.0
+25.0
+60.0
+49.0
+25.0
+33.0
+44.0
+10.0
+23.0
+22.0
+29.0
+22.0
+19.0
+58.0
+71.0
+33.0
+66.0
+31.0
+64.0
+28.0
+53.0
+33.0
+45.0
+21.0
+68.0
+37.0
+19.0
+48.0
+51.0
+35.0
+40.0
+53.0
+71.0
+51.0
+27.0
+60.0
+21.0
+45.0
+74.0
+46.0
+29.0
+44.0
+66.0
+55.0
+31.0
+48.0
+31.0
+39.0
+24.0
+55.0
+70.0
+24.0
+41.0
+68.0
+25.0
+39.0
+35.0
+39.0
+45.0
+29.0
+48.0
+65.0
+49.0
+34.0
+68.0
+43.0
+39.0
+32.0
+26.0
+47.0
+25.0
+37.0
+53.0
+37.0
+66.0
+57.0
+32.0
+44.0
+34.0
+47.0
+46.0
+24.0
+70.0
+56.0
+27.0
+28.0
+19.0
+38.0
+45.0
+35.0
+71.0
+53.0
+41.0
+65.0
+25.0
+46.0
+5.0
+14.0
+15.0
+35.0
+50.0
+30.0
+58.0
+46.0
+27.0
+45.0
+37.0
+50.0
+45.0
+40.0
+32.0
+20.0
+38.0
+46.0
+22.0
+50.0
+42.0
+66.0
+34.0
+32.0
+29.0
+30.0
+62.0
+66.0
+16.0
+63.0
+51.0
+25.0
+53.0
+41.0
+35.0
+37.0
+33.0
+25.0
+32.0
+55.0
+29.0
+36.0
+45.0
+14.0
+49.0
+32.0
+38.0
+37.0
+21.0
+29.0
+52.0
+46.0
+42.0
+57.0
+46.0
+58.0
+18.0
+48.0
+71.0
+41.0
+50.0
+43.0
+26.0
+30.0
+27.0
+47.0
+43.0
+38.0
+64.0
+36.0
+39.0
+74.0
+39.0
+56.0
+2.0
+39.0
+70.0
+34.0
+47.0
+59.0
+66.0
+36.0
+33.0
+41.0
+55.0
+66.0
+29.0
+77.0
+63.0
+68.0
+50.0
+50.0
+17.0
+24.0
+16.0
+59.0
+70.0
+16.0
+33.0
+41.0
+32.0
+40.0
+26.0
+7.0
+27.0
+29.0
+29.0
+19.0
+70.0
+23.0
+27.0
+36.0
+55.0
+21.0
+15.0
+14.0
+46.0
+37.0
+67.0
+29.0
+28.0
+33.0
+46.0
+55.0
+35.0
+49.0
+72.0
+20.0
+45.0
+69.0
+16.0
+53.0
+39.0
+38.0
+49.0
+25.0
+43.0
+28.0
+32.0
+67.0
+29.0
+70.0
+26.0
+48.0
+37.0
+16.0
+55.0
+30.0
+35.0
+43.0
diff --git a/your-code/main.ipynb b/your-code/main.ipynb
index a0a5b66..8e8dcf6 100644
--- a/your-code/main.ipynb
+++ b/your-code/main.ipynb
@@ -11,11 +11,15 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
- "# Libraries"
+ "from sklearn.linear_model import LinearRegression\n",
+ "import pandas as pd\n",
+ "import numpy as np\n",
+ "import matplotlib.pyplot as plt\n",
+ "from pylab import rcParams # library for plotting histograms"
]
},
{
@@ -29,11 +33,122 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
- "# your code here"
+ "import random\n",
+ "\n",
+ "def rolling_dice(times):\n",
+ " \n",
+ " dice = [1,2,3,4,5,6]\n",
+ " result = []\n",
+ " counter = 0\n",
+ " while counter < times:\n",
+ " result.append(random.choices(dice))\n",
+ " counter += 1\n",
+ " result_df = pd.DataFrame(result)\n",
+ " result_df[0] = result_df.rename(columns={0: \"value\"}, inplace=True)\n",
+ " result_df = result_df[[\"value\"]]\n",
+ " return result_df"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 5 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 6 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 7 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 8 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 9 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 6 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 7 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 8 \n",
+ " 4 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ " 5 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 9 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# plotting the value\n",
+ "result.plot(kind='bar', alpha=0.7)\n",
+ "plt.xlabel('roll(index)')\n",
+ "plt.ylabel('value')\n",
+ "plt.title('')\n",
+ "plt.show()"
]
},
{
@@ -61,22 +291,62 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 6,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "4 4\n",
+ "2 3\n",
+ "6 2\n",
+ "5 1\n",
+ "Name: value, dtype: int64"
+ ]
+ },
+ "execution_count": 6,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "# calculating FD\n",
+ "frequency_distribution = result['value'].value_counts()\n",
+ "frequency_distribution"
]
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 7,
+ "metadata": {
+ "scrolled": true
+ },
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAHACAYAAABeV0mSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA1oUlEQVR4nO3de1xVdb7/8fdOZCMKKigIeS3LC+YNnMLESyQmHnS8nJOTeUltxtHSJE7ngc2ccmpCJ3XMLNBJxUveipx0MEdLISs9ieKtTK1UyMBrgVJubuv3hw/3zy14YQvszfL1fDz2H+u7vmutz9qr4t13fddeFsMwDAEAAJjEXa4uAAAAoDIRbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQbgAAgKkQboAqlJycLIvFUu4nLi7O1eWZyrXftZeXl5o0aaI+ffooISFBp0+fLrPNyy+/LIvFUqHj/PLLL3r55ZeVlpZWoe3KO1bLli31H//xHxXaz82sXLlSc+fOLXedxWLRyy+/XKnHA9yRh6sLAO4ES5YsUdu2bR3agoODXVSNuV35rouKinT69Gl99tlnmjlzpmbNmqU1a9bo0UcftfcdP368HnvssQrt/5dfftH06dMlSb17977l7Zw5ljNWrlypgwcP6rnnniuzbseOHWratGmV1wC4GuEGqAYdOnRQWFjYLfUtKiqSxWKRhwf/ejrj2u966NChmjp1qnr06KEhQ4bo6NGjCgwMlCQ1bdq0yv/Y//LLL/L29q6WY93MQw895NLjA9WF21KAC6WlpclisWj58uV6/vnndffdd8tqterbb7+VJH388ceKjIyUr6+vvL299fDDD+uTTz4ps5/U1FR17txZVqtVrVq10qxZs8rcBjl+/LgsFouSk5PLbF/e7YqjR4/qiSeeUEBAgKxWq9q1a6e33nqr3PpXrVqlF198UcHBwfL19dWjjz6qw4cPlznOpk2bFBkZqfr168vb21vt2rVTQkKCJGn58uWyWCzasWNHme3+8pe/qHbt2vrxxx9v+p2Wp3nz5po9e7YuXLigBQsW2NvLu1W0detW9e7dW/7+/qpTp46aN2+uoUOH6pdfftHx48fVuHFjSdL06dPtt8DGjBnjsL89e/Zo2LBhatiwoe69997rHuuKdevWqWPHjvLy8tI999yjefPmOay/csvt+PHjDu1Xvv8rt8h69+6t1NRUnThxwuEW3RXlXeeDBw9q0KBBatiwoby8vNS5c2ctXbq03OPc6nUGXI1wA1SDkpISFRcXO3yuFh8fr6ysLCUlJWnDhg0KCAjQihUrFBUVJV9fXy1dulRr166Vn5+f+vXr5xBwPvnkEw0aNEg+Pj5avXq1Xn/9da1du1ZLlixxut6vv/5a3bp108GDBzV79mz961//0oABAzR58mT7LZmrTZs2TSdOnNA777yjhQsX6ujRo4qJiVFJSYm9z6JFixQdHa3S0lL7eU6ePFk//PCDJOnxxx9XkyZNygSo4uJiLViwQIMHD76tW3nR0dGqVauWPv300+v2OX78uAYMGCBPT08tXrxYmzZt0owZM1S3bl0VFhYqKChImzZtkiSNGzdOO3bs0I4dO/TnP//ZYT9DhgxR69at9d577ykpKemGde3du1fPPfecpk6dqnXr1ql79+6aMmWKZs2aVeFzfPvtt/Xwww+rSZMm9trKC4tXHD58WN27d9dXX32lefPm6YMPPlD79u01ZswY/e1vfyvT/1auM+AWDABVZsmSJYakcj9FRUXGtm3bDElGz549HbYrKCgw/Pz8jJiYGIf2kpISo1OnTsZvfvMbe9uDDz5oBAcHG7/++qu9LT8/3/Dz8zOu/lf82LFjhiRjyZIlZeqUZLz00kv25X79+hlNmzY18vLyHPo988wzhpeXl3H+/HnDMAx7/dHR0Q791q5da0gyduzYYRiGYVy4cMHw9fU1evToYZSWll73+3rppZcMT09P49SpU/a2NWvWGJKM9PT0625nGP//u961a9d1+wQGBhrt2rVzON7V39H7779vSDL27t173X2cOXOmzPd17f7+93//97rrrtaiRQvDYrGUOV7fvn0NX19fo6CgwOHcjh075tDvyve/bds2e9uAAQOMFi1alFv7tXUPHz7csFqtRlZWlkO//v37G97e3sbPP//scJybXWfAXTByA1SDZcuWadeuXQ6fq+fUDB061KH/F198ofPnz2v06NEOoz2lpaV67LHHtGvXLhUUFKigoEC7du3SkCFD5OXlZd/ex8dHMTExTtV66dIlffLJJxo8eLC8vb0djh8dHa1Lly5p586dDtsMHDjQYbljx46SpBMnTtjPJz8/XxMnTrzh00l//OMfJUn/+Mc/7G3z58/XAw88oJ49ezp1PlczDOOG6zt37ixPT0/9/ve/19KlS/X99987dZxrr+eNhISEqFOnTg5tTzzxhPLz87Vnzx6njn+rtm7dqsjISDVr1syhfcyYMfrll1/KjPrc7DoD7oJwA1SDdu3aKSwszOFztaCgIIflU6dOSZKGDRum2rVrO3xmzpwpwzB0/vx5/fTTTyotLVWTJk3KHLO8tltx7tw5FRcX68033yxz7OjoaEnS2bNnHbbx9/d3WLZarZKkX3/9VZJ05swZSbrphNrAwEA9/vjjWrBggUpKSrR//35t375dzzzzjFPncrWCggKdO3fuhre27r33Xn388ccKCAjQpEmTdO+99+ree+/VG2+8UaFjXXs9b+RG1+7cuXMVOm5FnTt3rtxar3xH1x7/ZtcZcBc8jgG4gWtHMxo1aiRJevPNN6/7hEtgYKD9yarc3Nwy669tuzKyY7PZHNqv/QPWsGFD1apVSyNHjtSkSZPKPXarVq1ucDZlXZmEe2V+zY1MmTJFy5cv14cffqhNmzapQYMGGjFiRIWOV57U1FSVlJTc9PHtiIgIRUREqKSkRBkZGXrzzTf13HPPKTAwUMOHD7+lY1Xkt3NudO2uhInrXbtrQ2ZF+fv7Kycnp0z7lYnbV/45BGoawg3ghh5++GE1aNBAX3/99Q1HLTw9PfWb3/xGH3zwgV5//XX7H8ELFy5ow4YNDn0DAwPl5eWl/fv3O7R/+OGHDsve3t7q06ePMjMz1bFjR3l6et72+XTv3l3169dXUlKShg8ffsM//qGhoerevbtmzpypgwcP6ve//73q1q17W8fPyspSXFyc6tevrz/84Q+3tE2tWrX04IMPqm3btnr33Xe1Z88eDR8+vNJHK7766ivt27fP4dbUypUr5ePjo65du0q6/GN/krR//361adPG3m/9+vVl9me1Wm+5tsjISK1bt04//vijw4jWsmXL5O3tzaPjqLEIN4Abqlevnt58802NHj1a58+f17BhwxQQEKAzZ85o3759OnPmjBITEyVJr7zyih577DH17dtXzz//vEpKSjRz5kzVrVtX58+ft+/TYrHoySef1OLFi3XvvfeqU6dO+vLLL7Vy5coyx3/jjTfUo0cPRURE6I9//KNatmypCxcu6Ntvv9WGDRu0devWCp/P7NmzNX78eD366KN6+umnFRgYqG+//Vb79u3T/PnzHfpPmTJFjz/+uCwWiyZOnFihYx08eNA+R+j06dPavn27lixZolq1amndunX2UaTyJCUlaevWrRowYICaN2+uS5cuafHixZJk//E/Hx8ftWjRQh9++KEiIyPl5+enRo0a2QNIRQUHB2vgwIF6+eWXFRQUpBUrVmjLli2aOXOmvL29JUndunVTmzZtFBcXp+LiYjVs2FDr1q3TZ599VmZ/DzzwgD744AMlJiYqNDRUd91113V/Y+mll17Sv/71L/Xp00f/+7//Kz8/P7377rtKTU3V3/72N9WvX9+pcwJcztUzmgEzu9kTPFeeQnnvvffKXZ+enm4MGDDA8PPzM2rXrm3cfffdxoABA8r0X79+vdGxY0fD09PTaN68uTFjxoxyn87Jy8szxo8fbwQGBhp169Y1YmJijOPHj5f79M+xY8eMsWPHGnfffbdRu3Zto3Hjxkb37t2NV1999ab1X+/JrI0bNxq9evUy6tata3h7exvt27c3Zs6cWea8bTabYbVajccee6zc76U81z6Z5unpaQQEBBi9evUyXnvtNeP06dNltrn2O9qxY4cxePBgo0WLFobVajX8/f2NXr16GevXr3fY7uOPPza6dOliWK1WQ5IxevRoh/2dOXPmpscyjMtPSw0YMMB4//33jZCQEMPT09No2bKlMWfOnDLbHzlyxIiKijJ8fX2Nxo0bG88++6yRmppa5mmp8+fPG8OGDTMaNGhgWCwWh2OWd50PHDhgxMTEGPXr1zc8PT2NTp06lbluFb3OgKtZDOMmjw8AqJFefvllTZ8+/aZPCLmjDRs2aODAgUpNTbVPYgaAW8VtKQBu4+uvv9aJEyf0/PPPq3Pnzurfv7+rSwJQA/EoOAC3MXHiRA0cOFANGzbUqlWrKvzGbgCQJG5LAQAAU2HkBgAAmArhBgAAmArhBgAAmMod97RUaWmpfvzxR/n4+DBZEQCAGsIwDF24cEHBwcG6664bj83cceHmxx9/LPMGXAAAUDNkZ2ff9CW8d1y48fHxkXT5y/H19XVxNQAA4Fbk5+erWbNm9r/jN3LHhZsrt6J8fX0JNwAA1DC3MqWECcUAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBUCDcAAMBU3CbcJCQkyGKx6Lnnnrthv/T0dIWGhsrLy0v33HOPkpKSqqdAAABQI7hFuNm1a5cWLlyojh073rDfsWPHFB0drYiICGVmZmratGmaPHmyUlJSqqlSAADg7lwebi5evKgRI0boH//4hxo2bHjDvklJSWrevLnmzp2rdu3aafz48Ro7dqxmzZpVTdUCAAB35/JwM2nSJA0YMECPPvroTfvu2LFDUVFRDm39+vVTRkaGioqKqqpEAABQg3i48uCrV6/Wnj17tGvXrlvqn5ubq8DAQIe2wMBAFRcX6+zZswoKCiqzjc1mk81msy/n5+ffXtEAAMCtuSzcZGdna8qUKdq8ebO8vLxueTuLxeKwbBhGue1XJCQkaPr06c4XWsnGJd9akDObRWO6uboEAMAdwmW3pXbv3q3Tp08rNDRUHh4e8vDwUHp6uubNmycPDw+VlJSU2aZJkybKzc11aDt9+rQ8PDzk7+9f7nHi4+OVl5dn/2RnZ1fJ+QAAAPfgspGbyMhIHThwwKHtqaeeUtu2bfU///M/qlWrVpltwsPDtWHDBoe2zZs3KywsTLVr1y73OFarVVartfIKBwAAbs1l4cbHx0cdOnRwaKtbt678/f3t7fHx8Tp58qSWLVsmSZowYYLmz5+v2NhYPf3009qxY4cWLVqkVatWVXv9AADAPbn8aakbycnJUVZWln25VatW2rhxo9LS0tS5c2e98sormjdvnoYOHerCKgEAgDtx6dNS10pLS3NYTk5OLtOnV69e2rNnT/UUBAAAahy3HrkBAACoKMINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFcINAAAwFZeGm8TERHXs2FG+vr7y9fVVeHi4Pvroo+v2T0tLk8ViKfP55ptvqrFqAADgzjxcefCmTZtqxowZat26tSRp6dKlGjRokDIzMxUSEnLd7Q4fPixfX1/7cuPGjau8VgAAUDO4NNzExMQ4LP/1r39VYmKidu7cecNwExAQoAYNGlRxdQAAoCZymzk3JSUlWr16tQoKChQeHn7Dvl26dFFQUJAiIyO1bdu2aqoQAADUBC4duZGkAwcOKDw8XJcuXVK9evW0bt06tW/fvty+QUFBWrhwoUJDQ2Wz2bR8+XJFRkYqLS1NPXv2LHcbm80mm81mX87Pz6+S8wAAAO7B5eGmTZs22rt3r37++WelpKRo9OjRSk9PLzfgtGnTRm3atLEvh4eHKzs7W7NmzbpuuElISND06dOrrH4AAOBeXH5bytPTU61bt1ZYWJgSEhLUqVMnvfHGG7e8/UMPPaSjR49ed318fLzy8vLsn+zs7MooGwAAuCmXj9xcyzAMh9tIN5OZmamgoKDrrrdarbJarZVRGgAAqAFcGm6mTZum/v37q1mzZrpw4YJWr16ttLQ0bdq0SdLlUZeTJ09q2bJlkqS5c+eqZcuWCgkJUWFhoVasWKGUlBSlpKS48jQAAIAbcWm4OXXqlEaOHKmcnBzVr19fHTt21KZNm9S3b19JUk5OjrKysuz9CwsLFRcXp5MnT6pOnToKCQlRamqqoqOjXXUKAADAzVgMwzBcXUR1ys/PV/369ZWXl+fwQ4DVZVzyrmo/pjtYNKabq0sAANRgFfn77fIJxQAAAJWJcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEzFpeEmMTFRHTt2lK+vr3x9fRUeHq6PPvrohtukp6crNDRUXl5euueee5SUlFRN1QIAgJrApeGmadOmmjFjhjIyMpSRkaFHHnlEgwYN0ldffVVu/2PHjik6OloRERHKzMzUtGnTNHnyZKWkpFRz5QAAwF15uPLgMTExDst//etflZiYqJ07dyokJKRM/6SkJDVv3lxz586VJLVr104ZGRmaNWuWhg4dWh0lAwAAN+c2c25KSkq0evVqFRQUKDw8vNw+O3bsUFRUlENbv379lJGRoaKiouooEwAAuDmXjtxI0oEDBxQeHq5Lly6pXr16Wrdundq3b19u39zcXAUGBjq0BQYGqri4WGfPnlVQUFCZbWw2m2w2m305Pz+/ck8AAAC4FZeHmzZt2mjv3r36+eeflZKSotGjRys9Pf26AcdisTgsG4ZRbvsVCQkJmj59euUWDdyiccm7XF2CSywa083VJQC4g7n8tpSnp6dat26tsLAwJSQkqFOnTnrjjTfK7dukSRPl5uY6tJ0+fVoeHh7y9/cvd5v4+Hjl5eXZP9nZ2ZV+DgAAwH24fOTmWoZhONxGulp4eLg2bNjg0LZ582aFhYWpdu3a5W5jtVpltVorvU4AAOCeXDpyM23aNG3fvl3Hjx/XgQMH9OKLLyotLU0jRoyQdHnUZdSoUfb+EyZM0IkTJxQbG6tDhw5p8eLFWrRokeLi4lx1CgAAwM24dOTm1KlTGjlypHJyclS/fn117NhRmzZtUt++fSVJOTk5ysrKsvdv1aqVNm7cqKlTp+qtt95ScHCw5s2bx2PgAADAzqXhZtGiRTdcn5ycXKatV69e2rNnTxVVBAAAajqXTygGAACoTIQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKoQbAABgKi4NNwkJCerWrZt8fHwUEBCg3/72tzp8+PANt0lLS5PFYinz+eabb6qpagAA4M5cGm7S09M1adIk7dy5U1u2bFFxcbGioqJUUFBw020PHz6snJwc++e+++6rhooBAIC783DlwTdt2uSwvGTJEgUEBGj37t3q2bPnDbcNCAhQgwYNqrA6AABQE7nVnJu8vDxJkp+f3037dunSRUFBQYqMjNS2bduqujQAAFBDuHTk5mqGYSg2NlY9evRQhw4drtsvKChICxcuVGhoqGw2m5YvX67IyEilpaWVO9pjs9lks9nsy/n5+VVSPwAAcA9uE26eeeYZ7d+/X5999tkN+7Vp00Zt2rSxL4eHhys7O1uzZs0qN9wkJCRo+vTplV4vAABwT25xW+rZZ5/V+vXrtW3bNjVt2rTC2z/00EM6evRouevi4+OVl5dn/2RnZ99uuQAAwI05NXJz7NgxtWrV6rYPbhiGnn32Wa1bt05paWlO7zMzM1NBQUHlrrNarbJarbdTJgAAqEGcCjetW7dWz549NW7cOA0bNkxeXl5OHXzSpElauXKlPvzwQ/n4+Cg3N1eSVL9+fdWpU0fS5ZGXkydPatmyZZKkuXPnqmXLlgoJCVFhYaFWrFihlJQUpaSkOFUDAAAwF6duS+3bt09dunTR888/ryZNmugPf/iDvvzyywrvJzExUXl5eerdu7eCgoLsnzVr1tj75OTkKCsry75cWFiouLg4dezYUREREfrss8+UmpqqIUOGOHMqAADAZCyGYRjOblxcXKwNGzYoOTlZH330ke677z6NGzdOI0eOVOPGjSuzzkqTn5+v+vXrKy8vT76+vtV+/HHJu6r9mO5g0Zhuri7BJbjeAFA5KvL3+7YmFHt4eGjw4MFau3atZs6cqe+++05xcXFq2rSpRo0apZycnNvZPQAAQIXdVrjJyMjQxIkTFRQUpDlz5iguLk7fffedtm7dqpMnT2rQoEGVVScAAMAtcWpC8Zw5c7RkyRIdPnxY0dHRWrZsmaKjo3XXXZezUqtWrbRgwQK1bdu2UosFAAC4GafCTWJiosaOHaunnnpKTZo0KbdP8+bNtWjRotsqDgAAoKKcCjfX+8G8q3l6emr06NHO7B4AAMBpTs25WbJkid57770y7e+9956WLl1620UBAAA4y6lwM2PGDDVq1KhMe0BAgF577bXbLgoAAMBZToWbEydOlPuqhBYtWjj84B4AAEB1cyrcBAQEaP/+/WXa9+3bJ39//9suCgAAwFlOhZvhw4dr8uTJ2rZtm0pKSlRSUqKtW7dqypQpGj58eGXXCAAAcMucelrq1Vdf1YkTJxQZGSkPj8u7KC0t1ahRo5hzAwAAXMqpcOPp6ak1a9bolVde0b59+1SnTh098MADatGiRWXXBwAAUCFOhZsr7r//ft1///2VVQsAAMBtcyrclJSUKDk5WZ988olOnz6t0tJSh/Vbt26tlOIAAAAqyqlwM2XKFCUnJ2vAgAHq0KGDLBZLZdcFAADgFKfCzerVq7V27VpFR0dXdj0AAAC3xalHwT09PdW6devKrgUAAOC2ORVunn/+eb3xxhsyDKOy6wEAALgtTt2W+uyzz7Rt2zZ99NFHCgkJUe3atR3Wf/DBB5VSHAAAQEU5FW4aNGigwYMHV3YtAAAAt82pcLNkyZLKrgMAAKBSODXnRpKKi4v18ccfa8GCBbpw4YIk6ccff9TFixcrrTgAAICKcmrk5sSJE3rssceUlZUlm82mvn37ysfHR3/729906dIlJSUlVXadAAAAt8SpkZspU6YoLCxMP/30k+rUqWNvHzx4sD755JNKKw4AAKCinH5a6vPPP5enp6dDe4sWLXTy5MlKKQwAAMAZTo3clJaWqqSkpEz7Dz/8IB8fn9suCgAAwFlOhZu+fftq7ty59mWLxaKLFy/qpZde4pUMAADApZy6LfX3v/9dffr0Ufv27XXp0iU98cQTOnr0qBo1aqRVq1ZVdo0AAAC3zKlwExwcrL1792rVqlXas2ePSktLNW7cOI0YMcJhgjEAAEB1cyrcSFKdOnU0duxYjR07tjLrAQAAuC1OhZtly5bdcP2oUaOcKgYAAOB2ORVupkyZ4rBcVFSkX375RZ6envL29ibcAAAAl3HqaamffvrJ4XPx4kUdPnxYPXr0YEIxAABwKaffLXWt++67TzNmzCgzqnMjCQkJ6tatm3x8fBQQEKDf/va3Onz48E23S09PV2hoqLy8vHTPPffwugcAAGBXaeFGkmrVqqUff/zxlvunp6dr0qRJ2rlzp7Zs2aLi4mJFRUWpoKDgutscO3ZM0dHRioiIUGZmpqZNm6bJkycrJSWlMk4BAADUcE7NuVm/fr3DsmEYysnJ0fz58/Xwww/f8n42bdrksLxkyRIFBARo9+7d6tmzZ7nbJCUlqXnz5vYfEWzXrp0yMjI0a9YsDR06tGInAgAATMepcPPb3/7WYdlisahx48Z65JFHNHv2bKeLycvLkyT5+fldt8+OHTsUFRXl0NavXz8tWrRIRUVFql27ttPHBwAANZ9T4aa0tLSy65BhGIqNjVWPHj3UoUOH6/bLzc1VYGCgQ1tgYKCKi4t19uxZBQUFOayz2Wyy2Wz25fz8/MotHAAAuBWnf8Svsj3zzDPav3+/Pvvss5v2tVgsDsuGYZTbLl2etDx9+vTKKRIAbmBc8i5Xl+ASi8Z0c3UJgAOnwk1sbOwt950zZ85N+zz77LNav369Pv30UzVt2vSGfZs0aaLc3FyHttOnT8vDw0P+/v5l+sfHxzvUm5+fr2bNmt1i9QAAoKZxKtxkZmZqz549Ki4uVps2bSRJR44cUa1atdS1a1d7v/JGUq5mGIaeffZZrVu3TmlpaWrVqtVNjx0eHq4NGzY4tG3evFlhYWHlzrexWq2yWq23cloAAMAEnAo3MTEx8vHx0dKlS9WwYUNJl3/Y76mnnlJERISef/75W9rPpEmTtHLlSn344Yfy8fGxj8jUr1/f/gLO+Ph4nTx50v7KhwkTJmj+/PmKjY3V008/rR07dmjRokX8eCAAAJDk5O/czJ49WwkJCfZgI0kNGzbUq6++WqGnpRITE5WXl6fevXsrKCjI/lmzZo29T05OjrKysuzLrVq10saNG5WWlqbOnTvrlVde0bx583gMHAAASHJy5CY/P1+nTp1SSEiIQ/vp06d14cKFW97PlYnAN5KcnFymrVevXtqzZ88tHwcAANw5nBq5GTx4sJ566im9//77+uGHH/TDDz/o/fff17hx4zRkyJDKrhEAAOCWOTVyk5SUpLi4OD355JMqKiq6vCMPD40bN06vv/56pRYIAABQEU6FG29vb7399tt6/fXX9d1338kwDLVu3Vp169at7PoAAAAq5LZenJmTk6OcnBzdf//9qlu37i3NoQEAAKhKToWbc+fOKTIyUvfff7+io6OVk5MjSRo/fvwtPwYOAABQFZwKN1OnTlXt2rWVlZUlb29ve/vjjz9e5k3fAAAA1cmpOTebN2/Wv//97zKvSrjvvvt04sSJSikMAADAGU6N3BQUFDiM2Fxx9uxZXnUAAABcyqlw07NnT/vrEKTL75AqLS3V66+/rj59+lRacQAAABXl1G2p119/Xb1791ZGRoYKCwv1wgsv6KuvvtL58+f1+eefV3aNAAAAt8ypkZv27dtr//79+s1vfqO+ffuqoKBAQ4YMUWZmpu69997KrhEAAOCWVXjkpqioSFFRUVqwYIGmT59eFTUBAAA4rcIjN7Vr19bBgwdlsViqoh4AAIDb4tRtqVGjRmnRokWVXQsAAMBtc2pCcWFhod555x1t2bJFYWFhZd4pNWfOnEopDgAAoKIqFG6+//57tWzZUgcPHlTXrl0lSUeOHHHow+0qAADgShUKN/fdd59ycnK0bds2SZdftzBv3jwFBgZWSXEAAAAVVaE5N9e+9fujjz5SQUFBpRYEAABwO5yaUHzFtWEHAADA1SoUbiwWS5k5NcyxAQAA7qRCc24Mw9CYMWPsL8e8dOmSJkyYUOZpqQ8++KDyKgQAAKiACoWb0aNHOyw/+eSTlVoMAADA7apQuFmyZElV1QEAAFApbmtCMQAAgLsh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFMh3AAAAFNxabj59NNPFRMTo+DgYFksFv3zn/+8Yf+0tDT7m8mv/nzzzTfVUzAAAHB7FXq3VGUrKChQp06d9NRTT2no0KG3vN3hw4fl6+trX27cuHFVlAcAAGogl4ab/v37q3///hXeLiAgQA0aNKj8ggAAQI1XI+fcdOnSRUFBQYqMjNS2bdtcXQ4AAHAjLh25qaigoCAtXLhQoaGhstlsWr58uSIjI5WWlqaePXuWu43NZpPNZrMv5+fnV1e5AADABWpUuGnTpo3atGljXw4PD1d2drZmzZp13XCTkJCg6dOnV1eJAADAxWrkbamrPfTQQzp69Oh118fHxysvL8/+yc7OrsbqAABAdatRIzflyczMVFBQ0HXXW61WWa3WaqwIAAC4kkvDzcWLF/Xtt9/al48dO6a9e/fKz89PzZs3V3x8vE6ePKlly5ZJkubOnauWLVsqJCREhYWFWrFihVJSUpSSkuKqUwAAAG7GpeEmIyNDffr0sS/HxsZKkkaPHq3k5GTl5OQoKyvLvr6wsFBxcXE6efKk6tSpo5CQEKWmpio6OrraawcAAO7JpeGmd+/eMgzjuuuTk5Mdll944QW98MILVVwVAACoyWr8hGIAAICrEW4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpuDTcfPrpp4qJiVFwcLAsFov++c9/3nSb9PR0hYaGysvLS/fcc4+SkpKqvlAAAFBjuDTcFBQUqFOnTpo/f/4t9T927Jiio6MVERGhzMxMTZs2TZMnT1ZKSkoVVwoAAGoKD1cevH///urfv/8t909KSlLz5s01d+5cSVK7du2UkZGhWbNmaejQoVVUJQAAqElq1JybHTt2KCoqyqGtX79+ysjIUFFRkYuqAgAA7sSlIzcVlZubq8DAQIe2wMBAFRcX6+zZswoKCiqzjc1mk81msy/n5+dXeZ0AAMB1alS4kSSLxeKwbBhGue1XJCQkaPr06VVeFwDgzjIueZerS3CJRWO6ubqEm6pRt6WaNGmi3Nxch7bTp0/Lw8ND/v7+5W4THx+vvLw8+yc7O7s6SgUAAC5So0ZuwsPDtWHDBoe2zZs3KywsTLVr1y53G6vVKqvVWh3lAQAAN+DSkZuLFy9q79692rt3r6TLj3rv3btXWVlZki6PuowaNcref8KECTpx4oRiY2N16NAhLV68WIsWLVJcXJwrygcAAG7IpSM3GRkZ6tOnj305NjZWkjR69GglJycrJyfHHnQkqVWrVtq4caOmTp2qt956S8HBwZo3bx6PgQMAADuXhpvevXvbJwSXJzk5uUxbr169tGfPniqsCgAA1GQ1akIxAADAzRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqRBuAACAqbg83Lz99ttq1aqVvLy8FBoaqu3bt1+3b1pamiwWS5nPN998U40VAwAAd+bScLNmzRo999xzevHFF5WZmamIiAj1799fWVlZN9zu8OHDysnJsX/uu+++aqoYAAC4O5eGmzlz5mjcuHEaP3682rVrp7lz56pZs2ZKTEy84XYBAQFq0qSJ/VOrVq1qqhgAALg7l4WbwsJC7d69W1FRUQ7tUVFR+uKLL264bZcuXRQUFKTIyEht27atKssEAAA1jIerDnz27FmVlJQoMDDQoT0wMFC5ubnlbhMUFKSFCxcqNDRUNptNy5cvV2RkpNLS0tSzZ89yt7HZbLLZbPbl/Pz8yjsJAADgdlwWbq6wWCwOy4ZhlGm7ok2bNmrTpo19OTw8XNnZ2Zo1a9Z1w01CQoKmT59eeQUDAAC35rLbUo0aNVKtWrXKjNKcPn26zGjOjTz00EM6evToddfHx8crLy/P/snOzna6ZgAA4P5cFm48PT0VGhqqLVu2OLRv2bJF3bt3v+X9ZGZmKigo6LrrrVarfH19HT4AAMC8XHpbKjY2ViNHjlRYWJjCw8O1cOFCZWVlacKECZIuj7qcPHlSy5YtkyTNnTtXLVu2VEhIiAoLC7VixQqlpKQoJSXFlacBAADciEvDzeOPP65z587pL3/5i3JyctShQwdt3LhRLVq0kCTl5OQ4/OZNYWGh4uLidPLkSdWpU0chISFKTU1VdHS0q04BAAC4GZdPKJ44caImTpxY7rrk5GSH5RdeeEEvvPBCNVQFAABqKpe/fgEAAKAyEW4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpEG4AAICpuDzcvP3222rVqpW8vLwUGhqq7du337B/enq6QkND5eXlpXvuuUdJSUnVVCkAAKgJXBpu1qxZo+eee04vvviiMjMzFRERof79+ysrK6vc/seOHVN0dLQiIiKUmZmpadOmafLkyUpJSanmygEAgLtyabiZM2eOxo0bp/Hjx6tdu3aaO3eumjVrpsTExHL7JyUlqXnz5po7d67atWun8ePHa+zYsZo1a1Y1Vw4AANyVy8JNYWGhdu/eraioKIf2qKgoffHFF+Vus2PHjjL9+/Xrp4yMDBUVFVVZrQAAoObwcNWBz549q5KSEgUGBjq0BwYGKjc3t9xtcnNzy+1fXFyss2fPKigoqMw2NptNNpvNvpyXlydJys/Pv91TcErhrxddclxXc9X37Wpc7zsL1/vOwvV2zXENw7hpX5eFmyssFovDsmEYZdpu1r+89isSEhI0ffr0Mu3NmjWraKm4DSsmuroCVCeu952F631ncfX1vnDhgurXr3/DPi4LN40aNVKtWrXKjNKcPn26zOjMFU2aNCm3v4eHh/z9/cvdJj4+XrGxsfbl0tJSnT9/Xv7+/jcMUWaTn5+vZs2aKTs7W76+vq4uB1WM631n4XrfWe7U620Yhi5cuKDg4OCb9nVZuPH09FRoaKi2bNmiwYMH29u3bNmiQYMGlbtNeHi4NmzY4NC2efNmhYWFqXbt2uVuY7VaZbVaHdoaNGhwe8XXYL6+vnfUvwx3Oq73nYXrfWe5E6/3zUZsrnDp01KxsbF65513tHjxYh06dEhTp05VVlaWJkyYIOnyqMuoUaPs/SdMmKATJ04oNjZWhw4d0uLFi7Vo0SLFxcW56hQAAICbcemcm8cff1znzp3TX/7yF+Xk5KhDhw7auHGjWrRoIUnKyclx+M2bVq1aaePGjZo6dareeustBQcHa968eRo6dKirTgEAALgZl08onjhxoiZOLH92UnJycpm2Xr16ac+ePVVclflYrVa99NJLZW7RwZy43ncWrvedhet9cxbjVp6pAgAAqCFc/m4pAACAykS4AQAApkK4AQAApkK4AQAApkK4uYMwdxwAcCcg3NxBrFarDh065OoyAACoUi7/nRtUvqvfpXW1kpISzZgxw/4erjlz5lRnWahCv/76q3bv3i0/Pz+1b9/eYd2lS5e0du1ah1/7Rs2VmZmpBg0aqFWrVpKkFStWKDExUVlZWWrRooWeeeYZDR8+3MVVoqr89NNPWrp0qY4ePaqgoCCNHj2aF0GXg9+5MaG77rpLnTp1KvMOrfT0dIWFhalu3bqyWCzaunWrawpEpTpy5IiioqKUlZUli8WiiIgIrVq1SkFBQZKkU6dOKTg4WCUlJS6uFJWha9eumj17tvr06aN33nlHkydP1tNPP6127drp8OHDeuedd/TGG29o7Nixri4VlSA4OFgHDhyQv7+/jh07pu7du0uSHnjgAR06dEgXLlzQzp071bZtWxdX6l4INyaUkJCgf/zjH3rnnXf0yCOP2Ntr166tffv2lfk/e9RsgwcPVnFxsZYsWaKff/5ZsbGxOnjwoNLS0tS8eXPCjcnUrVtXhw4dUvPmzdW1a1dNmDBBv//97+3rV65cqb/+9a/66quvXFglKstdd92l3NxcBQQE6He/+51yc3OVmpoqb29v2Ww2DRs2TF5eXnrvvfdcXapbYc6NCcXHx2vNmjX64x//qLi4OBUVFbm6JFShL774Qq+99poaNWqk1q1ba/369erfv78iIiL0/fffu7o8VLI6derozJkzkqSTJ0/qwQcfdFj/4IMP6tixY64oDVXs//7v//TnP/9Z3t7eki7Po/zTn/6knTt3urgy90O4Malu3bpp9+7dOnPmjMLCwnTgwAFZLBZXl4Uq8Ouvv8rDw3H63FtvvaWBAweqV69eOnLkiIsqQ1Xo37+/EhMTJV1+197777/vsH7t2rVq3bq1K0pDFbny326bzabAwECHdYGBgfawi/+PCcUmVq9ePS1dulSrV69W3759uS1hUm3btlVGRobatWvn0P7mm2/KMAwNHDjQRZWhKsycOVMPP/ywevXqpbCwMM2ePVtpaWn2OTc7d+7UunXrXF0mKlFkZKQ8PDyUn5+vI0eOKCQkxL4uKytLjRo1cmF17olwcwcYPny4evTood27d6tFixauLgeVbPDgwVq1apVGjhxZZt38+fNVWlqqpKQkF1SGqhAcHKzMzEzNmDFDGzZskGEY+vLLL5Wdna2HH35Yn3/+ucLCwlxdJirJSy+95LB85ZbUFRs2bFBERER1llQjMKEYAACYCnNuAACAqRBuAACAqRBuAACAqRBuAJhGy5YtNXfuXFeXAcDFCDcAAMBUCDcAAMBUCDcA3MKCBQt09913q7S01KF94MCBGj16tL777jsNGjRIgYGBqlevnrp166aPP/74uvs7fvy4LBaL9u7da2/7+eefZbFYlJaWZm/7+uuvFR0drXr16ikwMFAjR47U2bNnK/v0AFQjwg0At/Cf//mfOnv2rLZt22Zv++mnn/Tvf/9bI0aM0MWLFxUdHa2PP/5YmZmZ6tevn2JiYpSVleX0MXNyctSrVy917txZGRkZ2rRpk06dOqX/+q//qoxTAuAi/EIxALfg5+enxx57TCtXrlRkZKQk6b333pOfn58iIyNVq1YtderUyd7/1Vdf1bp167R+/Xo988wzTh0zMTFRXbt21WuvvWZvW7x4sZo1a6YjR47o/vvvv72TAuASjNwAcBsjRoxQSkqKbDabJOndd9/V8OHDVatWLRUUFOiFF15Q+/bt1aBBA9WrV0/ffPPNbY3c7N69W9u2bVO9evXsn7Zt20qSvvvuu0o5JwDVj5EbAG4jJiZGpaWlSk1NVbdu3bR9+3bNmTNHkvTf//3f+ve//61Zs2apdevWqlOnjoYNG6bCwsJy93XXXZf/3+3qN8wUFRU59CktLVVMTIxmzpxZZvugoKDKOi0A1YxwA8Bt1KlTR0OGDNG7776rb7/9Vvfff79CQ0MlSdu3b9eYMWM0ePBgSdLFixd1/Pjx6+6rcePGki7Pq+nSpYskOUwulqSuXbsqJSVFLVu2lIcH/zkEzILbUgDcyogRI5SamqrFixfrySeftLe3bt1aH3zwgfbu3at9+/bpiSeeKPNk1dXq1Kmjhx56SDNmzNDXX3+tTz/9VH/6058c+kyaNEnnz5/X7373O3355Zf6/vvvtXnzZo0dO1YlJSVVdo4AqhbhBoBbeeSRR+Tn56fDhw/riSeesLf//e9/V8OGDdW9e3fFxMSoX79+6tq16w33tXjxYhUVFSksLExTpkzRq6++6rA+ODhYn3/+uUpKStSvXz916NBBU6ZMUf369e23tQDUPBbj6hvSAAAANRz/awIAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEyFcAMAAEzl/wFwxbjt1Il4RgAAAABJRU5ErkJggg==",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# plotting FD\n",
+ "frequency_distribution.plot(kind='bar', alpha=0.7)\n",
+ "plt.xlabel('value')\n",
+ "plt.ylabel('Frequency')\n",
+ "plt.title('Frequency Distribution')\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "- Comment : Horizontal and vertical axes are interchanged. "
]
},
{
@@ -91,11 +361,16 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
- "# your code here"
+ "def calculate_mean(input_series):\n",
+ " count = 0\n",
+ " \n",
+ " for key,value in input_series.items():\n",
+ " count += key*value\n",
+ " return count/sum(input_series.values)"
]
},
{
@@ -107,11 +382,48 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 9,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "4 4\n",
+ "2 3\n",
+ "6 2\n",
+ "5 1\n",
+ "Name: value, dtype: int64"
+ ]
+ },
+ "execution_count": 9,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "# calculating FD\n",
+ "frequency_distribution = result['value'].value_counts()\n",
+ "frequency_distribution"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "3.9"
+ ]
+ },
+ "execution_count": 10,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_mean(frequency_distribution)"
]
},
{
@@ -124,11 +436,44 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 11,
"metadata": {},
"outputs": [],
"source": [
- "# your code here"
+ "def calculate_median(input_series):\n",
+ "\n",
+ " values_list = input_series.index.tolist()\n",
+ " values_list.sort()\n",
+ "\n",
+ " n = len(values_list)\n",
+ " if n % 2 == 0:\n",
+ " middle1 = values_list[n // 2 - 1]\n",
+ " middle2 = values_list[n // 2]\n",
+ " median = (middle1 + middle2) / 2\n",
+ " else:\n",
+ " median = values_list[n // 2]\n",
+ "\n",
+ " return median"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "4.5"
+ ]
+ },
+ "execution_count": 12,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_median(frequency_distribution)"
]
},
{
@@ -140,11 +485,37 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 13,
"metadata": {},
"outputs": [],
"source": [
- "# your code here"
+ "def calculate_quartiles(df, column):\n",
+ " \n",
+ " values_list = df[column].to_list()\n",
+ " values_list.sort()\n",
+ " \n",
+ " q75, q50, q25 = np.percentile(values_list, [75, 50, 25])\n",
+ " return q75, q50, q25"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "(4.75, 4.0, 2.5)"
+ ]
+ },
+ "execution_count": 14,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "calculate_quartiles(result, \"value\")"
]
},
{
@@ -158,22 +529,251 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 15,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " Unnamed: 0 \n",
+ " roll \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 0 \n",
+ " 0 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 1 \n",
+ " 1 \n",
+ " 2 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 2 \n",
+ " 2 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 3 \n",
+ " 3 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 4 \n",
+ " 4 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " Unnamed: 0 \n",
+ " roll \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 0 \n",
+ " 0 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 47 \n",
+ " 47 \n",
+ " 47 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 56 \n",
+ " 56 \n",
+ " 56 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 9 \n",
+ " 9 \n",
+ " 9 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 73 \n",
+ " 73 \n",
+ " 73 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 17 \n",
+ " 17 \n",
+ " 17 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 11 \n",
+ " 11 \n",
+ " 11 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 24 \n",
+ " 24 \n",
+ " 24 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 21 \n",
+ " 21 \n",
+ " 21 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 99 \n",
+ " 99 \n",
+ " 99 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
100 rows × 3 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# plotting the value\n",
+ "dice_hundred[\"value\"].plot(kind='bar', alpha=0.7)\n",
+ "plt.xlabel('roll(index)')\n",
+ "plt.ylabel('value')\n",
+ "plt.title('')\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : The number of times that each value appears (the frequency) differs by value (ex. 5 less vs. 4 and 6)"
]
},
{
@@ -185,11 +785,23 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 18,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "3.74"
+ ]
+ },
+ "execution_count": 18,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "frequency_distribution = dice_hundred['value'].value_counts()\n",
+ "calculate_mean(frequency_distribution)"
]
},
{
@@ -201,11 +813,29 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 19,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "6 23\n",
+ "4 22\n",
+ "2 17\n",
+ "3 14\n",
+ "1 12\n",
+ "5 12\n",
+ "Name: value, dtype: int64"
+ ]
+ },
+ "execution_count": 19,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "# already calculated in above cell\n",
+ "frequency_distribution"
]
},
{
@@ -217,22 +847,33 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 20,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAASIAAAEiCAYAAABdvt+2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAASuUlEQVR4nO3df2zU9R3H8de37e0rdXeYgv0Vuq6amrnVX6MOB05KtlYbY4L4K2u2gZlOY2HDi2EwQjymtgvJCFmITdgfqH808I8ocWh7ZlpUgkI3o6JTjFU7pGMg65XWHAf33R+kN2u7tle+X970fD6S/nHffve59312ffrttT0cz/M8AYChPOsBAIAQATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwFyB9QBflU6n9dlnnykcDstxHOtxAEyR53kaGBhQeXm58vLGv+Y570L02WefqaKiwnoMAD7p7e3VnDlzxj3nvAtROByWdGb4SCRiPI3/UqmUOjs71dDQoFAoZD1OzmF/g5XN/iYSCVVUVGS+psdz3oVo+NuxSCSSsyEqLCxUJBLhCyUA7G+wprK/k3mJhRerAZgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYO+9+sxqAP769+i++r+nme9rwA9+X5YoIgD1CBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5frMaOakm1qHkaX//OaqP/3Czr+vhf7giAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIC5rELU2tqqa6+9VuFwWMXFxVq8eLHef//9Eed4nqdYLKby8nLNmDFDdXV1OnDggK9DA8gtWYWoq6tLzc3N2rt3r+LxuE6dOqWGhgYNDg5mztmwYYM2btyozZs3a9++fSotLVV9fb0GBgZ8Hx5AbsjqHRpfeOGFEbe3bt2q4uJidXd364YbbpDnedq0aZPWrl2rJUuWSJKefPJJlZSUqL29Xffdd59/kwPIGWf1VrH9/f2SpKKiIklST0+P+vr61NDQkDnHdV0tXLhQe/bsGTNEyWRSyWQyczuRSEiSUqmUUqnU2Yx3Xhp+TLn42M4Hw/vq5nmBrT1duPn+78Hwvk5mL7LZrymHyPM8RaNRXX/99aqpqZEk9fX1SZJKSkpGnFtSUqJPPvlkzHVaW1u1fv36Ucc7OztVWFg41fHOe/F43HqEnPZIbdr3NXft2uX7mkHa8IPg1p7M83doaGjS6005RMuXL9dbb72lV199ddTnHGfkm5Z7njfq2LA1a9YoGo1mbicSCVVUVKihoUGRSGSq4523UqmU4vG46uvrFQqFrMfJOcP7u25/npJpf988/53Yjb6uF7SaWIfva7p5nh6pTU/q+Tv83c1kTClEK1as0M6dO7V7927NmTMnc7y0tFTSmSujsrKyzPEjR46Mukoa5rquXNcddTwUCuX0F2quPz5rybTj+7/iMd3+//L78X/ZZJ6/2exXVj818zxPy5cv19NPP62//vWvqqqqGvH5qqoqlZaWjrhsO3nypLq6ujR//vxs7grA10hWV0TNzc1qb2/Xs88+q3A4nHlNaObMmZoxY4Ycx9HKlSvV0tKi6upqVVdXq6WlRYWFhWpqagrkAQCY/rIKUVtbmySprq5uxPGtW7dq2bJlkqRVq1bpiy++0AMPPKDjx49r3rx56uzsVDgc9mVgALknqxB53sQ/DnQcR7FYTLFYbKozAfia4W/NAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMJfVP7AI/9TEOpQ87fi65sd/uNnX9YBzhSsiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJjLOkS7d+/WLbfcovLycjmOo2eeeWbE55ctWybHcUZ8XHfddX7NCyAHZR2iwcFBXXXVVdq8efP/Peemm27S4cOHMx+7du06qyEB5Las36GxsbFRjY2N457juq5KS0unPBSAr5dA3ir25ZdfVnFxsS666CItXLhQjz32mIqLi8c8N5lMKplMZm4nEglJUiqVUiqVCmI8U8OPyc3zAlv764z9/R833/89GN7XyexFNvvleJ435Wkdx9GOHTu0ePHizLHt27frm9/8piorK9XT06N169bp1KlT6u7uluu6o9aIxWJav379qOPt7e0qLCyc6mgAjA0NDampqUn9/f2KRCLjnut7iL7q8OHDqqys1LZt27RkyZJRnx/riqiiokJHjx6dcPjpKJVKKR6Pa93+PCXT/r55/juxG31dbzpif/+nJtbh+5punqdHatOqr69XKBQa99xEIqHZs2dPKkSB/yseZWVlqqys1MGDB8f8vOu6Y14phUKhCR/odJZMO77/Kx65vF/ZYn/l++P/ssl8fWazX4H/HtGxY8fU29ursrKyoO8KwDSV9RXRiRMn9OGHH2Zu9/T06M0331RRUZGKiooUi8V02223qaysTB9//LF+97vfafbs2br11lt9HRxA7sg6RPv379eiRYsyt6PRqCRp6dKlamtr09tvv62nnnpK//nPf1RWVqZFixZp+/btCofD/k0NIKdkHaK6ujqN9/p2R4f/L5AByG38rRkAc4QIgDlCBMAcIQJgLvBfaMT09+3Vfwlk3Y//cHMg62L64YoIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmBu2r9ndVDvpyzxnso4N4J8Dk8XXBEBMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwFzWIdq9e7duueUWlZeXy3EcPfPMMyM+73meYrGYysvLNWPGDNXV1enAgQN+zQsgB2UdosHBQV111VXavHnzmJ/fsGGDNm7cqM2bN2vfvn0qLS1VfX29BgYGznpYALkp63dobGxsVGNj45if8zxPmzZt0tq1a7VkyRJJ0pNPPqmSkhK1t7frvvvuO7tpAeQkX98qtqenR319fWpoaMgcc11XCxcu1J49e8YMUTKZVDKZzNxOJBKSpFQqpVQqNeF9uvmeD5OPbTL3P9U13Tz/5w5iXim4PWZ/zwjyOey34X2dzF5ks1+O53lT3gXHcbRjxw4tXrxYkrRnzx4tWLBAhw4dUnl5eea8X/3qV/rkk0/U0dExao1YLKb169ePOt7e3q7CwsKpjgbA2NDQkJqamtTf369IJDLuuYG8eb7jOCNue5436tiwNWvWKBqNZm4nEglVVFSooaFhwuElqSY2Om5+eSd2o+9rplIpxeNxrdufp2R67D2ZqiDmlYLbY/b3jCCfw35z8zw9UptWfX29QqHQuOcOf3czGb6GqLS0VJLU19ensrKyzPEjR46opKRkzP+N67pyXXfU8VAoNOEDlaTkaX+fbF+dISjJtOP77EHNG9Qes79nBPkcDspkvj6z2S9ff4+oqqpKpaWlisfjmWMnT55UV1eX5s+f7+ddAcghWV8RnThxQh9++GHmdk9Pj958800VFRXpW9/6llauXKmWlhZVV1erurpaLS0tKiwsVFNTk6+DA8gdWYdo//79WrRoUeb28Os7S5cu1RNPPKFVq1bpiy++0AMPPKDjx49r3rx56uzsVDgc9m9qADkl6xDV1dVpvB+0OY6jWCymWCx2NnMB+Brhb80AmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATDne4hisZgcxxnxUVpa6vfdAMghBUEs+r3vfU8vvvhi5nZ+fn4QdwMgRwQSooKCAq6CAExaICE6ePCgysvL5bqu5s2bp5aWFl1yySVjnptMJpVMJjO3E4mEJCmVSimVSk14X26+58/QY5jM/U91TTfP/7mDmFcKbo/Z3zOCfA77bXhfJ7MX2eyX43mer7vw/PPPa2hoSJdddpn+9a9/6dFHH9U//vEPHThwQLNmzRp1fiwW0/r160cdb29vV2FhoZ+jATiHhoaG1NTUpP7+fkUikXHP9T1EXzU4OKhLL71Uq1atUjQaHfX5sa6IKioqdPTo0QmHl6SaWIev837ZO7EbfV8zlUopHo9r3f48JdOOr2sHMa8U3B6zv2cE+Rz2m5vn6ZHatOrr6xUKhcY9N5FIaPbs2ZMKUSDfmn3ZhRdeqCuuuEIHDx4c8/Ou68p13VHHQ6HQhA9UkpKn/X2yfXWGoCTTju+zBzVvUHvM/p4R5HM4KJP5+sxmvwL/PaJkMqn33ntPZWVlQd8VgGnK9xA99NBD6urqUk9Pj15//XXdfvvtSiQSWrp0qd93BSBH+P6t2T//+U/99Kc/1dGjR3XxxRfruuuu0969e1VZWen3XQHIEb6HaNu2bX4vCSDH8bdmAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAOUIEwBwhAmCOEAEwR4gAmCNEAMwRIgDmCBEAc4QIgDlCBMAcIQJgjhABMEeIAJgjRADMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYCyxEjz/+uKqqqnTBBRdo7ty5euWVV4K6KwDTXCAh2r59u1auXKm1a9fq73//u370ox+psbFRn376aRB3B2CaCyREGzdu1C9/+Uvdc889uvzyy7Vp0yZVVFSora0tiLsDMM0V+L3gyZMn1d3drdWrV4843tDQoD179ow6P5lMKplMZm739/dLkj7//HOlUqkJ76/g1OBZTvz/HTt2zPc1U6mUhoaGVJDK0+m04+vaQcwrBbfH7O8ZQT6H/VaQ9jQ0lNaxY8cUCoXGPXdgYECS5HnexAt7Pjt06JAnyXvttddGHH/ssce8yy67bNT5Dz/8sCeJDz74yNGP3t7eCbvh+xXRMMcZ+V8jz/NGHZOkNWvWKBqNZm6n02l9/vnnmjVr1pjnT3eJREIVFRXq7e1VJBKxHifnsL/BymZ/Pc/TwMCAysvLJ1zX9xDNnj1b+fn56uvrG3H8yJEjKikpGXW+67pyXXfEsYsuusjvsc47kUiEL5QAsb/Bmuz+zpw5c1Lr+f5i9Te+8Q3NnTtX8Xh8xPF4PK758+f7fXcAckAg35pFo1H9/Oc/V21trX74wx9qy5Yt+vTTT3X//fcHcXcAprlAQnTXXXfp2LFj+v3vf6/Dhw+rpqZGu3btUmVlZRB3N624rquHH3541Lej8Af7G6yg9tfxvMn8bA0AgsPfmgEwR4gAmCNEAMwRIgDmCNE50NraqmuvvVbhcFjFxcVavHix3n//feuxclZra6scx9HKlSutR8kphw4d0s9+9jPNmjVLhYWFuvrqq9Xd3e3L2oToHOjq6lJzc7P27t2reDyuU6dOqaGhQYOD0+ePHaeLffv2acuWLbryyiutR8kpx48f14IFCxQKhfT888/r3Xff1R//+Eff/gqCH98b+Pe//63i4mJ1dXXphhtusB4nZ5w4cULf//739fjjj+vRRx/V1VdfrU2bNlmPlRNWr16t1157LbA3OOSKyMDwW50UFRUZT5JbmpubdfPNN+snP/mJ9Sg5Z+fOnaqtrdUdd9yh4uJiXXPNNfrzn//s2/qE6BzzPE/RaFTXX3+9ampqrMfJGdu2bdPf/vY3tba2Wo+Skz766CO1tbWpurpaHR0duv/++/XrX/9aTz31lC/rB/Y2IBjb8uXL9dZbb+nVV1+1HiVn9Pb26je/+Y06Ozt1wQUXWI+Tk9LptGpra9XS0iJJuuaaa3TgwAG1tbXpF7/4xVmvzxXRObRixQrt3LlTL730kubMmWM9Ts7o7u7WkSNHNHfuXBUUFKigoEBdXV3605/+pIKCAp0+fdp6xGmvrKxM3/3ud0ccu/zyy317H3quiM4Bz/O0YsUK7dixQy+//LKqqqqsR8opP/7xj/X222+POHb33XfrO9/5jn77298qPz/faLLcsWDBglG/cvLBBx/49ofshOgcaG5uVnt7u5599lmFw+HMm8bNnDlTM2bMMJ5u+guHw6Neb7vwwgs1a9YsXofzyYMPPqj58+erpaVFd955p9544w1t2bJFW7Zs8WV9fnx/Dvy/t7zdunWrli1bdm6H+Zqoq6vjx/c+e+6557RmzRodPHhQVVVVikajuvfee31ZmxABMMeL1QDMESIA5ggRAHOECIA5QgTAHCECYI4QATBHiACYI0QAzBEiAOYIEQBzhAiAuf8CkE34rsQ/UD0AAAAASUVORK5CYII=",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
"source": [
- "# your code here"
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "dice_hundred[\"value\"].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show() "
]
},
{
- "cell_type": "code",
- "execution_count": null,
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "- Comment : As larger numbers (4 and 6) appears more frequently compared to smaller values (2 and 3), mean is larger than the median (1 and 5 are found with same frequency, meaning they make no difference)."
]
},
{
@@ -244,22 +885,119 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 21,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " Unnamed: 0 \n",
+ " roll \n",
+ " value \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 0 \n",
+ " 0 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 1 \n",
+ " 1 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 2 \n",
+ " 2 \n",
+ " 1 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 3 \n",
+ " 3 \n",
+ " 6 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 4 \n",
+ " 4 \n",
+ " 5 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "dice_thousand[\"value\"].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show() "
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comments : The frequency distribution is more consistent for all values. It becomes more stable as the number of rolls increases."
]
},
{
@@ -274,11 +1012,142 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 23,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " observation \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 68.0 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 12.0 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 45.0 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 38.0 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 49.0 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 995 \n",
+ " 27.0 \n",
+ " \n",
+ " \n",
+ " 996 \n",
+ " 47.0 \n",
+ " \n",
+ " \n",
+ " 997 \n",
+ " 53.0 \n",
+ " \n",
+ " \n",
+ " 998 \n",
+ " 33.0 \n",
+ " \n",
+ " \n",
+ " 999 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
1000 rows × 1 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "population['observation'].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Assumption\n",
+ "- range of mean : 30-40\n",
+ "- std : 15 "
]
},
{
@@ -290,22 +1159,51 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 25,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "36.56"
+ ]
+ },
+ "execution_count": 25,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "np.mean(population['observation'])"
]
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 26,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "12.81008977329979"
+ ]
+ },
+ "execution_count": 26,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "np.std(population['observation'])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "Diffrence vs. assumption\n",
+ "- mean: falls inside the range assumed\n",
+ "- std: lower than assumption"
]
},
{
@@ -317,11 +1215,133 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 27,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " observation \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 25.0 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 29.0 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 29.0 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 995 \n",
+ " 26.0 \n",
+ " \n",
+ " \n",
+ " 996 \n",
+ " 22.0 \n",
+ " \n",
+ " \n",
+ " 997 \n",
+ " 21.0 \n",
+ " \n",
+ " \n",
+ " 998 \n",
+ " 19.0 \n",
+ " \n",
+ " \n",
+ " 999 \n",
+ " 28.0 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
1000 rows × 1 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "population_2['observation'].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show()"
]
},
{
@@ -332,14 +1352,10 @@
]
},
{
- "cell_type": "code",
- "execution_count": null,
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "- Comment : Distribution is centered in the range of 26-29 and std is assumed to be smaller compared to that of step 2 (meaning most of the values are found between smaller range of values)."
]
},
{
@@ -351,22 +1367,49 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 29,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "27.155"
+ ]
+ },
+ "execution_count": 29,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "np.mean(population_2['observation'])"
]
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 30,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "2.9683286543103704"
+ ]
+ },
+ "execution_count": 30,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "np.std(population_2['observation'])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : Same to above."
]
},
{
@@ -381,11 +1424,133 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 31,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/html": [
+ "\n",
+ "\n",
+ "
\n",
+ " \n",
+ " \n",
+ " observation \n",
+ " \n",
+ " \n",
+ " \n",
+ " \n",
+ " 0 \n",
+ " 21.0 \n",
+ " \n",
+ " \n",
+ " 1 \n",
+ " 21.0 \n",
+ " \n",
+ " \n",
+ " 2 \n",
+ " 24.0 \n",
+ " \n",
+ " \n",
+ " 3 \n",
+ " 31.0 \n",
+ " \n",
+ " \n",
+ " 4 \n",
+ " 54.0 \n",
+ " \n",
+ " \n",
+ " ... \n",
+ " ... \n",
+ " \n",
+ " \n",
+ " 995 \n",
+ " 16.0 \n",
+ " \n",
+ " \n",
+ " 996 \n",
+ " 55.0 \n",
+ " \n",
+ " \n",
+ " 997 \n",
+ " 30.0 \n",
+ " \n",
+ " \n",
+ " 998 \n",
+ " 35.0 \n",
+ " \n",
+ " \n",
+ " 999 \n",
+ " 43.0 \n",
+ " \n",
+ " \n",
+ "
\n",
+ "
1000 rows × 1 columns
\n",
+ "
"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# Plotting into histogram\n",
+ "rcParams['figure.figsize'] = 3, 3 \n",
+ "population_3['observation'].hist();\n",
+ "plt.tight_layout() \n",
+ "plt.show()"
]
},
{
@@ -397,22 +1562,49 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 33,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "41.989"
+ ]
+ },
+ "execution_count": 33,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "np.mean(population_3['observation'])"
]
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 34,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "16.136631587788084"
+ ]
+ },
+ "execution_count": 34,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "np.std(population_3['observation'])"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "- Comment : Mean is larger than mode and it is assume to be positively skewed. Distribution is varied in wider range resulting in larger std."
]
},
{
@@ -424,22 +1616,29 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 35,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "(53.0, 40.0, 30.0)"
+ ]
+ },
+ "execution_count": 35,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "calculate_quartiles(population_3, \"observation\")"
]
},
{
- "cell_type": "code",
- "execution_count": null,
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "- Comment : mean is slightly larger than median and the difference between median(q2) and q3 is larger than that of q1 and median, meaning it is positive distribution."
]
},
{
@@ -451,22 +1650,32 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 36,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "array([22., 67.])"
+ ]
+ },
+ "execution_count": 36,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
- "# your code here"
+ "values_list = population_3[\"observation\"].to_list()\n",
+ "values_list.sort()\n",
+ " \n",
+ "np.percentile(values_list, [10, 90])"
]
},
{
- "cell_type": "code",
- "execution_count": null,
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
- "\"\"\"\n",
- "your comments here\n",
- "\"\"\""
+ "- Comment : If we see the 10th and 90th percentile and compared the difference between the median, we find that the difference vs. 90th(67-40=23) is larger than that vs. 10th (40-22=18) , meaning it is positively skewed."
]
},
{
@@ -479,7 +1688,7 @@
},
{
"cell_type": "code",
- "execution_count": null,
+ "execution_count": 37,
"metadata": {},
"outputs": [],
"source": [
@@ -487,10 +1696,8 @@
]
},
{
- "cell_type": "code",
- "execution_count": null,
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
"\"\"\"\n",
"your comments here\n",
@@ -500,9 +1707,9 @@
],
"metadata": {
"kernelspec": {
- "display_name": "ironhack-3.7",
+ "display_name": "Python 3 (ipykernel)",
"language": "python",
- "name": "ironhack-3.7"
+ "name": "python3"
},
"language_info": {
"codemirror_mode": {
@@ -514,7 +1721,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
- "version": "3.7.3"
+ "version": "3.11.3"
}
},
"nbformat": 4,
diff --git a/your-code/roll_the_dice_hundred.csv b/your-code/roll_the_dice_hundred.csv
new file mode 100644
index 0000000..50975a2
--- /dev/null
+++ b/your-code/roll_the_dice_hundred.csv
@@ -0,0 +1,101 @@
+,roll,value
+0,0,1
+1,1,2
+2,2,6
+3,3,1
+4,4,6
+5,5,5
+6,6,2
+7,7,2
+8,8,4
+9,9,1
+10,10,5
+11,11,6
+12,12,5
+13,13,4
+14,14,5
+15,15,4
+16,16,4
+17,17,6
+18,18,2
+19,19,4
+20,20,4
+21,21,6
+22,22,3
+23,23,6
+24,24,6
+25,25,4
+26,26,1
+27,27,4
+28,28,4
+29,29,2
+30,30,6
+31,31,5
+32,32,5
+33,33,2
+34,34,3
+35,35,6
+36,36,6
+37,37,2
+38,38,3
+39,39,6
+40,40,6
+41,41,2
+42,42,5
+43,43,3
+44,44,4
+45,45,6
+46,46,2
+47,47,1
+48,48,4
+49,49,2
+50,50,3
+51,51,2
+52,52,2
+53,53,4
+54,54,6
+55,55,2
+56,56,1
+57,57,3
+58,58,2
+59,59,4
+60,60,4
+61,61,3
+62,62,4
+63,63,1
+64,64,3
+65,65,6
+66,66,3
+67,67,4
+68,68,4
+69,69,4
+70,70,2
+71,71,2
+72,72,5
+73,73,1
+74,74,5
+75,75,6
+76,76,2
+77,77,4
+78,78,6
+79,79,5
+80,80,6
+81,81,4
+82,82,1
+83,83,3
+84,84,3
+85,85,3
+86,86,5
+87,87,6
+88,88,5
+89,89,1
+90,90,6
+91,91,3
+92,92,6
+93,93,4
+94,94,1
+95,95,4
+96,96,6
+97,97,1
+98,98,3
+99,99,6
diff --git a/your-code/roll_the_dice_thousand.csv b/your-code/roll_the_dice_thousand.csv
new file mode 100644
index 0000000..f820dbb
--- /dev/null
+++ b/your-code/roll_the_dice_thousand.csv
@@ -0,0 +1,1001 @@
+,roll,value
+0,0,5
+1,1,6
+2,2,1
+3,3,6
+4,4,5
+5,5,2
+6,6,6
+7,7,5
+8,8,6
+9,9,6
+10,10,4
+11,11,3
+12,12,5
+13,13,6
+14,14,1
+15,15,3
+16,16,1
+17,17,1
+18,18,1
+19,19,1
+20,20,6
+21,21,2
+22,22,3
+23,23,4
+24,24,6
+25,25,5
+26,26,3
+27,27,2
+28,28,4
+29,29,1
+30,30,3
+31,31,4
+32,32,3
+33,33,3
+34,34,6
+35,35,2
+36,36,1
+37,37,2
+38,38,6
+39,39,4
+40,40,1
+41,41,4
+42,42,6
+43,43,1
+44,44,6
+45,45,3
+46,46,6
+47,47,4
+48,48,5
+49,49,1
+50,50,4
+51,51,4
+52,52,4
+53,53,6
+54,54,2
+55,55,6
+56,56,4
+57,57,6
+58,58,6
+59,59,2
+60,60,4
+61,61,1
+62,62,4
+63,63,3
+64,64,1
+65,65,4
+66,66,1
+67,67,2
+68,68,3
+69,69,3
+70,70,1
+71,71,6
+72,72,1
+73,73,6
+74,74,5
+75,75,2
+76,76,6
+77,77,5
+78,78,4
+79,79,6
+80,80,5
+81,81,1
+82,82,4
+83,83,4
+84,84,4
+85,85,4
+86,86,6
+87,87,2
+88,88,4
+89,89,5
+90,90,4
+91,91,3
+92,92,4
+93,93,3
+94,94,2
+95,95,5
+96,96,1
+97,97,6
+98,98,6
+99,99,1
+100,100,1
+101,101,3
+102,102,3
+103,103,2
+104,104,1
+105,105,3
+106,106,6
+107,107,2
+108,108,1
+109,109,4
+110,110,1
+111,111,2
+112,112,2
+113,113,4
+114,114,1
+115,115,2
+116,116,2
+117,117,4
+118,118,4
+119,119,1
+120,120,4
+121,121,5
+122,122,3
+123,123,2
+124,124,5
+125,125,5
+126,126,1
+127,127,4
+128,128,3
+129,129,6
+130,130,5
+131,131,1
+132,132,5
+133,133,1
+134,134,1
+135,135,2
+136,136,6
+137,137,6
+138,138,5
+139,139,5
+140,140,1
+141,141,5
+142,142,1
+143,143,1
+144,144,6
+145,145,6
+146,146,6
+147,147,6
+148,148,5
+149,149,1
+150,150,4
+151,151,6
+152,152,5
+153,153,4
+154,154,1
+155,155,4
+156,156,5
+157,157,5
+158,158,1
+159,159,2
+160,160,6
+161,161,5
+162,162,2
+163,163,3
+164,164,1
+165,165,2
+166,166,1
+167,167,6
+168,168,4
+169,169,5
+170,170,6
+171,171,4
+172,172,3
+173,173,1
+174,174,1
+175,175,5
+176,176,2
+177,177,5
+178,178,2
+179,179,3
+180,180,3
+181,181,4
+182,182,6
+183,183,3
+184,184,2
+185,185,4
+186,186,4
+187,187,3
+188,188,6
+189,189,2
+190,190,4
+191,191,2
+192,192,1
+193,193,4
+194,194,2
+195,195,3
+196,196,2
+197,197,1
+198,198,6
+199,199,4
+200,200,1
+201,201,5
+202,202,4
+203,203,3
+204,204,6
+205,205,3
+206,206,3
+207,207,6
+208,208,3
+209,209,2
+210,210,6
+211,211,2
+212,212,4
+213,213,1
+214,214,1
+215,215,2
+216,216,5
+217,217,5
+218,218,3
+219,219,5
+220,220,2
+221,221,1
+222,222,1
+223,223,2
+224,224,5
+225,225,4
+226,226,6
+227,227,5
+228,228,5
+229,229,5
+230,230,3
+231,231,2
+232,232,2
+233,233,5
+234,234,2
+235,235,3
+236,236,4
+237,237,2
+238,238,6
+239,239,2
+240,240,6
+241,241,3
+242,242,6
+243,243,5
+244,244,3
+245,245,2
+246,246,4
+247,247,6
+248,248,3
+249,249,6
+250,250,4
+251,251,5
+252,252,3
+253,253,5
+254,254,1
+255,255,2
+256,256,6
+257,257,1
+258,258,4
+259,259,4
+260,260,1
+261,261,1
+262,262,3
+263,263,3
+264,264,1
+265,265,3
+266,266,3
+267,267,4
+268,268,1
+269,269,2
+270,270,3
+271,271,3
+272,272,3
+273,273,1
+274,274,6
+275,275,6
+276,276,3
+277,277,3
+278,278,4
+279,279,4
+280,280,6
+281,281,5
+282,282,3
+283,283,6
+284,284,4
+285,285,3
+286,286,6
+287,287,6
+288,288,1
+289,289,5
+290,290,2
+291,291,1
+292,292,3
+293,293,6
+294,294,6
+295,295,6
+296,296,1
+297,297,3
+298,298,6
+299,299,5
+300,300,2
+301,301,2
+302,302,3
+303,303,4
+304,304,6
+305,305,4
+306,306,5
+307,307,5
+308,308,6
+309,309,6
+310,310,4
+311,311,4
+312,312,6
+313,313,2
+314,314,5
+315,315,5
+316,316,5
+317,317,5
+318,318,6
+319,319,1
+320,320,3
+321,321,3
+322,322,6
+323,323,6
+324,324,1
+325,325,6
+326,326,6
+327,327,4
+328,328,3
+329,329,4
+330,330,5
+331,331,6
+332,332,4
+333,333,1
+334,334,2
+335,335,3
+336,336,4
+337,337,3
+338,338,2
+339,339,3
+340,340,4
+341,341,3
+342,342,6
+343,343,3
+344,344,2
+345,345,6
+346,346,3
+347,347,3
+348,348,6
+349,349,2
+350,350,3
+351,351,1
+352,352,4
+353,353,2
+354,354,2
+355,355,6
+356,356,3
+357,357,2
+358,358,5
+359,359,4
+360,360,6
+361,361,2
+362,362,3
+363,363,3
+364,364,6
+365,365,5
+366,366,1
+367,367,3
+368,368,2
+369,369,5
+370,370,6
+371,371,6
+372,372,2
+373,373,3
+374,374,6
+375,375,2
+376,376,2
+377,377,3
+378,378,2
+379,379,1
+380,380,2
+381,381,1
+382,382,2
+383,383,3
+384,384,4
+385,385,2
+386,386,6
+387,387,4
+388,388,6
+389,389,3
+390,390,5
+391,391,1
+392,392,5
+393,393,2
+394,394,5
+395,395,3
+396,396,2
+397,397,2
+398,398,1
+399,399,4
+400,400,2
+401,401,2
+402,402,1
+403,403,5
+404,404,3
+405,405,2
+406,406,5
+407,407,6
+408,408,5
+409,409,5
+410,410,5
+411,411,4
+412,412,4
+413,413,1
+414,414,4
+415,415,4
+416,416,1
+417,417,2
+418,418,4
+419,419,5
+420,420,3
+421,421,6
+422,422,3
+423,423,1
+424,424,5
+425,425,1
+426,426,1
+427,427,5
+428,428,3
+429,429,4
+430,430,6
+431,431,3
+432,432,6
+433,433,6
+434,434,1
+435,435,2
+436,436,1
+437,437,1
+438,438,5
+439,439,1
+440,440,4
+441,441,1
+442,442,1
+443,443,2
+444,444,4
+445,445,2
+446,446,4
+447,447,6
+448,448,2
+449,449,2
+450,450,6
+451,451,2
+452,452,5
+453,453,5
+454,454,2
+455,455,1
+456,456,5
+457,457,5
+458,458,4
+459,459,1
+460,460,3
+461,461,1
+462,462,5
+463,463,5
+464,464,5
+465,465,1
+466,466,2
+467,467,2
+468,468,5
+469,469,5
+470,470,3
+471,471,6
+472,472,2
+473,473,6
+474,474,5
+475,475,3
+476,476,4
+477,477,6
+478,478,5
+479,479,5
+480,480,1
+481,481,3
+482,482,3
+483,483,2
+484,484,2
+485,485,1
+486,486,6
+487,487,1
+488,488,2
+489,489,3
+490,490,3
+491,491,2
+492,492,5
+493,493,3
+494,494,3
+495,495,1
+496,496,1
+497,497,5
+498,498,4
+499,499,3
+500,500,5
+501,501,3
+502,502,3
+503,503,1
+504,504,4
+505,505,5
+506,506,5
+507,507,5
+508,508,2
+509,509,6
+510,510,3
+511,511,1
+512,512,3
+513,513,4
+514,514,3
+515,515,3
+516,516,1
+517,517,4
+518,518,6
+519,519,1
+520,520,6
+521,521,3
+522,522,4
+523,523,1
+524,524,4
+525,525,4
+526,526,3
+527,527,1
+528,528,1
+529,529,1
+530,530,2
+531,531,2
+532,532,2
+533,533,1
+534,534,5
+535,535,5
+536,536,4
+537,537,5
+538,538,5
+539,539,2
+540,540,5
+541,541,2
+542,542,2
+543,543,3
+544,544,2
+545,545,3
+546,546,2
+547,547,4
+548,548,2
+549,549,1
+550,550,3
+551,551,1
+552,552,2
+553,553,6
+554,554,3
+555,555,5
+556,556,3
+557,557,2
+558,558,4
+559,559,2
+560,560,1
+561,561,4
+562,562,2
+563,563,2
+564,564,1
+565,565,2
+566,566,5
+567,567,2
+568,568,3
+569,569,6
+570,570,4
+571,571,2
+572,572,4
+573,573,5
+574,574,3
+575,575,4
+576,576,4
+577,577,5
+578,578,4
+579,579,5
+580,580,2
+581,581,1
+582,582,6
+583,583,3
+584,584,3
+585,585,6
+586,586,2
+587,587,1
+588,588,3
+589,589,1
+590,590,3
+591,591,5
+592,592,4
+593,593,5
+594,594,2
+595,595,3
+596,596,2
+597,597,6
+598,598,6
+599,599,2
+600,600,3
+601,601,1
+602,602,1
+603,603,4
+604,604,5
+605,605,3
+606,606,6
+607,607,3
+608,608,4
+609,609,5
+610,610,4
+611,611,6
+612,612,1
+613,613,2
+614,614,4
+615,615,5
+616,616,5
+617,617,3
+618,618,2
+619,619,2
+620,620,5
+621,621,6
+622,622,2
+623,623,1
+624,624,6
+625,625,3
+626,626,1
+627,627,2
+628,628,4
+629,629,2
+630,630,4
+631,631,1
+632,632,5
+633,633,4
+634,634,5
+635,635,6
+636,636,5
+637,637,1
+638,638,3
+639,639,6
+640,640,6
+641,641,5
+642,642,4
+643,643,2
+644,644,1
+645,645,3
+646,646,1
+647,647,6
+648,648,2
+649,649,4
+650,650,2
+651,651,4
+652,652,6
+653,653,1
+654,654,3
+655,655,1
+656,656,5
+657,657,4
+658,658,2
+659,659,6
+660,660,4
+661,661,3
+662,662,6
+663,663,2
+664,664,5
+665,665,2
+666,666,4
+667,667,1
+668,668,4
+669,669,6
+670,670,6
+671,671,5
+672,672,3
+673,673,4
+674,674,5
+675,675,2
+676,676,6
+677,677,6
+678,678,1
+679,679,6
+680,680,5
+681,681,4
+682,682,4
+683,683,3
+684,684,2
+685,685,3
+686,686,2
+687,687,2
+688,688,3
+689,689,1
+690,690,6
+691,691,2
+692,692,5
+693,693,6
+694,694,3
+695,695,6
+696,696,5
+697,697,2
+698,698,4
+699,699,3
+700,700,4
+701,701,4
+702,702,2
+703,703,2
+704,704,4
+705,705,3
+706,706,2
+707,707,2
+708,708,6
+709,709,6
+710,710,6
+711,711,4
+712,712,2
+713,713,6
+714,714,1
+715,715,2
+716,716,1
+717,717,2
+718,718,5
+719,719,4
+720,720,2
+721,721,1
+722,722,4
+723,723,1
+724,724,6
+725,725,4
+726,726,6
+727,727,3
+728,728,3
+729,729,3
+730,730,1
+731,731,2
+732,732,3
+733,733,1
+734,734,3
+735,735,6
+736,736,6
+737,737,4
+738,738,3
+739,739,1
+740,740,4
+741,741,3
+742,742,6
+743,743,3
+744,744,3
+745,745,6
+746,746,5
+747,747,4
+748,748,3
+749,749,2
+750,750,2
+751,751,1
+752,752,2
+753,753,1
+754,754,4
+755,755,6
+756,756,6
+757,757,3
+758,758,1
+759,759,1
+760,760,6
+761,761,4
+762,762,5
+763,763,1
+764,764,3
+765,765,4
+766,766,4
+767,767,5
+768,768,6
+769,769,5
+770,770,4
+771,771,5
+772,772,1
+773,773,4
+774,774,5
+775,775,2
+776,776,3
+777,777,6
+778,778,3
+779,779,1
+780,780,6
+781,781,5
+782,782,2
+783,783,2
+784,784,1
+785,785,6
+786,786,3
+787,787,6
+788,788,4
+789,789,2
+790,790,1
+791,791,1
+792,792,5
+793,793,2
+794,794,2
+795,795,4
+796,796,2
+797,797,4
+798,798,3
+799,799,1
+800,800,6
+801,801,2
+802,802,2
+803,803,3
+804,804,1
+805,805,5
+806,806,5
+807,807,1
+808,808,1
+809,809,6
+810,810,4
+811,811,3
+812,812,3
+813,813,2
+814,814,4
+815,815,1
+816,816,5
+817,817,5
+818,818,2
+819,819,3
+820,820,1
+821,821,6
+822,822,6
+823,823,5
+824,824,4
+825,825,2
+826,826,2
+827,827,3
+828,828,5
+829,829,3
+830,830,6
+831,831,1
+832,832,5
+833,833,1
+834,834,1
+835,835,3
+836,836,1
+837,837,4
+838,838,6
+839,839,4
+840,840,2
+841,841,1
+842,842,3
+843,843,1
+844,844,1
+845,845,1
+846,846,6
+847,847,4
+848,848,3
+849,849,5
+850,850,6
+851,851,1
+852,852,3
+853,853,2
+854,854,1
+855,855,6
+856,856,4
+857,857,6
+858,858,5
+859,859,1
+860,860,3
+861,861,1
+862,862,5
+863,863,6
+864,864,3
+865,865,6
+866,866,2
+867,867,6
+868,868,1
+869,869,3
+870,870,3
+871,871,1
+872,872,1
+873,873,6
+874,874,1
+875,875,3
+876,876,4
+877,877,3
+878,878,1
+879,879,4
+880,880,1
+881,881,4
+882,882,2
+883,883,1
+884,884,4
+885,885,3
+886,886,5
+887,887,1
+888,888,3
+889,889,3
+890,890,2
+891,891,6
+892,892,1
+893,893,6
+894,894,4
+895,895,5
+896,896,6
+897,897,5
+898,898,3
+899,899,3
+900,900,3
+901,901,6
+902,902,3
+903,903,4
+904,904,5
+905,905,5
+906,906,1
+907,907,1
+908,908,5
+909,909,6
+910,910,4
+911,911,1
+912,912,5
+913,913,4
+914,914,3
+915,915,3
+916,916,3
+917,917,6
+918,918,4
+919,919,4
+920,920,5
+921,921,6
+922,922,1
+923,923,3
+924,924,6
+925,925,4
+926,926,6
+927,927,1
+928,928,5
+929,929,4
+930,930,4
+931,931,5
+932,932,5
+933,933,2
+934,934,3
+935,935,2
+936,936,1
+937,937,5
+938,938,5
+939,939,2
+940,940,1
+941,941,5
+942,942,4
+943,943,5
+944,944,3
+945,945,6
+946,946,4
+947,947,2
+948,948,4
+949,949,1
+950,950,5
+951,951,2
+952,952,2
+953,953,4
+954,954,3
+955,955,6
+956,956,4
+957,957,2
+958,958,2
+959,959,3
+960,960,3
+961,961,4
+962,962,5
+963,963,3
+964,964,4
+965,965,5
+966,966,2
+967,967,3
+968,968,4
+969,969,1
+970,970,3
+971,971,4
+972,972,1
+973,973,1
+974,974,5
+975,975,4
+976,976,1
+977,977,6
+978,978,2
+979,979,2
+980,980,2
+981,981,5
+982,982,3
+983,983,4
+984,984,4
+985,985,4
+986,986,4
+987,987,1
+988,988,2
+989,989,6
+990,990,4
+991,991,6
+992,992,4
+993,993,1
+994,994,3
+995,995,1
+996,996,4
+997,997,4
+998,998,3
+999,999,6