2.本体系：线性回归.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "21b1f85c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>D1</th>\n",
       "      <th>D2</th>\n",
       "      <th>D3</th>\n",
       "      <th>D4</th>\n",
       "      <th>D5</th>\n",
       "      <th>D6</th>\n",
       "      <th>D7</th>\n",
       "      <th>D8</th>\n",
       "      <th>D9</th>\n",
       "      <th>D10</th>\n",
       "      <th>...</th>\n",
       "      <th>D20</th>\n",
       "      <th>D21</th>\n",
       "      <th>D22</th>\n",
       "      <th>D23</th>\n",
       "      <th>D24</th>\n",
       "      <th>D25</th>\n",
       "      <th>D26</th>\n",
       "      <th>D27</th>\n",
       "      <th>D28</th>\n",
       "      <th>y</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>233.762715</td>\n",
       "      <td>116.898844</td>\n",
       "      <td>116.863871</td>\n",
       "      <td>0.034973</td>\n",
       "      <td>30.917600</td>\n",
       "      <td>84706.25146</td>\n",
       "      <td>-84737.16906</td>\n",
       "      <td>-84706.25146</td>\n",
       "      <td>-338701.3354</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.040350</td>\n",
       "      <td>-45.693931</td>\n",
       "      <td>954.306069</td>\n",
       "      <td>-22.218856</td>\n",
       "      <td>5.469140</td>\n",
       "      <td>-338701.3354</td>\n",
       "      <td>-1.439090e+15</td>\n",
       "      <td>0.000013</td>\n",
       "      <td>-1.000127</td>\n",
       "      <td>-42.935489</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>233.754944</td>\n",
       "      <td>116.890477</td>\n",
       "      <td>116.864467</td>\n",
       "      <td>0.026010</td>\n",
       "      <td>30.913406</td>\n",
       "      <td>84703.02826</td>\n",
       "      <td>-84733.94167</td>\n",
       "      <td>-84703.02826</td>\n",
       "      <td>-338688.4594</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.040175</td>\n",
       "      <td>-45.815157</td>\n",
       "      <td>954.184843</td>\n",
       "      <td>-22.224564</td>\n",
       "      <td>5.432392</td>\n",
       "      <td>-338688.4594</td>\n",
       "      <td>-1.438920e+15</td>\n",
       "      <td>0.000009</td>\n",
       "      <td>-1.000118</td>\n",
       "      <td>-33.770706</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>233.756958</td>\n",
       "      <td>116.904940</td>\n",
       "      <td>116.852018</td>\n",
       "      <td>0.052922</td>\n",
       "      <td>30.843638</td>\n",
       "      <td>84702.00879</td>\n",
       "      <td>-84732.85243</td>\n",
       "      <td>-84702.00879</td>\n",
       "      <td>-338684.6606</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.037591</td>\n",
       "      <td>-45.967968</td>\n",
       "      <td>954.032032</td>\n",
       "      <td>-22.224253</td>\n",
       "      <td>5.414609</td>\n",
       "      <td>-338684.6606</td>\n",
       "      <td>-1.438870e+15</td>\n",
       "      <td>0.000006</td>\n",
       "      <td>-1.000107</td>\n",
       "      <td>-32.560866</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>233.757343</td>\n",
       "      <td>116.900775</td>\n",
       "      <td>116.856568</td>\n",
       "      <td>0.044207</td>\n",
       "      <td>30.849927</td>\n",
       "      <td>84703.15045</td>\n",
       "      <td>-84734.00037</td>\n",
       "      <td>-84703.15045</td>\n",
       "      <td>-338689.2021</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.037834</td>\n",
       "      <td>-45.913815</td>\n",
       "      <td>954.086185</td>\n",
       "      <td>-22.229211</td>\n",
       "      <td>5.374452</td>\n",
       "      <td>-338689.2021</td>\n",
       "      <td>-1.438930e+15</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>-1.000104</td>\n",
       "      <td>-23.811494</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>233.771601</td>\n",
       "      <td>116.910762</td>\n",
       "      <td>116.860839</td>\n",
       "      <td>0.049922</td>\n",
       "      <td>30.854660</td>\n",
       "      <td>84709.04221</td>\n",
       "      <td>-84739.89687</td>\n",
       "      <td>-84709.04221</td>\n",
       "      <td>-338712.7502</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.038212</td>\n",
       "      <td>-45.820667</td>\n",
       "      <td>954.179333</td>\n",
       "      <td>-22.221312</td>\n",
       "      <td>5.379066</td>\n",
       "      <td>-338712.7502</td>\n",
       "      <td>-1.439230e+15</td>\n",
       "      <td>0.000007</td>\n",
       "      <td>-1.000113</td>\n",
       "      <td>-27.856424</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>233.769472</td>\n",
       "      <td>116.904063</td>\n",
       "      <td>116.865409</td>\n",
       "      <td>0.038654</td>\n",
       "      <td>30.814714</td>\n",
       "      <td>84707.62008</td>\n",
       "      <td>-84738.43479</td>\n",
       "      <td>-84707.62008</td>\n",
       "      <td>-338707.2214</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.037772</td>\n",
       "      <td>-45.975482</td>\n",
       "      <td>954.024518</td>\n",
       "      <td>-22.211994</td>\n",
       "      <td>5.368161</td>\n",
       "      <td>-338707.2214</td>\n",
       "      <td>-1.439160e+15</td>\n",
       "      <td>0.000006</td>\n",
       "      <td>-1.000095</td>\n",
       "      <td>-32.063878</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>233.763229</td>\n",
       "      <td>116.902341</td>\n",
       "      <td>116.860888</td>\n",
       "      <td>0.041454</td>\n",
       "      <td>30.792690</td>\n",
       "      <td>84705.32736</td>\n",
       "      <td>-84736.12005</td>\n",
       "      <td>-84705.32736</td>\n",
       "      <td>-338698.1387</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.037038</td>\n",
       "      <td>-45.931454</td>\n",
       "      <td>954.068546</td>\n",
       "      <td>-22.218519</td>\n",
       "      <td>5.330169</td>\n",
       "      <td>-338698.1387</td>\n",
       "      <td>-1.439050e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000094</td>\n",
       "      <td>-29.522463</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>233.802668</td>\n",
       "      <td>116.925918</td>\n",
       "      <td>116.876749</td>\n",
       "      <td>0.049169</td>\n",
       "      <td>30.649775</td>\n",
       "      <td>84719.30084</td>\n",
       "      <td>-84749.95062</td>\n",
       "      <td>-84719.30084</td>\n",
       "      <td>-338754.6043</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.033831</td>\n",
       "      <td>-46.325084</td>\n",
       "      <td>953.674916</td>\n",
       "      <td>-22.192825</td>\n",
       "      <td>5.369025</td>\n",
       "      <td>-338754.6043</td>\n",
       "      <td>-1.439770e+15</td>\n",
       "      <td>0.000033</td>\n",
       "      <td>-1.000085</td>\n",
       "      <td>-33.272463</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>233.733487</td>\n",
       "      <td>116.876740</td>\n",
       "      <td>116.856747</td>\n",
       "      <td>0.019993</td>\n",
       "      <td>30.971601</td>\n",
       "      <td>84690.28825</td>\n",
       "      <td>-84721.25985</td>\n",
       "      <td>-84690.28825</td>\n",
       "      <td>-338637.2666</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.041815</td>\n",
       "      <td>-46.221557</td>\n",
       "      <td>953.778443</td>\n",
       "      <td>-22.197462</td>\n",
       "      <td>5.580163</td>\n",
       "      <td>-338637.2666</td>\n",
       "      <td>-1.438270e+15</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>-1.000105</td>\n",
       "      <td>-41.417543</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>233.693647</td>\n",
       "      <td>116.866452</td>\n",
       "      <td>116.827195</td>\n",
       "      <td>0.039257</td>\n",
       "      <td>31.075569</td>\n",
       "      <td>84678.89210</td>\n",
       "      <td>-84709.96767</td>\n",
       "      <td>-84678.89210</td>\n",
       "      <td>-338591.2661</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.039874</td>\n",
       "      <td>-45.725674</td>\n",
       "      <td>954.274326</td>\n",
       "      <td>-22.227629</td>\n",
       "      <td>5.248522</td>\n",
       "      <td>-338591.2661</td>\n",
       "      <td>-1.437680e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000139</td>\n",
       "      <td>-42.223893</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>233.770927</td>\n",
       "      <td>116.932869</td>\n",
       "      <td>116.838058</td>\n",
       "      <td>0.094811</td>\n",
       "      <td>30.917504</td>\n",
       "      <td>84708.85367</td>\n",
       "      <td>-84739.77117</td>\n",
       "      <td>-84708.85367</td>\n",
       "      <td>-338711.7447</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.040042</td>\n",
       "      <td>-45.823397</td>\n",
       "      <td>954.176603</td>\n",
       "      <td>-22.222212</td>\n",
       "      <td>5.454897</td>\n",
       "      <td>-338711.7447</td>\n",
       "      <td>-1.439220e+15</td>\n",
       "      <td>0.000013</td>\n",
       "      <td>-1.000126</td>\n",
       "      <td>-43.101152</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>233.684960</td>\n",
       "      <td>116.865398</td>\n",
       "      <td>116.819562</td>\n",
       "      <td>0.045836</td>\n",
       "      <td>31.080516</td>\n",
       "      <td>84674.96180</td>\n",
       "      <td>-84706.04231</td>\n",
       "      <td>-84674.96180</td>\n",
       "      <td>-338575.5251</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.038891</td>\n",
       "      <td>-45.789753</td>\n",
       "      <td>954.210247</td>\n",
       "      <td>-22.237178</td>\n",
       "      <td>5.174343</td>\n",
       "      <td>-338575.5251</td>\n",
       "      <td>-1.437480e+15</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>-1.000139</td>\n",
       "      <td>-25.204567</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>233.682678</td>\n",
       "      <td>116.861333</td>\n",
       "      <td>116.821345</td>\n",
       "      <td>0.039988</td>\n",
       "      <td>31.089534</td>\n",
       "      <td>84674.09623</td>\n",
       "      <td>-84705.18576</td>\n",
       "      <td>-84674.09623</td>\n",
       "      <td>-338572.0268</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.039063</td>\n",
       "      <td>-45.900897</td>\n",
       "      <td>954.099103</td>\n",
       "      <td>-22.241197</td>\n",
       "      <td>5.177862</td>\n",
       "      <td>-338572.0268</td>\n",
       "      <td>-1.437440e+15</td>\n",
       "      <td>0.000012</td>\n",
       "      <td>-1.000143</td>\n",
       "      <td>-24.566389</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>233.684299</td>\n",
       "      <td>116.863892</td>\n",
       "      <td>116.820407</td>\n",
       "      <td>0.043485</td>\n",
       "      <td>31.077070</td>\n",
       "      <td>84675.02508</td>\n",
       "      <td>-84706.10215</td>\n",
       "      <td>-84675.02508</td>\n",
       "      <td>-338575.7920</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.038881</td>\n",
       "      <td>-45.757424</td>\n",
       "      <td>954.242576</td>\n",
       "      <td>-22.228654</td>\n",
       "      <td>5.144039</td>\n",
       "      <td>-338575.7920</td>\n",
       "      <td>-1.437490e+15</td>\n",
       "      <td>0.000007</td>\n",
       "      <td>-1.000131</td>\n",
       "      <td>-36.915158</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>233.680232</td>\n",
       "      <td>116.857671</td>\n",
       "      <td>116.822562</td>\n",
       "      <td>0.035109</td>\n",
       "      <td>31.094235</td>\n",
       "      <td>84673.49629</td>\n",
       "      <td>-84704.59052</td>\n",
       "      <td>-84673.49629</td>\n",
       "      <td>-338569.6082</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.039223</td>\n",
       "      <td>-45.862397</td>\n",
       "      <td>954.137603</td>\n",
       "      <td>-22.233754</td>\n",
       "      <td>5.152476</td>\n",
       "      <td>-338569.6082</td>\n",
       "      <td>-1.437410e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000139</td>\n",
       "      <td>-33.082327</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>233.660056</td>\n",
       "      <td>116.844493</td>\n",
       "      <td>116.815563</td>\n",
       "      <td>0.028930</td>\n",
       "      <td>31.188605</td>\n",
       "      <td>84660.85502</td>\n",
       "      <td>-84692.04362</td>\n",
       "      <td>-84660.85502</td>\n",
       "      <td>-338518.6656</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.042624</td>\n",
       "      <td>-46.148498</td>\n",
       "      <td>953.851502</td>\n",
       "      <td>-22.207181</td>\n",
       "      <td>5.394659</td>\n",
       "      <td>-338518.6656</td>\n",
       "      <td>-1.436760e+15</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>-1.000131</td>\n",
       "      <td>-41.905118</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>233.659796</td>\n",
       "      <td>116.841464</td>\n",
       "      <td>116.818331</td>\n",
       "      <td>0.023133</td>\n",
       "      <td>31.189976</td>\n",
       "      <td>84660.74416</td>\n",
       "      <td>-84691.93414</td>\n",
       "      <td>-84660.74416</td>\n",
       "      <td>-338518.2167</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.042799</td>\n",
       "      <td>-46.267557</td>\n",
       "      <td>953.732443</td>\n",
       "      <td>-22.210997</td>\n",
       "      <td>5.406434</td>\n",
       "      <td>-338518.2167</td>\n",
       "      <td>-1.436750e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000133</td>\n",
       "      <td>-44.138426</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>233.766039</td>\n",
       "      <td>116.890806</td>\n",
       "      <td>116.875233</td>\n",
       "      <td>0.015573</td>\n",
       "      <td>31.072035</td>\n",
       "      <td>84703.46444</td>\n",
       "      <td>-84734.53647</td>\n",
       "      <td>-84703.46444</td>\n",
       "      <td>-338689.5696</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.046127</td>\n",
       "      <td>-48.576920</td>\n",
       "      <td>951.423080</td>\n",
       "      <td>-22.194475</td>\n",
       "      <td>5.665077</td>\n",
       "      <td>-338689.5696</td>\n",
       "      <td>-1.438940e+15</td>\n",
       "      <td>0.000037</td>\n",
       "      <td>-1.000149</td>\n",
       "      <td>-31.623366</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>233.989373</td>\n",
       "      <td>117.103265</td>\n",
       "      <td>116.886109</td>\n",
       "      <td>0.217156</td>\n",
       "      <td>30.653131</td>\n",
       "      <td>84791.85094</td>\n",
       "      <td>-84822.50407</td>\n",
       "      <td>-84791.85094</td>\n",
       "      <td>-339044.7912</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.034460</td>\n",
       "      <td>-48.262548</td>\n",
       "      <td>951.737452</td>\n",
       "      <td>-22.158039</td>\n",
       "      <td>5.418296</td>\n",
       "      <td>-339044.7912</td>\n",
       "      <td>-1.443470e+15</td>\n",
       "      <td>0.000023</td>\n",
       "      <td>-1.000057</td>\n",
       "      <td>7.783007</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>233.700905</td>\n",
       "      <td>116.850453</td>\n",
       "      <td>116.850453</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>31.175118</td>\n",
       "      <td>84677.15287</td>\n",
       "      <td>-84708.32799</td>\n",
       "      <td>-84677.15287</td>\n",
       "      <td>-338583.9110</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.044221</td>\n",
       "      <td>-47.126240</td>\n",
       "      <td>952.873760</td>\n",
       "      <td>-22.192904</td>\n",
       "      <td>5.482490</td>\n",
       "      <td>-338583.9110</td>\n",
       "      <td>-1.437590e+15</td>\n",
       "      <td>0.000025</td>\n",
       "      <td>-1.000147</td>\n",
       "      <td>-32.572789</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>233.692043</td>\n",
       "      <td>116.845771</td>\n",
       "      <td>116.846272</td>\n",
       "      <td>-0.000501</td>\n",
       "      <td>31.172002</td>\n",
       "      <td>84674.59514</td>\n",
       "      <td>-84705.76714</td>\n",
       "      <td>-84674.59514</td>\n",
       "      <td>-338573.6925</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.043952</td>\n",
       "      <td>-48.335758</td>\n",
       "      <td>951.664242</td>\n",
       "      <td>-22.182560</td>\n",
       "      <td>5.494526</td>\n",
       "      <td>-338573.6925</td>\n",
       "      <td>-1.437460e+15</td>\n",
       "      <td>0.000013</td>\n",
       "      <td>-1.000133</td>\n",
       "      <td>-18.205948</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>233.683895</td>\n",
       "      <td>116.841947</td>\n",
       "      <td>116.841947</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>31.119784</td>\n",
       "      <td>84675.26735</td>\n",
       "      <td>-84706.38713</td>\n",
       "      <td>-84675.26735</td>\n",
       "      <td>-338576.5903</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.040763</td>\n",
       "      <td>-51.289154</td>\n",
       "      <td>948.710846</td>\n",
       "      <td>-22.239280</td>\n",
       "      <td>5.159963</td>\n",
       "      <td>-338576.5903</td>\n",
       "      <td>-1.437500e+15</td>\n",
       "      <td>0.000021</td>\n",
       "      <td>-1.000150</td>\n",
       "      <td>-9.395707</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>233.706175</td>\n",
       "      <td>116.861443</td>\n",
       "      <td>116.844732</td>\n",
       "      <td>0.016712</td>\n",
       "      <td>31.108954</td>\n",
       "      <td>84683.77806</td>\n",
       "      <td>-84714.88701</td>\n",
       "      <td>-84683.77806</td>\n",
       "      <td>-338610.6764</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.041578</td>\n",
       "      <td>-51.498951</td>\n",
       "      <td>948.501049</td>\n",
       "      <td>-22.217949</td>\n",
       "      <td>5.259791</td>\n",
       "      <td>-338610.6764</td>\n",
       "      <td>-1.437930e+15</td>\n",
       "      <td>0.000023</td>\n",
       "      <td>-1.000152</td>\n",
       "      <td>-37.604792</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>233.677162</td>\n",
       "      <td>116.843316</td>\n",
       "      <td>116.833846</td>\n",
       "      <td>0.009470</td>\n",
       "      <td>31.087350</td>\n",
       "      <td>84671.94161</td>\n",
       "      <td>-84703.02896</td>\n",
       "      <td>-84671.94161</td>\n",
       "      <td>-338563.4170</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.039112</td>\n",
       "      <td>-51.550110</td>\n",
       "      <td>948.449890</td>\n",
       "      <td>-22.230283</td>\n",
       "      <td>5.164084</td>\n",
       "      <td>-338563.4170</td>\n",
       "      <td>-1.437330e+15</td>\n",
       "      <td>0.000005</td>\n",
       "      <td>-1.000135</td>\n",
       "      <td>-36.990459</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>233.662028</td>\n",
       "      <td>116.832662</td>\n",
       "      <td>116.829366</td>\n",
       "      <td>0.003297</td>\n",
       "      <td>31.196917</td>\n",
       "      <td>84661.95208</td>\n",
       "      <td>-84693.14900</td>\n",
       "      <td>-84661.95208</td>\n",
       "      <td>-338523.0207</td>\n",
       "      <td>0.999998</td>\n",
       "      <td>...</td>\n",
       "      <td>0.043453</td>\n",
       "      <td>-51.871485</td>\n",
       "      <td>948.128515</td>\n",
       "      <td>-22.217914</td>\n",
       "      <td>5.436756</td>\n",
       "      <td>-338523.0207</td>\n",
       "      <td>-1.436820e+15</td>\n",
       "      <td>0.000014</td>\n",
       "      <td>-1.000136</td>\n",
       "      <td>-45.003762</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>233.679314</td>\n",
       "      <td>116.829211</td>\n",
       "      <td>116.850103</td>\n",
       "      <td>-0.020892</td>\n",
       "      <td>31.174119</td>\n",
       "      <td>84669.04975</td>\n",
       "      <td>-84700.22387</td>\n",
       "      <td>-84669.04975</td>\n",
       "      <td>-338551.5025</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.043784</td>\n",
       "      <td>-61.664581</td>\n",
       "      <td>938.335419</td>\n",
       "      <td>-22.211691</td>\n",
       "      <td>5.487538</td>\n",
       "      <td>-338551.5025</td>\n",
       "      <td>-1.437180e+15</td>\n",
       "      <td>0.000017</td>\n",
       "      <td>-1.000137</td>\n",
       "      <td>-51.507276</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>233.712655</td>\n",
       "      <td>116.856327</td>\n",
       "      <td>116.856327</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>31.119080</td>\n",
       "      <td>84686.06121</td>\n",
       "      <td>-84717.18029</td>\n",
       "      <td>-84686.06121</td>\n",
       "      <td>-338619.7685</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.041903</td>\n",
       "      <td>-60.913220</td>\n",
       "      <td>939.086780</td>\n",
       "      <td>-22.217026</td>\n",
       "      <td>5.189871</td>\n",
       "      <td>-338619.7685</td>\n",
       "      <td>-1.438050e+15</td>\n",
       "      <td>0.000035</td>\n",
       "      <td>-1.000164</td>\n",
       "      <td>-17.178714</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>233.693456</td>\n",
       "      <td>116.846494</td>\n",
       "      <td>116.846963</td>\n",
       "      <td>-0.000469</td>\n",
       "      <td>31.184952</td>\n",
       "      <td>84674.16226</td>\n",
       "      <td>-84705.34721</td>\n",
       "      <td>-84674.16226</td>\n",
       "      <td>-338571.9092</td>\n",
       "      <td>0.999999</td>\n",
       "      <td>...</td>\n",
       "      <td>0.044217</td>\n",
       "      <td>-47.673505</td>\n",
       "      <td>952.326495</td>\n",
       "      <td>-22.185596</td>\n",
       "      <td>5.460037</td>\n",
       "      <td>-338571.9092</td>\n",
       "      <td>-1.437440e+15</td>\n",
       "      <td>0.000024</td>\n",
       "      <td>-1.000145</td>\n",
       "      <td>-30.992091</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>28 rows × 29 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            D1          D2          D3        D4         D5           D6  \\\n",
       "1   233.762715  116.898844  116.863871  0.034973  30.917600  84706.25146   \n",
       "2   233.754944  116.890477  116.864467  0.026010  30.913406  84703.02826   \n",
       "3   233.756958  116.904940  116.852018  0.052922  30.843638  84702.00879   \n",
       "4   233.757343  116.900775  116.856568  0.044207  30.849927  84703.15045   \n",
       "5   233.771601  116.910762  116.860839  0.049922  30.854660  84709.04221   \n",
       "6   233.769472  116.904063  116.865409  0.038654  30.814714  84707.62008   \n",
       "7   233.763229  116.902341  116.860888  0.041454  30.792690  84705.32736   \n",
       "8   233.802668  116.925918  116.876749  0.049169  30.649775  84719.30084   \n",
       "9   233.733487  116.876740  116.856747  0.019993  30.971601  84690.28825   \n",
       "10  233.693647  116.866452  116.827195  0.039257  31.075569  84678.89210   \n",
       "11  233.770927  116.932869  116.838058  0.094811  30.917504  84708.85367   \n",
       "12  233.684960  116.865398  116.819562  0.045836  31.080516  84674.96180   \n",
       "13  233.682678  116.861333  116.821345  0.039988  31.089534  84674.09623   \n",
       "14  233.684299  116.863892  116.820407  0.043485  31.077070  84675.02508   \n",
       "15  233.680232  116.857671  116.822562  0.035109  31.094235  84673.49629   \n",
       "16  233.660056  116.844493  116.815563  0.028930  31.188605  84660.85502   \n",
       "17  233.659796  116.841464  116.818331  0.023133  31.189976  84660.74416   \n",
       "18  233.766039  116.890806  116.875233  0.015573  31.072035  84703.46444   \n",
       "19  233.989373  117.103265  116.886109  0.217156  30.653131  84791.85094   \n",
       "20  233.700905  116.850453  116.850453  0.000000  31.175118  84677.15287   \n",
       "21  233.692043  116.845771  116.846272 -0.000501  31.172002  84674.59514   \n",
       "22  233.683895  116.841947  116.841947  0.000000  31.119784  84675.26735   \n",
       "23  233.706175  116.861443  116.844732  0.016712  31.108954  84683.77806   \n",
       "24  233.677162  116.843316  116.833846  0.009470  31.087350  84671.94161   \n",
       "25  233.662028  116.832662  116.829366  0.003297  31.196917  84661.95208   \n",
       "26  233.679314  116.829211  116.850103 -0.020892  31.174119  84669.04975   \n",
       "27  233.712655  116.856327  116.856327  0.000000  31.119080  84686.06121   \n",
       "28  233.693456  116.846494  116.846963 -0.000469  31.184952  84674.16226   \n",
       "\n",
       "             D7           D8           D9       D10  ...       D20        D21  \\\n",
       "1  -84737.16906 -84706.25146 -338701.3354  0.999999  ...  0.040350 -45.693931   \n",
       "2  -84733.94167 -84703.02826 -338688.4594  0.999999  ...  0.040175 -45.815157   \n",
       "3  -84732.85243 -84702.00879 -338684.6606  0.999999  ...  0.037591 -45.967968   \n",
       "4  -84734.00037 -84703.15045 -338689.2021  0.999999  ...  0.037834 -45.913815   \n",
       "5  -84739.89687 -84709.04221 -338712.7502  0.999999  ...  0.038212 -45.820667   \n",
       "6  -84738.43479 -84707.62008 -338707.2214  0.999999  ...  0.037772 -45.975482   \n",
       "7  -84736.12005 -84705.32736 -338698.1387  0.999999  ...  0.037038 -45.931454   \n",
       "8  -84749.95062 -84719.30084 -338754.6043  0.999999  ...  0.033831 -46.325084   \n",
       "9  -84721.25985 -84690.28825 -338637.2666  0.999999  ...  0.041815 -46.221557   \n",
       "10 -84709.96767 -84678.89210 -338591.2661  0.999999  ...  0.039874 -45.725674   \n",
       "11 -84739.77117 -84708.85367 -338711.7447  0.999999  ...  0.040042 -45.823397   \n",
       "12 -84706.04231 -84674.96180 -338575.5251  0.999999  ...  0.038891 -45.789753   \n",
       "13 -84705.18576 -84674.09623 -338572.0268  0.999999  ...  0.039063 -45.900897   \n",
       "14 -84706.10215 -84675.02508 -338575.7920  0.999999  ...  0.038881 -45.757424   \n",
       "15 -84704.59052 -84673.49629 -338569.6082  0.999999  ...  0.039223 -45.862397   \n",
       "16 -84692.04362 -84660.85502 -338518.6656  0.999999  ...  0.042624 -46.148498   \n",
       "17 -84691.93414 -84660.74416 -338518.2167  0.999999  ...  0.042799 -46.267557   \n",
       "18 -84734.53647 -84703.46444 -338689.5696  0.999999  ...  0.046127 -48.576920   \n",
       "19 -84822.50407 -84791.85094 -339044.7912  0.999999  ...  0.034460 -48.262548   \n",
       "20 -84708.32799 -84677.15287 -338583.9110  0.999999  ...  0.044221 -47.126240   \n",
       "21 -84705.76714 -84674.59514 -338573.6925  0.999999  ...  0.043952 -48.335758   \n",
       "22 -84706.38713 -84675.26735 -338576.5903  0.999999  ...  0.040763 -51.289154   \n",
       "23 -84714.88701 -84683.77806 -338610.6764  0.999999  ...  0.041578 -51.498951   \n",
       "24 -84703.02896 -84671.94161 -338563.4170  0.999999  ...  0.039112 -51.550110   \n",
       "25 -84693.14900 -84661.95208 -338523.0207  0.999998  ...  0.043453 -51.871485   \n",
       "26 -84700.22387 -84669.04975 -338551.5025  0.999999  ...  0.043784 -61.664581   \n",
       "27 -84717.18029 -84686.06121 -338619.7685  0.999999  ...  0.041903 -60.913220   \n",
       "28 -84705.34721 -84674.16226 -338571.9092  0.999999  ...  0.044217 -47.673505   \n",
       "\n",
       "           D22        D23       D24          D25           D26       D27  \\\n",
       "1   954.306069 -22.218856  5.469140 -338701.3354 -1.439090e+15  0.000013   \n",
       "2   954.184843 -22.224564  5.432392 -338688.4594 -1.438920e+15  0.000009   \n",
       "3   954.032032 -22.224253  5.414609 -338684.6606 -1.438870e+15  0.000006   \n",
       "4   954.086185 -22.229211  5.374452 -338689.2021 -1.438930e+15  0.000003   \n",
       "5   954.179333 -22.221312  5.379066 -338712.7502 -1.439230e+15  0.000007   \n",
       "6   954.024518 -22.211994  5.368161 -338707.2214 -1.439160e+15  0.000006   \n",
       "7   954.068546 -22.218519  5.330169 -338698.1387 -1.439050e+15  0.000011   \n",
       "8   953.674916 -22.192825  5.369025 -338754.6043 -1.439770e+15  0.000033   \n",
       "9   953.778443 -22.197462  5.580163 -338637.2666 -1.438270e+15  0.000003   \n",
       "10  954.274326 -22.227629  5.248522 -338591.2661 -1.437680e+15  0.000011   \n",
       "11  954.176603 -22.222212  5.454897 -338711.7447 -1.439220e+15  0.000013   \n",
       "12  954.210247 -22.237178  5.174343 -338575.5251 -1.437480e+15  0.000010   \n",
       "13  954.099103 -22.241197  5.177862 -338572.0268 -1.437440e+15  0.000012   \n",
       "14  954.242576 -22.228654  5.144039 -338575.7920 -1.437490e+15  0.000007   \n",
       "15  954.137603 -22.233754  5.152476 -338569.6082 -1.437410e+15  0.000011   \n",
       "16  953.851502 -22.207181  5.394659 -338518.6656 -1.436760e+15  0.000010   \n",
       "17  953.732443 -22.210997  5.406434 -338518.2167 -1.436750e+15  0.000011   \n",
       "18  951.423080 -22.194475  5.665077 -338689.5696 -1.438940e+15  0.000037   \n",
       "19  951.737452 -22.158039  5.418296 -339044.7912 -1.443470e+15  0.000023   \n",
       "20  952.873760 -22.192904  5.482490 -338583.9110 -1.437590e+15  0.000025   \n",
       "21  951.664242 -22.182560  5.494526 -338573.6925 -1.437460e+15  0.000013   \n",
       "22  948.710846 -22.239280  5.159963 -338576.5903 -1.437500e+15  0.000021   \n",
       "23  948.501049 -22.217949  5.259791 -338610.6764 -1.437930e+15  0.000023   \n",
       "24  948.449890 -22.230283  5.164084 -338563.4170 -1.437330e+15  0.000005   \n",
       "25  948.128515 -22.217914  5.436756 -338523.0207 -1.436820e+15  0.000014   \n",
       "26  938.335419 -22.211691  5.487538 -338551.5025 -1.437180e+15  0.000017   \n",
       "27  939.086780 -22.217026  5.189871 -338619.7685 -1.438050e+15  0.000035   \n",
       "28  952.326495 -22.185596  5.460037 -338571.9092 -1.437440e+15  0.000024   \n",
       "\n",
       "         D28          y  \n",
       "1  -1.000127 -42.935489  \n",
       "2  -1.000118 -33.770706  \n",
       "3  -1.000107 -32.560866  \n",
       "4  -1.000104 -23.811494  \n",
       "5  -1.000113 -27.856424  \n",
       "6  -1.000095 -32.063878  \n",
       "7  -1.000094 -29.522463  \n",
       "8  -1.000085 -33.272463  \n",
       "9  -1.000105 -41.417543  \n",
       "10 -1.000139 -42.223893  \n",
       "11 -1.000126 -43.101152  \n",
       "12 -1.000139 -25.204567  \n",
       "13 -1.000143 -24.566389  \n",
       "14 -1.000131 -36.915158  \n",
       "15 -1.000139 -33.082327  \n",
       "16 -1.000131 -41.905118  \n",
       "17 -1.000133 -44.138426  \n",
       "18 -1.000149 -31.623366  \n",
       "19 -1.000057   7.783007  \n",
       "20 -1.000147 -32.572789  \n",
       "21 -1.000133 -18.205948  \n",
       "22 -1.000150  -9.395707  \n",
       "23 -1.000152 -37.604792  \n",
       "24 -1.000135 -36.990459  \n",
       "25 -1.000136 -45.003762  \n",
       "26 -1.000137 -51.507276  \n",
       "27 -1.000164 -17.178714  \n",
       "28 -1.000145 -30.992091  \n",
       "\n",
       "[28 rows x 29 columns]"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#1. 输入库\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "\n",
    "#2. 打开d带中心数据\n",
    "with open('ML_DATA.csv', encoding = 'UTF-8') as f:\n",
    "    df = pd.read_csv(f, index_col = 0)  #index_col=0 告诉第一列是索引\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b17eeeb3",
   "metadata": {},
   "source": [
    "剔除描述符"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "109c51a9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>D3</th>\n",
       "      <th>D4</th>\n",
       "      <th>D11</th>\n",
       "      <th>D15</th>\n",
       "      <th>D16</th>\n",
       "      <th>D17</th>\n",
       "      <th>D19</th>\n",
       "      <th>D20</th>\n",
       "      <th>D22</th>\n",
       "      <th>D23</th>\n",
       "      <th>D24</th>\n",
       "      <th>D26</th>\n",
       "      <th>D27</th>\n",
       "      <th>D28</th>\n",
       "      <th>y</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>116.863871</td>\n",
       "      <td>0.034973</td>\n",
       "      <td>0.998788</td>\n",
       "      <td>-0.000199</td>\n",
       "      <td>3867.231391</td>\n",
       "      <td>-1.986133</td>\n",
       "      <td>0.000199</td>\n",
       "      <td>0.040350</td>\n",
       "      <td>954.306069</td>\n",
       "      <td>-22.218856</td>\n",
       "      <td>5.469140</td>\n",
       "      <td>-1.439090e+15</td>\n",
       "      <td>0.000013</td>\n",
       "      <td>-1.000127</td>\n",
       "      <td>-42.935489</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>116.864467</td>\n",
       "      <td>0.026010</td>\n",
       "      <td>0.998788</td>\n",
       "      <td>-0.000267</td>\n",
       "      <td>3903.249000</td>\n",
       "      <td>10.447788</td>\n",
       "      <td>0.000270</td>\n",
       "      <td>0.040175</td>\n",
       "      <td>954.184843</td>\n",
       "      <td>-22.224564</td>\n",
       "      <td>5.432392</td>\n",
       "      <td>-1.438920e+15</td>\n",
       "      <td>0.000009</td>\n",
       "      <td>-1.000118</td>\n",
       "      <td>-33.770706</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>116.852018</td>\n",
       "      <td>0.052922</td>\n",
       "      <td>0.998790</td>\n",
       "      <td>-0.000264</td>\n",
       "      <td>3953.777336</td>\n",
       "      <td>10.517678</td>\n",
       "      <td>0.000356</td>\n",
       "      <td>0.037591</td>\n",
       "      <td>954.032032</td>\n",
       "      <td>-22.224253</td>\n",
       "      <td>5.414609</td>\n",
       "      <td>-1.438870e+15</td>\n",
       "      <td>0.000006</td>\n",
       "      <td>-1.000107</td>\n",
       "      <td>-32.560866</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>116.856568</td>\n",
       "      <td>0.044207</td>\n",
       "      <td>0.998790</td>\n",
       "      <td>-0.000419</td>\n",
       "      <td>3940.268100</td>\n",
       "      <td>10.407253</td>\n",
       "      <td>0.000349</td>\n",
       "      <td>0.037834</td>\n",
       "      <td>954.086185</td>\n",
       "      <td>-22.229211</td>\n",
       "      <td>5.374452</td>\n",
       "      <td>-1.438930e+15</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>-1.000104</td>\n",
       "      <td>-23.811494</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>116.860839</td>\n",
       "      <td>0.049922</td>\n",
       "      <td>0.998790</td>\n",
       "      <td>-0.000091</td>\n",
       "      <td>3913.659032</td>\n",
       "      <td>10.472823</td>\n",
       "      <td>0.000302</td>\n",
       "      <td>0.038212</td>\n",
       "      <td>954.179333</td>\n",
       "      <td>-22.221312</td>\n",
       "      <td>5.379066</td>\n",
       "      <td>-1.439230e+15</td>\n",
       "      <td>0.000007</td>\n",
       "      <td>-1.000113</td>\n",
       "      <td>-27.856424</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>116.865409</td>\n",
       "      <td>0.038654</td>\n",
       "      <td>0.998792</td>\n",
       "      <td>-0.000154</td>\n",
       "      <td>3962.903026</td>\n",
       "      <td>10.520122</td>\n",
       "      <td>0.000407</td>\n",
       "      <td>0.037772</td>\n",
       "      <td>954.024518</td>\n",
       "      <td>-22.211994</td>\n",
       "      <td>5.368161</td>\n",
       "      <td>-1.439160e+15</td>\n",
       "      <td>0.000006</td>\n",
       "      <td>-1.000095</td>\n",
       "      <td>-32.063878</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>116.860888</td>\n",
       "      <td>0.041454</td>\n",
       "      <td>0.998792</td>\n",
       "      <td>-0.000054</td>\n",
       "      <td>3962.591314</td>\n",
       "      <td>10.518382</td>\n",
       "      <td>0.000403</td>\n",
       "      <td>0.037038</td>\n",
       "      <td>954.068546</td>\n",
       "      <td>-22.218519</td>\n",
       "      <td>5.330169</td>\n",
       "      <td>-1.439050e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000094</td>\n",
       "      <td>-29.522463</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>116.876749</td>\n",
       "      <td>0.049169</td>\n",
       "      <td>0.998798</td>\n",
       "      <td>-0.000030</td>\n",
       "      <td>4066.159834</td>\n",
       "      <td>10.525493</td>\n",
       "      <td>0.000789</td>\n",
       "      <td>0.033831</td>\n",
       "      <td>953.674916</td>\n",
       "      <td>-22.192825</td>\n",
       "      <td>5.369025</td>\n",
       "      <td>-1.439770e+15</td>\n",
       "      <td>0.000033</td>\n",
       "      <td>-1.000085</td>\n",
       "      <td>-33.272463</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>116.856747</td>\n",
       "      <td>0.019993</td>\n",
       "      <td>0.998785</td>\n",
       "      <td>-0.000203</td>\n",
       "      <td>3992.408767</td>\n",
       "      <td>10.522068</td>\n",
       "      <td>0.000376</td>\n",
       "      <td>0.041815</td>\n",
       "      <td>953.778443</td>\n",
       "      <td>-22.197462</td>\n",
       "      <td>5.580163</td>\n",
       "      <td>-1.438270e+15</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>-1.000105</td>\n",
       "      <td>-41.417543</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>116.827195</td>\n",
       "      <td>0.039257</td>\n",
       "      <td>0.998781</td>\n",
       "      <td>-0.000152</td>\n",
       "      <td>3776.169456</td>\n",
       "      <td>-0.704697</td>\n",
       "      <td>0.000438</td>\n",
       "      <td>0.039874</td>\n",
       "      <td>954.274326</td>\n",
       "      <td>-22.227629</td>\n",
       "      <td>5.248522</td>\n",
       "      <td>-1.437680e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000139</td>\n",
       "      <td>-42.223893</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>116.838058</td>\n",
       "      <td>0.094811</td>\n",
       "      <td>0.998788</td>\n",
       "      <td>-0.000190</td>\n",
       "      <td>3872.567480</td>\n",
       "      <td>8.163972</td>\n",
       "      <td>0.000217</td>\n",
       "      <td>0.040042</td>\n",
       "      <td>954.176603</td>\n",
       "      <td>-22.222212</td>\n",
       "      <td>5.454897</td>\n",
       "      <td>-1.439220e+15</td>\n",
       "      <td>0.000013</td>\n",
       "      <td>-1.000126</td>\n",
       "      <td>-43.101152</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>116.819562</td>\n",
       "      <td>0.045836</td>\n",
       "      <td>0.998781</td>\n",
       "      <td>-0.000353</td>\n",
       "      <td>3793.102605</td>\n",
       "      <td>-0.408854</td>\n",
       "      <td>0.000613</td>\n",
       "      <td>0.038891</td>\n",
       "      <td>954.210247</td>\n",
       "      <td>-22.237178</td>\n",
       "      <td>5.174343</td>\n",
       "      <td>-1.437480e+15</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>-1.000139</td>\n",
       "      <td>-25.204567</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>116.821345</td>\n",
       "      <td>0.039988</td>\n",
       "      <td>0.998780</td>\n",
       "      <td>-0.000445</td>\n",
       "      <td>3805.251289</td>\n",
       "      <td>-0.358081</td>\n",
       "      <td>0.000580</td>\n",
       "      <td>0.039063</td>\n",
       "      <td>954.099103</td>\n",
       "      <td>-22.241197</td>\n",
       "      <td>5.177862</td>\n",
       "      <td>-1.437440e+15</td>\n",
       "      <td>0.000012</td>\n",
       "      <td>-1.000143</td>\n",
       "      <td>-24.566389</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>116.820407</td>\n",
       "      <td>0.043485</td>\n",
       "      <td>0.998781</td>\n",
       "      <td>-0.000106</td>\n",
       "      <td>3783.130054</td>\n",
       "      <td>10.052937</td>\n",
       "      <td>0.000689</td>\n",
       "      <td>0.038881</td>\n",
       "      <td>954.242576</td>\n",
       "      <td>-22.228654</td>\n",
       "      <td>5.144039</td>\n",
       "      <td>-1.437490e+15</td>\n",
       "      <td>0.000007</td>\n",
       "      <td>-1.000131</td>\n",
       "      <td>-36.915158</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>116.822562</td>\n",
       "      <td>0.035109</td>\n",
       "      <td>0.998780</td>\n",
       "      <td>-0.000198</td>\n",
       "      <td>3792.298165</td>\n",
       "      <td>-2.174566</td>\n",
       "      <td>0.000623</td>\n",
       "      <td>0.039223</td>\n",
       "      <td>954.137603</td>\n",
       "      <td>-22.233754</td>\n",
       "      <td>5.152476</td>\n",
       "      <td>-1.437410e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000139</td>\n",
       "      <td>-33.082327</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>116.815563</td>\n",
       "      <td>0.028930</td>\n",
       "      <td>0.998776</td>\n",
       "      <td>-0.000091</td>\n",
       "      <td>3864.253177</td>\n",
       "      <td>10.510031</td>\n",
       "      <td>0.000679</td>\n",
       "      <td>0.042624</td>\n",
       "      <td>953.851502</td>\n",
       "      <td>-22.207181</td>\n",
       "      <td>5.394659</td>\n",
       "      <td>-1.436760e+15</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>-1.000131</td>\n",
       "      <td>-41.905118</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>116.818331</td>\n",
       "      <td>0.023133</td>\n",
       "      <td>0.998776</td>\n",
       "      <td>-0.000144</td>\n",
       "      <td>3882.868331</td>\n",
       "      <td>10.492924</td>\n",
       "      <td>0.000638</td>\n",
       "      <td>0.042799</td>\n",
       "      <td>953.732443</td>\n",
       "      <td>-22.210997</td>\n",
       "      <td>5.406434</td>\n",
       "      <td>-1.436750e+15</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>-1.000133</td>\n",
       "      <td>-44.138426</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>116.875233</td>\n",
       "      <td>0.015573</td>\n",
       "      <td>0.998782</td>\n",
       "      <td>-0.000191</td>\n",
       "      <td>4244.715508</td>\n",
       "      <td>10.154544</td>\n",
       "      <td>0.000155</td>\n",
       "      <td>0.046127</td>\n",
       "      <td>951.423080</td>\n",
       "      <td>-22.194475</td>\n",
       "      <td>5.665077</td>\n",
       "      <td>-1.438940e+15</td>\n",
       "      <td>0.000037</td>\n",
       "      <td>-1.000149</td>\n",
       "      <td>-31.623366</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>116.886109</td>\n",
       "      <td>0.217156</td>\n",
       "      <td>0.998800</td>\n",
       "      <td>-0.000382</td>\n",
       "      <td>3997.222613</td>\n",
       "      <td>8.089390</td>\n",
       "      <td>0.000705</td>\n",
       "      <td>0.034460</td>\n",
       "      <td>951.737452</td>\n",
       "      <td>-22.158039</td>\n",
       "      <td>5.418296</td>\n",
       "      <td>-1.443470e+15</td>\n",
       "      <td>0.000023</td>\n",
       "      <td>-1.000057</td>\n",
       "      <td>7.783007</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>116.850453</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.998777</td>\n",
       "      <td>-0.000270</td>\n",
       "      <td>3991.267727</td>\n",
       "      <td>10.131371</td>\n",
       "      <td>0.000687</td>\n",
       "      <td>0.044221</td>\n",
       "      <td>952.873760</td>\n",
       "      <td>-22.192904</td>\n",
       "      <td>5.482490</td>\n",
       "      <td>-1.437590e+15</td>\n",
       "      <td>0.000025</td>\n",
       "      <td>-1.000147</td>\n",
       "      <td>-32.572789</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>116.846272</td>\n",
       "      <td>-0.000501</td>\n",
       "      <td>0.998777</td>\n",
       "      <td>-0.000461</td>\n",
       "      <td>4130.662485</td>\n",
       "      <td>0.015884</td>\n",
       "      <td>0.000659</td>\n",
       "      <td>0.043952</td>\n",
       "      <td>951.664242</td>\n",
       "      <td>-22.182560</td>\n",
       "      <td>5.494526</td>\n",
       "      <td>-1.437460e+15</td>\n",
       "      <td>0.000013</td>\n",
       "      <td>-1.000133</td>\n",
       "      <td>-18.205948</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>116.841947</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.998779</td>\n",
       "      <td>-0.000726</td>\n",
       "      <td>2954.505165</td>\n",
       "      <td>7.459893</td>\n",
       "      <td>0.000657</td>\n",
       "      <td>0.040763</td>\n",
       "      <td>948.710846</td>\n",
       "      <td>-22.239280</td>\n",
       "      <td>5.159963</td>\n",
       "      <td>-1.437500e+15</td>\n",
       "      <td>0.000021</td>\n",
       "      <td>-1.000150</td>\n",
       "      <td>-9.395707</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>116.844732</td>\n",
       "      <td>0.016712</td>\n",
       "      <td>0.998780</td>\n",
       "      <td>-0.001043</td>\n",
       "      <td>2979.722916</td>\n",
       "      <td>10.348614</td>\n",
       "      <td>0.000511</td>\n",
       "      <td>0.041578</td>\n",
       "      <td>948.501049</td>\n",
       "      <td>-22.217949</td>\n",
       "      <td>5.259791</td>\n",
       "      <td>-1.437930e+15</td>\n",
       "      <td>0.000023</td>\n",
       "      <td>-1.000152</td>\n",
       "      <td>-37.604792</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>116.833846</td>\n",
       "      <td>0.009470</td>\n",
       "      <td>0.998780</td>\n",
       "      <td>-0.000574</td>\n",
       "      <td>2960.792158</td>\n",
       "      <td>-0.471472</td>\n",
       "      <td>0.001008</td>\n",
       "      <td>0.039112</td>\n",
       "      <td>948.449890</td>\n",
       "      <td>-22.230283</td>\n",
       "      <td>5.164084</td>\n",
       "      <td>-1.437330e+15</td>\n",
       "      <td>0.000005</td>\n",
       "      <td>-1.000135</td>\n",
       "      <td>-36.990459</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>116.829366</td>\n",
       "      <td>0.003297</td>\n",
       "      <td>0.998776</td>\n",
       "      <td>-0.000056</td>\n",
       "      <td>3007.712242</td>\n",
       "      <td>-1.949231</td>\n",
       "      <td>0.000718</td>\n",
       "      <td>0.043453</td>\n",
       "      <td>948.128515</td>\n",
       "      <td>-22.217914</td>\n",
       "      <td>5.436756</td>\n",
       "      <td>-1.436820e+15</td>\n",
       "      <td>0.000014</td>\n",
       "      <td>-1.000136</td>\n",
       "      <td>-45.003762</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>116.850103</td>\n",
       "      <td>-0.020892</td>\n",
       "      <td>0.998777</td>\n",
       "      <td>-0.000080</td>\n",
       "      <td>2438.494102</td>\n",
       "      <td>10.524549</td>\n",
       "      <td>0.000541</td>\n",
       "      <td>0.043784</td>\n",
       "      <td>938.335419</td>\n",
       "      <td>-22.211691</td>\n",
       "      <td>5.487538</td>\n",
       "      <td>-1.437180e+15</td>\n",
       "      <td>0.000017</td>\n",
       "      <td>-1.000137</td>\n",
       "      <td>-51.507276</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>116.856327</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.998779</td>\n",
       "      <td>-0.000013</td>\n",
       "      <td>2375.735107</td>\n",
       "      <td>10.507976</td>\n",
       "      <td>0.000531</td>\n",
       "      <td>0.041903</td>\n",
       "      <td>939.086780</td>\n",
       "      <td>-22.217026</td>\n",
       "      <td>5.189871</td>\n",
       "      <td>-1.438050e+15</td>\n",
       "      <td>0.000035</td>\n",
       "      <td>-1.000164</td>\n",
       "      <td>-17.178714</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>116.846963</td>\n",
       "      <td>-0.000469</td>\n",
       "      <td>0.998777</td>\n",
       "      <td>-0.000039</td>\n",
       "      <td>4016.004449</td>\n",
       "      <td>0.012835</td>\n",
       "      <td>0.000790</td>\n",
       "      <td>0.044217</td>\n",
       "      <td>952.326495</td>\n",
       "      <td>-22.185596</td>\n",
       "      <td>5.460037</td>\n",
       "      <td>-1.437440e+15</td>\n",
       "      <td>0.000024</td>\n",
       "      <td>-1.000145</td>\n",
       "      <td>-30.992091</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            D3        D4       D11       D15          D16        D17  \\\n",
       "1   116.863871  0.034973  0.998788 -0.000199  3867.231391  -1.986133   \n",
       "2   116.864467  0.026010  0.998788 -0.000267  3903.249000  10.447788   \n",
       "3   116.852018  0.052922  0.998790 -0.000264  3953.777336  10.517678   \n",
       "4   116.856568  0.044207  0.998790 -0.000419  3940.268100  10.407253   \n",
       "5   116.860839  0.049922  0.998790 -0.000091  3913.659032  10.472823   \n",
       "6   116.865409  0.038654  0.998792 -0.000154  3962.903026  10.520122   \n",
       "7   116.860888  0.041454  0.998792 -0.000054  3962.591314  10.518382   \n",
       "8   116.876749  0.049169  0.998798 -0.000030  4066.159834  10.525493   \n",
       "9   116.856747  0.019993  0.998785 -0.000203  3992.408767  10.522068   \n",
       "10  116.827195  0.039257  0.998781 -0.000152  3776.169456  -0.704697   \n",
       "11  116.838058  0.094811  0.998788 -0.000190  3872.567480   8.163972   \n",
       "12  116.819562  0.045836  0.998781 -0.000353  3793.102605  -0.408854   \n",
       "13  116.821345  0.039988  0.998780 -0.000445  3805.251289  -0.358081   \n",
       "14  116.820407  0.043485  0.998781 -0.000106  3783.130054  10.052937   \n",
       "15  116.822562  0.035109  0.998780 -0.000198  3792.298165  -2.174566   \n",
       "16  116.815563  0.028930  0.998776 -0.000091  3864.253177  10.510031   \n",
       "17  116.818331  0.023133  0.998776 -0.000144  3882.868331  10.492924   \n",
       "18  116.875233  0.015573  0.998782 -0.000191  4244.715508  10.154544   \n",
       "19  116.886109  0.217156  0.998800 -0.000382  3997.222613   8.089390   \n",
       "20  116.850453  0.000000  0.998777 -0.000270  3991.267727  10.131371   \n",
       "21  116.846272 -0.000501  0.998777 -0.000461  4130.662485   0.015884   \n",
       "22  116.841947  0.000000  0.998779 -0.000726  2954.505165   7.459893   \n",
       "23  116.844732  0.016712  0.998780 -0.001043  2979.722916  10.348614   \n",
       "24  116.833846  0.009470  0.998780 -0.000574  2960.792158  -0.471472   \n",
       "25  116.829366  0.003297  0.998776 -0.000056  3007.712242  -1.949231   \n",
       "26  116.850103 -0.020892  0.998777 -0.000080  2438.494102  10.524549   \n",
       "27  116.856327  0.000000  0.998779 -0.000013  2375.735107  10.507976   \n",
       "28  116.846963 -0.000469  0.998777 -0.000039  4016.004449   0.012835   \n",
       "\n",
       "         D19       D20         D22        D23       D24           D26  \\\n",
       "1   0.000199  0.040350  954.306069 -22.218856  5.469140 -1.439090e+15   \n",
       "2   0.000270  0.040175  954.184843 -22.224564  5.432392 -1.438920e+15   \n",
       "3   0.000356  0.037591  954.032032 -22.224253  5.414609 -1.438870e+15   \n",
       "4   0.000349  0.037834  954.086185 -22.229211  5.374452 -1.438930e+15   \n",
       "5   0.000302  0.038212  954.179333 -22.221312  5.379066 -1.439230e+15   \n",
       "6   0.000407  0.037772  954.024518 -22.211994  5.368161 -1.439160e+15   \n",
       "7   0.000403  0.037038  954.068546 -22.218519  5.330169 -1.439050e+15   \n",
       "8   0.000789  0.033831  953.674916 -22.192825  5.369025 -1.439770e+15   \n",
       "9   0.000376  0.041815  953.778443 -22.197462  5.580163 -1.438270e+15   \n",
       "10  0.000438  0.039874  954.274326 -22.227629  5.248522 -1.437680e+15   \n",
       "11  0.000217  0.040042  954.176603 -22.222212  5.454897 -1.439220e+15   \n",
       "12  0.000613  0.038891  954.210247 -22.237178  5.174343 -1.437480e+15   \n",
       "13  0.000580  0.039063  954.099103 -22.241197  5.177862 -1.437440e+15   \n",
       "14  0.000689  0.038881  954.242576 -22.228654  5.144039 -1.437490e+15   \n",
       "15  0.000623  0.039223  954.137603 -22.233754  5.152476 -1.437410e+15   \n",
       "16  0.000679  0.042624  953.851502 -22.207181  5.394659 -1.436760e+15   \n",
       "17  0.000638  0.042799  953.732443 -22.210997  5.406434 -1.436750e+15   \n",
       "18  0.000155  0.046127  951.423080 -22.194475  5.665077 -1.438940e+15   \n",
       "19  0.000705  0.034460  951.737452 -22.158039  5.418296 -1.443470e+15   \n",
       "20  0.000687  0.044221  952.873760 -22.192904  5.482490 -1.437590e+15   \n",
       "21  0.000659  0.043952  951.664242 -22.182560  5.494526 -1.437460e+15   \n",
       "22  0.000657  0.040763  948.710846 -22.239280  5.159963 -1.437500e+15   \n",
       "23  0.000511  0.041578  948.501049 -22.217949  5.259791 -1.437930e+15   \n",
       "24  0.001008  0.039112  948.449890 -22.230283  5.164084 -1.437330e+15   \n",
       "25  0.000718  0.043453  948.128515 -22.217914  5.436756 -1.436820e+15   \n",
       "26  0.000541  0.043784  938.335419 -22.211691  5.487538 -1.437180e+15   \n",
       "27  0.000531  0.041903  939.086780 -22.217026  5.189871 -1.438050e+15   \n",
       "28  0.000790  0.044217  952.326495 -22.185596  5.460037 -1.437440e+15   \n",
       "\n",
       "         D27       D28          y  \n",
       "1   0.000013 -1.000127 -42.935489  \n",
       "2   0.000009 -1.000118 -33.770706  \n",
       "3   0.000006 -1.000107 -32.560866  \n",
       "4   0.000003 -1.000104 -23.811494  \n",
       "5   0.000007 -1.000113 -27.856424  \n",
       "6   0.000006 -1.000095 -32.063878  \n",
       "7   0.000011 -1.000094 -29.522463  \n",
       "8   0.000033 -1.000085 -33.272463  \n",
       "9   0.000003 -1.000105 -41.417543  \n",
       "10  0.000011 -1.000139 -42.223893  \n",
       "11  0.000013 -1.000126 -43.101152  \n",
       "12  0.000010 -1.000139 -25.204567  \n",
       "13  0.000012 -1.000143 -24.566389  \n",
       "14  0.000007 -1.000131 -36.915158  \n",
       "15  0.000011 -1.000139 -33.082327  \n",
       "16  0.000010 -1.000131 -41.905118  \n",
       "17  0.000011 -1.000133 -44.138426  \n",
       "18  0.000037 -1.000149 -31.623366  \n",
       "19  0.000023 -1.000057   7.783007  \n",
       "20  0.000025 -1.000147 -32.572789  \n",
       "21  0.000013 -1.000133 -18.205948  \n",
       "22  0.000021 -1.000150  -9.395707  \n",
       "23  0.000023 -1.000152 -37.604792  \n",
       "24  0.000005 -1.000135 -36.990459  \n",
       "25  0.000014 -1.000136 -45.003762  \n",
       "26  0.000017 -1.000137 -51.507276  \n",
       "27  0.000035 -1.000164 -17.178714  \n",
       "28  0.000024 -1.000145 -30.992091  "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "excluded = ['D1','D2','D5','D6','D7','D8','D9','D10','D12','D13','D14','D18','D21','D25']\n",
    "df_new = df.drop(excluded, axis = 1)\n",
    "df_new"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1e15c948",
   "metadata": {},
   "source": [
    "设置自变量和因变量"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "2724489d",
   "metadata": {},
   "outputs": [],
   "source": [
    "y = df_new['y'].values\n",
    "X_1 = df_new.drop('y', axis = 1)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "14bb5b84",
   "metadata": {},
   "source": [
    "归一化"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "e88a94e4",
   "metadata": {},
   "outputs": [],
   "source": [
    "df = pd.DataFrame(X_1)\n",
    "X = df.apply(lambda x: (x - np.min(x)) / (np.max(x) - np.min(x)))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b8301586",
   "metadata": {},
   "source": [
    "切割数据"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "3cbad73e",
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.model_selection import train_test_split\n",
    "\n",
    "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size = 0.2,\n",
    "                                                   random_state = 1)    # 把样品分为测试和验证集，用split进行分割，测试集的比例为0.2，random_state意思是随机切割、抽取"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f8b95152",
   "metadata": {},
   "source": [
    "1. 进行线性回归算法【不交叉验证】"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "3163e939",
   "metadata": {},
   "outputs": [],
   "source": [
    "#进行线性回归\n",
    "\n",
    "from sklearn.linear_model import LinearRegression  #线性回归拟合\n",
    "\n",
    "lr = LinearRegression()\n",
    "lr.fit(X_train, y_train)\n",
    "\n",
    "y_pred_train_lr = lr.predict(X_train)\n",
    "y_pred_test_lr = lr.predict(X_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "d4c826a6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "training R2 = 0.974\n",
      "test R2 = -5.719\n",
      "RMSE(training)2.050\n",
      "RMSE(test)21.640\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAFzCAYAAADL1PXCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAhPUlEQVR4nO3df5BdZZ3n8fc3IQRQJkQC8qNpE50AoraALaLi7OiAgEsR6HWcDOOKNdRGGVit3YpKhlKxXGoZzTgu4+gM1rJLjTpUdHoW/B1j6eAMMDHB0BJ+mCjYNFGIscn4i0iS7/5xbpOb5J7uS9L3nnv7vl9VXffe55x77/dE7E8/53nOcyIzkSSpkVlVFyBJ6lyGhCSplCEhSSplSEiSShkSkqRShoQkqdQhVRcwHRYsWJALFy6sugxJ6irr16//WWYeM9k+MyIkFi5cyLp166ouQ5K6SkT8eKp9Kj3dFBE3R8QTEXFfXdvzIuIbEbGp9ji/yholqZdVPSbxf4EL9mm7BvhmZi4Gvll7LUmqQKUhkZl3AD/fp3kJcEvt+S3AJe2sSZK0RyeOSTw/M38CkJk/iYhjD+RDnn76acbGxnjqqaemt7oZ5rDDDqOvr485c+ZUXYqkDtSJIdGUiFgGLAPo7+/fb/vY2BhHHnkkCxcuJCLaXV5XyEy2bdvG2NgYixYtqrocSR2o6jGJRh6PiOMBao9PNNopM2/KzMHMHDzmmP1ncD311FMcffTRBsQkIoKjjz7a3pakUp0YErcDl9eeXw7cdqAfZEBMzX8jSZOpegrsPwB3AadExFhEXAHcAJwXEZuA82qvu9KTTz7JJz/5yWf9vje96U08+eSTk+7zgQ98gDVr1hxgZZLUnJgJNx0aHBzMfS+me+CBB3jxi19cUUWFRx55hIsuuoj77rtvr/Zdu3Yxe/bsiqraXyf8W0lqv4hYn5mDk+3TtQPX021kBIaHYXQU+vthaAgGBg7uM6+55hp++MMfcvrppzNnzhye+9zncvzxx7Nhwwbuv/9+LrnkEh599FGeeuop3v3ud7Ns2TJgzxXkv/zlL7nwwgs555xzuPPOOznxxBO57bbbOPzww3n729/ORRddxJvf/GYWLlzI5Zdfzhe/+EWefvppPv/5z3PqqaeydetWLrvsMrZt28YrX/lKvva1r7F+/XoWLFgwDf9ikqrUit9ZjXTimETbjYzAypUwPg59fcXjypVF+8G44YYbeNGLXsSGDRv46Ec/ytq1a7n++uu5//77Abj55ptZv34969at48Ybb2Tbtm37fcamTZu46qqr2LhxI0cddRT/+I//2PC7FixYwD333MOVV17JypUrAfjQhz7EG97wBu655x4uvfRSRkdHD+6AJHWEVv3OasSQoEjj+fOLn1mz9jwfHp7e7znrrLP2mmp644038vKXv5yzzz6bRx99lE2bNu33nkWLFnH66acD8IpXvIJHHnmk4WcPDQ3tt8+//Mu/sHTpUgAuuOAC5s93hRNpJmjX7ywwJICiuzZv3t5t8+YV7dPpOc95zjPPv/3tb7NmzRruuusu7r33Xs4444yGU1Hnzp37zPPZs2ezc+fOhp89sV/9PjNhvEnS/tr1OwsMCaA4n7d9+95t27cX7QfjyCOP5Be/+EXDbdu3b2f+/PkcccQRPPjgg9x9990H92UNnHPOOaxatQqA1atXMz4+Pu3fIan9WvU7qxFDgmLAZ3y8+Nm9e8/z2hmcA3b00Ufz2te+lpe+9KW85z3v2WvbBRdcwM6dOxkYGOD9738/Z5999sF9WQMf/OAHWb16NWeeeSZf/epXOf744znyyCOn/XsktVerfmc14hTYmnbNFGinHTt2MHv2bA455BDuuusurrzySjZs2LDffk6BlbrPdPzOcgrsszAw0P2hsK/R0VHe8pa3sHv3bg499FA+/elPV12SpGnSrt9ZhsQMtnjxYr73ve9VXYakLuaYhCSplCEhSSplSEiSShkSkqRShkQLHehS4QAf//jH+fWvfz3NFUnSs2NItJAhIanbOQV2QguupqtfKvy8887j2GOPZdWqVezYsYNLL72UD33oQ/zqV7/iLW95C2NjY+zatYv3v//9PP7442zZsoXXv/71LFiwgG9961vTdJCS9OwYErBn3d358/ded3f58oMKihtuuIH77ruPDRs2sHr1ar7whS+wdu1aMpOLL76YO+64g61bt3LCCSfw5S9/GSjWdJo3bx4f+9jH+Na3vuW9HyRVytNN0JZ1d1evXs3q1as544wzOPPMM3nwwQfZtGkTL3vZy1izZg3ve9/7+M53vsO8fZd2lKQK2ZOA4hRTX9/ebdO87m5msmLFCt7xjnfst239+vV85StfYcWKFbzxjW/kAx/4wLR9ryQdDHsS0LJ1d+uXCj///PO5+eab+eUvfwnAY489xhNPPMGWLVs44ogjeOtb38ry5cu555579nuvJFXFngQUg9S1W34yb14REOPjcMUVB/Wx9UuFX3jhhVx22WW8+tWvBuC5z30un/nMZ9i8eTPvec97mDVrFnPmzOFTn/oUAMuWLePCCy/k+OOPd+BaUmVcKnzCTFwrvEkuFS71JpcKfzZm4lrhknSQHJOQJJUyJCRJpWZ0SMyE8ZZW899I0mRmbEgcdthhbNu2zV+Ck8hMtm3bxmGHHVZ1KZI61IwduO7r62NsbIytW7dWXUpHO+yww+jb90JCSaqZsSExZ84cFi1aVHUZktTVZuzpJknSwTMkJEmlDAlJUqmODYmIuCAiHoqIzRFxTdX1SFIv6siQiIjZwN8AFwKnAX8cEadVW5Uk9Z6ODAngLGBzZv4oM38L3AosqbgmSeo5nRoSJwKP1r0eq7VJktqoU0MiGrTtdel0RCyLiHURsc4L5iSpNTo1JMaAk+pe9wFb6nfIzJsyczAzB4855pi2FidJvaJTQ+K7wOKIWBQRhwJLgdsrrkmSek5HLsuRmTsj4mrg68Bs4ObM3FhxWZLUczoyJAAy8yvAV6quQ5J6WaeebpIkdQBDQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTKkJAklTIkJEmlDAlJUilDQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTKkJAklTIkJEmlDAlJUilDQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUqCYmI+MOI2BgRuyNicJ9tKyJic0Q8FBHnV1GfJKlwSEXfex8wBPxdfWNEnAYsBV4CnACsiYiTM3NX+0uUJFUSEpn5AEBE7LtpCXBrZu4AHo6IzcBZwF3trVCVGhmB4WEYHYX+fhgagoGBqquSelKnjUmcCDxa93qs1qZeMTICK1fC+Dj09RWPK1cW7ZLarmU9iYhYAxzXYNO1mXlb2dsatGXJ5y8DlgH09/cfUI3qQMPDMH9+8QN7HoeH7U1IFWhZSGTmuQfwtjHgpLrXfcCWks+/CbgJYHBwsGGQqAuNjhY9iHrz5hXtktqu00433Q4sjYi5EbEIWAysrbgmtVN/P2zfvnfb9u1Fu6S2q2oK7KURMQa8GvhyRHwdIDM3AquA+4GvAVc5s6nHDA0V4xDj47B7957nQ0NVVyb1pMjs/jM1g4ODuW7duqrL0HRxdpPUFhGxPjMHJ9unquskpHIDA4aC1CE6bUxCktRB7ElIagnPGs4M9iQkTTuviZw5DAlJ067+mshZs/Y8Hx6uujI9W4aEpGk3OlpcA1nPayK7kyEhadp5TeTMYUhImnZeEzlzGBKSpt3AACxfXoxDjI0Vj8uXO7upGzkFVlJLeE3kzGBPQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTKkJAklXJZDqlLeec3tYM9CakLeec3tYshIXUh7/ymdjEkpC7knd/ULoaE1IW885vaxZCQupB3flO7GBJSF/LOb2oXp8BKXco7v6kd7ElIkkoZEpKkUoaEJKmUISFJKlVJSETERyPiwYgYiYh/ioij6ratiIjNEfFQRJxfRX2SpEJVPYlvAC/NzAHgB8AKgIg4DVgKvAS4APhkRMyuqEZJ6nmVhERmrs7MnbWXdwN9tedLgFszc0dmPgxsBs6qokZJUmeMSfwp8NXa8xOBR+u2jdXaJEkVaNnFdBGxBjiuwaZrM/O22j7XAjuBz068rcH+WfL5y4BlAP0uWCNJLdGykMjMcyfbHhGXAxcBf5CZE0EwBpxUt1sfsKXk828CbgIYHBxsGCSSpINT1eymC4D3ARdn5q/rNt0OLI2IuRGxCFgMrK2iRklSdWs3fQKYC3wjIgDuzsx3ZubGiFgF3E9xGuqqzNxVUY2S1PMqCYnM/N1Jtl0PXN/GciRJJTphdpMkqUMZEpKkUoaEJKmUISFJKmVISJJKGRKSpFJNT4GNiBOBF9S/JzPvaEVRkqTO0FRIRMRfAH9EcZHbxMVtCRgSmvFGRmB4GEZHob8fhoZgYKDqqqT2aLYncQlwSmbuaGEtUscZGYGVK2H+fOjrg/Hx4vXy5QaFekOzIfEjYA5gSGh6dfif6cPDRUDMn1+8nngcHu6oMqWWaXbg+tfAhoj4u4i4ceKnlYWpB0z8mT4+vvef6SMjVVf2jNFRmDdv77Z584p2qRc025O4vfYjTZ8u+DO9v7/IronSALZvL9qlXtBUSGTmLRFxKHByremhzHy6dWWpJ4yOFj2Ieh32Z/rQUNG5gaK07duL0LjiimrrktqlqdNNEfH7wCbgb4BPAj+IiN9rXVnqCf39xW/deh32Z/rAQDFIPX8+jI0Vjw5aq5c0e7rpL4E3ZuZDABFxMvAPwCtaVZh6QJf8mT4wYCiodzU7cD1nIiAAMvMHFLOdpAPnn+lSx2u2J7EuIv438Pe1138CrG9NSeop/pkudbRmQ+JK4CrgXUBQXGn9yVYVJUnqDM3ObtoBfKz2I0nqEZOGRESsysy3RMT3KdZq2ktmep5AkmawqXoS7649XtTqQiRJnWfSkMjMn9Se/llmvq9+W21l2Pft/y6pcxzo0lAdvqSU1DbNToE9r0HbhdNZiDTdDnRpqC5YUkpqm0lDIiKurI1HnBoRI3U/DwPfb0+J0oGpXxpq1qw9z4eHW/M+aSaaakzic8BXgf8JXFPX/ovM/HnLqpKmwYEuDdUFS0pJbTPVmMR2YHtE/C/g55n5C4CIODIiXpWZ/9aOIqVGpho3ONAVXF35Vdqj2TGJTwG/rHv9q1qbVIlmxg2Ghor28XHYvXvP86GhyT/7QN8nzUTNhkRk5jPXSWTmbpq/Wluads2MGxzo0lAuKSXt0fTtSyPiXezpPfwZxS1NpUo0O25woEtDTfo+58eqhzTbk3gn8BrgMWAMeBWwrFVFSVOp7FYUzo9Vj2kqJDLzicxcmpnHZubzM/OyzHyi1cVJZSobN3B+rHrMVGs3vTczPxIRf03jtZve1bLKpElMjBvUn/W54oo2nPVxfqx6zFRjEg/UHtdN55dGxIeBJcBu4Ang7Zm5pbZtBXAFsAt4V2Z+fTq/Wx1gms7pV3IrCufHqsdE3aSl9n1pxO9k5r/Xnr8LOC0z3xkRp1HcFvUs4ARgDXByZu6a7PMGBwdz3bppzTG1ysQ5/fnz975labdMH+r2+qU6EbE+Mwcn22eq001fpMFppgmZefGBFDYREDXPqfuOJcCttftXPBwRmykC464D+R51oPpz+vDM408/NczfPn+g8ycMVXaeS6rGVKebanepZwg4DvhM7fUfA48czBdHxPXA24DtwOtrzScCd9ftNlZr00zR4Jz+40/N4wffHGX8P+49Yahj/zj3lqvqIZPObsrMf87MfwbOyMw/yswv1n4uA86Z7L0RsSYi7mvws6T22ddm5knAZ4GrJ97WqIySz18WEesiYt3WrVunOk51igZzVx/esJ1fHd3vhCGpAzV7ncQxEfHCiRcRsQg4ZrI3ZOa5mfnSBj+37bPr54D/VHs+BpxUt60P2FLy+Tdl5mBmDh5zzKSlqJM0mLu6a9s4Pzx977mrThiSOkOzV1z/N+DbETFxlfVC4B0H+qURsTgzN9VeXgw8WHt+O/C5iPgYxcD1YmDtgX6POlCDc/r3nnsFm+YOUDdf6JkJQ17cLFWrqZDIzK9FxGLg1FrTg7XB5QN1Q0ScQjEF9scUV3STmRsjYhVwP7ATuGqqmU3qQvuc0z9nBO6ujX7VTxh63ev2TCTqirEKaQZqagpsRBwB/HfgBZn5X2qBcUpmfqnVBTbDKbDdr1GPYXh4/0sSJl5fd11lpUozxkFPga3zf4D1wKtrr8eAzwMdERLqfo0mDH38417cLFWt2YHrF2XmR4CnATLzNzSeiSRNm8oW8ZP0jGZD4rcRcTi16agR8SLgYMYkpCl58x+pes2GxAeBrwEnRcRngW8C721ZVRLe/EfqBFOOSUTELGA+xVXXZ1OcZnp3Zv6sxbVpBimbyjrVFFcvbpaq1ezspjsy8/faUM8BcXZTZytbE+/ii+H2210rT6pKM7Obmj3d9I2IWB4RJ0XE8yZ+pqFG9YCy+/R84hPev0fqdM1Ogf1TikHrP9un/YUN9pX2Ur+m3+OPwwMPwJNPFu2nnLL3vk5xlTpLsz2J04C/Ae4FNgB/DbykRTVphpmYyvr443DnnfCb38Chh8Lhh8Mdd8BPf7pnX6e4Sp2l2ZC4BXgxcCNFQLy41iZNaWIq6z33wNy5RduOHfCa1xTPv/c9p7hKnarZ002nZObL615/KyLubUVBPacHVrCbmMr6trcVr486Cs44A447Dp73PFi7tpji6v17pM7TbEh8LyLOzsy7ASLiVcC/tq6sHlE/7WeGr2A3MACXXLL/WkyHHQZLlrgWk9Spmj3d9Crgzoh4JCIeobid6H+IiO9HxEjLqpvpyqb9zNDpPV5BLXWfZnsSF7S0il7V4FaeM3l6j7eHlrpPs/eT+HGrC+lJ/f37n3/pwOk90zls4hXUUndp9nSTWqELzr9MDJuMj+89bDLiSUapJxgSVergFexGRorB5Le9DR56CH77254YNpG0j2bHJNQqHXj+pX7SFUBmcRHca14Dz3/+jB42kbQPexLaT/2kq6OOgohiquoDDxTbO3DYRFKLGBLaz+ho0VsAOPVUeOqpojfx5JMdOWwiqYUMCe2n/rahxx0Hr3510ZuI6KhhE0lt4JiE9jM0VIxJQNGjmDu3WK3VcJB6jz0J7aeDJ11JajN7EmqoAyddSaqAPQlJUil7Eh2mB1YOl9RFDIkO0qkrhxtcUu/ydFMH6cSVw127SepthkQHqb+IbULVS2B0YnBJah9DooPUX8Q2oeolMDoxuCS1jyHRQTpx5fBODC5J7WNIdJBOvIitE4NLUvtEZlb35RHLgY8Cx2Tmz2ptK4ArgF3AuzLz61N9zuDgYK5bt66ltfYyZzdJM1NErM/Mwcn2qWwKbEScBJwHjNa1nQYsBV4CnACsiYiTM3NXNVUKvPpa6mVVnm76K+C9QH1XZglwa2buyMyHgc3AWVUUJ0mqKCQi4mLgscy8d59NJwKP1r0eq7VJkirQstNNEbEGOK7BpmuBPwfe2OhtDdoaDppExDJgGUC/U20kqSVaFhKZeW6j9oh4GbAIuDciAPqAeyLiLIqew0l1u/cBW0o+/ybgJigGrqevcknShLafbsrM72fmsZm5MDMXUgTDmZn5U+B2YGlEzI2IRcBiYG27a5QkFTpqgb/M3BgRq4D7gZ3AVc5s2sOpqJLarfKL6Wo9ip/Vvb4+M1+Umadk5lerrK2TuNCepCpUHhJqjgvtSaqCIdElXGhPUhUMiS7hQnuSqmBIdIl9F9r7wQ/g29+GDRvguuscm5DUGoZEl6hfIXZkBDZuhJe+tGh3EFtSq3TUFFhNbmKhveuugxe8oAgM2PM4POyUWEnTy55EF3IQW1K72JPoEvUX0v3oR7BjByxeDI8/Dg88AE88AcceW+xnb0LSdLEn0QX2vZDuhBPgrrvgu9+Ff/1XePJJOOSQot2xCUnTyZ5EF6i/kA7g5JOLxzvvhCOOKHoQp54Kxx1XBIljE5KmiyHRBUZHix5Evd/93aIn8eY3F1dgT3BsQtJ08nRTFyi7kO7EE73ATlJrGRJdYN8L6SaeX3114/ahoaorljRTGBJdoP5CurGx4nH58uJUU6N2xyMkTRfHJLrExIV0zbZL0nSwJyFJKmVPYobwrnWSWsGexAzgXesktYohMQN41zpJrWJIzAAu+CepVQyJGcC71klqFUNiBii72M6L6iQdLENiBii72M7ZTZIOllNgZwgvqpPUCvYkJEmlDAlJUilDQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVqiQkIuK6iHgsIjbUft5Ut21FRGyOiIci4vwq6pMkFaq84vqvMnNlfUNEnAYsBV4CnACsiYiTM3NXFQVKUq/rtNNNS4BbM3NHZj4MbAbOqrgmSepZVYbE1RExEhE3R8T8WtuJwKN1+4zV2iRJFWhZSETEmoi4r8HPEuBTwIuA04GfAH858bYGH5Uln78sItZFxLqtW7e24hAkqee1bEwiM89tZr+I+DTwpdrLMeCkus19wJaSz78JuAlgcHCwYZBIkg5OVbObjq97eSlwX+357cDSiJgbEYuAxcDadtcnSSpUNbvpIxFxOsWppEeAdwBk5saIWAXcD+wErnJmkyRVp5KQyMz/PMm264Hr21HHyAgMD8PoaHE/6KEhb9wjSfU6bQps24yMwMqVxb2g+/qKx5Uri3ZJUqFnQ2J4uLgX9Pz5MGvWnufDw1VXJkmdo2dDYnQU5s3bu23evKJdklTo2ZDo74ft2/du2769aJckFXo2JIaGinGI8XHYvXvP86GhqiuTpM7RsyExMADLlxfjEGNjxePy5c5ukqR6Va4CW7mBAUNBkibTsz0JSdLUDAlJUilDQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTKkJAklTIkJEmlDAlJUilDQpJUypCQJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTKkJAklTIkJEmlKguJiPivEfFQRGyMiI/Uta+IiM21bedXVZ8kCQ6p4ksj4vXAEmAgM3dExLG19tOApcBLgBOANRFxcmbuqqJOSep1VfUkrgRuyMwdAJn5RK19CXBrZu7IzIeBzcBZFdUoST2vqpA4GXhdRPxbRPxzRLyy1n4i8GjdfmO1NklSBVp2uiki1gDHNdh0be175wNnA68EVkXEC4FosH+WfP4yYBlAf3//dJQsSdpHy0IiM88t2xYRVwLDmZnA2ojYDSyg6DmcVLdrH7Cl5PNvAm4CGBwcbBgkkqSDU9Xppv8HvAEgIk4GDgV+BtwOLI2IuRGxCFgMrK2oRknqeZXMbgJuBm6OiPuA3wKX13oVGyNiFXA/sBO4yplNklSdSkIiM38LvLVk2/XA9e2tSJLUiFdcS5JKGRKSpFKGhCSplCEhSSplSEiSShkSkqRShoQkqZQhIUkqZUhIkkoZEpKkUlWt3SSpWSMjMDwMo6PQ3w9DQzAwUHVV6hH2JKRONjICK1fC+Dj09RWPK1cW7VIbGBJSJxsehvnzi59Zs/Y8Hx6uujL1CENC6mSjozBv3t5t8+YV7VIbGBJSJ+vvh+3b927bvr1ol9rAkJA62dBQMQ4xPg67d+95PjRUdWXqEYaE1MkGBmD58mIcYmyseFy+3NlNahunwEqdbmDAUFBl7ElIkkoZEpKkUoaEJKmUISFJKmVISJJKGRKSpFKGhCSplCEhSSplSEiSShkSkqRSkZlV13DQImIr8OMWffwC4Gct+uwqzKTj8Vg610w6npl8LC/IzGMme8OMCIlWioh1mTlYdR3TZSYdj8fSuWbS8fT6sXi6SZJUypCQJJUyJKZ2U9UFTLOZdDweS+eaScfT08fimIQkqZQ9CUlSKUOiRERcFxGPRcSG2s+b6ratiIjNEfFQRJxfZZ3PRkQsj4iMiAV1bV11LBHx4YgYqf1vsjoiTqjb1lXHAhARH42IB2vH9E8RcVTdtq46noj4w4jYGBG7I2Jwn21ddSwAEXFBrd7NEXFN1fU8WxFxc0Q8ERH31bU9LyK+ERGbao/zp/ygzPSnwQ9wHbC8QftpwL3AXGAR8ENgdtX1NnE8JwFfp7ieZEG3HgvwO3XP3wX8bbceS63uNwKH1J7/BfAX3Xo8wIuBU4BvA4N17d14LLNrdb4QOLRW/2lV1/Usj+H3gDOB++raPgJcU3t+zcR/b5P92JN49pYAt2bmjsx8GNgMnFVxTc34K+C9QP0gVNcdS2b+e93L57DneLruWAAyc3Vm7qy9vBvoqz3vuuPJzAcy86EGm7ruWCjq25yZP8rM3wK3UhxH18jMO4Cf79O8BLil9vwW4JKpPseQmNzVtdMAN9d1y04EHq3bZ6zW1rEi4mLgscy8d59NXXcsABFxfUQ8CvwJ8IFac1ceyz7+FPhq7flMOJ4J3Xgs3VhzM56fmT8BqD0eO9UbDml5SR0sItYAxzXYdC3wKeDDFH+pfhj4S4r/E0eD/SufIjbFsfw5xWmN/d7WoK2jjyUzb8vMa4FrI2IFcDXwQTr0WGDq46ntcy2wE/jsxNsa7F/58TRzLI3e1qCt8mOZQjfW3BI9HRKZeW4z+0XEp4Ev1V6OUZzfn9AHbJnm0p61smOJiJdRnAe+NyKgqPeeiDiLLjuWBj4HfJkiJDryWGDq44mIy4GLgD/I2sliOvR4nsX/NvU68lim0I01N+PxiDg+M38SEccDT0z1Bk83laj9A064FJiYIXA7sDQi5kbEImAxsLbd9TUrM7+fmcdm5sLMXEjxH/+ZmflTuuxYACJicd3Li4EHa8+77ligmEEDvA+4ODN/XbepK4+nRDcey3eBxRGxKCIOBZZSHEe3ux24vPb8cqCs9/eMnu5JTOEjEXE6RRfzEeAdAJm5MSJWAfdTnB64KjN3VVXkwejSY7khIk4BdlPM1HondO2xAHyCYtbPN2o9vbsz853deDwRcSnw18AxwJcjYkNmnt+Nx5KZOyPiaooZgbOBmzNzY8VlPSsR8Q/A7wMLImKMosd9A7AqIq4ARoE/nPJz9vRuJUnam6ebJEmlDAlJUilDQpJUypCQJJUyJCRJpZwCKx2EiNgFfB+YQzG98xbg45m5OyJ+n2Ie+sO13X8GPAS8lmLRuEW11wD/IzO/0L7KpeYYEtLB+U1mng4QEcdSXAU+j2JOOsB3MvOifd8UEQuBL028V+pUnm6SpklmPgEso1gYstHaP1LXMSSkaZSZP6L4/9XE6pqvq7tx1bUVliYdEE83SdOvvhfR8HST1C3sSUjTKCJeCOyiidU1pW5gSEjTJCKOAf4W+ES6KJpmCE83SQfn8IjYwJ4psH8PfKzSiqRp5CqwkqRSnm6SJJUyJCRJpQwJSVIpQ0KSVMqQkCSVMiQkSaUMCUlSKUNCklTq/wMDxRyXmKSqLQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#作图\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "plt.figure(figsize = (6, 6))\n",
    "plt.scatter(y_train, y_pred_train_lr, alpha = 0.5, color = 'b', label = 'training')  #alpha设置透明度\n",
    "plt.scatter(y_test, y_pred_test_lr, alpha = 0.5, color = 'r', label = 'test')\n",
    "plt.legend()\n",
    "plt.xlabel('DFT')\n",
    "plt.ylabel('prediction')\n",
    "\n",
    "plt.savefig(r'C:/Users/Administrator.DESKTOP-CDML10A/Desktop/线性回归不交叉验证.tiff', dpi = 360)\n",
    "\n",
    "#算均方误差\n",
    "from sklearn.metrics import mean_squared_error #引用均方差\n",
    "\n",
    "print('training R2 = %.3f' %lr.score(X_train, y_train))\n",
    "print('test R2 = %.3f' %lr.score(X_test, y_test))\n",
    "rmse_tr_lr = mean_squared_error(y_train, y_pred_train_lr, squared = False) #squared = False 是表示对均方差求根，及得到均方根误差\n",
    "rmse_te_lr = mean_squared_error(y_test, y_pred_test_lr, squared = False)\n",
    "print('RMSE(training)%.3f'%rmse_tr_lr)\n",
    "print('RMSE(test)%.3f'%rmse_te_lr)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "31c4b711",
   "metadata": {},
   "source": [
    "2. 进行线性回归算法【7折交叉验证】"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "c81e3a92",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cross-validation results:\n",
      "Fold: 7, mean R2: -9.910\n",
      "Fold: 7, mean RMSE: 22.771\n"
     ]
    }
   ],
   "source": [
    "#进行线性回归\n",
    "\n",
    "from sklearn.linear_model import LinearRegression  #线性回归拟合\n",
    "\n",
    "lr = LinearRegression()\n",
    "lr.fit(X, y)\n",
    "\n",
    "from sklearn.model_selection import KFold, cross_val_score\n",
    "\n",
    "crossvalidation = KFold(n_splits = 7, shuffle = True)  #shuffle是进行打乱的意思，在交叉验证之前进行打乱\n",
    "r2_scores = cross_val_score(lr, X, y, scoring = 'r2', cv = crossvalidation)  #cv是分类方式，这里是交叉验证\n",
    "rmse_scores = cross_val_score(lr, X, y, scoring = 'neg_root_mean_squared_error',\n",
    "                             cv = crossvalidation)\n",
    "\n",
    "print('Cross-validation results:')\n",
    "print('Fold: %i, mean R2: %.3f' %(len(r2_scores), np.mean(r2_scores)))\n",
    "print('Fold: %i, mean RMSE: %.3f' %(len(rmse_scores), -np.mean(rmse_scores)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "3092789c",
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.model_selection import cross_val_predict\n",
    "y_cv = cross_val_predict(lr, X, y, cv = crossvalidation)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "36611bbc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUcAAAEvCAYAAADW7gNcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAoBUlEQVR4nO3deXhV5bn38e+diYQ5YZAhAqkgCOLrENGgVFFAxAkEW2q99NWecnDs0dap+Gqt5RynYhWpii09oFKKgK2tA4JDrRpEUFGgRUFQIlCBRMKwAxme949nBzZhJZBhZ+9k/z7XlYu919rDDeKPZ61nrfsx5xwiInKwpFgXICISjxSOIiIBFI4iIgEUjiIiARSOIiIBFI4iIgFSYl3AkerYsaPr1atXrMsQkWZm+fLl25xznapubzLh2KtXL5YtWxbrMkSkmTGzL4O267BaRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRJqHtWsb9OMUjiLS9D36KPTvDwUFDfaRTebeahGR/fbtg9mzoU8fOOMMGDsWkpIgK6vBvkIjRxFpOnbuhClT4Jhj4Oqr4Zln/PbsbLjxRmjZssG+SuEoIk3DI49Ajx7w059C797wyivwxBNR+zqFo4jEry++gNJS/7i8HM45B95/H958E0aOBLOofbXCUUTiz0cfwfjx/pzinDl+209/CvPnw6BBjVKCwlFE4oNz8PrrMGIEnHwyvPyyD8Rzz/X7ozhKDKLZahGJHz/5CWzbBv/zPzBxIrRvH7NSFI4iEhslJX62+emnYfFiaNsWFizwky7p6bGuTofVItLIduyABx6AnByYMAEqKmDTJr/v2GPjIhhBI0cRaUxbtkDfvlBcDMOHw7PP+hnoRj6feCQUjiISXZ99Bvn5cNVV0KUL3HorjBrlJ13imA6rRSQ6li71t/X16wc33QS7dvntd90V98EIGjlKPYVCIfLz8ykqKiIzM5O8vDwyMjJiXZbE0sqV/la+t97ys80//7kPx9atY11ZrSgcpU6cc8ybv4BpM57Ftc8mqVUmFbuLsIencv01VzBu7KVYHJ5HkigpK4Pt2+Goo6BVK1i/Hn79a/jxj6FNm1hXVycxC0czGwk8CiQDv3PO3R+rWqT25s1fwJSZC8gZOZGMtgc6oYSKC5kycxYAl40bG6vypLHs2QMzZvggPO44f+F2To6/7S+paZ+1i0n1ZpYMTAPOB/oDPzCz/rGoRWovFAoxbcaz5Ay/8qBgBMhom0XO8Cv57YznKCkpiVGFEnXbt8Mvfwk9e/pD6G7d/EXblZp4MELsJmQGAWudc1845/YBc4BLYlSL1FJ+fj6uffYhwVgpo20WLrM7+fn5jVyZNJqnnoJ77oHTToN//APefRcuvjjWVTWoWIVjd2BjxPOC8DZpAoqKikhqlVnja6xVJoWFhY1UkUTdqlX+Upz58/3za6+FTz+Fv/0NzjwztrVFSazCMehMvTvkRWYTzGyZmS3bunVrI5QlRyIz00++1MTtLiKrAbsyS4y88w5cdBEcfzzMm3dgGYLMTL+tGYtVOBYAR0c8zwY2VX2Rc266cy7XOZfbqVOnRitOapaXl4d9W0CoOHhkGCouxIq+Ji8vr5ErkwZ11VUwZIi/gPvee+Grr3xjiAQRq3D8AOhjZjlmlgaMB16MUS1SSxkZGVx/zRWsXzTrkIAMFReyYfEsrrvmh6THyT2ycoT27YNZs2D3bv/8ootg6lQfinffDR06xLa+RhaTS3mcc2VmdgOwEH8pzwzn3KpY1CJ1M27spQBMm/EkZGZjrTJxu4ugqICbw9c5ShOxc6fvjPPII/6w2Tk/ahw3LtaVxZQ5d8ipvriUm5vrli1bFusypIpQKMSSJUsoLCwkKyuLvLw8jRibirIy+MUvYNo0+PZbOPtsuP12OO+8uGwEES1mttw5l1t1u+6QkXrJyMhg6NChsS5DaqOoyE+opKTA22/7rji33eYvy5H9FI4iieLjj30fxb/+Fdau9R1yFi+GtLRYVxaXmv5l7CJSPefgjTf8ofJJJ8FLL8F11/lRIygYa6CRYwJRB50E9OWXMGwYdO4cF+uyNCUKxwSgDjoJpHJdltWr/exzr17wyitw1llxs/xAU6FwbOZCoRD3P/Agc19fSs9zr6BLjxySk5L9PnXQaT527IAnn4Tf/MYvRXDqqT4o09P9IbXUms45NlPOOZ6fN58Ro7/P47PmUzZgFF8WlfD+so/YtHkzDnXQaTYWLfIr9t1xBwwc6CdZ3n9fI8V6Ujg2U5X9FlP7nEGr3rm0634srTpmk96pJ18UbGHz5s2AOug0WZ99Bh9+6B+feKK/m2X5cnjtNTj33IS6TjFaFI7NUGS/RQOSWx9oAJGckkrLjtls2Pg15RUVgDroNCmR67L89Kd+W6dOfhW/JrAuS1OicGyGIvstprVsQ/mug4MvOSUVUtMpKvLb1UGnCXj7bRg61F+o/cYbfl2WP/0p1lU1awrHZiiy32Ln3gMp3/4VpbuqtBhLTqV0X6k66MSzsjIoLfWPP/0UPv/cL0fw1Vfwq1/5y3MkahSOzVBkv8WUtHT65o1g2z9mHxyQ5aVU7N2tDjrxaM8eePxx6NPHr88C8B//4ddlueWWJrtgVVOjS3maoby8POzhqYSKC8lom0WvU84GYM3i6SR36IFltKWk4F906tSCm//jSnXQiReFhb4JxGOPwbZtkJcHvXv7fS1axLa2BKSRYzNUtd+imZGTO5Rh/3kPffr2I/3r5Uy8KI/X/vwnLhs3VheAx4uxY33fxNNO8+cY333XzzxLTGjk2EzV1G/xV3f8l+6KiQerV8OUKXD//dCxo/+1ZUt/raLEnPo5NnPqtxiH3nkHHnzQd8dp2dIvWjVyZKyrSljq55ig1G8xjuzd65tAvPOOX3Lg3nvh+usTbvmBpkLhGKaONRIV+/bBe+/5LtstWsAJJ8D3vw/XXONHjRK3Ej4c1bFGoiJyXZZNm2DdOt8hZ9q0WFcmRyjhw7HyHuSckRPJaHvgLhF1rJE6KSrykyzTpvnHZ5/tQ7Jnz1hXJrWU0JfyRN6DHBmMoI41UktlZf7XkhJ4+GEfikuWwJtv+skWHX00OQkdjpH3IAdRxxo5rI8/hh/84EDPxK5d/e19CxZowaomLqHDMfIe5OqoY40cwjk/Ioxcl+Xkkw+MHjt1im190iAS+pxj5D3I1VHHGjnEs8/ClVfCUUfBf/83XHut1mVphhJ65JiXl4d9W0CoOHhkqI41AvjrE59+Gv78Z/98zBh46inYsAHuvFPB2EwldDhWvQc5Uqi4UB1rEt2OHX6d5169YMIEmDvXb2/d2j/X34tmLaEPq6Hme5BvDl/nKAnot7/1o8LiYhg+3B9Kn3NOrKuSRhS1cDSzh4CLgH3AOuBq59y34X13Aj8CyoGbnHMLo1XHEdTJZePGcuEFo3QPcqL77DN/HrFdO39L3/nnw223afmBBBW1xhNmNgJ4wzlXZmYPADjnbjez/sAfgUFAN2AxcKxzrrymz1PjCYmapUv94fMLL8DkyX7EKAmjusYTUTvn6Jx7zTkXvraBJUB2+PElwBzn3F7n3HpgLT4oRRrXwoWHrsvyox/FuiqJE411zvEaoHI1oO74sKxUEN4mTUiTbdTh3IG7VaZM8euyPPywn2DR8gMSoV7haGaLgS4BuyY55/4Sfs0koAx4rvJtAa8PPLY3swnABIAePXrUp1Sph8ggbN++PVu2/Junn3u+aTXq2LPHr8fy2GOwaJG/1/kPf/BNZtPSYl2dxKF6haNzblhN+83sKuBC4Fx34ORmAXB0xMuygU3VfP50YDr4c471qVVqL6hj0baCdRSs/pABw8Zx7JkX7Q/CuG3UUXVdltNPh2+/9eHYrVusq5M4Fs3Z6pHA7cBZzrk9EbteBGab2RT8hEwfYGm06pC6q9qxqLyinC3LPqL7gAtY/95c0jLakJPrG+keaNTxFBddeEF8zPTv2uUXqCoqglGj4I474Mwz1QRCjkg0LwJ/HGgDLDKzj83sSQDn3CpgLrAaeBW4/nAz1dL4qnYsKttXwmdL36T4y5Xs215A1unjWLNkEWWle/e/Jy4adaxeDQ895B+3bu1noT/5xN//PGSIglGOWNRGjs653jXsmwxMjtZ3S/1VdixKb5PJ+mVvsib/NSradqEsvT371q+k9IOXSUlJ5pvPP6Fb/1P3vy9mjTrefdcHYeW6LFdc4Tvk/PjHjV+LNAsJf4eMBKvsWLRh+Vus+fA9Og6bgEtrxa6SvSSnpVO2+1u2L3qKgpVLDgrHRm/U8fnncPXVPhw7dIBf/AJuuEHrski9KRwlUGZmJqU7trJ+nQ/G1NaZOBzsCYFzpLRqT9tTL2HL0rmUle4lJbVFYKOOqFzys28ffP015ORA585+SYLHHvPrsrRqVc/fuYincJRAeXl57Lzzbqz9d0ht7XteGkbLjHT2lISwpBSSklPJ6N6PrWs/pX3377Bh8SxuDjfqiMraPLt2+e44U6ZAVpZvNNuunf9V5xKlgSkcJVBGRgbDhuQxe/kWystKSU5JBaBFixZUlJVSvGU9LdLSCFkKX773InvapB/UqKNB1+b55huYOvXAuixnnQW3335gv4JRokDhmEBqe4g76vyRvLZqBiVbv4TUdEhOhfJSKC1hQJ8c0lu0YP3G9/i/Y4Zz3XXX7b98Z/9Md5VghDpe8vPyy/6e59GjfShq+QFpBArHBFDXQ9zBgwfTwR7nqGN7UlJulO4rJTUtlczMLJKTkggVF9LBQgcFI9RubZ6hQ4ce+oIVK/zM8+mnw003weWX+8f9+jXYn4nI4SgcE0BdD3ErmwFPmfmcv96xQ9eD3ht5jjFSndbmcQ7eesuH4sKF/hrFE0/0+9LSoF+/pns/tzRJCsdmrj6HuKFQiKzM9gzu153X59xP2579SW7TIbAZcGRwbdiwgdIdW2us65BLfm680Z9TDFiXJWjkW7ZzOzvvuo9zzziVUSPPY/Dgwc0qKPUPQewpHJu5uhziBoVR6+x+7PhiBcOG5HHB+DEMHjx4/6z08/PmHxJc/3z/DWjT8aD7ryuFigtJ2/YVZ/7zn9C3r7/H+Xvfg4ED4aqrDll+4Lk/zuH+6X/kqDPG0TqzEyV797Ll682UtTmWma/9mVeWrKJzm8fju/HFEYrKLL/UicKxmavLIW5Nh+HvvDqD1q3e3j+i2bTl30z7498OeW3qMafx8QtPAtB3yMUHvuzfGzlz5n1cvXkDqa+9CCUlcMst8N3v+p8Izjlmz/4jP/n5vWSNuI5NoST2FW6kZNcO2nTsStuuR9EyqyvbXn+agaMnMmXmHCDOGl/UUoPO8ku9KBybudouP1vdYbhzji2freCrTf/myZd20K+gjPLwCPH/XDGJ9Coj0159joMxE1nx7GSSdm8juU1Hrsl/hTFrVtC6rAw3fLifea5hXZZ58xcw+fHfk9H7NNr16I/DsdeKSW/TmT1Fm0hKTqFluw4kZx3Nzm8K4q/xRS01+Cy/1EtCrz6YCGq7/Gx1h+GVtxF2Pu86ss74Pu37n0GbXgPJ6H0aX39bwubNmw96vQE5fY7juyflMXpQH24YeSLDu2aSMXo0LF+OvfYanHtutdcoVgZFh+MGk9quMwCl+0ohKYmk1DTSsrqzq2grFa6C5NZZ7N1dHB+NL+qhNqdAJPoUjs1cbZefDToML9tXwpr81+g45HJ/t0xyKqX7Stm3Zyep7TrTsmM2GzZ+TXlFxf735Hyxiusev515L85ioBljx46ly8KFJD///BEtWFUZFK07dqV8l6+7IqKLt6WkYmkZ7Nu9k/JdhbRo1dZvj1XjiwZQp1l+iRodVseRaM1Q1mb52aDD8G/Wfkpyhx77byOk3F/v6Fq2oXzXGn/3TGo6RYXbOXvzWs5/eRb9/rWc3S3bMLv/yeT06uXfl3Tk/xZXBkXn3gP55PUFlO4qIimtlb/kJ8ySU9m3cxvlhRvp1PsaIAaNLxpQbU+BSHQpHONAtGcoa7P8bF5eHvbwVELFhfsP7/bt2Ulya/+4vMzfIZOZmYVr03p/cJGcSovib7nu8dvZ3aotf/r+T1h48lmsf+sZXjn//FrXXBkUKWnp9M0bwZp/zKbDkMuhImn/OjBlu7az55NXOX7w8GobXzQlQX/2kZr676+pUTjGgcaaoczIyAi+I6XKa/yF37P2N7pNC48Qy8tK2bOtgO8c3Z3kpCRSgVtbt+GY39/MLaNvZW/3vjx022/5qmdfdu3ZWe1F4kciMih6nXI2AGsWT/c9JdPaUrHnW0Jrl3Lq6KvpdcrZNV6U3lQE/dlXag6/v6YmautWN7Tmum51KBTi/LGX0y1ghhL8/xSbX32KVxbMbrT/KSJHsmRmU9GiNav+/hJZQ66gd99+9G7TknPemMe5i+fSdmcRKzt155K0dLoPOe+gi8TrO+p9ft58/49GRDfyb9atZOPa1Wz8YDHdevSkS/9BDfZ98aDqn33kKZDm8PuLR9WtW62RY4zV+z7kKAg6DP+kZwYLl37EqVlp/Pzen5O+N8QnJ5zBn88azavrP+SX4y+kW7euNR6y11Z150o7hgq44/676dblKIqKihrs++JBbU6BSHQpHGMsnmco9x+Gr17NpW3bcvyAATz99Exe7XEML54wmC9apMHad/dP6jT0iCaRg+JIToFIdCkcYyyuZyjfeQcefBD++ldswAAu+/TT/UF1QT2DqjYz8woKiQWFY4zF5QzlO+/4ZUyrrstiVu+g0r3D0lQoHGMsbmYo9+3zP61bw/btUFAQlXVZdO+wNBWarY4DMZ2hjFyX5aqr4Fe/gooKKC+H1NQG/ap4nJkX0Wx1HIvJxEPQuixnn+33JSXV6m6WIxWPM/Mi1VE4xpFGnXi48UaYO/fAuiynn35Eb6vPLY7xPDMvUpXCMVGsWOFnnu++2zeY/eUv4d57j3hdloaYSInrmXmRKhSOzVnQuixjxvhw7Nu3Vh/VEBMpcTkzL1KNqLcsM7OfmZkzs44R2+40s7VmtsbMzot2DQmpogKGDvXNZD/6yK/LsnEjjBtX64/a34S1ymw6RDZhfY6SkpIaP6e27dNEYimqI0czOxoYDnwVsa0/MB4YAHQDFpvZsc658mjWEo8avEXZ3r1+jecxY/yEyvDhMH68n4Wux+c25ERKbdqnxQMtdJW4on1Y/QhwG/CXiG2XAHOcc3uB9Wa2FhgEJEx74wa/EHrHDnjySfjNb2DLFli6FE49FSZNapB6G3IipancEqiL1SVq4WhmFwNfO+dWVPlL1B1YEvG8ILwtYTTYhdA7dvjD5SefhOJiGDYMnnkGcg+5ZKteojGREu+3BOpidanXOUczW2xmKwN+LgEmAXcHvS1gW+CV6GY2wcyWmdmyrVtrXge5qWiQ83ehkP81JQX+939h5EhYtgwWLfIB2cAjmtquQ9PUNdQ5Vmna6hWOzrlhzrnjq/4AXwA5wAoz2wBkAx+aWRf8SPHoiI/JBjZV8/nTnXO5zrncTp061afUuFGvRZQ++MBPqJx0kr+DpVUrWLcO/vQnOOWUqNWcaBMpWuhKIEqH1c65T4HOlc/DAZnrnNtmZi8Cs81sCn5Cpg+wNBp1xKMjPX+3ZcsW3njjDYoKCzlm3ToGvvIKyX//O7RvD9dd5ydfWrb0l+c0gqY2kVIfulhdIAbXOTrnVpnZXGA1UAZcn0gz1Yc7f+eA7RvXcf9v3qD9MSeRV7iNsQvnsLVFOpuuvJITpk7F2rZtvILDmspESkPQxeoCjRSOzrleVZ5PBiY3xnfHm8NdCL1p9QpGLH2dvkPHsmTElWyrKGd6z/680/dEPn9zDre8tiimEwHxPpHSEHSxuoDWrW501Z2/a7VrByPmTuWZKdfycOG/OXXdpwC4pGTezxtJalYXTQQ0kkQ7xyrBdPtgDFQ9f3d+wZfc9N6rZJSV8vfMTuRPuI/P+550yPvUtabxJNI5VgmmcIwBM+Oy/sdx4aMPsqSgAHv3XXakhXhz2DAeWrub7/Q7+aDXl1eUU1RYRGlpKbvKU9iyZUuMKk8ciXSOVYIpHBvbu+/6RhB//SsZN9zA0KlT/T3Qd91F+htvULHihf0vdcDmzZvZsPFrSE2H5FQK16/j/t98QEpqmu7SaASJcI5VgumcY2N5+WU480z/8957fl2We+456CVVL7bevHkzXxRsIb1TT1p1zCYtNY2U8hD9vnc7U2YuYN78BTH4jYgkBoVjNJWVHXg8f/6BdVm+/NIHY8eOB708ciJg17db2bDxa1p2zCY5JZXSXUVs+8ds+p4+nNYdNDkjEm06rI6GynVZHnkE5s2DQYPg17/2d7QcZl2WyhP9k6c8xPbSliSlt6KseBuEvmXAOWPpdcrZgCZnRKJNI8eGtHWr77Tdowfccgvk5By4z7l9+yNasMrMGDf2Ur476ET2frMeKspJ79aXjJ4nsPaDN9mw/C0qF0XTXRoi0aORY0MpK4MTT4RNm2q9LktV8+Yv4I2P19F5xH/Srkf//dtLdxWx5h+zAcjJHaq7NESiSCPH+lixAm691XfdTkmBxx+H1avhhRfqHIyVHWEGjJ5ISmoa5WWl+/elts6k45DLWbNkEbu2b9FdGiJRpHCsLefgzTd9m7ATT/S9FNes8fvGjIHjjqvXx1d2hGndvhO9ju7Onm0FhwSktevKyvmP6i4NkSjSYXVtVK7BsnQpdO4MkyfDtddCZs0dXGojsiNM165dAdiw8cv91zlSXkpFeRkjB5+kuzREokjheDh79/qR4QknQNeu0KYNPPFEvddlqU5kRxgDunXtylFHdaaoqIjSfaWkpqVS9E0mF14wSheAi0SRDqurs2OHv5OlVy8YMQJKSvx5xcWLYeLEqAQjBHfdTk5KpmOHjnTt2pVWqUkk79isc40iUaZwrGrLFrjjDn85zh13wPHHw7PPQosWjfL16ggjEh90WF3JOX9N4sqV8NBD/tzibbdFdfmB6qgjjEjsWeUFxfEuNzfXLVu2rOE/+IMPDhw+P/ywD8mvvoKePRv+u2opFAqpI4xIlJnZcufcIUt2JubI0Tm/Ut/99/vLctq189crgh89xkEwgjrCiMRSYp5zvOsuOO88Pwv90EN+pDhpUqyrEpE4khgjx1AI/vAHGDIEBg6EH/4QjjnG/9pIEy3SNIRCIfLz8ykqKiIzM5O8vDwyonRlgsS35h2ORUUwbZpvE7Z1qx8xDhwI/fv7H5Ew5xzz5i9g2oxnce2zSWrlrze1h6dyfXgSTNeVJpbmG4733OPbhO3eDaNG+UYQQ4bEuqoaadQSO/PmL2DKzAXkjJx40IqDoeJCpsycBRDTVR+l8TXfcCwt9fc633abHy3GMY1aYquy2UfVYATfN9M3Fn6Kiy68QFcLJJDmG46TJx/opRjnNGqJrcpmH0FrVIMaCyeq5jtb3USCcf+oZfiVNYxatBxCNEU2+6iOGgsnnuY7coxDQecUNWqJvchmH9VRY+HEo3BsBDWdUzxlQB+sZfsa369RS3Tl5eVhD08lVFwY+I9UqLhQjYUTkMKxEdR0TnHhgmlUlJdzzJkXV/t+jVqiq7LZx5SZsw45vVHZ7ONmNftIOFENRzO7EbgBKANecs7dFt5+J/AjoBy4yTm3MJp1xNLhZkIHjJ7I4kduZtf2LbTu0OXQ92vU0ijU7EOqilo4mtlQ4BLgBOfcXjPrHN7eHxgPDAC6AYvN7FjnXHm0aomlw51TbN2+E9375/LpvEc58Qe3atQSI2bGZePGcuEFo9TsQ4DojhyvBe53zu0FcM59E95+CTAnvH29ma0FBgH5UawlZo5kJrRD9jGcPKAby1/VqCXW1OxDKkUzHI8FhpjZZKAE+Jlz7gOgO7Ak4nUF4W2HMLMJwASAHj16RLHU6DmimdA9RVx4+aX86r7TNWoRiRP1CkczWwwceqIMJoU/OxM4HTgVmGtm38EvjVJVYFNJ59x0YDr4fo71qTVWajMTmp6erlGLSJyo10XgzrlhzrnjA37+gh8RLnDeUqAC6BjefnTEx2QDm+pTRzzTsgciTVM0D6v/DJwDvGVmxwJpwDbgRWC2mU3BT8j0AZZGsY6Y00yoSNMTzXCcAcwws5XAPuAq59dkWGVmc4HV+Et8rm+uM9WVNBMq0vRoDRkRSWjVrSHTfBtPiIjUg8JRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJIDCUUQkgMJRRCSAwlFEJEDUwtHMTjSzJWb2sZktM7NBEfvuNLO1ZrbGzM6LVg0iInWVEsXPfhC41zn3ipmNCj8/28z6A+OBAUA3YLGZHeucK49iLSIitRLNw2oHtA0/bgdsCj++BJjjnNvrnFsPrAUGBbxfRCRmojly/C9goZk9jA/hweHt3YElEa8rCG8TEYkb9QpHM1sMdAnYNQk4F7jZOTffzL4H/B4YBljA6101nz8BmADQo0eP+pQqIlIr9QpH59yw6vaZ2SzgJ+GnzwO/Cz8uAI6OeGk2Bw65q37+dGA6QG5ubmCAiohEQzTPOW4Czgo/Pgf4PPz4RWC8mbUwsxygD7A0inWIiNRaNM85/hh41MxSgBLCh8fOuVVmNhdYDZQB12umWkTiTdTC0Tn3DnBKNfsmA5Oj9d0iIvWlO2RERAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQCKBxFRAIoHEVEAigcRUQC1CsczewyM1tlZhVmlltl351mttbM1pjZeRHbTzGzT8P7HjMzq08NIiLRUN+R40rgUuDtyI1m1h8YDwwARgK/NbPk8O4ngAlAn/DPyHrWICLS4OoVjs65fzrn1gTsugSY45zb65xbD6wFBplZV6Ctcy7fOeeAWcDo+tQgIhIN0Trn2B3YGPG8ILyte/hx1e0iInEl5XAvMLPFQJeAXZOcc3+p7m0B21wN26v77gn4Q3B69OhxmEpFRBrOYcPROTesDp9bABwd8Twb2BTenh2wvbrvng5MB8jNza02REVEGlq0DqtfBMabWQszy8FPvCx1zm0GdprZ6eFZ6iuB6kafIiIxU99LecaYWQGQB7xkZgsBnHOrgLnAauBV4HrnXHn4bdcCv8NP0qwDXqlPDSIi0WB+0jj+5ebmumXLlsW6DBFpZsxsuXMut+p23SEjIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhJA4SgiEkDhKCISQOEoIhKgXuFoZpeZ2SozqzCz3Ijtw81suZl9Gv71nIh9p4S3rzWzx8zM6lODiEg01HfkuBK4FHi7yvZtwEXOuYHAVcAzEfueACYAfcI/I+tZg4hIg0upz5udc/8EqDr4c859FPF0FZBuZi2ALKCtcy4//L5ZwGjglfrUISLS0BrjnONY4CPn3F6gO1AQsa8gvC2QmU0ws2Vmtmzr1q1RLlNE5IDDjhzNbDHQJWDXJOfcXw7z3gHAA8CIyk0BL3PVvd85Nx2YDpCbm1vt60REGtphw9E5N6wuH2xm2cALwJXOuXXhzQVAdsTLsoFNdfl8EZFoisphtZm1B14C7nTOvVu53Tm3GdhpZqeHZ6mvBGocfYqIxEJ9L+UZY2YFQB7wkpktDO+6AegN/D8z+zj80zm871rgd8BaYB2ajBGROGTONY1Tebm5uW7ZsmWxLkNEmhkzW+6cy626vV6X8oiIRFsoFCI/P5+ioiIyMzPJy8sjIyMj6t+rcBSRuOScY978BUyb8SyufTZJrTKp2F2EPTyV66+5gnFjLz3kGuuGpHAUkbg0b/4CpsxcQM7IiWS0zdq/PVRcyJSZswC4bNzYqH2/Gk+ISNwJhUJMm/EsOcOvPCgYATLaZpEz/Ep+O+M5SkpKolaDwlFE4k5+fj6uffYhwVgpo20WLrM7+fn5UatB4SgicaeoqIikVpk1vsZaZVJYWBi1GhSOIhJ3MjP95EtN3O4isrKCR5YNQeEoInEnLy8P+7aAUHHwyDBUXIgVfU1eXl7UalA4ikjcycjI4PprrmD9olmHBGSouJANi2dx3TU/JD09PWo16FIeEYlL48ZeCsC0GU9CZjbWKhO3uwiKCrg5fJ1jNOn2QRGJa6FQiCVLllBYWEhWVhZ5eXkNOmLU7YMi0iRlZGQwdOjQRv9enXMUEQmgcBQRCaBwFBEJoHAUEQmgcBQRCaBwFBEJ0GSuczSzrcCXMSyhI7Atht9fE9VWN6qtbppbbT2dc52qbmwy4RhrZrYs6ELReKDa6ka11U2i1KbDahGRAApHEZEACscjNz3WBdRAtdWNaqubhKhN5xxFRAJo5CgiEkDhWAMz+4WZfW1mH4d/RkXsu9PM1prZGjM7L4Y1/szMnJl1jJfazOw+M/sk/Gf2mpl1i5fawjU8ZGb/Ctf4gpm1j5f6zOwyM1tlZhVmlltlX6xrGxn+7rVmdkdjf39APTPM7BszWxmxLcvMFpnZ5+Ffa16IpibOOf1U8wP8AvhZwPb+wAqgBZADrAOSY1Df0cBC/PWfHeOlNqBtxOObgCfjpbZwHSOAlPDjB4AH4qU+4DigL/AWkBsvf+eA5PB3fgdIC9fSv7H/21Wp6bvAycDKiG0PAneEH99R+d+2Lj8aOdbNJcAc59xe59x6YC0wKAZ1PALcBkSeOI55bc654oinrSLqi3lt4fpec86VhZ8uAbLjpT7n3D+dc2sCdsW6tkHAWufcF865fcCccE0x45x7G6i6yMwlwMzw45nA6Lp+vsLx8G4IH37NiBiidwc2RrymILyt0ZjZxcDXzrkVVXbFvDYAM5tsZhuBHwJ3x1NtVVwDvBJ+HI/1VYp1bbH+/iN1lHNuM0D41851/aCE7wRuZouBLgG7JgFPAPfhRz73Ab/G/89kAa9v8Gn/w9T2c/zh4SFvC9jWqLU55/7inJsETDKzO4EbgHsaq7YjqS/8mklAGfBc5dsao74jqS3obQHbGvNSk1h/f6NL+HB0zg07kteZ2dPA38JPC/Dn+yplA5sauLRqazOzgfjzTivMrPL7PzSzQbGuLcBs4CV8ODZKbXD4+szsKuBC4FwXPkHVWPXV4s8uUqP92cXp9x+pf5tZV+fcZjPrCnxT1w/SYXUNwn+4lcYAlbNiLwLjzayFmeUAfYCljVWXc+5T51xn51wv51wv/F/ck51zW2JdG4CZ9Yl4ejHwr/DjmNcWrm8kcDtwsXNuT8SuuKivGrGu7QOgj5nlmFkaMD5cU7x5Ebgq/PgqoLqR+GEl/MjxMB40sxPxhw8bgP8EcM6tMrO5wGr8Ydn1zrnyWBUZKU5qu9/M+gIV+Jn0iXFUG8Dj+FnfReGR9xLn3MR4qM/MxgBTgU7AS2b2sXPuvFjX5pwrM7Mb8FdHJAMznHOrGuv7g5jZH4GzgY5mVoA/OrkfmGtmPwK+Ai6r8+cfOKIQEZFKOqwWEQmgcBQRCaBwFBEJoHAUEQmgcBQRCaBwFBEJoHAUEQmgcBQRCfD/AZNO+YuTW/GLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize = (5, 5))\n",
    "plt.plot([-55, 10], [-55, 10], 'r--')\n",
    "plt.scatter(y, y_cv, s = 80, edgecolor = 'k', alpha = 0.7)\n",
    "\n",
    "plt.savefig(r'C:/Users/Administrator.DESKTOP-CDML10A/Desktop/线性回归7折交叉验证.tiff', dpi = 360)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "81ab5d70",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}