{ "cells": [ { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "###分时段预测\n", "###丰都站\n", "###ses算法" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "from pandasql import sqldf \n", "from sqlalchemy import create_engine\n", "import statsmodels.api as sm \n", "import numpy as np\n", "\n", "mysql_user = 'root'\n", "mysql_password = 'qx123456'\n", "mysql_host = 'algo-dev-internal.li-ai.com.cn'\n", "mysql_port = 3306\n", "mysql_db = 'algo'\n", "\n", "db_engine = create_engine(\"mysql+pymysql://{}:{}@{}:{}/{}?charset=utf8\".format(mysql_user, mysql_password, mysql_host, mysql_port, mysql_db),\n", " pool_recycle=7200,pool_size=2)" ] }, { "cell_type": "code", "execution_count": 61, "metadata": {}, "outputs": [], "source": [ "sql=\"select * from algo_dwd_station_change_static where organ_code='fengdu2' and algo='tes' and pre_time=1 and data_period=3\"\n", "\n", "change_data=pd.read_sql(sql, db_engine)" ] }, { "cell_type": "code", "execution_count": 62, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(16969, 9)" ] }, "execution_count": 62, "metadata": {}, "output_type": "execute_result" } ], "source": [ "change_data.shape" ] }, { "cell_type": "code", "execution_count": 63, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
iddayslotc_b_countprepre_timeorgan_codealgodata_period
01702832022-02-2422:30-22:59001fengdu2tes3
11702862022-02-2423:00-23:29001fengdu2tes3
21702892022-02-2423:30-23:59001fengdu2tes3
31702922022-02-2500:00-00:29001fengdu2tes3
41702952022-02-2500:30-00:59001fengdu2tes3
\n", "
" ], "text/plain": [ " id day slot c_b_count pre pre_time organ_code algo \\\n", "0 170283 2022-02-24 22:30-22:59 0 0 1 fengdu2 tes \n", "1 170286 2022-02-24 23:00-23:29 0 0 1 fengdu2 tes \n", "2 170289 2022-02-24 23:30-23:59 0 0 1 fengdu2 tes \n", "3 170292 2022-02-25 00:00-00:29 0 0 1 fengdu2 tes \n", "4 170295 2022-02-25 00:30-00:59 0 0 1 fengdu2 tes \n", "\n", " data_period \n", "0 3 \n", "1 3 \n", "2 3 \n", "3 3 \n", "4 3 " ] }, "execution_count": 63, "metadata": {}, "output_type": "execute_result" } ], "source": [ "change_data.head()" ] }, { "cell_type": "code", "execution_count": 64, "metadata": {}, "outputs": [], "source": [ "change_data[\"delta\"]=change_data[\"pre\"]-change_data[\"c_b_count\"]" ] }, { "cell_type": "code", "execution_count": 65, "metadata": {}, "outputs": [], "source": [ "change_data[\"delta_abs\"]=change_data[\"delta\"].abs()" ] }, { "cell_type": "code", "execution_count": 75, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 75, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA3tklEQVR4nO3de3TU9Z3/8VcSkgkBJhE0CVkCZksLBEEECoxXwJBUU49otltaillBqTR0DTkLyh5ELtoIyk2NoisSu4VVaCutQCEDFCgSLkZSuVhqKxVbnMluEUZBJkPy/f3RM98fwzUTZkg+w/NxzpzjfL/v7yef9+eT0ZdzycRZlmUJAADAIPEtPQEAAIBwEWAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMZp09ITiJbGxkYdOXJEHTp0UFxcXEtPBwAANIFlWfriiy+UlZWl+PgLP88SswHmyJEjys7ObulpAACAZvj000/VpUuXC56P2QDToUMHSf9YAKfTGbFxA4GAqqqqlJ+fr8TExIiN25rEeo/0Z75Y7zHW+5Niv0f6az6fz6fs7Gz7v+MXErMBJviykdPpjHiASUlJkdPpjMlfSin2e6Q/88V6j7HenxT7PdLf5bvU2z94Ey8AADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHHCCjANDQ164oknlJOTo7Zt2+prX/uaZs+eLcuy7BrLsjR9+nR17txZbdu2VV5enj766KOQcY4eParRo0fL6XQqLS1N48aN05dffhlS88EHH+i2225TcnKysrOzNXfu3MtoEwAAxJKwAsycOXP08ssv68UXX9SHH36oOXPmaO7cuXrhhRfsmrlz5+r555/X4sWLtXPnTrVr104FBQU6deqUXTN69Gjt379fbrdbq1ev1tatWzV+/Hj7vM/nU35+vrp166aamho9++yzmjFjhl599dUItAwAAEwX1h+y2759u+69914VFhZKkq6//nr9z//8j3bt2iXpH8++LFy4UNOmTdO9994rSfrpT3+qjIwMrVq1SqNGjdKHH36odevWaffu3Ro4cKAk6YUXXtDdd9+t5557TllZWVq2bJnq6+v1+uuvKykpSb1791Ztba3mz58fEnQAAMDVKawAc/PNN+vVV1/VH//4R33jG9/Q73//e23btk3z58+XJB06dEgej0d5eXn2NampqRo8eLCqq6s1atQoVVdXKy0tzQ4vkpSXl6f4+Hjt3LlT9913n6qrq3X77bcrKSnJrikoKNCcOXP0+eef65prrjlnbn6/X36/377v8/kk/eOvBQYCgXDavKjgWJEcs7WJ9R7pz3yx3mOs9yfFfo/0d/ljX0pYAebxxx+Xz+dTz549lZCQoIaGBj399NMaPXq0JMnj8UiSMjIyQq7LyMiwz3k8HqWnp4dOok0bdezYMaQmJyfnnDGC584XYMrLyzVz5sxzjldVVSklJSWcNpvE7XZHfMzWJtZ7pD/zxXqPsd6fFPs90l/4Tp482aS6sALMihUrtGzZMi1fvtx+Wae0tFRZWVkqLi5u1kQjZerUqSorK7PvB78MKj8/P+LfheR2uzVixIiY/H4LKfZ7pD/zxXqPsd6fFPs90l/zBV9BuZSwAszkyZP1+OOPa9SoUZKkPn366JNPPlF5ebmKi4uVmZkpSfJ6vercubN9ndfrVb9+/SRJmZmZqqurCxn39OnTOnr0qH19ZmamvF5vSE3wfrDmbA6HQw6H45zjiYmJUfnlida4rUms90h/5ov1HmO9Pyn2e6S/5o3ZFGF9CunkyZOKjw+9JCEhQY2NjZKknJwcZWZmauPGjfZ5n8+nnTt3yuVySZJcLpeOHTummpoau2bTpk1qbGzU4MGD7ZqtW7eGvA7mdrvVo0eP8758BAAAri5hPQNzzz336Omnn1bXrl3Vu3dv7dmzR/Pnz9fYsWMl/eOrr0tLS/XUU0/p61//unJycvTEE08oKytLI0eOlCT16tVL3/rWt/Twww9r8eLFCgQCmjhxokaNGqWsrCxJ0ve//33NnDlT48aN02OPPaZ9+/Zp0aJFWrBgQWS7vww3zFgvf8PFv+q7NfnLM4UtPQUAACImrADzwgsv6IknntCPfvQj1dXVKSsrSz/84Q81ffp0u2bKlCk6ceKExo8fr2PHjunWW2/VunXrlJycbNcsW7ZMEydO1J133qn4+HgVFRXp+eeft8+npqaqqqpKJSUlGjBggK699lpNnz6dj1ADAABJYQaYDh06aOHChVq4cOEFa+Li4jRr1izNmjXrgjUdO3bU8uXLL/qz+vbtq9/97nfhTA8AAFwl+C4kAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGCcsALM9ddfr7i4uHNuJSUlkqRTp06ppKREnTp1Uvv27VVUVCSv1xsyxuHDh1VYWKiUlBSlp6dr8uTJOn36dEjN5s2b1b9/fzkcDnXv3l2VlZWX1yUAAIgpYQWY3bt367PPPrNvbrdbkvSd73xHkjRp0iS98847WrlypbZs2aIjR47o/vvvt69vaGhQYWGh6uvrtX37dr3xxhuqrKzU9OnT7ZpDhw6psLBQw4YNU21trUpLS/XQQw9p/fr1kegXAADEgDbhFF933XUh95955hl97Wtf0x133KHjx49ryZIlWr58uYYPHy5JWrp0qXr16qUdO3ZoyJAhqqqq0oEDB7RhwwZlZGSoX79+mj17th577DHNmDFDSUlJWrx4sXJycjRv3jxJUq9evbRt2zYtWLBABQUFEWobAACYrNnvgamvr9fPfvYzjR07VnFxcaqpqVEgEFBeXp5d07NnT3Xt2lXV1dWSpOrqavXp00cZGRl2TUFBgXw+n/bv32/XnDlGsCY4BgAAQFjPwJxp1apVOnbsmP7t3/5NkuTxeJSUlKS0tLSQuoyMDHk8HrvmzPASPB88d7Ean8+nr776Sm3btj3vfPx+v/x+v33f5/NJkgKBgAKBQPOaPI/gWI54K2JjXgnhrEGwNpLr1prQn/livcdY70+K/R7p7/LHvpRmB5glS5borrvuUlZWVnOHiKjy8nLNnDnznONVVVVKSUmJ+M+bPbAx4mNG09q1a8O+Jvgep1hFf+aL9R5jvT8p9nukv/CdPHmySXXNCjCffPKJNmzYoF/+8pf2sczMTNXX1+vYsWMhz8J4vV5lZmbaNbt27QoZK/gppTNrzv7kktfrldPpvOCzL5I0depUlZWV2fd9Pp+ys7OVn58vp9PZnDbPKxAIyO1264n34uVvjIvYuNG2b0bT3z8U7HHEiBFKTEyM4qxaBv2ZL9Z7jPX+pNjvkf6aL/gKyqU0K8AsXbpU6enpKiwstI8NGDBAiYmJ2rhxo4qKiiRJBw8e1OHDh+VyuSRJLpdLTz/9tOrq6pSeni7pH+nN6XQqNzfXrjn72QK3222PcSEOh0MOh+Oc44mJiVH55fE3xsnfYE6Aac4aRGvtWgv6M1+s9xjr/Umx3yP9NW/Mpgj7TbyNjY1aunSpiouL1abN/88/qampGjdunMrKyvTb3/5WNTU1evDBB+VyuTRkyBBJUn5+vnJzczVmzBj9/ve/1/r16zVt2jSVlJTY4eORRx7Rxx9/rClTpugPf/iDXnrpJa1YsUKTJk0Kd6oAACBGhf0MzIYNG3T48GGNHTv2nHMLFixQfHy8ioqK5Pf7VVBQoJdeesk+n5CQoNWrV2vChAlyuVxq166diouLNWvWLLsmJydHa9as0aRJk7Ro0SJ16dJFr732Gh+hBgAAtrADTH5+vizr/J/ASU5OVkVFhSoqKi54fbdu3S75htKhQ4dqz5494U4NAABcJfguJAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOM3+NmrgSrj+8TURH9ORYGnuIOmGGeuj8n1Wf3mm8NJFAIDLwjMwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOGEHmL/97W/6wQ9+oE6dOqlt27bq06eP3nvvPfu8ZVmaPn26OnfurLZt2yovL08fffRRyBhHjx7V6NGj5XQ6lZaWpnHjxunLL78Mqfnggw902223KTk5WdnZ2Zo7d24zWwQAALEmrADz+eef65ZbblFiYqJ+85vf6MCBA5o3b56uueYau2bu3Ll6/vnntXjxYu3cuVPt2rVTQUGBTp06ZdeMHj1a+/fvl9vt1urVq7V161aNHz/ePu/z+ZSfn69u3bqppqZGzz77rGbMmKFXX301Ai0DAADTtQmneM6cOcrOztbSpUvtYzk5OfY/W5alhQsXatq0abr33nslST/96U+VkZGhVatWadSoUfrwww+1bt067d69WwMHDpQkvfDCC7r77rv13HPPKSsrS8uWLVN9fb1ef/11JSUlqXfv3qqtrdX8+fNDgg4AALg6hRVgfv3rX6ugoEDf+c53tGXLFv3TP/2TfvSjH+nhhx+WJB06dEgej0d5eXn2NampqRo8eLCqq6s1atQoVVdXKy0tzQ4vkpSXl6f4+Hjt3LlT9913n6qrq3X77bcrKSnJrikoKNCcOXP0+eefhzzjE+T3++X3++37Pp9PkhQIBBQIBMJp86KCYznirYiNeSWEswbB2kiuW3M5EiK/zsG9i9YetvS6tab9i5ZY7zHW+5Niv0f6u/yxLyWsAPPxxx/r5ZdfVllZmf7zP/9Tu3fv1r//+78rKSlJxcXF8ng8kqSMjIyQ6zIyMuxzHo9H6enpoZNo00YdO3YMqTnzmZ0zx/R4POcNMOXl5Zo5c+Y5x6uqqpSSkhJOm00ye2BjxMeMprVr14Z9jdvtjsJMwjN3UPTGjtYeNmeto6E17F+0xXqPsd6fFPs90l/4Tp482aS6sAJMY2OjBg4cqJ/85CeSpJtuukn79u3T4sWLVVxcHP4sI2jq1KkqKyuz7/t8PmVnZys/P19OpzNiPycQCMjtduuJ9+Llb4yL2LjRtm9GQZNrgz2OGDFCiYmJUZzVpd0wY33Ex3TEW5o9sDFqexjOWkdDa9q/aIn1HmO9Pyn2e6S/5gu+gnIpYQWYzp07Kzc3N+RYr1699Itf/EKSlJmZKUnyer3q3LmzXeP1etWvXz+7pq6uLmSM06dP6+jRo/b1mZmZ8nq9ITXB+8GaszkcDjkcjnOOJyYmRuWXx98YJ3+DOQGmOWsQrbULRzTXOFp72NJrFtQa9i/aYr3HWO9Piv0e6a95YzZFWJ9CuuWWW3Tw4MGQY3/84x/VrVs3Sf94Q29mZqY2btxon/f5fNq5c6dcLpckyeVy6dixY6qpqbFrNm3apMbGRg0ePNiu2bp1a8jrYG63Wz169Djvy0cAAODqElaAmTRpknbs2KGf/OQn+tOf/qTly5fr1VdfVUlJiSQpLi5OpaWleuqpp/TrX/9ae/fu1QMPPKCsrCyNHDlS0j+esfnWt76lhx9+WLt27dK7776riRMnatSoUcrKypIkff/731dSUpLGjRun/fv366233tKiRYtCXiICAABXr7BeQvrmN7+pt99+W1OnTtWsWbOUk5OjhQsXavTo0XbNlClTdOLECY0fP17Hjh3TrbfeqnXr1ik5OdmuWbZsmSZOnKg777xT8fHxKioq0vPPP2+fT01NVVVVlUpKSjRgwABde+21mj59Oh+hBgAAksIMMJL07W9/W9/+9rcveD4uLk6zZs3SrFmzLljTsWNHLV++/KI/p2/fvvrd734X7vQAAMBVgO9CAgAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjhBVgZsyYobi4uJBbz5497fOnTp1SSUmJOnXqpPbt26uoqEherzdkjMOHD6uwsFApKSlKT0/X5MmTdfr06ZCazZs3q3///nI4HOrevbsqKyub3yEAAIg5YT8D07t3b3322Wf2bdu2bfa5SZMm6Z133tHKlSu1ZcsWHTlyRPfff799vqGhQYWFhaqvr9f27dv1xhtvqLKyUtOnT7drDh06pMLCQg0bNky1tbUqLS3VQw89pPXr119mqwAAIFa0CfuCNm2UmZl5zvHjx49ryZIlWr58uYYPHy5JWrp0qXr16qUdO3ZoyJAhqqqq0oEDB7RhwwZlZGSoX79+mj17th577DHNmDFDSUlJWrx4sXJycjRv3jxJUq9evbRt2zYtWLBABQUFl9kuAACIBWEHmI8++khZWVlKTk6Wy+VSeXm5unbtqpqaGgUCAeXl5dm1PXv2VNeuXVVdXa0hQ4aourpaffr0UUZGhl1TUFCgCRMmaP/+/brppptUXV0dMkawprS09KLz8vv98vv99n2fzydJCgQCCgQC4bZ5QcGxHPFWxMa8EsJZg2BtJNetuRwJkV/n4N5Faw9bet1a0/5FS6z3GOv9SbHfI/1d/tiXElaAGTx4sCorK9WjRw999tlnmjlzpm677Tbt27dPHo9HSUlJSktLC7kmIyNDHo9HkuTxeELCS/B88NzFanw+n7766iu1bdv2vHMrLy/XzJkzzzleVVWllJSUcNpsktkDGyM+ZjStXbs27GvcbncUZhKeuYOiN3a09rA5ax0NrWH/oi3We4z1/qTY75H+wnfy5Mkm1YUVYO666y77n/v27avBgwerW7duWrFixQWDxZUydepUlZWV2fd9Pp+ys7OVn58vp9MZsZ8TCATkdrv1xHvx8jfGRWzcaNs3o+kvvwV7HDFihBITE6M4q0u7YUbk3/vkiLc0e2Bj1PYwnLWOhta0f9ES6z3Gen9S7PdIf80XfAXlUsJ+CelMaWlp+sY3vqE//elPGjFihOrr63Xs2LGQZ2G8Xq/9npnMzEzt2rUrZIzgp5TOrDn7k0ter1dOp/OiIcnhcMjhcJxzPDExMSq/PP7GOPkbzAkwzVmDaK1dOKK5xtHaw5Zes6DWsH/RFus9xnp/Uuz3SH/NG7MpLuvvwHz55Zf685//rM6dO2vAgAFKTEzUxo0b7fMHDx7U4cOH5XK5JEkul0t79+5VXV2dXeN2u+V0OpWbm2vXnDlGsCY4BgAAQFgB5j/+4z+0ZcsW/eUvf9H27dt13333KSEhQd/73veUmpqqcePGqaysTL/97W9VU1OjBx98UC6XS0OGDJEk5efnKzc3V2PGjNHvf/97rV+/XtOmTVNJSYn97Mkjjzyijz/+WFOmTNEf/vAHvfTSS1qxYoUmTZoU+e4BAICRwnoJ6a9//au+973v6e9//7uuu+463XrrrdqxY4euu+46SdKCBQsUHx+voqIi+f1+FRQU6KWXXrKvT0hI0OrVqzVhwgS5XC61a9dOxcXFmjVrll2Tk5OjNWvWaNKkSVq0aJG6dOmi1157jY9QAwAAW1gB5s0337zo+eTkZFVUVKiiouKCNd26dbvkpzSGDh2qPXv2hDM1AABwFeG7kAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcS4rwDzzzDOKi4tTaWmpfezUqVMqKSlRp06d1L59exUVFcnr9YZcd/jwYRUWFiolJUXp6emaPHmyTp8+HVKzefNm9e/fXw6HQ927d1dlZeXlTBUAAMSQZgeY3bt365VXXlHfvn1Djk+aNEnvvPOOVq5cqS1btujIkSO6//777fMNDQ0qLCxUfX29tm/frjfeeEOVlZWaPn26XXPo0CEVFhZq2LBhqq2tVWlpqR566CGtX7++udMFAAAxpFkB5ssvv9To0aP1X//1X7rmmmvs48ePH9eSJUs0f/58DR8+XAMGDNDSpUu1fft27dixQ5JUVVWlAwcO6Gc/+5n69eunu+66S7Nnz1ZFRYXq6+slSYsXL1ZOTo7mzZunXr16aeLEifqXf/kXLViwIAItAwAA07VpzkUlJSUqLCxUXl6ennrqKft4TU2NAoGA8vLy7GM9e/ZU165dVV1drSFDhqi6ulp9+vRRRkaGXVNQUKAJEyZo//79uummm1RdXR0yRrDmzJeqzub3++X3++37Pp9PkhQIBBQIBJrT5nkFx3LEWxEb80oIZw2CtZFct+ZyJER+nYN7F609bOl1a037Fy2x3mOs9yfFfo/0d/ljX0rYAebNN9/U+++/r927d59zzuPxKCkpSWlpaSHHMzIy5PF47Jozw0vwfPDcxWp8Pp+++uortW3b9pyfXV5erpkzZ55zvKqqSikpKU1vsIlmD2yM+JjRtHbt2rCvcbvdUZhJeOYOit7Y0drD5qx1NLSG/Yu2WO8x1vuTYr9H+gvfyZMnm1QXVoD59NNP9eijj8rtdis5OblZE4uWqVOnqqyszL7v8/mUnZ2t/Px8OZ3OiP2cQCAgt9utJ96Ll78xLmLjRtu+GQVNrg32OGLECCUmJkZxVpd2w4zIv+/JEW9p9sDGqO1hOGsdDa1p/6Il1nuM9f6k2O+R/pov+ArKpYQVYGpqalRXV6f+/fvbxxoaGrR161a9+OKLWr9+verr63Xs2LGQZ2G8Xq8yMzMlSZmZmdq1a1fIuMFPKZ1Zc/Ynl7xer5xO53mffZEkh8Mhh8NxzvHExMSo/PL4G+PkbzAnwDRnDaK1duGI5hpHaw9bes2CWsP+RVus9xjr/Umx3yP9NW/MpgjrTbx33nmn9u7dq9raWvs2cOBAjR492v7nxMREbdy40b7m4MGDOnz4sFwulyTJ5XJp7969qqurs2vcbrecTqdyc3PtmjPHCNYExwAAAFe3sJ6B6dChg2644YaQY+3atVOnTp3s4+PGjVNZWZk6duwop9OpH//4x3K5XBoyZIgkKT8/X7m5uRozZozmzp0rj8ejadOmqaSkxH4G5ZFHHtGLL76oKVOmaOzYsdq0aZNWrFihNWvWRKJnAABguGZ9CuliFixYoPj4eBUVFcnv96ugoEAvvfSSfT4hIUGrV6/WhAkT5HK51K5dOxUXF2vWrFl2TU5OjtasWaNJkyZp0aJF6tKli1577TUVFLTsewsAAEDrcNkBZvPmzSH3k5OTVVFRoYqKigte061bt0t+UmPo0KHas2fP5U4PAADEIL4LCQAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACME1aAefnll9W3b185nU45nU65XC795je/sc+fOnVKJSUl6tSpk9q3b6+ioiJ5vd6QMQ4fPqzCwkKlpKQoPT1dkydP1unTp0NqNm/erP79+8vhcKh79+6qrKxsfocAACDmhBVgunTpomeeeUY1NTV67733NHz4cN17773av3+/JGnSpEl65513tHLlSm3ZskVHjhzR/fffb1/f0NCgwsJC1dfXa/v27XrjjTdUWVmp6dOn2zWHDh1SYWGhhg0bptraWpWWluqhhx7S+vXrI9QyAAAwXZtwiu+5556Q+08//bRefvll7dixQ126dNGSJUu0fPlyDR8+XJK0dOlS9erVSzt27NCQIUNUVVWlAwcOaMOGDcrIyFC/fv00e/ZsPfbYY5oxY4aSkpK0ePFi5eTkaN68eZKkXr16adu2bVqwYIEKCgoi1DYAADBZWAHmTA0NDVq5cqVOnDghl8ulmpoaBQIB5eXl2TU9e/ZU165dVV1drSFDhqi6ulp9+vRRRkaGXVNQUKAJEyZo//79uummm1RdXR0yRrCmtLT0ovPx+/3y+/32fZ/PJ0kKBAIKBALNbfMcwbEc8VbExrwSwlmDYG0k1625HAmRX+fg3kVrD1t63VrT/kVLrPcY6/1Jsd8j/V3+2JcSdoDZu3evXC6XTp06pfbt2+vtt99Wbm6uamtrlZSUpLS0tJD6jIwMeTweSZLH4wkJL8HzwXMXq/H5fPrqq6/Utm3b886rvLxcM2fOPOd4VVWVUlJSwm3zkmYPbIz4mNG0du3asK9xu91RmEl45g6K3tjR2sPmrHU0tIb9i7ZY7zHW+5Niv0f6C9/JkyebVBd2gOnRo4dqa2t1/Phx/fznP1dxcbG2bNkS9gQjberUqSorK7Pv+3w+ZWdnKz8/X06nM2I/JxAIyO1264n34uVvjIvYuNG2b0bTX34L9jhixAglJiZGcVaXdsOMyL/3yRFvafbAxqjtYThrHQ2taf+iJdZ7jPX+pNjvkf6aL/gKyqWEHWCSkpLUvXt3SdKAAQO0e/duLVq0SN/97ndVX1+vY8eOhTwL4/V6lZmZKUnKzMzUrl27QsYLfkrpzJqzP7nk9XrldDov+OyLJDkcDjkcjnOOJyYmRuWXx98YJ3+DOQGmOWsQrbULRzTXOFp72NJrFtQa9i/aYr3HWO9Piv0e6a95YzbFZf8dmMbGRvn9fg0YMECJiYnauHGjfe7gwYM6fPiwXC6XJMnlcmnv3r2qq6uza9xut5xOp3Jzc+2aM8cI1gTHAAAACOsZmKlTp+quu+5S165d9cUXX2j58uXavHmz1q9fr9TUVI0bN05lZWXq2LGjnE6nfvzjH8vlcmnIkCGSpPz8fOXm5mrMmDGaO3euPB6Ppk2bppKSEvvZk0ceeUQvvviipkyZorFjx2rTpk1asWKF1qxZE/nuAQCAkcIKMHV1dXrggQf02WefKTU1VX379tX69es1YsQISdKCBQsUHx+voqIi+f1+FRQU6KWXXrKvT0hI0OrVqzVhwgS5XC61a9dOxcXFmjVrll2Tk5OjNWvWaNKkSVq0aJG6dOmi1157jY9QAwAAW1gBZsmSJRc9n5ycrIqKClVUVFywplu3bpf8lMbQoUO1Z8+ecKYGAACuInwXEgAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwTlgBpry8XN/85jfVoUMHpaena+TIkTp48GBIzalTp1RSUqJOnTqpffv2KioqktfrDak5fPiwCgsLlZKSovT0dE2ePFmnT58Oqdm8ebP69+8vh8Oh7t27q7KysnkdAgCAmBNWgNmyZYtKSkq0Y8cOud1uBQIB5efn68SJE3bNpEmT9M4772jlypXasmWLjhw5ovvvv98+39DQoMLCQtXX12v79u164403VFlZqenTp9s1hw4dUmFhoYYNG6ba2lqVlpbqoYce0vr16yPQMgAAMF2bcIrXrVsXcr+yslLp6emqqanR7bffruPHj2vJkiVavny5hg8fLklaunSpevXqpR07dmjIkCGqqqrSgQMHtGHDBmVkZKhfv36aPXu2HnvsMc2YMUNJSUlavHixcnJyNG/ePElSr169tG3bNi1YsEAFBQURah0AAJgqrABztuPHj0uSOnbsKEmqqalRIBBQXl6eXdOzZ0917dpV1dXVGjJkiKqrq9WnTx9lZGTYNQUFBZowYYL279+vm266SdXV1SFjBGtKS0svOBe/3y+/32/f9/l8kqRAIKBAIHA5bYYIjuWItyI25pUQzhoEayO5bs3lSIj8Ogf3Llp72NLr1pr2L1pivcdY70+K/R7p7/LHvpRmB5jGxkaVlpbqlltu0Q033CBJ8ng8SkpKUlpaWkhtRkaGPB6PXXNmeAmeD567WI3P59NXX32ltm3bnjOf8vJyzZw585zjVVVVSklJaV6TFzF7YGPEx4ymtWvXhn2N2+2OwkzCM3dQ9MaO1h42Z62joTXsX7TFeo+x3p8U+z3SX/hOnjzZpLpmB5iSkhLt27dP27Zta+4QETV16lSVlZXZ930+n7Kzs5Wfny+n0xmxnxMIBOR2u/XEe/HyN8ZFbNxo2zej6S+9BXscMWKEEhMTozirS7thRuTf9+SItzR7YGPU9jCctY6G1rR/0RLrj8OraQ9jtUf6a77gKyiX0qwAM3HiRK1evVpbt25Vly5d7OOZmZmqr6/XsWPHQp6F8Xq9yszMtGt27doVMl7wU0pn1pz9ySWv1yun03neZ18kyeFwyOFwnHM8MTExKr88/sY4+RvM+Rdnc9YgWmsXjmiucbT2sKXXLKg17F+0xfrj8GrYw1jvkf6aN2ZThPUpJMuyNHHiRL399tvatGmTcnJyQs4PGDBAiYmJ2rhxo33s4MGDOnz4sFwulyTJ5XJp7969qqurs2vcbrecTqdyc3PtmjPHCNYExwAAAFe3sJ6BKSkp0fLly/WrX/1KHTp0sN+zkpqaqrZt2yo1NVXjxo1TWVmZOnbsKKfTqR//+MdyuVwaMmSIJCk/P1+5ubkaM2aM5s6dK4/Ho2nTpqmkpMR+BuWRRx7Riy++qClTpmjs2LHatGmTVqxYoTVr1kS4fQAAYKKwnoF5+eWXdfz4cQ0dOlSdO3e2b2+99ZZds2DBAn37299WUVGRbr/9dmVmZuqXv/ylfT4hIUGrV69WQkKCXC6XfvCDH+iBBx7QrFmz7JqcnBytWbNGbrdbN954o+bNm6fXXnuNj1ADAABJYT4DY1mX/thpcnKyKioqVFFRccGabt26XfKTGkOHDtWePXvCmR4AALhK8F1IAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOAQYAABgHAIMAAAwDgEGAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGCcsAPM1q1bdc899ygrK0txcXFatWpVyHnLsjR9+nR17txZbdu2VV5enj766KOQmqNHj2r06NFyOp1KS0vTuHHj9OWXX4bUfPDBB7rtttuUnJys7OxszZ07N/zuAABATAo7wJw4cUI33nijKioqznt+7ty5ev7557V48WLt3LlT7dq1U0FBgU6dOmXXjB49Wvv375fb7dbq1au1detWjR8/3j7v8/mUn5+vbt26qaamRs8++6xmzJihV199tRktAgCAWNMm3Avuuusu3XXXXec9Z1mWFi5cqGnTpunee++VJP30pz9VRkaGVq1apVGjRunDDz/UunXrtHv3bg0cOFCS9MILL+juu+/Wc889p6ysLC1btkz19fV6/fXXlZSUpN69e6u2tlbz588PCToAAODqFHaAuZhDhw7J4/EoLy/PPpaamqrBgwerurpao0aNUnV1tdLS0uzwIkl5eXmKj4/Xzp07dd9996m6ulq33367kpKS7JqCggLNmTNHn3/+ua655ppzfrbf75ff77fv+3w+SVIgEFAgEIhYj8GxHPFWxMa8EsJZg2BtJNetuRwJkV/n4N5Faw9bet1a0/5FS6w/Dq+mPYzVHunv8se+lIgGGI/HI0nKyMgIOZ6RkWGf83g8Sk9PD51Emzbq2LFjSE1OTs45YwTPnS/AlJeXa+bMmeccr6qqUkpKSjM7urDZAxsjPmY0rV27Nuxr3G53FGYSnrmDojd2tPawOWsdDa1h/6It1h+HV8MexnqP9Be+kydPNqkuogGmJU2dOlVlZWX2fZ/Pp+zsbOXn58vpdEbs5wQCAbndbj3xXrz8jXERGzfa9s0oaHJtsMcRI0YoMTExirO6tBtmrI/4mI54S7MHNkZtD8NZ62hoTfsXLbH+OLya9jBWe6S/5gu+gnIpEQ0wmZmZkiSv16vOnTvbx71er/r162fX1NXVhVx3+vRpHT161L4+MzNTXq83pCZ4P1hzNofDIYfDcc7xxMTEqPzy+Bvj5G8w51+czVmDaK1dOKK5xtHaw5Zes6DWsH/RFuuPw6thD2O9R/pr3phNEdG/A5OTk6PMzExt3LjRPubz+bRz5065XC5Jksvl0rFjx1RTU2PXbNq0SY2NjRo8eLBds3Xr1pDXwdxut3r06HHel48AAMDVJewA8+WXX6q2tla1tbWS/vHG3draWh0+fFhxcXEqLS3VU089pV//+tfau3evHnjgAWVlZWnkyJGSpF69eulb3/qWHn74Ye3atUvvvvuuJk6cqFGjRikrK0uS9P3vf19JSUkaN26c9u/fr7feekuLFi0KeYkIAABcvcJ+Cem9997TsGHD7PvBUFFcXKzKykpNmTJFJ06c0Pjx43Xs2DHdeuutWrdunZKTk+1rli1bpokTJ+rOO+9UfHy8ioqK9Pzzz9vnU1NTVVVVpZKSEg0YMEDXXnutpk+fzkeoAQCApGYEmKFDh8qyLvzRxbi4OM2aNUuzZs26YE3Hjh21fPnyi/6cvn376ne/+1240wMAAFcBvgsJAAAYhwADAACMQ4ABAADGIcAAAADjEGAAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMA4BBgAAGIcAAwAAjEOAAQAAxiHAAAAA4xBgAACAcQgwAADAOG1aegIAgNbr+sfXRGVcR4KluYOkG2asl78hLqJj/+WZwoiOh9aJZ2AAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDitOsBUVFTo+uuvV3JysgYPHqxdu3a19JQAAEAr0GoDzFtvvaWysjI9+eSTev/993XjjTeqoKBAdXV1LT01AADQwlptgJk/f74efvhhPfjgg8rNzdXixYuVkpKi119/vaWnBgAAWlir/Dbq+vp61dTUaOrUqfax+Ph45eXlqbq6+rzX+P1++f1++/7x48clSUePHlUgEIjY3AKBgE6ePKk2gXg1NEb2G1Sj6e9//3uTa4M9/v3vf1diYmIUZ3VpbU6fiPyYjZZOnmyM2h6Gs9bR0Jr2L1pi/XHYmvYwGo9BKbqPw5Z+DErh7+Hg8o1XYFaR44i3NO2mxqj8jn7xxReSJMuyLl5otUJ/+9vfLEnW9u3bQ45PnjzZGjRo0HmvefLJJy1J3Lhx48aNG7cYuH366acXzQqt8hmY5pg6darKysrs+42NjTp69Kg6deqkuLjIpXufz6fs7Gx9+umncjqdERu3NYn1HunPfLHeY6z3J8V+j/TXfJZl6YsvvlBWVtZF61plgLn22muVkJAgr9cbctzr9SozM/O81zgcDjkcjpBjaWlp0ZqinE5nTP5SninWe6Q/88V6j7HenxT7PdJf86Smpl6yplW+iTcpKUkDBgzQxo3//zXBxsZGbdy4US6XqwVnBgAAWoNW+QyMJJWVlam4uFgDBw7UoEGDtHDhQp04cUIPPvhgS08NAAC0sFYbYL773e/qf//3fzV9+nR5PB7169dP69atU0ZGRovOy+Fw6Mknnzzn5apYEus90p/5Yr3HWO9Piv0e6S/64izrUp9TAgAAaF1a5XtgAAAALoYAAwAAjEOAAQAAxiHAAAAA4xBgzqOiokLXX3+9kpOTNXjwYO3ateui9StXrlTPnj2VnJysPn36aO3atVdops0XTo+VlZWKi4sLuSUnJ1/B2YZn69atuueee5SVlaW4uDitWrXqktds3rxZ/fv3l8PhUPfu3VVZWRn1eTZXuP1t3rz5nP2Li4uTx+O5MhMOU3l5ub75zW+qQ4cOSk9P18iRI3Xw4MFLXmfK47A5/Zn2GHz55ZfVt29f+4+cuVwu/eY3v7noNabsnxR+f6bt39meeeYZxcXFqbS09KJ1V3oPCTBneeutt1RWVqYnn3xS77//vm688UYVFBSorq7uvPXbt2/X9773PY0bN0579uzRyJEjNXLkSO3bt+8Kz7zpwu1R+sdfW/zss8/s2yeffHIFZxyeEydO6MYbb1RFRUWT6g8dOqTCwkINGzZMtbW1Ki0t1UMPPaT169dHeabNE25/QQcPHgzZw/T09CjN8PJs2bJFJSUl2rFjh9xutwKBgPLz83XixIW/VNCkx2Fz+pPMegx26dJFzzzzjGpqavTee+9p+PDhuvfee7V///7z1pu0f1L4/Ulm7d+Zdu/erVdeeUV9+/a9aF2L7GFkvn4xdgwaNMgqKSmx7zc0NFhZWVlWeXn5eev/9V//1SosLAw5NnjwYOuHP/xhVOd5OcLtcenSpVZqauoVml1kSbLefvvti9ZMmTLF6t27d8ix7373u1ZBQUEUZxYZTenvt7/9rSXJ+vzzz6/InCKtrq7OkmRt2bLlgjUmPg6DmtKfyY/BoGuuucZ67bXXznvO5P0Lulh/pu7fF198YX3961+33G63dccdd1iPPvroBWtbYg95BuYM9fX1qqmpUV5enn0sPj5eeXl5qq6uPu811dXVIfWSVFBQcMH6ltacHiXpyy+/VLdu3ZSdnX3J/9MwjWl72Fz9+vVT586dNWLECL377rstPZ0mO378uCSpY8eOF6wxeQ+b0p9k7mOwoaFBb775pk6cOHHBr4Ixef+a0p9k5v6VlJSosLDwnL05n5bYQwLMGf7v//5PDQ0N5/y134yMjAu+X8Dj8YRV39Ka02OPHj30+uuv61e/+pV+9rOfqbGxUTfffLP++te/XokpR92F9tDn8+mrr75qoVlFTufOnbV48WL94he/0C9+8QtlZ2dr6NChev/991t6apfU2Nio0tJS3XLLLbrhhhsuWGfa4zCoqf2Z+Bjcu3ev2rdvL4fDoUceeURvv/22cnNzz1tr4v6F05+J+/fmm2/q/fffV3l5eZPqW2IPW+1XCaD1cLlcIf9ncfPNN6tXr1565ZVXNHv27BacGZqiR48e6tGjh33/5ptv1p///GctWLBA//3f/92CM7u0kpIS7du3T9u2bWvpqURFU/sz8THYo0cP1dbW6vjx4/r5z3+u4uJibdmy5YL/kTdNOP2Ztn+ffvqpHn30Ubnd7lb9ZmMCzBmuvfZaJSQkyOv1hhz3er3KzMw87zWZmZlh1be05vR4tsTERN10003605/+FI0pXnEX2kOn06m2bdu20Kyia9CgQa0+FEycOFGrV6/W1q1b1aVLl4vWmvY4lMLr72wmPAaTkpLUvXt3SdKAAQO0e/duLVq0SK+88so5tSbuXzj9na21719NTY3q6urUv39/+1hDQ4O2bt2qF198UX6/XwkJCSHXtMQe8hLSGZKSkjRgwABt3LjRPtbY2KiNGzde8LVNl8sVUi9Jbrf7oq+FtqTm9Hi2hoYG7d27V507d47WNK8o0/YwEmpra1vt/lmWpYkTJ+rtt9/Wpk2blJOTc8lrTNrD5vR3NhMfg42NjfL7/ec9Z9L+XcjF+jtba9+/O++8U3v37lVtba19GzhwoEaPHq3a2tpzwovUQnsYtbcHG+rNN9+0HA6HVVlZaR04cMAaP368lZaWZnk8HsuyLGvMmDHW448/bte/++67Vps2baznnnvO+vDDD60nn3zSSkxMtPbu3dtSLVxSuD3OnDnTWr9+vfXnP//ZqqmpsUaNGmUlJydb+/fvb6kWLuqLL76w9uzZY+3Zs8eSZM2fP9/as2eP9cknn1iWZVmPP/64NWbMGLv+448/tlJSUqzJkydbH374oVVRUWElJCRY69ata6kWLirc/hYsWGCtWrXK+uijj6y9e/dajz76qBUfH29t2LChpVq4qAkTJlipqanW5s2brc8++8y+nTx50q4x+XHYnP5Meww+/vjj1pYtW6xDhw5ZH3zwgfX4449bcXFxVlVVlWVZZu+fZYXfn2n7dz5nfwqpNewhAeY8XnjhBatr165WUlKSNWjQIGvHjh32uTvuuMMqLi4OqV+xYoX1jW98w0pKSrJ69+5trVmz5grPOHzh9FhaWmrXZmRkWHfffbf1/vvvt8Csmyb4seGzb8GeiouLrTvuuOOca/r162clJSVZ//zP/2wtXbr0is+7qcLtb86cOdbXvvY1Kzk52erYsaM1dOhQa9OmTS0z+SY4X2+SQvbE5Mdhc/oz7TE4duxYq1u3blZSUpJ13XXXWXfeeaf9H3fLMnv/LCv8/kzbv/M5O8C0hj2MsyzLit7zOwAAAJHHe2AAAIBxCDAAAMA4BBgAAGAcAgwAADAOAQYAABiHAAMAAIxDgAEAAMYhwAAAAOMQYAAAgHEIMAAAwDgEGAAAYBwCDAAAMM7/A46FLYxtp9qIAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "change_data[\"delta_abs\"].hist()" ] }, { "cell_type": "code", "execution_count": 71, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1.0" ] }, "execution_count": 71, "metadata": {}, "output_type": "execute_result" } ], "source": [ "change_data[\"delta_abs\"].quantile(0.85)" ] }, { "cell_type": "code", "execution_count": 72, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.697742943013731" ] }, "execution_count": 72, "metadata": {}, "output_type": "execute_result" } ], "source": [ "change_data[\"delta_abs\"].mean()" ] }, { "cell_type": "code", "execution_count": 73, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "4" ] }, "execution_count": 73, "metadata": {}, "output_type": "execute_result" } ], "source": [ "change_data[\"delta_abs\"].max()" ] }, { "cell_type": "code", "execution_count": 74, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1.0" ] }, "execution_count": 74, "metadata": {}, "output_type": "execute_result" } ], "source": [ "change_data[\"delta_abs\"].median()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "py3916", "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.9.16" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }