Commit 71c31986 authored by Waad ALMASRI's avatar Waad ALMASRI
parents 3485b5d6 b6dc6f87
...@@ -40,7 +40,7 @@ ...@@ -40,7 +40,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 1, "execution_count": 79,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
...@@ -287,7 +287,7 @@ ...@@ -287,7 +287,7 @@
"22 1/12/86 6 58 200 1" "22 1/12/86 6 58 200 1"
] ]
}, },
"execution_count": 1, "execution_count": 79,
"metadata": {}, "metadata": {},
"output_type": "execute_result" "output_type": "execute_result"
} }
...@@ -317,12 +317,13 @@ ...@@ -317,12 +317,13 @@
"Les vols où aucun incident n'est relevé n'apportant aucun information\n", "Les vols où aucun incident n'est relevé n'apportant aucun information\n",
"sur l'influence de la température ou de la pression sur les\n", "sur l'influence de la température ou de la pression sur les\n",
"dysfonctionnements, nous nous concentrons sur les expériences où au\n", "dysfonctionnements, nous nous concentrons sur les expériences où au\n",
"moins un joint a été défectueux.\n" "moins un joint a été défectueux.\n",
"\n"
] ]
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 2, "execution_count": 80,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
...@@ -425,7 +426,7 @@ ...@@ -425,7 +426,7 @@
"22 1/12/86 6 58 200 1" "22 1/12/86 6 58 200 1"
] ]
}, },
"execution_count": 2, "execution_count": 80,
"metadata": {}, "metadata": {},
"output_type": "execute_result" "output_type": "execute_result"
} }
...@@ -448,12 +449,12 @@ ...@@ -448,12 +449,12 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 3, "execution_count": 81,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAFaNJREFUeJzt3X2QZXV95/H3p2cGGASFwGZiMSAQWFdKCWALGtxkiMRCqxzWwgfYSjRGnWwJlTImRuK6hLCmaiUxJlaIOroaYUuRh1Vnd3ERNK3REmHUCY/BzCJCgwHFUWkY5oH+7h/3zvFOd0/37aHPvUz3+1XVNfec+zvnfvvL4X76PNxzU1VIkgQwMuwCJElPH4aCJKlhKEiSGoaCJKlhKEiSGoaCJKnRWigk+XiSh5Pcvofnk+SDSTYnuTXJKW3VIknqT5t7Cn8PnDXL868Aju/+rAM+1GItkqQ+tBYKVfVV4MezDDkbuLw6bgIOSfLstuqRJM1t+RBf+wjg/p7p8e68H0wdmGQdnb0JVq5c+cIjjzxyIAU+VZOTk4yMeNqmlz2Zzp5MZ09m9lT68t3vfvdHVfVv5ho3zFDIDPNmvOdGVa0H1gOMjo7Wxo0b26xrwYyNjbFmzZphl/G0Yk+msyfT2ZOZPZW+JPl+P+OGGcXjQO+f/KuBB4dUiySJ4YbCBuAN3auQXgz8tKqmHTqSJA1Oa4ePknwaWAMcnmQc+FNgBUBVfRi4DnglsBl4HHhTW7VIkvrTWihU1XlzPF/A+W29viRp/jy9L0lqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqtBoKSc5KcneSzUkunOH5o5L8Q5LvJLk1ySvbrEeSNLvWQiHJMuAy4BXACcB5SU6YMuw9wFVVdTJwLvB3bdUjSZpbm3sKpwKbq+qeqtoOXAmcPWVMAc/sPn4W8GCL9UiS5pCqamfFyWuAs6rqLd3p3wZOq6oLesY8G/gicCjwDODMqvrWDOtaB6wDWLVq1QuvvPLKVmpeaBMTExx00EHDLuNpxZ5MZ0+msyczeyp9OeOMM75VVaNzjVu+V2vvT2aYNzWBzgP+vqren+QlwBVJnl9Vk7stVLUeWA8wOjpaa9asaaPeBTc2Nsa+Uuug2JPp7Ml09mRmg+hLm4ePxoEje6ZXM/3w0JuBqwCq6hvAAcDhLdYkSZpFm6FwC3B8kmOS7EfnRPKGKWPuA14GkOR5dELhhy3WJEmaRWuhUFU7gQuA64G76FxldEeSS5Ks7Q77Q+CtSf4J+DTwO9XWSQ5J0pzaPKdAVV0HXDdl3kU9j+8ETm+zBklS//xEsySpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqGgiSpYShIkhqthkKSs5LcnWRzkgv3MOZ1Se5MckeST7VZjyRpdsv7GZTk+VV1+3xWnGQZcBnwm8A4cEuSDVV1Z8+Y44E/AU6vqi1JfnE+ryFJWlj97il8OMnNSd6W5JA+lzkV2FxV91TVduBK4OwpY94KXFZVWwCq6uE+1y1JakFfewpV9dLuX/W/C2xMcjPwiaq6YZbFjgDu75keB06bMubfAiT5OrAMuLiq/u/UFSVZB6wDWLVqFWNjY/2UPXQTExP7TK2DYk+msyfT2ZOZDaIvfYUCQFX9S5L3ABuBDwInJwnw7qr6nzMskplWM8PrHw+sAVYD/9g9VPWTKa+9HlgPMDo6WmvWrOm37KEaGxtjX6l1UOzJdPZkOnsys0H0pa/DR0lOTPIB4C7gN4BXVdXzuo8/sIfFxoEje6ZXAw/OMObzVbWjqr4H3E0nJCRJQ9DvOYW/Bb4N/EpVnV9V3waoqgeB9+xhmVuA45Mck2Q/4Fxgw5QxnwPOAEhyOJ3DSffM71eQJC2Ufg8fvRLYWlVPAiQZAQ6oqser6oqZFqiqnUkuAK6nc77g41V1R5JLgI1VtaH73MuT3Ak8Cbyzqh55ir+TJGkv9RsKNwJnAhPd6QOBLwK/OttCVXUdcN2UeRf1PC7gHd0fSdKQ9Xv46ICq2hUIdB8f2E5JkqRh6TcUHktyyq6JJC8EtrZTkiRpWPo9fPR24Ooku64eejbw+nZKkiQNS78fXrslyb8Dnkvn8wf/XFU7Wq1MkjRwfX94DXgRcHR3mZOTUFWXt1KVJGko+r0h3hXALwOb6Fw6Cp1PJxsKkrSI9LunMAqc0L2EVJK0SPV79dHtwC+1WYgkafj63VM4HLize3fUbbtmVtXaVqqSJA1Fv6FwcZtFSJKeHvq9JPUrSZ4DHF9VNyY5kM79jCRJi0i/t85+K3AN8JHurCPo3OFUkrSI9Hui+XzgdOBn0PnCHcDvU5akRabfUNjW/Z5lAJIsZ/q3qEmS9nH9hsJXkrwbWJnkN4Grgf/VXlmSpGHoNxQuBH4I3Ab8Hp3vSNjTN65JkvZR/V59NAl8tPsjSVqk+r330feY4RxCVR274BVJkoZmPvc+2uUA4LXALyx8OZKkYerrnEJVPdLz80BV/TXwGy3XJkkasH4PH53SMzlCZ8/h4FYqkiQNTb+Hj97f83gncC/wugWvRpI0VP1efXRG24VIkoav38NH75jt+ar6q4UpR5I0TPO5+uhFwIbu9KuArwL3t1GUJGk45vMlO6dU1aMASS4Grq6qt7RVmCRp8Pq9zcVRwPae6e3A0QtejSRpqPrdU7gCuDnJZ+l8svnVwOWtVSVJGop+rz768yRfAP59d9abquo77ZUlSRqGfg8fARwI/Kyq/gYYT3JMSzVJkoak36/j/FPgXcCfdGetAP5HW0VJkoaj3z2FVwNrgccAqupBvM2FJC06/YbC9qoqurfPTvKM9kqSJA1Lv6FwVZKPAIckeStwI37hjiQtOv1effSX3e9m/hnwXOCiqrqh1cokSQM3555CkmVJbqyqG6rqnVX1R/0GQpKzktydZHOSC2cZ95oklWR0T2MkSe2bMxSq6kng8STPms+KkywDLgNeAZwAnJfkhBnGHQz8PvDN+axfkrTw+v1E8xPAbUluoHsFEkBV/f4sy5wKbK6qewCSXAmcDdw5Zdx/BS4F/qjfoiVJ7eg3FP5P92c+jmD3u6iOA6f1DkhyMnBkVf3vJHsMhSTrgHUAq1atYmxsbJ6lDMfExMQ+U+ug2JPp7Ml09mRmg+jLrKGQ5Kiquq+qPrkX684M86pn3SPAB4DfmWtFVbUeWA8wOjpaa9as2YtyBm9sbIx9pdZBsSfT2ZPp7MnMBtGXuc4pfG7XgyTXznPd48CRPdOrgQd7pg8Gng+MJbkXeDGwwZPNkjQ8c4VC71/7x85z3bcAxyc5Jsl+wLn8/Et6qKqfVtXhVXV0VR0N3ASsraqN83wdSdICmSsUag+P51RVO4ELgOuBu4CrquqOJJckWTu/MiVJgzDXieZfSfIzOnsMK7uP6U5XVT1ztoWr6jrguinzLtrD2DV9VSxJas2soVBVywZViCRp+ObzfQqSpEXOUJAkNQwFSVLDUJAkNZZMKDwysY1/uv8nPDKxbdilSNK8PTKxja07nmz9PWxJhMLnNz3A6e/7Mr/1sW9y+vu+zIZNDwy7JEnq2673sO/98LHW38MWfSg8MrGNd117K0/smOTRbTt5Ysckf3ztre4xSNon9L6HPVnV+nvYog+F8S1bWTGy+6+5YmSE8S1bh1SRJPVv0O9hiz4UVh+6kh2Tk7vN2zE5yepDVw6pIknq36DfwxZ9KBx20P5ces6JHLBihIP3X84BK0a49JwTOeyg/YddmiTNqfc9bFnS+ntYv1+ys09be9IRnH7c4Yxv2crqQ1caCJL2Kbvew27+xtf4+tqXtvoetiRCATppaxhI2lcddtD+rFyxrPX3sUV/+EiS1D9DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSQ1DQZLUMBQkSY1WQyHJWUnuTrI5yYUzPP+OJHcmuTXJl5I8p816JEmzay0UkiwDLgNeAZwAnJfkhCnDvgOMVtWJwDXApW3VI0maW5t7CqcCm6vqnqraDlwJnN07oKr+oaoe707eBKxusR5J0hyWt7juI4D7e6bHgdNmGf9m4AszPZFkHbAOYNWqVYyNjS1Qie2amJjYZ2odFHsynT2Zzp7MbBB9aTMUMsO8mnFg8lvAKPDrMz1fVeuB9QCjo6O1Zs2aBSqxXWNjY+wrtQ6KPZnOnkxnT2Y2iL60GQrjwJE906uBB6cOSnIm8J+BX6+qbS3WI0maQ5vnFG4Bjk9yTJL9gHOBDb0DkpwMfARYW1UPt1iLJKkPrYVCVe0ELgCuB+4CrqqqO5JckmRtd9hfAAcBVyfZlGTDHlYnSRqANg8fUVXXAddNmXdRz+Mz23z9peSRiW2Mb9nK6kNXcthB+7e+3GJmT4Zr80OPsuXxHWx+6FGOW3XwsMtZcloNBQ3G5zc9wLuuvZUVIyPsmJzk0nNOZO1JR7S23GJmT4bros/dxuU33ccfvmAnf/CBr/KGlxzFJWe/YNhlLSne5mIf98jENt517a08sWOSR7ft5Ikdk/zxtbfyyMTs5+z3drnFzJ4M1+aHHuXym+7bbd7l37iPzQ89OqSKliZDYR83vmUrK0Z2/8+4YmSE8S1bW1luMbMnw7Xp/p/Ma77aYSjs41YfupIdk5O7zdsxOcnqQ1e2stxiZk+G66QjD5nXfLXDUNjHHXbQ/lx6zokcsGKEg/dfzgErRrj0nBPnPEG6t8stZvZkuI5bdTBveMlRu817w0uO8mTzgHmieRFYe9IRnH7c4fO+YmZvl1vM7MlwXXL2C3jDi4/mtm/dxI1/8GIDYQgMhUXisIP236s3sL1dbjGzJ8N13KqDGT9whYEwJB4+kiQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUsNQkCQ1DAVJUqPVUEhyVpK7k2xOcuEMz++f5DPd57+Z5Og265Ekza61UEiyDLgMeAVwAnBekhOmDHszsKWqjgM+ALyvrXokSXNrc0/hVGBzVd1TVduBK4Gzp4w5G/hk9/E1wMuSpMWaJEmzWN7iuo8A7u+ZHgdO29OYqtqZ5KfAYcCPegclWQes605OJLm7lYoX3uFM+V1kT2ZgT6azJzN7Kn15Tj+D2gyFmf7ir70YQ1WtB9YvRFGDlGRjVY0Ou46nE3synT2Zzp7MbBB9afPw0ThwZM/0auDBPY1Jshx4FvDjFmuSJM2izVC4BTg+yTFJ9gPOBTZMGbMBeGP38WuAL1fVtD0FSdJgtHb4qHuO4ALgemAZ8PGquiPJJcDGqtoA/HfgiiSb6ewhnNtWPUOyzx3yGgB7Mp09mc6ezKz1vsQ/zCVJu/iJZklSw1CQJDUMhQWS5N4ktyXZlGRjd97FSR7oztuU5JXDrnPQkhyS5Jok/5zkriQvSfILSW5I8i/dfw8ddp2DtIeeLNltJclze37vTUl+luTtS3k7maUnrW8nnlNYIEnuBUar6kc98y4GJqrqL4dV17Al+STwj1X1se5VaAcC7wZ+XFX/rXtPrEOr6l1DLXSA9tCTt7PEtxVobo/zAJ0Pup7PEt5OdpnSkzfR8nbinoJak+SZwK/RucqMqtpeVT9h99ubfBL4D8OpcPBm6Yk6Xgb8v6r6Pkt4O5mityetMxQWTgFfTPKt7m05drkgya1JPr6Udn+7jgV+CHwiyXeSfCzJM4BVVfUDgO6/vzjMIgdsTz2Bpb2t7HIu8Onu46W8nfTq7Qm0vJ0YCgvn9Ko6hc5dYc9P8mvAh4BfBk4CfgC8f4j1DcNy4BTgQ1V1MvAYMO0W6kvMnnqy1LcVuofS1gJXD7uWp4sZetL6dmIoLJCqerD778PAZ4FTq+qhqnqyqiaBj9K5c+xSMg6MV9U3u9PX0HlDfCjJswG6/z48pPqGYcaeuK0AnT+ovl1VD3Wnl/J2sstuPRnEdmIoLIAkz0hy8K7HwMuB23dt0F2vBm4fRn3DUlX/Ctyf5LndWS8D7mT325u8Efj8EMobij31ZKlvK13nsfthkiW7nfTYrSeD2E68+mgBJDmWzt4BdA4PfKqq/jzJFXR28wq4F/i9XcdIl4okJwEfA/YD7qFz9cQIcBVwFHAf8NqqWjI3QtxDTz7IEt5WkhxI5zb6x1bVT7vzDmNpbycz9aT19xRDQZLU8PCRJKlhKEiSGoaCJKlhKEiSGoaCJKnR2jevSYPWvYTxS93JXwKepHNLCeh8mHD7UAqbRZLfBa7rfn5BGjovSdWi9HS6Q22SZVX15B6e+xpwQVVtmsf6llfVzgUrUOrh4SMtCUnemOTm7j3o/y7JSJLlSX6S5C+SfDvJ9UlOS/KVJPfsuld9krck+Wz3+buTvKfP9b43yc3AqUn+LMktSW5P8uF0vJ7OB5E+011+vyTjSQ7prvvFSW7sPn5vko8kuYHOzfSWJ/mr7mvfmuQtg++qFiNDQYtekufTuSXAr1bVSXQOm57bffpZwBe7NzPcDlxM59YTrwUu6VnNqd1lTgH+Y5KT+ljvt6vq1Kr6BvA3VfUi4AXd586qqs8Am4DXV9VJfRzeOhl4VVX9NrAOeLiqTgVeROcmjEftTX+kXp5T0FJwJp03zo1JAFbSuX0AwNaquqH7+Dbgp1W1M8ltwNE967i+qrYAJPkc8FI6///sab3b+fmtTwBeluSdwAHA4cC3gC/M8/f4fFU90X38cuB5SXpD6Hg6t4OQ9pqhoKUgwMer6r/sNjNZTufNe5dJYFvP497/P6aefKs51ru1uifsuvew+Vs6d0N9IMl76YTDTHby8z34qWMem/I7va2qvoS0gDx8pKXgRuB1SQ6HzlVKe3Go5eXpfLfygXS+Eezr81jvSjoh86Pu3XTP6XnuUeDgnul7gRd2H/eOm+p64G3dANr1nb4r5/k7SdO4p6BFr6puS/JnwI1JRoAdwH8CHpzHar4GfIrOF5xcsetqoX7WW1WPpPO9zLcD3we+2fP0J4CPJdlK57zFxcBHk/wrcPMs9XyEzt1DN3UPXT1MJ6ykp8RLUqU5dK/seX5VvX3YtUht8/CRJKnhnoIkqeGegiSpYShIkhqGgiSpYShIkhqGgiSp8f8B+Q9eu+sB8EwAAAAASUVORK5CYII=\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFYRJREFUeJzt3XuQpXV95/H3Zy7AIBMhsJm4MxBBCFlKAXG4GEx2IokLbgmxiBHcDS5ZMqGE3TK7m8BariHGVEWM2WiJjiOLCqmERFEgu+MiJNUaExCQTIaLgcwiQjMGBFFoHObW3/3jnHlyprun5/TQzzlM9/tV1TXnufa3vz6cj8/l/E6qCkmSABYMuwBJ0kuHoSBJahgKkqSGoSBJahgKkqSGoSBJarQWCkmuSfJkkvt2szxJPppkY5INSU5qqxZJUn/aPFP4DHDmNMvPAo7p/qwGPtFiLZKkPrQWClX1VeB706xyDnBtddwBHJzkFW3VI0nas0VD/N3Lgcd6pke7874zccUkq+mcTbBkyZLXHX744QMp8MUaHx9nwQJv2/SyJ5PZk6nZl8leTE8eeuihp6rqX+xpvWGGQqaYN+WYG1W1FlgLsHLlyrr77rvbrGvWjIyMsGrVqmGX8ZJiTyazJ1OzL5O9mJ4k+XY/6w0zhkeB3v/LvwLYNKRaJEkMNxRuBi7oPoV0GvCDqpp06UiSNDitXT5K8qfAKuCwJKPAbwOLAapqDbAOeDOwEfghcGFbtUiS+tNaKFTV+XtYXsAlbf1+SdLMeWtfktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktQwFCRJDUNBktRoNRSSnJnkwSQbk1w+xfKXJ/mLJH+f5P4kF7ZZjyRpeq2FQpKFwFXAWcBxwPlJjpuw2iXAA1V1ArAK+HCS/dqqSZI0vTbPFE4BNlbVw1W1FbgeOGfCOgUsTRLgIOB7wPYWa5IkTWNRi/teDjzWMz0KnDphnY8BNwObgKXA26tqfOKOkqwGVgMsW7aMkZGRNuqddWNjY/tMrYNiTyazJ1OzL5MNoidthkKmmFcTpv8NsB54I/Aq4NYkf11Vz+6yUdVaYC3AypUra9WqVbNfbQtGRkbYV2odFHsymT2Zmn2ZbBA9afPy0ShweM/0CjpnBL0uBL5QHRuBbwE/1WJNkqRptBkKdwHHJDmye/P4PDqXino9CpwBkGQZcCzwcIs1SZKm0drlo6ranuRS4BZgIXBNVd2f5OLu8jXA7wKfSXIvnctNl1XVU23VJEmaXpv3FKiqdcC6CfPW9LzeBLypzRokSf3zE82SpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqGAqSpIahIElqtBoKSc5M8mCSjUku3806q5KsT3J/kq+0WY8kaXqL+lkpyaur6r6Z7DjJQuAq4BeAUeCuJDdX1QM96xwMfBw4s6oeTfJjM/kdkqTZ1e+ZwpokdyZ5V/eNvB+nABur6uGq2gpcD5wzYZ13AF+oqkcBqurJPvctSWpBX2cKVfWGJMcAvwrcneRO4NNVdes0my0HHuuZHgVOnbDOTwKLk4wAS4GPVNW1E3eUZDWwGmDZsmWMjIz0U/bQjY2N7TO1Doo9mcyeTM2+TDaInvQVCgBV9Y9J3gvcDXwUeG2SAO+pqi9MsUmm2s0Uv/91wBnAEuD2JHdU1UMTfvdaYC3AypUra9WqVf2WPVQjIyPsK7UOij2ZzJ5Mzb5MNoie9HtP4XjgQuDfArcCb6mqe5L8S+B2YKpQGAUO75leAWyaYp2nqup54PkkXwVOAB5CkjRw/d5T+BhwD3BCVV1SVfcAVNUm4L272eYu4JgkRybZDzgPuHnCOjcBP5NkUZID6Vxe+uZM/whJ0uzo9/LRm4HNVbUDIMkC4ICq+mFVXTfVBlW1PcmlwC3AQuCaqro/ycXd5Wuq6ptJ/i+wARgHrp7pU06SpNnTbyjcBvw8MNadPhD4MvDT021UVeuAdRPmrZkw/SHgQ33WIUlqUb+Xjw6oqp2BQPf1ge2UJEkaln5D4fkkJ+2cSPI6YHM7JUmShqXfy0fvBj6XZOfTQ68A3t5OSZKkYen3w2t3Jfkp4Fg6nz/4h6ra1mplkqSB6/vDa8DJwCu727w2CVN9+liStO/q98Nr1wGvAtYDO7qzCzAUJGkO6fdMYSVwXFVNHKZCkjSH9Pv00X3Aj7dZiCRp+Po9UzgMeKA7OuqWnTOr6uxWqpIkDUW/oXBFm0VIkl4a+n0k9StJfgI4pqpu6w5et7Dd0iRJg9bXPYUkvwZ8Hvhkd9Zy4Ma2ipIkDUe/N5ovAU4HnoXOF+4Afp+yJM0x/YbClu73LAOQZBGTv0VNkrSP6zcUvpLkPcCSJL8AfA74i/bKkiQNQ7+hcDnwXeBe4NfpfEfC7r5xTZK0j+r36aNx4FPdH0nSHNXv2EffYop7CFV11KxXJEkampmMfbTTAcDbgB+d/XIkScPU1z2Fqnq65+fxqvoj4I0t1yZJGrB+Lx+d1DO5gM6Zw9JWKpIkDU2/l48+3PN6O/AI8MuzXo0kaaj6ffro59ouRJI0fP1ePvov0y2vqj+cnXIkScM0k6ePTgZu7k6/Bfgq8FgbRUmShmMmX7JzUlU9B5DkCuBzVXVRW4VJkgav32EujgC29kxvBV4569VIkoaq3zOF64A7k3yRzieb3wpc21pVkqSh6Pfpo99L8iXgZ7qzLqyqv2uvLEnSMPR7+QjgQODZqvoIMJrkyJZqkiQNSb9fx/nbwGXAf+/OWgz8cVtFSZKGo98zhbcCZwPPA1TVJhzmQpLmnH5DYWtVFd3hs5O8rL2SJEnD0m8o/HmSTwIHJ/k14Db8wh1JmnP6ffroD7rfzfwscCzwvqq6tdXKJEkDt8czhSQLk9xWVbdW1W9W1X/rNxCSnJnkwSQbk1w+zXonJ9mR5JdmUrwkaXbtMRSqagfwwyQvn8mOkywErgLOAo4Dzk9y3G7W+yBwy0z2L0maff1+ovkF4N4kt9J9Agmgqv7zNNucAmysqocBklwPnAM8MGG9/wTcQGfAPUnSEPUbCv+n+zMTy9l1FNVR4NTeFZIsp/O46xuZJhSSrAZWAyxbtoyRkZEZljIcY2Nj+0ytg2JPJrMnU7Mvkw2iJ9OGQpIjqurRqvrsXuw7U8yrCdN/BFxWVTuSqVbvblS1FlgLsHLlylq1atVelDN4IyMj7Cu1Doo9mcyeTM2+TDaInuzpnsKNO18kuWGG+x4FDu+ZXgFsmrDOSuD6JI8AvwR8PMkvzvD3SJJmyZ4uH/X+3/ejZrjvu4BjumMkPQ6cB7yjd4WqasZPSvIZ4H9X1Y1IkoZiT6FQu3m9R1W1PcmldJ4qWghcU1X3J7m4u3zNjCqVJLVuT6FwQpJn6ZwxLOm+pjtdVfUj021cVeuAdRPmTRkGVfUf+qpYktSaaUOhqhYOqhBJ0vDN5PsUJElznKEgSWoYCpKkhqEgSWrMq1B4emwLf//Y93l6bMuwS5GkGXl6bAubt+1o/f1r3oTCTesf5/QP/hX//uqvc/oH/4qb1z8+7JIkqS8737++9d3nW3//mheh8PTYFi67YQMvbBvnuS3beWHbOL91wwbPGCS95PW+f+2oav39a16Ewugzm1m8YNc/dfGCBYw+s3lIFUlSfwb9/jUvQmHFIUvYNj6+y7xt4+OsOGTJkCqSpP4M+v1rXoTCoQftz5XnHs8BixewdP9FHLB4AVeeezyHHrT/sEuTpGn1vn8tTFp//+r3S3b2eWefuJzTjz6M0Wc2s+KQJQaCpH3GzvevO2//Gn9z9htaff+aN6EAncQ1DCTtiw49aH+WLF7Y+nvYvLh8JEnqj6EgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkhqEgSWoYCpKkRquhkOTMJA8m2Zjk8imW/7skG7o/f5vkhDbrkSRNr7VQSLIQuAo4CzgOOD/JcRNW+xbwr6vqeOB3gbVt1SNJ2rM2zxROATZW1cNVtRW4Hjind4Wq+tuqeqY7eQewosV6JEl7sKjFfS8HHuuZHgVOnWb9/wh8aaoFSVYDqwGWLVvGyMjILJXYrrGxsX2m1kGxJ5PZk6nZl8kG0ZM2QyFTzKspV0x+jk4ovGGq5VW1lu6lpZUrV9aqVatmqcR2jYyMsK/UOij2ZDJ7MjX7MtkgetJmKIwCh/dMrwA2TVwpyfHA1cBZVfV0i/VIkvagzXsKdwHHJDkyyX7AecDNvSskOQL4AvArVfVQi7VIkvrQ2plCVW1PcilwC7AQuKaq7k9ycXf5GuB9wKHAx5MAbK+qlW3VJEmaXpuXj6iqdcC6CfPW9Ly+CLiozRrmi6fHtjD6zGZWHLKEQw/av/Xt5jJ7Mnwbn3iOZ364jY1PPMfRy5YOu5x5pdVQ0GDctP5xLrthA4sXLGDb+DhXnns8Z5+4vLXt5jJ7Mnzvu/Ferr3jUf7ra7bzG//zq1zw+iN4/zmvGXZZ84bDXOzjnh7bwmU3bOCFbeM8t2U7L2wb57du2MDTY1ta2W4usyfDt/GJ57j2jkd3mXft7Y+y8YnnhlTR/GMo7ONGn9nM4gW7/s+4eMECRp/Z3Mp2c5k9Gb71j31/RvM1+wyFfdyKQ5awbXx8l3nbxsdZcciSVraby+zJ8J14+MEzmq/ZZyjs4w49aH+uPPd4Dli8gKX7L+KAxQu48tzj93iDdG+3m8vsyfAdvWwpF7z+iF3mXfD6I7zZPEDeaJ4Dzj5xOacffdiMn5jZ2+3mMnsyfO8/5zVccNorufcbd3Dbb5xmIAyYoTBHHHrQ/nv1Bra3281l9mT4jl62lNEDFxsIQ+DlI0lSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDUMBUlSw1CQJDVaDYUkZyZ5MMnGJJdPsTxJPtpdviHJSW3WI0maXmuhkGQhcBVwFnAccH6S4yasdhZwTPdnNfCJtuqRJO1Zm2cKpwAbq+rhqtoKXA+cM2Gdc4Brq+MO4OAkr2ixJknSNBa1uO/lwGM906PAqX2ssxz4Tu9KSVbTOZMAGEvy4OyW2prDgKeGXcRLjD2ZzJ5Mzb5M9mJ68hP9rNRmKGSKebUX61BVa4G1s1HUICW5u6pWDruOlxJ7Mpk9mZp9mWwQPWnz8tEocHjP9Apg016sI0kakDZD4S7gmCRHJtkPOA+4ecI6NwMXdJ9COg34QVV9Z+KOJEmD0drlo6ranuRS4BZgIXBNVd2f5OLu8jXAOuDNwEbgh8CFbdUzJPvcJa8BsCeT2ZOp2ZfJWu9JqiZdwpckzVN+olmS1DAUJEkNQ2EWJXkkyb1J1ie5uzvviiSPd+etT/LmYdc5SEkOTvL5JP+Q5JtJXp/kR5PcmuQfu/8eMuw6B2k3PZm3x0mSY3v+7vVJnk3y7vl8nEzTk9aPE+8pzKIkjwArq+qpnnlXAGNV9QfDqmuYknwW+Ouqurr7FNqBwHuA71XV73fHxDqkqi4baqEDtJuevJt5fJzs1B0e53E6H3S9hHl8nOw0oScX0vJx4pmCWpPkR4CfBf4XQFVtrarv0xne5LPd1T4L/OJwKhy8aXqijjOA/1dV32YeHycT9PakdYbC7Crgy0m+0R2aY6dLu6PAXjOfToGBo4DvAp9O8ndJrk7yMmDZzs+jdP/9sWEWOWC76wnM3+Ok13nAn3Zfz+fjpFdvT6Dl48RQmF2nV9VJdEZ/vSTJz9IZ+fVVwIl0xnT68BDrG7RFwEnAJ6rqtcDzwKQh1OeZ3fVkPh8nAHQvpZ0NfG7YtbxUTNGT1o8TQ2EWVdWm7r9PAl8ETqmqJ6pqR1WNA5+iM3rsfDEKjFbV17vTn6fzhvjEztFwu/8+OaT6hmHKnszz42Sns4B7quqJ7vR8Pk522qUngzhODIVZkuRlSZbufA28CbhvwlDgbwXuG0Z9w1BV/wQ8luTY7qwzgAfoDG/yzu68dwI3DaG8odhdT+bzcdLjfHa9TDJvj5Meu/RkEMeJTx/NkiRH0Tk7gM4lgj+pqt9Lch2dU70CHgF+fT6N75TkROBqYD/gYTpPTywA/hw4AngUeFtVfW9oRQ7YbnryUeb3cXIgnWH0j6qqH3TnHcr8Pk6m6knr7yeGgiSp4eUjSVLDUJAkNQwFSVLDUJAkNQwFSVKjtW9ekwat+wjjX3YnfxzYQWdICeh8kHDrUAqbRpJfBdZ1P78gDZ2PpGpOeimNTptkYVXt2M2yrwGXVtX6GexvUVVtn7UCpR5ePtK8kOSdSe7sjkH/8SQLkixK8v0kH0pyT5Jbkpya5CtJHt45Vn2Si5J8sbv8wSTv7XO/H0hyJ3BKkt9JcleS+5KsScfb6XwQ6c+62++XZDTJwd19n5bktu7rDyT5ZJJb6QymtyjJH3Z/94YkFw2+q5qLDAXNeUleTWdIgJ+uqhPpXDY9r7v45cCXuwMZbgWuoDP0xNuA9/fs5pTuNicB70hyYh/7vaeqTqmq24GPVNXJwGu6y86sqj8D1gNvr6oT+7i89VrgLVX1K8Bq4MmqOgU4mc4AjEfsTX+kXt5T0Hzw83TeOO9OArCEzvABAJur6tbu63uBH1TV9iT3Aq/s2cctVfUMQJIbgTfQ+e9nd/vdyj8PewJwRpLfBA4ADgO+AXxphn/HTVX1Qvf1m4B/laQ3hI6hMxyEtNcMBc0HAa6pqv+xy8xkEZ03753GgS09r3v/+5h48632sN/N1b1h1x3D5mN0RkN9PMkH6ITDVLbzz2fwE9d5fsLf9K6q+kukWeTlI80HtwG/nOQw6DyltBeXWt6UzncrH0jnG8H+Zgb7XUInZJ7qjqR7bs+y54ClPdOPAK/rvu5db6JbgHd1A2jnd/oumeHfJE3imYLmvKq6N8nvALclWQBsAy4GNs1gN18D/oTOF5xct/NpoX72W1VPp/O9zPcB3wa+3rP408DVSTbTuW9xBfCpJP8E3DlNPZ+kM3ro+u6lqyfphJX0ovhIqrQH3Sd7Xl1V7x52LVLbvHwkSWp4piBJanimIElqGAqSpIahIElqGAqSpIahIElq/H/IxmFZztFAcQAAAABJRU5ErkJggg==\n",
"text/plain": [ "text/plain": [
"<Figure size 432x288 with 1 Axes>" "<Figure size 432x288 with 1 Axes>"
] ]
...@@ -500,7 +501,7 @@ ...@@ -500,7 +501,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 4, "execution_count": 82,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
...@@ -524,10 +525,10 @@ ...@@ -524,10 +525,10 @@
" <th>Method:</th> <td>IRLS</td> <th> Log-Likelihood: </th> <td> -2.5250</td> \n", " <th>Method:</th> <td>IRLS</td> <th> Log-Likelihood: </th> <td> -2.5250</td> \n",
"</tr>\n", "</tr>\n",
"<tr>\n", "<tr>\n",
" <th>Date:</th> <td>Sat, 13 Apr 2019</td> <th> Deviance: </th> <td> 0.22231</td> \n", " <th>Date:</th> <td>Wed, 26 Aug 2020</td> <th> Deviance: </th> <td> 0.22231</td> \n",
"</tr>\n", "</tr>\n",
"<tr>\n", "<tr>\n",
" <th>Time:</th> <td>19:11:24</td> <th> Pearson chi2: </th> <td> 0.236</td> \n", " <th>Time:</th> <td>16:53:28</td> <th> Pearson chi2: </th> <td> 0.236</td> \n",
"</tr>\n", "</tr>\n",
"<tr>\n", "<tr>\n",
" <th>No. Iterations:</th> <td>4</td> <th> Covariance Type: </th> <td>nonrobust</td>\n", " <th>No. Iterations:</th> <td>4</td> <th> Covariance Type: </th> <td>nonrobust</td>\n",
...@@ -555,8 +556,8 @@ ...@@ -555,8 +556,8 @@
"Model Family: Binomial Df Model: 1\n", "Model Family: Binomial Df Model: 1\n",
"Link Function: logit Scale: 1.0000\n", "Link Function: logit Scale: 1.0000\n",
"Method: IRLS Log-Likelihood: -2.5250\n", "Method: IRLS Log-Likelihood: -2.5250\n",
"Date: Sat, 13 Apr 2019 Deviance: 0.22231\n", "Date: Wed, 26 Aug 2020 Deviance: 0.22231\n",
"Time: 19:11:24 Pearson chi2: 0.236\n", "Time: 16:53:28 Pearson chi2: 0.236\n",
"No. Iterations: 4 Covariance Type: nonrobust\n", "No. Iterations: 4 Covariance Type: nonrobust\n",
"===============================================================================\n", "===============================================================================\n",
" coef std err z P>|z| [0.025 0.975]\n", " coef std err z P>|z| [0.025 0.975]\n",
...@@ -567,7 +568,7 @@ ...@@ -567,7 +568,7 @@
"\"\"\"" "\"\"\""
] ]
}, },
"execution_count": 4, "execution_count": 82,
"metadata": {}, "metadata": {},
"output_type": "execute_result" "output_type": "execute_result"
} }
...@@ -578,6 +579,7 @@ ...@@ -578,6 +579,7 @@
"data[\"Success\"]=data.Count-data.Malfunction\n", "data[\"Success\"]=data.Count-data.Malfunction\n",
"data[\"Intercept\"]=1\n", "data[\"Intercept\"]=1\n",
"\n", "\n",
"data[\"Frequency\"]=data.Malfunction/data.Count\n",
"logmodel=sm.GLM(data['Frequency'], data[['Intercept','Temperature']], family=sm.families.Binomial(sm.families.links.logit)).fit()\n", "logmodel=sm.GLM(data['Frequency'], data[['Intercept','Temperature']], family=sm.families.Binomial(sm.families.links.logit)).fit()\n",
"\n", "\n",
"logmodel.summary()" "logmodel.summary()"
...@@ -605,12 +607,12 @@ ...@@ -605,12 +607,12 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 5, "execution_count": 16,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEKCAYAAAAcgp5RAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAGyFJREFUeJzt3X2UVPWd5/H3pxuQBhEjkhkFDWSWtHF9ABRQWZ3WqGhORLPrE2vGMRNCdmeMk83Knng2E43Rc2YHd2I26zgy6jgxiUo8iiSHCahjT2Y8PoCCILAIY4g2JEGND7Q2Snd/9497u6kuqunqpvqhfnxe5/Tpurd+de/3V7fvp27fuvUrRQRmZpaumsEuwMzM+peD3swscQ56M7PEOejNzBLnoDczS5yD3swscT0GvaR7Je2U9HI390vS/5G0VdI6SdMrX6aZmfVVOUf09wEX7Of+C4Ep+c8C4M4DL8vMzCqlx6CPiF8Av9tPk4uBH0TmWeBwSUdVqkAzMzswwyqwjAnA6wXTTfm8Xxc3lLSA7KifkSNHnnLsscdWYPVDU3t7OzU16b4FknL/Uu4buH/V7pVXXnkzIsb35jGVCHqVmFdyXIWIWAwsBqivr4/NmzdXYPVDU2NjIw0NDYNdRr9JuX8p9w3cv2on6Ve9fUwlXvaagGMKpicCOyqwXDMzq4BKBP0y4Or86pvTgHcjYp/TNmZmNjh6PHUj6QGgAThSUhNwIzAcICL+FlgOfBbYCnwAfLG/ijUzs97rMegjYl4P9wfwZxWryMyqwp49e2hqamL37t2DXUoXY8eOZdOmTYNdxgEbOXIkEydOZPjw4Qe8rEq8GWtmB6GmpibGjBnDpEmTkEpdkzE4du3axZgxYwa7jAMSEbz11ls0NTUxefLkA15eutcgmVm/2r17N+PGjRtSIZ8KSYwbN65i/y056M2szxzy/aeSz62D3swscT5Hb2ZVq7a2lhNPPLFzeunSpYwbN24QKxqaHPRmVrXq6upYu3Ztl3m7du3qvN3a2sqwYY45n7oxs6T86Ec/4rLLLuOiiy7i/PPPB2DRokXMmDGDk046iRtvvLGz7a233kp9fT3nnnsu8+bN47bbbgOgoaGB1atXA/Dmm28yadIkANra2li4cGHnsu666y5g77ALl156KccddxxXXXUV2ZXnsGrVKs444wxOPvlkZs6cya5duzjzzDO7vEDNnj2bdevW9dtz4pc6Mztg3/7pBjbueK+iyzz+6MO48aJ/v982LS0tTJ06FYDJkyfz6KOPAvDMM8+wbt06jjjiCFauXMmWLVt4/vnniQjmzp3LL37xC0aPHs2DDz7ImjVraG1tZfr06Zxyyin7Xd8999zD2LFjWbVqFR9++CGzZ8/ufDFZs2YNGzZs4Oijj2b27Nk8/fTTzJw5kyuuuIKHHnqIGTNm8N5771FXV8f8+fO57777uP3223nllVf48MMPOemkkyrwrJXmoDezqlXq1A3AeeedxxFHHAHAypUrWblyJdOmTQOgubmZLVu2sGvXLj7/+c8zatQoAObOndvj+lauXMm6det4+OGHAXj33XfZsmULI0aMYObMmUycOBGAqVOnsm3bNsaOHctRRx3FjBkzADjssMMAuOyyy/jOd77DokWLuPfee7nmmmsO7InogYPezA5YT0feA2306NGdtyOCG264ga985Std2tx+++3dXsI4bNgw2tvbAbpcyx4RfP/732fOnDld2jc2NnLIIYd0TtfW1tLa2kpElFzHqFGjOO+883jsscdYsmRJ52mi/uJz9GaWtDlz5nDvvffS3NwMwPbt29m5cydnnXUWjz76KC0tLezatYuf/vSnnY+ZNGkSL7zwAkDn0XvHsu6880727NkDwCuvvML777/f7bqPO+44duzYwapVq4DsjeLW1lYA5s+fz3XXXceMGTM6//voLz6iN7OknX/++WzatInTTz8dgEMPPZQf/vCHTJ8+nSuuuIKpU6fyiU98gjPPPLPzMddffz2XX345999/P+ecc07n/Pnz57Nt2zamT59ORDB+/HiWLl3a7bpHjBjBQw89xFe/+lVaWlqoq6vjiSee4NBDD+WUU07hsMMO44tfHIBxICNiUH4+9alPRcqeeuqpwS6hX6Xcv5T7FlG5/m3cuLEiy6m09957r0+Pu/HGG2PRokUVrqZ727dvjylTpkRbW1u3bUo9x8Dq6GXe+tSNmdkA+8EPfsCsWbO49dZbB+RrD33qxswMuOmmmwZsXVdffTVXX331gK3PR/Rm1mcRJb8e2iqgks+tg97M+mTkyJG89dZbDvt+EPl49CNHjqzI8nzqxsz6ZOLEiTQ1NfHGG28Mdild7N69u2IBOZg6vmGqEhz0ZtYnw4cPr8i3H1VaY2Nj56dgLeNTN2ZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeLKCnpJF0jaLGmrpG+UuP9YSU9JWiNpnaTPVr5UMzPrix6DXlItcAdwIXA8ME/S8UXNvgksiYhpwJXA31S6UDMz65tyjuhnAlsj4tWI+Ah4ELi4qE0Ah+W3xwI7KleimZkdCPX0De6SLgUuiIj5+fQfAbMi4tqCNkcBK4GPAaOBcyPihRLLWgAsABg/fvwpS5YsqVQ/hpzm5mYOPfTQwS6j36Tcv5T7Bu5ftTv77LNfiIhTe/OYcr4cXCXmFb86zAPui4j/Lel04H5JJ0REe5cHRSwGFgPU19dHQ0NDb2qtKo2Njbh/1SnlvoH7dzAq59RNE3BMwfRE9j018yVgCUBEPAOMBI6sRIFmZnZgygn6VcAUSZMljSB7s3VZUZvXgM8ASPo0WdC/UclCzcysb3oM+ohoBa4FVgCbyK6u2SDpZklz82b/HfiypJeAB4BroqeT/2ZmNiDKOUdPRCwHlhfN+1bB7Y3A7MqWZmZmleBPxpqZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWuLKCXtIFkjZL2irpG920uVzSRkkbJP24smWamVlfDeupgaRa4A7gPKAJWCVpWURsLGgzBbgBmB0Rb0v6eH8VbGZmvVPOEf1MYGtEvBoRHwEPAhcXtfkycEdEvA0QETsrW6aZmfVVj0f0wATg9YLpJmBWUZtPAUh6GqgFboqInxcvSNICYAHA+PHjaWxs7EPJ1aG5udn9q1Ip9w3cv4NROUGvEvOixHKmAA3AROBfJJ0QEe90eVDEYmAxQH19fTQ0NPS23qrR2NiI+1edUu4buH8Ho3JO3TQBxxRMTwR2lGjzWETsiYhfApvJgt/MzAZZOUG/CpgiabKkEcCVwLKiNkuBswEkHUl2KufVShZqZmZ902PQR0QrcC2wAtgELImIDZJuljQ3b7YCeEvSRuApYGFEvNVfRZuZWfnKOUdPRCwHlhfN+1bB7QC+nv+YmdkQ4k/GmpklzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeIc9GZmiXPQm5klzkFvZpY4B72ZWeLKCnpJF0jaLGmrpG/sp92lkkLSqZUr0czMDkSPQS+pFrgDuBA4Hpgn6fgS7cYA1wHPVbpIMzPru3KO6GcCWyPi1Yj4CHgQuLhEu+8AfwXsrmB9ZmZ2gIaV0WYC8HrBdBMwq7CBpGnAMRHxM0nXd7cgSQuABQDjx4+nsbGx1wVXi+bmZvevSqXcN3D/DkblBL1KzIvOO6Ua4LvANT0tKCIWA4sB6uvro6Ghoawiq1FjYyPuX3VKuW/g/h2Myjl10wQcUzA9EdhRMD0GOAFolLQNOA1Y5jdkzcyGhnKCfhUwRdJkSSOAK4FlHXdGxLsRcWRETIqIScCzwNyIWN0vFZuZWa/0GPQR0QpcC6wANgFLImKDpJslze3vAs3M7MCUc46eiFgOLC+a961u2jYceFlmZlYp/mSsmVniHPRmZolz0JuZJc5Bb2aWOAe9mVniyrrqxqxSlq7ZzqIVm9nxTgtHH17Hwjn1XDJtwmCXZf3A23rocNDbgFm6Zjs3PLKelj1tAGx/p4UbHlkP4ABIjLf10OJTNzZgFq3Y3Lnjd2jZ08aiFZsHqSLrL97WQ4uD3gbMjndaejXfqpe39dDioLcBc/Thdb2ab9XL23pocdDbgFk4p5664bVd5tUNr2XhnPpBqsj6i7f10OI3Y23AdLwJ5ysx0udtPbQ46G1AXTJtgnf2g4S39dDhUzdmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOAe9mVniHPRmZolz0JuZJc5Bb2aWOH9nrJmVJSJoaw/aA9oj8h9oa49u72tvL9Eugvb20svYp13RMtoj8mloi6L15u02vbaH15/9Vba89qAtStS+T61BW3vpPnZMd3df5DV31Fiq9s52HevrspyCdlG0vKK+tkfftp2D3gZN4Y7dJRwiiPaCnaI4GPLbe3f0vfdFqeX1cud8eUcrv3uxqcdgiYKAKL1Tlg6P0kGwd3n7C7G2gvW2t5doV+r5KgqSDz7YzSHPPLnf9RSHS1tfE2awbHy57KYS1ErUSNTUQI1ErYQENTUdt0WNoLZmb7u9j8nuqylaRo0K5teI2hoxvOPxUr4sAFFbky1b+bq7PK5omWv68HSUFfSSLgC+B9QCd0fEXxbd/3VgPtAKvAH8SUT8qg/19Nk+r/qljgAKjhpKvcruEzzFO0vnEUL2uP0dtaz/TSvvvbRjnzCIotvdBl3Rq3rbfu7b35HBPvOLQ6rLkcq+Ry7dHYE1N3/AyFVPdd7X2ZcSwVNYb9fneyD/Qnpp3UsH9PDC8FBhQOThURwEtTXZziyU7/D7Pr62RojC8MkeO6y2hkOG7Q2cjmDqWEbn8vL17vztbzn6qCP3CaWOdvvczutUPr+wPqlrqGX9UNewLFFDTan11AhREKI1RcuVgL3PRZdw7QjfGvHsM88we/YZBbUL1RQ9FwXPrfLlVotb+vCYHoNeUi1wB3Ae0ASskrQsIjYWNFsDnBoRH0j6r8BfAVfsb7nbm9s596//uXQQlArtff716Ro8Q9Lavrz2ltbxh93lFb9oR+g86ug4AqnpLmz2hkoWGnuPHrJ1lA6O2jwUJHizpoWjfv/wLjt58Tpqa9Sl9pJHPwV9KWyX9TN7LHmfawtrrum6jNqa7LkR+4ZNjYAuobdvUBTWvnrV85x+2qx9gqXrkV/R9ig6EhzK4dHY2EhDw8mDXUa/+djIGj4+ZuRglzGklHNEPxPYGhGvAkh6ELgY6Az6iHiqoP2zwBd6WujwGlH/e2M6d8CSQdBNKOzdWfNQoGvoFR/xFP6bpOLbRUc1xUddxQGyN/T2hlvx0UiNxIurVzNr1oyCYN7fkR37HAWpKJiGmiwspg12Gf2iaXQNnxg3erDLMKuYcoJ+AvB6wXQTMGs/7b8E/GOpOyQtABYAjB8/nssmvFdmmX0Q+U8FteU/5RirD2ja+EJlCxhCmpubaWxsHOwy+kXKfQP372BUTtCXOpwsGaGSvgCcCvxhqfsjYjGwGKC+vj4aGhrKq7IKZUe8DYNdRr9JuX8p9w3cv4NROUHfBBxTMD0R2FHcSNK5wP8E/jAiPqxMeWZmdqDK+cDUKmCKpMmSRgBXAssKG0iaBtwFzI2InZUv08zM+qrHoI+IVuBaYAWwCVgSERsk3Sxpbt5sEXAo8BNJayUt62ZxZmY2wMq6jj4ilgPLi+Z9q+D2uRWuy6zXlq7ZzqIVm9nxTgtHH17Hwjn1APvMu2TahAGtoT/X1xvfXLqeB557na+dsIcv3bCcebOO4ZZLThzssmwA+JOxloSla7ZzwyPradmTXRe1/Z0WFv7kJRDsaYvOeTc8sh6gX8K3VA39ub7e+ObS9fzw2dc6p9siOqcd9unzoGaWhEUrNncGbIc97dEZ8h1a9rSxaMXmAauhP9fXGw8893qv5ltaHPSWhB3vtPRL20rU0F/r6422bsab6G6+pcVBb0k4+vC6fmlbiRr6a329UdvNp6u7m29pcdBbEhbOqadueG2XecNrxPDarkFWN7y2803agaihP9fXG/NmHdOr+ZYWvxlrSeh4s3Mwr7rprobBfiMW9r7h2nFOvlbyVTcHEQe9JeOSaRNKhupABm13NQwFt1xyIrdcciKNjY3821UNg12ODSCfujEzS5yD3swscQ56M7PEOejNzBLnoDczS5yD3swscQ56M7PEOejNzBLnoDczS5yD3swscQ56M7PEOejNzBLnoDczS5yD3swscQ56M7PEOejNzBLnoDczS5yD3swscQ56M7PEOejNzBLnoDczS5yD3swscQ56M7PEOejNzBLnoDczS5yD3swscQ56M7PEOejNzBJXVtBLukDSZklbJX2jxP2HSHoov/85SZMqXaiZmfVNj0EvqRa4A7gQOB6YJ+n4omZfAt6OiH8HfBf4X5Uu1MzM+qacI/qZwNaIeDUiPgIeBC4uanMx8A/57YeBz0hS5co0M7O+GlZGmwnA6wXTTcCs7tpERKukd4FxwJuFjSQtABbkkx9KerkvRVeJIynqf2JS7l/KfQP3r9rV9/YB5QR9qSPz6EMbImIxsBhA0uqIOLWM9Vcl9696pdw3cP+qnaTVvX1MOadumoBjCqYnAju6ayNpGDAW+F1vizEzs8orJ+hXAVMkTZY0ArgSWFbUZhnwx/ntS4F/ioh9jujNzGzg9XjqJj/nfi2wAqgF7o2IDZJuBlZHxDLgHuB+SVvJjuSvLGPdiw+g7mrg/lWvlPsG7l+163X/5ANvM7O0+ZOxZmaJc9CbmSVuQIJe0khJz0t6SdIGSd/O50/Oh0zYkg+hMGIg6ukPkmolrZH0s3w6pb5tk7Re0tqOS7skHSHp8bx/j0v62GDX2VeSDpf0sKT/J2mTpNNT6Z+k+ny7dfy8J+lrCfXvv+WZ8rKkB/KsSWnf+/O8bxskfS2f1+ttN1BH9B8C50TEycBU4AJJp5ENlfDdiJgCvE02lEK1+nNgU8F0Sn0DODsiphZcn/wN4Mm8f0/m09Xqe8DPI+I44GSy7ZhE/yJic77dpgKnAB8Aj5JA/yRNAK4DTo2IE8guFrmSRPY9SScAXyYbneBk4HOSptCXbRcRA/oDjAJeJPt07ZvAsHz+6cCKga6nQn2amD/h5wA/I/sAWRJ9y+vfBhxZNG8zcFR++yhg82DX2ce+HQb8kvzChNT6V9Sn84GnU+kfez+RfwTZFYQ/A+aksu8BlwF3F0z/BfA/+rLtBuwcfX5qYy2wE3gc+DfgnYhozZs0kW24anQ72QZoz6fHkU7fIPuU80pJL+TDWAD8XkT8GiD//fFBq+7AfBJ4A/j7/NTb3ZJGk07/Cl0JPJDfrvr+RcR24DbgNeDXwLvAC6Sz770MnCVpnKRRwGfJPpja6203YEEfEW2R/fs4kexfkU+XajZQ9VSKpM8BOyPihcLZJZpWXd8KzI6I6WQjmP6ZpLMGu6AKGgZMB+6MiGnA+1ThaYye5Oep5wI/GexaKiU/N30xMBk4GhhN9jdarCr3vYjYRHYa6nHg58BLQOt+H9SNAb/qJiLeARqB04DD8yEToPTQCtVgNjBX0jaykT3PITvCT6FvAETEjvz3TrLzuzOB30o6CiD/vXPwKjwgTUBTRDyXTz9MFvyp9K/DhcCLEfHbfDqF/p0L/DIi3oiIPcAjwBmkte/dExHTI+Issg+jbqEP226grroZL+nw/HYd2QbaBDxFNmQCZEMoPDYQ9VRSRNwQERMjYhLZv8b/FBFXkUDfACSNljSm4zbZed6X6TrsRdX2LyJ+A7wuqWNEwM8AG0mkfwXmsfe0DaTRv9eA0ySNyodF79h2Sex7AJI+nv8+FviPZNuw19tuQD4ZK+kksvHqa8leXJZExM2SPkl2FHwEsAb4QkR82O8F9RNJDcD1EfG5VPqW9+PRfHIY8OOIuFXSOGAJcCzZDndZRFTlQHaSpgJ3AyOAV4Evkv+dkkb/RpG9afnJiHg3n5fE9ssv1b6C7JTGGmA+2Tn5qt/3ACT9C9l7fnuAr0fEk33Zdh4Cwcwscf5krJlZ4hz0ZmaJc9CbmSXOQW9mljgHvZlZ4sr5cnCzAZVfPvZkPvn7QBvZMAUAMyPio0EpbD8k/QmwPL8u32xI8eWVNqRJuglojojbhkAttRHR1s19/wpcGxFre7G8YQVjspj1G5+6saoi6Y+VfbfBWkl/I6lG0jBJ70haJOlFSSskzZL0z5JelfTZ/LHzJT2a379Z0jfLXO4tkp4HZkr6tqRV+Rjhf6vMFWTDbz+UP36EpKaCT4OfJumJ/PYtku6S9DjZQGrDJP11vu51kuYP/LNqqXPQW9XIx+f+PHBGPkDeMPZ+Ef1YYGU++NpHwE1kH4m/DLi5YDEz88dMB/6zpKllLPfFiJgZEc8A34uIGcCJ+X0XRMRDwFrgisjGfu/p1NI04KKI+CNgAdmgeDOBGWSDxh3bl+fHrDs+R2/V5FyyMFydDW1CHdlH+wFaIuLx/PZ64N2IaJW0HphUsIwVEfE2gKSlwH8g2w+6W+5H7B0CAuAzkhYCI4EjyYbF/cde9uOxiNid3z4f+LSkwheWKWQfbTerCAe9VRMB90bEX3SZmY1UWHgU3U72rWYdtwv/zovflIoeltsS+RtZ+Zgx/xeYHhHbJd1CFviltLL3P+biNu8X9elPI+JJzPqJT91YNXkCuFzSkZBdndOH0xznK/uO2FFkY5k/3Yvl1pG9cLyZj+j5nwru2wWMKZjeRvbVfRS1K7YC+NOOYXWVfcdrXS/7ZLZfPqK3qhER6/PRCp+QVEM2ot9/oXfjjf8r8GPgD4D7O66SKWe5EfGWpH8gG6b5V8BzBXf/PXC3pBay9wFuAv5O0m+A5/dTz11koxCuzU8b7SR7ATKrGF9eaQeN/IqWEyLia4Ndi9lA8qkbM7PE+YjezCxxPqI3M0ucg97MLHEOejOzxDnozcwS56A3M0vc/wcowwoTqhaBUgAAAABJRU5ErkJggg==\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGzdJREFUeJzt3X+UVOWd5/H3t6tBGhohoGGAJoHM4cA6UX41jUriNkYBc+KvWQ2io4k7LHEnJJPdIxs5J7OaWT1n57S7h0zWiIwyTOLR1nEVNcsG1E3HiauxQRAEhh9riDadBDGj0Noo3f3dP+6t6qrqbrq6qO6qevy8zulD3VvPfe7z7aI+dfupW7fM3RERkbBUFHsAIiJSeAp3EZEAKdxFRAKkcBcRCZDCXUQkQAp3EZEA9RvuZrbBzI6a2Rt93G9m9rdmdsjMdpnZ3MIPU0REBiKXI/eNwNLT3H8FMD3+WQncf+bDEhGRM9FvuLv7i8AfTtPkauDHHnkFGGtmEws1QBERGbjKAvQxGXg7bbklXvfb7IZmtpLo6J6qqqp5U6ZMyWuHXV1dVFSE8XaBailNodQSSh2gWpIOHDhwzN3P7a9dIcLdelnX6zUN3H09sB6gtrbWt23bltcOm5qaqK+vz2vbUqNaSlMotYRSB6iWJDP7TS7tCvEy2AKkH4LXAK0F6FdERPJUiHB/BrglPmvmQuB9d+8xJSMiIkOn32kZM3sUqAfOMbMW4E5gGIC7rwM2A18GDgEfArcO1mBFRCQ3/Ya7uy/v534HvlmwEYlIWTh16hQtLS2cPHlySPY3ZswY9u3bNyT7Gmy51DJixAhqamoYNmxYXvsoxBuqIvIJ1NLSwujRo5k6dSpmvZ1XUVgnTpxg9OjRg76fodBfLe7Ou+++S0tLC9OmTctrH2GcVyQiQ+7kyZOMHz9+SIL9k8bMGD9+/Bn9VaRwF5G8KdgHz5n+bhXuIiIB0py7iJStRCLB+eefn1retGkTU6dOLd6ASojCXUTKVlVVFTt37uzz/o6ODiorP5kxp2kZEQnKxo0buf7667nyyitZvHgxAA0NDcyfP58LLriAO++8M9X2nnvuYcaMGVx22WUsX76ce++9F4D6+nqSl0c5duxY6q+Bzs5OVq9enerrgQceALovJ3Ddddcxc+ZMbrrpJqKzxKG5uZmLL76YWbNmUVdXx4kTJ1iyZEnGi9LChQvZtWtXQX8Pn8yXNBEpqO8/u4e9rccL2ud5k87mziv/5LRt2tvbmT17NgDTpk3jqaeeAuDll19m165djBs3jq1bt3Lw4EFeffVV3J2rrrqKF198kVGjRtHY2MiOHTvo6Ohg7ty5zJs377T7e+ihhxgzZgzNzc189NFHLFy4MPUCsmPHDvbs2cOkSZNYuHAhL730EnV1dSxbtozHHnuM+fPnc/z4caqqqrjlllvYuHEja9eu5cCBA3z00UdccMEFBfitdVO4i0jZ6mta5vLLL2fcuHEAbN26la1btzJnzhwA2traOHjwICdOnODaa69l5MiRAFx11VX97m/r1q3s2rWLJ554AoD333+fgwcPMnz4cOrq6qipqQFg9uzZHD58mDFjxjBx4kTmz58PwNlnnw3Atddey8KFC2loaGDDhg18/etfP7NfRC8U7iJyxvo7wh5qo0aNSt12d9asWcM3vvGNjDZr167t83TDyspKurq6ADLONXd3fvjDH7JkyZKM9k1NTZx11lmp5UQiQUdHB+7e6z5GjhzJ5ZdfztNPP83jjz9OvlfIPR3NuYtI0JYsWcKGDRtoa2sD4MiRIxw9epRLLrmEp556ivb2dk6cOMGzzz6b2mbq1Kls374dIHWUnuzr/vvv59SpUwAcOHCADz74oM99z5w5k9bWVpqbm4Hok6kdHR0ArFixgm9/+9vMnz8/9VdGIenIXUSCtnjxYvbt28dFF10EQHV1NQ8//DBz585l2bJlzJ49m89+9rN88YtfTG1z++2389WvfpWf/OQnXHrppan1K1as4PDhw8ydOxd359xzz2XTpk197nv48OE89thjfOtb36K9vZ2qqiqef/55AObNm8fZZ5/NrbcO0rUW3b0oP/PmzfN8/fznP89721KjWkpTKLUMZh179+4dtL57c/z48UHt/8477/SGhoZB3UfS8ePH/ciRIz59+nTv7Ozss11vv2Ngm+eQsZqWEREZYo888ggLFizgnnvuGbSvDtS0jIgIcNdddw3Zvm688cYeb/AWmo7cRSRv7r1+XbIUwJn+bhXuIpKXESNG8O677yrgB4HH13MfMWJE3n1oWkZE8lJTU0NLSwvvvPPOkOzv5MmTZxR2pSSXWpLfxJQvhbuI5GXYsGF5f0tQPpqamlKfMi13Q1GLpmVERAKkcBcRCZDCXUQkQAp3EZEAKdxFRAKkcBcRCZDCXUQkQAp3EZEAKdxFRAKkcBcRCZDCXUQkQAp3EZEAKdxFRAKkcBcRCZDCXUQkQAp3EZEA5RTuZrbUzPab2SEzu6OX+8eY2bNm9rqZ7TGzWws/VBERyVW/4W5mCeA+4ArgPGC5mZ2X1eybwF53nwXUA//NzIYXeKwiIpKjXI7c64BD7v6mu38MNAJXZ7VxYLSZGVAN/AHoKOhIRUQkZ9bfN5eb2XXAUndfES/fDCxw91VpbUYDzwAzgdHAMnf/X730tRJYCTBhwoR5jY2NeQ26ra2N6urqvLYtNaqlNIVSSyh1gGpJWrRo0XZ3r+2vXS5fkG29rMt+RVgC7AQuBf4YeM7M/sndj2ds5L4eWA9QW1vr9fX1Oey+p6amJvLdttSoltIUSi2h1AGqZaBymZZpAaakLdcArVltbgWe9Mgh4NdER/EiIlIEuYR7MzDdzKbFb5LeQDQFk+4t4EsAZjYBmAG8WciBiohI7vqdlnH3DjNbBWwBEsAGd99jZrfF968D/guw0cx2E03jfNfdjw3iuEVE5DRymXPH3TcDm7PWrUu73QosLuzQREQkX/qEqohIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIByinczWypme03s0NmdkcfberNbKeZ7TGzXxR2mCIiMhCV/TUwswRwH3A50AI0m9kz7r43rc1Y4EfAUnd/y8w+PVgDFhGR/uVy5F4HHHL3N939Y6ARuDqrzY3Ak+7+FoC7Hy3sMEVEZCDM3U/fwOw6oiPyFfHyzcACd1+V1mYtMAz4E2A08AN3/3Evfa0EVgJMmDBhXmNjY16Dbmtro7q6Oq9tS41qKU2h1BJKHaBakhYtWrTd3Wv7a9fvtAxgvazLfkWoBOYBXwKqgJfN7BV3P5Cxkft6YD1AbW2t19fX57D7npqamsh321KjWkpTKLWEUgeoloHKJdxbgClpyzVAay9tjrn7B8AHZvYiMAs4gIiIDLlc5tybgelmNs3MhgM3AM9ktXka+KKZVZrZSGABsK+wQxURkVz1e+Tu7h1mtgrYAiSADe6+x8xui+9f5+77zOxnwC6gC3jQ3d8YzIGLiEjfcpmWwd03A5uz1q3LWm4AGgo3NBERyZc+oSoiEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIByinczWypme03s0Nmdsdp2s03s04zu65wQxQRkYHqN9zNLAHcB1wBnAcsN7Pz+mj3N8CWQg9SREQGJpcj9zrgkLu/6e4fA43A1b20+xbwP4GjBRyfiIjkwdz99A2iKZal7r4iXr4ZWODuq9LaTAYeAS4FHgJ+6u5P9NLXSmAlwIQJE+Y1NjbmNei2tjaqq6vz2rbUqJbSFEotodQBqiVp0aJF2929tr92lTn0Zb2sy35FWAt81907zXprHm/kvh5YD1BbW+v19fU57L6npqYm8t221KiW0hRKLaHUAaploHIJ9xZgStpyDdCa1aYWaIyD/Rzgy2bW4e6bCjJKEREZkFzCvRmYbmbTgCPADcCN6Q3cfVrytpltJJqWUbCLiBRJv+Hu7h1mtoroLJgEsMHd95jZbfH96wZ5jCIiMkC5HLnj7puBzVnreg11d//6mQ9LRETOhD6hKiISIIW7iEiAFO4iIgFSuIuIBEjhLiISoJzOlhEZLJt2HKFhy35a32tn0tgqVi+ZwTVzJhd7WJIjPX6lS+EuRbNpxxHWPLmb9lOdABx5r501T+4GUECUAT1+pU3TMlI0DVv2p4Ihqf1UJw1b9hdpRDIQevxKm8Jdiqb1vfYBrZfSosevtCncpWgmja0a0HopLXr8SpvCXYpm9ZIZVA1LZKyrGpZg9ZIZRRqRDIQev9KmN1SlaJJvuulsi/Kkx6+0KdylqK6ZM1lhUMb0+JUuTcuIiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgBTuIiIBUriLiARI4S4iEiCFu4hIgPQdqiKSoavL6XSns8txhy6Plru6onVdyXVdTpdHbTq7utt0xcvpbaKf7vVdXaS173ubzi5Sbfa9dYq3X/lNarm7PXGfcT/ueHx/97hJG0e0nN0muW16ne6Zv4vM9k6nkzEeT+/D02qNl5P7rJ8E9fWD+zgq3KVkZD8hkqHicRBET55kG9KehJlP8vQnWXKb04ZH6knXHQxvHDnFse0tWWPKDLfsfWYHTPc2PYOlOxy6t+m+Py18soIlo413B1SX9+yz050PPzzJWS+/EPeT2Wdv++3yYv8v6MfeN3JqVmFQYUZFhZEwi5YrjES8bPG6RIVRYfH6CsMMEvFytH3cT2pdtDy8siJubyTifsBIVGT22b0tGcuf+uh3g/t7IsdwN7OlwA+ABPCgu//XrPtvAr4bL7YB/97dXy/kQAdD8kmf+YSnO0RSr7ZpbdJDp7cnd3/bpx21uDu7ftfB8ddb+36SZh0Z9NwHWUcSWUHQ1d1PakwZYdMzJLt6CZLMcfUMlU532to+ZETzz+Px9AzAjPBLO5JJLpec3fn9F06GQPTETwuF7Cd7HDLpwVIRbxMFUWawGFEflRUVnFVpcX+k2ifDKn2/R3//eyZNPKc76LL67G2/6X2mB1oqsNLWRTWQ6is7BLtDMnus3WNJtkn1k94mDthEhfHKyy/zhYULM+usyBpDfNvMCvt/ocCamo4N+j76DXczSwD3AZcDLUCzmT3j7nvTmv0a+Nfu/i9mdgWwHlgwGAP+xYF3+N4vP6TqtV9kBFnGEVp6qPU4ius+qvJSyZOdO864i+wndnqYdB8xdD+Zsp8MmUcZWdvHbXo7WrG0J+CxinYm/tHYVD+ZT0LSjpjiJ3bqqCrtyCrtiZ1ZU+YTurfwSMT9ZuwjOzxS9VtakKSHV7T9tldf5eKLLszpd5PcPjn2UtLU1ER9/axiD6MgPjWignNHn1XsYZSNXI7c64BD7v4mgJk1AlcDqXB39/+b1v4VoKaQg0xXfVaCCaMqmPDp6h5BlgqKrCd1X0FD9vanOUrp7jvz6Cc9ELqDI3P/GUdyWcG0fds2Llwwv8f22WPv3nd32Gb3XWxRkMwp9jAK4u1RFUwZN7LYwxDJm3k/h69mdh2w1N1XxMs3AwvcfVUf7W8HZibbZ923ElgJMGHChHmNjY15DbqtrY3q6uq8ti01qqU0hVJLKHWAaklatGjRdnev7a9dLkfuvR0S9vqKYGaLgD8HvtDb/e6+nmjKhtraWq/P8+3i6Agxv21LjWopTaHUEkodoFoGKpdwbwGmpC3XAK3ZjczsAuBB4Ap3f7cwwxMRkXzk8iGmZmC6mU0zs+HADcAz6Q3M7DPAk8DN7n6g8MMUEZGB6PfI3d07zGwVsIXoVMgN7r7HzG6L718H/GdgPPCj+I29jlzmhEREZHDkdJ67u28GNmetW5d2ewXQ4w1UkaG2accRGrbsp/W9diaNrWL1khkAPdZdM2fykOx7MPaTi+9t2s2jv3qb73z+FH++ZjPLF0zh7mvOL8pYpDj0CVUJxqYdR1jz5G7aT3UCcOS9dlb/4+tgcKrTU+vWPLkboKDB29u+B2M/ufjept08/MpbqeVO99SyAv6TQxcOk2A0bNmfCtekU12eCvak9lOdNGzZP+j7Hoz95OLRX709oPUSJoW7BKP1vfZBaXsm/RV6P7no7OOzK32tlzAp3CUYk8ZWDUrbM+mv0PvJRaKPTyv3tV7CpHCXYKxeMoOqYYmMdcMqjGGJzFCrGpZIvdE6mPsejP3kYvmCKQNaL2HSG6oSjOQbl8U4W6avfRfjbJnkm6bJOfaEmc6W+QRSuEtQrpkzuddAHYqQ7WvfxXD3Nedz9zXn09TUxP+7qb7Yw5Ei0LSMiEiAFO4iIgFSuIuIBEjhLiISIIW7iEiAFO4iIgFSuIuIBEjhLiISIIW7iEiAFO4iIgFSuIuIBEjhLiISIIW7iEiAFO4iIgFSuIuIBEjhLiISIIW7iEiAFO4iIgFSuIuIBEjhLiISIIW7iEiAFO4iIgFSuIuIBEjhLiISIIW7iEiAFO4iIgFSuIuIBEjhLiISoJzC3cyWmtl+MztkZnf0cr+Z2d/G9+8ys7mFH6qIiOSq33A3swRwH3AFcB6w3MzOy2p2BTA9/lkJ3F/gcYqIyADkcuReBxxy9zfd/WOgEbg6q83VwI898gow1swmFnisIiKSo8oc2kwG3k5bbgEW5NBmMvDb9EZmtpLoyB6gzcz2D2i03c4BjuW5balRLaUplFpCqQNUS9Jnc2mUS7hbL+s8jza4+3pgfQ77PP2AzLa5e+2Z9lMKVEtpCqWWUOoA1TJQuUzLtABT0pZrgNY82oiIyBDJJdybgelmNs3MhgM3AM9ktXkGuCU+a+ZC4H13/212RyIiMjT6nZZx9w4zWwVsARLABnffY2a3xfevAzYDXwYOAR8Ctw7ekIECTO2UENVSmkKpJZQ6QLUMiLn3mBoXEZEyp0+oiogESOEuIhKgkg93MxthZq+a2etmtsfMvh+vH2dmz5nZwfjfTxV7rLkws4SZ7TCzn8bL5VrHYTPbbWY7zWxbvK5caxlrZk+Y2T+b2T4zu6gcazGzGfHjkfw5bmbfKdNa/kP8fH/DzB6Nc6Ds6gAws7+M69hjZt+J1w16LSUf7sBHwKXuPguYDSyNz8i5A3jB3acDL8TL5eAvgX1py+VaB8Aid5+ddr5uudbyA+Bn7j4TmEX0+JRdLe6+P348ZgPziE5ueIoyq8XMJgPfBmrd/fNEJ3LcQJnVAWBmnwf+HdEn/WcBXzGz6QxFLe5eNj/ASOA1ok/I7gcmxusnAvuLPb4cxl8TP5CXAj+N15VdHfFYDwPnZK0ru1qAs4FfE59cUM61ZI1/MfBSOdZC9yfexxGd0ffTuJ6yqiMe5/XAg2nLfwX8p6GopRyO3JNTGTuBo8Bz7v4rYILH59LH/366mGPM0VqiB7YrbV051gHRJ5C3mtn2+LISUJ61fA54B/j7eLrsQTMbRXnWku4G4NH4dlnV4u5HgHuBt4guYfK+u2+lzOqIvQFcYmbjzWwk0SnjUxiCWsoi3N2906M/NWuAuvhPnbJiZl8Bjrr79mKPpUAWuvtcoiuCftPMLin2gPJUCcwF7nf3OcAHlMGf+6cTf9jwKuAfiz2WfMTzz1cD04BJwCgz+7Pijio/7r4P+BvgOeBnwOtAx1DsuyzCPcnd3wOagKXA75NXnoz/PVrEoeViIXCVmR0murLmpWb2MOVXBwDu3hr/e5RoXreO8qylBWiJ/xoEeIIo7MuxlqQrgNfc/ffxcrnVchnwa3d/x91PAU8CF1N+dQDg7g+5+1x3vwT4A3CQIail5MPdzM41s7Hx7SqiB/6fiS558LW42deAp4szwty4+xp3r3H3qUR/Mv8fd/8zyqwOADMbZWajk7eJ5kPfoAxrcfffAW+b2Yx41ZeAvZRhLWmW0z0lA+VXy1vAhWY20syM6DHZR/nVAYCZfTr+9zPAnxI9NoNeS8l/QtXMLgD+gegd8wrgcXf/azMbDzwOfIboP8P17v6H4o00d2ZWD9zu7l8pxzrM7HNER+sQTWs84u73lGMtAGY2G3gQGA68SXT5jArKs5aRRG9Gfs7d34/Xld3jEp/yvIxoCmMHsAKopszqADCzfwLGA6eA/+juLwzFY1Ly4S4iIgNX8tMyIiIycAp3EZEAKdxFRAKkcBcRCZDCXUQkQLl8QbbIkIpPE3shXvwjoJPoEgEAde7+cVEGdhpm9m+BzfF58yJFp1MhpaSZ2V1Am7vfWwJjSbh7Zx/3/RJY5e47B9BfpbsPyUfR5ZNH0zJSVszsaxZd33+nmf3IzCrMrNLM3jOzBjN7zcy2mNkCM/uFmb1pZl+Ot11hZk/F9+83s+/l2O/dZvYq0XWNvm9mzfH1uddZZBnR5agfi7cfbmYtaZ+svtDMno9v321mD5jZc0QXK6s0s/8e73uXma0Y+t+qhEjhLmUjvmDctcDF8YXkKoku5QAwBtgaX8zsY+Auoo+tXw/8dVo3dfE2c4EbzWx2Dv2+5u517v4y8AN3nw+cH9+31N0fA3YCyzy6nnp/00ZzgCvd/WZgJdEF5eqA+UQXYftMPr8fkXSac5dychlRAG6LLjlCFdFH7QHa3f25+PZuosvEdpjZbmBqWh9b3P1fAMxsE/AFoudBX/1+TPelFgC+ZGargRHAOcB24H8PsI6n3f1kfHsx8K/MLP3FZDrRR9JF8qZwl3JiwAZ3/6uMlWaVRCGc1EX0DV7J2+n/z7PfZPJ++m33+I2p+Lot/wOY6+5HzOxuopDvTQfdfxlnt/kgq6a/cPcXECkgTctIOXke+KqZnQPRWTV5TGEstug7U0cSXTP8pQH0W0X0YnEsvirmv0m77wQwOm35MNFX3ZHVLtsW4C/iF5Lk96BWDbAmkR505C5lw913x1cLfN7MKoiusncb0DqAbn4JPAL8MfCT5NktufTr7u+a2T8QXd74N8Cv0u7+e+BBM2snmte/C/g7M/sd8OppxvMA0ZUBd8ZTQkeJXnREzohOhZRPjPhMlM+7+3eKPRaRwaZpGRGRAOnIXUQkQDpyFxEJkMJdRCRACncRkQAp3EVEAqRwFxEJ0P8HfLcy7/zjy3oAAAAASUVORK5CYII=\n",
"text/plain": [ "text/plain": [
"<Figure size 432x288 with 1 Axes>" "<Figure size 432x288 with 1 Axes>"
] ]
...@@ -648,7 +650,7 @@ ...@@ -648,7 +650,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 6, "execution_count": 17,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
...@@ -664,6 +666,27 @@ ...@@ -664,6 +666,27 @@
"print(np.sum(data.Malfunction)/np.sum(data.Count))" "print(np.sum(data.Malfunction)/np.sum(data.Count))"
] ]
}, },
{
"cell_type": "code",
"execution_count": 57,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.01270572944054793"
]
},
"execution_count": 57,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"p = np.sum(data.Malfunction)/np.sum(data.Count)\n",
"1-(1-p**2)**3"
]
},
{ {
"cell_type": "markdown", "cell_type": "markdown",
"metadata": {}, "metadata": {},
...@@ -672,7 +695,7 @@ ...@@ -672,7 +695,7 @@
"un joint primaire un joint secondaire sur chacune des trois parties du\n", "un joint primaire un joint secondaire sur chacune des trois parties du\n",
"lançeur, la probabilité de défaillance des deux joints d'un lançeur\n", "lançeur, la probabilité de défaillance des deux joints d'un lançeur\n",
"est de $p^2 \\approx 0.00425$. La probabilité de défaillance d'un des\n", "est de $p^2 \\approx 0.00425$. La probabilité de défaillance d'un des\n",
"lançeur est donc de $1-(1-p^2)^3 \\approx 1.2%$. Ça serait vraiment\n", "lançeur est donc de $1-(1-p^2)^3 \\approx$ 1.2\\%. Ça serait vraiment\n",
"pas de chance... Tout est sous contrôle, le décollage peut donc avoir\n", "pas de chance... Tout est sous contrôle, le décollage peut donc avoir\n",
"lieu demain comme prévu.\n", "lieu demain comme prévu.\n",
"\n", "\n",
...@@ -686,6 +709,305 @@ ...@@ -686,6 +709,305 @@
"analyse et de regarder ce jeu de données sous tous les angles afin\n", "analyse et de regarder ce jeu de données sous tous les angles afin\n",
"d'expliquer ce qui ne va pas." "d'expliquer ce qui ne va pas."
] ]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Nouvelle perspective"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"data = pd.read_csv(\"shuttle.csv\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Tout d'abord, nous allons inspecter graphiquement la relation entre le dysfonctionnement des joints toriques et la température."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/opt/conda/lib/python3.6/site-packages/scipy/stats/stats.py:1713: FutureWarning: Using a non-tuple sequence for multidimensional indexing is deprecated; use `arr[tuple(seq)]` instead of `arr[seq]`. In the future this will be interpreted as an array index, `arr[np.array(seq)]`, which will result either in an error or a different result.\n",
" return np.add.reduce(sorted[indexer] * weights, axis=axis) / sumval\n",
"/opt/conda/lib/python3.6/site-packages/matplotlib/axes/_axes.py:6571: UserWarning: The 'normed' kwarg is deprecated, and has been replaced by the 'density' kwarg.\n",
" warnings.warn(\"The 'normed' kwarg is deprecated, and has been \"\n"
]
},
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fce15975198>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import seaborn as sns\n",
"sns.distplot(data.Temperature.tolist())"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/opt/conda/lib/python3.6/site-packages/matplotlib/axes/_axes.py:6571: UserWarning: The 'normed' kwarg is deprecated, and has been replaced by the 'density' kwarg.\n",
" warnings.warn(\"The 'normed' kwarg is deprecated, and has been \"\n"
]
},
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fce13783710>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl0XOWd5vHvr0qlzdp3WYst7zYGb7IxS4JNWAxJxyQBAoSkQyfNkIRM5/RMT2d6m17SMz2nl0mnQ4cmDIfhJEBnxyEGx4DZAsaWd8vyIsuLZO2Wtdra3/lDMq0IY5Xtkm7p6vmco4OrdKl6KJvHr9773veacw4REfGXgNcBREQk8lTuIiI+pHIXEfEhlbuIiA+p3EVEfEjlLiLiQyp3EREfUrmLiPiQyl1ExIdivHrjrKwsN3PmTK/eXkRkUtqxY0ezcy57rOM8K/eZM2dSVlbm1duLiExKZnYinOM0LSMi4kMqdxERH1K5i4j4kMpdRMSHVO4iIj6kchcR8SGVu4iID6ncRUR8SOUuIuJDnl2hKt569r2TXke4Ig9cW+x1BJGoppG7iIgPqdxFRHxI5S4i4kMqdxERHxqz3M3sKTNrNLP9H/J9M7PvmFmlme01s+WRjykiIpcinJH708C6i3z/DmDu8NfDwPeuPJaIiFyJMcvdOfcm0HKRQ9YDz7ghW4E0M8uPVEAREbl0kZhzLwCqRzyuGX7uA8zsYTMrM7OypqamCLy1iIhcSCTK3S7wnLvQgc65J5xzpc650uzsMW8BKCIilykS5V4DFI14XAjURuB1RUTkMkWi3DcAXxheNbMaaHPO1UXgdUVE5DKNubeMmT0HrAGyzKwG+B9ACMA59ziwEbgTqATOAg+NV1gREQnPmOXunLt/jO874GsRSyQiIldMV6iKiPiQyl1ExIdU7iIiPqRyFxHxIZW7iIgPqdxFRHxI5S4i4kMqdxERH1K5i4j4kMpdRMSHVO4iIj6kchcR8SGVu4iID6ncRUR8SOUuIuJDKncRER9SuYuI+JDKXUTEh1TuIiI+pHIXEfEhlbuIiA+p3EVEfEjlLiLiQyp3EREfUrmLiPiQyl1ExIdU7iIiPqRyFxHxIZW7iIgPqdxFRHworHI3s3VmdsjMKs3smxf4fqqZ/dLM9phZuZk9FPmoIiISrjHL3cyCwGPAHcAi4H4zWzTqsK8BB5xzS4A1wD+aWWyEs4qISJjCGbmvAiqdc1XOuV7geWD9qGMckGxmBiQBLUB/RJOKiEjYwin3AqB6xOOa4edG+i6wEKgF9gF/4JwbjEhCERG5ZOGUu13gOTfq8e3AbmA6sBT4rpmlfOCFzB42szIzK2tqarrksCIiEp5wyr0GKBrxuJChEfpIDwE/c0MqgWPAgtEv5Jx7wjlX6pwrzc7OvtzMIiIyhnDKfTsw18xKhk+S3gdsGHXMSeBjAGaWC8wHqiIZVEREwhcz1gHOuX4zexTYBASBp5xz5Wb2yPD3Hwf+BnjazPYxNI3zx8655nHMLSIiFzFmuQM45zYCG0c99/iIX9cCt0U2moiIXC5doSoi4kMqdxERH1K5i4j4kMpdRMSHwjqhKv5R23qO6paz7KluJS0xRFFGIgG70HVqIjKZqdyniOqWs/yvlyrYuK/+t55PSwixrDiNNfNzCAX1g5yIX6jcp4Cn3j7G3710kGDA+M83z2FlSQZlx89Q23qOvTVtbDnUxOGGTh5YVUz6NG3mKeIHKnefe/KtKr71qwpuWZjLt+5aTF5qPADVLefITYlnWXE6FXXt/Kismu9uqeR3r5tBceY0j1OLyJXSz+E+9v/eOc63flXBx6/O5/EHl79f7KMtzE/h0bVzSIwN8szWE7R09U5wUhGJNJW7T71ztJn/saGc2xbl8u37lhIzxnx6ZlIcv3vdTJyDZ949TnffwMQEFZFxoXL3ofbuPv7ox3uZlTWNf75vWdgnSrOS43jg2mKaO3v4cVk1zo3e2VlEJguVuw/99S8PUNd2jn+8dwkJscFL+ndnZyexbnE+FfUd7KlpG6eEIjLeVO4+89rBBn6yo4avrZ3DsuL0y3qN62dnUpSewIt7a+nq0d0SRSYjlbuP9A8M8j83HmRW9jS+fvPcy36dgBmfWl5Id98AG/fVRTChiEwUlbuP/HRnDZWNnfy32xcQG3Nlv7V5KfHcNC+bXdWtVDV3RiihiEwUlbtPdPcN8H82H2FZcRq3X5UbkddcMz+HlPgYNu2v18lVkUlG5e4TT79znPr2br65bgEWob1iQsEAH1uQS/WZcxys74jIa4rIxFC5+8C53gH+7Y2jrJ2fzbWzMiP62stnpJM5LZZfH6hnUKN3kUlD5e4DP9lZw5mzfXx17ZyIv3YwYNy6KJeG9h72VLdG/PVFZHyo3Ce5wUHHU28fY0lRGqUzLm/p41gWF6SSnxrPlkNNGr2LTBIq90nulYoGjjV38fsfKYnYXPtoATM+Oi+b5s4eDta1j8t7iEhkqdwnuSffOkZBWgLrrsob1/dZPD2V9MQQbxxu0soZkUlA5T6J7atpY9vxFn7vxpIxNwa7UsGAcePcbKrPnOP46bPj+l4icuVU7pPYs9tOEh8KcE9p4YS834ridBJjg7x1pGlC3k9ELp/KfZLq6ulnw+5TfPzq6aTEhybkPWNjAlw3K5OD9R00dfRMyHuKyOVRuU9Sv9pbR1fvAPevKprQ911VkkHQjK3HTk/o+4rIpVG5T1LPbT/JnJwkVozT8scPkxwfYnFBCjtPnKGnXzf0EIlWKvdJ6FB9B7tOtnLfyqJxW/54MdfNyqSnf5DduqhJJGqp3Cehf99eTShofHr5xJxIHa0oI5HpqfG8e/S0lkWKRCmV+yTTPzDIhj2n+NiCXDKmxXqSwcxYPSuTxo4ejp3u8iSDiFycyn2S+c3R0zR39nLXsume5lhSlEZ8KMD2Yy2e5hCRCwur3M1snZkdMrNKM/vmhxyzxsx2m1m5mb0R2Zhy3gu7TpEcH8Oa+Tme5ggFAywtSqO8tp1zvTqxKhJtxix3MwsCjwF3AIuA+81s0ahj0oB/BT7pnLsKuGccsk5553oH2FRez52L84kPXdqNr8dD6YwM+gcdu6vPeB1FREYJZ+S+Cqh0zlU553qB54H1o455APiZc+4kgHOuMbIxBYY2CevqHWC9x1My501PS2B6WjxlJ87oxKpIlAmn3AuA6hGPa4afG2kekG5mr5vZDjP7QqQCyn94Yfcp8lLiWV0S2RtyXInSGRnUtXVT29rtdRQRGSGccr/QQurRw7QYYAXwceB24M/NbN4HXsjsYTMrM7OypibtT3IpWs/28sbhJj65dDqBwMSvbf8wSwrTiAkY20/oxKpINAmn3GuAkde4FwK1FzjmZedcl3OuGXgTWDL6hZxzTzjnSp1zpdnZ2ZebeUr69YEG+gYcn7gm3+sovyUhNsjVBansqW6lt3/Q6zgiMiycct8OzDWzEjOLBe4DNow65gXgI2YWY2aJwLVARWSjTm0v7aujMD2BqwtSvY7yAStmptPTP8j+2javo4jIsDHL3TnXDzwKbGKosH/knCs3s0fM7JHhYyqAl4G9wDbgSefc/vGLPbW0nevj7cpm7rw635PtBsZSkjmNzGmxlB3XqhmRaBETzkHOuY3AxlHPPT7q8d8Dfx+5aHLe5uEpmTuvjq4pmfPMjNIZ6Ww60EBzRw9ZyXFeRxKZ8nSF6iSwcV8dBWkJLCmMvimZ85bNSCdgUHZCo3eRaKByj3Lt3X28daSJOxbnReWUzHkp8SHm5yaz8+QZBga15l3Eayr3KPdqxfCUTJStkrmQ0pkZdPb0c6i+w+soIlOeyj3KbdrfQG5KHEsL07yOMqZ5uckkxcWw46SmZkS8pnKPYt19A7x5pIlbFuZG1YVLHyYYMJYVp3Govp2O7j6v44hMaSr3KPbO0WbO9g5w66Jcr6OEbUVxOoMO3aVJxGMq9yj26/IGkuJiuG529OwlM5aclHiK0hPYoc3ERDylco9Sg4OOVyoauWl+NnEx3m/veymWz0insaOHU63nvI4iMmWp3KPUrupWmjt7uG0STcmcd34zMa15F/GOyj1K/fpAPTEB8/yOS5cjPhRkcUEqe2ta6RvQZmIiXlC5R6nNBxpYPSuT1ISQ11Euy/LidLr7BjlQ2+51FJEpSeUehSobO6lq6ppUq2RGm5U9jbTEEDs0NSPiCZV7FNp8oAGAWyZxuQfMWF6cztGmTlrP9nodR2TKUblHoc0H6llckEJBWoLXUa7I8uJ0HLBTV6yKTDiVe5Rp7OhmV3Urty7M8zrKFcuYFsusrGnsOHGGQa15F5lQKvco82pFI84xqefbRyqdmc6Zs31UNXV5HUVkSlG5R5nNBxooSEtgYX6y11Ei4qrpqcSHAmw/rhtoi0wklXsU6erp5+3KZm67Kjeq926/FKFggGVF6Ryoa6erp9/rOCJThso9irx1pIne/kFuXeiPKZnzVs7MYGDQsUsnVkUmjMo9imw+0EhKfAwrSzK8jhJRealDm4lt12ZiIhNG5R4lBgYdWw41snZBDqGg/35bVs7MoKmjhxOnz3odRWRK8F+LTFK7Tp6hpauXW3w2JXPeNYVpxMUE2KYTqyITQuUeJTZXNBATMG6an+11lHERGxNgeXE6+0610akTqyLjTuUeJV4Z3igsJX5ybhQWjmtLhk6slmn0LjLuVO5RoKqpk6NNXdyycPJt73spclLimZU9jW3HWnTFqsg4U7lHgVcrGgH4mE/n20daXZJJ67k+DtV3eB1FxNdU7lFgc0UDC/KSKcpI9DrKuFuYn0JKfAxbq057HUXE11TuHjvT1UvZ8RbfrpIZLRgwVpVkcqSxk4b2bq/jiPiWyt1jWw41Mugm997tl+rakgxCQeM3lc1eRxHxLZW7x16paCA7OY5rClK9jjJhpsXFsKw4nd3VrVoWKTJOVO4e6ukf4M3DzdyyMIdAwB8bhYXrhtlZ9A863tPcu8i4CKvczWydmR0ys0oz++ZFjltpZgNmdnfkIvrXe1UtdPb0T5n59pGyk+NYkJfM1qrT9A0Meh1HxHfGLHczCwKPAXcAi4D7zWzRhxz3v4FNkQ7pV69UNBAfCnDDnCyvo3jixjlZdPUO6CbaIuMgnJH7KqDSOVflnOsFngfWX+C4rwM/BRojmM+3nHO8cqCBj8zNJj4U9DqOJ0qyplGckcgbh5voH9ToXSSSwin3AqB6xOOa4efeZ2YFwKeAxy/2Qmb2sJmVmVlZU1PTpWb1lQN17dS2dftu7/ZLYWbcvCCHtnN97DrZ6nUcEV8Jp9wvdKZv9LXj3wb+2Dk3cLEXcs494Zwrdc6VZmf7c4OscL1yoBEzWLvA31sOjGVuThIFaQm8cbiJgUFtSSASKeGUew1QNOJxIVA76phS4HkzOw7cDfyrmd0VkYQ+9UpFA8uK0shOjvM6iqfOj95bunrZU6PRu0ikhFPu24G5ZlZiZrHAfcCGkQc450qcczOdczOBnwBfdc79IuJpfaKu7Rz7TrVNib1kwrEgL5npqfG8WtFAv1bOiETEmOXunOsHHmVoFUwF8CPnXLmZPWJmj4x3QD86v1HYrVPoqtSLMTNuX5zHmbN9vHdM2wGLREJMOAc55zYCG0c9d8GTp865L155LH97taKB4oxE5uYkeR0laszNSWZOThKvHWxkeXE6CbFTcwWRSKToCtUJ1tXTz2+OnuaWhbmYTa2rUsey7qo8zvUN8OaRqb2SSiQSVO4T7K0jTfT2D3LLoqm9SuZCpqclsLQojd9UNnO6s8frOCKTmsp9gr28v560xBArZ2Z4HSUq3X5VHsGAsWFPLU53axK5bCr3CdTTP8CrFY3cujCXUFAf/YWkJoS4dVEuRxo72Xeqzes4IpOWGmYC/aaymY6efu68Ot/rKFFt9axMCtIS+NXeOs71XvS6OBH5ECr3CfTSvnqS42K4fk6m11GiWsCMu5YW0NnTz4t7R18vJyLhULlPkL6BQTZXNHDLolziYrTMbywF6QmsXZDDrupWdldr10iRS6VynyBbq07TeraPdYvzvI4yaaydn8OMjERe2F1LS1ev13FEJhWV+wR5aX89ibFBbpo3tTdMuxTBgHHvyiLM4LltJ+nt19YEIuEK6wpVuTIDg45fl9ezdkHOlN27/XKlJ8Zyz4oifrD1BD/eUc39q4oJmPHseye9jnZFHri22OsI4nMauU+A7cdbaO7s5c7FWiVzORbmp3Dn1fmU17azqbze6zgik4JG7hPg5f31xMUEWDNfUzKX6/rZmZzu6uGtI83ExQRYOz9H2zeIXITKfZwNDjpe2l/HmvnZTIvTx325zIxPXDOdnr5BXqlopG/Acdsi7c8j8mHUNuNsV3UrDe093KEpmSsWMOMzKwoJBQO8cbiJrp5+PrlkOjG62lfkA1Tu4+zl/XXEBgPcvFAbhUVCwIz1S6eTGBvk9cNN1Ld388CqYtISY72OJhJVNOQZR845Nu6r58a5WaTEh7yO4xtmxm1X5fHgtcU0dfTwL69Vsv1YC4PaaEzkfSr3cbS3po1Tred04dI4WTQ9lUfXziEvNZ6f7z7FE29WceJ0l9exRKKCyn0cbdhTS2wwwO1XqdzHS2ZSHF++sYS7lxdyurOHf3uziiffruJwQ4dG8jKlac59nAwMOn65p5a1C7JJTdCUzHgyM5bPSGdxQSrvHTvNW0eaefqd46QlhlhenM5V01PIS4nXyhqZUlTu42Rr1WkaO3pYv7TA6yhTRmxMgI/MzWb1rEwO1LWz4/gZXjvYyGsHG0lNCDE/N5kFecnMyk4iNkY/tIq/qdzHyS92nSI5LoabF2iVzEQLBQMsKUxjSWEa7d19HK7v4GB9B7urW9l2vIWYgDEjM5E5wzflzk+NJ6BRvfiMyn0cdPcN8PL+em5fnKe9ZDyWEh+idGYGpTMz6B8Y5NjpLg7Xd1DZ1Mmm8no2lUNibJB5uckszE9hXk4Scfo9Ex9QuY+D1w810tHTz12akokqMcEAc3OSmZuTDEB7dx9HGzs50tjJ4YahkX1MwFiQn8KyojTm5SYTDGhEL5OTyn0c/HTnKbKT47hutu64FM1S4kMsK05nWXE6A4OOky1n2XeqjX01rew/1UZyfAyrZmawsiRD1ynIpKNyj7Dmzh62HGzkSzeWaNQ3iQQDRknWNEqypvHxq/M53NDBe8dO8+rBRl4/3ETpjHRumpetK2Fl0lC5R9gvdp2if9Bx94pCr6PIZQoGjIX5KSzMT6G5c2gnyrLjZyg7fobVszJYuyCHxFj9ryPRTX9CI8g5x4/LalhSlMbc3GSv40gEZCXF8allBaydn82WQ428c/Q0O0+28rGFOayelalVNhK1tNg3gvadauNQQwf3aNTuO2mJsXxqWSFf/9hcCtITeHFvHY+/cZS6tnNeRxO5IJV7BP24rIa4mAC/s2S611FknOSlxPPQ9TP5bGkRZ8728diWSl4/1KitDiTqaFomQrr7Btiwp5bbr8rTdgM+Z2bDU29JvLC7ll8faOBIYyf3rCjUCVeJGmGN3M1snZkdMrNKM/vmBb7/OTPbO/z1jpktiXzU6Pbi3jrazvVx36oir6PIBEmMjeG+lUXcvbyQU2fO8S+vVbL/VJvXsUSAMMrdzILAY8AdwCLgfjNbNOqwY8BNzrlrgL8Bnoh00Gj3g60nmJ09jetmaW37VHJ+07JHb55DxrRYnt12kp/vOkXfwKDX0WSKC2fkvgqodM5VOed6geeB9SMPcM6945w7M/xwKzClzijuP9XG7upWHlw9QzsPTlFZSXH8p5tm8dG52Ww/3sL336qi7Vyf17FkCgun3AuA6hGPa4af+zBfAl66klCTzQ+2niAhFOTTy6fU32kySkwgwLrFeTywqpjGjh4e21LJ8WbdPES8EU65X2goesGlAWa2lqFy/+MP+f7DZlZmZmVNTU3hp4xi7d19vLC7lk8uma4TqQLA4oJUvnLTbOJiAjz5dhVbq07jtJpGJlg45V4DjDxLWAjUjj7IzK4BngTWO+dOX+iFnHNPOOdKnXOl2dnZl5M36vy4rIZzfQM8uHqG11EkiuSmxPPVNXOYm5PMhj21/HzXKfo1Dy8TKJxy3w7MNbMSM4sF7gM2jDzAzIqBnwGfd84djnzM6NQ/MMhTbx9j5cx0ri5M9TqORJmE2CCfv24Ga+ZnU3bijObhZUKNWe7OuX7gUWATUAH8yDlXbmaPmNkjw4f9BZAJ/KuZ7TazsnFLHEVeLq/nVOs5vvyRWV5HkSgVMOO2RUPz8A3tmoeXiRPWRUzOuY3AxlHPPT7i118GvhzZaNHNOcf33zrGzMxEblmY63UciXKLC1LJTo7jB1tP8OTbVRRlJGh1lYwrbT9wmcpOnGFPdau29pWwjZyH//MXyvmjn+ylu2/A61jiU9p+4DI98WYVaYkh7l6hK1IlfOfn4Rs7evjOq0c43NDB4w+uYHpagtfRxGc0cr8M5bVtbD7QwO9eN5OEWN1vUy5NwIw/vHUe3/9CKVVNXXz8O2/xyoEGr2OJz6jcL8N3Xj1CcnwMv3djiddRZBK7dVEuGx69gfzUBL78TBl/8cJ+TdNIxKjcL9GB2nY2lTfw0A0lumhJrtis7CR+/rXr+fKNJTzz7gnWf/c3HKrv8DqW+IDK/RJ959UjJMfF8KUbNGqXyIiLCfJnn1jE0w+t5HRXD5/87ts89fYxBgZ1VatcPpX7Jdh/qo2Xy+t56MYSUhM1apfIWjM/h5f+4KNcPzuTv37xAJ/53jscrG/3OpZMUir3MDnn+NtfVZCeGOJLmmuXcZKdHMdTX1zJtz+7lJMtZ/n4d97mLzeU03q21+toMsmo3MP0akUj71ad5hu3zNNcu4wrM+OuZQW88oc3cd/KIp559zg3/f3rPP7GUbp6+r2OJ5OEebVbXWlpqSsrmxy7FPQNDHL7t98EYNM3PkooOPn/Tnz2vZNeR5Aw1bd183J5HYcbOkmMDXLDnCxWzswgKc67y1QeuLbYs/ee6sxsh3OudKzjdBFTGJ597yRVTV18/wulvih2mVzyUuP54vUlnGw5y2sHG9h8oIEtBxu5pjCVZcXplGRNI6BtDGQUlfsYGtq7+YdNh7hxTha3LMzxOo5MYcUZiXzx+hIa2rvZWnWaXdWt7DzZSnJ8DIvyU5ifm0xJ9jTiYnRhnajcx/SXG8rpHRjkW3ct1iZPEhVyU+JZv7SAOxbnc7C+nT01bew8eYb3jrUQsKGRflF6Ijkp8WROiyVzWixpibHaA2mKUblfxOYDDby0v54/un0+M7OmeR1H5LfExgS4pjCNawrT6B8Y5ETLWY42dlJ95iy7q1vp6f+Pm4MEDFITQkyLiyEhFCQhNjj0z1CQuFCQuJjA8FeQuNAHfx0bDGhwM8mo3D9E29k+/uKF/czPTebhj2q/doluMcEAs7OTmJ2dBAwt3e3o6aels5fTXb20dPVw5mwfZ3v7Ods7QEtXL+f6BujuGyCca6VCQSMlPkRqYojU+BC1reeYk5PEgvxkZmcn6VxUFFK5X4Bzjm/+bC9NHT382+dX6A+uTDpmQ2WcEh+66E+dzjn6Bx3dfQP09g/S3T9IT/8APX2D9Iz4dWdPP23n+mg710dVcxd7T7W9fwVtbDDAnJwkFuancO2sDG6Yk0WBdrn0nMr9Ap7fXs1L++v55h0LuKYwzes4IuPGzAgF7ZIHMPeUFlLV1EVFXTsV9e1U1HWw5VAjP91ZA8CMzESun53FR+ZmsWZ+NomxqpqJpk98lCMNHfzVL8u5YU4mD+v2eSIXFAoGmJ+XzPy8ZO6iABj6KeBQQwfvVJ7mnaPNvLinlue2nSQ+FGDNvBzuvCafmxfkeLo+fyrRpzzCma5evvxMGUlxMfzTvUsJaHWBSNjMjAV5KSzIS+H3biyhf2CQshNneGlfHS/tr+fl8nriYgLcvCCHe0uL+MjcLGI05TluVO7D+gYG+coPd1DX2s1zD68mNyXe60gik1pMMMDqWZmsnpXJX/zOVew4cYaN++rYsKeWl/bXk5sSx90rCrlnRZFWo40DlTtDP07+2c/3s7WqhX+6dwkrZqR7HUnEV4IBY1VJBqtKMviTOxfy2sEGflRWw/deP8pjW46yelYGD66ewW2L8oiN0Wg+EqZ8uTvn+KtfHuDfy6r5+s1z+PTyQq8jifhabEyAdYvzWbc4n4b2bn6yo4bntp3k0Wd3kZUUy72lRdy/qpiijESvo05qU7rcnXP83UsHefqd43zpxhL+8NZ5XkcSmVJyU+L52to5fOWm2bx5pIkfbD3J428c5XtvHOWmedl87toZ3LwgR1fXXoYpW+79A4P8+QvlPLftJA+uLubPPr5QV+CJeCQQMNbMz2HN/BxqW8/x/PZqnt92kt9/pozpqfHct6qYz64s0rmwSzAlt/zt6unn0Wd3suVQE19dM5v/etv8KbcyRlv+ypWYiC1/+wYGebWikR++d4K3jjQTDBi3Lszlc6uLuWF21pT7f/Y8bfn7ISrq2nn02Z0ca+7ibz+1mM9dO8PrSCJyAaFggHWL81i3OI/jzV08t+0kPyqr5uXyemZmJvLAtcXcvaKIjGmxXkeNSlNm5D446Pjheyf4m19VkJoQ4tufXcoNc7Im7P2jjUbuciW8ullHd98Am8rr+cHWE2w/fobYYIA7r87jc6tnUDojfUpMrWrkPkJFXTt/+vN97DzZykfnZfOP9ywhOznO61giconiQ0HWLy1g/dICDtV38Ox7J/jZzlP8YnctxRmJfOKafH5nyXQW5CVPiaK/GF+P3KtbzvLd1yr5yc4aUhNC/MmdC/nM8oIp/5sOGrnLlYmm2+yd7e3nxb11/HJPLe8cPc3AoGN29jQ+cc10bl2Uy6L8FF/Nz0/pkfvemlaeefcEv9h1ikDA+PzqGXzjlrmkJWpuTsRvEmNjuLe0iHtLizjd2cNL++v55Z5avvPaEf751SNkJ8exZl42Ny/I4Ya5WaTET40b3Pum3KtbzvLy/no27Kll36k2EmODfO7aYr6yZg55qVo+JTIVZCbF8eDqGTy4egbNnT28caiJLYca2VRez4931BAwWJifwsqZGZTOTGdU0LRnAAAE9ElEQVTlzAzfLq8Mq9zNbB3wz0AQeNI593ejvm/D378TOAt80Tm3M8JZ3+eco66tm701bWytOs27R09zqKEDgMUFKfz1+qv41LICkqfI39Ai8kFZSXF8ZkUhn1lRSP/AIDtPtvJ2ZTM7TrTw79urefqd4wDkp8azIC+ZBfkpLMhLZmF+CjMyEyf9vWjHLHczCwKPAbcCNcB2M9vgnDsw4rA7gLnDX9cC3xv+Z8RtPtDAf//ZXpo7ewGIDwVYOTODz6woYN1V+RRn6pJlEfltMcHA+3vbwNAa+oq6drYda6G8tp2Kunbermymb2DoHKQZ5CbHU5ieQFFGIkXpCWQlx5GeGEvGtFjSE2NJSwwN3Y4wFCQ2GCAUNMwM5xy9A4P0DTj6+gdp7+6jubOHpo5emjt7aO7sYXlxOh+dlz2+/81hHLMKqHTOVQ39R9vzwHpgZLmvB55xQ2dnt5pZmpnlO+fqIh24IC2BNfNzuKYwlcUFqVw1PWXS/w0rIhMrFPyP+8+e19s/SFVzJwfrOjh+uovqlnPUnDnLtmMtvLD73Ji3IzSDoBn9Ydy38CtrZkdFuRcA1SMe1/DBUfmFjikAIl7ui6an8A/3LIn0y4rIFBcbE3h/P/rR+gcGaT3Xx5muXlq6ejlztpfWs33v34qwt3/otoQDg45QMEDs8E3FQ0EjKT5EVlIsWUlxZCXFkTEtdkJ2vgyn3C+0hmj0X03hHIOZPQw8PPyw08wOhfH+lysLaB7H15/s9PlcnD6fi/icPp+xjOfnE9Zl9eGUew1QNOJxIVB7GcfgnHsCeCKcYFfKzMrCWQs6VenzuTh9Phenz+fiouHzCedng+3AXDMrMbNY4D5gw6hjNgBfsCGrgbbxmG8XEZHwjDlyd871m9mjwCaGlkI+5ZwrN7NHhr//OLCRoWWQlQwthXxo/CKLiMhYwlrn7pzbyFCBj3zu8RG/dsDXIhvtik3I9M8kps/n4vT5XJw+n4vz/PPxbG8ZEREZP7oTrYiID/m23M3sHjMrN7NBM9NZ/WFmts7MDplZpZl90+s80cbMnjKzRjPb73WWaGNmRWa2xcwqhv/f+gOvM0UTM4s3s21mtmf48/krL/P4ttyB/cCngTe9DhItRmwlcQewCLjfzBZ5myrqPA2s8zpElOoH/otzbiGwGvia/vz8lh7gZufcEmApsG549aAnfFvuzrkK59x4XiQ1Gb2/lYRzrhc4v5WEDHPOvQm0eJ0jGjnn6s5vCOic6wAqGLoSXRhaWOKc6xx+GBr+8uykpm/LXS7ow7aJELkkZjYTWAa8522S6GJmQTPbDTQCm51znn0+k3o/dzN7Bci7wLf+1Dn3wkTnmQTC2iZC5GLMLAn4KfAN51y713miiXNuAFhqZmnAz81ssXPOk/M3k7rcnXO3eJ1hkglrmwiRD2NmIYaK/YfOuZ95nSdaOedazex1hs7feFLumpaZWsLZSkLkgoZvyvN/gQrn3D95nSfamFn28IgdM0sAbgEOepXHt+VuZp8ysxrgOuBXZrbJ60xec871A+e3kqgAfuScK/c2VXQxs+eAd4H5ZlZjZl/yOlMUuQH4PHCzme0e/rrT61BRJB/YYmZ7GRpIbXbOvehVGF2hKiLiQ74duYuITGUqdxERH1K5i4j4kMpdRMSHVO4iIj6kchcR8SGVu4iID6ncRUR86P8D7VKw1Jfd4zUAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sns.distplot(data.Malfunction.tolist())"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PathCollection at 0x7fce13693c88>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFNNJREFUeJzt3W+MXXd95/H3pxNHNWm2KXigxLFxQJa1dAshuorJpoJE3SQOKusgtVpHLFSo1EtFpKVaWZtUCLq7XWkld/cBu2kik3oDW3DULYnxA8DhwWrDQhN5TNI4CbhrQsDjYclAmlDAUmLvdx/MmezN+I7nzMwdz5/zfklXc8/v/M6539+cuR8fnzl3fqkqJEnd8QvLXYAk6cIy+CWpYwx+SeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjrlouQsYZMOGDbVly5blLkOSVo2jR4/+qKpG2/RdkcG/ZcsWxsbGlrsMSVo1knyvbV8v9UhSxxj8ktQxBr8kdYzBL0kdY/BLUscY/JLUMXMGf5JNSf5Hkm8leSrJvxzQJ0k+leREkieSXN23bkeS4826O4Y9AEnS/LS5j/8M8K+q6ptJLgWOJvlqVT3d1+cWYGvz2A7cDWxPMgLcBdwIjANHkhyase1QHHzsFHsPH2fihdNcftl69ty8jVvfsXHYLyNpyHzvXnhzBn9V/QD4QfP875N8C9gI9If3TuCzNTWB7yNJLkvyRmALcKKqngFIcn/Td6jBf/CxU9z5wDFOv3wWgFMvnObOB44B+AMkrWC+d5fHvK7xJ9kCvAN4dMaqjcDJvuXxpm229qHae/j4Kz84006/fJa9h48P+6UkDZHv3eXROviT/BLwBeBjVfWTmasHbFLnaR+0/91JxpKMTU5Oti0LgIkXTs+rXdLK4Ht3ebQK/iTrmAr9z1XVAwO6jAOb+pavACbO036OqtpXVb2q6o2Otvo7Q6+4/LL182qXtDL43l0ebe7qCfDnwLeq6j/N0u0Q8MHm7p53Ai82vxs4AmxNcmWSi4FdTd+h2nPzNtavG3lV2/p1I+y5eduwX0rSEPneXR5t7uq5DvgAcCzJ403bHwGbAarqHuBLwHuAE8DPgQ81684kuR04DIwA+6vqqaGOgP//SyDvDJBWF9+7yyNTN+KsLL1er/yzzJLUXpKjVdVr09dP7kpSxxj8ktQxBr8kdYzBL0kdY/BLUscY/JLUMQa/JHWMwS9JHWPwS1LHGPyS1DEGvyR1jMEvSR1j8EtSxxj8ktQxBr8kdcycE7Ek2Q/8FvBcVf2jAev3AO/v298/BEar6vkkzwJ/D5wFzrT9W9GSpKXT5oz/PmDHbCuram9VXVVVVwF3Av+zqp7v63JDs97Ql6QVYM7gr6qHgefn6te4DTiwqIokSUtqaNf4k7yGqf8ZfKGvuYCHkhxNsntYryVJWrg2k6239V7g6zMu81xXVRNJXg98Ncm3m/9BnKP5h2E3wObNm4dYliSp3zDv6tnFjMs8VTXRfH0OeBC4ZraNq2pfVfWqqjc6OjrEsiRJ/YYS/El+GXg38MW+tkuSXDr9HLgJeHIYrydJWrg2t3MeAK4HNiQZBz4JrAOoqnuabu8DHqqqn/Vt+gbgwSTTr/P5qvrK8EqXJC3EnMFfVbe16HMfU7d99rc9A7x9oYVJkpaGn9yVpI4x+CWpYwx+SeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4JaljDH5J6hiDX5I6xuCXpI4x+CWpYwx+SeoYg1+SOmbO4E+yP8lzSQZOm5jk+iQvJnm8eXyib92OJMeTnEhyxzALlyQtTJsz/vuAHXP0+VpVXdU8/i1AkhHgLuAW4K3AbUneuphiJUmLN2fwV9XDwPML2Pc1wImqeqaqXgLuB3YuYD+SpCEa1jX+a5P8TZIvJ/m1pm0jcLKvz3jTNlCS3UnGkoxNTk4OqSxJ0kzDCP5vAm+qqrcD/xk42LRnQN+abSdVta+qelXVGx0dHUJZkqRBFh38VfWTqvpp8/xLwLokG5g6w9/U1/UKYGKxrydJWpxFB3+SX02S5vk1zT5/DBwBtia5MsnFwC7g0GJfT5K0OBfN1SHJAeB6YEOSceCTwDqAqroH+G3gD5KcAU4Du6qqgDNJbgcOAyPA/qp6aklGIUlqLVMZvbL0er0aGxtb7jIkadVIcrSqem36+sldSeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4JaljDH5J6hiDX5I6xuCXpI4x+CWpYwx+SeoYg1+SOsbgl6SOmTP4k+xP8lySJ2dZ//4kTzSPbyR5e9+6Z5McS/J4Ev/AviStAG3O+O8Ddpxn/XeBd1fV24B/B+ybsf6Gqrqq7QQBkqSlNefUi1X1cJIt51n/jb7FR5iaVF2StEIN+xr/7wFf7lsu4KEkR5PsHvJrSZIWYM4z/raS3MBU8P9GX/N1VTWR5PXAV5N8u6oenmX73cBugM2bNw+rLEnSDEM540/yNuBeYGdV/Xi6vaommq/PAQ8C18y2j6raV1W9quqNjo4OoyxJ0gCLDv4km4EHgA9U1d/2tV+S5NLp58BNwMA7gyRJF86cl3qSHACuBzYkGQc+CawDqKp7gE8ArwP+LAnAmeYOnjcADzZtFwGfr6qvLMEYJEnz0OauntvmWP9h4MMD2p8B3n7uFpKk5eQndyWpYwx+SeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4JaljDH5J6hiDX5I6xuCXpI4x+CWpYwx+SeoYg1+SOmbO4E+yP8lzSQbOnpUpn0pyIskTSa7uW7cjyfFm3R3DLFyStDBtJlu/D/gvwGdnWX8LsLV5bAfuBrYnGQHuAm4ExoEjSQ5V1dOLLXq1O/jYKfYePs7EC6e5/LL17Ll5G7e+Y+Oi+642a3lsS+HjB49x4NGTnK1iJOG27Zv4k1t/fbnL0irUZgauh5NsOU+XncBnq6qAR5JcluSNwBbgRDMTF0nub/p2OvgPPnaKOx84xumXzwJw6oXT3PnAMYBzQm8+fVebtTy2pfDxg8f4i0e+/8ry2apXlg1/zdcwrvFvBE72LY83bbO1d9rew8dfCbtpp18+y97DxxfVd7VZy2NbCgcePTmvdul8hhH8GdBW52kfvJNkd5KxJGOTk5NDKGtlmnjhdOv2+fRdbdby2JbC2Rr81pmtXTqfYQT/OLCpb/kKYOI87QNV1b6q6lVVb3R0dAhlrUyXX7a+dft8+q42a3lsS2Ekg86jZm+XzmcYwX8I+GBzd887gRer6gfAEWBrkiuTXAzsavp22p6bt7F+3cir2tavG2HPzdsW1Xe1WctjWwq3bd80r3bpfOb85W6SA8D1wIYk48AngXUAVXUP8CXgPcAJ4OfAh5p1Z5LcDhwGRoD9VfXUEoxhVZn+xWWbu1nm03e1WctjWwrTv8D1rh4NQ2oFXiPs9Xo1Nja23GVI0qqR5GhV9dr09ZO7ktQxBr8kdYzBL0kdY/BLUscY/JLUMQa/JHWMwS9JHWPwS1LHGPyS1DEGvyR1jMEvSR1j8EtSxxj8ktQxBr8kdYzBL0kdY/BLUse0Cv4kO5IcT3IiyR0D1u9J8njzeDLJ2SSvbdY9m+RYs87ZVSRpmbWZenEEuAu4kakJ1I8kOVRVT0/3qaq9wN6m/3uBP6yq5/t2c0NV/WiolUuSFqTNGf81wImqeqaqXgLuB3aep/9twIFhFCdJGr42wb8RONm3PN60nSPJa4AdwBf6mgt4KMnRJLtne5Eku5OMJRmbnJxsUZYkaSHaBH8GtM02Q/t7ga/PuMxzXVVdDdwCfDTJuwZtWFX7qqpXVb3R0dEWZUmSFqJN8I8Dm/qWrwAmZum7ixmXeapqovn6HPAgU5eOJEnLpE3wHwG2JrkyycVMhfuhmZ2S/DLwbuCLfW2XJLl0+jlwE/DkMAqXJC3MnHf1VNWZJLcDh4ERYH9VPZXkI836e5qu7wMeqqqf9W3+BuDBJNOv9fmq+sowByBJmp9UzXa5fvn0er0aG/OWf0lqK8nRquq16esndyWpYwx+SeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4JaljDH5J6hiDX5I6xuCXpI4x+CWpYwx+SeoYg1+SOqZV8CfZkeR4khNJ7hiw/vokLyZ5vHl8ou22kqQLa84ZuJKMAHcBNzI1/+6RJIeq6ukZXb9WVb+1wG0lSRdImzP+a4ATVfVMVb0E3A/sbLn/xWwrSVoCbYJ/I3Cyb3m8aZvp2iR/k+TLSX5tnttKki6QOS/1ABnQNnOi3m8Cb6qqnyZ5D3AQ2Npy26kXSXYDuwE2b97coixJ0kK0OeMfBzb1LV8BTPR3qKqfVNVPm+dfAtYl2dBm27597KuqXlX1RkdH5zEESdJ8tAn+I8DWJFcmuRjYBRzq75DkV5OkeX5Ns98ft9lWknRhzXmpp6rOJLkdOAyMAPur6qkkH2nW3wP8NvAHSc4Ap4FdVVXAwG2XaCySpBYylc8rS6/Xq7GxseUuQ5JWjSRHq6rXpq+f3JWkjjH4JaljDH5J6hiDX5I6xuCXpI4x+CWpYwx+SeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4JaljDH5J6hiDX5I6plXwJ9mR5HiSE0nuGLD+/UmeaB7fSPL2vnXPJjmW5PEkzq4iSctszqkXk4wAdwE3MjV5+pEkh6rq6b5u3wXeXVV/l+QWYB+wvW/9DVX1oyHWLUlaoDZn/NcAJ6rqmap6Cbgf2Nnfoaq+UVV/1yw+Alwx3DIlScPSJvg3Aif7lsebttn8HvDlvuUCHkpyNMnu2TZKsjvJWJKxycnJFmVJkhZizks9QAa0DZyhPckNTAX/b/Q1X1dVE0leD3w1yber6uFzdli1j6lLRPR6vZU3A7wkrRFtzvjHgU19y1cAEzM7JXkbcC+ws6p+PN1eVRPN1+eAB5m6dCRJWiZtgv8IsDXJlUkuBnYBh/o7JNkMPAB8oKr+tq/9kiSXTj8HbgKeHFbxkqT5m/NST1WdSXI7cBgYAfZX1VNJPtKsvwf4BPA64M+SAJypqh7wBuDBpu0i4PNV9ZUlGYkkqZVUrbzL6b1er8bGvOVfktpKcrQ54Z6Tn9yVpI4x+CWpYwx+SeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4JaljDH5J6hiDX5I6xuCXpI4x+CWpYwx+SeqYVsGfZEeS40lOJLljwPok+VSz/okkV7fdVpJ0Yc05A1eSEeAu4Eam5t89kuRQVT3d1+0WYGvz2A7cDWxvua20LA4+doq9h48z8cJpLr9sPXtu3sat79g4sO/7P/3XfP07z7+yfN1bXsvnfv/aRe93Pn0/fvAYBx49ydkqRhJu276JP7n11+c56gtT61LV4PdgOOYMfqYmRz9RVc8AJLkf2An0h/dO4LM1NZ3XI0kuS/JGYEuLbaUL7uBjp7jzgWOcfvksAKdeOM2dDxwDOOcNNzP0Ab7+ned5/6f/+pzwn89+59P34weP8RePfP+V5bNVrywvNPiWqtalqsHvwfC0udSzETjZtzzetLXp02Zb6YLbe/j4K2+0aadfPsvew8fP6Tsz9M/XPp/9zqfvgUdPntN2vvY2lqrWparB78HwtAn+DGibOVHvbH3abDu1g2R3krEkY5OTky3KkhZu4oXT82pfiv3Op+/ZWebGnq29jaWqdalq8HswPG2CfxzY1Ld8BTDRsk+bbQGoqn1V1auq3ujoaIuypIW7/LL182pfiv3Op+9IBp1Dzd7exlLVulQ1+D0YnjbBfwTYmuTKJBcDu4BDM/ocAj7Y3N3zTuDFqvpBy22lC27PzdtYv27kVW3r142w5+Zt5/S97i2vHbiPQe3z2e98+t62fdM5bedrb2Opal2qGvweDM+cv9ytqjNJbgcOAyPA/qp6KslHmvX3AF8C3gOcAH4OfOh82y7JSKR5mP6lWZs7KT73+9e2vqtnPvudT9/pX14O846Wpap1qWrwezA8qUVcH1sqvV6vxsbGlrsMSVo1khytql6bvn5yV5I6xuCXpI4x+CWpYwx+SeoYg1+SOsbgl6SOWZG3cyaZBL633HXMsAH40XIXsQQc1+qzVsfmuBbnTVXV6s8erMjgX4mSjLW9R3Y1cVyrz1odm+O6cLzUI0kdY/BLUscY/O3tW+4ClojjWn3W6tgc1wXiNX5J6hjP+CWpYwz+AZI8m+RYkseTjDVtf5zkVNP2eJL3LHed89XMhfxXSb6d5FtJrk3y2iRfTfK/m6+/stx1LsQsY1vVxyzJtr7aH0/ykyQfW+3H7DzjWtXHCyDJHyZ5KsmTSQ4k+cWVeLy81DNAkmeBXlX9qK/tj4GfVtWfLlddi5XkM8DXqureZmKc1wB/BDxfVf8hyR3Ar1TVv17WQhdglrF9jFV+zKYlGQFOAduBj7IGjhmcM64PsYqPV5KNwP8C3lpVp5P8JVNzlbyVFXa8POPviCT/AHgX8OcAVfVSVb0A7AQ+03T7DHDr8lS4cOcZ21rym8B3qup7rIFj1qd/XGvBRcD6JBcxdfIxwQo8Xgb/YAU8lORokt197bcneSLJ/pXw37V5ejMwCfzXJI8luTfJJcAbmmkyab6+fjmLXKDZxgar+5j12wUcaJ6vhWM2rX9csIqPV1WdAv4U+D7wA6amoH2IFXi8DP7Brquqq4FbgI8meRdwN/AW4CqmDup/XMb6FuIi4Grg7qp6B/Az4I7lLWloZhvbaj9mADSXrv4p8N+Xu5ZhGjCuVX28mn+odgJXApcDlyT558tb1WAG/wBVNdF8fQ54ELimqn5YVWer6v8CnwauWc4aF2AcGK+qR5vlv2IqLH+Y5I0Azdfnlqm+xRg4tjVwzKbdAnyzqn7YLK+FYwYzxrUGjtc/Ab5bVZNV9TLwAPCPWYHHy+CfIcklSS6dfg7cBDw5feAa7wOeXI76Fqqq/g9wMsm2puk3gaeBQ8DvNm2/C3xxGcpblNnGttqPWZ/bePXlkFV/zBqvGtcaOF7fB96Z5DVJwtTP4bdYgcfLu3pmSPJmps7yYeoSwuer6t8n+W9M/Re0gGeBfzF93W61SHIVcC9wMfAMU3dR/ALwl8Bmpn5wf6eqnl+2IhdolrF9itV/zF4DnATeXFUvNm2vY5Ufs1nGtRbeY/8G+GfAGeAx4MPAL7HCjpfBL0kd46UeSeoYg1+SOsbgl6SOMfglqWMMfknqGINfkjrG4JekjjH4Jalj/h9PxBfh4dyEbAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"import matplotlib.pyplot as plt\n",
"plt.scatter(x=data.Temperature, y=data.Malfunction)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nous remarquons que des températures:\n",
" - inférieures à 65 F causent toujours un dysfonctionnement\n",
" - entre 65 et 75 F causent rarement un dysfonctionnement\n",
" - supérieures à 75 F ne causent pas de dysfonctionnement\n",
" \n",
"Essayons de transformer le problème en binaire, ç.à.d ajouter une nouvelle colonne *Malfunction2* telle que:\n",
" - Malfunction2 = 1 si Malfunction >1\n",
" - 0 autrement"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PathCollection at 0x7fce13689550>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEFBJREFUeJzt3W2MXGd5xvH/xTpRHV5qIAsijt0ElLq1xGunMZSWUqWQlw91qFo1BgpELW4kgko/RIQqaqkAlSogAUrAMmkaKDQRpalJqxTTLy1tIVHWJI3zUlMTIF6bEqdpqAqWEpu7H3Zsjde73pn1jMfz+P+TrN3zzJ0z9zPP7JWzZ87spKqQJLXlGeNuQJI0fIa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUErxnXHZ599dp133nnjuntJmkg7dux4vKqml6obW7ifd955zMzMjOvuJWkiJfluP3WelpGkBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUFLhnuSm5M8luSBRW5Pkk8k2Z3k/iSvGn6bkqRB9PMmpluAG4DPLnL7pcAF3X8bgE91v47Etnv3cv32Xex78gDnrFrJNRev4/JXrh7V3Z2QQXqdpHkNotV5jcp123Zy6917OFTFVMKmDWv44OUvHXdbQ+Fz4eRaMtyr6qtJzjtOyUbgszX3Sdt3JVmV5EVV9b0h9XjEtnv38r7bd3Lg6UMA7H3yAO+7fSfAKfckGaTXSZrXIFqd16hct20nn7vr0SPbh6qObE96wPtcOPmGcc59NbCnZ3u2OzZ012/fdeTJcdiBpw9x/fZdo7i7EzJIr5M0r0G0Oq9RufXuPQONTxKfCyffMMI9C4zVgoXJ5iQzSWb2798/8B3te/LAQOPjNEivkzSvQbQ6r1E5VAv+2Cw6Pkl8Lpx8wwj3WWBNz/a5wL6FCqtqa1V1qqozPb3kHzU7xjmrVg40Pk6D9DpJ8xpEq/MalaksdJy0+Pgk8blw8g0j3O8A3ta9aubVwA9Gcb4d4JqL17HyjKmjxlaeMcU1F68bxd2dkEF6naR5DaLVeY3Kpg1rBhqfJD4XTr4lX1BNcivweuDsJLPAHwNnAFTVFuBO4DJgN/Aj4MpRNXv4hZdJeMV9kF4naV6DaHVeo3L4RdMWr5bxuXDypcZ0Pq/T6ZR/z12SBpNkR1V1lqrzHaqS1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQX2Fe5JLkuxKsjvJtQvc/pNJ/i7Jvyd5MMmVw29VktSvJcM9yRRwI3ApsB7YlGT9vLJ3AQ9V1cuB1wMfTXLmkHuVJPWpnyP3C4HdVfVIVT0F3AZsnFdTwLOTBHgW8ARwcKidSpL61k+4rwb29GzPdsd63QD8LLAP2An8flX9eP6OkmxOMpNkZv/+/ctsWZK0lH7CPQuM1bzti4H7gHOAVwA3JHnOMf9R1daq6lRVZ3p6euBmJUn96SfcZ4E1PdvnMneE3utK4Paasxv4NvAzw2lRkjSofsL9HuCCJOd3XyS9ArhjXs2jwEUASV4IrAMeGWajkqT+rViqoKoOJrka2A5MATdX1YNJrurevgX4AHBLkp3MncZ5b1U9PsK+JUnHsWS4A1TVncCd88a29Hy/D3jjcFuTJC2X71CVpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWpQX+Ge5JIku5LsTnLtIjWvT3JfkgeT/PNw25QkDWLFUgVJpoAbgTcAs8A9Se6oqod6alYBnwQuqapHk7xgVA1LkpbWz5H7hcDuqnqkqp4CbgM2zqt5M3B7VT0KUFWPDbdNSdIg+gn31cCenu3Z7livnwaem+SfkuxI8rZhNShJGtySp2WALDBWC+zn54CLgJXA15PcVVXfPGpHyWZgM8DatWsH71aS1Jd+jtxngTU92+cC+xao+XJV/bCqHge+Crx8/o6qamtVdaqqMz09vdyeJUlL6Cfc7wEuSHJ+kjOBK4A75tV8CfilJCuSnAVsAB4ebquSpH4teVqmqg4muRrYDkwBN1fVg0mu6t6+paoeTvJl4H7gx8BNVfXAKBuXJC0uVfNPn58cnU6nZmZmxnLfkjSpkuyoqs5Sdb5DVZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDWor3BPckmSXUl2J7n2OHU/n+RQkt8YXouSpEEtGe5JpoAbgUuB9cCmJOsXqfszYPuwm5QkDaafI/cLgd1V9UhVPQXcBmxcoO7dwN8Ajw2xP0nSMvQT7quBPT3bs92xI5KsBt4EbDnejpJsTjKTZGb//v2D9ipJ6lM/4Z4Fxmre9seA91bVoePtqKq2VlWnqjrT09P99ihJGtCKPmpmgTU92+cC++bVdIDbkgCcDVyW5GBVbRtKl5KkgfQT7vcAFyQ5H9gLXAG8ubegqs4//H2SW4C/N9glaXyWDPeqOpjkauaugpkCbq6qB5Nc1b39uOfZJUknXz9H7lTVncCd88YWDPWqeseJtyVJOhG+Q1WSGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUF9hXuSS5LsSrI7ybUL3P6WJPd3/30tycuH36okqV9LhnuSKeBG4FJgPbApyfp5Zd8GfrmqXgZ8ANg67EYlSf3r58j9QmB3VT1SVU8BtwEbewuq6mtV9T/dzbuAc4fbpiRpEP2E+2pgT8/2bHdsMb8D/MOJNCVJOjEr+qjJAmO1YGHyK8yF+y8ucvtmYDPA2rVr+2xRkjSofo7cZ4E1PdvnAvvmFyV5GXATsLGq/nuhHVXV1qrqVFVnenp6Of1KkvrQT7jfA1yQ5PwkZwJXAHf0FiRZC9wO/HZVfXP4bUqSBrHkaZmqOpjkamA7MAXcXFUPJrmqe/sW4I+A5wOfTAJwsKo6o2tbknQ8qVrw9PnIdTqdmpmZGct9S9KkSrKjn4Nn36EqSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGrSin6IklwAfB6aAm6rqw/NuT/f2y4AfAe+oqm8MuVdpWbbdu5frt+9i35MHOGfVSq65eB2Xv3L1grVv+fTX+bdvPXFk+7UveR6ff+drTmifg9Ret20nt969h0NVTCVs2rCGD17+0mXM+uT0O4oefAyGI1V1/IJkCvgm8AZgFrgH2FRVD/XUXAa8m7lw3wB8vKo2HG+/nU6nZmZmTqx7aQnb7t3L+27fyYGnDx0ZW3nGFH/66y895odqfrAfNj/gB9nnILXXbdvJ5+569Jj7f+ur155QuI2q31H04GOwtCQ7qqqzVF0/p2UuBHZX1SNV9RRwG7BxXs1G4LM15y5gVZIX9d2tNCLXb9911A8TwIGnD3H99l3H1C4U7AuND7LPQWpvvXvPgve/2Hi/RtXvKHrwMRiefsJ9NdD7yM52xwatIcnmJDNJZvbv3z9or9LA9j15YKDxYe9zkNpDi/wWvdh4v0bV7yh68DEYnn7CPQuMzX+k+6mhqrZWVaeqOtPT0/30J52Qc1atHGh82PscpHYqC/0YLT7er1H1O4oefAyGp59wnwXW9GyfC+xbRo100l1z8TpWnjF11NjKM6a45uJ1x9S+9iXPW3Af88cH2ecgtZs2rDlm7Hjj/RpVv6PowcdgePq5WuYe4IIk5wN7gSuAN8+ruQO4OsltzL2g+oOq+t5QO5WW4fALVf1cofD5d76mr6tlBtnnILWHXzAc9pUio+p3FD34GAzPklfLwJGrYT7G3KWQN1fVh5JcBVBVW7qXQt4AXMLcpZBXVtVxL4XxahlJGly/V8v0dZ17Vd0J3DlvbEvP9wW8a9AmJUmj4TtUJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqUF9vYhrJHSf7ge+O5c4Xdzbw+LibGJFW5+a8Jk+rcztZ8/qpqlryj3ONLdxPRUlm+nnn1yRqdW7Oa/K0OrdTbV6elpGkBhnuktQgw/1oW8fdwAi1OjfnNXlandspNS/PuUtSgzxyl6QGndbhnuQ7SXYmuS/JTHfs/Un2dsfu6/4t+4mSZFWSLyb5jyQPJ3lNkucl+cck/9n9+txx9zmoRebVwnqt6+n/viT/m+Q9k75mx5lXC2v2B0keTPJAkluT/MSptl6n9WmZJN8BOlX1eM/Y+4H/q6qPjKuvE5XkM8C/VNVNSc4EzgL+EHiiqj6c5FrguVX13rE2OqBF5vUeJny9eiWZYu4TzzYw9xkJE71mh82b15VM8JolWQ38K7C+qg4k+QJzn3exnlNovU7rI/cWJXkO8DrgzwGq6qmqehLYCHymW/YZ4PLxdLg8x5lXay4CvlVV32XC12ye3nm1YAWwMskK5g4y9nGKrdfpHu4FfCXJjiSbe8avTnJ/kpvH/avVMrwY2A/8RZJ7k9yU5JnACw9/rm336wvG2eQyLDYvmOz1mu8K4Nbu95O+Zr165wUTvGZVtRf4CPAo8D3mPjP6K5xi63W6h/trq+pVwKXAu5K8DvgU8BLgFcwt3EfH2N9yrABeBXyqql4J/BC4drwtDcVi85r09Tqie6rp14C/Hncvw7TAvCZ6zbr/M9oInA+cAzwzyVvH29WxTutwr6p93a+PAX8LXFhV36+qQ1X1Y+DTwIXj7HEZZoHZqrq7u/1F5kLx+0leBND9+tiY+luuBefVwHr1uhT4RlV9v7s96Wt22FHzamDNfhX4dlXtr6qngduBX+AUW6/TNtyTPDPJsw9/D7wReODw4nS9CXhgHP0tV1X9F7Anybru0EXAQ8AdwNu7Y28HvjSG9pZtsXlN+nrNs4mjT11M9Jr1OGpeDazZo8Crk5yVJMw9Fx/mFFuv0/ZqmSQvZu5oHeZ+5f+rqvpQkr9k7tfFAr4D/N7h82iTIskrgJuAM4FHmLs64RnAF4C1zD05f7Oqnhhbk8uwyLw+wYSvF0CSs4A9wIur6gfdsecz+Wu20Lxa+Bn7E+C3gIPAvcDvAs/iFFqv0zbcJallp+1pGUlqmeEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KD/h/HPwqjVlmBwQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"data['Malfunction2'] = [1 if m>=1 else 0 for m in data.Malfunction]\n",
"plt.scatter(x=data.Temperature, y=data.Malfunction2)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"simpletable\">\n",
"<caption>Generalized Linear Model Regression Results</caption>\n",
"<tr>\n",
" <th>Dep. Variable:</th> <td>Malfunction2</td> <th> No. Observations: </th> <td> 23</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Model:</th> <td>GLM</td> <th> Df Residuals: </th> <td> 21</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Model Family:</th> <td>Binomial</td> <th> Df Model: </th> <td> 1</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Link Function:</th> <td>logit</td> <th> Scale: </th> <td> 1.0000</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Method:</th> <td>IRLS</td> <th> Log-Likelihood: </th> <td> -10.158</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Date:</th> <td>Wed, 26 Aug 2020</td> <th> Deviance: </th> <td> 20.315</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Time:</th> <td>16:57:48</td> <th> Pearson chi2: </th> <td> 23.2</td> \n",
"</tr>\n",
"<tr>\n",
" <th>No. Iterations:</th> <td>5</td> <th> Covariance Type: </th> <td>nonrobust</td>\n",
"</tr>\n",
"</table>\n",
"<table class=\"simpletable\">\n",
"<tr>\n",
" <td></td> <th>coef</th> <th>std err</th> <th>z</th> <th>P>|z|</th> <th>[0.025</th> <th>0.975]</th> \n",
"</tr>\n",
"<tr>\n",
" <th>Intercept</th> <td> 15.0429</td> <td> 7.379</td> <td> 2.039</td> <td> 0.041</td> <td> 0.581</td> <td> 29.505</td>\n",
"</tr>\n",
"<tr>\n",
" <th>Temperature</th> <td> -0.2322</td> <td> 0.108</td> <td> -2.145</td> <td> 0.032</td> <td> -0.444</td> <td> -0.020</td>\n",
"</tr>\n",
"</table>"
],
"text/plain": [
"<class 'statsmodels.iolib.summary.Summary'>\n",
"\"\"\"\n",
" Generalized Linear Model Regression Results \n",
"==============================================================================\n",
"Dep. Variable: Malfunction2 No. Observations: 23\n",
"Model: GLM Df Residuals: 21\n",
"Model Family: Binomial Df Model: 1\n",
"Link Function: logit Scale: 1.0000\n",
"Method: IRLS Log-Likelihood: -10.158\n",
"Date: Wed, 26 Aug 2020 Deviance: 20.315\n",
"Time: 16:57:48 Pearson chi2: 23.2\n",
"No. Iterations: 5 Covariance Type: nonrobust\n",
"===============================================================================\n",
" coef std err z P>|z| [0.025 0.975]\n",
"-------------------------------------------------------------------------------\n",
"Intercept 15.0429 7.379 2.039 0.041 0.581 29.505\n",
"Temperature -0.2322 0.108 -2.145 0.032 -0.444 -0.020\n",
"===============================================================================\n",
"\"\"\""
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import statsmodels.api as sm\n",
"data['Intercept']=1\n",
"logmodel=sm.GLM(data['Malfunction2'], data[['Intercept','Temperature']], family=sm.families.Binomial(sm.families.links.logit)).fit()\n",
"\n",
"logmodel.summary()"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl8VPW9//HXZ2ayJwRIIAIhBJBFNtlR3EAtBcWlrljFhVp+1Nrrcmur7W21Xr2tV3tb27rUYrVqBalVaxVFRSJuKAioIATZCVvY14Rs398fM2AISybJJGdm8n4+HvPIzJnvOefzZch7Ts7yPeacQ0RE4ovP6wJERCTyFO4iInFI4S4iEocU7iIicUjhLiIShxTuIiJxqNZwN7O/mlmxmS06xvtmZn8ws+Vm9oWZDYx8mSIiUhfhbLk/DYw+zvtjgG6hx0TgsYaXJSIiDVFruDvnZgPbj9PkIuAZFzQHaGlm7SJVoIiI1F0gAsvoAKyr9rooNG1jzYZmNpHg1j0pKSmDOnbsWOeV7St3bCmJ7qtq7eBPq/a82uua031mR7zns9BzC34DmwXb+ULvHXxEk6qqKny+8A/jlJRXHvO9lAR/ndd/vOXVdR1e96UhqtcSMKio9usSyVqaus91/UyiWUP6smzZsq3OuTa1tYtEuB8tYo6avs65J4AnAAYPHuzmzZtX55XtKinn1XfeZ9CgwbjQag6OoOAcOBxVDpwLvht87+A0qHIu2M6FpuGorDrY1lFZFWxTVRV8vzL0vLLKffPchV5XOSoqv3ldXllFZZWjrLKKikpHRWUV5VWhn5WOsooqDlRUUVZZRVlFJQcqqti6YzeJySmUlldxoKKS0vIq9pdVUHWM7y8HVIYeAOlJAVokB8hMTaRVagKt0hJpnZpIVnoiWelJtElPok1GEjktkmibkUxioPF+OQoKChgxYkTY7U/7zbus31lyxPQOLVP48M6z67z+Yy3vaGpbh9d9aYiud02nMvRLcVvfCn77ZfDX3G9G4a/Pi9h6mrrPdf1MollD+mJma8JpF4lwLwKqb4LnAhsisNyjykxJoGOGj17tWzTWKprU0T5k54JfECVllewvq2R/WQV7SivYd6CSvQfK2VNaceixq6Q89Chjx/5ylmzYzbZ9ZewqKT/q+rLTk+jQMpn2LVPo2DqVjq1SyMtKo3NWGh1apeBvwj8H7vh2D+566cvDtgBTEvzc8e0eEVtegi/4p1B55Tfflg1ZR13W3RjrCcdVwzry3Jy1R50eSdHUZzlSJML9VeBmM5sKDAN2OeeO2CUj4TMzkgJ+kgJ+WqbWbxnllVXs2FdG8Z4DbNlzgM27S9m0u5RNu0pZv7OEws17mLm0mLKKqkPzJPiNztlpdGubQbecdE5q14Je7VqQ2yoFs8iH/sUDOgDw4IxCNuwsoX3LFO74do9D0yO1vEiuo67rjvR6wnHfxX0BmPJJcG+p34yrhnU8ND1SoqnPcqRaw93MpgAjgGwzKwLuBhIAnHOPA9OB84DlwH7ghsYqVsKX4PfRtkUybVskH7NNVZVjy94DrN66j9Xb9rFy6z5WFO9l0YZdTF+08dDursyUBPrlZtIvN5NBnVoxKK81makJEanz4gEdIhoGx1peUwROpPvSEPdd3Jf7Lu5LQUEBK64e0WjriaY+y+FqDXfn3FW1vO+AH0asImkyPp+R0yKZnBbJDOuSddh7+8sqKNy0h6827mbR+l18UbSLx99bSWXoYEDPEzI4tWsWw7tmM7xrFmlJkfgjULxUXl5OUVERpaWlXpdyVJmZmSxZssTrMiIinL4kJyeTm5tLQkL9NqT0GylHlZoYYEBeKwbktTo0bX9ZBZ+v28Vna7YzZ+V2nv9kLU99uJoEvzEkvzWdEso5ccd+clvVc1+SeKqoqIiMjAzy8/MbZTdcQ+3Zs4eMjAyvy4iI2vrinGPbtm0UFRXRuXPneq1D4S5hS00McGrXLE7tmsXNZ8OBiko+W7OD9wq3MKuwmI82lzFl6Sz6dGjBBf3ac2H/9rTLTPG6bAlTaWlp1AZ7c2NmZGVlsWXLlnovQ+Eu9ZYU8Id2y2Rz13knMW36u+xM78TrX27i128s5TdvLmV41yzGDcljVO8ckgJNe7631J2CPXo09LNQuEvEtE31ccWZXZl4ZldWbd3Hvxau5x/zivjRlAW0TkvkqqEdGX9KPidkHvsgr4hERnxc7iVRp3N2Gree2533fzKSZyYMZXCnVjxasILTH3iX219YyPLivV6XKFHIzBg/fvyh1xUVFbRp04axY8ced76CgoJDbQ4cOMC5555L//79eeGFFyJW2yuvvMJXX3116PUvf/lL3nnnnXot64UXXqBfv37069eP4cOH8/nnn0eqzEO05S6NyuczzuzehjO7t2Httv089dEqpn66jpcXrue8vu247dzunNg23esyJUqkpaWxaNEiSkpKSElJ4e2336ZDh7qdarlgwQLKy8tZuHBhRGt75ZVXGDt2LL169QLg3nvvrfey8vPzee+992jVqhVvvPEGEydO5JNPPolUqYC23KUJ5WWlcvcFvfngpyO5aURX3ivcwqjfvcddL33B5t3RefqdNL0xY8bw+uuvAzBlyhSuuuqbs7E//fRThg8fzoABAxg+fDiFhYWHzVtcXMw111zDwoUL6d+/PytWrCA/P5+tW7cCMG/evENXhN9zzz1MmDCBESNG0KVLF/7whz8cWs4zzzxDv379OPnkkxk/fjwfffQRr776Knfccceh5V5//fW8+OKLAMycOZMBAwbQt29fJkyYwIEDB4BgiN99990MHDiQvn37snTpUgCGDRtGq1bBM9FOOeUUioqKIv7vqC13aXJZ6Unc8e2eTDitM3+atZzn5qzhXws38KOzu/G90zs36vg3Ep5f/XsxX23YHdFl9mrfgrsv6F1ru3HjxnHvvfcyduxYvvjiCyZMmMD7778PQM+ePZk9ezaBQIB33nmHn/3sZ/zzn/88NG/btm2ZPHkyDz30EK+99lqt61q6dCmzZs1iz5499OjRgx/84AcsW7aM+++/nw8//JDs7Gy2b99O69atufDCCxk7diyXXXbZYcsoLS3l+uuvZ+bMmXTv3p1rr72Wxx57jFtvvRWA7Oxs5s+fz6OPPspDDz3E5MmTD5v/ySefZMyYMbXWWlf6LRLPZKUncfcFvZl5+wiGd83mgTeXMvr3s5mzcpvXpYmH+vXrx+rVq5kyZQrnnXf4QGe7du3i8ssvp0+fPtx2220sXry4Qes6//zzSUpKIjs7m7Zt27J582beffddLrvsMrKzswFo3br1cZdRWFhI586d6d69OwDXXXcds2fPPvT+JZdcAsCgQYNYvXr1YfPOmjWLJ598kgceeKBB/TgabbmL5/KyUpl83WBmLS3ml68uYtwTc7j21E78dHRPXfnqkXC2sBvThRdeyI9//GMKCgrYtu2bL/tf/OIXjBw5kpdffpnVq1eHNbJiIBCgqio4hlLNq2+TkpIOPff7/VRUVOCcq9NpiM4dfwjyg+s4uPyDvvjiC2688UbeeOMNsrKyjjV7vWnLXaLGyJ5tmXHrmdxwWj7PzlnD6Idns2DtDq/LEg9MmDCBX/7yl/Tte/hgZ7t27Tp0gPXpp58Oa1n5+fl89tlnAIftwjmWc845h2nTph36Utm+PXivooyMDPbs2XNE+549e7J69WqWL18OwLPPPstZZ5113HWsW7eOSy65hGefffbQFn+kKdwlqqQmBrj7gt68MPFUqqrg8sc/5rGCFVQda4B7iUu5ubnccsstR0z/yU9+wl133cVpp51GZWV4N2e5++67ueWWWzjjjDPw+2u/kK537978/Oc/56yzzuLkk0/m9ttvB4LHAh588EEGDBjAihUrDrVPTk7mqaee4vLLL6dv3774fD4mTZp03HU88MADbNu2jZtuuon+/fszePDgsPpSJ845Tx6DBg1y9TVr1qx6zxtt1Jdj27m/zP3guXmu009fczc89anbVVIW0eUfT7x8LnXpx1dffdV4hUTA7t27vS4hYsLty9E+E2CeCyNjteUuUSszJYFHvjuQey/qzexlW7j4Tx+yvPjIP4tF5EgKd4lqZsa1p+bz9xuHsbu0nO888hEfLt/qdVkiUU/hLjFhWJcsXr35dNq3TOH6pz7l5QWRv+hDaj/zQ5pOQz8LhbvEjPYtU5g26VQGd2rNbS98zl9mr/S6pLiSnJzMtm3bFPBRwIXGc09Orv8gezqJWGJKZkoCT08Ywu3TPuf+6UsoLa/kR+d087qsuJCbm0tRUVGDxhBvTKWlpQ0Ku2gSTl8O3ompvhTuEnOSAn7+MG4ASQEfv317GaUVlfx4VA+NRd5ACQkJ9b7rT1MoKChgwIABXpcREU3RF4W7xCS/z3jospNJCvh5ZNYKEv1+bjlXW/AiByncJWb5fMb9F/ehvLKK372zjLQkPzee0cXrskSigsJdYprPZ/zmkr7sO1DBfa8voUVKAlcM7uh1WSKe09kyEvMCfh8PjxvAGd2yueulL5m9LDoPCIo0JYW7xIXEgI9Hrx5It7bp3PT3+SzdFNmxyEVijcJd4kZGcgJ/vX4IaUl+bnhqLsV7dHcnab4U7hJX2rdM4cnrhrBjfxk//Pt8yiqqvC5JxBMKd4k7fTpk8sCl/Zi7egf3v/5V7TOIxCGdLSNx6aL+HVi0fhd/eX8VfXNbctmg+l/pJxKLtOUuceuno3syvGsW//XKlxoqWJodhbvErYDfx+/H9Sc9KcAP/76A0vLw7twjEg8U7hLX2mYk89sr+lO4eQ/3vqb979J8KNwl7p3VvQ2TzurK85+sZcbiTV6XI9IkFO7SLPznqO70bt+Cn7/8Jdv2HvC6HJFGp3CXZiHB7+P/rujP7pIKfv7yIt2QQuKewl2ajR4nZHD7qO68uXgT/1q4wetyRBpVWOFuZqPNrNDMlpvZnUd5P9PM/m1mn5vZYjO7IfKlijTc98/owsC8ltzz78Vs1e4ZiWO1hruZ+YFHgDFAL+AqM+tVo9kPga+ccycDI4DfmllihGsVaTC/z3jg0n7sO1DBf+vsGYlj4Wy5DwWWO+dWOufKgKnARTXaOCDDgvc5Swe2AxURrVQkQrrlZHDTiBP518INFBQWe12OSKOw2g4smdllwGjn3I2h1+OBYc65m6u1yQBeBXoCGcCVzrnXj7KsicBEgJycnEFTp06tV9F79+4lPT29XvNGG/XFG+VVjl98WEJ5JfzP6SkkBQ6//2os9eV44qUfoL4cNHLkyM+cc4NrbeicO+4DuByYXO31eOCPNdpcBvwOMOBEYBXQ4njLHTRokKuvWbNm1XveaKO+eOeTldtcp5++5h54Y8kR78VaX44lXvrhnPpyEDDP1ZLbzrmwdssUAdXvW5YL1DzV4AbgpdC6l4fCvWcYyxbxzNDOrblkQAf+8v5KVm3d53U5IhEVTrjPBbqZWefQQdJxBHfBVLcWOAfAzHKAHsDKSBYq0hjuHNOTpICfX/17sc59l7hSa7g75yqAm4EZwBJgmnNusZlNMrNJoWb/DQw3sy+BmcBPnXNbG6tokUhp2yKZW8/tRkHhFmYu0cFViR9hjefunJsOTK8x7fFqzzcAoyJbmkjTuG54Pi/MXce9r33FGd2zSQr4vS5JpMF0hao0ewl+Hz8//yTWbt/Psx+v8bockYhQuIsAI3q05Yxu2fzx3eXs3F/mdTkiDaZwFwn52Xknsbu0nD+9u9zrUkQaTOEuEnJSuxZcMagjf/t4NcX7q7wuR6RBFO4i1dw+qjt+n/Hy19o1I7FN4S5STU6LZK4f3pk5GytZumm31+WI1JvCXaSGSWd1ITkAv31rmdeliNSbwl2khpapiYzpnMDbX21m/todXpcjUi8Kd5GjGNUpgez0RB6aUeh1KSL1onAXOYrkgPGDESfy0YptfLJym9fliNSZwl3kGK4elkd2ehIPz/za61JE6kzhLnIMyQl+Jp3VhY9WbGPu6u1elyNSJwp3keO4elgnstMTefgdbb1LbFG4ixxHSqKf/3dmVz5YvpV52nqXGKJwF6nF1afk0TotkUdmacwZiR0Kd5FapCYGmHBaPrMKt7B4wy6vyxEJi8JdJAzjT80nPSnAYwUrvC5FJCwKd5EwZKYkcM0pnZj+5UbdTFtigsJdJEwTTs8n4Pfx5/e09S7RT+EuEqa2GclcMTiXf84vYvPuUq/LETkuhbtIHXz/jC5UVDn+9tFqr0sROS6Fu0gddMpKY3TvE3huzhr2HqjwuhyRY1K4i9TR98/swu7SCqbNXed1KSLHpHAXqaOBea0Ykt+KJz9YRUWl7rUq0UnhLlIP3z+jC+t3ljB90SavSxE5KoW7SD2ce1IOnbPT+OsHq7wuReSoFO4i9eDzGTecls/CdTt1Kz6JSgp3kXq6dGAuGckBbb1LVFK4i9RTWlKAq4bm8caiTazfWeJ1OSKHUbiLNMC1p3bCOcczH6/2uhSRwyjcRRogt1UqY/q0Y8ona9lfpouaJHoo3EUa6PrT8tldWsErCzZ4XYrIIQp3kQYa3KkVvdq14G8frcY553U5IoDCXaTBzIzrh+dTuHkPc1bqPqsSHcIKdzMbbWaFZrbczO48RpsRZrbQzBab2XuRLVMkul3Yvz0tUxM0WqREjVrD3cz8wCPAGKAXcJWZ9arRpiXwKHChc643cHkj1CoStZIT/Fw5pCNvfaXTIiU6hLPlPhRY7pxb6ZwrA6YCF9Vo813gJefcWgDnXHFkyxSJfuNP6QTA3+es8bgSEbDaDgCZ2WXAaOfcjaHX44Fhzrmbq7X5PZAA9AYygIedc88cZVkTgYkAOTk5g6ZOnVqvovfu3Ut6enq95o026kt0qm9fHp5fyvKdlfzfiFQSfNYIldWNPpPo1JC+jBw58jPn3ODa2gXCWNbR/ofW/EYIAIOAc4AU4GMzm+OcW3bYTM49ATwBMHjwYDdixIgwVn+kgoIC6jtvtFFfolN9++Jrv4Vr//op+1t356L+HSJfWB3pM4lOTdGXcHbLFAEdq73OBWqe0FsEvOmc2+ec2wrMBk6OTIkiseP0E7PpnJ3GMx9r14x4K5xwnwt0M7POZpYIjANerdHmX8AZZhYws1RgGLAksqWKRD+fz7h6WB6frdnB4g27vC5HmrFaw905VwHcDMwgGNjTnHOLzWySmU0KtVkCvAl8AXwKTHbOLWq8skWi1+WDOpKc4OO5OWu9LkWasXD2ueOcmw5MrzHt8RqvHwQejFxpIrEpMzWBC09uzysL1nPnmJ5kpiR4XZI0Q7pCVaQRjD8ln5LySl5ZsN7rUqSZUriLNIK+uZn0y83kuTlrNN6MeELhLtJIrhnWia+L9zJ3tW7DJ01P4S7SSMae3I6M5ADP6YpV8YDCXaSRpCYGuHRgLm8s2sjWvQe8LkeaGYW7SCO6elge5ZWOf8wr8roUaWYU7iKNqFtOBkM7t2bKp2upqtKBVWk6CneRRnb1sDzWbt/PB8u3el2KNCMKd5FGNrrPCbROS+T5T3TFqjQdhbtII0sK+LlsUC5vL9nM5t2lXpcjzYTCXaQJXDU0j8oqx7S567wuRZoJhbtIE+icncZpJ2Yxde46KnVgVZqAwl2kiXx3aCfW7yxh9rItXpcizYDCXaSJfKtXDtnpifxdB1alCSjcRZpIYsDH5YM78u7SzWzcVeJ1ORLnFO4iTeiqIXlUOXhBB1alkSncRZpQXlYqZ3TL5oW566iorPK6HIljCneRJnb1sDw27iqloFAHVqXxKNxFmtg5J+XQNiOJ5z/VgVVpPAp3kSaW4Pdx5ZCOzCospmjHfq/LkTilcBfxwJVDOgI6sCqNR+Eu4oHcVqmM6N6GF+auo1wHVqURKNxFPPLdYZ0o3nOAmUuKvS5F4pDCXcQjI3u04YQWyTqwKo1C4S7ikYDfx7ihHXn/6y2s3aYDqxJZCncRD105pCM+M6bM1da7RJbCXcRD7TJTOKdnW6bNXUdZhQ6sSuQo3EU8dvUpndi2r4wZizd5XYrEEYW7iMfOODGbjq1TeG7OGq9LkTiicBfxmM9nfHdoJz5ZtZ3lxXu8LkfihMJdJApcMTiXBL/x3BwdWJXIULiLRIGs9CTG9GnHP+cXsb+swutyJA4o3EWixPhTO7GntIJ/f77B61IkDijcRaLE4E6t6JGTwTMfr8E553U5EuMU7iJRwsy45tROLN6wm4XrdnpdjsS4sMLdzEabWaGZLTezO4/TboiZVZrZZZErUaT5+M6ADqQl+nlWp0VKA9Ua7mbmBx4BxgC9gKvMrNcx2j0AzIh0kSLNRXpSgEsG5vLaFxvZvq/M63IkhoWz5T4UWO6cW+mcKwOmAhcdpd2PgH8CGr9UpAGuOaUTZRVVupGHNEggjDYdgOr/y4qAYdUbmFkH4DvA2cCQYy3IzCYCEwFycnIoKCioY7lBe/furfe80UZ9iU5e96VHKx+TCwrp4dbiM6v3crzuRySpL3UTTrgf7X9WzUP5vwd+6pyrtOP8R3TOPQE8ATB48GA3YsSIMMs8XEFBAfWdN9qoL9HJ676UZG3kB3+fT2VOL87ulVPv5Xjdj0hSX+omnN0yRUDHaq9zgZon4g4GpprZauAy4FEzuzgiFYo0Q9/qlUO7zGT+9tFqr0uRGBVOuM8FuplZZzNLBMYBr1Zv4Jzr7JzLd87lAy8CNznnXol4tSLNRMDv4+pheXywfKvGm5F6qTXcnXMVwM0Ez4JZAkxzzi02s0lmNqmxCxRprsYNzSPR7+PZj3VapNRdOPvccc5NB6bXmPb4Mdpe3/CyRCQ7PYmx/drx4mdF/Oe3e9AiOcHrkiSG6ApVkSh2w2md2VdWyT/mFXldisQYhbtIFOubm8mQ/FY8/dEqKqs03oyET+EuEuVuOK0z67aX8M6SzV6XIjFE4S4S5Ub1yqFDyxSe+nCV16VIDFG4i0S5gN/HdcM7MWfldhZv2OV1ORIjFO4iMeDKwXmkJvp58gNtvUt4FO4iMSAzNYErBnfk1YUb2LSr1OtyJAYo3EVixPdO70yVczytIQkkDAp3kRjRsXUqY/q04++frGHvAd1EW45P4S4SQ248ozN7SiuYprHepRYKd5EYMiCvFYM7teLJD1ZRXlnldTkSxRTuIjHm/53VlfU7S3j9i41elyJRTOEuEmPO6dmW7jnpPFawgioNSSDHoHAXiTE+nzHprK4Ubt7DrELdsliOTuEuEoMuOLk9HVqm8FjBCq9LkSilcBeJQQl+HxPP7MK8NTv4dNV2r8uRKKRwF4lRVwzuSFZaIn9892uvS5EopHAXiVEpiX4mntmF97/eyvy1O7wuR6KMwl0khl1zSidapSbwx5naepfDKdxFYlhaUoAbz+jCrMItfFG00+tyJIoo3EVi3LWndiIzJYE/aOtdqlG4i8S4jOQEbjy9M+8sKdbWuxyicBeJA9eflk+r1AQeemuZ16VIlFC4i8SBjOQEfjCiK7OXbeGTldu8LkeigMJdJE5ce2o+bTOSeOitQpzTmDPNncJdJE4kJ/j50dknMnf1Dt5btsXrcsRjCneROHLlkDw6tk7hgTcLqdSIkc2awl0kjiQGfPx4VA+WbNzNywvWe12OeEjhLhJnLujXnn65mfz2rULKKrX13lwp3EXijM9n/Oy8k9i4q5S3Vpd7XY54ROEuEodO6ZLFuSe15bWV5Wzde8DrcsQDCneROHXXeSdRXgUPvlnodSniAYW7SJzq2iadb3UKMO2zdRqWoBlSuIvEsYtOTCQrLYl7Xl2sm2k3M2GFu5mNNrNCM1tuZnce5f2rzeyL0OMjMzs58qWKSF2lBIw7x/Rk/tqdOjWymak13M3MDzwCjAF6AVeZWa8azVYBZznn+gH/DTwR6UJFpH4uGdCBAXkt+Z/pS9i5v8zrcqSJhLPlPhRY7pxb6ZwrA6YCF1Vv4Jz7yDl38D5fc4DcyJYpIvXl8xn/852+7Cwp59fTl3pdjjQRq22AITO7DBjtnLsx9Ho8MMw5d/Mx2v8Y6HmwfY33JgITAXJycgZNnTq1XkXv3buX9PT0es0bbdSX6BQvfanej2mFZUxfVc5dQ5Pp0drvcWV1Fy+fCTSsLyNHjvzMOTe41obOueM+gMuBydVejwf+eIy2I4ElQFZtyx00aJCrr1mzZtV73mijvkSneOlL9X7sO1DuTvvNTHf2Q7NcaXmFd0XVU7x8Js41rC/APFdLvjrnwtotUwR0rPY6F9hQs5GZ9QMmAxc55zSgtEiUSU0McN/FfVixZZ9uydcMhBPuc4FuZtbZzBKBccCr1RuYWR7wEjDeOadbwYhEqRE92nLF4FweK1jBwnU69z2e1RruzrkK4GZgBsFdLtOcc4vNbJKZTQo1+yWQBTxqZgvNbF6jVSwiDfJfY3uR0yKZ/5y2kNLySq/LkUYS1nnuzrnpzrnuzrmuzrn7Q9Med849Hnp+o3OulXOuf+hR+85+EfFEi+QEHri0Hyu27OOhGRqaIF7pClWRZujM7m245pQ8Jn+witm6a1NcUriLNFP/dX4vuuekc/u0z9myRyNHxhuFu0gzlZzg549XDWRPaTn/+Y/PNfZMnFG4izRjPU7I4BdjezF72RYee2+F1+VIBCncRZq5q4flcVH/9jz0ViHvaf973FC4izRzZsavL+lLj5wMbpm6gHXb93tdkkSAwl1ESE0M8Pg1g6iscnz/mXnsPVDhdUnSQAp3EQEgPzuNP313IF8X7+U/piygUgdYY5rCXUQOOat7G+65sDfvLi3mvte/8rocaYCA1wWISHQZf0onVm/dx5MfrKJ9ZgrfP7OL1yVJPSjcReQIPzvvJDbtKuX+6UvITEngiiEda59JoorCXUSO4PcZv7uyP7tLy7nzpS/ISA4wpm87r8uSOtA+dxE5qsSAjz+PH8SAvFb8aMoC3ly00euSpA4U7iJyTKmJAZ66YQj9cjP54fMLeP0LBXysULiLyHG1SE7gme8NY2BeS340ZT4vflbkdUkSBoW7iNQqPSnA0zcMZXjXbH78j895tGD5wfsmS5RSuItIWNKSAvz1+iFceHJ7/vfNQu55dTEVlVVelyXHoLNlRCRsiQEfv7+yPzkmf/gSAAAM7ElEQVQtkvjL+6tYuXUff7xqAC1TE70uTWrQlruI1InPZ/z8/F7876X9mLNyGxc/8iGFm/Z4XZbUoHAXkXq5YkhHpk48hX1llVz4pw+Y8ula7YePIgp3Eam3QZ1aM/0/zmBIfmvueulLfjRlATv3l3ldlqBwF5EGapORxDMThnLHt3vw5qJNfOt3s5m5ZLPXZTV7CncRaTCfz/jhyBN55YenkZWWyPf+No9bpi6geHep16U1Wwp3EYmYPh0yefXm07nlnG688eUmzvnte/z1g1WU65TJJqdwF5GISgz4uO1b3Zlx25n0z2vJva99xajfzebNRRt1wLUJKdxFpFF0zk7jmQlDefK6wfh9xqTn5nPxox8xa2mxQr4J6CImEWk0ZsY5J+VwVvc2vPhZEX+atZwbnp5Lv9xMJp7ZhdG9TyDg1zZmY1C4i0ijC/h9jBuax6WDcnlpfhGPFazg5ucX0KFlCtee2onLBuWSlZ7kdZlxReEuIk0mwe/jyiF5XD6oIzOXFjP5/ZX8+o2l/PatZYzqncOlg3I548Rsbc1HgMJdRJqcz2d8q1cO3+qVw9eb9/D8p2t5af56XvtiI9npiZzftx2j+7RjaOfW+H3mdbkxSeEuIp7qlpPB3Rf05s4xPXmvcAsvL1jP1Lnr+NvHa8hKS2REj7aM7NkGK9dB2LpQuItIVEgK+BnV+wRG9T6BfQcqKCjcwozFm3hnyWb+Ob8IA/5c+AHDu2YzrHNrBnZqRWZKgtdlRy2Fu4hEnbSkAOf3a8f5/dpRUVnFwnU7efbteWyo9DH5/ZU8/t4KzKBb23T65bbk5I4t6dO+BT1PaEFKot/r8qOCwl1EolrA72Nwfmv2dktkxIjh7C+rYOG6ncxdtYOF63bw7tLiQ7f+8xnkZ6XRLSed7jkZdG2TTufsNDq3SaNFcvPayle4i0hMSU0MMLxrNsO7ZgPgnKNoRwlLNu7mq427WbpxD8uK9/DOkmIqq77ZT98qNYG81ql0aJVCh5YptG+ZwgktkjkhM5mcFslkpyeRGIifs3TCCnczGw08DPiByc6539R430LvnwfsB653zs2PcK0iceuVBet5cEYhG3aW0L5lCnd8uwf/mLeWD1dsP9TmtK6tuXxw3hHtgCOmzVuznSmfrOPWPuV8767pXDWsI/dd3Des9V48oMMxp4cz/8F1VzqH36xO6z5aX8Jd78UDOjCq9wmH2hyoqGTttv2s3LqPVVv3sW77ftZu38/SjXuYuaSYAxVHjneTmZJAVnoi2WlJtE5LpFVaAi1TE2mZkkBm6JGRnEBGcoCM5ADpyQHSkwKkJPgJxmD0qDXczcwPPAJ8CygC5prZq865r6o1GwN0Cz2GAY+FfopILV5ZsJ67XvqSkvJKANbvLOHWFxYe0e7DFdsPC/v1O0u448XPwUF5aAt1/c4Sbn9hIdVjq9I5npuzFuCwkD3aeu966UvmrdnOPz9bf8R04LCgPdr8DVn3Hf/4HAzKK7/pS7jrPVq7pICfbjkZdMvJOOLf0jnHtn1lbNpVyubdpRTvOcCW0GPbvgNs3VvGii172bGmnJ37y6ioOv6ZOmaQmuAnJTFAaqKflAQ/yYl+UhJ8JCf4SQp88zMx4KP1gQpGHHeJDRfOlvtQYLlzbiWAmU0FLgKqh/tFwDMuOGDEHDNraWbtnHMbI16xSJx5cEbhoaCqq4NBWN2xxl+c8sm6wwL2aOstKa88tNVdc/qDMwoPC8+jzd+QdZcfJUDDXe/R2h2PmZGdnkR2ehJ9OmQet61zjv1llewsKWd3STl7SivYXVLOvrIK9pRWsPdABfvLKtkX+llSFvxZWlFFaVkl2/eVUVpeyYGKKsoqqjhQUcVZ7Rr/tE6rbQAfM7sMGO2cuzH0ejwwzDl3c7U2rwG/cc59EHo9E/ipc25ejWVNBCaGXvYACutZdzawtZ7zRhv1JTo1WV8STzhxUGMtu3L/Lvyp34RX2ablnzVkvQ2ZPwLzZgNbjzdv9XVEuYb8/+rknGtTW6NwttyPtiOp5jdCOG1wzj0BPBHGOo9fkNk859zghi4nGqgv0Sle+mJm8yp2Fcd8PyB+PhNomr6Ec2i4COhY7XUusKEebUREpImEE+5zgW5m1tnMEoFxwKs12rwKXGtBpwC7tL9dRMQ7te6Wcc5VmNnNwAyCp0L+1Tm32Mwmhd5/HJhO8DTI5QRPhbyh8UoGIrBrJ4qoL9EpXvoSL/0A9aVOaj2gKiIisSd+LscSEZFDFO4iInEo6sPdzJLN7FMz+9zMFpvZr0LTW5vZ22b2dehnK69rDYeZ+c1sQejagFjux2oz+9LMFprZvNC0WO1LSzN70cyWmtkSMzs1FvtiZj1Cn8fBx24zuzVG+3Jb6Pd9kZlNCeVAzPUDwMxuCfVjsZndGprW6H2J+nAHDgBnO+dOBvoDo0Nn5NwJzHTOdQNmhl7HgluAJdVex2o/AEY65/pXO183VvvyMPCmc64ncDLBzyfm+uKcKwx9Hv2BQQRPbniZGOuLmXUA/gMY7JzrQ/BEjnHEWD8AzKwP8H2CV/qfDIw1s240RV+cczHzAFKB+QTHrSkE2oWmtwMKva4vjPpzQx/k2cBroWkx149QrauB7BrTYq4vQAtgFaGTC2K5LzXqHwV8GIt9AToA64DWBM/oey3Un5jqR6jOywkOtnjw9S+AnzRFX2Jhy/3groyFQDHwtnPuEyDHhc6lD/1s62WNYfo9wQ+2+hAcsdgPCF6B/JaZfRYaVgJisy9dgC3AU6HdZZPNLI3Y7Et144Apoecx1Rfn3HrgIWAtsJHgdTNvEWP9CFkEnGlmWWaWSvCU8Y40QV9iItydc5Uu+KdmLjA09KdOTDGzsUCxcy5Wxr6ozWnOuYEERwT9oZmd6XVB9RQABgKPOecGAPuIgT/3jyd0seGFwD+8rqU+QvufLwI6A+2BNDO7xtuq6sc5twR4AHgbeBP4HKhoinXHRLgf5JzbCRQAo4HNZtYOIPSz2MPSwnEacKGZrQamAmeb2XPEXj8AcM5tCP0sJrhfdyix2ZcioCj01yDAiwTDPhb7ctAYYL5zbnPodaz15VxglXNui3OuHHgJGE7s9QMA59yTzrmBzrkzge3A1zRBX6I+3M2sjZm1DD1PIfjBLyU45MF1oWbXAf/ypsLwOOfucs7lOufyCf7J/K5z7hpirB8AZpZmZhkHnxPcH7qIGOyLc24TsM7MeoQmnUNwOOuY60s1V/HNLhmIvb6sBU4xs1QzM4KfyRJirx8AmFnb0M884BKCn02j9yXqr1A1s37A3wgeMfcB05xz95pZFjANyCP4n+Fy59z2Yy8pepjZCODHzrmxsdgPM+tCcGsdgrs1nnfO3R+LfQEws/7AZCARWElw+AwfsdmXVIIHI7s453aFpsXc5xI65flKgrswFgA3AunEWD8AzOx9IAsoB253zs1sis8k6sNdRETqLup3y4iISN0p3EVE4pDCXUQkDincRUTikMJdRCQOhXODbJEmFTpNbGbo5QlAJcEhAgCGOufKPCnsOMxsAjA9dN68iOd0KqRENTO7B9jrnHsoCmrxO+cqj/HeB8DNzrmFdVhewDnXJJeiS/Oj3TISU8zsOguO77/QzB41M5+ZBcxsp5k9aGbzzWyGmQ0zs/fMbKWZnRea90Yzezn0fqGZ/VeYy73PzD4lOK7Rr8xsbmh87sct6EqCw1G/EJo/0cyKql1ZfYqZvRN6fp+Z/dnM3iY4WFnAzP4vtO4vzOzGpv9XlXikcJeYERow7jvA8NBAcgGCQzkAZAJvhQYzKwPuIXjZ+uXAvdUWMzQ0z0Dgu2bWP4zlznfODXXOfQw87JwbAvQNvTfaOfcCsBC40gXHU69tt9EA4ALn3HhgIsEB5YYCQwgOwpZXn38fkeq0z11iybkEA3BecMgRUgheag9Q4px7O/T8S4LDxFaY2ZdAfrVlzHDO7QAws1eA0wn+HhxruWV8M9QCwDlmdgeQDGQDnwFv1LEf/3LOlYaejwJOMrPqXybdCF6SLlJvCneJJQb81Tn3i8MmmgUIhvBBVQTv4HXwefX/5zUPMrlallviQgemQuO2/AkY6Jxbb2b3EQz5o6ngm7+Ma7bZV6NPNznnZiISQdotI7HkHeAKM8uG4Fk19diFMcqC90xNJThm+Id1WG4KwS+LraFRMS+t9t4eIKPa69UEb3VHjXY1zQBuCn2RHLwPakod+yRyBG25S8xwzn0ZGi3wHTPzERxlbxKwoQ6L+QB4HugKPHvw7JZwluuc22ZmfyM4vPEa4JNqbz8FTDazEoL79e8B/mJmm4BPj1PPnwmODLgwtEuomOCXjkiD6FRIaTZCZ6L0cc7d6nUtIo1Nu2VEROKQttxFROKQttxFROKQwl1EJA4p3EVE4pDCXUQkDincRUTi0P8H1Tp8BrA8PtUAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"import numpy as np\n",
"data_pred = pd.DataFrame({'Temperature': np.linspace(start=30, stop=90, num=121), 'Intercept': 1})\n",
"data_pred['Malfunction2'] = logmodel.predict(data_pred[['Intercept','Temperature']])\n",
"data_pred.plot(x=\"Temperature\",y=\"Malfunction2\",kind=\"line\",ylim=[0,1])\n",
"plt.scatter(x=data[\"Temperature\"],y=data[\"Malfunction2\"])\n",
"plt.grid(True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Conclusion\n",
"Avec cette nouvelle perspective, on peut déduire qu'à une température de 31 F, un dysfonctionnement est extrèmement probable."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
} }
], ],
"metadata": { "metadata": {
...@@ -705,7 +1027,7 @@ ...@@ -705,7 +1027,7 @@
"name": "python", "name": "python",
"nbconvert_exporter": "python", "nbconvert_exporter": "python",
"pygments_lexer": "ipython3", "pygments_lexer": "ipython3",
"version": "3.7.3" "version": "3.6.4"
} }
}, },
"nbformat": 4, "nbformat": 4,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment