2,2/2.2 homwork .ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Задание\n",
    "Построить модель машинного обучения, позволяющую предсказать, является ли цветок цветком класса Ирис-Вирджиния или нет.\n",
    "\n",
    "План решения:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "1. Из библиотеки scikit-learn подгрузите данные. Мы уже работали с ними в рамках Темы 1."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "\n",
    "from sklearn.model_selection import train_test_split\n",
    "from sklearn.datasets import load_iris\n",
    "from sklearn.metrics import ConfusionMatrixDisplay\n",
    "from sklearn.linear_model import LogisticRegression\n",
    "from sklearn.metrics import confusion_matrix"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>sepal length (cm)</th>\n",
       "      <th>sepal width (cm)</th>\n",
       "      <th>petal length (cm)</th>\n",
       "      <th>petal width (cm)</th>\n",
       "      <th>target</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>5.1</td>\n",
       "      <td>3.5</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>4.9</td>\n",
       "      <td>3.0</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>4.7</td>\n",
       "      <td>3.2</td>\n",
       "      <td>1.3</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4.6</td>\n",
       "      <td>3.1</td>\n",
       "      <td>1.5</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5.0</td>\n",
       "      <td>3.6</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   sepal length (cm)  sepal width (cm)  petal length (cm)  petal width (cm)  \\\n",
       "0                5.1               3.5                1.4               0.2   \n",
       "1                4.9               3.0                1.4               0.2   \n",
       "2                4.7               3.2                1.3               0.2   \n",
       "3                4.6               3.1                1.5               0.2   \n",
       "4                5.0               3.6                1.4               0.2   \n",
       "\n",
       "   target  \n",
       "0       0  \n",
       "1       0  \n",
       "2       0  \n",
       "3       0  \n",
       "4       0  "
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = load_iris(as_frame=True).frame\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Теперь создадим новое поле target_virginica и заполним его значениями, где 0 - цветок не является классом Ирис-Вирджиния, 1 - цветок принадлежит классу Ирис-Вирджиния:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>sepal length (cm)</th>\n",
       "      <th>sepal width (cm)</th>\n",
       "      <th>petal length (cm)</th>\n",
       "      <th>petal width (cm)</th>\n",
       "      <th>target</th>\n",
       "      <th>target_virginica</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>5.1</td>\n",
       "      <td>3.5</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>4.9</td>\n",
       "      <td>3.0</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>4.7</td>\n",
       "      <td>3.2</td>\n",
       "      <td>1.3</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4.6</td>\n",
       "      <td>3.1</td>\n",
       "      <td>1.5</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5.0</td>\n",
       "      <td>3.6</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   sepal length (cm)  sepal width (cm)  petal length (cm)  petal width (cm)  \\\n",
       "0                5.1               3.5                1.4               0.2   \n",
       "1                4.9               3.0                1.4               0.2   \n",
       "2                4.7               3.2                1.3               0.2   \n",
       "3                4.6               3.1                1.5               0.2   \n",
       "4                5.0               3.6                1.4               0.2   \n",
       "\n",
       "   target  target_virginica  \n",
       "0       0                 0  \n",
       "1       0                 0  \n",
       "2       0                 0  \n",
       "3       0                 0  \n",
       "4       0                 0  "
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['target_virginica'] = np.where(df['target'] == 2, 1, 0)\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "target_virginica  target\n",
       "0                 0         50\n",
       "                  1         50\n",
       "1                 2         50\n",
       "Name: count, dtype: int64"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(df.groupby('target_virginica')['target'].value_counts())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "2. Обучите библиотечную модель логистической регрессии на обучающей части данных."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [],
   "source": [
    "X = df.drop(['target', 'target_virginica'], axis=1)\n",
    "y = df['target_virginica']\n",
    "\n",
    "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42)\n",
    "\n",
    "model = LogisticRegression()\n",
    "model.fit(X_train, y_train)\n",
    "\n",
    "y_pred = model.predict(X_test)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "3. Оцените качество обученной на предыдущем шаге модели, выведя матрицу ошибок для прогнозов обученной моделью на тестовой части данных."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<sklearn.metrics._plot.confusion_matrix.ConfusionMatrixDisplay at 0x2069ffc2600>"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf8AAAGwCAYAAACn/2wHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAzlElEQVR4nO3deXRV5b3/8c9OIAlDEghDQiAiiEwKQbFNQ0XgJxrSLmSoQ7n0EhBw1YpDUxS8FQRR06tXUQoFJwi2UsBWI6KXFlGmgvYyxIpCSkIgQRIEGUKCJOGc/fsD2fGYBHKyz8lJzn6/1nrWunt49v4eL+XL832evbdhmqYpAADgGCGBDgAAADQskj8AAA5D8gcAwGFI/gAAOAzJHwAAhyH5AwDgMCR/AAAcplmgA7DD7XbryJEjioyMlGEYgQ4HAOAl0zR15swZxcfHKyTEf+PRc+fOqaKiwvZ1wsLCFBER4YOIAqtJJ/8jR44oISEh0GEAAGwqLCxUly5d/HLtc+fOqVvX1ir+ymX7WnFxccrPz2/y/wBo0sk/MjJSknRo15WKas0MBoLTmJ79Ah0C4DfnVamtet/6+9wfKioqVPyVS4d2XqmoyPrnipIzbnUdeFAVFRUk/0C6WOqPah1i6/+hQGPWzGge6BAA//n2BfMNMXXbOtJQ68j638et4JlebtLJHwCAunKZbrlsfM3GZbp9F0yAkfwBAI7glim36p/97fRtbKiVAwDgMIz8AQCO4JZbdgr39no3LiR/AIAjuExTLrP+pXs7fRsbyv4AADgMI38AgCOw4K8KyR8A4AhumXKR/CVR9gcAwHEY+QMAHIGyfxWSPwDAEVjtX4WyPwAADsPIHwDgCO5vm53+wYLkDwBwBJfN1f52+jY2JH8AgCO4TNn8qp/vYgk05vwBAHAYRv4AAEdgzr8KyR8A4AhuGXLJsNU/WFD2BwDAYRj5AwAcwW1eaHb6BwuSPwDAEVw2y/52+jY2lP0BAHAYRv4AAEdg5F+FkT8AwBHcpmG7eWPz5s0aOXKk4uPjZRiGsrKyPI4bhlFje/bZZ2u95pw5c6qd37t3b6//W5D8AQDwg7KyMiUmJmrRokU1Hi8qKvJoS5culWEY+tnPfnbJ615zzTUe/bZu3ep1bJT9AQCO4Kuyf0lJicf+8PBwhYeHVzs/NTVVqamptV4vLi7OY/udd97RsGHD1L1790vG0axZs2p9vcXIHwDgCC6F2G6SlJCQoOjoaKtlZGTYju3o0aN67733NHny5Mueu3//fsXHx6t79+4aP368CgoKvL4fI38AgCOY9Zi3/35/SSosLFRUVJS1v6ZRv7eWL1+uyMhIjR079pLnJSUlKTMzU7169VJRUZHmzp2rwYMHa8+ePYqMjKzz/Uj+AAB4ISoqyiP5+8LSpUs1fvx4RUREXPK8704j9O/fX0lJSeratatWr15dp6rBRSR/AIAjNNZH/bZs2aKcnBytWrXK675t2rRRz549lZub61U/5vwBAI7gMkNsN3947bXXNHDgQCUmJnrdt7S0VHl5eerUqZNX/Uj+AAD4QWlpqbKzs5WdnS1Jys/PV3Z2tscCvZKSEr355puaMmVKjde4+eabtXDhQmt7+vTp2rRpkw4ePKht27ZpzJgxCg0N1bhx47yKjbI/AMAR3DLktjHmdcu7L/vs2LFDw4YNs7bT09MlSWlpacrMzJQkrVy5UqZp1pq88/LydPz4cWv78OHDGjdunL7++mt16NBBN954oz7++GN16NDBq9gM0zSb7HeKSkpKFB0drZP/7q6oSIoYCE4p8QMCHQLgN+fNSm3UOzp9+rTPF9FddDFXrPnXVWoVGVrv65Sdcem2/nl+jbWhkDEBAHAYyv4AAEewu2jP1XQL5dWQ/AEAjnBhzr/+j+vZ6dvYUPYHAMBhGPkDABzB/Z3389evP2V/AACaFOb8q5D8AQCO4FZIgz7n35gx5w8AgMMw8gcAOILLNOSy8UlfO30bG5I/AMARXDYX/Lko+wMAgKaKkT8AwBHcZojcNlb7u1ntDwBA00LZvwplfwAAHIaRPwDAEdyyt2Lf7btQAo7kDwBwBPsv+QmeYnnw/BIAAFAnjPwBAI5g/93+wTNeJvkDABzBLUNu2Znz5w1/AAA0KYz8qwTPLwEAAHXCyB8A4Aj2X/ITPONlkj8AwBHcpiG3nef8g+irfsHzzxgAAFAnjPwBAI7gtln2D6aX/JD8AQCOYP+rfsGT/IPnlwAAgDph5A8AcASXDLlsvKjHTt/GhuQPAHAEyv5VgueXAACAOmHkDwBwBJfsle5dvgsl4Ej+AABHoOxfheQPAHAEPuxTJXh+CQAAqBNG/gAARzBlyG1jzt/kUT8AAJoWyv5VgueXAACAOiH5AwAc4eInfe00b2zevFkjR45UfHy8DMNQVlaWx/GJEyfKMAyPNmLEiMted9GiRbryyisVERGhpKQk/fOf//QqLonkDwBwCNe3X/Wz07xRVlamxMRELVq0qNZzRowYoaKiIqv9+c9/vuQ1V61apfT0dD3++OPatWuXEhMTlZKSoq+++sqr2JjzBwDAD1JTU5WamnrJc8LDwxUXF1fnaz7//POaOnWqJk2aJElasmSJ3nvvPS1dulQzZ86s83UY+QMAHMFXZf+SkhKPVl5eXu+YNm7cqI4dO6pXr16699579fXXX9d6bkVFhXbu3Knhw4db+0JCQjR8+HBt377dq/uS/AEAjuBWiO0mSQkJCYqOjrZaRkZGveIZMWKEXn/9dW3YsEH//d//rU2bNik1NVUuV80vEj5+/LhcLpdiY2M99sfGxqq4uNire1P2BwDAC4WFhYqKirK2w8PD63Wdn//859b/3a9fP/Xv319XXXWVNm7cqJtvvtl2nJfCyB8A4Agu07DdJCkqKsqj1Tf5f1/37t3Vvn175ebm1ni8ffv2Cg0N1dGjRz32Hz161Kt1AxLJHwDgEA39qJ+3Dh8+rK+//lqdOnWq8XhYWJgGDhyoDRs2VP0mt1sbNmxQcnKyV/ei7A8AcATT5lf9TC/7lpaWeozi8/PzlZ2drZiYGMXExGju3Ln62c9+pri4OOXl5emRRx5Rjx49lJKSYvW5+eabNWbMGE2bNk2SlJ6errS0NN1www364Q9/qBdeeEFlZWXW6v+6IvkDAOAHO3bs0LBhw6zt9PR0SVJaWpoWL16sf/3rX1q+fLlOnTql+Ph43XrrrZo3b57HNEJeXp6OHz9ubd911106duyYZs+ereLiYg0YMEDr1q2rtgjwckj+AABHcMmQy8bHebztO3ToUJmmWevxv/3tb5e9xsGDB6vtmzZtmlUJqC+SPwDAEdymbM3bu2vP400OC/4AAHAYRv6o5rOPW+nNP3TU/s9a6sTR5nr8tXwNSj1tHT95rJleeypeOzdFqux0qK79Uanue/KwOnevCGDUgH0jJx7X7fd+pZgO53Xgixb6w2OdlZPdMtBhwUfcNhf82enb2ATPL4HPnDsbou7XfKNpTx+udsw0pbl3d1PRoTDNWXZAi/6eo9guFZp5Vw+dO8sfJzRdQ247qXseP6I3no/TfSk9deCLCD214oCi21UGOjT4iFuG7RYsGsXf1r74PCF85wf/74wmzijWj78z2r/oywPh2ruzle7/3WH1GvCNEnqU6/7fHVb5OUMfvd2m4YMFfGTsPce1bkWM/r4qRgX7I7RgRheVf2MoZdyJQIcG+FzAk7+vPk+IhlFZceFfvmHhbmtfSIjUPMzU5//XOlBhAbY0a+7W1f3PateWSGufaRravSVSfQeeDWBk8CVfveEvGAQ8+X/384R9+/bVkiVL1LJlSy1dujTQoaEGCT3OqWPnCi3N6KQzp0JVWWFo1cKOOl4UphNHWUKCpikqxqXQZtKpY55/hk8eb6a2Hc4HKCr42sU5fzstWAT0l3j7ecLy8vJqn1JEw2rWXJr9Wr6+zIvQ7X376bar+uvTba31g/9XIiN4/ncBAEEtoEO1S32ecN++fdXOz8jI0Ny5cxsqPNTi6v7faPEHOSorCVFlpaE27Vx64KdXq2d/yqNomkpOhMp1XmrzvVF+2/bndfIYFa1g4Za99/Oz4C9AHn30UZ0+fdpqhYWFgQ7J0VpFudWmnUtfHgjT/k9bKjmFSgyapvOVIdr/r5a67sYz1j7DMDXgxlJ9sZNH/YKFaXOlvxlEyT+g/6T19vOE4eHhPvt0Imr3TVmIjuRX/XcuLgxT3p4WimxzXh27VGrzu9GKbudSx84Vyt8boSWzuyh5xGkNHHrmElcFGre3Xm6v6S8U6t+ftlTO7pYaM/WYIlq69feVMYEODT5i98t8/v6qX0MKaPL/7ucJR48eLanq84R231uM+vv3py31yO09rO2X5nSWJN1y5wlNf6FAJ44210tzOuvU8WaK6Xhew+84of946GhtlwOahE1r2iq6nUsTHi5W2w7ndeDzFvrt+G46dbx5oEMDfC7gk1m++jwhfCdxUKn+diS71uOjpxzX6CnHaz0ONFVrlrXXmmXtAx0G/IQ3/FUJePL31ecJAQC4FMr+VQKe/CXffJ4QAADUTaNI/gAA+Jvd9/MH06N+JH8AgCNQ9q8SPKsXAABAnTDyBwA4AiP/KiR/AIAjkPyrUPYHAMBhGPkDAByBkX8Vkj8AwBFM2Xtcz/RdKAFH8gcAOAIj/yrM+QMA4DCM/AEAjsDIvwrJHwDgCCT/KpT9AQBwGEb+AABHYORfheQPAHAE0zRk2kjgdvo2NpT9AQBwGEb+AABHcMuw9ZIfO30bG5I/AMARmPOvQtkfAACHYeQPAHAEFvxVIfkDAByBsn8Vkj8AwBEY+Vdhzh8AAIch+QMAHMH8tuxf3+btyH/z5s0aOXKk4uPjZRiGsrKyrGOVlZWaMWOG+vXrp1atWik+Pl4TJkzQkSNHLnnNOXPmyDAMj9a7d2+v/1uQ/AEAjmBKMk0bzcv7lZWVKTExUYsWLap27OzZs9q1a5dmzZqlXbt26a233lJOTo5uu+22y173mmuuUVFRkdW2bt3qZWTM+QMA4BepqalKTU2t8Vh0dLTWr1/vsW/hwoX64Q9/qIKCAl1xxRW1XrdZs2aKi4uzFRsjfwCAI1x8w5+dJkklJSUerby83CfxnT59WoZhqE2bNpc8b//+/YqPj1f37t01fvx4FRQUeH0vkj8AwBEurva30yQpISFB0dHRVsvIyLAd27lz5zRjxgyNGzdOUVFRtZ6XlJSkzMxMrVu3TosXL1Z+fr4GDx6sM2fOeHU/yv4AAHihsLDQI0GHh4fbul5lZaXuvPNOmaapxYsXX/Lc704j9O/fX0lJSeratatWr16tyZMn1/meJH8AgCO4TUOGD17yExUVdcnRuTcuJv5Dhw7pww8/9Pq6bdq0Uc+ePZWbm+tVP8r+AABHsLXS/9vmSxcT//79+/XBBx+oXbt2Xl+jtLRUeXl56tSpk1f9SP4AAPhBaWmpsrOzlZ2dLUnKz89Xdna2CgoKVFlZqdtvv107duzQG2+8IZfLpeLiYhUXF6uiosK6xs0336yFCxda29OnT9emTZt08OBBbdu2TWPGjFFoaKjGjRvnVWyU/QEAjtDQr/fdsWOHhg0bZm2np6dLktLS0jRnzhytWbNGkjRgwACPfh999JGGDh0qScrLy9Px48etY4cPH9a4ceP09ddfq0OHDrrxxhv18ccfq0OHDl7FRvIHADhCQyf/oUOHyrzEXMGljl108OBBj+2VK1d6FUNtSP4AAEfw1YK/YMCcPwAADsPIHwDgCHZX7Pt6tX8gkfwBAI5wIfnbmfP3YTABRtkfAACHYeQPAHCEhl7t35iR/AEAjmB+2+z0DxaU/QEAcBhG/gAAR6DsX4XkDwBwBur+FpI/AMAZbI78FUQjf+b8AQBwGEb+AABH4A1/VUj+AABHYMFfFcr+AAA4DCN/AIAzmIa9RXtBNPIn+QMAHIE5/yqU/QEAcBhG/gAAZ+AlPxaSPwDAEVjtX6VOyX/NmjV1vuBtt91W72AAAID/1Sn5jx49uk4XMwxDLpfLTjwAAPhPEJXu7ahT8ne73f6OAwAAv6LsX8XWav9z5875Kg4AAPzL9EELEl4nf5fLpXnz5qlz585q3bq1Dhw4IEmaNWuWXnvtNZ8HCAAAfMvr5P/UU08pMzNTzzzzjMLCwqz91157rV599VWfBgcAgO8YPmjBwevk//rrr+vll1/W+PHjFRoaau1PTEzUvn37fBocAAA+Q9nf4nXy//LLL9WjR49q+91utyorK30SFAAA8B+vk3/fvn21ZcuWavv/8pe/6LrrrvNJUAAA+Bwjf4vXb/ibPXu20tLS9OWXX8rtduutt95STk6OXn/9da1du9YfMQIAYB9f9bN4PfIfNWqU3n33XX3wwQdq1aqVZs+erb179+rdd9/VLbfc4o8YAQCAD9Xr3f6DBw/W+vXrfR0LAAB+wyd9q9T7wz47duzQ3r17JV1YBzBw4ECfBQUAgM/xVT+L18n/8OHDGjdunP7xj3+oTZs2kqRTp05p0KBBWrlypbp06eLrGAEAgA95Pec/ZcoUVVZWau/evTpx4oROnDihvXv3yu12a8qUKf6IEQAA+y4u+LPTgoTXI/9NmzZp27Zt6tWrl7WvV69e+v3vf6/Bgwf7NDgAAHzFMC80O/2DhdfJPyEhocaX+bhcLsXHx/skKAAAfI45f4vXZf9nn31W999/v3bs2GHt27Fjhx588EH9z//8j0+DAwAAvlen5N+2bVvFxMQoJiZGkyZNUnZ2tpKSkhQeHq7w8HAlJSVp165duvvuu/0dLwAA9dPAc/6bN2/WyJEjFR8fL8MwlJWV5RmOaWr27Nnq1KmTWrRooeHDh2v//v2Xve6iRYt05ZVXKiIiQklJSfrnP//pVVxSHcv+L7zwgtcXBgCgUWngsn9ZWZkSExN19913a+zYsdWOP/PMM1qwYIGWL1+ubt26adasWUpJSdEXX3yhiIiIGq+5atUqpaena8mSJUpKStILL7yglJQU5eTkqGPHjnWOrU7JPy0trc4XBAAgmJWUlHhsX6yCf19qaqpSU1NrvIZpmnrhhRf02GOPadSoUZIufDU3NjZWWVlZ+vnPf15jv+eff15Tp07VpEmTJElLlizRe++9p6VLl2rmzJl1/g1ez/l/17lz51RSUuLRAABolHz0YZ+EhARFR0dbLSMjw+tQ8vPzVVxcrOHDh1v7oqOjlZSUpO3bt9fYp6KiQjt37vToExISouHDh9fapzZer/YvKyvTjBkztHr1an399dfVjrtcLm8vCQCA//mo7F9YWKioqChrd02j/sspLi6WJMXGxnrsj42NtY593/Hjx+VyuWrss2/fPq/u7/XI/5FHHtGHH36oxYsXKzw8XK+++qrmzp2r+Ph4vf76695eDgCAJiUqKsqj1Sf5B5rXyf/dd9/VH/7wB/3sZz9Ts2bNNHjwYD322GN6+umn9cYbb/gjRgAA7GtEb/iLi4uTJB09etRj/9GjR61j39e+fXuFhoZ61ac2Xif/EydOqHv37pIu/OvnxIkTkqQbb7xRmzdv9vZyAAA0iItv+LPTfKVbt26Ki4vThg0brH0lJSX65JNPlJycXGOfsLAwDRw40KOP2+3Whg0bau1TG6+Tf/fu3ZWfny9J6t27t1avXi3pQkXg4od+AABwutLSUmVnZys7O1vShUV+2dnZKigokGEYeuihh/Tkk09qzZo1+uyzzzRhwgTFx8dr9OjR1jVuvvlmLVy40NpOT0/XK6+8ouXLl2vv3r269957VVZWZq3+ryuvF/xNmjRJn376qYYMGaKZM2dq5MiRWrhwoSorK/X88897ezkAABpGAz/nv2PHDg0bNszaTk9Pl3Th8fnMzEw98sgjKisr0z333KNTp07pxhtv1Lp16zye8c/Ly9Px48et7bvuukvHjh3T7NmzVVxcrAEDBmjdunXVFgFejmGapq1CxqFDh7Rz50716NFD/fv3t3Mpr5WUlCg6Olon/91dUZG2nloEGq2U+AGBDgHwm/NmpTbqHZ0+fdpjBb0vXcwVV/z3kwppUfPLc+rC/c05Fcx4zK+xNhSvR/7f17VrV3Xt2tUXsQAA4DeGbH7Vz2eRBF6dkv+CBQvqfMEHHnig3sEAAAD/q1Pynz9/fp0uZhhGQJL/mJ791Mxo3uD3BRrC/tevD3QIgN+4vzkn3fNOw9zM7uN6PnzUL9DqlPwvru4HAKDJauAFf40Zq+QAAHAY2wv+AABoEhj5W0j+AABHsPuWPl++4S/QKPsDAOAwjPwBAM5A2d9Sr5H/li1b9Itf/ELJycn68ssvJUl//OMftXXrVp8GBwCAz5g+aEHC6+T/17/+VSkpKWrRooV2796t8vJySdLp06f19NNP+zxAAADgW14n/yeffFJLlizRK6+8oubNq16s8+Mf/1i7du3yaXAAAPhKY/qkb6B5Peefk5Ojm266qdr+6OhonTp1yhcxAQDge7zhz+L1yD8uLk65ubnV9m/dulXdu3f3SVAAAPgcc/4Wr5P/1KlT9eCDD+qTTz6RYRg6cuSI3njjDU2fPl333nuvP2IEAAA+5HXZf+bMmXK73br55pt19uxZ3XTTTQoPD9f06dN1//33+yNGAABs4yU/VbxO/oZh6Le//a0efvhh5ebmqrS0VH379lXr1q39ER8AAL7Bc/6Wer/kJywsTH379vVlLAAAoAF4nfyHDRsmw6h9xeOHH35oKyAAAPzC7uN6Th75DxgwwGO7srJS2dnZ2rNnj9LS0nwVFwAAvkXZ3+J18p8/f36N++fMmaPS0lLbAQEAAP/y2Vf9fvGLX2jp0qW+uhwAAL7Fc/4Wn33Vb/v27YqIiPDV5QAA8Cke9avidfIfO3asx7ZpmioqKtKOHTs0a9YsnwUGAAD8w+vkHx0d7bEdEhKiXr166YknntCtt97qs8AAAIB/eJX8XS6XJk2apH79+qlt27b+igkAAN9jtb/FqwV/oaGhuvXWW/l6HwCgyeGTvlW8Xu1/7bXX6sCBA/6IBQAANACvk/+TTz6p6dOna+3atSoqKlJJSYlHAwCg0eIxP0lezPk/8cQT+s1vfqOf/OQnkqTbbrvN4zW/pmnKMAy5XC7fRwkAgF3M+VvqnPznzp2rX/7yl/roo4/8GQ8AAPCzOid/07zwT54hQ4b4LRgAAPyFl/xU8epRv0t9zQ8AgEaNsr/Fq+Tfs2fPy/4D4MSJE7YCAgAA/uVV8p87d261N/wBANAUUPav4lXy//nPf66OHTv6KxYAAPyHsr+lzs/5M98PAEBwqHPyv7jaHwCAJsnOC37qUTW48sorZRhGtXbffffVeH5mZma1cyMiIurxQy+vzmV/t9vtlwAAAGgIDT3n/3//938eL77bs2ePbrnlFt1xxx219omKilJOTk7VPf1Udff6k74AADRJDTzn36FDB4/t3/3ud7rqqqsu+b4cwzAUFxdXn+i84vW7/QEAcLLvf9OmvLz8sn0qKir0pz/9SXffffclR/OlpaXq2rWrEhISNGrUKH3++ee+DN1C8gcAOIOP5vwTEhIUHR1ttYyMjMveOisrS6dOndLEiRNrPadXr15aunSp3nnnHf3pT3+S2+3WoEGDdPjw4Xr+4NpR9gcAOIKv5vwLCwsVFRVl7Q8PD79s39dee02pqamKj4+v9Zzk5GQlJydb24MGDVKfPn300ksvad68efUPvAYkfwAAvBAVFeWR/C/n0KFD+uCDD/TWW295dZ/mzZvruuuuU25urrchXhZlfwCAMzTwo34XLVu2TB07dtRPf/pTr/q5XC599tln6tSpU/1ufAmM/AEAjhCI1/u63W4tW7ZMaWlpatbMM+VOmDBBnTt3ttYMPPHEE/rRj36kHj166NSpU3r22Wd16NAhTZkypf5B14LkDwCAn3zwwQcqKCjQ3XffXe1YQUGBQkKqCvAnT57U1KlTVVxcrLZt22rgwIHatm2b+vbt6/O4SP4AAGcIwLv9b7311lrfkLtx40aP7fnz52v+/Pn1CMx7JH8AgDPwYR8LC/4AAHAYRv4AAEcwvm12+gcLkj8AwBko+1tI/gAARwjEo36NFXP+AAA4DCN/AIAzUPa3kPwBAM4RRAncDsr+AAA4DCN/AIAjsOCvCskfAOAMzPlbKPsDAOAwjPwBAI5A2b8KyR8A4AyU/S2U/QEAcBhG/gAAR6DsX4XkDwBwBsr+FpI/AMAZSP4W5vwBAHAYRv4AAEdgzr8KyR8A4AyU/S2U/QEAcBhG/gAARzBMU4ZZ/+G7nb6NDckfAOAMlP0tlP0BAHAYRv4AAEdgtX8Vkj8AwBko+1so+wMA4DCM/AEAjkDZvwrJHwDgDJT9LSR/AIAjMPKvwpw/AAAOw8gfAOAMlP0tJH8AgGMEU+neDsr+AAA4DCN/AIAzmOaFZqd/kCD5AwAcgdX+VSj7AwDgMCR/AIAzmD5oXpgzZ44Mw/BovXv3vmSfN998U71791ZERIT69eun999/37ub1hHJHwDgCIbbfvPWNddco6KiIqtt3bq11nO3bdumcePGafLkydq9e7dGjx6t0aNHa8+ePTZ+dc2Y8wcAwAslJSUe2+Hh4QoPD6/x3GbNmikuLq5O133xxRc1YsQIPfzww5KkefPmaf369Vq4cKGWLFliL+jvYeSPOhs58biWf/KF3j3wL724dr96DTgb6JCAeonYd0adns9Vtwc+09UTdqnVzlMex1v930nFP7Nf3e/9VFdP2KWwQ/xZDwo+KvsnJCQoOjraahkZGbXecv/+/YqPj1f37t01fvx4FRQU1Hru9u3bNXz4cI99KSkp2r59e71+7qWQ/FEnQ247qXseP6I3no/TfSk9deCLCD214oCi21UGOjTAayHlblVc0VJfTUio+XiFW+d6ttbxuzo3cGTwp4ur/e00SSosLNTp06et9uijj9Z4v6SkJGVmZmrdunVavHix8vPzNXjwYJ05c6bG84uLixUbG+uxLzY2VsXFxT797yAFOPlv3rxZI0eOVHx8vAzDUFZWViDDwSWMvee41q2I0d9Xxahgf4QWzOii8m8MpYw7EejQAK+dTYzW17fHq+yGNjUeP/PjdjoxupPOXhPZsIHBvy4+52+nSYqKivJotZX8U1NTdccdd6h///5KSUnR+++/r1OnTmn16tUN+atrFNDkX1ZWpsTERC1atCiQYeAymjV36+r+Z7VrS9VfhKZpaPeWSPUdSDkUAOqiTZs26tmzp3Jzc2s8HhcXp6NHj3rsO3r0aJ3XDHgjoMk/NTVVTz75pMaMGVOn88vLy1VSUuLR4H9RMS6FNpNOHfNcH3ryeDO17XA+QFEBgHd8Vfavr9LSUuXl5alTp041Hk9OTtaGDRs89q1fv17Jycn2blyDJjXnn5GR4bHIIiGh5vk6AACqaeDn/KdPn65Nmzbp4MGD2rZtm8aMGaPQ0FCNGzdOkjRhwgSP9QIPPvig1q1bp+eee0779u3TnDlztGPHDk2bNs3Or65Rk0r+jz76qMcii8LCwkCH5AglJ0LlOi+1+d4ov2378zp5jKdFAaAmhw8f1rhx49SrVy/deeedateunT7++GN16NBBklRQUKCioiLr/EGDBmnFihV6+eWXlZiYqL/85S/KysrStdde6/PYmtTf3Jd6lhL+c74yRPv/1VLX3XhG29dFS5IMw9SAG0u1JrNdgKMDgLpp6Hf7r1y58pLHN27cWG3fHXfcoTvuuMO7G9VDk0r+CJy3Xm6v6S8U6t+ftlTO7pYaM/WYIlq69feVMYEODfCacc6l5kfLre3mx8oVduis3K2a6Xz7MIWUnlezryvU7NSFR1nDis5JklzRzeVq0zwgMcMH+KqfheSPOtm0pq2i27k04eFite1wXgc+b6Hfju+mU8f5ixBNT0T+WXXJ2G9td1jxpSSp5MYYHb3nSrXafVpxrxyyjnf6w0FJ0tej43RibHyDxgr4Q0CTf2lpqccjD/n5+crOzlZMTIyuuOKKAEaGmqxZ1l5rlrUPdBiAbd/0idT+16+v9fiZwe10ZjBTWsGGT/pWCWjy37Fjh4YNG2Ztp6enS5LS0tKUmZkZoKgAAEGpHiv2q/UPEgFN/kOHDpUZRHMoAAA0Bcz5AwAcgbJ/FZI/AMAZ3OaFZqd/kCD5AwCcgTl/S5N6wx8AALCPkT8AwBEM2Zzz91kkgUfyBwA4A2/4s1D2BwDAYRj5AwAcgUf9qpD8AQDOwGp/C2V/AAAchpE/AMARDNOUYWPRnp2+jQ3JHwDgDO5vm53+QYKyPwAADsPIHwDgCJT9q5D8AQDOwGp/C8kfAOAMvOHPwpw/AAAOw8gfAOAIvOGvCskfAOAMlP0tlP0BAHAYRv4AAEcw3Beanf7BguQPAHAGyv4Wyv4AADgMI38AgDPwkh8LyR8A4Ai83rcKZX8AAByGkT8AwBlY8Gch+QMAnMGUZOdxveDJ/SR/AIAzMOdfhTl/AAAchpE/AMAZTNmc8/dZJAFH8gcAOAML/iyU/QEA8IOMjAz94Ac/UGRkpDp27KjRo0crJyfnkn0yMzNlGIZHi4iI8HlsJH8AgDO4fdC8sGnTJt133336+OOPtX79elVWVurWW29VWVnZJftFRUWpqKjIaocOHfLuxnVA2R8A4AgNvdp/3bp1HtuZmZnq2LGjdu7cqZtuuqn2+xiG4uLi6hVjXTHyBwDACyUlJR6tvLy8Tv1Onz4tSYqJibnkeaWlperatasSEhI0atQoff7557Zj/j6SPwDAGS4u+LPTJCUkJCg6OtpqGRkZl7212+3WQw89pB//+Me69tpraz2vV69eWrp0qd555x396U9/ktvt1qBBg3T48GGf/WeQKPsDAJzCR6v9CwsLFRUVZe0ODw+/bNf77rtPe/bs0datWy95XnJyspKTk63tQYMGqU+fPnrppZc0b968egZeHckfAAAvREVFeST/y5k2bZrWrl2rzZs3q0uXLl7dq3nz5rruuuuUm5vrbZiXRNkfAOAMPir71/12pqZNm6a3335bH374obp16+Z1yC6XS5999pk6derkdd9LYeQPAHAGtyTDZn8v3HfffVqxYoXeeecdRUZGqri4WJIUHR2tFi1aSJImTJigzp07W+sGnnjiCf3oRz9Sjx49dOrUKT377LM6dOiQpkyZYiPw6kj+AABHaOhH/RYvXixJGjp0qMf+ZcuWaeLEiZKkgoIChYRUFeFPnjypqVOnqri4WG3bttXAgQO1bds29e3bt95x14TkDwCAH5h1+MfCxo0bPbbnz5+v+fPn+ymiKiR/AIAz8G5/C8kfAOAMblMybCRwd/Akf1b7AwDgMIz8AQDOQNnfQvIHADiEzeSv4En+lP0BAHAYRv4AAGeg7G8h+QMAnMFtylbpntX+AACgqWLkDwBwBtN9odnpHyRI/gAAZ2DO30LyBwA4A3P+Fub8AQBwGEb+AABnoOxvIfkDAJzBlM3k77NIAo6yPwAADsPIHwDgDJT9LSR/AIAzuN2SbDyr7w6e5/wp+wMA4DCM/AEAzkDZ30LyBwA4A8nfQtkfAACHYeQPAHAGXu9rIfkDABzBNN0ybXyZz07fxobkDwBwBtO0N3pnzh8AADRVjPwBAM5g2pzzD6KRP8kfAOAMbrdk2Ji3D6I5f8r+AAA4DCN/AIAzUPa3kPwBAI5gut0ybZT9g+lRP8r+AAA4DCN/AIAzUPa3kPwBAM7gNiWD5C9R9gcAwHEY+QMAnME0Jdl5zj94Rv4kfwCAI5huU6aNsr9J8gcAoIkx3bI38udRPwAAUAeLFi3SlVdeqYiICCUlJemf//znJc9/88031bt3b0VERKhfv356//33fR4TyR8A4Aim27TdvLVq1Sqlp6fr8ccf165du5SYmKiUlBR99dVXNZ6/bds2jRs3TpMnT9bu3bs1evRojR49Wnv27LH78z2Q/AEAzmC67TcvPf/885o6daomTZqkvn37asmSJWrZsqWWLl1a4/kvvviiRowYoYcfflh9+vTRvHnzdP3112vhwoV2f72HJj3nf3HxxXlV2npvA9CYub85F+gQAL9xf1MuqWEW09nNFedVKUkqKSnx2B8eHq7w8PBq51dUVGjnzp169NFHrX0hISEaPny4tm/fXuM9tm/frvT0dI99KSkpysrKqn/gNWjSyf/MmTOSpK3y/XwI0Gjc806gIwD87syZM4qOjvbLtcPCwhQXF6etxfZzRevWrZWQkOCx7/HHH9ecOXOqnXv8+HG5XC7FxsZ67I+NjdW+fftqvH5xcXGN5xcXF9sL/HuadPKPj49XYWGhIiMjZRhGoMNxhJKSEiUkJKiwsFBRUVGBDgfwKf58NzzTNHXmzBnFx8f77R4RERHKz89XRUWF7WuZplkt39Q06m/smnTyDwkJUZcuXQIdhiNFRUXxlyOCFn++G5a/RvzfFRERoYiICL/f57vat2+v0NBQHT161GP/0aNHFRcXV2OfuLg4r86vLxb8AQDgB2FhYRo4cKA2bNhg7XO73dqwYYOSk5Nr7JOcnOxxviStX7++1vPrq0mP/AEAaMzS09OVlpamG264QT/84Q/1wgsvqKysTJMmTZIkTZgwQZ07d1ZGRoYk6cEHH9SQIUP03HPP6ac//alWrlypHTt26OWXX/ZpXCR/eCU8PFyPP/54k5zjAi6HP9/wtbvuukvHjh3T7NmzVVxcrAEDBmjdunXWor6CggKFhFQV4QcNGqQVK1boscce03/913/p6quvVlZWlq699lqfxmWYwfSyYgAAcFnM+QMA4DAkfwAAHIbkDwCAw5D8AQBwGJI/6szbz1ICTcXmzZs1cuRIxcfHyzAMn79HHWhsSP6oE28/Swk0JWVlZUpMTNSiRYsCHQrQIHjUD3WSlJSkH/zgB9ZnJd1utxISEnT//fdr5syZAY4O8B3DMPT2229r9OjRgQ4F8BtG/risi5+lHD58uLXvcp+lBAA0XiR/XNalPkvp689MAgD8j+QPAIDDkPxxWfX5LCUAoPEi+eOy6vNZSgBA48VX/VAnl/ssJdCUlZaWKjc319rOz89Xdna2YmJidMUVVwQwMsA/eNQPdbZw4UI9++yz1mcpFyxYoKSkpECHBdi2ceNGDRs2rNr+tLQ0ZWZmNnxAgJ+R/AEAcBjm/AEAcBiSPwAADkPyBwDAYUj+AAA4DMkfAACHIfkDAOAwJH8AAByG5A8AgMOQ/AGbJk6cqNGjR1vbQ4cO1UMPPdTgcWzcuFGGYejUqVO1nmMYhrKysup8zTlz5mjAgAG24jp48KAMw1B2drat6wDwHZI/gtLEiRNlGIYMw1BYWJh69OihJ554QufPn/f7vd966y3NmzevTufWJWEDgK/xYR8ErREjRmjZsmUqLy/X+++/r/vuu0/NmzfXo48+Wu3ciooKhYWF+eS+MTExPrkOAPgLI38ErfDwcMXFxalr16669957NXz4cK1Zs0ZSVan+qaeeUnx8vHr16iVJKiws1J133qk2bdooJiZGo0aN0sGDB61rulwupaenq02bNmrXrp0eeeQRff/zGN8v+5eXl2vGjBlKSEhQeHi4evTooddee00HDx60PibTtm1bGYahiRMnSrrwyeSMjAx169ZNLVq0UGJiov7yl7943Of9999Xz5491aJFCw0bNswjzrqaMWOGevbsqZYtW6p79+6aNWuWKisrq5330ksvKSEhQS1bttSdd96p06dPexx/9dVX1adPH0VERKh37976wx/+4HUsABoOyR+O0aJFC1VUVFjbGzZsUE5OjtavX6+1a9eqsrJSKSkpioyM1JYtW/SPf/xDrVu31ogRI6x+zz33nDIzM7V06VJt3bpVJ06c0Ntvv33J+06YMEF//vOftWDBAu3du1cvvfSSWrdurYSEBP31r3+VJOXk5KioqEgvvviiJCkjI0Ovv/66lixZos8//1y//vWv9Ytf/EKbNm2SdOEfKWPHjtXIkSOVnZ2tKVOmaObMmV7/N4mMjFRmZqa++OILvfjii3rllVc0f/58j3Nyc3O1evVqvfvuu1q3bp12796tX/3qV9bxN954Q7Nnz9ZTTz2lvXv36umnn9asWbO0fPlyr+MB0EBMIAilpaWZo0aNMk3TNN1ut7l+/XozPDzcnD59unU8NjbWLC8vt/r88Y9/NHv16mW63W5rX3l5udmiRQvzb3/7m2maptmpUyfzmWeesY5XVlaaXbp0se5lmqY5ZMgQ88EHHzRN0zRzcnJMSeb69etrjPOjjz4yJZknT5609p07d85s2bKluW3bNo9zJ0+ebI4bN840TdN89NFHzb59+3ocnzFjRrVrfZ8k8+233671+LPPPmsOHDjQ2n788cfN0NBQ8/Dhw9a+//3f/zVDQkLMoqIi0zRN86qrrjJXrFjhcZ158+aZycnJpmmaZn5+vinJ3L17d633BdCwmPNH0Fq7dq1at26tyspKud1u/cd//IfmzJljHe/Xr5/HPP+nn36q3NxcRUZGelzn3LlzysvL0+nTp1VUVKSkpCTrWLNmzXTDDTdUK/1flJ2drdDQUA0ZMqTOcefm5urs2bO65ZZbPPZXVFTouuuukyTt3bvXIw5JSk5OrvM9Llq1apUWLFigvLw8lZaW6vz584qKivI454orrlDnzp097uN2u5WTk6PIyEjl5eVp8uTJmjp1qnXO+fPnFR0d7XU8ABoGyR9Ba9iwYVq8eLHCwsIUHx+vZs08/7i3atXKY7u0tFQDBw7UG2+8Ue1aHTp0qFcMLVq08LpPaWmpJOm9997zSLrShXUMvrJ9+3aNHz9ec+fOVUpKiqKjo7Vy5Uo999xzXsf6yiuvVPvHSGhoqM9iBeBbJH8ErVatWqlHjx51Pv/666/XqlWr1LFjx2qj34s6deqkTz75RDfddJOkCyPcnTt36vrrr6/x/H79+sntdmvTpk0aPnx4teMXKw8ul8va17dvX4WHh6ugoKDWikGfPn2sxYsXffzxx5f/kd+xbds2de3aVb/97W+tfYcOHap2XkFBgY4cOaL4+HjrPiEhIerVq5diY2MVHx+vAwcOaPz48V7dH0DgsOAP+Nb48ePVvn17jRo1Slu2bFF+fr42btyoBx54QIcPH5YkPfjgg/rd736nrKws7du3T7/61a8u+Yz+lVdeqbS0NN19993Kysqyrrl69WpJUteuXWUYhtauXatjx46ptLRUkZGRmj59un79619r+fLlysvL065du/T73//eWkT3y1/+Uvv379fDDz+snJwcrVixQpmZmV793quvvloFBQVauXKl8vLytGDBghoXL0ZERCgtLU2ffvqptmzZogceeEB33nmn4uLiJElz585VRkaGFixYoH//+9/67LPPtGzZMj3//PNexQOg4ZD8gW+1bNlSmzdv1hVXXKGxY8eqT58+mjx5ss6dO2dVAn7zm9/oP//zP5WWlqbk5GRFRkZqzJgxl7zu4sWLdfvtt+tXv/qVevfuralTp6qsrEyS1LlzZ82dO1czZ85UbGyspk2bJkmaN2+eZs2apYyMDPXp00cjRozQe++9p27dukm6MA//17/+VVlZWUpMTNSSJUv09NNPe/V7b7vtNv3617/WtGnTNGDAAG3btk2zZs2qdl6PHj00duxY/eQnP9Gtt96q/v37ezzKN2XKFL366qtatmyZ+vXrpyFDhigzM9OKFUDjY5i1rVQCAABBiZE/AAAOQ/IHAMBhSP4AADgMyR8AAIch+QMA4DAkfwAAHIbkDwCAw5D8AQBwGJI/AAAOQ/IHAMBhSP4AADjM/weaS5sLoNazdAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ConfusionMatrixDisplay(confusion_matrix(y_test, y_pred)).plot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Вывод:\n",
    "\n",
    "Обученная модель прошла оценку по матрице ошибок со следующими показателями качества:\n",
    "- полнота = 1\n",
    "- точность = 1\n",
    "- F1-мера = 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "4. Сконструируйте свой класс MyLogisticRegression, решающий задачу логистической регрессии методом градиентного спуска. Он должен:\n",
    "- инициализироваться величиной шага градиентного спуска и количеством итераций градиентного спуска\n",
    "- содержать метод fit для обучения модели, на вход которой будут подаваться обучающие данные\n",
    "- чтобы реализовать функцию sigmoid, можно воспользоваться библиотечной функцией np.exp\n",
    "- для расчета градиента необходимо использовать формулу, приведенную в курсе в блоке «Логистическая регрессия для решения задачи бинарной классификации».\n",
    "- должен содержаться метод predict, который будет возвращать лейбл 0/1.\n",
    "- лейбл должен рассчитываться из вероятности принадлежности классу 1.\n",
    "- вероятность должна рассчитываться как сигмоида от результата произведения матрицы фичей X и вектора весов модели w."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "class MyLogisticRegression:\n",
    "    def __init__(self, learning_rate=0.01, n_iterations=1000):\n",
    "        self.learning_rate = learning_rate\n",
    "        self.n_iterations = n_iterations\n",
    "        self.weights = None\n",
    "\n",
    "    def sigmoid(self, x):\n",
    "        return 1 / (1 + np.exp(-x))\n",
    "\n",
    "    def fit(self, X, y):\n",
    "        n_samples, n_features = X.shape\n",
    "        self.weights = np.zeros(n_features)\n",
    "\n",
    "        for _ in range(self.n_iterations):\n",
    "            linear_model = np.dot(X, self.weights)\n",
    "            predictions = self.sigmoid(linear_model)\n",
    "\n",
    "            gradient = np.dot(X.T, (predictions - y)) / n_samples\n",
    "            self.weights -= self.learning_rate * gradient\n",
    "\n",
    "    def predict(self, X):\n",
    "        linear_model = np.dot(X, self.weights)\n",
    "        predictions = self.sigmoid(linear_model)\n",
    "        return np.where(predictions >= 0.5, 1, 0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "5. Обучите модель с помощью реализованного класса с произвольными параметрами градиентного спуска."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [],
   "source": [
    "my_model = MyLogisticRegression(learning_rate=0.1, n_iterations=1000)\n",
    "my_model.fit(X_train, y_train)\n",
    "\n",
    "my_predictions = my_model.predict(X_test)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "6. Оцените качество модели, обученной с помощью класса MyLogisticRegression, выведя матрицу ошибок для прогнозов обученной моделью на тестовой части данных. Подберите такие параметры градиентного спуска, чтобы разделение получилось не хуже, чем у библиотечной модели."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<sklearn.metrics._plot.confusion_matrix.ConfusionMatrixDisplay at 0x2069feac4a0>"
      ]
     },
     "execution_count": 48,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf8AAAGwCAYAAACn/2wHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAzlElEQVR4nO3deXRV5b3/8c9OIAlDEghDQiAiiEwKQbFNQ0XgJxrSLmSoQ7n0EhBw1YpDUxS8FQRR06tXUQoFJwi2UsBWI6KXFlGmgvYyxIpCSkIgQRIEGUKCJOGc/fsD2fGYBHKyz8lJzn6/1nrWunt49v4eL+XL832evbdhmqYpAADgGCGBDgAAADQskj8AAA5D8gcAwGFI/gAAOAzJHwAAhyH5AwDgMCR/AAAcplmgA7DD7XbryJEjioyMlGEYgQ4HAOAl0zR15swZxcfHKyTEf+PRc+fOqaKiwvZ1wsLCFBER4YOIAqtJJ/8jR44oISEh0GEAAGwqLCxUly5d/HLtc+fOqVvX1ir+ymX7WnFxccrPz2/y/wBo0sk/MjJSknRo15WKas0MBoLTmJ79Ah0C4DfnVamtet/6+9wfKioqVPyVS4d2XqmoyPrnipIzbnUdeFAVFRUk/0C6WOqPah1i6/+hQGPWzGge6BAA//n2BfMNMXXbOtJQ68j638et4JlebtLJHwCAunKZbrlsfM3GZbp9F0yAkfwBAI7glim36p/97fRtbKiVAwDgMIz8AQCO4JZbdgr39no3LiR/AIAjuExTLrP+pXs7fRsbyv4AADgMI38AgCOw4K8KyR8A4AhumXKR/CVR9gcAwHEY+QMAHIGyfxWSPwDAEVjtX4WyPwAADsPIHwDgCO5vm53+wYLkDwBwBJfN1f52+jY2JH8AgCO4TNn8qp/vYgk05vwBAHAYRv4AAEdgzr8KyR8A4AhuGXLJsNU/WFD2BwDAYRj5AwAcwW1eaHb6BwuSPwDAEVw2y/52+jY2lP0BAHAYRv4AAEdg5F+FkT8AwBHcpmG7eWPz5s0aOXKk4uPjZRiGsrKyPI4bhlFje/bZZ2u95pw5c6qd37t3b6//W5D8AQDwg7KyMiUmJmrRokU1Hi8qKvJoS5culWEY+tnPfnbJ615zzTUe/bZu3ep1bJT9AQCO4Kuyf0lJicf+8PBwhYeHVzs/NTVVqamptV4vLi7OY/udd97RsGHD1L1790vG0axZs2p9vcXIHwDgCC6F2G6SlJCQoOjoaKtlZGTYju3o0aN67733NHny5Mueu3//fsXHx6t79+4aP368CgoKvL4fI38AgCOY9Zi3/35/SSosLFRUVJS1v6ZRv7eWL1+uyMhIjR079pLnJSUlKTMzU7169VJRUZHmzp2rwYMHa8+ePYqMjKzz/Uj+AAB4ISoqyiP5+8LSpUs1fvx4RUREXPK8704j9O/fX0lJSeratatWr15dp6rBRSR/AIAjNNZH/bZs2aKcnBytWrXK675t2rRRz549lZub61U/5vwBAI7gMkNsN3947bXXNHDgQCUmJnrdt7S0VHl5eerUqZNX/Uj+AAD4QWlpqbKzs5WdnS1Jys/PV3Z2tscCvZKSEr355puaMmVKjde4+eabtXDhQmt7+vTp2rRpkw4ePKht27ZpzJgxCg0N1bhx47yKjbI/AMAR3DLktjHmdcu7L/vs2LFDw4YNs7bT09MlSWlpacrMzJQkrVy5UqZp1pq88/LydPz4cWv78OHDGjdunL7++mt16NBBN954oz7++GN16NDBq9gM0zSb7HeKSkpKFB0drZP/7q6oSIoYCE4p8QMCHQLgN+fNSm3UOzp9+rTPF9FddDFXrPnXVWoVGVrv65Sdcem2/nl+jbWhkDEBAHAYyv4AAEewu2jP1XQL5dWQ/AEAjnBhzr/+j+vZ6dvYUPYHAMBhGPkDABzB/Z3389evP2V/AACaFOb8q5D8AQCO4FZIgz7n35gx5w8AgMMw8gcAOILLNOSy8UlfO30bG5I/AMARXDYX/Lko+wMAgKaKkT8AwBHcZojcNlb7u1ntDwBA00LZvwplfwAAHIaRPwDAEdyyt2Lf7btQAo7kDwBwBPsv+QmeYnnw/BIAAFAnjPwBAI5g/93+wTNeJvkDABzBLUNu2Znz5w1/AAA0KYz8qwTPLwEAAHXCyB8A4Aj2X/ITPONlkj8AwBHcpiG3nef8g+irfsHzzxgAAFAnjPwBAI7gtln2D6aX/JD8AQCOYP+rfsGT/IPnlwAAgDph5A8AcASXDLlsvKjHTt/GhuQPAHAEyv5VgueXAACAOmHkDwBwBJfsle5dvgsl4Ej+AABHoOxfheQPAHAEPuxTJXh+CQAAqBNG/gAARzBlyG1jzt/kUT8AAJoWyv5VgueXAACAOiH5AwAc4eInfe00b2zevFkjR45UfHy8DMNQVlaWx/GJEyfKMAyPNmLEiMted9GiRbryyisVERGhpKQk/fOf//QqLonkDwBwCNe3X/Wz07xRVlamxMRELVq0qNZzRowYoaKiIqv9+c9/vuQ1V61apfT0dD3++OPatWuXEhMTlZKSoq+++sqr2JjzBwDAD1JTU5WamnrJc8LDwxUXF1fnaz7//POaOnWqJk2aJElasmSJ3nvvPS1dulQzZ86s83UY+QMAHMFXZf+SkhKPVl5eXu+YNm7cqI4dO6pXr16699579fXXX9d6bkVFhXbu3Knhw4db+0JCQjR8+HBt377dq/uS/AEAjuBWiO0mSQkJCYqOjrZaRkZGveIZMWKEXn/9dW3YsEH//d//rU2bNik1NVUuV80vEj5+/LhcLpdiY2M99sfGxqq4uNire1P2BwDAC4WFhYqKirK2w8PD63Wdn//859b/3a9fP/Xv319XXXWVNm7cqJtvvtl2nJfCyB8A4Agu07DdJCkqKsqj1Tf5f1/37t3Vvn175ebm1ni8ffv2Cg0N1dGjRz32Hz161Kt1AxLJHwDgEA39qJ+3Dh8+rK+//lqdOnWq8XhYWJgGDhyoDRs2VP0mt1sbNmxQcnKyV/ei7A8AcATT5lf9TC/7lpaWeozi8/PzlZ2drZiYGMXExGju3Ln62c9+pri4OOXl5emRRx5Rjx49lJKSYvW5+eabNWbMGE2bNk2SlJ6errS0NN1www364Q9/qBdeeEFlZWXW6v+6IvkDAOAHO3bs0LBhw6zt9PR0SVJaWpoWL16sf/3rX1q+fLlOnTql+Ph43XrrrZo3b57HNEJeXp6OHz9ubd911106duyYZs+ereLiYg0YMEDr1q2rtgjwckj+AABHcMmQy8bHebztO3ToUJmmWevxv/3tb5e9xsGDB6vtmzZtmlUJqC+SPwDAEdymbM3bu2vP400OC/4AAHAYRv6o5rOPW+nNP3TU/s9a6sTR5nr8tXwNSj1tHT95rJleeypeOzdFqux0qK79Uanue/KwOnevCGDUgH0jJx7X7fd+pZgO53Xgixb6w2OdlZPdMtBhwUfcNhf82enb2ATPL4HPnDsbou7XfKNpTx+udsw0pbl3d1PRoTDNWXZAi/6eo9guFZp5Vw+dO8sfJzRdQ247qXseP6I3no/TfSk9deCLCD214oCi21UGOjT4iFuG7RYsGsXf1r74PCF85wf/74wmzijWj78z2r/oywPh2ruzle7/3WH1GvCNEnqU6/7fHVb5OUMfvd2m4YMFfGTsPce1bkWM/r4qRgX7I7RgRheVf2MoZdyJQIcG+FzAk7+vPk+IhlFZceFfvmHhbmtfSIjUPMzU5//XOlBhAbY0a+7W1f3PateWSGufaRravSVSfQeeDWBk8CVfveEvGAQ8+X/384R9+/bVkiVL1LJlSy1dujTQoaEGCT3OqWPnCi3N6KQzp0JVWWFo1cKOOl4UphNHWUKCpikqxqXQZtKpY55/hk8eb6a2Hc4HKCr42sU5fzstWAT0l3j7ecLy8vJqn1JEw2rWXJr9Wr6+zIvQ7X376bar+uvTba31g/9XIiN4/ncBAEEtoEO1S32ecN++fdXOz8jI0Ny5cxsqPNTi6v7faPEHOSorCVFlpaE27Vx64KdXq2d/yqNomkpOhMp1XmrzvVF+2/bndfIYFa1g4Za99/Oz4C9AHn30UZ0+fdpqhYWFgQ7J0VpFudWmnUtfHgjT/k9bKjmFSgyapvOVIdr/r5a67sYz1j7DMDXgxlJ9sZNH/YKFaXOlvxlEyT+g/6T19vOE4eHhPvt0Imr3TVmIjuRX/XcuLgxT3p4WimxzXh27VGrzu9GKbudSx84Vyt8boSWzuyh5xGkNHHrmElcFGre3Xm6v6S8U6t+ftlTO7pYaM/WYIlq69feVMYEODT5i98t8/v6qX0MKaPL/7ucJR48eLanq84R231uM+vv3py31yO09rO2X5nSWJN1y5wlNf6FAJ44210tzOuvU8WaK6Xhew+84of946GhtlwOahE1r2iq6nUsTHi5W2w7ndeDzFvrt+G46dbx5oEMDfC7gk1m++jwhfCdxUKn+diS71uOjpxzX6CnHaz0ONFVrlrXXmmXtAx0G/IQ3/FUJePL31ecJAQC4FMr+VQKe/CXffJ4QAADUTaNI/gAA+Jvd9/MH06N+JH8AgCNQ9q8SPKsXAABAnTDyBwA4AiP/KiR/AIAjkPyrUPYHAMBhGPkDAByBkX8Vkj8AwBFM2Xtcz/RdKAFH8gcAOAIj/yrM+QMA4DCM/AEAjsDIvwrJHwDgCCT/KpT9AQBwGEb+AABHYORfheQPAHAE0zRk2kjgdvo2NpT9AQBwGEb+AABHcMuw9ZIfO30bG5I/AMARmPOvQtkfAACHYeQPAHAEFvxVIfkDAByBsn8Vkj8AwBEY+Vdhzh8AAIch+QMAHMH8tuxf3+btyH/z5s0aOXKk4uPjZRiGsrKyrGOVlZWaMWOG+vXrp1atWik+Pl4TJkzQkSNHLnnNOXPmyDAMj9a7d2+v/1uQ/AEAjmBKMk0bzcv7lZWVKTExUYsWLap27OzZs9q1a5dmzZqlXbt26a233lJOTo5uu+22y173mmuuUVFRkdW2bt3qZWTM+QMA4BepqalKTU2t8Vh0dLTWr1/vsW/hwoX64Q9/qIKCAl1xxRW1XrdZs2aKi4uzFRsjfwCAI1x8w5+dJkklJSUerby83CfxnT59WoZhqE2bNpc8b//+/YqPj1f37t01fvx4FRQUeH0vkj8AwBEurva30yQpISFB0dHRVsvIyLAd27lz5zRjxgyNGzdOUVFRtZ6XlJSkzMxMrVu3TosXL1Z+fr4GDx6sM2fOeHU/yv4AAHihsLDQI0GHh4fbul5lZaXuvPNOmaapxYsXX/Lc704j9O/fX0lJSeratatWr16tyZMn1/meJH8AgCO4TUOGD17yExUVdcnRuTcuJv5Dhw7pww8/9Pq6bdq0Uc+ePZWbm+tVP8r+AABHsLXS/9vmSxcT//79+/XBBx+oXbt2Xl+jtLRUeXl56tSpk1f9SP4AAPhBaWmpsrOzlZ2dLUnKz89Xdna2CgoKVFlZqdtvv107duzQG2+8IZfLpeLiYhUXF6uiosK6xs0336yFCxda29OnT9emTZt08OBBbdu2TWPGjFFoaKjGjRvnVWyU/QEAjtDQr/fdsWOHhg0bZm2np6dLktLS0jRnzhytWbNGkjRgwACPfh999JGGDh0qScrLy9Px48etY4cPH9a4ceP09ddfq0OHDrrxxhv18ccfq0OHDl7FRvIHADhCQyf/oUOHyrzEXMGljl108OBBj+2VK1d6FUNtSP4AAEfw1YK/YMCcPwAADsPIHwDgCHZX7Pt6tX8gkfwBAI5wIfnbmfP3YTABRtkfAACHYeQPAHCEhl7t35iR/AEAjmB+2+z0DxaU/QEAcBhG/gAAR6DsX4XkDwBwBur+FpI/AMAZbI78FUQjf+b8AQBwGEb+AABH4A1/VUj+AABHYMFfFcr+AAA4DCN/AIAzmIa9RXtBNPIn+QMAHIE5/yqU/QEAcBhG/gAAZ+AlPxaSPwDAEVjtX6VOyX/NmjV1vuBtt91W72AAAID/1Sn5jx49uk4XMwxDLpfLTjwAAPhPEJXu7ahT8ne73f6OAwAAv6LsX8XWav9z5875Kg4AAPzL9EELEl4nf5fLpXnz5qlz585q3bq1Dhw4IEmaNWuWXnvtNZ8HCAAAfMvr5P/UU08pMzNTzzzzjMLCwqz91157rV599VWfBgcAgO8YPmjBwevk//rrr+vll1/W+PHjFRoaau1PTEzUvn37fBocAAA+Q9nf4nXy//LLL9WjR49q+91utyorK30SFAAA8B+vk3/fvn21ZcuWavv/8pe/6LrrrvNJUAAA+Bwjf4vXb/ibPXu20tLS9OWXX8rtduutt95STk6OXn/9da1du9YfMQIAYB9f9bN4PfIfNWqU3n33XX3wwQdq1aqVZs+erb179+rdd9/VLbfc4o8YAQCAD9Xr3f6DBw/W+vXrfR0LAAB+wyd9q9T7wz47duzQ3r17JV1YBzBw4ECfBQUAgM/xVT+L18n/8OHDGjdunP7xj3+oTZs2kqRTp05p0KBBWrlypbp06eLrGAEAgA95Pec/ZcoUVVZWau/evTpx4oROnDihvXv3yu12a8qUKf6IEQAA+y4u+LPTgoTXI/9NmzZp27Zt6tWrl7WvV69e+v3vf6/Bgwf7NDgAAHzFMC80O/2DhdfJPyEhocaX+bhcLsXHx/skKAAAfI45f4vXZf9nn31W999/v3bs2GHt27Fjhx588EH9z//8j0+DAwAAvlen5N+2bVvFxMQoJiZGkyZNUnZ2tpKSkhQeHq7w8HAlJSVp165duvvuu/0dLwAA9dPAc/6bN2/WyJEjFR8fL8MwlJWV5RmOaWr27Nnq1KmTWrRooeHDh2v//v2Xve6iRYt05ZVXKiIiQklJSfrnP//pVVxSHcv+L7zwgtcXBgCgUWngsn9ZWZkSExN19913a+zYsdWOP/PMM1qwYIGWL1+ubt26adasWUpJSdEXX3yhiIiIGq+5atUqpaena8mSJUpKStILL7yglJQU5eTkqGPHjnWOrU7JPy0trc4XBAAgmJWUlHhsX6yCf19qaqpSU1NrvIZpmnrhhRf02GOPadSoUZIufDU3NjZWWVlZ+vnPf15jv+eff15Tp07VpEmTJElLlizRe++9p6VLl2rmzJl1/g1ez/l/17lz51RSUuLRAABolHz0YZ+EhARFR0dbLSMjw+tQ8vPzVVxcrOHDh1v7oqOjlZSUpO3bt9fYp6KiQjt37vToExISouHDh9fapzZer/YvKyvTjBkztHr1an399dfVjrtcLm8vCQCA//mo7F9YWKioqChrd02j/sspLi6WJMXGxnrsj42NtY593/Hjx+VyuWrss2/fPq/u7/XI/5FHHtGHH36oxYsXKzw8XK+++qrmzp2r+Ph4vf76695eDgCAJiUqKsqj1Sf5B5rXyf/dd9/VH/7wB/3sZz9Ts2bNNHjwYD322GN6+umn9cYbb/gjRgAA7GtEb/iLi4uTJB09etRj/9GjR61j39e+fXuFhoZ61ac2Xif/EydOqHv37pIu/OvnxIkTkqQbb7xRmzdv9vZyAAA0iItv+LPTfKVbt26Ki4vThg0brH0lJSX65JNPlJycXGOfsLAwDRw40KOP2+3Whg0bau1TG6+Tf/fu3ZWfny9J6t27t1avXi3pQkXg4od+AABwutLSUmVnZys7O1vShUV+2dnZKigokGEYeuihh/Tkk09qzZo1+uyzzzRhwgTFx8dr9OjR1jVuvvlmLVy40NpOT0/XK6+8ouXLl2vv3r269957VVZWZq3+ryuvF/xNmjRJn376qYYMGaKZM2dq5MiRWrhwoSorK/X88897ezkAABpGAz/nv2PHDg0bNszaTk9Pl3Th8fnMzEw98sgjKisr0z333KNTp07pxhtv1Lp16zye8c/Ly9Px48et7bvuukvHjh3T7NmzVVxcrAEDBmjdunXVFgFejmGapq1CxqFDh7Rz50716NFD/fv3t3Mpr5WUlCg6Olon/91dUZG2nloEGq2U+AGBDgHwm/NmpTbqHZ0+fdpjBb0vXcwVV/z3kwppUfPLc+rC/c05Fcx4zK+xNhSvR/7f17VrV3Xt2tUXsQAA4DeGbH7Vz2eRBF6dkv+CBQvqfMEHHnig3sEAAAD/q1Pynz9/fp0uZhhGQJL/mJ791Mxo3uD3BRrC/tevD3QIgN+4vzkn3fNOw9zM7uN6PnzUL9DqlPwvru4HAKDJauAFf40Zq+QAAHAY2wv+AABoEhj5W0j+AABHsPuWPl++4S/QKPsDAOAwjPwBAM5A2d9Sr5H/li1b9Itf/ELJycn68ssvJUl//OMftXXrVp8GBwCAz5g+aEHC6+T/17/+VSkpKWrRooV2796t8vJySdLp06f19NNP+zxAAADgW14n/yeffFJLlizRK6+8oubNq16s8+Mf/1i7du3yaXAAAPhKY/qkb6B5Peefk5Ojm266qdr+6OhonTp1yhcxAQDge7zhz+L1yD8uLk65ubnV9m/dulXdu3f3SVAAAPgcc/4Wr5P/1KlT9eCDD+qTTz6RYRg6cuSI3njjDU2fPl333nuvP2IEAAA+5HXZf+bMmXK73br55pt19uxZ3XTTTQoPD9f06dN1//33+yNGAABs4yU/VbxO/oZh6Le//a0efvhh5ebmqrS0VH379lXr1q39ER8AAL7Bc/6Wer/kJywsTH379vVlLAAAoAF4nfyHDRsmw6h9xeOHH35oKyAAAPzC7uN6Th75DxgwwGO7srJS2dnZ2rNnj9LS0nwVFwAAvkXZ3+J18p8/f36N++fMmaPS0lLbAQEAAP/y2Vf9fvGLX2jp0qW+uhwAAL7Fc/4Wn33Vb/v27YqIiPDV5QAA8Cke9avidfIfO3asx7ZpmioqKtKOHTs0a9YsnwUGAAD8w+vkHx0d7bEdEhKiXr166YknntCtt97qs8AAAIB/eJX8XS6XJk2apH79+qlt27b+igkAAN9jtb/FqwV/oaGhuvXWW/l6HwCgyeGTvlW8Xu1/7bXX6sCBA/6IBQAANACvk/+TTz6p6dOna+3atSoqKlJJSYlHAwCg0eIxP0lezPk/8cQT+s1vfqOf/OQnkqTbbrvN4zW/pmnKMAy5XC7fRwkAgF3M+VvqnPznzp2rX/7yl/roo4/8GQ8AAPCzOid/07zwT54hQ4b4LRgAAPyFl/xU8epRv0t9zQ8AgEaNsr/Fq+Tfs2fPy/4D4MSJE7YCAgAA/uVV8p87d261N/wBANAUUPav4lXy//nPf66OHTv6KxYAAPyHsr+lzs/5M98PAEBwqHPyv7jaHwCAJsnOC37qUTW48sorZRhGtXbffffVeH5mZma1cyMiIurxQy+vzmV/t9vtlwAAAGgIDT3n/3//938eL77bs2ePbrnlFt1xxx219omKilJOTk7VPf1Udff6k74AADRJDTzn36FDB4/t3/3ud7rqqqsu+b4cwzAUFxdXn+i84vW7/QEAcLLvf9OmvLz8sn0qKir0pz/9SXffffclR/OlpaXq2rWrEhISNGrUKH3++ee+DN1C8gcAOIOP5vwTEhIUHR1ttYyMjMveOisrS6dOndLEiRNrPadXr15aunSp3nnnHf3pT3+S2+3WoEGDdPjw4Xr+4NpR9gcAOIKv5vwLCwsVFRVl7Q8PD79s39dee02pqamKj4+v9Zzk5GQlJydb24MGDVKfPn300ksvad68efUPvAYkfwAAvBAVFeWR/C/n0KFD+uCDD/TWW295dZ/mzZvruuuuU25urrchXhZlfwCAMzTwo34XLVu2TB07dtRPf/pTr/q5XC599tln6tSpU/1ufAmM/AEAjhCI1/u63W4tW7ZMaWlpatbMM+VOmDBBnTt3ttYMPPHEE/rRj36kHj166NSpU3r22Wd16NAhTZkypf5B14LkDwCAn3zwwQcqKCjQ3XffXe1YQUGBQkKqCvAnT57U1KlTVVxcrLZt22rgwIHatm2b+vbt6/O4SP4AAGcIwLv9b7311lrfkLtx40aP7fnz52v+/Pn1CMx7JH8AgDPwYR8LC/4AAHAYRv4AAEcwvm12+gcLkj8AwBko+1tI/gAARwjEo36NFXP+AAA4DCN/AIAzUPa3kPwBAM4RRAncDsr+AAA4DCN/AIAjsOCvCskfAOAMzPlbKPsDAOAwjPwBAI5A2b8KyR8A4AyU/S2U/QEAcBhG/gAAR6DsX4XkDwBwBsr+FpI/AMAZSP4W5vwBAHAYRv4AAEdgzr8KyR8A4AyU/S2U/QEAcBhG/gAARzBMU4ZZ/+G7nb6NDckfAOAMlP0tlP0BAHAYRv4AAEdgtX8Vkj8AwBko+1so+wMA4DCM/AEAjkDZvwrJHwDgDJT9LSR/AIAjMPKvwpw/AAAOw8gfAOAMlP0tJH8AgGMEU+neDsr+AAA4DCN/AIAzmOaFZqd/kCD5AwAcgdX+VSj7AwDgMCR/AIAzmD5oXpgzZ44Mw/BovXv3vmSfN998U71791ZERIT69eun999/37ub1hHJHwDgCIbbfvPWNddco6KiIqtt3bq11nO3bdumcePGafLkydq9e7dGjx6t0aNHa8+ePTZ+dc2Y8wcAwAslJSUe2+Hh4QoPD6/x3GbNmikuLq5O133xxRc1YsQIPfzww5KkefPmaf369Vq4cKGWLFliL+jvYeSPOhs58biWf/KF3j3wL724dr96DTgb6JCAeonYd0adns9Vtwc+09UTdqnVzlMex1v930nFP7Nf3e/9VFdP2KWwQ/xZDwo+KvsnJCQoOjraahkZGbXecv/+/YqPj1f37t01fvx4FRQU1Hru9u3bNXz4cI99KSkp2r59e71+7qWQ/FEnQ247qXseP6I3no/TfSk9deCLCD214oCi21UGOjTAayHlblVc0VJfTUio+XiFW+d6ttbxuzo3cGTwp4ur/e00SSosLNTp06et9uijj9Z4v6SkJGVmZmrdunVavHix8vPzNXjwYJ05c6bG84uLixUbG+uxLzY2VsXFxT797yAFOPlv3rxZI0eOVHx8vAzDUFZWViDDwSWMvee41q2I0d9Xxahgf4QWzOii8m8MpYw7EejQAK+dTYzW17fHq+yGNjUeP/PjdjoxupPOXhPZsIHBvy4+52+nSYqKivJotZX8U1NTdccdd6h///5KSUnR+++/r1OnTmn16tUN+atrFNDkX1ZWpsTERC1atCiQYeAymjV36+r+Z7VrS9VfhKZpaPeWSPUdSDkUAOqiTZs26tmzp3Jzc2s8HhcXp6NHj3rsO3r0aJ3XDHgjoMk/NTVVTz75pMaMGVOn88vLy1VSUuLR4H9RMS6FNpNOHfNcH3ryeDO17XA+QFEBgHd8Vfavr9LSUuXl5alTp041Hk9OTtaGDRs89q1fv17Jycn2blyDJjXnn5GR4bHIIiGh5vk6AACqaeDn/KdPn65Nmzbp4MGD2rZtm8aMGaPQ0FCNGzdOkjRhwgSP9QIPPvig1q1bp+eee0779u3TnDlztGPHDk2bNs3Or65Rk0r+jz76qMcii8LCwkCH5AglJ0LlOi+1+d4ov2378zp5jKdFAaAmhw8f1rhx49SrVy/deeedateunT7++GN16NBBklRQUKCioiLr/EGDBmnFihV6+eWXlZiYqL/85S/KysrStdde6/PYmtTf3Jd6lhL+c74yRPv/1VLX3XhG29dFS5IMw9SAG0u1JrNdgKMDgLpp6Hf7r1y58pLHN27cWG3fHXfcoTvuuMO7G9VDk0r+CJy3Xm6v6S8U6t+ftlTO7pYaM/WYIlq69feVMYEODfCacc6l5kfLre3mx8oVduis3K2a6Xz7MIWUnlezryvU7NSFR1nDis5JklzRzeVq0zwgMcMH+KqfheSPOtm0pq2i27k04eFite1wXgc+b6Hfju+mU8f5ixBNT0T+WXXJ2G9td1jxpSSp5MYYHb3nSrXafVpxrxyyjnf6w0FJ0tej43RibHyDxgr4Q0CTf2lpqccjD/n5+crOzlZMTIyuuOKKAEaGmqxZ1l5rlrUPdBiAbd/0idT+16+v9fiZwe10ZjBTWsGGT/pWCWjy37Fjh4YNG2Ztp6enS5LS0tKUmZkZoKgAAEGpHiv2q/UPEgFN/kOHDpUZRHMoAAA0Bcz5AwAcgbJ/FZI/AMAZ3OaFZqd/kCD5AwCcgTl/S5N6wx8AALCPkT8AwBEM2Zzz91kkgUfyBwA4A2/4s1D2BwDAYRj5AwAcgUf9qpD8AQDOwGp/C2V/AAAchpE/AMARDNOUYWPRnp2+jQ3JHwDgDO5vm53+QYKyPwAADsPIHwDgCJT9q5D8AQDOwGp/C8kfAOAMvOHPwpw/AAAOw8gfAOAIvOGvCskfAOAMlP0tlP0BAHAYRv4AAEcw3Beanf7BguQPAHAGyv4Wyv4AADgMI38AgDPwkh8LyR8A4Ai83rcKZX8AAByGkT8AwBlY8Gch+QMAnMGUZOdxveDJ/SR/AIAzMOdfhTl/AAAchpE/AMAZTNmc8/dZJAFH8gcAOAML/iyU/QEA8IOMjAz94Ac/UGRkpDp27KjRo0crJyfnkn0yMzNlGIZHi4iI8HlsJH8AgDO4fdC8sGnTJt133336+OOPtX79elVWVurWW29VWVnZJftFRUWpqKjIaocOHfLuxnVA2R8A4AgNvdp/3bp1HtuZmZnq2LGjdu7cqZtuuqn2+xiG4uLi6hVjXTHyBwDACyUlJR6tvLy8Tv1Onz4tSYqJibnkeaWlperatasSEhI0atQoff7557Zj/j6SPwDAGS4u+LPTJCUkJCg6OtpqGRkZl7212+3WQw89pB//+Me69tpraz2vV69eWrp0qd555x396U9/ktvt1qBBg3T48GGf/WeQKPsDAJzCR6v9CwsLFRUVZe0ODw+/bNf77rtPe/bs0datWy95XnJyspKTk63tQYMGqU+fPnrppZc0b968egZeHckfAAAvREVFeST/y5k2bZrWrl2rzZs3q0uXLl7dq3nz5rruuuuUm5vrbZiXRNkfAOAMPir71/12pqZNm6a3335bH374obp16+Z1yC6XS5999pk6derkdd9LYeQPAHAGtyTDZn8v3HfffVqxYoXeeecdRUZGqri4WJIUHR2tFi1aSJImTJigzp07W+sGnnjiCf3oRz9Sjx49dOrUKT377LM6dOiQpkyZYiPw6kj+AABHaOhH/RYvXixJGjp0qMf+ZcuWaeLEiZKkgoIChYRUFeFPnjypqVOnqri4WG3bttXAgQO1bds29e3bt95x14TkDwCAH5h1+MfCxo0bPbbnz5+v+fPn+ymiKiR/AIAz8G5/C8kfAOAMblMybCRwd/Akf1b7AwDgMIz8AQDOQNnfQvIHADiEzeSv4En+lP0BAHAYRv4AAGeg7G8h+QMAnMFtylbpntX+AACgqWLkDwBwBtN9odnpHyRI/gAAZ2DO30LyBwA4A3P+Fub8AQBwGEb+AABnoOxvIfkDAJzBlM3k77NIAo6yPwAADsPIHwDgDJT9LSR/AIAzuN2SbDyr7w6e5/wp+wMA4DCM/AEAzkDZ30LyBwA4A8nfQtkfAACHYeQPAHAGXu9rIfkDABzBNN0ybXyZz07fxobkDwBwBtO0N3pnzh8AADRVjPwBAM5g2pzzD6KRP8kfAOAMbrdk2Ji3D6I5f8r+AAA4DCN/AIAzUPa3kPwBAI5gut0ybZT9g+lRP8r+AAA4DCN/AIAzUPa3kPwBAM7gNiWD5C9R9gcAwHEY+QMAnME0Jdl5zj94Rv4kfwCAI5huU6aNsr9J8gcAoIkx3bI38udRPwAAUAeLFi3SlVdeqYiICCUlJemf//znJc9/88031bt3b0VERKhfv356//33fR4TyR8A4Aim27TdvLVq1Sqlp6fr8ccf165du5SYmKiUlBR99dVXNZ6/bds2jRs3TpMnT9bu3bs1evRojR49Wnv27LH78z2Q/AEAzmC67TcvPf/885o6daomTZqkvn37asmSJWrZsqWWLl1a4/kvvviiRowYoYcfflh9+vTRvHnzdP3112vhwoV2f72HJj3nf3HxxXlV2npvA9CYub85F+gQAL9xf1MuqWEW09nNFedVKUkqKSnx2B8eHq7w8PBq51dUVGjnzp169NFHrX0hISEaPny4tm/fXuM9tm/frvT0dI99KSkpysrKqn/gNWjSyf/MmTOSpK3y/XwI0Gjc806gIwD87syZM4qOjvbLtcPCwhQXF6etxfZzRevWrZWQkOCx7/HHH9ecOXOqnXv8+HG5XC7FxsZ67I+NjdW+fftqvH5xcXGN5xcXF9sL/HuadPKPj49XYWGhIiMjZRhGoMNxhJKSEiUkJKiwsFBRUVGBDgfwKf58NzzTNHXmzBnFx8f77R4RERHKz89XRUWF7WuZplkt39Q06m/smnTyDwkJUZcuXQIdhiNFRUXxlyOCFn++G5a/RvzfFRERoYiICL/f57vat2+v0NBQHT161GP/0aNHFRcXV2OfuLg4r86vLxb8AQDgB2FhYRo4cKA2bNhg7XO73dqwYYOSk5Nr7JOcnOxxviStX7++1vPrq0mP/AEAaMzS09OVlpamG264QT/84Q/1wgsvqKysTJMmTZIkTZgwQZ07d1ZGRoYk6cEHH9SQIUP03HPP6ac//alWrlypHTt26OWXX/ZpXCR/eCU8PFyPP/54k5zjAi6HP9/wtbvuukvHjh3T7NmzVVxcrAEDBmjdunXWor6CggKFhFQV4QcNGqQVK1boscce03/913/p6quvVlZWlq699lqfxmWYwfSyYgAAcFnM+QMA4DAkfwAAHIbkDwCAw5D8AQBwGJI/6szbz1ICTcXmzZs1cuRIxcfHyzAMn79HHWhsSP6oE28/Swk0JWVlZUpMTNSiRYsCHQrQIHjUD3WSlJSkH/zgB9ZnJd1utxISEnT//fdr5syZAY4O8B3DMPT2229r9OjRgQ4F8BtG/risi5+lHD58uLXvcp+lBAA0XiR/XNalPkvp689MAgD8j+QPAIDDkPxxWfX5LCUAoPEi+eOy6vNZSgBA48VX/VAnl/ssJdCUlZaWKjc319rOz89Xdna2YmJidMUVVwQwMsA/eNQPdbZw4UI9++yz1mcpFyxYoKSkpECHBdi2ceNGDRs2rNr+tLQ0ZWZmNnxAgJ+R/AEAcBjm/AEAcBiSPwAADkPyBwDAYUj+AAA4DMkfAACHIfkDAOAwJH8AAByG5A8AgMOQ/AGbJk6cqNGjR1vbQ4cO1UMPPdTgcWzcuFGGYejUqVO1nmMYhrKysup8zTlz5mjAgAG24jp48KAMw1B2drat6wDwHZI/gtLEiRNlGIYMw1BYWJh69OihJ554QufPn/f7vd966y3NmzevTufWJWEDgK/xYR8ErREjRmjZsmUqLy/X+++/r/vuu0/NmzfXo48+Wu3ciooKhYWF+eS+MTExPrkOAPgLI38ErfDwcMXFxalr16669957NXz4cK1Zs0ZSVan+qaeeUnx8vHr16iVJKiws1J133qk2bdooJiZGo0aN0sGDB61rulwupaenq02bNmrXrp0eeeQRff/zGN8v+5eXl2vGjBlKSEhQeHi4evTooddee00HDx60PibTtm1bGYahiRMnSrrwyeSMjAx169ZNLVq0UGJiov7yl7943Of9999Xz5491aJFCw0bNswjzrqaMWOGevbsqZYtW6p79+6aNWuWKisrq5330ksvKSEhQS1bttSdd96p06dPexx/9dVX1adPH0VERKh37976wx/+4HUsABoOyR+O0aJFC1VUVFjbGzZsUE5OjtavX6+1a9eqsrJSKSkpioyM1JYtW/SPf/xDrVu31ogRI6x+zz33nDIzM7V06VJt3bpVJ06c0Ntvv33J+06YMEF//vOftWDBAu3du1cvvfSSWrdurYSEBP31r3+VJOXk5KioqEgvvviiJCkjI0Ovv/66lixZos8//1y//vWv9Ytf/EKbNm2SdOEfKWPHjtXIkSOVnZ2tKVOmaObMmV7/N4mMjFRmZqa++OILvfjii3rllVc0f/58j3Nyc3O1evVqvfvuu1q3bp12796tX/3qV9bxN954Q7Nnz9ZTTz2lvXv36umnn9asWbO0fPlyr+MB0EBMIAilpaWZo0aNMk3TNN1ut7l+/XozPDzcnD59unU8NjbWLC8vt/r88Y9/NHv16mW63W5rX3l5udmiRQvzb3/7m2maptmpUyfzmWeesY5XVlaaXbp0se5lmqY5ZMgQ88EHHzRN0zRzcnJMSeb69etrjPOjjz4yJZknT5609p07d85s2bKluW3bNo9zJ0+ebI4bN840TdN89NFHzb59+3ocnzFjRrVrfZ8k8+233671+LPPPmsOHDjQ2n788cfN0NBQ8/Dhw9a+//3f/zVDQkLMoqIi0zRN86qrrjJXrFjhcZ158+aZycnJpmmaZn5+vinJ3L17d633BdCwmPNH0Fq7dq1at26tyspKud1u/cd//IfmzJljHe/Xr5/HPP+nn36q3NxcRUZGelzn3LlzysvL0+nTp1VUVKSkpCTrWLNmzXTDDTdUK/1flJ2drdDQUA0ZMqTOcefm5urs2bO65ZZbPPZXVFTouuuukyTt3bvXIw5JSk5OrvM9Llq1apUWLFigvLw8lZaW6vz584qKivI454orrlDnzp097uN2u5WTk6PIyEjl5eVp8uTJmjp1qnXO+fPnFR0d7XU8ABoGyR9Ba9iwYVq8eLHCwsIUHx+vZs08/7i3atXKY7u0tFQDBw7UG2+8Ue1aHTp0qFcMLVq08LpPaWmpJOm9997zSLrShXUMvrJ9+3aNHz9ec+fOVUpKiqKjo7Vy5Uo999xzXsf6yiuvVPvHSGhoqM9iBeBbJH8ErVatWqlHjx51Pv/666/XqlWr1LFjx2qj34s6deqkTz75RDfddJOkCyPcnTt36vrrr6/x/H79+sntdmvTpk0aPnx4teMXKw8ul8va17dvX4WHh6ugoKDWikGfPn2sxYsXffzxx5f/kd+xbds2de3aVb/97W+tfYcOHap2XkFBgY4cOaL4+HjrPiEhIerVq5diY2MVHx+vAwcOaPz48V7dH0DgsOAP+Nb48ePVvn17jRo1Slu2bFF+fr42btyoBx54QIcPH5YkPfjgg/rd736nrKws7du3T7/61a8u+Yz+lVdeqbS0NN19993Kysqyrrl69WpJUteuXWUYhtauXatjx46ptLRUkZGRmj59un79619r+fLlysvL065du/T73//eWkT3y1/+Uvv379fDDz+snJwcrVixQpmZmV793quvvloFBQVauXKl8vLytGDBghoXL0ZERCgtLU2ffvqptmzZogceeEB33nmn4uLiJElz585VRkaGFixYoH//+9/67LPPtGzZMj3//PNexQOg4ZD8gW+1bNlSmzdv1hVXXKGxY8eqT58+mjx5ss6dO2dVAn7zm9/oP//zP5WWlqbk5GRFRkZqzJgxl7zu4sWLdfvtt+tXv/qVevfuralTp6qsrEyS1LlzZ82dO1czZ85UbGyspk2bJkmaN2+eZs2apYyMDPXp00cjRozQe++9p27dukm6MA//17/+VVlZWUpMTNSSJUv09NNPe/V7b7vtNv3617/WtGnTNGDAAG3btk2zZs2qdl6PHj00duxY/eQnP9Gtt96q/v37ezzKN2XKFL366qtatmyZ+vXrpyFDhigzM9OKFUDjY5i1rVQCAABBiZE/AAAOQ/IHAMBhSP4AADgMyR8AAIch+QMA4DAkfwAAHIbkDwCAw5D8AQBwGJI/AAAOQ/IHAMBhSP4AADjM/weaS5sLoNazdAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ConfusionMatrixDisplay(confusion_matrix(y_test, my_predictions)).plot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Вывод:\n",
    "\n",
    "Качество модели, обученной с помощью класса MyLogisticRegression, получилось не хуже:\n",
    "- полнота = 1\n",
    "- точность = 1\n",
    "- F1-мера = 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Евклидово расстояние между векторами: 577.3214009544423\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "vector1 = np.arange(1, 101)\n",
    "vector2 = np.arange(100, 0, -1)\n",
    "\n",
    "euclidean_distance = np.linalg.norm(vector1 - vector2)\n",
    "\n",
    "print(\"Евклидово расстояние между векторами:\", euclidean_distance)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Iteration 1: x = 8.0, loss = -4.0\n",
      "Iteration 2: x = 9.6, loss = -0.8000000000000007\n",
      "Iteration 3: x = 9.92, loss = -0.16000000000000014\n",
      "Iteration 4: x = 9.984, loss = -0.03200000000000003\n"
     ]
    }
   ],
   "source": [
    "def gradient_descent(x, learning_rate, iterations):\n",
    "    losses = []\n",
    "    for i in range(iterations):\n",
    "        gradient = 2*(x-10)\n",
    "        x = x - learning_rate * gradient\n",
    "        loss = (x-10)*2\n",
    "        losses.append(loss)\n",
    "        print(f\"Iteration {i+1}: x = {x}, loss = {loss}\")\n",
    "    return x, losses\n",
    "\n",
    "x = 0\n",
    "learning_rate = 0.4\n",
    "iterations = 4\n",
    "\n",
    "final_x, losses = gradient_descent(x, learning_rate, iterations)"
   ]
  }
 ],
 "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.12.3"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}