4_1_LogisticRegression.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "LlS4gyYqJVpn"
      },
      "source": [
        "# 로지스틱 회귀"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "tQBDVUjBbIsL"
      },
      "source": [
        "## 럭키백의 확률"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ILi_LPl9JVpw"
      },
      "source": [
        "### 데이터 준비하기"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 206
        },
        "id": "Mba6QeEmLn3r",
        "outputId": "46d85238-2159-4f68-ad30-71be88b0d91d"
      },
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "  Species  Weight  Length  Diagonal   Height   Width\n",
              "0   Bream   242.0    25.4      30.0  11.5200  4.0200\n",
              "1   Bream   290.0    26.3      31.2  12.4800  4.3056\n",
              "2   Bream   340.0    26.5      31.1  12.3778  4.6961\n",
              "3   Bream   363.0    29.0      33.5  12.7300  4.4555\n",
              "4   Bream   430.0    29.0      34.0  12.4440  5.1340"
            ],
            "text/html": [
              "\n",
              "  <div id=\"df-0feb6e8a-b7ff-4ef4-9d81-95bf525845d6\" class=\"colab-df-container\">\n",
              "    <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>Species</th>\n",
              "      <th>Weight</th>\n",
              "      <th>Length</th>\n",
              "      <th>Diagonal</th>\n",
              "      <th>Height</th>\n",
              "      <th>Width</th>\n",
              "    </tr>\n",
              "  </thead>\n",
              "  <tbody>\n",
              "    <tr>\n",
              "      <th>0</th>\n",
              "      <td>Bream</td>\n",
              "      <td>242.0</td>\n",
              "      <td>25.4</td>\n",
              "      <td>30.0</td>\n",
              "      <td>11.5200</td>\n",
              "      <td>4.0200</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>1</th>\n",
              "      <td>Bream</td>\n",
              "      <td>290.0</td>\n",
              "      <td>26.3</td>\n",
              "      <td>31.2</td>\n",
              "      <td>12.4800</td>\n",
              "      <td>4.3056</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>2</th>\n",
              "      <td>Bream</td>\n",
              "      <td>340.0</td>\n",
              "      <td>26.5</td>\n",
              "      <td>31.1</td>\n",
              "      <td>12.3778</td>\n",
              "      <td>4.6961</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>3</th>\n",
              "      <td>Bream</td>\n",
              "      <td>363.0</td>\n",
              "      <td>29.0</td>\n",
              "      <td>33.5</td>\n",
              "      <td>12.7300</td>\n",
              "      <td>4.4555</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>4</th>\n",
              "      <td>Bream</td>\n",
              "      <td>430.0</td>\n",
              "      <td>29.0</td>\n",
              "      <td>34.0</td>\n",
              "      <td>12.4440</td>\n",
              "      <td>5.1340</td>\n",
              "    </tr>\n",
              "  </tbody>\n",
              "</table>\n",
              "</div>\n",
              "    <div class=\"colab-df-buttons\">\n",
              "\n",
              "  <div class=\"colab-df-container\">\n",
              "    <button class=\"colab-df-convert\" onclick=\"convertToInteractive('df-0feb6e8a-b7ff-4ef4-9d81-95bf525845d6')\"\n",
              "            title=\"Convert this dataframe to an interactive table.\"\n",
              "            style=\"display:none;\">\n",
              "\n",
              "  <svg xmlns=\"http://www.w3.org/2000/svg\" height=\"24px\" viewBox=\"0 -960 960 960\">\n",
              "    <path d=\"M120-120v-720h720v720H120Zm60-500h600v-160H180v160Zm220 220h160v-160H400v160Zm0 220h160v-160H400v160ZM180-400h160v-160H180v160Zm440 0h160v-160H620v160ZM180-180h160v-160H180v160Zm440 0h160v-160H620v160Z\"/>\n",
              "  </svg>\n",
              "    </button>\n",
              "\n",
              "  <style>\n",
              "    .colab-df-container {\n",
              "      display:flex;\n",
              "      gap: 12px;\n",
              "    }\n",
              "\n",
              "    .colab-df-convert {\n",
              "      background-color: #E8F0FE;\n",
              "      border: none;\n",
              "      border-radius: 50%;\n",
              "      cursor: pointer;\n",
              "      display: none;\n",
              "      fill: #1967D2;\n",
              "      height: 32px;\n",
              "      padding: 0 0 0 0;\n",
              "      width: 32px;\n",
              "    }\n",
              "\n",
              "    .colab-df-convert:hover {\n",
              "      background-color: #E2EBFA;\n",
              "      box-shadow: 0px 1px 2px rgba(60, 64, 67, 0.3), 0px 1px 3px 1px rgba(60, 64, 67, 0.15);\n",
              "      fill: #174EA6;\n",
              "    }\n",
              "\n",
              "    .colab-df-buttons div {\n",
              "      margin-bottom: 4px;\n",
              "    }\n",
              "\n",
              "    [theme=dark] .colab-df-convert {\n",
              "      background-color: #3B4455;\n",
              "      fill: #D2E3FC;\n",
              "    }\n",
              "\n",
              "    [theme=dark] .colab-df-convert:hover {\n",
              "      background-color: #434B5C;\n",
              "      box-shadow: 0px 1px 3px 1px rgba(0, 0, 0, 0.15);\n",
              "      filter: drop-shadow(0px 1px 2px rgba(0, 0, 0, 0.3));\n",
              "      fill: #FFFFFF;\n",
              "    }\n",
              "  </style>\n",
              "\n",
              "    <script>\n",
              "      const buttonEl =\n",
              "        document.querySelector('#df-0feb6e8a-b7ff-4ef4-9d81-95bf525845d6 button.colab-df-convert');\n",
              "      buttonEl.style.display =\n",
              "        google.colab.kernel.accessAllowed ? 'block' : 'none';\n",
              "\n",
              "      async function convertToInteractive(key) {\n",
              "        const element = document.querySelector('#df-0feb6e8a-b7ff-4ef4-9d81-95bf525845d6');\n",
              "        const dataTable =\n",
              "          await google.colab.kernel.invokeFunction('convertToInteractive',\n",
              "                                                    [key], {});\n",
              "        if (!dataTable) return;\n",
              "\n",
              "        const docLinkHtml = 'Like what you see? Visit the ' +\n",
              "          '<a target=\"_blank\" href=https://colab.research.google.com/notebooks/data_table.ipynb>data table notebook</a>'\n",
              "          + ' to learn more about interactive tables.';\n",
              "        element.innerHTML = '';\n",
              "        dataTable['output_type'] = 'display_data';\n",
              "        await google.colab.output.renderOutput(dataTable, element);\n",
              "        const docLink = document.createElement('div');\n",
              "        docLink.innerHTML = docLinkHtml;\n",
              "        element.appendChild(docLink);\n",
              "      }\n",
              "    </script>\n",
              "  </div>\n",
              "\n",
              "\n",
              "<div id=\"df-7e811f85-fe7c-4ab0-abb9-4ee3d3e75901\">\n",
              "  <button class=\"colab-df-quickchart\" onclick=\"quickchart('df-7e811f85-fe7c-4ab0-abb9-4ee3d3e75901')\"\n",
              "            title=\"Suggest charts\"\n",
              "            style=\"display:none;\">\n",
              "\n",
              "<svg xmlns=\"http://www.w3.org/2000/svg\" height=\"24px\"viewBox=\"0 0 24 24\"\n",
              "     width=\"24px\">\n",
              "    <g>\n",
              "        <path d=\"M19 3H5c-1.1 0-2 .9-2 2v14c0 1.1.9 2 2 2h14c1.1 0 2-.9 2-2V5c0-1.1-.9-2-2-2zM9 17H7v-7h2v7zm4 0h-2V7h2v10zm4 0h-2v-4h2v4z\"/>\n",
              "    </g>\n",
              "</svg>\n",
              "  </button>\n",
              "\n",
              "<style>\n",
              "  .colab-df-quickchart {\n",
              "      --bg-color: #E8F0FE;\n",
              "      --fill-color: #1967D2;\n",
              "      --hover-bg-color: #E2EBFA;\n",
              "      --hover-fill-color: #174EA6;\n",
              "      --disabled-fill-color: #AAA;\n",
              "      --disabled-bg-color: #DDD;\n",
              "  }\n",
              "\n",
              "  [theme=dark] .colab-df-quickchart {\n",
              "      --bg-color: #3B4455;\n",
              "      --fill-color: #D2E3FC;\n",
              "      --hover-bg-color: #434B5C;\n",
              "      --hover-fill-color: #FFFFFF;\n",
              "      --disabled-bg-color: #3B4455;\n",
              "      --disabled-fill-color: #666;\n",
              "  }\n",
              "\n",
              "  .colab-df-quickchart {\n",
              "    background-color: var(--bg-color);\n",
              "    border: none;\n",
              "    border-radius: 50%;\n",
              "    cursor: pointer;\n",
              "    display: none;\n",
              "    fill: var(--fill-color);\n",
              "    height: 32px;\n",
              "    padding: 0;\n",
              "    width: 32px;\n",
              "  }\n",
              "\n",
              "  .colab-df-quickchart:hover {\n",
              "    background-color: var(--hover-bg-color);\n",
              "    box-shadow: 0 1px 2px rgba(60, 64, 67, 0.3), 0 1px 3px 1px rgba(60, 64, 67, 0.15);\n",
              "    fill: var(--button-hover-fill-color);\n",
              "  }\n",
              "\n",
              "  .colab-df-quickchart-complete:disabled,\n",
              "  .colab-df-quickchart-complete:disabled:hover {\n",
              "    background-color: var(--disabled-bg-color);\n",
              "    fill: var(--disabled-fill-color);\n",
              "    box-shadow: none;\n",
              "  }\n",
              "\n",
              "  .colab-df-spinner {\n",
              "    border: 2px solid var(--fill-color);\n",
              "    border-color: transparent;\n",
              "    border-bottom-color: var(--fill-color);\n",
              "    animation:\n",
              "      spin 1s steps(1) infinite;\n",
              "  }\n",
              "\n",
              "  @keyframes spin {\n",
              "    0% {\n",
              "      border-color: transparent;\n",
              "      border-bottom-color: var(--fill-color);\n",
              "      border-left-color: var(--fill-color);\n",
              "    }\n",
              "    20% {\n",
              "      border-color: transparent;\n",
              "      border-left-color: var(--fill-color);\n",
              "      border-top-color: var(--fill-color);\n",
              "    }\n",
              "    30% {\n",
              "      border-color: transparent;\n",
              "      border-left-color: var(--fill-color);\n",
              "      border-top-color: var(--fill-color);\n",
              "      border-right-color: var(--fill-color);\n",
              "    }\n",
              "    40% {\n",
              "      border-color: transparent;\n",
              "      border-right-color: var(--fill-color);\n",
              "      border-top-color: var(--fill-color);\n",
              "    }\n",
              "    60% {\n",
              "      border-color: transparent;\n",
              "      border-right-color: var(--fill-color);\n",
              "    }\n",
              "    80% {\n",
              "      border-color: transparent;\n",
              "      border-right-color: var(--fill-color);\n",
              "      border-bottom-color: var(--fill-color);\n",
              "    }\n",
              "    90% {\n",
              "      border-color: transparent;\n",
              "      border-bottom-color: var(--fill-color);\n",
              "    }\n",
              "  }\n",
              "</style>\n",
              "\n",
              "  <script>\n",
              "    async function quickchart(key) {\n",
              "      const quickchartButtonEl =\n",
              "        document.querySelector('#' + key + ' button');\n",
              "      quickchartButtonEl.disabled = true;  // To prevent multiple clicks.\n",
              "      quickchartButtonEl.classList.add('colab-df-spinner');\n",
              "      try {\n",
              "        const charts = await google.colab.kernel.invokeFunction(\n",
              "            'suggestCharts', [key], {});\n",
              "      } catch (error) {\n",
              "        console.error('Error during call to suggestCharts:', error);\n",
              "      }\n",
              "      quickchartButtonEl.classList.remove('colab-df-spinner');\n",
              "      quickchartButtonEl.classList.add('colab-df-quickchart-complete');\n",
              "    }\n",
              "    (() => {\n",
              "      let quickchartButtonEl =\n",
              "        document.querySelector('#df-7e811f85-fe7c-4ab0-abb9-4ee3d3e75901 button');\n",
              "      quickchartButtonEl.style.display =\n",
              "        google.colab.kernel.accessAllowed ? 'block' : 'none';\n",
              "    })();\n",
              "  </script>\n",
              "</div>\n",
              "    </div>\n",
              "  </div>\n"
            ],
            "application/vnd.google.colaboratory.intrinsic+json": {
              "type": "dataframe",
              "variable_name": "fish",
              "summary": "{\n  \"name\": \"fish\",\n  \"rows\": 159,\n  \"fields\": [\n    {\n      \"column\": \"Species\",\n      \"properties\": {\n        \"dtype\": \"category\",\n        \"num_unique_values\": 7,\n        \"samples\": [\n          \"Bream\",\n          \"Roach\",\n          \"Pike\"\n        ],\n        \"semantic_type\": \"\",\n        \"description\": \"\"\n      }\n    },\n    {\n      \"column\": \"Weight\",\n      \"properties\": {\n        \"dtype\": \"number\",\n        \"std\": 357.9783165508931,\n        \"min\": 0.0,\n        \"max\": 1650.0,\n        \"num_unique_values\": 101,\n        \"samples\": [\n          770.0,\n          51.5,\n          197.0\n        ],\n        \"semantic_type\": \"\",\n        \"description\": \"\"\n      }\n    },\n    {\n      \"column\": \"Length\",\n      \"properties\": {\n        \"dtype\": \"number\",\n        \"std\": 10.716328098884247,\n        \"min\": 8.4,\n        \"max\": 63.4,\n        \"num_unique_values\": 93,\n        \"samples\": [\n          14.7,\n          18.8,\n          19.6\n        ],\n        \"semantic_type\": \"\",\n        \"description\": \"\"\n      }\n    },\n    {\n      \"column\": \"Diagonal\",\n      \"properties\": {\n        \"dtype\": \"number\",\n        \"std\": 11.610245832690964,\n        \"min\": 8.8,\n        \"max\": 68.0,\n        \"num_unique_values\": 124,\n        \"samples\": [\n          39.2,\n          27.2,\n          23.1\n        ],\n        \"semantic_type\": \"\",\n        \"description\": \"\"\n      }\n    },\n    {\n      \"column\": \"Height\",\n      \"properties\": {\n        \"dtype\": \"number\",\n        \"std\": 4.2862076199688675,\n        \"min\": 1.7284,\n        \"max\": 18.957,\n        \"num_unique_values\": 154,\n        \"samples\": [\n          15.438,\n          7.293,\n          2.8728\n        ],\n        \"semantic_type\": \"\",\n        \"description\": \"\"\n      }\n    },\n    {\n      \"column\": \"Width\",\n      \"properties\": {\n        \"dtype\": \"number\",\n        \"std\": 1.685803869992167,\n        \"min\": 1.0476,\n        \"max\": 8.142,\n        \"num_unique_values\": 152,\n        \"samples\": [\n          3.1571,\n          1.3936,\n          3.6835\n        ],\n        \"semantic_type\": \"\",\n        \"description\": \"\"\n      }\n    }\n  ]\n}"
            }
          },
          "metadata": {},
          "execution_count": 2
        }
      ],
      "source": [
        "import pandas as pd\n",
        "\n",
        "fish = pd.read_csv('https://bit.ly/fish_csv_data')\n",
        "fish.head()  # dataframe의 처음 5개 행만 출력"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "UWJWlRCHVWUg",
        "outputId": "451a9987-5f9b-4cb6-e04a-df7abaad8ac0"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Bream' 'Roach' 'Whitefish' 'Parkki' 'Perch' 'Pike' 'Smelt']\n"
          ]
        }
      ],
      "source": [
        "# Species 열에서 고유한 값만 출력\n",
        "print(pd.unique(fish['Species']))"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "# dataframe에서 target을 제외한 5개 열 선택하여 2차원 numpy array로 만들기\n",
        "fish_input = fish[['Weight','Length','Diagonal','Height','Width']].to_numpy()\n",
        "\n",
        "# dataframe에서 target 열 선택하여 1차원 numpy array로 만들기\n",
        "fish_target = fish['Species'].to_numpy()"
      ],
      "metadata": {
        "id": "jaaNAlmXc_UB"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# train set, test set 나누기\n",
        "from sklearn.model_selection import train_test_split\n",
        "train_input, test_input, train_target, test_target = train_test_split(fish_input, fish_target, random_state=42)"
      ],
      "metadata": {
        "id": "R9pwGWR1dJHO"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "# 표준화\n",
        "from sklearn.preprocessing import StandardScaler\n",
        "ss = StandardScaler()\n",
        "ss.fit(train_input)\n",
        "\n",
        "train_scaled = ss.transform(train_input)\n",
        "test_scaled = ss.transform(test_input)"
      ],
      "metadata": {
        "id": "porDOI-geW9a"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oAxk-V5kQcgc"
      },
      "source": [
        "### k-최근접 이웃 분류기의 확률 예측"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "# k-최근접 이웃 분류기\n",
        "from sklearn.neighbors import KNeighborsClassifier\n",
        "\n",
        "# 이웃 개수 3\n",
        "kn = KNeighborsClassifier(n_neighbors=3)\n",
        "\n",
        "# train\n",
        "kn.fit(train_scaled, train_target)\n",
        "\n",
        "# score\n",
        "print(kn.score(train_scaled, train_target))\n",
        "print(kn.score(test_scaled, test_target))"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "yPzkbp0dg_e1",
        "outputId": "d6bbba65-c932-4db7-be22-c8179fa4ec99"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "0.8907563025210085\n",
            "0.85\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "a879-O42RhFO",
        "outputId": "a9a77272-d6bd-4940-be68-6cf168ca682d"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Bream' 'Parkki' 'Perch' 'Pike' 'Roach' 'Smelt' 'Whitefish']\n"
          ]
        }
      ],
      "source": [
        "print(kn.classes_)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "EucmtF8HVOS_",
        "outputId": "564c99ea-37f1-46d4-fff1-db1b7c614a0b"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Perch' 'Smelt' 'Pike' 'Perch' 'Perch']\n"
          ]
        }
      ],
      "source": [
        "# test set의 처음 5개 샘플의 target 값 예측\n",
        "print(kn.predict(test_scaled[:5]))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "OSDr8WSKXbUa",
        "outputId": "0683a443-bc0c-4987-f8e2-e1bae01df230"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[0.     0.     1.     0.     0.     0.     0.    ]\n",
            " [0.     0.     0.     0.     0.     1.     0.    ]\n",
            " [0.     0.     0.     1.     0.     0.     0.    ]\n",
            " [0.     0.     0.6667 0.     0.3333 0.     0.    ]\n",
            " [0.     0.     0.6667 0.     0.3333 0.     0.    ]]\n"
          ]
        }
      ],
      "source": [
        "import numpy as np\n",
        "\n",
        "# 처음 5개 샘플의 클래스별 확률값\n",
        "proba = kn.predict_proba(test_scaled[:5])\n",
        "# 소수점 4자리까지 출력 (반올림)\n",
        "print(np.round(proba, decimals=4))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "Hk-ywsfKkf7t",
        "outputId": "2ec40948-c3f8-4f03-f41e-9e56773bbd49"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[['Roach' 'Perch' 'Perch']]\n"
          ]
        }
      ],
      "source": [
        "# test_scaled[3:4]: 4번째 샘플을 2차원 배열로 만듦\n",
        "# kneighbors 메서드의 인자는 항상 2차원 배열\n",
        "distances, indexes = kn.kneighbors(test_scaled[3:4])\n",
        "print(train_target[indexes])"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Q9_wuI_0tEqL"
      },
      "source": [
        "## 로지스틱 회귀"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 449
        },
        "id": "8rdDSaZ5uji2",
        "outputId": "4cfb37b9-2f96-4665-ba0c-102513f99f42"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 640x480 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA9y0lEQVR4nO3deXRU5eHG8Wcmy2QPhGxAAmHfIbIFXMFG44Z1p2gFqdJqkYr5tQouUK0aV6QqlYpraym4VK2CKKK4FAQE2WQNW0IgO8lknUlm7u+P1GjKYoAkd5bv55w5SW7uJc/MgcnDve99X4thGIYAAAB8hNXsAAAAAC2JcgMAAHwK5QYAAPgUyg0AAPAplBsAAOBTKDcAAMCnUG4AAIBPCTQ7QFtzu906dOiQIiMjZbFYzI4DAACawTAMVVRUqFOnTrJaT3xuxu/KzaFDh5ScnGx2DAAAcApyc3OVlJR0wn38rtxERkZKanhxoqKiTE4DAACaw263Kzk5ufH3+In4Xbn5/lJUVFQU5QYAAC/TnCElDCgGAAA+hXIDAAB8CuUGAAD4FMoNAADwKZQbAADgUyg3AADAp1BuAACAT6HcAAAAn2Jqufniiy80btw4derUSRaLRe++++5PHrNy5UoNHTpUNptNPXv21KuvvtrqOQEAgPcwtdxUVVVpyJAhmjdvXrP237dvny699FKNHTtWGzdu1PTp03XLLbfoo48+auWkAADAW5i6/MLFF1+siy++uNn7z58/X926ddNTTz0lSerXr5+++uorPf3008rIyDjmMQ6HQw6Ho/Fru91+eqEBAIBH86oxN6tXr1Z6enqTbRkZGVq9evVxj8nKylJ0dHTjgxXBAQDwbV5VbvLz85WQkNBkW0JCgux2u2pqao55zMyZM1VeXt74yM3NbYuoAAD4FbfbUKWjXofLa3So7Ni/k9uKz68KbrPZZLPZzI4BAIDHM4yGglJWXafSKqeOVDtVXlOnsur/PmqcstfUq7ymTvbaOtlrGh4VjnpVOuplGA1/Tlq3GC3+zWjTnodXlZvExEQVFBQ02VZQUKCoqCiFhoaalAoAAM/lchsqrXKqqMKhwopaFVU4VFTpUEmlUyWVDpVUOVX838+PVDtV5zJO6+cFWi2yWFoo/KlmMPfHn5zRo0dr6dKlTbYtX75co0eb1w4BADCLy22owF6rQ2U1OlReq8NlDZeE8u21KrA7VGCvVWGFQy73yRWW0KAAxYQHKzo0SO3DgxQdGqTo0Iav24UFKSokSFGhgf/9GKTIkMCGhy1IIUFWWUxuN6aWm8rKSmVnZzd+vW/fPm3cuFExMTHq0qWLZs6cqby8PP3tb3+TJN1666167rnndNddd+lXv/qVPv30U73xxhtasmSJWU8BAIBWYxiGyqrrdKC0WgdKqnSgpFoHj1Qrt7RGB8uqdbisVvXNKC4Wi9QhPFhxkSGKi7QpLsKm2IhgxUbY1CEiWDHhweoQ3vB5+7BghQYHtMGzaz2mlptvvvlGY8eObfw6MzNTkjRp0iS9+uqrOnz4sHJychq/361bNy1ZskR33nmn/vznPyspKUkvvvjicW8DBwDAGzjqXdpXXKW9RVXaU1ipvcVV2ltUqX3FVbLX1p/w2KAAixKjQ9QxOlSdokPUsV3Dx/ioECVEhSgxKkSxEcEKDPCqe4hOi8UwjNO7uOZl7Ha7oqOjVV5erqioKLPjAAD8iNtt6EBptbYftmtHfoV2F1RoV0GF9pdUn/DSUUKUTV1jwtWlQ5i6xIQpOSZUSe3DlNw+TPGRNlmtJg9yaQMn8/vbq8bcAADgLepcbu0qqNCWg+XaklfeWGiqna5j7h9pC1SP+Ah1jwtXj7gIdY8NV7e4cHWNCff6y0RtjXIDAMBpMgxDB4/UaEPOEX2bU6aNuWXadtguZ737qH1tgVb1SYxU38RI9U6IVK+ESPVOiFBiVIjpA3F9BeUGAICT5HIb2nbIrjX7SrRuf6k25JSpqMJx1H6RIYEanBStgZ2jNaBTtPp3jFRKh3C/Gv9iBsoNAAA/we02tO2wXf/JLtbXe0v0zf4jqnA0HegbFGBR/07RGtqlnVKT22lIUjt17RDG2RgTUG4AADiG/PJafbGrSF9mF2tVdrFKqpxNvh9pC9SIbjEa2S1Gw7u218DO0QoJYmyMJ6DcAACghktNmw6W6bMdhVqxvVDbDtubfD88OECjunfQ6B4dNKp7B/XrGKUAP7hLyRtRbgAAfstZ79bqvSVatvWwlm8rUHHlD2dnLBZpSFI7ndsrVmf3ilNqcjsFBzJWxhtQbgAAfsVZ79aXu4u0ZPNhfbK9oMkkeZG2QJ3bO07n943XmD5x6hDBwsveiHIDAPB5brehbw4c0Xsb87Rky2GVVdc1fi82wqaMAQm6aGCiRnXvoCDuZPJ6lBsAgM/KLa3Wm+sP6u31B5VXVtO4PS7SpssGd9QlgzpqaJf2jJ3xMZQbAIBPqa1z6eNtBXpjXa7+s6dY3y8yFGEL1EUDE3VFameN7tGBQuPDKDcAAJ+QV1ajf3x9QIvW5ar0R7dtn9Wzg64bnqyMAYncqu0nKDcAAK9lGIZW7ynRa6v3a/m2An2/9mRiVIiuG56ka4cnKzkmzNyQaHOUGwCA16l3ubVky2HN/3yvtv9oPpoze3TQxNEpSu8XzxIHfoxyAwDwGtXOer2xLlcLvtzXOEA4NChA1wxL0sTRXdUrIdLkhPAElBsAgMerctTr718f0Atf7G0cT9MhPFg3nZmiX47qqvbhwSYnhCeh3AAAPFa1s15/X31Af/1RqenaIUxTzumua4YlMUAYx0S5AQB4HGe9W/9Yc0DPfZrduGBl1w5h+t35vfTz1E6Mp8EJUW4AAB7DMAwt2XJYjy/bqZzSaklSl5gwTTu/p648ozOlBs1CuQEAeIQ1e0v0yIc7tCm3TFLDLMLT03vpuuHJLImAk0K5AQCY6nB5jR5ZukPvbzokSQoLDtBvzu2hW87ppnAbv6Zw8vhbAwAwhbPerZe+2qdnP92taqdLVov0i5FddGd6b8VFsho3Th3lBgDQ5v6TXaz7392qvcVVkqRhXdvrgcsHaGDnaJOTwRdQbgAAbaa8uk4PLdmmN9cflCTFRth0zyV9deUZnWWxsJAlWgblBgDQJj7ccliz/v2diiocslikG0d11e8z+igqJMjsaPAxlBsAQKsqrnTovne2atl3+ZKkHnHhevyawRrWNcbkZPBVlBsAQKv5ZFuB7n57s0qqnAq0WnTbmB6aOrYnMwujVVFuAAAtrspRr4eWbNc/1+ZIkvomRmrOdanq3ynK5GTwB5QbAECL2phbpumLvtX+kmpZLNItZ3fT/13Yh7M1aDOUGwBAizAMQy99tU+PfrhD9W5DnaJD9OR1Q3Rmj1izo8HPUG4AAKetvKZOd721SR99VyBJunRQRz1y1SBFh3InFNoe5QYAcFq2HCzXbxeuV25pjYICLLr/sv66cVRX5q2BaSg3AIBTtnhdju5/9zs5XW4ltQ/VX24YqsFJ7cyOBT9HuQEAnLQ6l1sPL9muV1ftlyRd0D9BT14zRNFhXIaC+Sg3AICTcqTKqakLN2jVnhJJ0p3pvTXt/J6yWrkMBc9AuQEANNvO/Ard8rd1yi2tUVhwgOZcl6qLBiaaHQtognIDAGiWL3cX6bbXN6jSUa/kmFAtmDhcfROZlA+eh3IDAPhJb36Tq5n/2qJ6t6G0bjGa/8thah8ebHYs4JgoNwCA4zIMQ8+syNbTn+ySJP08tZMev2awbIHMNgzPRbkBABxTncut+97ZqsXf5EqSbhvTQ3+4sA8Dh+HxKDcAgKPU1rl0+8IN+mR7oawW6cGfD9QvR3U1OxbQLJQbAEATVY56TfnbN1q1p0S2QKvmXT9U6f0TzI4FNBvlBgDQqLy6Tje9ulbf5pQpPDhAL900QqO6dzA7FnBSKDcAAElSUYVDN760RjvyKxQdGqTXfjVSqcntzI4FnDTKDQBABfZaTXjha+0trlJshE2v3zKSOWzgtSg3AODnCu21mrCgodh0bheq129JU7fYcLNjAaeMcgMAfqyowtFQbIoais2iX49SckyY2bGA02I1OwAAwBzFlQ5dv+Br7SmqUsfoEP1zCsUGvoFyAwB+qKTSoRsWrNHuwkolRjUUmy4dKDbwDZQbAPAzFbV1mvTKWu0sqFB8pE3//PUopTDGBj6EcgMAfqS2zqVbXvtGW/Ps6hAerIVTRjF4GD6HcgMAfqLe5dbtC7/Vmn2lirAF6rVfjVTP+AizYwEtjnIDAH7A7TZ099tb9Mn2AgUHWvXipOEa2Dna7FhAq6DcAICPMwxDDy/drrc3HFSA1aJ51w9lSQX4NMoNAPi4l77ap5e+2idJevzqwbqARTDh4yg3AODDPtxyWA8v3S5JuueSvrp6WJLJiYDWR7kBAB+1IeeIpi/eKMOQbhzVVVPO6W52JKBNmF5u5s2bp5SUFIWEhCgtLU1r16494f5z585Vnz59FBoaquTkZN15552qra1to7QA4B0OlFRpymvfyFHv1vl94zV7XH9ZLBazYwFtwtRys3jxYmVmZmr27NnasGGDhgwZooyMDBUWFh5z/4ULF2rGjBmaPXu2tm/frpdeekmLFy/WPffc08bJAcBzHalyavIr61RS5dTAzlF6dsIZCgww/f+yQJsx9W/7nDlzNGXKFE2ePFn9+/fX/PnzFRYWppdffvmY+69atUpnnXWWrr/+eqWkpOjCCy/UhAkTfvJsDwD4izqXW7e+vr5xhe+XJ41QuI01kuFfTCs3TqdT69evV3p6+g9hrFalp6dr9erVxzzmzDPP1Pr16xvLzN69e7V06VJdcsklx/05DodDdru9yQMAfNUD73/XOEnfyzeNUHxUiNmRgDZnWp0vLi6Wy+VSQkLTWxITEhK0Y8eOYx5z/fXXq7i4WGeffbYMw1B9fb1uvfXWE16WysrK0gMPPNCi2QHAE/1jzQG9/nWOLBbpz79IVZ/ESLMjAabwqouwK1eu1COPPKK//OUv2rBhg/71r39pyZIl+tOf/nTcY2bOnKny8vLGR25ubhsmBoC2sWZviWa/950k6fcX9tHP+jGXDfyXaWduYmNjFRAQoIKCgibbCwoKlJiYeMxj7r//ft1444265ZZbJEmDBg1SVVWVfv3rX+vee++V1Xp0V7PZbLLZbC3/BADAQ+SWVuu2f2xQvdvQuCGd9NsxPcyOBJjKtDM3wcHBGjZsmFasWNG4ze12a8WKFRo9evQxj6murj6qwAQEBEhqmF4cAPxNtbNev/77epX+986ox68ezC3f8HumDqHPzMzUpEmTNHz4cI0cOVJz585VVVWVJk+eLEmaOHGiOnfurKysLEnSuHHjNGfOHJ1xxhlKS0tTdna27r//fo0bN66x5ACAvzAMQ/e+s1XbD9sVGxGsF24crtBg3gsBU8vN+PHjVVRUpFmzZik/P1+pqalatmxZ4yDjnJycJmdq7rvvPlksFt13333Ky8tTXFycxo0bp4cfftispwAApnl9TY7e+TavcTHMTu1CzY4EeASL4WfXc+x2u6Kjo1VeXq6oqCiz4wDAKdmYW6Zr569SncvQPZf01a/PZZwNfNvJ/P72qrulAABSaZVTU/+xQXUuQxcNSGTNKOB/UG4AwIu43IamL96ovLIadYsN1+PXMoAY+F+UGwDwIs9+ultf7CpSSJBVz/9yqKJCgsyOBHgcyg0AeInVe0r0zIrdkqRHrhykvomMGwSOhXIDAF6gtMqp6Yu/lduQrhuepKuGJpkdCfBYlBsA8HCGYeiutzapwO5Qj7hw/fHyAWZHAjwa5QYAPNyrq/brk+2FCg606tkJQxUWbOoUZYDHo9wAgAfbmleurKU7JEn3XtJP/Tsxzgb4KZQbAPBQVY56Tfvnt3K63Lqgf4Imju5qdiTAK1BuAMBD/emDbdpXXKWO0SF64hrmswGai3IDAB5o+bYCLVqXK4tFmnNdqtqFBZsdCfAalBsA8DDFlQ7NeHuzJGnKOd01ukcHkxMB3oVyAwAexDAMzXh7s0qqnOqbGKn/u7C32ZEAr0O5AQAPsnhdbsNt3wFWPT0+VbbAALMjAV6HcgMAHuJASZUe/GCbJOn3Gb3VryO3fQOngnIDAB7A5TaU+cYmVTtdSusWo5vP7m52JMBrUW4AwAO88p99Wn/giCJsgXrquiEKsHLbN3CqKDcAYLK9RZV64qOdkqT7Lu2npPZhJicCvBvlBgBM5HIb+sNbm+Wod+ucXrEaPyLZ7EiA16PcAICJfnw56tGrmYUYaAmUGwAwyb7iqsbLUfdc0k+d24WanAjwDZQbADCBy23oD29ukqPerbN7xmrCSC5HAS2FcgMAJvjb6v365sARhQcHKOuqQVyOAloQ5QYA2lheWU3j5agZF/dVcgx3RwEtiXIDAG3IMAzd984WVTtdGt61vW5I62p2JMDnUG4AoA29v/mwPttZpOAAq7KuGiQrk/UBLY5yAwBtpKzaqQff/06S9NuxPdQrIdLkRIBvotwAQBt5eMl2FVc61TM+QreN6WF2HMBnUW4AoA2syi7Wm+sPSpIevWqQbIEBJicCfBflBgBaWW2dS/e8s0WS9MtRXTQ8JcbkRIBvo9wAQCv7y8o92l9SrYQom+66qK/ZcQCfR7kBgFa0t6hS81fukSTNumyAokKCTE4E+D7KDQC0EsMwNOu97+R0uXVu7zhdMijR7EiAX6DcAEAreX/zYX2VXazgQKsevHwASywAbYRyAwCtwF5bpz99sE2SdPvYnkqJDTc5EeA/KDcA0ArmfLxLRRUOdY8N12/O6252HMCvUG4AoIVtzSvX31bvlyT96YqBzGkDtDHKDQC0ILfb0H3vbpXbkC4f0kln9Yw1OxLgdyg3ANCC3tpwUBtzyxQeHKD7Lu1ndhzAL1FuAKCFlNfU6bEPd0iSpqf3VnxUiMmJAP9EuQGAFjL3k10qqWpYGPOms1LMjgP4LcoNALSAHfl2/W31AUnSH8cNUFAAb6+AWfjXBwCnyTAMzX7vO7nchi4emKizezGIGDAT5QYATtP7mw9rzb5ShQRZdS+DiAHTUW4A4DRUOer1yJLtkqSpY3oqqX2YyYkAUG4A4DQ8v3KP8u216hITpinnMhMx4AkoNwBwinJLq/XCl3slSfde2k8hQcxEDHgCyg0AnKKsD7fLWe/WmT066ML+CWbHAfBflBsAOAVf7y3R0i35slqkWeP6y2KxmB0JwH9RbgDgJLnchh54f5sk6fq0LuqbGGVyIgA/RrkBgJO0eF2uth+2KyokUJkX9DE7DoD/QbkBgJNQXlOnpz7eKalh/aiY8GCTEwH4X5QbADgJz67YrZIqp3rEhevG0V3NjgPgGCg3ANBM+4ur9Nrq/ZKk+y7rz/pRgIfiXyYANNNjy3aozmXo3N5xGtsn3uw4AI6DcgMAzbB2X6k+3Npw6/e9l7B+FODJKDcA8BPcbkMPLWm49fsXI7uoT2KkyYkAnIjp5WbevHlKSUlRSEiI0tLStHbt2hPuX1ZWpqlTp6pjx46y2Wzq3bu3li5d2kZpAfij9zblafPBckXYAnVnem+z4wD4CYFm/vDFixcrMzNT8+fPV1pamubOnauMjAzt3LlT8fFHX892Op264IILFB8fr7feekudO3fWgQMH1K5du7YPD8Av1DhdenxZw63fvx3bQ3GRNpMTAfgpppabOXPmaMqUKZo8ebIkaf78+VqyZIlefvllzZgx46j9X375ZZWWlmrVqlUKCgqSJKWkpLRlZAB+5sUv9+pwea06twvVr87qZnYcAM1g2mUpp9Op9evXKz09/YcwVqvS09O1evXqYx7z73//W6NHj9bUqVOVkJCggQMH6pFHHpHL5Truz3E4HLLb7U0eANAchfZaPf/5HknS3Rf3ZdVvwEuYVm6Ki4vlcrmUkNB0Jd2EhATl5+cf85i9e/fqrbfeksvl0tKlS3X//ffrqaee0kMPPXTcn5OVlaXo6OjGR3Jycos+DwC+6+lPdqna6VJqcjuNG9zR7DgAmsn0AcUnw+12Kz4+Xi+88IKGDRum8ePH695779X8+fOPe8zMmTNVXl7e+MjNzW3DxAC81a6CCi1e1/B+cd+l/Vj1G/Aipo25iY2NVUBAgAoKCppsLygoUGJi4jGP6dixo4KCghQQ8MOp4X79+ik/P19Op1PBwUev8WKz2WSzMQAQwMl59MMdchtSxoAEDU+JMTsOgJNg2pmb4OBgDRs2TCtWrGjc5na7tWLFCo0ePfqYx5x11lnKzs6W2+1u3LZr1y517NjxmMUGAE7Fqj3F+nRHoQKtFt19UV+z4wA4SaZelsrMzNSCBQv02muvafv27brttttUVVXVePfUxIkTNXPmzMb9b7vtNpWWluqOO+7Qrl27tGTJEj3yyCOaOnWqWU8BgI9xuw09snS7JOmGtC7qHhdhciIAJ8vUW8HHjx+voqIizZo1S/n5+UpNTdWyZcsaBxnn5OTIav2hfyUnJ+ujjz7SnXfeqcGDB6tz58664447dPfdd5v1FAD4mPc25Wlrnl2RtkD97me9zI4D4BRYDMMwzA7Rlux2u6Kjo1VeXq6oqCiz4wDwILV1Lv3sqc+VV1ajP2T00dSxPc2OBOC/Tub3t1fdLQUArenVVfuVV1ajjtEhuvlsJuwDvBXlBgAkHalyat5n2ZKk/7uwDxP2AV6McgMAkuZ9lq2K2nr1TYzUlWd0NjsOgNNAuQHg93JLq/W31QckSTMu7qsAKxP2Ad6McgPA781ZvktOl1tn9uig83rHmR0HwGmi3ADwa98dKte7G/MkSTMvZpkFwBdQbgD4tUc/3CHDkMYN6aRBSdFmxwHQAig3APzWV7uL9eXuYgUFWPSHC/uYHQdAC6HcAPBLbrehrA+/X2ahq7p0CDM5EYCWQrkB4Jfe33xI3x2yK8IWqGnnMxMx4EsoNwD8jrPerSc/3ilJ+s253dUhwmZyIgAtiXIDwO8sXHNAuaU1iou06eZzWGYB8DWUGwB+pdJRr2c/bVhm4Y6f9VJYcKDJiQC0NMoNAL+y4Iu9KqlyqltsuMaPSDY7DoBWQLkB4DeKKhxa8OVeSdIfMvooKIC3QMAX8S8bgN949tPdqna6NCQpWhcPTDQ7DoBWQrkB4BcOlFRp4ZocSdLdF/dlmQXAh1FuAPiFJz/epXq3ofN6x+nMHrFmxwHQiig3AHze1rxyvb/pkCTprotYZgHwdZQbAD7vsWU7JEk/T+2kAZ1YHBPwdZQbAD5tVfYPi2P+3wWctQH8QbNnr9q8ebMGDhwoq9WqzZs3n3DfwYMHn3YwADhdhmE0nrW5fmQXFscE/ESzy01qaqry8/MVHx+v1NRUWSwWGYbR+P3vv7ZYLHK5XK0SFgBOxodb87XpYLnCggN0+/m9zI4DoI00u9zs27dPcXFxjZ8DgCerd7n15EcNi2Peck53xUWyOCbgL5pdbrp27XrMzwHAE73xzUHtLa5STHiwprA4JuBXTnnFuN27d+uzzz5TYWGh3G53k+/NmjXrtIMBwKmqcbo095NdkqTbx/ZUZEiQyYkAtKVTKjcLFizQbbfdptjYWCUmJjaZ6dNisVBuAJjqlVX7VFjhUFL7UN0wqovZcQC0sVMqNw899JAefvhh3X333S2dBwBOS1m1U8+v3CNJyrygt2yBASYnAtDWTmmemyNHjujaa69t6SwAcNqeX7lHFbX16psYqZ+ndjY7DgATnFK5ufbaa/Xxxx+3dBYAOC2Hy2v06qr9khqWWQiwsjgm4I+afVnqmWeeafy8Z8+euv/++/X1119r0KBBCgpqOljvd7/7XcslBIBmmrt8txz1bo1MidHYPvFmxwFgEovx45n4TqBbt+bdSmmxWLR3797TCtWa7Ha7oqOjVV5erqioKLPjAGgh2YUVuvDpL+Q2pLdvO1PDurY3OxKAFnQyv79PahK/Y/m+G/34jikAaGtPfLRTbkO6oH8CxQbwc6e8cOZLL72kgQMHKiQkRCEhIRo4cKBefPHFlswGAM2yIeeIPvquQFaLdFcGi2MC/u6UbgWfNWuW5syZo2nTpmn06NGSpNWrV+vOO+9UTk6OHnzwwRYNCQDHYxiGHvuwYXHMq4cmqVdCpMmJAJit2WNufiwuLk7PPPOMJkyY0GT7P//5T02bNk3FxcUtFrClMeYG8C2f7SzU5FfWKTjQqs9+P0ad24WaHQlAKziZ39+ndFmqrq5Ow4cPP2r7sGHDVF9ffyp/JACcNLf7h7M2k0Z3pdgAkHSK5ebGG2/U888/f9T2F154QTfccMNphwKA5nhvU5525FcoMiRQvx3T0+w4ADzEKS+c+dJLL+njjz/WqFGjJElr1qxRTk6OJk6cqMzMzMb95syZc/opAeB/OOpdeurjhsUxbxvTQ+3Dg01OBMBTnFK52bp1q4YOHSpJ2rOnYQ2X2NhYxcbGauvWrY37cXs4gNbyj69zdPBIjRKibJp8ZvPm4QLgH06p3Hz22WctnQMAmq2itk7PfZYtSZqe3luhwSyOCeAHpzzPDQCYZcEXe1Va5VT3uHBdOyzJ7DgAPAzlBoBXKayo1YIvG2ZMvyujjwIDeBsD0BTvCgC8yjMrdqumzqXU5HbKGJBodhwAHohyA8Br7C2q1D/X5kqSZlzcl5sWABwT5QaA13jio51yuQ39rG+8RnXvYHYcAB6KcgPAK2zIOaIPt+bLapHuvriv2XEAeDDKDQCPZxiGHl3asMzCNcOS1JvFMQGcAOUGgMf7ZHuh1u4vlS3Qqjsv6G12HAAejnIDwKPVu9x6bFnDWZtfnd1NHaNZHBPAiVFuAHi0t9YfVHZhpdqFBenW83qYHQeAF6DcAPBYNU6Xnv6kYXHM28f2VHRokMmJAHgDyg0Aj/Xil3tVYHcoqX2obhzd1ew4ALwE5QaARyqqcGj+53skSX/I6CNbIItjAmgeyg0Aj/TnFbtU5XRpcFK0xg3uZHYcAF6EcgPA42QX/rDMwj2X9JPVyjILAJqPcgPA4zy2bIdcbkPp/VhmAcDJ84hyM2/ePKWkpCgkJERpaWlau3Zts45btGiRLBaLrrjiitYNCKDNrN1XquXbChRgtWgGyywAOAWml5vFixcrMzNTs2fP1oYNGzRkyBBlZGSosLDwhMft379fv//973XOOee0UVIArc0wDD28dLsk6RcjktUznmUWAJw808vNnDlzNGXKFE2ePFn9+/fX/PnzFRYWppdffvm4x7hcLt1www164IEH1L179zZMC6A1fbD5sDbllik8OEDT01lmAcCpMbXcOJ1OrV+/Xunp6Y3brFar0tPTtXr16uMe9+CDDyo+Pl4333zzT/4Mh8Mhu93e5AHA89TWuRqXWfjNeT0UF2kzOREAb2VquSkuLpbL5VJCQkKT7QkJCcrPzz/mMV999ZVeeuklLViwoFk/IysrS9HR0Y2P5OTk084NoOW9umq/Dh6pUWJUiG45p5vZcQB4MdMvS52MiooK3XjjjVqwYIFiY2ObdczMmTNVXl7e+MjNzW3llABOVkmlQ/M+zZbUMGFfWHCgyYkAeDNT30FiY2MVEBCggoKCJtsLCgqUmJh41P579uzR/v37NW7cuMZtbrdbkhQYGKidO3eqR4+mC+vZbDbZbJzeBjzZ3E92q8JRr4Gdo3TlGZ3NjgPAy5l65iY4OFjDhg3TihUrGre53W6tWLFCo0ePPmr/vn37asuWLdq4cWPj4/LLL9fYsWO1ceNGLjkBXmh3QYUWrs2RJN13aX8m7ANw2kw/95uZmalJkyZp+PDhGjlypObOnauqqipNnjxZkjRx4kR17txZWVlZCgkJ0cCBA5sc365dO0k6ajsA7/DI0u1yuQ1d2D+BCfsAtAjTy8348eNVVFSkWbNmKT8/X6mpqVq2bFnjIOOcnBxZrV41NAhAM32xq0if7SxSoNWimZf0MzsOAB9hMQzDMDtEW7Lb7YqOjlZ5ebmioqLMjgP4LZfb0CV//lI7Cyr0q7O6ada4/mZHAuDBTub3N6dEAJhi0boc7SyoUHRokH73s55mxwHgQyg3ANpceU2dnvp4lyRpenovtQsLNjkRAF9CuQHQ5p5ZsVulVU71jI/QL0d1NTsOAB9DuQHQpvYUVeq1VfslSfdf1l9BAbwNAWhZvKsAaFMPfbBN9W5DP+sbr/N6x5kdB4APotwAaDOf7SzUZzuLFBRg0b2Xcus3gNZBuQHQJupcbj30wTZJ0k1npqh7XITJiQD4KsoNgDbx99UHtKeoSh3CgzXtZ73MjgPAh1FuALS64kqHnv6k4dbv32f0UVRIkMmJAPgyyg2AVvf4sh2qqG1Y9fu64SxwC6B1UW4AtKpvc47ojW8OSpIeuHygAlj1G0Aro9wAaDVut6HZ//5OknT10CQN69re5EQA/AHlBkCreeObXG0+WK5IW6DuvriP2XEA+AnKDYBWUVbt1GPLdkiSpl/QW/GRISYnAuAvKDcAWsWc5bt0pLpOvRMiNHE060cBaDuUGwAt7rtD5Xr96wOSpD9ePoD1owC0Kd5xALQot9vQfe9ulduQLh3cUWf2iDU7EgA/Q7kB0KIWrcvVtzllCg8O0P2X9jc7DgA/RLkB0GKKKx2Ng4gzL+yjxGgGEQNoe5QbAC0ma+kOldfUqX/HKE1iEDEAk1BuALSINXtL9PaGg7JYpIevHKhABhEDMAnvPgBOm7Perfve3SpJmjCyi87owkzEAMxDuQFw2l76ap92F1aqQ3iw7s7oa3YcAH6OcgPgtOSUVOvPK3ZJku65pJ+iw4JMTgTA31FuAJwywzB077tbVFvn1ujuHXTV0M5mRwIAyg2AU/fOt3n6cnexbIFWZV01SBaLxexIAEC5AXBqSiod+tMH2yRJd6T3UkpsuMmJAKAB5QbAKfnTB9t0pLpO/TpGaco53c2OAwCNKDcATtrnu4r07sZDslqkR68axMKYADwK70gATkq1s173vrNFknTTmd00JLmduYEA4H9QbgCclCc/2qWDR2rUuV2o/u/C3mbHAYCjUG4ANNu6/aV6ZdU+SQ1LLITbAk1OBABHo9wAaJYap0t/eHOTDEO6bniSxvSJNzsSABwT5QZAszz58U7tL6lWYlSI7r20v9lxAOC4KDcAftI3+0v18n8aLkdlXT1I0aEssQDAc1FuAJxQbZ1Lf3hrswxDunZYksZyOQqAh6PcADihpz7eqX3FVUqIsum+y7gcBcDzUW4AHNeavSV68av/Xo66istRALwD5QbAMVXU1inzjU2Nl6PO75tgdiQAaBbKDYBjeuD9bcorq1FyTKhmXz7A7DgA0GyUGwBHWbb1sN5af1BWizTnulRFMFkfAC9CuQHQRGFFrWb+q2HtqFvP66ERKTEmJwKAk0O5AdDIMAzd9dZmHamuU/+OUZqeztpRALwP5QZAo9fX5GjlziIFB1o19xepCg7kLQKA9+GdC4AkaUe+XQ99sE2SNOOivuqdEGlyIgA4NZQbAKpxunT7wm/lqHdrbJ84TT4rxexIAHDKKDcA9OAH3ym7sFLxkTY9ee0QWSwWsyMBwCmj3AB+7oPNh/TPtbmyWKS541PVIcJmdiQAOC2UG8CP5ZZWa+bbDbd9Tx3TU2f2jDU5EQCcPsoN4Kec9W79btG3qnDUa1jX9pqe3svsSADQIig3gJ/K+nC7vs0pU1RIoP78i1QFBvB2AMA38G4G+KEPNh/SK//ZL0l66rpUJbUPMzcQALQgyg3gZ7ILK3T3W5slSbeN6aEL+rPaNwDfQrkB/EiVo163vr5BVU6XRnfvoP+7gOUVAPgeyg3gJwzD0Ix/bWmcz+aZCWcwzgaAT+KdDfATr67ar/c3HVKg1aK/3DBUcZHMZwPAN1FuAD+wKrtYDy3ZLkmacXFfDU+JMTkRALQejyg38+bNU0pKikJCQpSWlqa1a9ced98FCxbonHPOUfv27dW+fXulp6efcH/A3+WUVOu3CzfI5TZ05RmddfPZ3cyOBACtyvRys3jxYmVmZmr27NnasGGDhgwZooyMDBUWFh5z/5UrV2rChAn67LPPtHr1aiUnJ+vCCy9UXl5eGycHPF+lo163/G2dyqrrNCQpWllXDWLdKAA+z2IYhmFmgLS0NI0YMULPPfecJMntdis5OVnTpk3TjBkzfvJ4l8ul9u3b67nnntPEiRN/cn+73a7o6GiVl5crKirqtPMDnsrtNvSb19dr+bYCxUfa9O/bz1ZidIjZsQDglJzM729Tz9w4nU6tX79e6enpjdusVqvS09O1evXqZv0Z1dXVqqurU0zMsccQOBwO2e32Jg/AH8z9ZJeWbytQcIBV828cRrEB4DdMLTfFxcVyuVxKSGg6iVhCQoLy8/Ob9Wfcfffd6tSpU5OC9GNZWVmKjo5ufCQnJ592bsDTvbcxT898mi1JeuSqQRrapb3JiQCg7Zg+5uZ0PProo1q0aJHeeecdhYQc+3+lM2fOVHl5eeMjNze3jVMCbWvN3hL94c2GGYh/fW53XTMsyeREANC2As384bGxsQoICFBBQUGT7QUFBUpMTDzhsU8++aQeffRRffLJJxo8ePBx97PZbLLZmM8D/mFPUaV+/ff1crrcunhgomZc1NfsSADQ5kw9cxMcHKxhw4ZpxYoVjdvcbrdWrFih0aNHH/e4xx9/XH/605+0bNkyDR8+vC2iAh6vuNKhm15Zq/KaOp3RpZ2eHp8qq5U7owD4H1PP3EhSZmamJk2apOHDh2vkyJGaO3euqqqqNHnyZEnSxIkT1blzZ2VlZUmSHnvsMc2aNUsLFy5USkpK49iciIgIRUREmPY8ADPVOF265bVvlFtaoy4xYXpx4nCFBAWYHQsATGF6uRk/fryKioo0a9Ys5efnKzU1VcuWLWscZJyTkyOr9YcTTM8//7ycTqeuueaaJn/O7Nmz9cc//rEtowMeod7l1h2LvtXG3DK1CwvSq5NHqEMEl2IB+C/T57lpa8xzA19iGIbuemuz3lx/UMGBVr1+c5pGdmNpBQC+x2vmuQFw6gzD0CNLt+vN9QdltUjPTjiDYgMAotwAXuv5z/dowZf7JEmPXT1YGQNOfIchAPgLyg3ghRauydHjy3ZKku67tJ+uHc7klADwPcoN4GXe25ine9/dIkmaOraHbjmnu8mJAMCzUG4AL/L+pkO6c/FGGYZ0fVoX/f7CPmZHAgCPQ7kBvMSSzYc1ffFGuQ1p/PBkPfTzgbJYmKQPAP4X5QbwAh9uOazfLfpWLreha4YlKeuqQcw+DADHQbkBPNyyrfma9s+GYnPVGZ312NWDKTYAcAKmz1AM4Pje25inzDc2yeU2dEVqJz1x7RAFUGwA4IQoN4CHWrgmR/e+u0WGIV11Rmc9fs1gig0ANAPlBvBAL3yxR48s3SFJunFUVz1w+QAuRQFAM1FuAA9iGIaeXr5Lz3yaLUm6bUwP3ZXRh7uiAOAkUG4AD1HvcuuP73+n17/OkSTddVEf/XZMT5NTAYD3odwAHqDaWa9pC7/Vih2FslikBy8foBtHp5gdCwC8EuUGMFlRhUO3vLZOmw6WyxZo1Z9/cYYuGsgimABwqig3gIn2FFXqplfWKre0Ru3DgvTipBEa1rW92bEAwKtRbgCTfLW7WFMXblB5TZ26xITp1ckj1D0uwuxYAOD1KDdAGzMMQ6/8Z78eXrpdLreh1OR2enHScMVG2MyOBgA+gXIDtCFHvUv3vbNVb64/KEm6amhnPXLlIIUEBZicDAB8B+UGaCOF9lrd+vp6bcgpk9Ui3XNJP918djfmsAGAFka5AdrAquxi/W7RtyqudCoqJFDPXT9U5/aOMzsWAPgkyg3QilxuQ899mq25K3bJMKS+iZF6/pfD1C023OxoAOCzKDdAKymudGj6oo36KrtYkjR+eLL+ePkAhQYzvgYAWhPlBmgFX+wq0u/f3KTCCodCgwL00BUDdfWwJLNjAYBfoNwALai2zqVHP9yhV1ftlyT1jI/QX24Yqt4JkeYGAwA/QrkBWsjWvHJNX7xR2YWVkqSJo7tq5sX9uAwFAG2McgOcJme9W3/9fI+e+XS36lyG4iJteuKawRrTJ97saADglyg3wGnYlFumu9/erB35FZKkiwYk6pGrBikmPNjkZADgvyg3wCmodtZrzse79PJ/9sltSDHhwZp1WX/9PLUTk/IBgMkoN8BJ+nRHgWb/+zvlltZIkn6e2kmzLuuvDqwNBQAegXIDNNOBkio9+P42rdhRKEnqFB2ih64cqPP7JpicDADwY5Qb4CfUOF36y8ps/fWLvXLWuxVotejmc7pp2vm9FGHjnxAAeBremYHjcLkNvb3hoOZ8vEv59lpJ0jm9YjV73AD1jI8wOR0A4HgoN8D/MAxDK3cV6dGlO7SzoOEuqM7tQnX/Zf2UMSCRAcMA4OEoN8CPrD9wRE99vFOr9pRIkqJCAjXt/F66cXRXhQQxGR8AeAPKDSDp25wjevqT3fpiV5EkKTjAqpvOStFvx/RQuzDmrAEAb0K5gV/bkHNEz6zYrZU7G0pNgNWia4Ym6fbzeyo5JszkdACAU0G5gd9xuw2t3FWo+Z/v1dp9pZIaSs3VQzvr9rG91KUDpQYAvBnlBn7DUe/S+5sO64Uv9mhXQcPilkEBFl2R2lm3n99TXTuEm5wQANASKDfweYfLa7RwTY7+uTZHxZVOSVKELVA3pHXR5LO6KTE6xOSEAICWRLmBTzIMQ1/vLdXfv96vj74rkMttSJISo0I06cwU3TCqi6JCgkxOCQBoDZQb+JT88lq9veGg3vgmVwdKqhu3p3WL0aQzU3RB/wQFBVhNTAgAaG2UG3i92jqXPtleoH9tyNPKnYX670kaRdgCdXlqJ00c3VV9E6PMDQkAaDOUG3ilepdbq/aU6N2Nefr4uwJVOuobvzcipb2uG56sSwd3VFgwf8UBwN/wzg+v4ax3a/XeEi3bmq/l2/IbBwdLDcsj/Dy1k64elqQecaz7BAD+jHIDj1ZRW6evdhdr+bYCfbK9QPbaH87QtA8L0mWDO+nnqZ00rGt71nwCAEii3MDDGIahfcVV+mxnkT7dUaC1+0pV5zIavx8bEawLByTq4oGJGtW9A4ODAQBHodzAdKVVTv0nu1hf7S7WV9nFyiurafL9brHhGtsnXhcNTNSwru0VYOUMDQDg+Cg3aHOlVU6t3Veir/eWas2+Uu3It8v44eSMggIsGpESo/P7xuv8vvHqzhgaAMBJoNygVbndhvYWV2rDgTJtyDmi9QeOaHdh5VH79UmI1Nm9YnV2r1ildYvhLicAwCnjNwhajGEYyrfXavPBcm05WK7NeeXalFum8pq6o/btnRChtG4dlNY9RiO7xSg+kiUQAAAtg3KDU1LvcmtvcZW2HbJr+2G7th1u+Pjj27O/FxJk1eCkdhrapb2Gdmmn4SkxigkPNiE1AMAfUG5wQs56t3JKq5RdWKXdBRXaVVipXfkV2ltc2eQupu8FWC3qkxCpwUnRGpQUrcGd26lvx0juagIAtBnKDVTncivvSI0OlFbrQEmVDpRUa19xlfYWVSr3SE3jopP/Kzw4QH07Rql/xyj16xilfh0j1a9jlEKCAtr4GQAA8APKjR+orXMpv7xWh8prdLisVgeP1OjgkWrlHqnWwSM1Olxee9wCIzWUmO5xEeqVEKHeCZHqnRChXvGR6twuVFZuywYAeBjKjRerc7lVUulUUYVDRZW1KqpwKL/coYKKWhWU1yrfXqsCe+0xx8H8L1ugVV07hKlLTLhSOoQpJTZc3ePC1SMuQvGRNmb/BQB4DcqNhzAMQ5WOepVV16m8pk5Hqp0qrXLqSJVTpdV1DR+rnCqudKikyqmSSoeOVB99F9LxhAYFqGO7EHWKDlVS++8fYY0f4yNtnIUBAPgEyk0L+f4sSqWjTpUOl6oc9aqorVelo16VtXWqqK1XxX+32WvrZK+pk722vuFjTZ3KaupOeGnoeAKtFsVG2BQXaVNsRLASo0MUHxmixOgQJUTZlBDVUGjahQVx9gUA4Bc8otzMmzdPTzzxhPLz8zVkyBA9++yzGjly5HH3f/PNN3X//fdr//796tWrlx577DFdcsklbZj4aOsPHNEvXvj6tP+c4ECr2ocFqV1osGLCGx7tw4PUPixYHcKD1SHCpg4RwYqNsKlDeLDahwVzxgUAgB8xvdwsXrxYmZmZmj9/vtLS0jR37lxlZGRo586dio+PP2r/VatWacKECcrKytJll12mhQsX6oorrtCGDRs0cOBAE55BgwhboAKtFkWEBCo8OFARtsCGz22BigwJVFRIw7bIkCBFhQQqKjRIUSFBDR9DA9UuNFjtwoK40wgAgNNkMQzj5K+FtKC0tDSNGDFCzz33nCTJ7XYrOTlZ06ZN04wZM47af/z48aqqqtIHH3zQuG3UqFFKTU3V/Pnzj9rf4XDI4XA0fm2325WcnKzy8nJFRUW12PP4/mXk0g8AAC3PbrcrOjq6Wb+/TZ1Zzel0av369UpPT2/cZrValZ6ertWrVx/zmNWrVzfZX5IyMjKOu39WVpaio6MbH8nJyS33BH7EYrFQbAAA8ACmlpvi4mK5XC4lJCQ02Z6QkKD8/PxjHpOfn39S+8+cOVPl5eWNj9zc3JYJDwAAPJLpY25am81mk81mMzsGAABoI6aeuYmNjVVAQIAKCgqabC8oKFBiYuIxj0lMTDyp/QEAgH8xtdwEBwdr2LBhWrFiReM2t9utFStWaPTo0cc8ZvTo0U32l6Tly5cfd38AAOBfTL8slZmZqUmTJmn48OEaOXKk5s6dq6qqKk2ePFmSNHHiRHXu3FlZWVmSpDvuuEPnnXeennrqKV166aVatGiRvvnmG73wwgtmPg0AAOAhTC8348ePV1FRkWbNmqX8/HylpqZq2bJljYOGc3JyZLX+cILpzDPP1MKFC3XffffpnnvuUa9evfTuu++aOscNAADwHKbPc9PWTuY+eQAA4Bm8Zp4bAACAlka5AQAAPoVyAwAAfArlBgAA+BTKDQAA8CmUGwAA4FNMn+emrX1/57vdbjc5CQAAaK7vf283ZwYbvys3FRUVkqTk5GSTkwAAgJNVUVGh6OjoE+7jd5P4ud1uHTp0SJGRkbJYLGbHMZ3dbldycrJyc3OZ1LCV8Vq3LV7vtsNr3Xb8+bU2DEMVFRXq1KlTk5ULjsXvztxYrVYlJSWZHcPjREVF+d0/FLPwWrctXu+2w2vddvz1tf6pMzbfY0AxAADwKZQbAADgUyg3fs5ms2n27Nmy2WxmR/F5vNZti9e77fBatx1e6+bxuwHFAADAt3HmBgAA+BTKDQAA8CmUGwAA4FMoNwAAwKdQbnBMDodDqampslgs2rhxo9lxfM7+/ft18803q1u3bgoNDVWPHj00e/ZsOZ1Os6P5hHnz5iklJUUhISFKS0vT2rVrzY7kc7KysjRixAhFRkYqPj5eV1xxhXbu3Gl2LL/w6KOPymKxaPr06WZH8ViUGxzTXXfdpU6dOpkdw2ft2LFDbrdbf/3rX/Xdd9/p6aef1vz583XPPfeYHc3rLV68WJmZmZo9e7Y2bNigIUOGKCMjQ4WFhWZH8ymff/65pk6dqq+//lrLly9XXV2dLrzwQlVVVZkdzaetW7dOf/3rXzV48GCzo3g0bgXHUT788ENlZmbq7bff1oABA/Ttt98qNTXV7Fg+74knntDzzz+vvXv3mh3Fq6WlpWnEiBF67rnnJDWsJ5ecnKxp06ZpxowZJqfzXUVFRYqPj9fnn3+uc8891+w4PqmyslJDhw7VX/7yFz300ENKTU3V3LlzzY7lkThzgyYKCgo0ZcoU/f3vf1dYWJjZcfxKeXm5YmJizI7h1ZxOp9avX6/09PTGbVarVenp6Vq9erWJyXxfeXm5JPF3uBVNnTpVl156aZO/3zg2v1s4E8dnGIZuuukm3XrrrRo+fLj2799vdiS/kZ2drWeffVZPPvmk2VG8WnFxsVwulxISEppsT0hI0I4dO0xK5fvcbremT5+us846SwMHDjQ7jk9atGiRNmzYoHXr1pkdxStw5sYPzJgxQxaL5YSPHTt26Nlnn1VFRYVmzpxpdmSv1dzX+sfy8vJ00UUX6dprr9WUKVNMSg6cuqlTp2rr1q1atGiR2VF8Um5uru644w794x//UEhIiNlxvAJjbvxAUVGRSkpKTrhP9+7ddd111+n999+XxWJp3O5yuRQQEKAbbrhBr732WmtH9XrNfa2Dg4MlSYcOHdKYMWM0atQovfrqq7Ja+f/G6XA6nQoLC9Nbb72lK664onH7pEmTVFZWpvfee8+8cD7q9ttv13vvvacvvvhC3bp1MzuOT3r33Xd15ZVXKiAgoHGby+WSxWKR1WqVw+Fo8j1QbvAjOTk5stvtjV8fOnRIGRkZeuutt5SWlqakpCQT0/mevLw8jR07VsOGDdPrr7/Om1MLSUtL08iRI/Xss89Karhk0qVLF91+++0MKG5BhmFo2rRpeuedd7Ry5Ur16tXL7Eg+q6KiQgcOHGiybfLkyerbt6/uvvtuLgUeA2Nu0KhLly5Nvo6IiJAk9ejRg2LTwvLy8jRmzBh17dpVTz75pIqKihq/l5iYaGIy75eZmalJkyZp+PDhGjlypObOnauqqipNnjzZ7Gg+ZerUqVq4cKHee+89RUZGKj8/X5IUHR2t0NBQk9P5lsjIyKMKTHh4uDp06ECxOQ7KDWCC5cuXKzs7W9nZ2UcVR06mnp7x48erqKhIs2bNUn5+vlJTU7Vs2bKjBhnj9Dz//POSpDFjxjTZ/sorr+imm25q+0DAj3BZCgAA+BRGLwIAAJ9CuQEAAD6FcgMAAHwK5QYAAPgUyg0AAPAplBsAAOBTKDcAAMCnUG4AAIBPodwAAACfQrkBAAA+hXIDAAB8CuUGgNfbv3+/LBbLUY//XdQRgH9gVXAAXi85OVmHDx9u/Do/P1/p6ek699xzTUwFwCysCg7Ap9TW1mrMmDGKi4vTe++9J6uVE9SAv+HMDQCf8qtf/UoVFRVavnw5xQbwU5QbAD7joYce0kcffaS1a9cqMjLS7DgATMJlKQA+4e2339aECRP04Ycf6mc/+5nZcQCYiHIDwOtt3bpVaWlpyszM1NSpUxu3BwcHKyYmxsRkAMxAuQHg9V599VVNnjz5qO3nnXeeVq5c2faBAJiKcgMAAHwKtxIAAACfQrkBAAA+hXIDAAB8CuUGAAD4FMoNAADwKZQbAADgUyg3AADAp1BuAACAT6HcAAAAn0K5AQAAPoVyAwAAfMr/A/JLZgEUSZtJAAAAAElFTkSuQmCC\n"
          },
          "metadata": {}
        }
      ],
      "source": [
        "# 로지스틱 함수 예시\n",
        "\n",
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "from scipy.special import expit\n",
        "\n",
        "# -5와 5 사이 0.1 간격으로 z 값 배열 생성\n",
        "z = np.arange(-5, 5, 0.1)\n",
        "\n",
        "# 각 z마다 phi 계산\n",
        "phi = 1 / (1 + np.exp(-z))\n",
        "\n",
        "# scipy 라이브러리로도 시그모이드 함수 계산 가능\n",
        "# phi = expit(z)\n",
        "\n",
        "# 그래프 그리기\n",
        "plt.plot(z, phi)\n",
        "plt.xlabel('z')\n",
        "plt.ylabel('phi')\n",
        "plt.show()"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "1J6LGKpUJbFE"
      },
      "source": [
        "### 로지스틱 회귀로 이진 분류 수행하기"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "Khxh-3t5-2Tk"
      },
      "outputs": [],
      "source": [
        "# train_target == 'Bream'  ; target 값이 Bream인 것만 True, 나머지는 False인 배열 반환\n",
        "# (train_target == 'Bream') | (train_target == 'Smelt')  ; Bream과 Smelt 행만 True인 배열 반환\n",
        "bream_smelt_indexes = (train_target == 'Bream') | (train_target == 'Smelt')\n",
        "\n",
        "# 넘파이 boolean 인덱싱으로 원하는 행 추출\n",
        "train_bream_smelt = train_scaled[bream_smelt_indexes]\n",
        "target_bream_smelt = train_target[bream_smelt_indexes]"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 74
        },
        "id": "jEzP0aeXANra",
        "outputId": "9bf8491d-319b-4fa2-9438-8649f81c6499"
      },
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "LogisticRegression()"
            ],
            "text/html": [
              "<style>#sk-container-id-1 {color: black;background-color: white;}#sk-container-id-1 pre{padding: 0;}#sk-container-id-1 div.sk-toggleable {background-color: white;}#sk-container-id-1 label.sk-toggleable__label {cursor: pointer;display: block;width: 100%;margin-bottom: 0;padding: 0.3em;box-sizing: border-box;text-align: center;}#sk-container-id-1 label.sk-toggleable__label-arrow:before {content: \"▸\";float: left;margin-right: 0.25em;color: #696969;}#sk-container-id-1 label.sk-toggleable__label-arrow:hover:before {color: black;}#sk-container-id-1 div.sk-estimator:hover label.sk-toggleable__label-arrow:before {color: black;}#sk-container-id-1 div.sk-toggleable__content {max-height: 0;max-width: 0;overflow: hidden;text-align: left;background-color: #f0f8ff;}#sk-container-id-1 div.sk-toggleable__content pre {margin: 0.2em;color: black;border-radius: 0.25em;background-color: #f0f8ff;}#sk-container-id-1 input.sk-toggleable__control:checked~div.sk-toggleable__content {max-height: 200px;max-width: 100%;overflow: auto;}#sk-container-id-1 input.sk-toggleable__control:checked~label.sk-toggleable__label-arrow:before {content: \"▾\";}#sk-container-id-1 div.sk-estimator input.sk-toggleable__control:checked~label.sk-toggleable__label {background-color: #d4ebff;}#sk-container-id-1 div.sk-label input.sk-toggleable__control:checked~label.sk-toggleable__label {background-color: #d4ebff;}#sk-container-id-1 input.sk-hidden--visually {border: 0;clip: rect(1px 1px 1px 1px);clip: rect(1px, 1px, 1px, 1px);height: 1px;margin: -1px;overflow: hidden;padding: 0;position: absolute;width: 1px;}#sk-container-id-1 div.sk-estimator {font-family: monospace;background-color: #f0f8ff;border: 1px dotted black;border-radius: 0.25em;box-sizing: border-box;margin-bottom: 0.5em;}#sk-container-id-1 div.sk-estimator:hover {background-color: #d4ebff;}#sk-container-id-1 div.sk-parallel-item::after {content: \"\";width: 100%;border-bottom: 1px solid gray;flex-grow: 1;}#sk-container-id-1 div.sk-label:hover label.sk-toggleable__label {background-color: #d4ebff;}#sk-container-id-1 div.sk-serial::before {content: \"\";position: absolute;border-left: 1px solid gray;box-sizing: border-box;top: 0;bottom: 0;left: 50%;z-index: 0;}#sk-container-id-1 div.sk-serial {display: flex;flex-direction: column;align-items: center;background-color: white;padding-right: 0.2em;padding-left: 0.2em;position: relative;}#sk-container-id-1 div.sk-item {position: relative;z-index: 1;}#sk-container-id-1 div.sk-parallel {display: flex;align-items: stretch;justify-content: center;background-color: white;position: relative;}#sk-container-id-1 div.sk-item::before, #sk-container-id-1 div.sk-parallel-item::before {content: \"\";position: absolute;border-left: 1px solid gray;box-sizing: border-box;top: 0;bottom: 0;left: 50%;z-index: -1;}#sk-container-id-1 div.sk-parallel-item {display: flex;flex-direction: column;z-index: 1;position: relative;background-color: white;}#sk-container-id-1 div.sk-parallel-item:first-child::after {align-self: flex-end;width: 50%;}#sk-container-id-1 div.sk-parallel-item:last-child::after {align-self: flex-start;width: 50%;}#sk-container-id-1 div.sk-parallel-item:only-child::after {width: 0;}#sk-container-id-1 div.sk-dashed-wrapped {border: 1px dashed gray;margin: 0 0.4em 0.5em 0.4em;box-sizing: border-box;padding-bottom: 0.4em;background-color: white;}#sk-container-id-1 div.sk-label label {font-family: monospace;font-weight: bold;display: inline-block;line-height: 1.2em;}#sk-container-id-1 div.sk-label-container {text-align: center;}#sk-container-id-1 div.sk-container {/* jupyter's `normalize.less` sets `[hidden] { display: none; }` but bootstrap.min.css set `[hidden] { display: none !important; }` so we also need the `!important` here to be able to override the default hidden behavior on the sphinx rendered scikit-learn.org. See: https://github.com/scikit-learn/scikit-learn/issues/21755 */display: inline-block !important;position: relative;}#sk-container-id-1 div.sk-text-repr-fallback {display: none;}</style><div id=\"sk-container-id-1\" class=\"sk-top-container\"><div class=\"sk-text-repr-fallback\"><pre>LogisticRegression()</pre><b>In a Jupyter environment, please rerun this cell to show the HTML representation or trust the notebook. <br />On GitHub, the HTML representation is unable to render, please try loading this page with nbviewer.org.</b></div><div class=\"sk-container\" hidden><div class=\"sk-item\"><div class=\"sk-estimator sk-toggleable\"><input class=\"sk-toggleable__control sk-hidden--visually\" id=\"sk-estimator-id-1\" type=\"checkbox\" checked><label for=\"sk-estimator-id-1\" class=\"sk-toggleable__label sk-toggleable__label-arrow\">LogisticRegression</label><div class=\"sk-toggleable__content\"><pre>LogisticRegression()</pre></div></div></div></div></div>"
            ]
          },
          "metadata": {},
          "execution_count": 14
        }
      ],
      "source": [
        "from sklearn.linear_model import LogisticRegression\n",
        "\n",
        "lr = LogisticRegression()\n",
        "# train\n",
        "lr.fit(train_bream_smelt, target_bream_smelt)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "Mm60bpr7EQKU",
        "outputId": "e1d292e1-3fc0-4b4a-ba40-d0d0c80b83a2"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Bream' 'Smelt']\n"
          ]
        }
      ],
      "source": [
        "print(lr.classes_)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "3H_qieV-_CTt",
        "outputId": "475b55bb-992d-4c92-fa90-cd48344c551f"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[0.99759855 0.00240145]\n",
            " [0.02735183 0.97264817]\n",
            " [0.99486072 0.00513928]\n",
            " [0.98584202 0.01415798]\n",
            " [0.99767269 0.00232731]]\n"
          ]
        }
      ],
      "source": [
        "# 처음 5개 샘플의 클래스별 확률 출력\n",
        "print(lr.predict_proba(train_bream_smelt[:5]))"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "print(lr.predict(train_bream_smelt[:5]))"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "rxKmlZJR7P8_",
        "outputId": "5fe60215-cd1c-49e5-caa9-f656f47931fe"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Bream' 'Smelt' 'Bream' 'Bream' 'Bream']\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "1mvoYhUVQmFY",
        "outputId": "34c0c944-7b96-451c-e0fa-14c8814a5275"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[-0.4037798  -0.57620209 -0.66280298 -1.01290277 -0.73168947]] [-2.16155132]\n"
          ]
        }
      ],
      "source": [
        "# 계수 출력\n",
        "print(lr.coef_, lr.intercept_)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "SxrRy9m8A5Hy",
        "outputId": "7319dfa2-1cc9-4d35-931d-633550cf1321"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[-6.02927744  3.57123907 -5.26568906 -4.24321775 -6.0607117 ]\n"
          ]
        }
      ],
      "source": [
        "# z 값 출력\n",
        "decisions = lr.decision_function(train_bream_smelt[:5])\n",
        "print(decisions)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "SeuhSRuiA9yZ",
        "outputId": "a05c6327-fb35-4dac-bf6b-40252a6d9d04"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[0.00240145 0.97264817 0.00513928 0.01415798 0.00232731]\n"
          ]
        }
      ],
      "source": [
        "# phi 값 출력\n",
        "print(expit(decisions))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "6ee-s4l7EuVo"
      },
      "source": [
        "### 로지스틱 회귀로 다중 분류 수행하기"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "7QugsbD2X8bf",
        "outputId": "d825bac5-521a-4619-83d1-6608711b85e3"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "0.9327731092436975\n",
            "0.925\n"
          ]
        }
      ],
      "source": [
        "# max_iter로 반복 횟수 지정\n",
        "# C 매개변수 값을 크게 하여 규제 완화\n",
        "lr = LogisticRegression(C=20, max_iter=1000)\n",
        "lr.fit(train_scaled, train_target)\n",
        "\n",
        "print(lr.score(train_scaled, train_target))\n",
        "print(lr.score(test_scaled, test_target))"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "print(lr.classes_)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "KHN5OOhO_zuY",
        "outputId": "2dcafda8-b4c9-43c4-c28f-69906e7c3945"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Bream' 'Parkki' 'Perch' 'Pike' 'Roach' 'Smelt' 'Whitefish']\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "pqZosYezZOi3",
        "outputId": "a31583bb-b5f0-4bdf-c5bd-777e0a18f4c9"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[0.    0.014 0.841 0.    0.136 0.007 0.003]\n",
            " [0.    0.003 0.044 0.    0.007 0.946 0.   ]\n",
            " [0.    0.    0.034 0.935 0.015 0.016 0.   ]\n",
            " [0.011 0.034 0.306 0.007 0.567 0.    0.076]\n",
            " [0.    0.    0.904 0.002 0.089 0.002 0.001]]\n"
          ]
        }
      ],
      "source": [
        "# 처음 5개 샘플의 클래스별 확률 출력\n",
        "proba = lr.predict_proba(test_scaled[:5])\n",
        "# 소수점 3번째 자리까지 출력 (반올림)\n",
        "print(np.round(proba, decimals=3))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "0taO0XnF9dha",
        "outputId": "2d6fcc7d-216e-4f3e-9fe5-c48f29468801"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "['Perch' 'Smelt' 'Pike' 'Roach' 'Perch']\n"
          ]
        }
      ],
      "source": [
        "print(lr.predict(test_scaled[:5]))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "1swPv6ZOZTjg",
        "outputId": "8db0b77d-27e3-4661-e616-515e21ce9fb5"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "(7, 5) (7,)\n"
          ]
        }
      ],
      "source": [
        "print(lr.coef_.shape, lr.intercept_.shape)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "s9iRz1iAd7Oe",
        "outputId": "252c9024-f981-4ecd-fcfb-f9e1bf2d822d"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[ -6.5    1.03   5.16  -2.73   3.34   0.33  -0.63]\n",
            " [-10.86   1.93   4.77  -2.4    2.98   7.84  -4.26]\n",
            " [ -4.34  -6.23   3.17   6.49   2.36   2.42  -3.87]\n",
            " [ -0.68   0.45   2.65  -1.19   3.26  -5.75   1.26]\n",
            " [ -6.4   -1.99   5.82  -0.11   3.5   -0.11  -0.71]]\n"
          ]
        }
      ],
      "source": [
        "decision = lr.decision_function(test_scaled[:5])\n",
        "# 소수점 2번째 자리까지 출력 (반올림)\n",
        "print(np.round(decision, decimals=2))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "49CcsDHZeJma",
        "outputId": "3a823853-f328-4874-97fb-6c9559219c3a"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[0.    0.014 0.841 0.    0.136 0.007 0.003]\n",
            " [0.    0.003 0.044 0.    0.007 0.946 0.   ]\n",
            " [0.    0.    0.034 0.935 0.015 0.016 0.   ]\n",
            " [0.011 0.034 0.306 0.007 0.567 0.    0.076]\n",
            " [0.    0.    0.904 0.002 0.089 0.002 0.001]]\n"
          ]
        }
      ],
      "source": [
        "from scipy.special import softmax\n",
        "\n",
        "# axis=1은 소프트맥스를 계산할 축이 각 행인 것 = 각 샘플에 대해 계산\n",
        "proba = softmax(decision, axis=1)\n",
        "print(np.round(proba, decimals=3))"
      ]
    }
  ],
  "metadata": {
    "colab": {
      "provenance": []
    },
    "kernelspec": {
      "display_name": "default:Python",
      "language": "python",
      "name": "conda-env-default-py"
    },
    "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.9.10"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}