{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Расчет площадей лесных пожаров" ] }, { "cell_type": "code", "execution_count": 119, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import math\n", "import numpy as np\n", "import pandas as pd\n", "import random\n", "import sklearn\n", "\n", "# importing sklearn libraries\n", "from sklearn import neural_network, linear_model, preprocessing, svm, tree\n", "from sklearn.ensemble import RandomForestRegressor, RandomForestClassifier\n", "from sklearn.metrics import accuracy_score, mean_squared_error, r2_score\n", "from sklearn.model_selection import train_test_split\n", "from sklearn.linear_model import LassoLarsCV\n", "from sklearn.linear_model import Lasso\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Загрузка данных" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "В данные входят следующие величины:\n", "\n", " X - Пространственная координата оси X на карте парка Монтезиньо: от 1 до 9\n", " Y - Пространственная координата оси X на карте парка Монтезиньо: от 2 до 9\n", " month - название месяца: \"jan\" до \"dec\"\n", " day - день недели: \"mon\" до \"sun\"\n", " FFMC - FFMC-индекс FWI-системы: от 18.7 до 96.20\n", " DMC - DMC-индекс FWI-системы: от 1.1 до 291.3\n", " DC - DC-индекс FWI-системы: от 7.9 до 860.6\n", " ISI - ISI-индекс FWI-системы: от 0.0 до 56.10\n", " temp - температура воздуха в Цельсиях: от 2.2 до 33.30\n", " RH - относительная влажность в%: от 15.0 до 100\n", " wind - скорочть ветра в км/ч: от 0.40 до 9.40\n", " rain - количество осадков в мм/м2: от 0.0 до 6.4\n", " area - размер территории леса, охваченной огнем(в гектарах): от 0.00 до 1090.84\n" ] }, { "cell_type": "code", "execution_count": 120, "metadata": {}, "outputs": [], "source": [ "forest_fires = pd.read_csv('forestfires.csv')" ] }, { "cell_type": "code", "execution_count": 121, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
XYmonthdayFFMCDMCDCISItempRHwindrainarea
075marfri86.226.294.35.18.2516.70.00.0
174octtue90.635.4669.16.718.0330.90.00.0
274octsat90.643.7686.96.714.6331.30.00.0
386marfri91.733.377.59.08.3974.00.20.0
486marsun89.351.3102.29.611.4991.80.00.0
\n", "
" ], "text/plain": [ " X Y month day FFMC DMC DC ISI temp RH wind rain area\n", "0 7 5 mar fri 86.2 26.2 94.3 5.1 8.2 51 6.7 0.0 0.0\n", "1 7 4 oct tue 90.6 35.4 669.1 6.7 18.0 33 0.9 0.0 0.0\n", "2 7 4 oct sat 90.6 43.7 686.9 6.7 14.6 33 1.3 0.0 0.0\n", "3 8 6 mar fri 91.7 33.3 77.5 9.0 8.3 97 4.0 0.2 0.0\n", "4 8 6 mar sun 89.3 51.3 102.2 9.6 11.4 99 1.8 0.0 0.0" ] }, "execution_count": 121, "metadata": {}, "output_type": "execute_result" } ], "source": [ "forest_fires.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Статистический анализ данных" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Рассмотрим значение стандартных статистических величин у каждого из данных признаков, а так же попарную корелляцию" ] }, { "cell_type": "code", "execution_count": 122, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
XYFFMCDMCDCISItempRHwindrainarea
count517.000000517.000000517.000000517.000000517.000000517.000000517.000000517.000000517.000000517.000000517.000000
mean4.6692464.29980790.644681110.872340547.9400399.02166318.88916844.2882014.0176020.02166312.847292
std2.3137781.2299005.52011164.046482248.0661924.5594775.80662516.3174691.7916530.29595963.655818
min1.0000002.00000018.7000001.1000007.9000000.0000002.20000015.0000000.4000000.0000000.000000
25%3.0000004.00000090.20000068.600000437.7000006.50000015.50000033.0000002.7000000.0000000.000000
50%4.0000004.00000091.600000108.300000664.2000008.40000019.30000042.0000004.0000000.0000000.520000
75%7.0000005.00000092.900000142.400000713.90000010.80000022.80000053.0000004.9000000.0000006.570000
max9.0000009.00000096.200000291.300000860.60000056.10000033.300000100.0000009.4000006.4000001090.840000
\n", "
" ], "text/plain": [ " X Y FFMC DMC DC ISI \\\n", "count 517.000000 517.000000 517.000000 517.000000 517.000000 517.000000 \n", "mean 4.669246 4.299807 90.644681 110.872340 547.940039 9.021663 \n", "std 2.313778 1.229900 5.520111 64.046482 248.066192 4.559477 \n", "min 1.000000 2.000000 18.700000 1.100000 7.900000 0.000000 \n", "25% 3.000000 4.000000 90.200000 68.600000 437.700000 6.500000 \n", "50% 4.000000 4.000000 91.600000 108.300000 664.200000 8.400000 \n", "75% 7.000000 5.000000 92.900000 142.400000 713.900000 10.800000 \n", "max 9.000000 9.000000 96.200000 291.300000 860.600000 56.100000 \n", "\n", " temp RH wind rain area \n", "count 517.000000 517.000000 517.000000 517.000000 517.000000 \n", "mean 18.889168 44.288201 4.017602 0.021663 12.847292 \n", "std 5.806625 16.317469 1.791653 0.295959 63.655818 \n", "min 2.200000 15.000000 0.400000 0.000000 0.000000 \n", "25% 15.500000 33.000000 2.700000 0.000000 0.000000 \n", "50% 19.300000 42.000000 4.000000 0.000000 0.520000 \n", "75% 22.800000 53.000000 4.900000 0.000000 6.570000 \n", "max 33.300000 100.000000 9.400000 6.400000 1090.840000 " ] }, "execution_count": 122, "metadata": {}, "output_type": "execute_result" } ], "source": [ "forest_fires.describe()" ] }, { "cell_type": "code", "execution_count": 123, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
XYFFMCDMCDCISItempRHwindrainarea
X1.0000000.539548-0.021039-0.048384-0.0859160.006210-0.0512580.0852230.0187980.0653870.063385
Y0.5395481.000000-0.0463080.007782-0.101178-0.024488-0.0241030.062221-0.0203410.0332340.044873
FFMC-0.021039-0.0463081.0000000.3826190.3305120.5318050.431532-0.300995-0.0284850.0567020.040122
DMC-0.0483840.0077820.3826191.0000000.6821920.3051280.4695940.073795-0.1053420.0747900.072994
DC-0.085916-0.1011780.3305120.6821921.0000000.2291540.496208-0.039192-0.2034660.0358610.049383
ISI0.006210-0.0244880.5318050.3051280.2291541.0000000.394287-0.1325170.1068260.0676680.008258
temp-0.051258-0.0241030.4315320.4695940.4962080.3942871.000000-0.527390-0.2271160.0694910.097844
RH0.0852230.062221-0.3009950.073795-0.039192-0.132517-0.5273901.0000000.0694100.099751-0.075519
wind0.018798-0.020341-0.028485-0.105342-0.2034660.106826-0.2271160.0694101.0000000.0611190.012317
rain0.0653870.0332340.0567020.0747900.0358610.0676680.0694910.0997510.0611191.000000-0.007366
area0.0633850.0448730.0401220.0729940.0493830.0082580.097844-0.0755190.012317-0.0073661.000000
\n", "
" ], "text/plain": [ " X Y FFMC DMC DC ISI temp \\\n", "X 1.000000 0.539548 -0.021039 -0.048384 -0.085916 0.006210 -0.051258 \n", "Y 0.539548 1.000000 -0.046308 0.007782 -0.101178 -0.024488 -0.024103 \n", "FFMC -0.021039 -0.046308 1.000000 0.382619 0.330512 0.531805 0.431532 \n", "DMC -0.048384 0.007782 0.382619 1.000000 0.682192 0.305128 0.469594 \n", "DC -0.085916 -0.101178 0.330512 0.682192 1.000000 0.229154 0.496208 \n", "ISI 0.006210 -0.024488 0.531805 0.305128 0.229154 1.000000 0.394287 \n", "temp -0.051258 -0.024103 0.431532 0.469594 0.496208 0.394287 1.000000 \n", "RH 0.085223 0.062221 -0.300995 0.073795 -0.039192 -0.132517 -0.527390 \n", "wind 0.018798 -0.020341 -0.028485 -0.105342 -0.203466 0.106826 -0.227116 \n", "rain 0.065387 0.033234 0.056702 0.074790 0.035861 0.067668 0.069491 \n", "area 0.063385 0.044873 0.040122 0.072994 0.049383 0.008258 0.097844 \n", "\n", " RH wind rain area \n", "X 0.085223 0.018798 0.065387 0.063385 \n", "Y 0.062221 -0.020341 0.033234 0.044873 \n", "FFMC -0.300995 -0.028485 0.056702 0.040122 \n", "DMC 0.073795 -0.105342 0.074790 0.072994 \n", "DC -0.039192 -0.203466 0.035861 0.049383 \n", "ISI -0.132517 0.106826 0.067668 0.008258 \n", "temp -0.527390 -0.227116 0.069491 0.097844 \n", "RH 1.000000 0.069410 0.099751 -0.075519 \n", "wind 0.069410 1.000000 0.061119 0.012317 \n", "rain 0.099751 0.061119 1.000000 -0.007366 \n", "area -0.075519 0.012317 -0.007366 1.000000 " ] }, "execution_count": 123, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Попарная корелляция между признаками\n", "forest_fires.corr()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Линейная модель" ] }, { "cell_type": "code", "execution_count": 124, "metadata": {}, "outputs": [], "source": [ "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.33, random_state=42)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Значение весов и параметр регуляризации" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Для начала расчитаем коэффициенты регрессии в нашей модели и поймем все ли признаки имеют некоторый вес в нашей модели" ] }, { "cell_type": "code", "execution_count": 125, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([0. , 0. , 0. , 0.01194366, 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. ])" ] }, "execution_count": 125, "metadata": {}, "output_type": "execute_result" } ], "source": [ "est = LassoLarsCV(cv=10, precompute=False).fit(X_train, y_train)\n", "est.coef_" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Как мы видим в данной модели только 1 признак имеет некоторый вес. Этот признак - DMC. Построим график зависимости веса от параметра регуляризации" ] }, { "cell_type": "code", "execution_count": 126, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl4VPXd/vH3JzuELUDYEiAZqrZq3ZrHarWLVq2VJAQREFBBUURQREVcUBYXFHFDFhWVRUUWAVmCa63V2ro0qK1rtUwCJAQSdsKa5fv7I/H5UR6RkJnkzEzu13XlYjI5mXPPV7w9npn5HHPOISIi4S/K6wAiIhIcKnQRkQihQhcRiRAqdBGRCKFCFxGJECp0EZEIoUIXEYkQKnQRkQihQhcRiRAxDbmztm3burS0tIbcpYhI2Fu9evVm51zykbZr0EJPS0sjLy+vIXcpIhL2zGxtbbbTKRcRkQihQhcRiRAqdBGRCKFCFxGJECp0EZEIoUIXEYkQKnQRkQgRFoW+++OP2Tp3Lq6y0usoIiIhKywKfdfrr7PpgQdZ238A+//zH6/jiIiEpLAo9PZ3302nyQ9xYO1a8nteTOmMGbgDB7yOJSISUsKi0M2MlllZ+Fbl0vz889n8xFTye/dh7+dfeB1NRCRkhEWhfy+mTRtSHn2E1BnTqdy2jYK+fdk0eTJVe/d6HU1ExHNhVejfa37uufhyV9KqVy+2PjcLf04Ouz/+2OtYIiKeCstCB4hu0YKO995Dlzmzocqx7oqBFI8fT2VZmdfRREQ8EbaF/r3EM87At3wZrQcNYvuil/FnZrHrL3/xOpaISIM7YqGb2SwzKzGz//MKpJmNMjNnZm3rJ17tRDVtSvvbbyNt/ktEN29G4dDrKLp1NBXbtnkZS0SkQdXmCH0OcOGhd5pZZ+B8YF2QM9VZk5NPJn3JEtoOH87O11/Hf1F3dqxahXPO62giIvXuiIXunHsP2PoDP3oMGA2EVFtaXBzJN1xP+uLFxKamsuGWURQOG075pk1eRxMRqVd1OoduZtlAkXPun0HOEzQJxx1L2oL5tBs9mt0ffIC/eybbFi3S0bqIRKyjLnQzawqMAcbWcvshZpZnZnmlpaVHu7uAWHQ0ba66Et/yZSQcfzwbx45j3aArObAuZM4SiYgETV2O0LsB6cA/zawASAU+MbMOP7Sxc26mcy7DOZeRnHzEi1bXi7iuXekyZzYd7pnAvi+/xJ/dgy2z52jYl4hElKMudOfc5865ds65NOdcGlAInOac2xj0dEFkUVEk9emDb1UuiWeeScmkSRT068++b7/1OpqISFDU5m2L84EPgOPMrNDMBtd/rPoT2749qTOm0+mRhylfv578XpdQOnWahn2JSNizhnyRMCMjw+Xl5TXY/o6kYts2Nt0/kZ25ucQf8xM63n8/TU46yetYIiL/xcxWO+cyjrRd2H9SNBAxSUmkPDyZ1CdnULlzFwWX9mPTg5M07EtEwlKjLvTvNT/nnOphX717s3XOHPzZPdj94UdexxIROSoq9BrRzZvTccJ4usydC1HGukGDKL57LJW7dnkdTUSkVlToh0j85en4li2j9eCr2L5kCf7umez685+9jiUickQq9B8Q1aQJ7W+9lbSFC4lu1YrCYcMpuvlmKrZs8TqaiMhhqdB/RJOfn0j64pdpO+IGdr71J/zdM9mxcqXGB4hISFKhH4HFxZE8bBi+pUuI7dqFDbeOpnDodZQXF3sdTUTkv6jQayn+mGNIe+kl2t9xO7s//hh/ZhbbFizAVVV5HU1EBFChHxWLjqb1wIH4Viwn4aSfs3H8BNYNHMSBggKvo4mIqNDrIq5zZ7rMmkXH++9j3zff4O+Rw5bnnsNVVHgdTUQaMRV6HZkZrXr1wpebS+LZZ1My+WEK+l7Kvm++8TqaiDRSKvQAxbZvR+q0qaQ8/hjlGzeSf0lvSqZMoUrDvkSkganQg8DMaHHhhfhyV9Ky+0VsefIp8ntezJ5PP/U6mog0Iir0IIpJSqLTpEl0nvk0VXv2sLb/ADZOnEjVnj1eRxORRkCFXg+a/eY3+FauJKlfP7Y9/wL+rGx2//3vXscSkQinQq8n0c0S6TD2brq++AIWE8O6qwazYcwYKnfs8DqaiEQoFXo9a5qRQfryZbS55hp2LFvOmsxMdr71ltexRCQCqdAbQFRCAu1uuZm0hQuJadOWohtGUDjyJio2b/Y6mohEkNpcU3SWmZWY2RcH3TfZzL4xs3+Z2Stm1qp+Y0aGJieeQPrLi0geOZKyt99mTfdMti9bpmFfIhIUtTlCnwNceMh9bwEnOudOAr4F7ghyrohlsbG0HXot6cuXEe/zUXz7Hawfci3lGzZ4HU1EwtwRC9059x6w9ZD73nTOff859w+B1HrIFtHifT66znuR9mPGsGf1avyZWWydN0/DvkSkzoJxDv0q4LUgPE6jY1FRtL78MnwrVtDklFPYdO99rL38Cvb7872OJiJhKKBCN7MxQAUw70e2GWJmeWaWV1paGsjuIlZcagqdn3uWjhMnsv+778jPyWHzzGdw5eVeRxORMFLnQjezgUAmMMD9yKt6zrmZzrkM51xGcnJyXXcX8cyMVhf3pNuqXJr99reUPvoo+X37su+rr7yOJiJhok6FbmYXArcB2c45fa49iGKSk0md+gQpU6ZQUVJKfu8+lDz2OFX793sdTURCXG3etjgf+AA4zswKzWwwMA1oDrxlZp+Z2VP1nLPRafGHC+iWu5KW2dlsefpp8nN6sueTT7yOJSIhzBryPdAZGRkuLy+vwfYXKcre/xsbx46lvLiYpAEDaHfTSKISE72OJSINxMxWO+cyjrSdPikaBpqdfRa+lStIGjCAbfPmsSYri7K/vu91LBEJMSr0MBGVmEiHu8bQdd6LRMUnsP6aa9hw+x1Ubt/udTQRCREq9DDT9LTTSF/2Cm2GXsuOlStZk5nFzjfe9DqWiIQAFXoYioqPp93IkaQvfpmYdskU3XgjhTeMoLykxOtoIuIhFXoYS/jZz0hftIjkW26m7N138WdmsX3pKxr2JdJIqdDDnMXE0Paaa0hftoz4Y46h+M47WT/4ag4UFnkdTUQamAo9QsT70un6wvO0H3s3ez/7DH92NluffwFXWel1NBFpICr0CGJRUbTu3x/fyhU0/cUv2DRxImsvu5z9a9Z4HU1EGoAKPQLFpqTQeebTdJr0IAf8fvJzerL5qac07EskwqnQI5SZ0bJHD3yrcml23u8pfXwK+b37sPfLL72OJiL1RIUe4WLatiX1scdInTaVii2bKejTl5JHHqFq3z6vo4lIkKnQG4nm551Ht9xcWvbMYcszz1YP+9JcHZGIokJvRKJbtqTTfffRZdZzuPJy1l52ORvvuYfKsjKvo4lIEKjQG6HEX/0K38oVtB54BdvmL8CflU3Ze+95HUtEAqRCb6Simjal/R130PWleUQlNmX9kGvZcNttVGzb5nU0EakjFXoj1/TUU0lfupS2w65jx6pX8XfPZOdrr2l8gEgYUqELUXFxJI8YQfqSxcR27EjRTTdTeP0NlG/SsC+RcKJCl/+VcNxxpC1cQLtbR7H7/ffxZ2ayffFiHa2LhInaXFN0lpmVmNkXB93X2szeMrPvav5Mqt+Y0lAsJoY2gwfjW76MhOOOo/iuu1l35VUcWL/e62gicgS1OUKfA1x4yH23A287544B3q75XiJIXFoaXZ6fS4fx49n3+ef4s3uwde5cDfsSCWFHLHTn3HvA1kPu7gHMrbk9F8gJci4JARYVRdKlffGtyiXx9NPZ9MCDFPTvz/7vvvM6moj8gLqeQ2/vnCsGqPmz3eE2NLMhZpZnZnmlpaV13J14KbZDB1KfepJOkydTvnYd/ot7UTpjBu7AAa+jichB6v1FUefcTOdchnMuIzk5ub53J/XEzGiZlYnv1VW0OP98Nj8xlfxLerP388+9jiYiNepa6JvMrCNAzZ96f1sjEdO6NSmPPkLqjOlUbt9OQd9L2fTQZKr27vU6mkijV9dCXwEMrLk9EFgenDgSLpqfey6+Vbm06tWLrbNm4c/JYfdHH3sdS6RRq83bFucDHwDHmVmhmQ0GHgTON7PvgPNrvpdGJrp5czreew9d5syGKse6gQMpHjeeyl27vI4m0ihZQ35oJCMjw+VpZGtEqtq7l9InprJ17lxikpPpMGE8zX/3O69jiUQEM1vtnMs40nb6pKgERVSTJrS/bTRpC+YT3aIFhUOvo2jUrVRsPfQdryJSX1ToElRNTjqJ9CWLaXv99ex84w383TPZkbtK4wNEGoAKXYLO4uJIvn549bCvzp3ZMGoUhdcNo3zjRq+jiUQ0FbrUm4RjjyVt/ku0u+02dn/4If7MLLYtXISrqvI6mkhEUqFLvbLoaNpcOQjfiuUknHACG8eNY92gKzmwdq3X0UQijgpdGkRcly50mTObDvfew76vvsLfI4cts2Zr2JdIEKnQpcGYGUm9e1cP+zrzTEoeeoiCS/ux79tvvY4mEhFU6NLgYtu3J3XGdFIefYTyoiLye11C6dRpVGnYl0hAVOjiCTOjxUUX4VuVS4sLL2Tz9OkU9OrF3n/+0+toImFLhS6eiklKImXyQ6Q+9SSVu8oouLQfmx54kKo9e7yOJhJ2VOgSEpr/7nf4clfS6tK+bJ07F3+PHHZ/+KHXsUTCigpdQkZ0s2Z0HDeOLs/PhShj3aArKb77bip37vQ6mkhYUKFLyEk8/XR8y5fT5urBbF+yFH9mFrv+/GevY4mEPBW6hKSohATajRpF2sKFRLdqReGw4RTdfDMVW7Z4HU0kZKnQJaQ1+fmJpC9+meQbR7DrrT/hv6g7O1as0LAvkR+gQpeQZ3FxtL3uOtJfWUpcWhobRt/G+qFDKS8u9jqaSEhRoUvYiP/JT+j60jza33kHez7+R/Wwr/nzNexLpEZAhW5mN5nZl2b2hZnNN7OEYAUT+SEWHU3rK67At3IFTU4+iY0T7mHdFQM5UFDgdTQRz9W50M0sBRgBZDjnTgSigUuDFUzkx8SlptL5uefoeP997Pv3v6uHfT37LK6iwutoIp4J9JRLDNDEzGKApsCGwCOJ1I6Z0apXL3y5uST++mxKHn6Egj592ffNN15HE/FEnQvdOVcEPAysA4qBHc65N4MVTKS2Ytu3I3XqVFIef5zyTZvIv6Q3JVOmaNiXNDqBnHJJAnoA6UAnINHMLvuB7YaYWZ6Z5ZWWltY9qciPMDNaXPgHfLkradm9O1uefIr8nhez59NPvY4m0mACOeVyHpDvnCt1zpUDS4FfHbqRc26mcy7DOZeRnJwcwO5EjiwmKYlOkx6k8zMzqdq7h7X9B7Bx4kSqdu/2OppIvQuk0NcBZ5hZUzMz4PfA18GJJRKYZr/+Nb4VK0nq149tz7+AP7sHZX/7m9exROpVIOfQPwIWA58An9c81swg5RIJWHSzRDqMvZuuL76AxcayfvDVbBgzhsodO7yOJlIvrCE/Qp2RkeHy8vIabH8i36vav5/N06azZdYsolsn0WHsWFqcf77XsURqxcxWO+cyjrSdPikqjUJUfDztbrmZtEULiWmbTNENIyi8cSQVeqFeIogKXRqVJiecQPqihSTfdBNl77zDmswsti9bpmFfEhFU6NLoWGwsba8dQvqyV4j3+Si+/Q7WXzOE8qIir6OJBESFLo1WvM9H13kv0v6uu9jzySf4s7LZOm+ehn1J2FKhS6NmUVG0vmwAvhUraHLqqWy69z7WXnY5+/35XkcTOWoqdBEgLjWFzs8+Q8cHHmD/mjXk5+Sw+emZuPJyr6OJ1JoKXaSGmdGqZw7dclfS7JxzKH3sMfL79mXfV195HU2kVlToIoeISU4mdcrjpDwxhYqSUvJ796Hk0ceo2r/f62giP0qFLnIYLS64gG6rcmnZowdbZs4kP6cnez75xOtYIoelQhf5EdEtW9Jp4v10fvZZ3P79rB1wGRvvvY/KMg37ktCjQhephWZnn4Vv5QqSLruMbS+9hD87i7K/vu91LJH/okIXqaWoxEQ6jLmTrvPmEZXQhPXXXMOG226ncvt2r6OJACp0kaPW9LRTSX9lKW2GXsuOVatY0z2Tna+/4XUsERW6SF1ExcfTbuRI0l9eRGz79hSNHEnhDSMoLynxOpo0Yip0kQAk/OxnpC1aSPItN1P27rv4M7PYvmSphn2JJ1ToIgGymBjaXnMN6cuXEX/sMRSPGcP6wVdzoLDQ62jSyKjQRYIkPj2drs8/T4dxY9n72WfVw76efwFXWel1NGkkVOgiQWRRUST164cvdyVN/yeDTRMnsnbAZexfs8braNIIBFToZtbKzBab2Tdm9rWZnRmsYCLhLLZTJzo//TSdHprEgfx88nN6svmppzTsS+pVoEfoU4DXnXM/BU4Gvg48kkhkMDNaZmfje3UVzc77PaWPTyH/kt7s/eJLr6NJhKpzoZtZC+A3wHMAzrkDzjl9wkLkEDFt2pD62GOkTptK5datFPTtS8kjj1C1b5/X0STCBHKE7gNKgdlm9qmZPWtmiUHKJRJxmp93Hr5VubTsmcOWZ54lv0cOe/7xD69jSQQJpNBjgNOAJ51zpwK7gdsP3cjMhphZnpnlleoK69LIRbdoQaf77qPL7Fm4ykrWXn4FxRMmUFlW5nU0iQCBFHohUOic+6jm+8VUF/x/cc7NdM5lOOcykpOTA9idSORIPPNMfCuW03rgQLYvWIg/K5uyd9/1OpaEuToXunNuI7DezI6ruev3gC7tIlJLUU2b0v6O20mb/xJRiU1Zf+1QikaPpmLbNq+jSZgK9F0uNwDzzOxfwCnAxMAjiTQuTU45hfSlS2k7bBg7X30Nf/dMdr72msYHyFELqNCdc5/VnE45yTmX45zToYVIHUTFxZE84gbSlywmtlMnim66mcLrb6B8k4Z9Se3pk6IiISThuONIWzCfdrfeyu7338efmcm2l1/W0brUigpdJMRYTAxtBl+Fb8VyEn76UzbePZZ1V17FgfXrvY4mIU6FLhKi4rp2pcvcOXSYMIF9n3+OPyubLXPmaNiXHJYKXSSEWVQUSX374FuVS+Ivf0nJg5Mo6N+f/d9953U0CUEqdJEwENuhA6lPPUmnhx+mfN16/Bf3onT6dNyBA15HkxCiQhcJE2ZGy8zu+Fbl0uKCC9g8dRr5vS5h7+efex1NQoQKXSTMxLRuTcojD5M6YwaVO3dS0PdSNk16iKq9e72OJh5ToYuEqebnnoMvdyWtevdm6+zZ+HvksPujj72OJR5SoYuEsejmzek4YTxd5swBYN3AgRSPHUflrl3eBhNPqNBFIkDiGb/Et3wZra+6iu2LF+PPzGLXO+94HUsamApdJEJENWlC+9G3krZgPtEtWlB43TCKbhlFxdatXkeTBqJCF4kwTU46ifQli2l7w/XsfPNN/N0z2ZG7SuMDGgEVukgEsrg4kocPx7d0CbFdOrNh1CgKrxtG+caNXkeTeqRCF4lg8cccQ9pLL9Hu9tvY/eGH+Ltnsm3BQlxVldfRpB6o0EUinEVH02bQIHwrV5Dw85+zcfx41g26kgNr13odTYJMhS7SSMR17kyX2bPocO897PvqK/zZPdjy3CxcRYXX0SRIVOgijYiZkdS7d/Wwr7POomTyZAr69Wffv7/1OpoEgQpdpBGKbd+e1OnTSHnsUcqLisjv1YvSJ6ZSpWFfYS3gQjezaDP71MxygxFIRBqGmdHij3+sHvZ10R/ZPGMG+RdfzN7PPvM6mtRRMI7QbwS+DsLjiIgHYpKSSHnoITo//RRVZbsp6NefTQ88SNWePV5Hk6MUUKGbWSrQHXg2OHFExCvNfvvb6mFfl/Zl69y5+LN7sPuDD7yOJUch0CP0x4HRwGHf1GpmQ8wsz8zySktLA9ydiNSn6GbN6DhuHF1feB6LjmbdlVdRfPfdVO7c6XU0qYU6F7qZZQIlzrnVP7adc26mcy7DOZeRnJxc192JSANq+j//Q/ryZbS5ejDblyzF3z2TXW+/7XUsOYJAjtDPArLNrABYAJxrZi8GJZWIeC4qIYF2o0aRtnAh0a1bUzj8egpvuomKzZu9jiaHUedCd87d4ZxLdc6lAZcCf3bOXRa0ZCISEpr8/ETSF79M8sgbKfvT29XDvlas0LCvEKT3oYvIEVlsLG2HDiV92SvEpaezYfRtrL/2Wso3bPA6mhwkKIXunPuLcy4zGI8lIqErvls3us57kfZ33smef+Thz8xi2/z5GvYVInSELiJHxaKjaX3F5fhWrqDJKSezccI9rL3iCvbn53sdrdFToYtIncSlptL5uefoeP/97P/2O/JzerLl2Wc17MtDKnQRqTMzo1Wvi/HlrqTZb35NycOPUNCnL/u++cbraI2SCl1EAhbbrh2pU6eSMmUK5SUl5F/Sm5LHH6dq/36vozUqKnQRCZoWf7iAbrkraZmZyZannia/58Xs+eRTr2M1Gip0EQmq6Fat6PTgA3R+5hmq9u1l7YABbLx/IlW7d3sdLeKp0EWkXjT79dn4VqwkqX9/tr3wAv7sHpT97W9ex4poKnQRqTfRzRLpcPdddJ33IhYXx/rBV7PhzjFU7tjhdbSIpEIXkXrX9Be/IH3ZK7QZMoQdy5ezJjOTnW++6XWsiKNCF5EGERUfT7ubbyL95UXEtE2maMSNFI64kQqN1Q4aFbqINKiE448nfdFCkm+6ibK//IU1mVlsf2WZhn0FgQpdRBqcxcbS9tohpC97hfhu3Si+4w7WXzOE8qIir6OFNRW6iHgm3uej64sv0P7uu9j7ySesycpm64vzNOyrjlToIuIpi4qi9YAB+FauoOlpp7HpvvtYe9nl7Pf7vY4WdlToIhISYlNS6PzMTDo++AD716whv0cOm596Glde7nW0sKFCF5GQYWa0ysmh26pcmp17LqWPP05+n77s++orr6OFBRW6iIScmLZtSZ3yOClPTKFicyn5vftQ8sijGvZ1BCp0EQlZLS64gG65ubTs0YMtzzxDfo8c9qxe7XWskFXnQjezzmb2jpl9bWZfmtmNwQwmIgIQ3bIlnSbeT+fnnsWVl7N2wGVsvOdeKss07OtQgRyhVwC3OOd+BpwBDDez44MTS0TkvzU76yx8K5aTdMXlbJs/H392FmV//avXsUJKnQvdOVfsnPuk5vYu4GsgJVjBREQOFZWYSIc776TrS/OIatKU9dcMYcNtt1OxbZvX0UJCUM6hm1kacCrw0Q/8bIiZ5ZlZXqlmNohIEDQ99VTSX1lKm+uGsmPVKvyZWex8/Y1GPz4g4EI3s2bAEmCkc27noT93zs10zmU45zKSk5MD3Z2ICABRcXG0u/FG0he/TGyHDhSNHEnRiBGUl5R4Hc0zARW6mcVSXebznHNLgxNJRKT2En76U9IWLqDdqFsoe++v+DOz2L5kaaM8Wg/kXS4GPAd87Zx7NHiRRESOjsXE0Obqq0lf9goJxx5L8ZgxrB88mAOFhV5Ha1CBHKGfBVwOnGtmn9V8XRSkXCIiRy0+PZ0uz8+lw/hx7P3nv/BnZbP1+edxlZVeR2sQ1pD/W5KRkeHy8vIabH8i0niVFxdTPG4cu9/7K01OOYWO991L/E9+4nWsOjGz1c65jCNtp0+KikhEiu3Ykc5PP02nyQ9xoKCA/J4Xs/nJJyN62JcKXUQilpnRMisL36pcmp9/HqVTniD/kt7s/eJLr6PVCxW6iES8mDZtSHn0UVKnT6Ny2zYK+vSh5OGHqdq3z+toQaVCF5FGo/nvf48vdyWtevViy7PPkd8jh90ff+x1rKBRoYtIoxLdogUd772HLnNm46qqWHfFQIrHj6eyrMzraAFToYtIo5R4xhn4li+j9aBBbF/0Mv7MLMrefdfrWAFRoYtIoxXVtCntb7+NtPkvEdUskfXXDqXo1tFhO+xLhS4ijV6Tk08mfelS2g4fzs7XXsPfPZOdr74aduMDVOgiIlQP+0q+4XrSlywhNiWFoptvoXD49ZRv2uR1tFpToYuIHCThuGNJWzCfdqNHs/vvf8ffPZNtixaFxdG6Cl1E5BAWHU2bq67Et3wZCccfz8ax41g36EoOrFvndbQfpUIXETmMuK5d6TJnNh0mTGDfl1/iz+7BltlzQnbYlwpdRORHWFQUSX374MtdSeIZZ1AyaRIF/fqz79tvvY72f6jQRURqIbZDB1KfnEGnRx6mfP168ntdQum06bgDB7yO9r9U6CIitWRmtOzeHd+rq2jxhz+wedo08ntdwt5//cvraIAKXUTkqMUkJZHy8GRSn5xB5c6dFFzaj02THqJq715Pc6nQRUTqqPk551QP++rdm62zZ+PP7sHuDz/yLE+gF4m+0Mz+bWb/MbPbgxVKRCRcRDdvTscJ4+kydy6YsW7QIIrHjqNy164GzxLIRaKjgenAH4HjgX5mdnywgomIhJPEX55ePezrqqvYvngx/u6Z7PrzOw2aIZAj9NOB/zjn/M65A8ACoEdwYomIhJ+oJk1oP/pW0hYuILpVKwqHDaPollFUbN3aIPuPCeB3U4D1B31fCPwysDg/7C/DnqJlXOv6eGgRkfrhGwA+2FcJJff8iR0HtvK7GUPrdZeBHKHbD9z3f4YdmNkQM8szs7zS0tIAdiciIj8mkCP0QqDzQd+nAhsO3cg5NxOYCZCRkVGn6Tb1/V81EZFIEMgR+j+AY8ws3czigEuBFcGJJSIiR6vOR+jOuQozux54A4gGZjnnvgxaMhEROSqBnHLBOfcq8GqQsoiISAD0SVERkQihQhcRiRAqdBGRCKFCFxGJECp0EZEIYQ15JWszKwXWBvAQbYHNQYoTabQ2h6e1OTytzeGF0tp0dc4lH2mjBi30QJlZnnMuw+scoUhrc3ham8PT2hxeOK6NTrmIiEQIFbqISIQIt0Kf6XWAEKa1OTytzeFpbQ4v7NYmrM6hi4jI4YXbEbqIiBxGSBT6kS42bWY3m9lXZvYvM3vbzLoe9LOHzOxLM/vazJ4wsx+68EbYCnBtJpnZFzVffRs2ecOoxfoMNbPPzewzM3v/4OvemtkdNb/3bzP7Q8Mmr391XRsza2Nm75hZmZlNa/jk9S+AtTnfzFbX/Gy1mZ3b8Ol/hHOqoykaAAADBklEQVTO0y+qR++uAXxAHPBP4PhDtjkHaFpz+zpgYc3tXwF/q3mMaOAD4HdeP6cQWZvuwFtUT9RMBPKAFl4/Jw/Wp8VBt7OB12tuH1+zfTyQXvM40V4/pxBZm0TgbGAoMM3r5xJia3Mq0Knm9olAkdfP5+CvUDhCP+LFpp1z7zjn9tR8+yHVV0eC6kveJVD9DyUeiAU2NUjqhhHI2hwPvOucq3DO7ab6L+2FDZS7odRmfXYe9G0i//8yiT2ABc65/c65fOA/NY8XKeq8Ns653c6594F9DRW2gQWyNp86576/MtuXQIKZxTdA5loJhUL/oYtNp/zI9oOB1wCccx8A7wDFNV9vOOe+rqecXqjz2lBd4H80s6Zm1pbqI/nOh/3N8FSr9TGz4Wa2BngIGHE0vxvGAlmbSBestekFfOqc218vKesgFAq9VhebBjCzy4AMYHLN9z8Bfkb1UWkKcK6Z/aaecnqhzmvjnHuT6ouP/B2YT/XpqIr6iemZWq2Pc266c64bcBtw19H8bhgLZG0iXcBrY2YnAJOAa+slYR2FQqHX6mLTZnYeMAbIPui/iD2BD51zZc65MqqPTs+o57wNKZC1wTl3v3PuFOfc+VT/Jf6unvM2tFqtz0EWADl1/N1wE8jaRLqA1sbMUoFXgCucc2vqJWFdeX0Sn+oX7fxUvzD1/QsUJxyyzalUv4hxzCH39wX+VPMYscDbQJbXzylE1iYaaFNz+yTgCyDG6+fkwfocc9DtLCCv5vYJ/PeLon4i60XROq/NQfcNIjJfFA3k702rmu17ef08fvC5eR2gZpEuAr6tKaYxNffdQ/URJzWlvQn4rOZrRc390cDTwNfAV8CjXj+XEFqbhJo1+YrqF0tP8fq5eLQ+U6h+8eozql9vOeGg3x1T83v/Bv7o9XMJsbUpALYCZVQf0R7f0PlDcW2oPvWy+6B/3z4D2nn9fL7/0idFRUQiRCicQxcRkSBQoYuIRAgVuohIhFChi4hECBW6iEiEUKGLiEQIFbqISIRQoYuIRIj/B+EBGGgFsSMvAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "est.alphas_\n", "alphas = est.alphas_\n", "plt.plot(alphas, est.coef_path_.T);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Зависимости общей ошибки от параметра регуляризации" ] }, { "cell_type": "code", "execution_count": 127, "metadata": {}, "outputs": [], "source": [ "X = pd.get_dummies(X, drop_first=True)" ] }, { "cell_type": "code", "execution_count": 128, "metadata": {}, "outputs": [], "source": [ "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.33, random_state=42)" ] }, { "cell_type": "code", "execution_count": 129, "metadata": {}, "outputs": [], "source": [ "errors = []\n", "for i in np.linspace(1e-1, 1e+2, 30):\n", " est = sklearn.linear_model.Lasso(alpha=i)\n", " errors.append(sklearn.model_selection.cross_val_score(est, X, y).mean()**2)" ] }, { "cell_type": "code", "execution_count": 130, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X2QXfV93/H35z7svZJ2hXhYu0QilRKrjTFtZFvFatzpuOABQTIR6cAUJg2qRx2lHjx1mrQx+B/iB6bxTBMSpjYdYhSEJzVmiFM0rhyqATxJZmxgCYQHyx424BgZGS2WkPWAdvfe++0f53d3z67uo7Srhb2f18yde873/M7D1dXoo3N+556fIgIzM7NeFJb6AMzM7J3DoWFmZj1zaJiZWc8cGmZm1jOHhpmZ9cyhYWZmPXNomJlZzxwaZmbWM4eGmZn1rLTUB7DQLrrooli/fv1SH4aZ2TvK008//UZEjHZrt+xCY/369YyNjS31YZiZvaNI+ode2vnylJmZ9cyhYWZmPXNomJlZzxwaZmbWM4eGmZn1zKFhZmY9c2iYmVnPHBrJo/tf50vfGl/qwzAze1tzaCR//dIb/K9v/f1SH4aZ2duaQyMZqZY4PlkjIpb6UMzM3rYcGslwpUQj4ORUfakPxczsbcuhkQxXs8dwHZ+sLfGRmJm9fTk0kpFqGYBjpxwaZmbtODSSkUp2pnHs1PQSH4mZ2duXQyMZ8eUpM7OuHBpJs0/Dl6fMzNrrOTQkFSU9I+kbaX6DpCckvSTpa5KGUr2S5sfT8vW5bdyW6t+XdHWuvjXVxiXdmqu33MdiGE6Xp447NMzM2urnTOOTwP7c/BeAOyNiI3AE2JHqO4AjEfEe4M7UDkmXAjcC7wO2Al9KQVQEvghcA1wK3JTadtrHgpvpCPflKTOztnoKDUnrgF8GvpzmBVwBPJSa7AauS9Pb0jxp+ZWp/TbggYiYjIhXgHHg8vQaj4iXI2IKeADY1mUfC27YHeFmZl31eqbxR8DvAo00fyHwZkQ0/1t+AFibptcCrwKk5UdT+5n6vHXa1TvtY8EVC2LVUNGXp8zMOugaGpJ+BTgUEU/nyy2aRpdlC1VvdYw7JY1JGpuYmGjVpCfD1ZI7ws3MOujlTOPDwK9K+gHZpaMryM481kgqpTbrgNfS9AHgEoC0/DzgcL4+b5129Tc67GOOiLgnIjZHxObR0dEePlJrw5WSb7k1M+uga2hExG0RsS4i1pN1ZD8WEb8OPA5cn5ptBx5O03vSPGn5Y5E9BXAPcGO6u2oDsBF4EngK2JjulBpK+9iT1mm3j0UxUi27I9zMrIOz+Z3Gp4DfljRO1v9wb6rfC1yY6r8N3AoQES8CDwLfBf4SuCUi6qnP4hPAI2R3Zz2Y2nbax6IYqZbcEW5m1kGpe5NZEfEt4Ftp+mWyO5/mtzkF3NBm/TuAO1rU9wJ7W9Rb7mOxDFdK/PjoqXO1OzOzdxz/IjxnxB3hZmYdOTRyhitld4SbmXXg0Mhpjt7XaHj0PjOzVhwaOTNPup3y2YaZWSsOjRw/tNDMrDOHRo5H7zMz68yhkTM7Trh/q2Fm1opDI2fEAzGZmXXk0MiZHSfcoWFm1opDI2fY44SbmXXk0MiZ7Qh3n4aZWSsOjZyV5SKSb7k1M2vHoZFTKIjhSsmPRzcza8OhMc9IxQ8tNDNrx6Exz3C15MtTZmZtODTmGan6SbdmZu04NOYZrnj0PjOzdhwa84xU3RFuZtZO19CQVJX0pKS/k/SipM+k+n2SXpH0bHptSnVJukvSuKTnJH0gt63tkl5Kr+25+gclPZ/WuUuSUv0CSftS+32Szl/4P4K5PHqfmVl7vZxpTAJXRMQvApuArZK2pGX/LSI2pdezqXYNsDG9dgJ3QxYAwO3Ah8jG/b49FwJ3p7bN9bam+q3AoxGxEXg0zS+q4Yo7ws3M2ukaGpE5nmbL6dVpaLttwP1pve8AayRdDFwN7IuIwxFxBNhHFkAXA6sj4tsREcD9wHW5be1O07tz9UUzUi3z1nSdWr2x2LsyM3vH6alPQ1JR0rPAIbJ/+J9Ii+5Il6DulFRJtbXAq7nVD6Rap/qBFnWAd0fEQYD0/q6eP9kZmhmIyf0aZman6Sk0IqIeEZuAdcDlki4DbgN+AfgXwAXAp1JztdrEGdR7JmmnpDFJYxMTE/2sepphPx7dzKytvu6eiog3gW8BWyPiYLoENQn8KVk/BWRnCpfkVlsHvNalvq5FHeD1dPmK9H6ozXHdExGbI2Lz6OhoPx/pNKsdGmZmbfVy99SopDVpegXwUeB7uX/MRdbX8EJaZQ9wc7qLagtwNF1aegS4StL5qQP8KuCRtOyYpC1pWzcDD+e21bzLanuuvmiGK9mTbn15yszsdKUe2lwM7JZUJAuZByPiG5IekzRKdnnpWeA/pfZ7gWuBceAk8DGAiDgs6XPAU6ndZyPicJr+OHAfsAL4ZnoB/D7woKQdwA+BG870g/ZqxEO+mpm11TU0IuI54P0t6le0aR/ALW2W7QJ2taiPAZe1qP8EuLLbMS4k92mYmbXnX4TP4yFfzczac2jMMzt6n0PDzGw+h8Y81XKBYkHu0zAza8GhMY8kRjymhplZSw6NFoY9ep+ZWUsOjRY8TriZWWsOjRZWV8seiMnMrAWHRgvD1ZJ/EW5m1oJDowV3hJuZtebQaMEd4WZmrTk0Whj2OOFmZi05NFpYXS0zVWswWasv9aGYmb2tODRamBm9z5eozMzmcGi0MPt4dIeGmVmeQ6OFYT/p1sysJYdGCx5Tw8ysNYdGC6tnHo/uX4WbmeU5NFqY6Qh3n4aZ2RwOjRbcEW5m1lrX0JBUlfSkpL+T9KKkz6T6BklPSHpJ0tckDaV6Jc2Pp+Xrc9u6LdW/L+nqXH1rqo1LujVXb7mPxeY+DTOz1no505gEroiIXwQ2AVslbQG+ANwZERuBI8CO1H4HcCQi3gPcmdoh6VLgRuB9wFbgS5KKkorAF4FrgEuBm1JbOuxjUVVKRYaKBYeGmdk8XUMjMsfTbDm9ArgCeCjVdwPXpeltaZ60/EpJSvUHImIyIl4BxoHL02s8Il6OiCngAWBbWqfdPhbdSLXkjnAzs3l66tNIZwTPAoeAfcDfA29GRPO/4geAtWl6LfAqQFp+FLgwX5+3Trv6hR32Mf/4dkoakzQ2MTHRy0fqyo9HNzM7XU+hERH1iNgErCM7M3hvq2bpXW2WLVS91fHdExGbI2Lz6OhoqyZ9G6748ehmZvP1dfdURLwJfAvYAqyRVEqL1gGvpekDwCUAafl5wOF8fd467epvdNjHossuTzk0zMzyerl7alTSmjS9AvgosB94HLg+NdsOPJym96R50vLHIiJS/cZ0d9UGYCPwJPAUsDHdKTVE1lm+J63Tbh+LbrhS9uPRzczmKXVvwsXA7nSXUwF4MCK+Iem7wAOSPg88A9yb2t8LfEXSONkZxo0AEfGipAeB7wI14JaIqANI+gTwCFAEdkXEi2lbn2qzj0W3ulrie+4INzObo2toRMRzwPtb1F8m69+YXz8F3NBmW3cAd7So7wX29rqPc8Ed4WZmp/MvwttodoRnV8nMzAwcGm2NVMvUGsGp6cZSH4qZ2duGQ6ONmUeJTLpfw8ysyaHRxmo/f8rM7DQOjTY8TriZ2ekcGm14TA0zs9M5NNoY8eh9ZmancWi0MeI+DTOz0zg02nBomJmdzqHRxir3aZiZncah0Ua5WKBaLjg0zMxyHBodjFTL7gg3M8txaHQwUvGYGmZmeQ6NDjwQk5nZXA6NDvx4dDOzuRwaHXiccDOzuRwaHbgj3MxsLodGB8OVkscJNzPL6Roaki6R9Lik/ZJelPTJVP89ST+S9Gx6XZtb5zZJ45K+L+nqXH1rqo1LujVX3yDpCUkvSfqapKFUr6T58bR8/UJ++G5GUp9Go+HR+8zMoLczjRrwOxHxXmALcIukS9OyOyNiU3rtBUjLbgTeB2wFviSpKKkIfBG4BrgUuCm3nS+kbW0EjgA7Un0HcCQi3gPcmdqdMyPVEhFwcrp+LndrZva21TU0IuJgRPxtmj4G7AfWdlhlG/BARExGxCvAOHB5eo1HxMsRMQU8AGyTJOAK4KG0/m7guty2dqfph4ArU/tzYriSPenWneFmZpm++jTS5aH3A0+k0ickPSdpl6TzU20t8GputQOp1q5+IfBmRNTm1edsKy0/mtqfE7MPLXRnuJkZ9BEakoaBPwd+KyJ+CtwN/DywCTgI/EGzaYvV4wzqnbY1/9h2ShqTNDYxMdHxc/Rjdpxwn2mYmUGPoSGpTBYYfxYRXweIiNcjoh4RDeBPyC4/QXamcElu9XXAax3qbwBrJJXm1edsKy0/Dzg8//gi4p6I2BwRm0dHR3v5SD0Zqfjx6GZmeb3cPSXgXmB/RPxhrn5xrtmvAS+k6T3AjenOpw3ARuBJ4ClgY7pTaoiss3xPRATwOHB9Wn878HBuW9vT9PXAY6n9OdEcvc99GmZmmVL3JnwY+A3geUnPptqnye5+2kR2uegHwG8CRMSLkh4Evkt259UtEVEHkPQJ4BGgCOyKiBfT9j4FPCDp88AzZCFFev+KpHGyM4wbz+Kz9q15eer4pPs0zMygh9CIiL+hdd/C3g7r3AHc0aK+t9V6EfEys5e38vVTwA3djnGxePQ+M7O5/IvwDlYNOTTMzPIcGh0UC2LVUNGhYWaWODS6GKmW3adhZpY4NLrwmBpmZrMcGl149D4zs1kOjS6GPU64mdkMh0YX2ZmG+zTMzMCh0dVIpew+DTOzxKHRxXDV44SbmTU5NLoYqZY4MVWn7tH7zMwcGt0MV5rPn/LZhpmZQ6MLD8RkZjbLodHFzOPRfaZhZubQ6Gbm8pQ7w83MHBrdDPvx6GZmMxwaXaz2OOFmZjMcGl0MV7I+DXeEm5k5NLpq3j3lPg0zM4dGVyuHiki+e8rMDHoIDUmXSHpc0n5JL0r6ZKpfIGmfpJfS+/mpLkl3SRqX9JykD+S2tT21f0nS9lz9g5KeT+vcJUmd9nEuSfKTbs3Mkl7ONGrA70TEe4EtwC2SLgVuBR6NiI3Ao2ke4BpgY3rtBO6GLACA24EPAZcDt+dC4O7Utrne1lRvt49zanW17NAwM6OH0IiIgxHxt2n6GLAfWAtsA3anZruB69L0NuD+yHwHWCPpYuBqYF9EHI6II8A+YGtatjoivh0RAdw/b1ut9nFODVdKHvLVzIw++zQkrQfeDzwBvDsiDkIWLMC7UrO1wKu51Q6kWqf6gRZ1Ouxj/nHtlDQmaWxiYqKfj9QTj95nZpbpOTQkDQN/DvxWRPy0U9MWtTiDes8i4p6I2BwRm0dHR/tZtSceJ9zMLNNTaEgqkwXGn0XE11P59XRpifR+KNUPAJfkVl8HvNalvq5FvdM+zil3hJuZZXq5e0rAvcD+iPjD3KI9QPMOqO3Aw7n6zekuqi3A0XRp6RHgKknnpw7wq4BH0rJjkrakfd08b1ut9nFOjbgj3MwMgFIPbT4M/AbwvKRnU+3TwO8DD0raAfwQuCEt2wtcC4wDJ4GPAUTEYUmfA55K7T4bEYfT9MeB+4AVwDfTiw77OKdGqu4INzODHkIjIv6G1v0OAFe2aB/ALW22tQvY1aI+BlzWov6TVvs410YqJU5NN5iuNygX/XtIMxtc/hewB8N+lIiZGeDQ6ElzTA33a5jZoHNo9KA5et8x92uY2YBzaPTAT7o1M8s4NHow4tH7zMwAh0ZPZsYJ96/CzWzAOTR6MDtOuPs0zGywOTR6sHqmI9xnGmY22BwaPaiUCpQKcke4mQ08h0YPJPnx6GZmODR65sejm5k5NHo2XCm7I9zMBp5Do0e+PGVm5tDo2UjFl6fMzBwaPRr2mYaZmUOjVyPuCDczc2j0qtkRno0xZWY2mBwaPRqplpiuB5O1xlIfipnZkukaGpJ2STok6YVc7fck/UjSs+l1bW7ZbZLGJX1f0tW5+tZUG5d0a66+QdITkl6S9DVJQ6leSfPjafn6hfrQZ2Lm8ei+RGVmA6yXM437gK0t6ndGxKb02gsg6VLgRuB9aZ0vSSpKKgJfBK4BLgVuSm0BvpC2tRE4AuxI9R3AkYh4D3BnardkPHqfmVkPoRERfwUc7nF724AHImIyIl4BxoHL02s8Il6OiCngAWCbJAFXAA+l9XcD1+W2tTtNPwRcmdoviebofX7+lJkNsrPp0/iEpOfS5avzU20t8GquzYFUa1e/EHgzImrz6nO2lZYfTe2XxOyZhn8VbmaD60xD427g54FNwEHgD1K91ZlAnEG907ZOI2mnpDFJYxMTE52O+4zNjN7nPg0zG2BnFBoR8XpE1COiAfwJ2eUnyM4ULsk1XQe81qH+BrBGUmlefc620vLzaHOZLCLuiYjNEbF5dHT0TD5SVx4n3MzsDEND0sW52V8DmndW7QFuTHc+bQA2Ak8CTwEb051SQ2Sd5Xsi+9HD48D1af3twMO5bW1P09cDj8US/kjCl6fMzKDUrYGkrwIfAS6SdAC4HfiIpE1kl4t+APwmQES8KOlB4LtADbglIuppO58AHgGKwK6IeDHt4lPAA5I+DzwD3Jvq9wJfkTROdoZx41l/2rMw7Ftuzcy6h0ZE3NSifG+LWrP9HcAdLep7gb0t6i8ze3krXz8F3NDt+M6VSqnIUKngW27NbKD5F+F9WF0tuSPczAaaQ6MPw5WSO8LNbKA5NPqQPR7dHeFmNrgcGn0YqZTdEW5mA82h0QcPxGRmg86h0QePE25mg86h0QePE25mg86h0YfhNOSrR+8zs0Hl0OjDSLVMvRG8NV1f6kMxM1sSDo0+eCAmMxt0Do0+zDwe3aFhZgPKodEHjxNuZoPOodGH4Uo25Kt/FW5mg8qh0QcPxGRmg86h0Qd3hJvZoHNo9MHjhJvZoHNo9KF5puHLU2Y2qBwafSgVC6woF90RbmYDq2toSNol6ZCkF3K1CyTtk/RSej8/1SXpLknjkp6T9IHcOttT+5ckbc/VPyjp+bTOXZLUaR9LbaTq50+Z2eDq5UzjPmDrvNqtwKMRsRF4NM0DXANsTK+dwN2QBQBwO/AhsvHAb8+FwN2pbXO9rV32saT8eHQzG2RdQyMi/go4PK+8DdidpncD1+Xq90fmO8AaSRcDVwP7IuJwRBwB9gFb07LVEfHtyJ4CeP+8bbXax5IaqXiccDMbXGfap/HuiDgIkN7fleprgVdz7Q6kWqf6gRb1TvtYUiPVMsfdp2FmA2qhO8LVohZnUO9vp9JOSWOSxiYmJvpdvS/DFV+eMrPBdaah8Xq6tER6P5TqB4BLcu3WAa91qa9rUe+0j9NExD0RsTkiNo+Ojp7hR+qNO8LNbJCdaWjsAZp3QG0HHs7Vb053UW0BjqZLS48AV0k6P3WAXwU8kpYdk7Ql3TV187xttdrHknJHuJkNslK3BpK+CnwEuEjSAbK7oH4feFDSDuCHwA2p+V7gWmAcOAl8DCAiDkv6HPBUavfZiGh2rn+c7A6tFcA304sO+1hSzSFfG42gUGh1dc3MbPnqGhoRcVObRVe2aBvALW22swvY1aI+BlzWov6TVvtYaiPV7Em3J6ZqM9NmZoPCvwjv07AHYjKzAebQ6JMHYjKzQebQ6NPs49H9Ww0zGzwOjT55nHAzG2QOjT41O799ecrMBlHXu6dsrublqV1/8wo/PnqKD224kEt/ZjVF335rZgPAodGnf7S6yn/4pfU89r1DfP7/7geyINm8/nwu33ABH9pwAf9s7RqGSj6JM7PlR9lPK5aPzZs3x9jY2DnZ18Gjb/HkK4d54pXDPPnKYcYPHQegWi7wgZ/NQuSfvHuECGhE0Iig3she2TTUI2ik2lvTdU5M1jg5Vef4ZI2TUzVOTGa1E1PNZTXqjWDlUIlVlSIrh0oMV0qsHCqyqpLVVg2VWJVqI9USw5Vy9l4tMVJJ79UyK8tF/0DRzACQ9HREbO7azqGxcN44PsnYD7IQeeLlw+z/8U/p94+3WBCrhopZEFSyf/xXDTXDocjKSomixImpGicn65yYqs0Lmux9qtboui8JhodKM4EyXCkxXC1nwVKZrY1U5843p1elwFpVKfnMyuwdrtfQ8OWpBXTRcIWtl13M1ssuBuDoW9McOHKSYkEUJQrpvViYnS4UmKmtGCoyVCyQBi88K9P1Bicn6xyfqnHs1DTHT9U4Nlnj2Kkax0/VOD45zbFTaX4ya3Niss7Rt6b50ZGTHJ/M2p2Yqve0v6FSYSZQVlWygFvVnB46vdY8IxqemZ89c1o1VKRUdAiZvR05NBbReSvKnLfivCXZd7lY4LyVBc5bWSZ7rNeZqTdi5mzm+KkaPz2VTZ+YzEKoWT8+lUJmssbxyTrHJ6f5yfEpfviTLICal9h6VSkVZi6x5S+/rRzKwmZlCp3mZboVzXpql59fmc7UquWFCWSzQebQsI6KBbG6WmZ1tQxnmX+NXL9NFiT1XKBkl9ZOpPrJqVoKq7n9PK//9BQnp+ozbSd7uAzXJMHKcpEVc8IkC5QVQ0VWlLP5FXOCZrbdinJxpt3892q5SKXkULLlz6Fh50yhoJlLUQs1DGOt3uDkdJ23pmbDpfmevWan32oGU266GWI/OTE1W5uqc3K6Tr3RX4dUQcwESbU8GyrVUpHqUJEV5cKc5c021VSvzKvlpyul5jpZvezLd7ZEHBr2jlYqFlhdLGRnQgsoIpiqN3grBUszTN6abv9+qtku1U6l2lvTWV/R60dn2zXfp+tndiNKsSCqpUIWNKUsSCrpbKcZLJVUr5aKVObVKrl1h0pZKFXKhaxemt1Oc7rZZqhU8G+SBpxDw6wFSekfzCJrFnE/0/VGCpfGTMicmm5wqpYLnlqDU1N1TtXqTDbb1WbXmazNbmOylq1z+MRU1ra5Tu79bG+YLBY0EyRDxQKVcvY+NCdgmrXZdqdN52rl4ux65WK2rDzTRgwVi5RLml1WLFAuanbdYsG3j58jDg2zJVRO/wCOVM/N/iKC6XrMCaCpeoPJFDiTtQaTtQZTtTQ/3ZgJpan6bH0q1645PZlb1rzte6rWmFlvpm2aX2jFgigXRbmQBU65KEqFZihl0+VSgXIhC59SUQyl91IxX8+3z7ZXKopSIWtXSu2a+yvNLM/ey0VRLGTtmm2a883t5JfPvhcoFjWn/nbsI3NomA0QSdn/3EsFOEdB1UozvKZTgEzXs9CZrmchM10Lpur1VAumag1q9dkAaq6bbz8znV61eqT5oDZTn13v5HRWr9WD6Ub2nu0jqKX5Zts+u7cWTEFQKhQoFNK7skuyBc2GS/NVEPz3f/vPuXzDBYt6TA4NMzvn8uG1qrLUR9NdoxHUGlmYNEOo3gimG81Ayp7qMJ3qtcZsm1qjxbJGI3siRKMxs7yWtpHNN3LrzX2aRLNN85gaMTu/qlJc9D8Lh4aZWReFghgqiCE/GPzs/gQk/UDS85KelTSWahdI2ifppfR+fqpL0l2SxiU9J+kDue1sT+1fkrQ9V/9g2v54Wvftd4HPzGyALERs/puI2JR7ZsmtwKMRsRF4NM0DXANsTK+dwN2QhQxwO/Ah4HLg9mbQpDY7c+ttXYDjNTOzM7QY51rbgN1pejdwXa5+f2S+A6yRdDFwNbAvIg5HxBFgH7A1LVsdEd+O7KmK9+e2ZWZmS+BsQyOA/yfpaUk7U+3dEXEQIL03f/y7Fng1t+6BVOtUP9CifhpJOyWNSRqbmJg4y49kZmbtnG1H+Icj4jVJ7wL2Sfpeh7at+iPiDOqnFyPuAe6B7NHonQ/ZzMzO1FmdaUTEa+n9EPAXZH0Sr6dLS6T3Q6n5AeCS3OrrgNe61Ne1qJuZ2RI549CQtErSSHMauAp4AdgDNO+A2g48nKb3ADenu6i2AEfT5atHgKsknZ86wK8CHknLjknaku6aujm3LTMzWwJnc3nq3cBfpLtgS8D/joi/lPQU8KCkHcAPgRtS+73AtcA4cBL4GEBEHJb0OeCp1O6zEXE4TX8cuI9sQIhvppeZmS2RZTfcq6QJ4B/OcPWLgDcW8HDeCfyZB4M/82A4m8/8jyNitFujZRcaZ0PSWC9j5C4n/syDwZ95MJyLz+zfxJuZWc8cGmZm1jOHxlz3LPUBLAF/5sHgzzwYFv0zu0/DzMx65jMNMzPrmUMDkLRV0vfTI9hv7b7GO4+kSyQ9Lmm/pBclfTLVWz7KfjmRVJT0jKRvpPkNkp5In/lrkoaW+hgXkqQ1kh6S9L30ff/L5f49S/ov6e/1C5K+Kqm63L5nSbskHZL0Qq7W91AUZ2vgQ0NSEfgi2aPbLwVuknTp0h7VoqgBvxMR7wW2ALekz9nuUfbLySeB/bn5LwB3ps98BNixJEe1eP4Y+MuI+AXgF8k++7L9niWtBf4zsDkiLgOKwI0sv+/5Pk4fHqKvoSgWwsCHBtnzssYj4uWImAIeIHuM+7ISEQcj4m/T9DGyf0jW0v5R9suCpHXALwNfTvMCrgAeSk2W1WeWtBr418C9ABExFRFvssy/Z7KnUqyQVAJWAgdZZt9zRPwVcHheud+hKM6aQ6P9o9mXLUnrgfcDT9D+UfbLxR8Bvws00vyFwJsRUUvzy+37/jlgAvjTdEnuy+nZcMv2e46IHwH/g+yxRQeBo8DTLO/vuanfoSjOmkOjj0ewLweShoE/B34rIn661MezmCT9CnAoIp7Ol1s0XU7fdwn4AHB3RLwfOMEyuhTVSrqOvw3YAPwMsIrs8sx8y+l77mbR/p47NNo/mn3ZkVQmC4w/i4ivp3K7R9kvBx8GflXSD8guO15BduaxJl3GgOX3fR8ADkTEE2n+IbIQWc7f80eBVyJiIiKmga8Dv8Ty/p6b+h2K4qw5NLKn625Md1oMkXWg7VniY1pw6Vr+vcD+iPjD3KJ2j7J/x4uI2yJiXUSsJ/teH4uIXwceB65PzZbbZ/4x8Kqkf5pKVwLfZRl/z2SXpbZIWpn+njc/87L9nnP6HYrirPnHfYCka8n+B1oEdkXEHUt8SAtO0r8C/hp4ntnr+58m69d4EPhZ0qPsc4+mXzYkfQT4rxHxK5J+juzM4wLgGeDfR8TkUh7fQpK0iazjfwh4mWwYggLL+HuW9Bng35HL4AZBAAAAaElEQVTdJfgM8B/JruEvm+9Z0leBj5A9yfZ14Hbg/9Die03h+T/J7rY6CXwsIsYW5DgcGmZm1itfnjIzs545NMzMrGcODTMz65lDw8zMeubQMDOznjk0zMysZw4NMzPrmUPDzMx69v8BY50P7zcnoB4AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(np.linspace(1e-1, 1e+2, 30), errors);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Зависимость ошибки от числа признаков" ] }, { "cell_type": "code", "execution_count": 131, "metadata": {}, "outputs": [], "source": [ "mses1 = []\n", "mses2 = []\n", "for i in range(1,13):\n", " est=Lasso(alpha=10)\n", " est.fit(X_train.iloc[:,0:i],y_train)\n", " predictions = est.predict(X_test.iloc[:,0:i])\n", " mses1 = np.append(mses1, [mean_squared_error(y_test, predictions)])" ] }, { "cell_type": "code", "execution_count": 132, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xt4VPd95/H3V3cJ0GWQAAECDTbBxgiEgQGcxMZ2zDX1LdnGbmLTW9x6k7awbTZxdzd5kvTJNrHTJtm0ybqJGzt17BYXxc4ax3YujuNwB4HBxhgwAiSQEAiQQAh0+e0fc2TEVRoxM2cun9fzzOOZn8+c8z02zGfme875HXPOISIi6SfD7wJERMQfCgARkTSlABARSVMKABGRNKUAEBFJUwoAEZE0pQAQEUlTCgARkTSlABARSVNZfhdwJaWlpa6ystLvMkREksqmTZuOOOfK+lsuoQOgsrKSjRs3+l2GiEhSMbN9A1lOLSARkTSlABARSVMKABGRNKUAEBFJUwoAEZE0pQAQEUlTCgARkTSV0NcBSPo5cvIMT6/dT3dPT9y2OfeaUuZeMzxu2xNJFAoASSjff20PP3hjL2bx2Z5zsGJTPb/7/G1kZMRpoyIJQgEgCaOru4fntx7kjskj+ZcHZ8Zlmys31/Pf/mMrtQeOM2N8SVy2KZIodAxAEsbv9hylue0M904fE7dtfmTySHIyM1i17VDctimSKBQAkjB+WttAYV4Wt10/Im7bLMzL5sMTS3lp2yF6elzctiuSCBQAkhBOneni59sbWTJ1NLlZmXHd9uKqcg6e6GBL/fG4blfEbwoASQgvv9XI6c5u7r0xfu2fXh+ZPJLsTGPVm2oDSXpRAEhCqKltoCKQz0wfDsQW5Wfz4YllvLS9EefUBpL0oQAQ3zW1dvC73Ue4p3oMFq/zPy+wuKqchuOn2XJAbSBJHwoA8d3zWxrocXDPjWN9q+GO3jaQzgaSNKIAEN+t3NxAdUUxwdIhvtXQ2wZatU1tIEkfCgDx1Y5DrbzT2ObLwd8L9baBttaf8LsUkbhQAIivamobyMowPjp1tN+lqA0kaUcBIL7p7nE8v6WBeZNGEBiS43c5FOVn86FrS3nxzUNqA0laUACIb1bvOUJT65mEaP/06m0Dvak2kKQBBYD4pmZzA8PysrjtuvhN/dCf+ZNHqQ0kaUMBIL5oP9vFz99q5KNTy8nLju/UD1dSVJDNB68t5cVtagNJ6lMAiC9efquR9rPd3DPdv3P/L2dxVTn1x06zrUFtIEltCgDxxcrNDYwt8Wfqh/7MnzySrAzjRc0NJClOASBxd7h36ofpYxLyLlzFBTlqA0laUABI3D2/5WB46oc43vglUkvUBpI0oACQuFtZ28C0imImlA31u5TLmn+D1wbS2UCSwhQAElfvNLay41BrXG/7OBjFBTncdG0pq9QGkhTWbwCY2SQz29Ln0Wpmy8ys2szWemMbzSx0wftmmVm3mX28z9hSM9vlPZbGYocksdVsDk/98HvT/J/6oT9LqkZxoOU02xta/S5FJCb6DQDn3E7nXLVzrhqYAbQDNcA3gC9741/0XgNgZpnA14GX+4wFgC8Bs4EQ8CUzS7xTQCRmunscP93SwLxJZQkx9UN/5k8eRabaQJLCIm0B3Q7scc7tAxxQ6I0XAQf7LPcXwH8Ch/uMLQBedc61OOeOAa8CCwdVtSSlNXuO0tR6hrsTvP3Tq2RIDjddM1xtIElZkQbAfcAz3vNlwKNmdgB4DHgEwMzGAPcA37/gvWOAA31e13tj5zGzh7yW0sbm5uYIy5NEtrK2nmG5WXzk+pF+lzJgS6rK2d/SzlsH1QaS1DPgADCzHOBOYIU39DCw3DlXASwHfuiNfwv4vHOu+8JVXGK1F32tcs497pyb6ZybWVZWNtDyJMG1n+3i59sbWVyVWFM/9Gf+DWoDSeqK5BfAImCzc67Je70UWOk9X0G4rw8wE3jWzOqAjwP/bGZ3E/7GX9FnfWM5v20kKeyVt5rCUz8k0MyfAxFQG0hSWCQBcD/n2j8Q/vC+xXt+G7ALwDkXdM5VOucqgeeA/+qc+ynhA8LzzazEO/g7nz4HiSW1raxtYExxPqHKgN+lRGxxVTn7jqoNJKlnQAFgZgXAHZz7xg/waeCbZrYV+Brw0JXW4ZxrAb4KbPAeX/HGJMUdbu3gjV3N3D19dEJO/dCfBV4bSFNES6rJGshCzrl2YPgFY28QPi30Su/7wwtePwE8EVmJkuxe2No79UPizfw5EH3bQJ9bMAmz5AsxkUvRlcAScys3NzB1bBHXjkjcqR/6s7iqnLqj7bx9SG0gSR0KAImpnY1tvH2oNaEnfhsItYEkFSkAJKZW1taTmSRTP1xJYEgOcycMZ9W2Rp0NJClDASAx093jeL72ILd8oIzSobl+l3PVFleVs/fIKXYcavO7FJGoUABIzKx97yiNrR1J3/7pteCGkWoDSUpRAEjMrNzcwLDcLO6YnDxTP1zJ8KG5zJkQ0EVhkjIUABITp8928/Pth1hUNSqppn7oz+Kqct47cop3GtUGkuSnAJCYeOXtRk6d7U7ac/8vZ8ENo8gw1AaSlKAAkJhYuTk89cPsYPJN/XAlpUNzmTNhuG4YLylBASBRd7itg9/uauau6uSc+qE/i6vKea/5FDub1AaS5KYAkKj72dZD9Di4N8lm/hyohVPCbaAX31QbSJKbAkCirqa2nqoxRVw7YpjfpcRE6dBcZgfVBpLkpwCQqNrV1Mb2huSf+qE/i6eqDSTJTwEgUbWytoHMDOPO6uSe+qE/C3vPBlIbSJKYAkCipqfH8XxtAzdPLE2JqR+upGxYLqFgQG0gSWoKAImatXuPcvBEB/fcmFrn/l/Okqpy9jSf4t2mk36XIjIoCgCJmprNDQzNzWJ+ikz90J8FU0Zhhm4YL0lLASBRcfpsNy9tb2TRlNSa+uFKRgzLI1QZ0FXBkrQUABIVr+5o4uSZLu5J0XP/L2fJ1HJ2Hz7JuzobSJKQAkCiomZzPaOL8pgTHN7/wilkYW8bSGcDSRJSAMhVa247w+u7jnDX9DEpOfXDlagNJMlMASBX7WdbD9Ld47g3xS/+upwlU8vZdfgku9QGkiSjAJCrVlPbwJQxhUwcmZpTP/Rnoc4GkiSlAJCrsvtwG9saTnB3dXp++4dwG2iW2kCShBQAclVWbm4gw0j5qR/6s6SqnHebTrL7sNpAkjwUADJoPT2O57cc5MMTyxgxLM/vcny16P2zgRr9LkVkwBQAMmjr9rbQcPx0ys77H4kRhXnMGq82kCQXBYAMWk1tPUNyMpk/eZTfpSSExVWj2NnUxu7DmhtIkoMCQAalo7Obl7Y1snBKOfk56TH1Q38WVZVjumG8JBEFgAzKq2830XamS+2fPkYW5jFzfIkCQJKGAkAGpaa2gVGFecyZkF5TP/RncVU57zS2sadZbSBJfAoAidiRk2f4zbvN3DV9NJlpNvVDfxZNKQd0pzBJDgoAidi5qR/S48YvkRhVFG4D6apgSQYKAIlYTW0Dk8sLmTQqPad+6I/aQJIsFAASkd2HT/Jm/Qkd/L2CRVXh02LVBpJEpwCQiNTU1oenfpiW3lM/XEl5UT4z1AaSJNBvAJjZJDPb0ufRambLzKzazNZ6YxvNLOQtf5eZvdln/EN91rXUzHZ5j6Wx3DGJvp4ex09rD/KhiWWMKEzvqR/609sGek9tIElg/QaAc26nc67aOVcNzADagRrgG8CXvfEveq8BfglM88b/GPgBgJkFgC8Bs4EQ8CUzK4ny/kgMra/zpn5I03n/I7G4tw2kXwGSwCJtAd0O7HHO7QMcUOiNFwEHAZxzJ51zzhsf4i0HsAB41TnX4pw7BrwKLLya4iW+ajY3UJCTyfwbRvpdSsIrL8rnxnHFvLhNk8NJ4oo0AO4DnvGeLwMeNbMDwGPAI70Lmdk9ZvYO8CLhXwEAY4ADfdZV742dx8we8lpHG5ubmyMsT2Klo7ObVdsOsXDKKApysvwuJyksripnx6FW9h455XcpIpc04AAwsxzgTmCFN/QwsNw5VwEsB37Yu6xzrsY5dx1wN/DV3lVcYrXuogHnHnfOzXTOzSwrKxtoeRJjv9jhTf2gc/8HbHGVd1GY2kCSoCL5BbAI2Oyca/JeLwVWes9XEO7rn8c59zpwjZmVEv7GX9HnX4/FaxtJ4qvZ3MDIwlzmXqOpHwZqdLHXBtLpoJKgIgmA+znX/oHwh/ct3vPbgF0AZnatmZn3/EYgBzgKvAzMN7MS7+DvfG9MEtxRb+qHu6vHaOqHCC2uKuftQ63UqQ0kCWhAAWBmBcAdnPvGD/Bp4JtmthX4GvCQN/4xYLuZbQH+CfiEC2sh3A7a4D2+4o1JgvvZ1oN09Tju0cVfEettA+maAElEAzqa55xrB4ZfMPYG4dNCL1z268DXL7OeJ4AnIi9T/FRT28D15YVcN6qw/4XlPKOL85k+rphV2w7xmVuv9bsckfPoSmC5op2NbWytP6Fz/6/Ckqpy3jrYyr6jagNJYlEAyBU9taaO3KwMPj5DZ/8M1iK1gSRBKQDksk6c7mTl5gbunDaakiE5fpeTtMYU51NdUazTQSXhKADksp7bVM/pzm6W3lTpdylJb0lVOdsbWtl/tN3vUkTepwCQS+rpcfx4TR0zxpcwZUyR3+Ukvd4potUGkkSiAJBLen1XM3VH23lw7ni/S0kJY0sKmKY2kCQYBYBc0lNr9lE6NPf9e9zK1VtSNYptDSfUBpKEoQCQi+w7eopf7zzMH8weR06W/ohES2+YvvK2ZgiVxKC/3XKRf1u7j0wzPjl7nN+lpJSKQAGVwwtYt1cXwEtiUADIedrPdvHvGw6wYMooRuquX1EXCgbYUNdCT89FE+GKxJ0CQM7z/JaDtHZ0sXRupd+lpKRQcDjH2zvZdVi3ihT/KQDkfc45nlxdx/Xlhcyq1N06Y2F2MADA+r1Hfa5ERAEgfWyoO8Y7jW0snTseb0ZvibKxJfmUF+XpOIAkBAWAvO/JNXUU5mVxV7UmfosVMyMUDLB+bwvnbp0t4g8FgADQeKKDl7c38olZFeTnZPpdTkoLBQMcbjvD/hZdDyD+UgAIAD9Zv59u5/jUHF35G2u9xwHUBhK/KQCEs109/GTdfm6dNILxw4f4XU7Ku6ZsKIEhOaxXAIjPFADCS9sPceTkGc37EydmxqzKEgWA+E4BIDy5uo5g6RBunljmdylpIxQczv6Wdg6dOO13KZLGFABpblv9CTbvP84Dc8aTkaFTP+Pl3PUA+hUg/lEApLmn1tRRkJPJx3TLx7i6vryQoblZCgDxlQIgjR07dZbntx7knuljKMrP9ructJKZYczUcQDxmQIgjf37xgOc7erhQc3744tQMMCuwyc5evKM36VImlIApKnuHseP1+xjzoQAk0YN87uctNR7HGBD3TGfK5F0pQBIU7/c0UTD8dOa9dNHVWOKyc3KUBtIfKMASFNPrdlHeVEed0we6XcpaSsnK4Mbx5Wwvk4zg4o/FABpaPfhNt7YfYRPzRlPVqb+CPgpFAzw9sFWWjs6/S5F0pD+9qehH6/ZR05mBp+YVeF3KWlvdjBAj4NN+3QcQOJPAZBm2jo6eW5TPR+dWk7p0Fy/y0l708eVkJVhbNBxAPGBAiDN1NQ2cOpsNw/eVOl3KQLk52QydWyRDgSLLxQAaaT3lo/TxhZRXVHsdzniCQWHs7X+OB2d3X6XImlGAZBGfrf7KHuaT+nCrwQTCpbQ2e2o3X/c71IkzSgA0siTa+oYPiSHJVPL/S5F+pgxPoCZJoaT+FMApIkDLe38ckcT94UqyMvWLR8TSVF+NtePKtT1ABJ3CoA08fS6/QB8crZu+pKIQsEAm/Yd42xXj9+lSBpRAKSBjs5u/n3DfuZPHsXo4ny/y5FLmB0M0NHZw/aDJ/wuRdJIvwFgZpPMbEufR6uZLTOzajNb641tNLOQt/wnzexN77HazKb1WddCM9tpZrvN7Aux3DE552dbD3KsvZMHb9K3/0Q1SzeIER/0GwDOuZ3OuWrnXDUwA2gHaoBvAF/2xr/ovQbYC9zinJsKfBV4HMDMMoF/AhYBk4H7zWxylPdHLuCc48k1dUwcMZS5E4b7XY5cRunQXK4pG6IAkLiKtAV0O7DHObcPcEChN14EHARwzq12zvVe174W6L3VVAjY7Zx7zzl3FngWuOtqipf+bd5/nO0NrTx4UyVmuuVjIgsFh7OhroXuHud3KZImIg2A+4BnvOfLgEfN7ADwGPDIJZb/E+Al7/kY4ECff1fvjZ3HzB7yWkobm5ubIyxPLvTUmjqG5WZx7/SL/lNLgpkdDNDW0cU7ja1+lyJpYsABYGY5wJ3ACm/oYWC5c64CWA788ILlbyUcAJ/vHbrEai/6quOce9w5N9M5N7OsrGyg5cklHG7rYNW2Q3x85liG5Gb5XY70I6TjABJnkfwCWARsds41ea+XAiu95ysIt3gAMLOpwA+Au5xzvSc31wN9p58ci9c2kth4dv0BOrsdD8zRwd9kMLo4n7El+QoAiZtIAuB+zrV/IPzhfYv3/DZgF4CZjSMcDA84597ts/wGYKKZBb1fE/cBLwy2cLmyzu4enl63j5s/UMaEsqF+lyMDFAoGWL+3Bed0HEBib0ABYGYFwB2c+8YP8Gngm2a2Ffga8JA3/kVgOPDPvaeIAjjnuoDPAi8DO4D/cM69FZW9kIu88lYTTa1nWDpX3/6TyexggKOnzvLekVN+lyJpYECNYedcO+EP9b5jbxA+LfTCZf8U+NPLrGcVsCryMiVST66uoyKQz7xJI/wuRSIwq/LccYBr9MtNYkxXAqegtw+2sr6uhQfmjCczQ6d+JpNg6RBKh+bqOIDEhQIgBf14bR152Rn8/kzd8jHZmBmzveMAIrGmAEgxJ9o7qalt4O7qMRQX5PhdjgxCKBig4fhp6o+1+12KpDgFQIpZsekAHZ09PKCDv0lL1wNIvCgAUkhPj+OpNfuYVVnCDaOL/C5HBmnSyGEU5mUpACTmUvLy0KMnzzD3f/8qrtsszM/mz2+ZwANzx5Ob5c8NV37zbjP7W9r53IJJvmxfoiMjw96/HkAkllIyAPKyM/mTDwfjus3tDSf4uxd38KPVdXxuwSR+b+poMuJ8Bs6PVtcxYlguC24YFdftSvSFggF+seMwh9s6GDEsz+9yJEWlZAAMyc3i8wuvi/t2f7urma+teoe/enYLP/jtXh5ZfB03XVMal23vPXKK37zbzLKPTCQnS529ZBcKhi+72bD3mO7hLDGjT4oo+vDEMl78iw/xD78/jaMnz/AH/7KOP/7RBt5taov5tn+8Zh/ZmcYfzB4X821J7N0wupCCnEzW79V9giV2FABRlpFh3HvjWH71N/P4wqLr2FDXwsJvvc7nn3uTxhMdMdnmqTNdrNh0gEVTytUuSBHZmRnMGF/COh0HkBhSAMRIXnYmf37LNbz+uVv5ow8GWVlbz7zHfs1jL++kraMzqtv66ZYG2jq6WKpbPqaUUGWAnU1tHG8/63cpkqIUADFWMiSH//XRyfzqr+cxf/Iovvvr3cx79DWeWlNHZ3fPVa/fOcdTq/dxw+hCbhxXcvUFS8IIBQM4BxvrjvW/sMggKADipCJQwHfun84Ln/0gE0cO5YvPv8X8f3ydl7Yduqqpf9ftbWFnUxtL5+qWj6lmWkUxOZkZbKhTG0hiQwEQZ1PHFvPMp+fwr384i+xM4+GnN3Pv91azcZB/yZ9cXUdxQTZ3Vo+OcqXit7zsTKZVFOk4gMSMAsAHZsat141g1V9+mK9/rIqGY6f5+PfX8Gc/3sie5pMDXs/B46d55e0mPjGzgrxsfy4+k9gKBQNsbzjBqTNdfpciKUgB4KOszAw+MWscr31uHn8z/wP8bvdR5v/j6/zPn26jue1Mv+//ybr99DjHp3TLx5QVCg6nq8dRu/+436VIClIAJICCnCw+e9tEXvvcPD45exzPrj/AvEd/zbd/seuy3/zOdHXzzPr93H7dSCoCBXGuWOJlxvgSMgxdDyAxoQBIIKVDc/nKXVN4ZfnN3PyBMv7xF+8y77HX+Mm6/XRdcMbQqm2HOHrqrE79THFDc7OYMkbHASQ2FAAJaELZUL73qRn858M3MS5QwN/WbGPht3/LL95uev+MoSdX72NC2RA+GKepJsQ/ocoAtQeOc6ar2+9SJMUoABLYjPElPPfnc/m/D8ygp8fxp09t5BOPr+Un6/az5cBxHpwzPu4Tzkn8hYIBznb18Gb9Cb9LkRSjAEhwZsaCG0bx8vKb+erdU3iv+SR/W7ONITmZfGzGWL/Lkzjoe6N4kWhKydlAU1F2ZgYPzBnPPdPH8OTqOsqL8hiWl+13WRIHJUNymDRyGOv2tvCZW/2uRlKJAiDJDM3N4jO3Xut3GRJnoWCAlZvr6eruIStTP9wlOvQnSSQJhIIBTp3t5u1DrX6XIilEASCSBHSjeIkFBYBIEhhZmEfl8AJdDyBRpQAQSRKzKgNsqGuhp2fws8eK9KUAEEkSoWCA4+2d7I5gwkCRK1EAiCSJ2d6N4tUGkmhRAIgkiYpAPqMK83QgWKJGASCSJMyMUDDA+r1Hr+ouciK9FAAiSSQUDNDUeob9Le1+lyIpQAEgkkRme9cD6DiARIMCQCSJXDtiKIEhOToOIFGhABBJImbGrMoSBYBEhQJAJMmEgsPZ39LOoROn/S5Fkly/AWBmk8xsS59Hq5ktM7NqM1vrjW00s5C3/HVmtsbMzpjZ31ywroVmttPMdpvZF2K1UyKpbLbmBZIo6TcAnHM7nXPVzrlqYAbQDtQA3wC+7I1/0XsN0AL8JfBY3/WYWSbwT8AiYDJwv5lNjtaOiKSL68sLGZqbpQCQqxZpC+h2YI9zbh/ggEJvvAg4COCcO+yc2wB0XvDeELDbOfeec+4s8Cxw16ArF0lTmRnGTB0HkCiINADuA57xni8DHjWzA4S/7T/Sz3vHAAf6vK73xs5jZg95LaWNzc3NEZYnkh5mVQbYdfgkR0+e8bsUSWIDDgAzywHuBFZ4Qw8Dy51zFcBy4If9reISYxddzuice9w5N9M5N7OsrGyg5Ymkld7jABvqjvlciSSzSH4BLAI2O+eavNdLgZXe8xWEWzxXUg9U9Hk9Fq9tJCKRqRpbRG5WhtpAclUiCYD7Odf+gfCH9y3e89uAXf28fwMw0cyC3q+J+4AXIti+iHhyszKZPq6YDXUKABm8Ad0U3swKgDuAP+sz/Gng22aWBXQAD3nLjgI2Ej5A3GNmy4DJzrlWM/ss8DKQCTzhnHsransikmZCweF891e7aOvoZFhett/lSBIaUAA459qB4ReMvUH4tNALl20k3N651HpWAasiL1NELjQ7GOA7DjbtO8a8SSP8LkeSkK4EFklS08cVk5VhOg4gg6YAEElSBTlZVI0tUgDIoCkARJJYKBhga/1xOjq7/S5FkpACQCSJzQ4G6Ox21O4/7ncpkoQUACJJbMb4AGaaGE4GRwEgksSK8rO5flQh6+uO+l2KJCEFgEiSCwUDbNp3jLNdPX6XIklGASCS5ELBAB2dPWw/eMLvUiTJKABEktysSt0gRgZHASCS5MqG5TKhbIgCQCKmABBJAbODATbUtdDdc9EM6yKXpQAQSQGhYIC2ji52Nrb5XYokEQWASAoIBcNzNa7fq9NBZeAUACIpYExxPmOK81mv+wNIBBQAIilidjDA+r0tOKfjADIwCgCRFBEKBjhy8izvHTnldymSJBQAIikiFNT1ABIZBYBIigiWDqF0aK4CQAZMASCSIszs/eMAIgOhABBJIaFggIbjp6k/1u53KZIEFAAiKUTHASQSCgCRFDJp5DAK87IUADIgCgCRFJKRYcyq1HEAGRgFgEiKCQUDvHfkFIfbOvwuRRKcAkAkxfQeB9iw95jPlUiiUwCIpJgpY4rIz87UxHDSLwWASIrJzsxgxvgS1tfpF4BcmQJAJAWFggHeaWzlRHun36VIAlMAiKSgUDCAc7Bxn84GkstTAIikoOqKYnIyM3Q6qFyRAkAkBeVlZzKtooh1CgC5giy/CxCR2AgFA3zvtT3c8Q+/8bsUGYTrygv5P/dPj+k2FAAiKeq/zKjgQMtpunp6/C5FBqGiJD/m21AAiKSoytIhfCfG3yAluekYgIhImlIAiIikqX4DwMwmmdmWPo9WM1tmZtVmttYb22hmIW95M7PvmNluM3vTzG7ss66lZrbLeyyN5Y6JiMiV9XsMwDm3E6gGMLNMoAGoAf4F+LJz7iUzWwx8A5gHLAImeo/ZwPeA2WYWAL4EzAQcsMnMXnDO6Xp1EREfRNoCuh3Y45zbR/hDvNAbLwIOes/vAp5yYWuBYjMrBxYArzrnWrwP/VeBhVe9ByIiMiiRngV0H/CM93wZ8LKZPUY4SG7yxscAB/q8p94bu9z4eczsIeAhgHHjxkVYnoiIDNSAfwGYWQ5wJ7DCG3oYWO6cqwCWAz/sXfQSb3dXGD9/wLnHnXMznXMzy8rKBlqeiIhEKJIW0CJgs3OuyXu9FFjpPV8BhLzn9UBFn/eNJdweuty4iIj4IJIW0P2ca/9A+MP7FuA14DZglzf+AvBZM3uW8EHgE865Q2b2MvA1MyvxlpsPPHKlDW7atOmIme2LoEY/lQJH/C4ihlJ5/7RvySuV9+9q9m38QBYaUACYWQFwB/BnfYY/DXzbzLKADry+PbAKWAzsBtqBPwJwzrWY2VeBDd5yX3HOXXGmKudc0vSAzGyjc26m33XESirvn/YteaXy/sVj3wYUAM65dmD4BWNvADMusawDPnOZ9TwBPBF5mSIiEm26ElhEJE0pAKLncb8LiLFU3j/tW/JK5f2L+b5ZuGMjIiLpRr8ARETSlALgKplZhZn92sx2mNlbZvZXftcUbWaWaWa1Zvb//K4l2sys2MyeM7N3vP+Hc/2uKVrMbLn3Z3K7mT1jZnl+13Q1zOwJMztsZtv7jAXM7FVvgslX+5xmnlQus2+Pen8u3zSzGjMrjvZ2FQBXrwv4a+fc9cAc4DNmNtnnmqLtr4AdfhcRI98Gfu6cuw6YRorsp5mNAf4SmOmcmwJkEp7KJZn9iIvnD/sdAwxCAAACVklEQVQC8Evn3ETgl97rZPQjLt63V4EpzrmpwLv0c93UYCgArpJz7pBzbrP3vI3wB8hFcxwlKzMbCywBfuB3LdFmZoXAzXjTmDjnzjrnjvtbVVRlAfnetToFJPmV986514ELrx26C3jSe/4kcHdci4qSS+2bc+4V51yX93It4dkTokoBEEVmVglMB9b5W0lUfQv470Aq3lh2AtAM/KvX4vqBmQ3xu6hocM41AI8B+4FDhK/If8XfqmJipHPuEIS/jAEjfK4nVv4YeCnaK1UARImZDQX+E1jmnGv1u55oMLOPAoedc5v8riVGsoAbge8556YDp0jeFsJ5vF74XUAQGA0MMbNP+VuVDIaZ/Q/Creano71uBUAUmFk24Q//p51zK/tbPol8ELjTzOqAZ4HbzOzf/C0pquqBeudc7y+25wgHQir4CLDXOdfsnOskPHHjTf28Jxk1efcbwfvnYZ/riSrvzokfBT7pYnDOvgLgKpmZEe4h73DO/YPf9USTc+4R59xY51wl4QOIv3LOpcy3SOdcI3DAzCZ5Q7cDb/tYUjTtB+aYWYH3Z/R2UuQA9wVeIDwzMd4/n/exlqgys4XA54E7vel4ok4BcPU+CDxA+Ntx732TF/tdlAzYXwBPm9mbhG99+jWf64kK71fNc8BmYBvhv+tJfdWsmT0DrAEmmVm9mf0J8PfAHWa2i/CElX/vZ42DdZl9+y4wDHjV+1z5ftS3qyuBRUTSk34BiIikKQWAiEiaUgCIiKQpBYCISJpSAIiIpCkFgIhImlIAiIikKQWAiEia+v/tf99OV+SEDwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure()\n", "plt.plot(range(1,13), mses1);" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "# Вывод" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Как говорилось уже в самом описании к этим данным - такого количества недостатосно для описания полноценной линейной модели. Как мы видим она чаще ошибается чем оказывается права, что не является хорошим показателем данной модели. Однако, как мы видим, она выдает постоянный уровень ошибки при 8 и более признаках." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.0" } }, "nbformat": 4, "nbformat_minor": 2 }