{ "cells": [ { "cell_type": "markdown", "metadata": { "id": "1k2vhsMVv0Pk" }, "source": [ "# Tutorial 4 Assignment - Logistic Regression\n", "\n", "We have provided you with a preprocessed dataset, the first cell will load and set everything up for you.\n", "The objectives for you to complete are as follows:\n", "1. Code up the commented functions on your own.\n", "2. Every step that you must code are explained as comments, use them as hints.\n", "\n", "The last cell has the code set up for training the model. We expect each one to have trained the model, and note down the best accuracy that they can achieve, and the conditions required to do the same." ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "nj3rHkttqucf", "outputId": "0ff8d346-148a-4dd9-df8c-4696d2a002c3" }, "outputs": [], "source": [ "# !wget -O rainDataset.csv \"https://docs.google.com/spreadsheets/d/1RNtDIvisrnOmjJxS7aPm-45NtOH3qd5-mgd2bHeSOGA/export?format=csv&gid=1727131321\"\n", "import pandas as pd\n", "from sklearn.preprocessing import StandardScaler\n", "from sklearn.model_selection import train_test_split\n", "import numpy as np\n", "import math\n", "import matplotlib.pyplot as plt\n", "df=pd.read_csv('rainDataset.csv')\n", "%matplotlib inline\n", "df.head()\n", "X = df.drop(['RainTomorrow'], axis=1)\n", "y = df['RainTomorrow']\n", "scaler = StandardScaler()\n", "X_scaled = scaler.fit_transform(X)\n", "X_train, X_test, y_train, y_test = train_test_split(X_scaled, y, test_size=0.2, random_state=42, stratify=y)\n" ] }, { "cell_type": "markdown", "metadata": { "id": "kiZ4TBA6wFea" }, "source": [ "# CODE BELOW" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "id": "2A-RJ6Rscyoa" }, "outputs": [], "source": [ "import numpy as np\n", "\n", "def sigmoid(z):\n", " \"\"\"\n", " Compute the sigmoid function.\n", "\n", " Parameters:\n", " z : numpy array\n", " Linear combination of weights and input features.\n", "\n", " Returns:\n", " numpy array\n", " Sigmoid of input z.\n", " \"\"\"\n", " return 1 / (1 + math.e**(-z))\n", "\n", "def initialize_weights(n_features):\n", " \"\"\"\n", " Initialize weights and bias to zero.\n", "\n", " Parameters:\n", " n_features : int\n", " Number of features in the dataset.\n", "\n", " Returns:\n", " tuple\n", " Initialized weights and bias.\n", " \"\"\"\n", " # initialize the weights and bias to zero (hint: make sure dimentions are correct)\n", " w = np.zeros(n_features)\n", " b = np.array([0])\n", " return w, b\n", "\n", "def compute_cost(y, y_pred):\n", " \"\"\"\n", " Compute the cost function for logistic regression.\n", "\n", " Parameters:\n", " y : numpy array\n", " Actual labels.\n", " y_pred : numpy array\n", " Predicted probabilities.\n", "\n", " Returns:\n", " float\n", " The cost value.\n", " \"\"\"\n", " # compute the cost\n", " # Cross Entropy Loss:\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", " return cost\n", "\n", "def compute_gradients(X, y, y_pred):\n", " \"\"\"\n", " Compute the gradients for weights and bias.\n", "\n", " Parameters:\n", " X : numpy array\n", " Feature matrix.\n", " y : numpy array\n", " Actual labels.\n", " y_pred : numpy array\n", " Predicted probabilities.\n", "\n", " Returns:\n", " tuple\n", " Gradients of weights and bias.\n", " \"\"\"\n", " m = X.shape[0]\n", "\n", " # compute dw\n", " dw = np.mean(np.dot(X.T, (y-y_pred)))\n", " # compute db\n", " db = np.mean(y-y_pred)\n", " return dw, db\n", "\n", "\n", "def optimize(X, y, weights, bias, learning_rate, num_iterations):\n", " \"\"\"\n", " Perform gradient descent to optimize weights and bias.\n", "\n", " Parameters:\n", " X : numpy array\n", " Feature matrix.\n", " y : numpy array\n", " Actual labels.\n", " weights : numpy array\n", " Weights of the model.\n", " bias : float\n", " Bias of the model.\n", " learning_rate : float\n", " Learning rate for gradient descent.\n", " num_iterations : int\n", " Number of iterations for gradient descent.\n", "\n", " Returns:\n", " tuple\n", " Optimized weights, bias, and the list of costs.\n", " \"\"\"\n", " costs = []\n", "\n", " for i in range(num_iterations):\n", " # Compute linear model\n", " logits = computePred(X, weights, bias)\n", " # Apply sigmoid function\n", " yPred = sigmoid(logits)\n", " # Compute cost\n", " cost = compute_cost(y, yPred)\n", " costs.append(cost)\n", " # Compute gradients\n", " dw, db = compute_gradients(X, y, yPred)\n", " # Update weights and bias\n", " weights = weights - learning_rate * dw\n", " bias = bias - learning_rate * db\n", " return weights, bias, costs\n", "\n", "def predict(X, weights, bias):\n", " \"\"\"\n", " Predict the binary labels for a dataset.\n", "\n", " Parameters:\n", " X : numpy array\n", " Feature matrix.\n", " weights : numpy array\n", " Weights of the model.\n", " bias : float\n", " Bias of the model.\n", "\n", " Returns:\n", " numpy array\n", " Predicted binary labels (0 or 1).\n", " \"\"\"\n", " z = computePred(X, weights, bias)\n", " y_pred = sigmoid(z)\n", " predictions = [1 if i > 0.5 else 0 for i in y_pred]\n", " return np.array(predictions)\n", "\n", "def computePred(X, weights, bias):\n", " return np.dot(X, weights) + bias" ] }, { "cell_type": "markdown", "metadata": { "id": "v9rwH83Rwrfk" }, "source": [ "# COMPUTE ACCURACY" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "id": "eyGEV4mWB-rW", "scrolled": true }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Accuracy: 0.41052779633601744\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n", "/tmp/ipykernel_2559765/2227838026.py:50: RuntimeWarning: divide by zero encountered in log\n", " cost = -np.mean(y * np.log(y_pred) + (1-y) * np.log((1-y_pred)))\n", "/tmp/ipykernel_2559765/2227838026.py:15: RuntimeWarning: overflow encountered in power\n", " return 1 / (1 + math.e**(-z))\n" ] } ], "source": [ "weights, bias = initialize_weights(len(X_train[0]))\n", "weights, bias,costs = optimize(X_train, y_train, weights, bias, 0.01, 100)\n", "y_pred = predict(X_test,weights,bias)\n", "matches = np.sum(y_test == y_pred)\n", "mismatches = np.sum(y_test != y_pred)\n", "print(f\"Accuracy: {matches/(matches+mismatches)}\")" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGdCAYAAAAxCSikAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAok0lEQVR4nO3df3DU9Z3H8VeyJBsukk0lzeaHi6GtnBwicQLZRjpTHBejR0/peBgYaWgqZ8UcIvE8krEkrZWkFr2mlZwZmLRyej1SONoyQoMavI5oJBJOjxwkgVMJQjaQYna5VLN093v/HGtXksACIZ+E52PmOw6f7/vzyefzmei+/O73+yXGsixLAAAABosd6QkAAACcD4EFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGC8cSM9gcslFArp+PHjmjBhgmJiYkZ6OgAA4AJYlqXTp08rIyNDsbGDX0cZM4Hl+PHjcrlcIz0NAABwEY4eParrrrtu0PNjJrBMmDBB+v8FJyUljfR0AADABfD7/XK5XOHP8cFcVGCpqanR2rVr5fV6NWPGDD333HPKzc0dsHbOnDn6/e9/f077X//1X2v79u2SpO9///vatGmTjh49qvj4eOXk5GjNmjVyu90XPKezXwMlJSURWAAAGGXOdztH1Dfd1tfXq6SkRBUVFdq3b59mzJih/Px8nThxYsD6rVu3qqurK3y0trbKZrNpwYIF4ZopU6Zo3bp12r9/v3bv3q2srCzdcccdOnnyZLTTAwAAY1BMtH9bs9vt1qxZs7Ru3Trp/292dblcWr58uUpLS8/bv7q6WuXl5erq6lJiYuKANX6/Xw6HQ6+99ppuv/32C5rX2T4+n48rLAAAjBIX+vkd1RWWQCCglpYWeTyezwaIjZXH41FTU9MFjVFXV6eFCxcOGlYCgYDWr18vh8OhGTNmDDpOf3+//H5/xAEAAMamqAJLT0+PgsGgnE5nRLvT6ZTX6z1v/+bmZrW2tmrp0qXnnHv55Zd1zTXXKCEhQT/5yU/06quvKiUlZdCxqqqq5HA4wgdPCAEAMHZd0RfH1dXVafr06QPeoHvbbbfp3Xff1VtvvaU777xT991336D3xUhSWVmZfD5f+Dh69Ogwzx4AAIyUqAJLSkqKbDaburu7I9q7u7uVlpY2ZN++vj5t2rRJDzzwwIDnExMT9ZWvfEVf/epXVVdXp3Hjxqmurm7Q8ex2e/iJIJ4MAgBgbIsqsJx95LixsTHcFgqF1NjYqLy8vCH7bt68Wf39/Vq8ePEF/axQKKT+/v5opgcAAMaoqN/DUlJSoiVLlmjmzJnKzc1VdXW1+vr6VFRUJEkqLCxUZmamqqqqIvrV1dVp/vz5mjhxYkR7X1+f1qxZo7vvvlvp6enq6elRTU2Njh07FvHoMwAAuHpFHVgKCgp08uRJlZeXy+v1Kjs7Ww0NDeEbcTs7O8/5uwDa29u1e/duvfLKK+eMZ7PZ1NbWpo0bN6qnp0cTJ07UrFmz9MYbb2jatGmXsjYAADBGRP0eFlPxHhYAAEafYXkPCwAAwEggsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgvIsKLDU1NcrKylJCQoLcbream5sHrZ0zZ45iYmLOOebNmydJOnPmjFatWqXp06crMTFRGRkZKiws1PHjxy9+VQAAYEyJOrDU19erpKREFRUV2rdvn2bMmKH8/HydOHFiwPqtW7eqq6srfLS2tspms2nBggWSpD/+8Y/at2+fVq9erX379mnr1q1qb2/X3XfffemrAwAAY0KMZVlWNB3cbrdmzZqldevWSZJCoZBcLpeWL1+u0tLS8/avrq5WeXm5urq6lJiYOGDNO++8o9zcXB05ckSTJk26oHn5/X45HA75fD4lJSVFsyQAADBCLvTzO6orLIFAQC0tLfJ4PJ8NEBsrj8ejpqamCxqjrq5OCxcuHDSsSJLP51NMTIySk5MHrenv75ff7484AADA2BRVYOnp6VEwGJTT6Yxodzqd8nq95+3f3Nys1tZWLV26dNCaTz/9VKtWrdKiRYuGTFpVVVVyOBzhw+VyRbMUAAAwilzRp4Tq6uo0ffp05ebmDnj+zJkzuu+++2RZlp5//vkhxyorK5PP5wsfR48eHaZZAwCAkTYumuKUlBTZbDZ1d3dHtHd3dystLW3Ivn19fdq0aZOefPLJAc+fDStHjhzRrl27znsfit1ul91uj2b6AABglIrqCkt8fLxycnLU2NgYbguFQmpsbFReXt6QfTdv3qz+/n4tXrz4nHNnw8qhQ4f02muvaeLEidFMCwAAjHFRXWGRpJKSEi1ZskQzZ85Ubm6uqqur1dfXp6KiIklSYWGhMjMzVVVVFdGvrq5O8+fPPyeMnDlzRn/7t3+rffv26eWXX1YwGAzfD3PttdcqPj7+0lYIAABGvagDS0FBgU6ePKny8nJ5vV5lZ2eroaEhfCNuZ2enYmMjL9y0t7dr9+7deuWVV84Z79ixY9q2bZskKTs7O+Lc66+/rjlz5kQ7RQAAMMZE/R4WU/EeFgAARp9heQ8LAADASCCwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADDeRQWWmpoaZWVlKSEhQW63W83NzYPWzpkzRzExMecc8+bNC9ds3bpVd9xxhyZOnKiYmBi9++67F7caAAAwJkUdWOrr61VSUqKKigrt27dPM2bMUH5+vk6cODFg/datW9XV1RU+WltbZbPZtGDBgnBNX1+fvva1r+npp5++tNUAAIAxKcayLCuaDm63W7NmzdK6deskSaFQSC6XS8uXL1dpael5+1dXV6u8vFxdXV1KTEyMOPfhhx9q8uTJ+s///E9lZ2dHtRC/3y+HwyGfz6ekpKSo+gIAgJFxoZ/fUV1hCQQCamlpkcfj+WyA2Fh5PB41NTVd0Bh1dXVauHDhOWElWv39/fL7/REHAAAYm6IKLD09PQoGg3I6nRHtTqdTXq/3vP2bm5vV2tqqpUuXRj/Tz6mqqpLD4QgfLpfrkscEAABmuqJPCdXV1Wn69OnKzc295LHKysrk8/nCx9GjRy/LHAEAgHnGRVOckpIim82m7u7uiPbu7m6lpaUN2bevr0+bNm3Sk08+eXEz/Ry73S673X5ZxgIAAGaL6gpLfHy8cnJy1NjYGG4LhUJqbGxUXl7ekH03b96s/v5+LV68+OJnCwAArkpRXWGRpJKSEi1ZskQzZ85Ubm6uqqur1dfXp6KiIklSYWGhMjMzVVVVFdGvrq5O8+fP18SJE88Z89SpU+rs7NTx48clSe3t7ZKktLS08165AQAAY1/UgaWgoEAnT55UeXm5vF6vsrOz1dDQEL4Rt7OzU7GxkRdu2tvbtXv3br3yyisDjrlt27Zw4JGkhQsXSpIqKir0/e9/P9opAgCAMSbq97CYivewAAAw+gzLe1gAAABGAoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAONdVGCpqalRVlaWEhIS5Ha71dzcPGjtnDlzFBMTc84xb968cI1lWSovL1d6errGjx8vj8ejQ4cOXdyKAADAmBN1YKmvr1dJSYkqKiq0b98+zZgxQ/n5+Tpx4sSA9Vu3blVXV1f4aG1tlc1m04IFC8I1P/7xj/Wzn/1MtbW12rNnjxITE5Wfn69PP/300lYHAADGhBjLsqxoOrjdbs2aNUvr1q2TJIVCIblcLi1fvlylpaXn7V9dXa3y8nJ1dXUpMTFRlmUpIyNDjz32mP7hH/5BkuTz+eR0OvXCCy9o4cKFFzQvv98vh8Mhn8+npKSkaJYEAABGyIV+fkd1hSUQCKilpUUej+ezAWJj5fF41NTUdEFj1NXVaeHChUpMTJQkffDBB/J6vRFjOhwOud3uIcfs7++X3++POAAAwNgUVWDp6elRMBiU0+mMaHc6nfJ6veft39zcrNbWVi1dujTcdrZftGNWVVXJ4XCED5fLFc1SAADAKHJFnxKqq6vT9OnTlZube8ljlZWVyefzhY+jR49eljkCAADzRBVYUlJSZLPZ1N3dHdHe3d2ttLS0Ifv29fVp06ZNeuCBByLaz/aLdky73a6kpKSIAwAAjE1RBZb4+Hjl5OSosbEx3BYKhdTY2Ki8vLwh+27evFn9/f1avHhxRPvkyZOVlpYWMabf79eePXvOOyYAALg6jIu2Q0lJiZYsWaKZM2cqNzdX1dXV6uvrU1FRkSSpsLBQmZmZqqqqiuhXV1en+fPna+LEiRHtMTExevTRR/XUU0/phhtu0OTJk7V69WplZGRo/vz5l7o+AAAwBkQdWAoKCnTy5EmVl5fL6/UqOztbDQ0N4ZtmOzs7FRsbeeGmvb1du3fv1iuvvDLgmP/4j/+ovr4+Pfjgg+rt7dXXvvY1NTQ0KCEh4WLXBQAAxpCo38NiKt7DAgDA6DMs72EBAAAYCQQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMa7qMBSU1OjrKwsJSQkyO12q7m5ecj63t5eFRcXKz09XXa7XVOmTNGOHTvC50+fPq1HH31U119/vcaPH69bb71V77zzzsVMDQAAjEFRB5b6+nqVlJSooqJC+/bt04wZM5Sfn68TJ04MWB8IBDR37lx9+OGH2rJli9rb27VhwwZlZmaGa5YuXapXX31VL774ovbv36877rhDHo9Hx44du7TVAQCAMSHGsiwrmg5ut1uzZs3SunXrJEmhUEgul0vLly9XaWnpOfW1tbVau3at2traFBcXd875Tz75RBMmTNBvf/tbzZs3L9yek5Oju+66S0899dQFzcvv98vhcMjn8ykpKSmaJQEAgBFyoZ/fUV1hCQQCamlpkcfj+WyA2Fh5PB41NTUN2Gfbtm3Ky8tTcXGxnE6nbrrpJlVWVioYDEqS/vSnPykYDCohISGi3/jx47V79+5B59Lf3y+/3x9xAACAsSmqwNLT06NgMCin0xnR7nQ65fV6B+zz/vvva8uWLQoGg9qxY4dWr16tZ599NnzlZMKECcrLy9MPf/hDHT9+XMFgUC+99JKamprU1dU16FyqqqrkcDjCh8vlimYpAABgFBn2p4RCoZBSU1O1fv165eTkqKCgQE888YRqa2vDNS+++KIsy1JmZqbsdrt+9rOfadGiRYqNHXx6ZWVl8vl84ePo0aPDvRQAADBCxkVTnJKSIpvNpu7u7oj27u5upaWlDdgnPT1dcXFxstls4bapU6fK6/UqEAgoPj5eX/7yl/X73/9efX198vv9Sk9PV0FBgb70pS8NOhe73S673R7N9AEAwCgV1RWW+Ph45eTkqLGxMdwWCoXU2NiovLy8AfvMnj1bhw8fVigUCrd1dHQoPT1d8fHxEbWJiYlKT0/Xxx9/rJ07d+qee+6JfkUAAGDMiforoZKSEm3YsEEbN27UwYMHtWzZMvX19amoqEiSVFhYqLKysnD9smXLdOrUKa1YsUIdHR3avn27KisrVVxcHK7ZuXOnGhoa9MEHH+jVV1/VbbfdphtvvDE8JgAAuLpF9ZWQJBUUFOjkyZMqLy+X1+tVdna2GhoawjfidnZ2Rtx74nK5tHPnTq1cuVI333yzMjMztWLFCq1atSpc4/P5VFZWpo8++kjXXnut7r33Xq1Zs2bAx6ABAMDVJ+r3sJiK97AAADD6DMt7WAAAAEYCgQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA411UYKmpqVFWVpYSEhLkdrvV3Nw8ZH1vb6+Ki4uVnp4uu92uKVOmaMeOHeHzwWBQq1ev1uTJkzV+/Hh9+ctf1g9/+ENZlnUx0wMAAGPMuGg71NfXq6SkRLW1tXK73aqurlZ+fr7a29uVmpp6Tn0gENDcuXOVmpqqLVu2KDMzU0eOHFFycnK45umnn9bzzz+vjRs3atq0adq7d6+KiorkcDj0yCOPXPoqAQDAqBZjRXkZw+12a9asWVq3bp0kKRQKyeVyafny5SotLT2nvra2VmvXrlVbW5vi4uIGHPMb3/iGnE6n6urqwm333nuvxo8fr5deeumC5uX3++VwOOTz+ZSUlBTNkgAAwAi50M/vqL4SCgQCamlpkcfj+WyA2Fh5PB41NTUN2Gfbtm3Ky8tTcXGxnE6nbrrpJlVWVioYDIZrbr31VjU2Nqqjo0OS9N5772n37t266667Bp1Lf3+//H5/xAEAAMamqL4S6unpUTAYlNPpjGh3Op1qa2sbsM/777+vXbt26f7779eOHTt0+PBhPfzwwzpz5owqKiokSaWlpfL7/brxxhtls9kUDAa1Zs0a3X///YPOpaqqSj/4wQ+imT4AABilhv0poVAopNTUVK1fv145OTkqKCjQE088odra2nDNr371K/3rv/6rfvnLX2rfvn3auHGjnnnmGW3cuHHQccvKyuTz+cLH0aNHh3spAABghER1hSUlJUU2m03d3d0R7d3d3UpLSxuwT3p6uuLi4mSz2cJtU6dOldfrVSAQUHx8vB5//HGVlpZq4cKFkqTp06fryJEjqqqq0pIlSwYc1263y263RzN9AAAwSkV1hSU+Pl45OTlqbGwMt4VCITU2NiovL2/APrNnz9bhw4cVCoXCbR0dHUpPT1d8fLwk6Y9//KNiYyOnYrPZIvoAAICrV9RfCZWUlGjDhg3auHGjDh48qGXLlqmvr09FRUWSpMLCQpWVlYXrly1bplOnTmnFihXq6OjQ9u3bVVlZqeLi4nDN3/zN32jNmjXavn27PvzwQ/3617/WP/3TP+mb3/zm5VonAAAYxaJ+D0tBQYFOnjyp8vJyeb1eZWdnq6GhIXwjbmdnZ8TVEpfLpZ07d2rlypW6+eablZmZqRUrVmjVqlXhmueee06rV6/Www8/rBMnTigjI0Pf/e53VV5efrnWCQAARrGo38NiKt7DAgDA6DMs72EBAAAYCQQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMa7qMBSU1OjrKwsJSQkyO12q7m5ecj63t5eFRcXKz09XXa7XVOmTNGOHTvC57OyshQTE3POUVxcfDHTAwAAY8y4aDvU19erpKREtbW1crvdqq6uVn5+vtrb25WamnpOfSAQ0Ny5c5WamqotW7YoMzNTR44cUXJycrjmnXfeUTAYDP+5tbVVc+fO1YIFCy5lbQAAYIyIsSzLiqaD2+3WrFmztG7dOklSKBSSy+XS8uXLVVpaek59bW2t1q5dq7a2NsXFxV3Qz3j00Uf18ssv69ChQ4qJibmgPn6/Xw6HQz6fT0lJSdEsCQAAjJAL/fyO6iuhQCCglpYWeTyezwaIjZXH41FTU9OAfbZt26a8vDwVFxfL6XTqpptuUmVlZcQVlc//jJdeeknf+c53hgwr/f398vv9EQcAABibogosPT09CgaDcjqdEe1Op1Ner3fAPu+//762bNmiYDCoHTt2aPXq1Xr22Wf11FNPDVj/m9/8Rr29vfr2t7895FyqqqrkcDjCh8vlimYpAABgFBn2p4RCoZBSU1O1fv165eTkqKCgQE888YRqa2sHrK+rq9Ndd92ljIyMIcctKyuTz+cLH0ePHh2mFQAAgJEW1U23KSkpstls6u7ujmjv7u5WWlragH3S09MVFxcnm80Wbps6daq8Xq8CgYDi4+PD7UeOHNFrr72mrVu3nncudrtddrs9mukDAIBRKqorLPHx8crJyVFjY2O4LRQKqbGxUXl5eQP2mT17tg4fPqxQKBRu6+joUHp6ekRYkaRf/OIXSk1N1bx586JfCQAAGLOi/kqopKREGzZs0MaNG3Xw4EEtW7ZMfX19KioqkiQVFhaqrKwsXL9s2TKdOnVKK1asUEdHh7Zv367Kyspz3rESCoX0i1/8QkuWLNG4cVE/bQ0AAMawqJNBQUGBTp48qfLycnm9XmVnZ6uhoSF8I25nZ6diYz/LQS6XSzt37tTKlSt18803KzMzUytWrNCqVasixn3ttdfU2dmp73znO5djXQAAYAyJ+j0spuI9LAAAjD7D8h4WAACAkUBgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADAegQUAABiPwAIAAIxHYAEAAMYjsAAAAOMRWAAAgPEILAAAwHgEFgAAYDwCCwAAMB6BBQAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4FxVYampqlJWVpYSEBLndbjU3Nw9Z39vbq+LiYqWnp8tut2vKlCnasWNHRM2xY8e0ePFiTZw4UePHj9f06dO1d+/ei5keAAAYY8ZF26G+vl4lJSWqra2V2+1WdXW18vPz1d7ertTU1HPqA4GA5s6dq9TUVG3ZskWZmZk6cuSIkpOTwzUff/yxZs+erdtuu02/+93v9MUvflGHDh3SF77whUtfIQAAGPViLMuyoungdrs1a9YsrVu3TpIUCoXkcrm0fPlylZaWnlNfW1urtWvXqq2tTXFxcQOOWVpaqjfffFNvvPHGxa5Dfr9fDodDPp9PSUlJFz0OAAC4ci708zuqr4QCgYBaWlrk8Xg+GyA2Vh6PR01NTQP22bZtm/Ly8lRcXCyn06mbbrpJlZWVCgaDETUzZ87UggULlJqaqltuuUUbNmwYci79/f3y+/0RBwAAGJuiCiw9PT0KBoNyOp0R7U6nU16vd8A+77//vrZs2aJgMKgdO3Zo9erVevbZZ/XUU09F1Dz//PO64YYbtHPnTi1btkyPPPKINm7cOOhcqqqq5HA4wofL5YpmKQAAYBSJ+h6WaIVCIaWmpmr9+vWy2WzKycnRsWPHtHbtWlVUVIRrZs6cqcrKSknSLbfcotbWVtXW1mrJkiUDjltWVqaSkpLwn/1+P6EFAIAxKqrAkpKSIpvNpu7u7oj27u5upaWlDdgnPT1dcXFxstls4bapU6fK6/UqEAgoPj5e6enp+qu/+quIflOnTtW///u/DzoXu90uu90ezfQBAMAoFdVXQvHx8crJyVFjY2O4LRQKqbGxUXl5eQP2mT17tg4fPqxQKBRu6+joUHp6uuLj48M17e3tEf06Ojp0/fXXR7seAAAwBkX9HpaSkhJt2LBBGzdu1MGDB7Vs2TL19fWpqKhIklRYWKiysrJw/bJly3Tq1CmtWLFCHR0d2r59uyorK1VcXByuWblypd5++21VVlbq8OHD+uUvf6n169dH1AAAgKtX1PewFBQU6OTJkyovL5fX61V2drYaGhrCN+J2dnYqNvazHORyubRz506tXLlSN998szIzM7VixQqtWrUqXDNr1iz9+te/VllZmZ588klNnjxZ1dXVuv/++y/XOgEAwCgW9XtYTMV7WAAAGH2G5T0sAAAAI4HAAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADjEVgAAIDxCCwAAMB4BBYAAGA8AgsAADDeuJGewOViWZYkye/3j/RUAADABTr7uX32c3wwYyawnD59WpLkcrlGeioAACBKp0+flsPhGPR8jHW+SDNKhEIhHT9+XBMmTFBMTMxIT2dE+f1+uVwuHT16VElJSSM9nTGLfb5y2Osrg32+MtjnSJZl6fTp08rIyFBs7OB3qoyZKyyxsbG67rrrRnoaRklKSuJfhiuAfb5y2Osrg32+Mtjnzwx1ZeUsbroFAADGI7AAAADjEVjGILvdroqKCtnt9pGeypjGPl857PWVwT5fGezzxRkzN90CAICxiyssAADAeAQWAABgPAILAAAwHoEFAAAYj8AySp06dUr333+/kpKSlJycrAceeED/+7//O2SfTz/9VMXFxZo4caKuueYa3Xvvveru7h6w9g9/+IOuu+46xcTEqLe3d5hWYb7h2Of33ntPixYtksvl0vjx4zV16lT99Kc/vQKrMUdNTY2ysrKUkJAgt9ut5ubmIes3b96sG2+8UQkJCZo+fbp27NgRcd6yLJWXlys9PV3jx4+Xx+PRoUOHhnkV5ruc+3zmzBmtWrVK06dPV2JiojIyMlRYWKjjx49fgZWY73L/Tv+5hx56SDExMaqurh6GmY8iFkalO++805oxY4b19ttvW2+88Yb1la98xVq0aNGQfR566CHL5XJZjY2N1t69e62vfvWr1q233jpg7T333GPdddddliTr448/HqZVmG849rmurs565JFHrP/4j/+w/ud//sd68cUXrfHjx1vPPffcFVjRyNu0aZMVHx9v/fznP7f++7//2/q7v/s7Kzk52eru7h6w/s0337RsNpv14x//2Dpw4ID1ve99z4qLi7P2798frvnRj35kORwO6ze/+Y313nvvWXfffbc1efJk65NPPrmCKzPL5d7n3t5ey+PxWPX19VZbW5vV1NRk5ebmWjk5OVd4ZeYZjt/ps7Zu3WrNmDHDysjIsH7yk59cgdWYi8AyCh04cMCSZL3zzjvhtt/97ndWTEyMdezYsQH79Pb2WnFxcdbmzZvDbQcPHrQkWU1NTRG1//zP/2x9/etftxobG6/qwDLc+/znHn74Yeu22267zCswU25urlVcXBz+czAYtDIyMqyqqqoB6++77z5r3rx5EW1ut9v67ne/a1mWZYVCISstLc1au3Zt+Hxvb69lt9utf/u3fxu2dZjucu/zQJqbmy1J1pEjRy7jzEef4drrjz76yMrMzLRaW1ut66+//qoPLHwlNAo1NTUpOTlZM2fODLd5PB7FxsZqz549A/ZpaWnRmTNn5PF4wm033nijJk2apKampnDbgQMH9OSTT+pf/uVfhvxLqK4Gw7nPn+fz+XTttdde5hWYJxAIqKWlJWJ/YmNj5fF4Bt2fpqamiHpJys/PD9d/8MEH8nq9ETUOh0Nut3vIPR/LhmOfB+Lz+RQTE6Pk5OTLOPvRZbj2OhQK6Vvf+pYef/xxTZs2bRhXMHpc3Z9Io5TX61VqampE27hx43TttdfK6/UO2ic+Pv6c/7A4nc5wn/7+fi1atEhr167VpEmThnEFo8Nw7fPnvfXWW6qvr9eDDz54GWdvpp6eHgWDQTmdzoj2ofbH6/UOWX/2n9GMOdYNxz5/3qeffqpVq1Zp0aJFV/Vf4Ddce/30009r3LhxeuSRR4Zp5qMPgcUgpaWliomJGfJoa2sbtp9fVlamqVOnavHixcP2M0ww0vv851pbW3XPPfeooqJCd9xxxxX5mcClOnPmjO677z5ZlqXnn39+pKcz5rS0tOinP/2pXnjhBcXExIz0dIwxbqQngM889thj+va3vz1kzZe+9CWlpaXpxIkTEe1/+tOfdOrUKaWlpQ3YLy0tTYFAQL29vRH/99/d3R3us2vXLu3fv19btmyR/v/JC0lKSUnRE088oR/84AeXvEYTjPQ+n3XgwAHdfvvtevDBB/W9733vktY0WqSkpMhms53zdNpA+3NWWlrakPVn/9nd3a309PSImuzs7GFYhfmGY5/POhtWjhw5ol27dl3VV1c0THv9xhtv6MSJExFXuoPBoB577DFVV1frww8/HJa1GG+kb6JB9M7eDLp3795w286dOy/oZtAtW7aE29ra2iJuBj18+LC1f//+8PHzn//ckmS99dZbg97tPpYN1z5blmW1trZaqamp1uOPPz7MqzBPbm6u9fd///fhPweDQSszM3PIGxS/8Y1vRLTl5eWdc9PtM888Ez7v8/m46fYy77NlWVYgELDmz59vTZs2zTpx4sQwzn50udx73dPTE/Hf4v3791sZGRnWqlWrrLa2tmFejbkILKPUnXfead1yyy3Wnj17rN27d1s33HBDxOO2H330kfWXf/mX1p49e8JtDz30kDVp0iRr165d1t69e628vDwrLy9v0J/x+uuvX9VPCVnDtM/79++3vvjFL1qLFy+2urq6wsfV8gGwadMmy263Wy+88IJ14MAB68EHH7SSk5Mtr9drWZZlfetb37JKS0vD9W+++aY1btw465lnnrEOHjxoVVRUDPhYc3JysvXb3/7W+q//+i/rnnvu4bHmy7zPgUDAuvvuu63rrrvOevfddyN+d/v7+0dsnSYYjt/pz+MpIQLLqPWHP/zBWrRokXXNNddYSUlJVlFRkXX69Onw+Q8++MCSZL3++uvhtk8++cR6+OGHrS984QvWX/zFX1jf/OY3ra6urkF/BoFlePa5oqLCknTOcf3111/x9Y2U5557zpo0aZIVHx9v5ebmWm+//Xb43Ne//nVryZIlEfW/+tWvrClTpljx8fHWtGnTrO3bt0ecD4VC1urVqy2n02nZ7Xbr9ttvt9rb26/Yekx1Off57O/6QMef//5frS737/TnEVgsK8Y6e6MCAACAoXhKCAAAGI/AAgAAjEdgAQAAxiOwAAAA4xFYAACA8QgsAADAeAQWAABgPAILAAAwHoEFAAAYj8ACAACMR2ABAADGI7AAAADj/R8QIyCsWgIl/gAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(costs)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[0.6931471805599452,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf,\n", " inf]" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ "costs" ] } ], "metadata": { "colab": { "provenance": [] }, "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.9" } }, "nbformat": 4, "nbformat_minor": 4 }