Skip to content

M
mooc-rr
Commit 941fd0a0 authored by cb1e5ba91280d02583d0604166e71c1c's avatar cb1e5ba91280d02583d0604166e71c1c
Browse files
Options

no commit message

parent 9a5c8704
Show whitespace changes
Inline Side-by-side
Showing with 1276 additions and 60 deletions
module3/exo3/exercice.ipynb
View file @ 941fd0a0
......@@ -9,7 +9,7 @@
},
{
"cell_type": "code",
"execution_count": 1,
"execution_count": 55,
"metadata": {},
"outputs": [],
"source": [
......@@ -23,7 +23,7 @@
"InteractiveShell.ast_node_interactivity = \"last_expr\" \n",
"# ‘all’|’last’|’last_expr’|’none’\n",
"#https://jupyter-console.readthedocs.io/en/4.0.1/config_options.html\n",
"pd.options.display.max_rows = 10"
"pd.options.display.max_rows = 15"
]
},
{
......@@ -35,7 +35,7 @@
},
{
"cell_type": "code",
"execution_count": 2,
"execution_count": 56,
"metadata": {},
"outputs": [
{
......@@ -145,6 +145,32 @@
" <td>314.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>1958</td>\n",
" <td>6</td>\n",
" <td>21351</td>\n",
" <td>1958.4548</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>1958</td>\n",
" <td>7</td>\n",
" <td>21381</td>\n",
" <td>1958.5370</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" <td>315.86</td>\n",
" <td>315.22</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
......@@ -158,6 +184,32 @@
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>749</th>\n",
" <td>2020</td>\n",
" <td>6</td>\n",
" <td>43997</td>\n",
" <td>2020.4563</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" </tr>\n",
" <tr>\n",
" <th>750</th>\n",
" <td>2020</td>\n",
" <td>7</td>\n",
" <td>44027</td>\n",
" <td>2020.5383</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" </tr>\n",
" <tr>\n",
" <th>751</th>\n",
" <td>2020</td>\n",
" <td>8</td>\n",
......@@ -234,7 +286,11 @@
"2 1958 3 21259 1958.2027 315.70 314.44 316.19 \n",
"3 1958 4 21290 1958.2877 317.45 315.16 317.30 \n",
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 \n",
"5 1958 6 21351 1958.4548 -99.99 -99.99 317.24 \n",
"6 1958 7 21381 1958.5370 315.86 315.19 315.86 \n",
".. ... ... ... ... ... ... ... \n",
"749 2020 6 43997 2020.4563 -99.99 -99.99 -99.99 \n",
"750 2020 7 44027 2020.5383 -99.99 -99.99 -99.99 \n",
"751 2020 8 44058 2020.6230 -99.99 -99.99 -99.99 \n",
"752 2020 9 44089 2020.7077 -99.99 -99.99 -99.99 \n",
"753 2020 10 44119 2020.7896 -99.99 -99.99 -99.99 \n",
......@@ -247,7 +303,11 @@
"2 314.91 315.70 314.44 \n",
"3 314.99 317.45 315.16 \n",
"4 315.06 317.51 314.71 \n",
"5 315.14 317.24 315.14 \n",
"6 315.22 315.86 315.19 \n",
".. ... ... ... \n",
"749 -99.99 -99.99 -99.99 \n",
"750 -99.99 -99.99 -99.99 \n",
"751 -99.99 -99.99 -99.99 \n",
"752 -99.99 -99.99 -99.99 \n",
"753 -99.99 -99.99 -99.99 \n",
......@@ -257,7 +317,7 @@
"[756 rows x 10 columns]"
]
},
"execution_count": 2,
"execution_count": 56,
"metadata": {},
"output_type": "execute_result"
}
......@@ -289,7 +349,7 @@
},
{
"cell_type": "code",
"execution_count": 3,
"execution_count": 57,
"metadata": {},
"outputs": [
{
......@@ -310,7 +370,7 @@
},
{
"cell_type": "code",
"execution_count": 4,
"execution_count": 58,
"metadata": {},
"outputs": [
{
......@@ -425,7 +485,7 @@
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 317.51 314.71"
]
},
"execution_count": 4,
"execution_count": 58,
"metadata": {},
"output_type": "execute_result"
}
......@@ -443,7 +503,7 @@
},
{
"cell_type": "code",
"execution_count": 5,
"execution_count": 59,
"metadata": {},
"outputs": [],
"source": [
......@@ -486,7 +546,7 @@
},
{
"cell_type": "code",
"execution_count": 6,
"execution_count": 60,
"metadata": {},
"outputs": [
{
......@@ -601,7 +661,7 @@
"75 1964 4 23482 1964.2896 NaN NaN 321.83 319.45 321.83 319.45"
]
},
"execution_count": 6,
"execution_count": 60,
"metadata": {},
"output_type": "execute_result"
}
......@@ -612,7 +672,7 @@
},
{
"cell_type": "code",
"execution_count": 7,
"execution_count": 61,
"metadata": {},
"outputs": [
{
......@@ -727,7 +787,7 @@
"6 1958 7 21381 1958.5370 315.86 315.19 315.86 315.22 315.86 315.19"
]
},
"execution_count": 7,
"execution_count": 61,
"metadata": {},
"output_type": "execute_result"
}
......@@ -746,7 +806,7 @@
},
{
"cell_type": "code",
"execution_count": 8,
"execution_count": 90,
"metadata": {
"scrolled": true
},
......@@ -870,6 +930,146 @@
" <td>314.93</td>\n",
" <td>316.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-09-01</th>\n",
" <td>8</td>\n",
" <td>1958</td>\n",
" <td>9</td>\n",
" <td>21443</td>\n",
" <td>1958.7068</td>\n",
" <td>313.21</td>\n",
" <td>316.08</td>\n",
" <td>312.46</td>\n",
" <td>315.35</td>\n",
" <td>313.21</td>\n",
" <td>316.08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-11-01</th>\n",
" <td>10</td>\n",
" <td>1958</td>\n",
" <td>11</td>\n",
" <td>21504</td>\n",
" <td>1958.8740</td>\n",
" <td>313.33</td>\n",
" <td>315.20</td>\n",
" <td>313.62</td>\n",
" <td>315.46</td>\n",
" <td>313.33</td>\n",
" <td>315.20</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-12-01</th>\n",
" <td>11</td>\n",
" <td>1958</td>\n",
" <td>12</td>\n",
" <td>21534</td>\n",
" <td>1958.9562</td>\n",
" <td>314.67</td>\n",
" <td>315.43</td>\n",
" <td>314.77</td>\n",
" <td>315.51</td>\n",
" <td>314.67</td>\n",
" <td>315.43</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-01-01</th>\n",
" <td>12</td>\n",
" <td>1959</td>\n",
" <td>1</td>\n",
" <td>21565</td>\n",
" <td>1959.0411</td>\n",
" <td>315.58</td>\n",
" <td>315.54</td>\n",
" <td>315.62</td>\n",
" <td>315.57</td>\n",
" <td>315.58</td>\n",
" <td>315.54</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-02-01</th>\n",
" <td>13</td>\n",
" <td>1959</td>\n",
" <td>2</td>\n",
" <td>21596</td>\n",
" <td>1959.1260</td>\n",
" <td>316.49</td>\n",
" <td>315.86</td>\n",
" <td>316.27</td>\n",
" <td>315.63</td>\n",
" <td>316.49</td>\n",
" <td>315.86</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-03-01</th>\n",
" <td>14</td>\n",
" <td>1959</td>\n",
" <td>3</td>\n",
" <td>21624</td>\n",
" <td>1959.2027</td>\n",
" <td>316.65</td>\n",
" <td>315.38</td>\n",
" <td>316.98</td>\n",
" <td>315.69</td>\n",
" <td>316.65</td>\n",
" <td>315.38</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-04-01</th>\n",
" <td>15</td>\n",
" <td>1959</td>\n",
" <td>4</td>\n",
" <td>21655</td>\n",
" <td>1959.2877</td>\n",
" <td>317.72</td>\n",
" <td>315.42</td>\n",
" <td>318.09</td>\n",
" <td>315.77</td>\n",
" <td>317.72</td>\n",
" <td>315.42</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-05-01</th>\n",
" <td>16</td>\n",
" <td>1959</td>\n",
" <td>5</td>\n",
" <td>21685</td>\n",
" <td>1959.3699</td>\n",
" <td>318.29</td>\n",
" <td>315.49</td>\n",
" <td>318.65</td>\n",
" <td>315.85</td>\n",
" <td>318.29</td>\n",
" <td>315.49</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-06-01</th>\n",
" <td>17</td>\n",
" <td>1959</td>\n",
" <td>6</td>\n",
" <td>21716</td>\n",
" <td>1959.4548</td>\n",
" <td>318.15</td>\n",
" <td>316.03</td>\n",
" <td>318.04</td>\n",
" <td>315.94</td>\n",
" <td>318.15</td>\n",
" <td>316.03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959-07-01</th>\n",
" <td>18</td>\n",
" <td>1959</td>\n",
" <td>7</td>\n",
" <td>21746</td>\n",
" <td>1959.5370</td>\n",
" <td>316.54</td>\n",
" <td>315.86</td>\n",
" <td>316.67</td>\n",
" <td>316.03</td>\n",
" <td>316.54</td>\n",
" <td>315.86</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
......@@ -882,6 +1082,16 @@
"1958-05-01 4 1958 5 21320 1958.3699 317.51 314.71 317.86 \n",
"1958-07-01 6 1958 7 21381 1958.5370 315.86 315.19 315.86 \n",
"1958-08-01 7 1958 8 21412 1958.6219 314.93 316.19 314.00 \n",
"1958-09-01 8 1958 9 21443 1958.7068 313.21 316.08 312.46 \n",
"1958-11-01 10 1958 11 21504 1958.8740 313.33 315.20 313.62 \n",
"1958-12-01 11 1958 12 21534 1958.9562 314.67 315.43 314.77 \n",
"1959-01-01 12 1959 1 21565 1959.0411 315.58 315.54 315.62 \n",
"1959-02-01 13 1959 2 21596 1959.1260 316.49 315.86 316.27 \n",
"1959-03-01 14 1959 3 21624 1959.2027 316.65 315.38 316.98 \n",
"1959-04-01 15 1959 4 21655 1959.2877 317.72 315.42 318.09 \n",
"1959-05-01 16 1959 5 21685 1959.3699 318.29 315.49 318.65 \n",
"1959-06-01 17 1959 6 21716 1959.4548 318.15 316.03 318.04 \n",
"1959-07-01 18 1959 7 21746 1959.5370 316.54 315.86 316.67 \n",
"\n",
" s4 s5 s6 \n",
"period \n",
......@@ -889,10 +1099,20 @@
"1958-04-01 314.99 317.45 315.16 \n",
"1958-05-01 315.06 317.51 314.71 \n",
"1958-07-01 315.22 315.86 315.19 \n",
"1958-08-01 315.29 314.93 316.19 "
"1958-08-01 315.29 314.93 316.19 \n",
"1958-09-01 315.35 313.21 316.08 \n",
"1958-11-01 315.46 313.33 315.20 \n",
"1958-12-01 315.51 314.67 315.43 \n",
"1959-01-01 315.57 315.58 315.54 \n",
"1959-02-01 315.63 316.49 315.86 \n",
"1959-03-01 315.69 316.65 315.38 \n",
"1959-04-01 315.77 317.72 315.42 \n",
"1959-05-01 315.85 318.29 315.49 \n",
"1959-06-01 315.94 318.15 316.03 \n",
"1959-07-01 316.03 316.54 315.86 "
]
},
"execution_count": 8,
"execution_count": 90,
"metadata": {},
"output_type": "execute_result"
}
......@@ -906,7 +1126,7 @@
"#print(period.shape)\n",
"df = pd.concat([df,period],axis=1)\n",
"df = df.set_index('period') \n",
"df.head()"
"df[0:15]"
]
},
{
......@@ -918,7 +1138,7 @@
},
{
"cell_type": "code",
"execution_count": 9,
"execution_count": 63,
"metadata": {},
"outputs": [
{
......@@ -949,7 +1169,7 @@
},
{
"cell_type": "code",
"execution_count": 10,
"execution_count": 64,
"metadata": {},
"outputs": [
{
......@@ -1001,7 +1221,7 @@
},
{
"cell_type": "code",
"execution_count": 11,
"execution_count": 65,
"metadata": {},
"outputs": [
{
......@@ -1058,7 +1278,7 @@
},
{
"cell_type": "code",
"execution_count": 111,
"execution_count": 66,
"metadata": {},
"outputs": [
{
......@@ -1202,7 +1422,7 @@
"1964-04-01 321.83 319.45 "
]
},
"execution_count": 111,
"execution_count": 66,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1218,7 +1438,7 @@
},
{
"cell_type": "code",
"execution_count": 112,
"execution_count": 67,
"metadata": {},
"outputs": [
{
......@@ -1362,7 +1582,7 @@
"2020-01-01 413.05 413.37 413.33 "
]
},
"execution_count": 112,
"execution_count": 67,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1380,7 +1600,7 @@
},
{
"cell_type": "code",
"execution_count": 113,
"execution_count": 68,
"metadata": {},
"outputs": [
{
......@@ -1503,6 +1723,34 @@
" <td>315.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-08-01</th>\n",
" <td>7</td>\n",
" <td>1958</td>\n",
" <td>8</td>\n",
" <td>21412</td>\n",
" <td>1958.6219</td>\n",
" <td>314.930</td>\n",
" <td>316.19</td>\n",
" <td>314.00</td>\n",
" <td>315.29</td>\n",
" <td>314.93</td>\n",
" <td>316.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-09-01</th>\n",
" <td>8</td>\n",
" <td>1958</td>\n",
" <td>9</td>\n",
" <td>21443</td>\n",
" <td>1958.7068</td>\n",
" <td>313.210</td>\n",
" <td>316.08</td>\n",
" <td>312.46</td>\n",
" <td>315.35</td>\n",
" <td>313.21</td>\n",
" <td>316.08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
......@@ -1517,6 +1765,34 @@
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-07-01</th>\n",
" <td>738</td>\n",
" <td>2019</td>\n",
" <td>7</td>\n",
" <td>43661</td>\n",
" <td>2019.5370</td>\n",
" <td>411.850</td>\n",
" <td>411.03</td>\n",
" <td>412.36</td>\n",
" <td>411.57</td>\n",
" <td>411.85</td>\n",
" <td>411.03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-08-01</th>\n",
" <td>739</td>\n",
" <td>2019</td>\n",
" <td>8</td>\n",
" <td>43692</td>\n",
" <td>2019.6219</td>\n",
" <td>410.080</td>\n",
" <td>411.62</td>\n",
" <td>410.23</td>\n",
" <td>411.81</td>\n",
" <td>410.08</td>\n",
" <td>411.62</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019-09-01</th>\n",
" <td>740</td>\n",
" <td>2019</td>\n",
......@@ -1599,7 +1875,11 @@
"1958-05-01 4 1958 5 21320 1958.3699 317.510 314.71 317.86 \n",
"1958-06-01 5 1958 6 21351 1958.4548 316.685 314.95 317.24 \n",
"1958-07-01 6 1958 7 21381 1958.5370 315.860 315.19 315.86 \n",
"1958-08-01 7 1958 8 21412 1958.6219 314.930 316.19 314.00 \n",
"1958-09-01 8 1958 9 21443 1958.7068 313.210 316.08 312.46 \n",
"... ... ... .. ... ... ... ... ... \n",
"2019-07-01 738 2019 7 43661 2019.5370 411.850 411.03 412.36 \n",
"2019-08-01 739 2019 8 43692 2019.6219 410.080 411.62 410.23 \n",
"2019-09-01 740 2019 9 43723 2019.7068 408.550 412.06 408.52 \n",
"2019-10-01 741 2019 10 43753 2019.7890 408.430 412.07 408.67 \n",
"2019-11-01 742 2019 11 43784 2019.8740 410.280 412.56 410.29 \n",
......@@ -1613,7 +1893,11 @@
"1958-05-01 315.06 317.51 314.71 \n",
"1958-06-01 315.14 317.24 315.14 \n",
"1958-07-01 315.22 315.86 315.19 \n",
"1958-08-01 315.29 314.93 316.19 \n",
"1958-09-01 315.35 313.21 316.08 \n",
"... ... ... ... \n",
"2019-07-01 411.57 411.85 411.03 \n",
"2019-08-01 411.81 410.08 411.62 \n",
"2019-09-01 412.05 408.55 412.06 \n",
"2019-10-01 412.29 408.43 412.07 \n",
"2019-11-01 412.54 410.28 412.56 \n",
......@@ -1623,7 +1907,7 @@
"[743 rows x 11 columns]"
]
},
"execution_count": 113,
"execution_count": 68,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1635,7 +1919,7 @@
},
{
"cell_type": "code",
"execution_count": 114,
"execution_count": 69,
"metadata": {},
"outputs": [
{
......@@ -1758,18 +2042,32 @@
" <td>319.40</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <th>1964-06-01</th>\n",
" <td>77</td>\n",
" <td>1964</td>\n",
" <td>6</td>\n",
" <td>23543</td>\n",
" <td>1964.4563</td>\n",
" <td>321.8900</td>\n",
" <td>319.7500</td>\n",
" <td>321.64</td>\n",
" <td>319.52</td>\n",
" <td>321.89</td>\n",
" <td>319.75</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-07-01</th>\n",
" <td>78</td>\n",
" <td>1964</td>\n",
" <td>7</td>\n",
" <td>23573</td>\n",
" <td>1964.5383</td>\n",
" <td>320.4400</td>\n",
" <td>319.7800</td>\n",
" <td>320.19</td>\n",
" <td>319.55</td>\n",
" <td>320.44</td>\n",
" <td>319.78</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964-08-01</th>\n",
......@@ -1843,7 +2141,6 @@
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>12 rows × 11 columns</p>\n",
"</div>"
],
"text/plain": [
......@@ -1854,7 +2151,8 @@
"1964-03-01 74 1964 3 23451 1964.2049 320.9150 319.4650 320.74 \n",
"1964-04-01 75 1964 4 23482 1964.2896 321.5875 319.4325 321.83 \n",
"1964-05-01 76 1964 5 23512 1964.3716 322.2600 319.4000 322.34 \n",
"... ... ... .. ... ... ... ... ... \n",
"1964-06-01 77 1964 6 23543 1964.4563 321.8900 319.7500 321.64 \n",
"1964-07-01 78 1964 7 23573 1964.5383 320.4400 319.7800 320.19 \n",
"1964-08-01 79 1964 8 23604 1964.6230 318.6900 320.0000 318.24 \n",
"1964-09-01 80 1964 9 23635 1964.7077 316.7100 319.6500 316.64 \n",
"1964-10-01 81 1964 10 23665 1964.7896 316.8700 319.9100 316.59 \n",
......@@ -1868,17 +2166,16 @@
"1964-03-01 319.41 320.74 319.41 \n",
"1964-04-01 319.45 321.83 319.45 \n",
"1964-05-01 319.49 322.26 319.40 \n",
"... ... ... ... \n",
"1964-06-01 319.52 321.89 319.75 \n",
"1964-07-01 319.55 320.44 319.78 \n",
"1964-08-01 319.58 318.69 320.00 \n",
"1964-09-01 319.60 316.71 319.65 \n",
"1964-10-01 319.62 316.87 319.91 \n",
"1964-11-01 319.64 317.68 319.58 \n",
"1964-12-01 319.66 318.71 319.49 \n",
"\n",
"[12 rows x 11 columns]"
"1964-12-01 319.66 318.71 319.49 "
]
},
"execution_count": 114,
"execution_count": 69,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1890,7 +2187,7 @@
},
{
"cell_type": "code",
"execution_count": 115,
"execution_count": 70,
"metadata": {},
"outputs": [],
"source": [
......@@ -1915,7 +2212,7 @@
},
{
"cell_type": "code",
"execution_count": 116,
"execution_count": 71,
"metadata": {},
"outputs": [
{
......@@ -1965,7 +2262,7 @@
},
{
"cell_type": "code",
"execution_count": 117,
"execution_count": 72,
"metadata": {},
"outputs": [
{
......@@ -1988,7 +2285,7 @@
},
{
"cell_type": "code",
"execution_count": 118,
"execution_count": 96,
"metadata": {
"scrolled": true
},
......@@ -2013,22 +2310,23 @@
"\n",
"p0 = [f0,ampl0,phase0,offset0]\n",
"x = np.array(data_cube['Date 2'])\n",
"popt, pcov = curve_fit(my_sin, x , s1, p0 = p0)\n",
"print(popt)\n",
"ylent = my_sin(x,popt[0],popt[1],popt[2],popt[3])\n",
"popt_l, pcov_l = curve_fit(my_sin, x , s1, p0 = p0)\n",
"print(popt_l)\n",
"ylent = my_sin(x,popt_l[0],popt_l[1],popt_l[2],popt_l[3])\n",
"fl = p0[0]\n",
"\n",
"f0 = f_rapide \n",
"ampl0= 3.5\n",
"p0 = [f0,ampl0,phase0,offset0]\n",
"popt, pcov = curve_fit(my_sin, x , s1, p0 = p0)\n",
"print(popt)\n",
"yrapide = my_sin(x,popt[0],popt[1],popt[2],popt[3]) \n",
" "
"popt_r, pcov_r = curve_fit(my_sin, x , s1, p0 = p0)\n",
"print(popt_r)\n",
"yrapide = my_sin(x,popt_r[0],popt_r[1],popt_r[2],popt_r[3]) \n",
"fr = p0[0]"
]
},
{
"cell_type": "code",
"execution_count": 119,
"execution_count": 97,
"metadata": {
"scrolled": true
},
......@@ -2036,11 +2334,11 @@
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f8e88a52898>,\n",
" <matplotlib.lines.Line2D at 0x7f8e889fbef0>]"
"[<matplotlib.lines.Line2D at 0x7fbbaa92a9b0>,\n",
" <matplotlib.lines.Line2D at 0x7fbbaa947c88>]"
]
},
"execution_count": 119,
"execution_count": 97,
"metadata": {},
"output_type": "execute_result"
},
......@@ -2069,6 +2367,924 @@
"ax2.set(title = 'Variation concentration CO2 rapide',xlabel='Période',ylabel='Concentration (ppm)')\n",
"ax2.plot(x,yrapide,x,s1)"
]
},
{
"cell_type": "code",
"execution_count": 75,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" Période lente 30.958333333333332\n",
" Période rapide 0.9986559139784947\n"
]
}
],
"source": [
"print(f' Période lente {1/fl}\\n Période rapide {1/fr}')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nous avons pour la variation lente une fréquence de 0.0294 Hz et pour la variation rapide de 1.0004 Hz. Sachant que 1 Hz \n",
"représente 1 an, nous avons un peu plus d'un an pour la variation rapide, et pour la variation lente une période d'environ 34 ans."
]
},
{
"cell_type": "code",
"execution_count": 76,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"34.013605442176875"
]
},
"execution_count": 76,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"1/0.0294"
]
},
{
"cell_type": "code",
"execution_count": 92,
"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>Mn</th>\n",
" <th>1</th>\n",
" <th>2</th>\n",
" <th>3</th>\n",
" <th>4</th>\n",
" <th>5</th>\n",
" <th>6</th>\n",
" <th>7</th>\n",
" <th>8</th>\n",
" <th>9</th>\n",
" <th>10</th>\n",
" <th>11</th>\n",
" <th>12</th>\n",
" </tr>\n",
" <tr>\n",
" <th>Yr</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1958</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>315.70</td>\n",
" <td>317.45</td>\n",
" <td>317.51</td>\n",
" <td>NaN</td>\n",
" <td>315.86</td>\n",
" <td>314.93</td>\n",
" <td>313.21</td>\n",
" <td>NaN</td>\n",
" <td>313.33</td>\n",
" <td>314.67</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959</th>\n",
" <td>315.58</td>\n",
" <td>316.49</td>\n",
" <td>316.65</td>\n",
" <td>317.72</td>\n",
" <td>318.29</td>\n",
" <td>318.15</td>\n",
" <td>316.54</td>\n",
" <td>314.80</td>\n",
" <td>313.84</td>\n",
" <td>313.33</td>\n",
" <td>314.81</td>\n",
" <td>315.58</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1960</th>\n",
" <td>316.43</td>\n",
" <td>316.98</td>\n",
" <td>317.58</td>\n",
" <td>319.03</td>\n",
" <td>320.04</td>\n",
" <td>319.59</td>\n",
" <td>318.18</td>\n",
" <td>315.90</td>\n",
" <td>314.17</td>\n",
" <td>313.83</td>\n",
" <td>315.00</td>\n",
" <td>316.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1961</th>\n",
" <td>316.89</td>\n",
" <td>317.70</td>\n",
" <td>318.54</td>\n",
" <td>319.48</td>\n",
" <td>320.58</td>\n",
" <td>319.77</td>\n",
" <td>318.57</td>\n",
" <td>316.79</td>\n",
" <td>314.99</td>\n",
" <td>315.31</td>\n",
" <td>316.10</td>\n",
" <td>317.01</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1962</th>\n",
" <td>317.94</td>\n",
" <td>318.55</td>\n",
" <td>319.68</td>\n",
" <td>320.57</td>\n",
" <td>321.02</td>\n",
" <td>320.62</td>\n",
" <td>319.61</td>\n",
" <td>317.40</td>\n",
" <td>316.25</td>\n",
" <td>315.42</td>\n",
" <td>316.69</td>\n",
" <td>317.70</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1963</th>\n",
" <td>318.74</td>\n",
" <td>319.07</td>\n",
" <td>319.86</td>\n",
" <td>321.38</td>\n",
" <td>322.25</td>\n",
" <td>321.48</td>\n",
" <td>319.74</td>\n",
" <td>317.77</td>\n",
" <td>316.21</td>\n",
" <td>315.99</td>\n",
" <td>317.07</td>\n",
" <td>318.35</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1964</th>\n",
" <td>319.57</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>322.26</td>\n",
" <td>321.89</td>\n",
" <td>320.44</td>\n",
" <td>318.69</td>\n",
" <td>316.71</td>\n",
" <td>316.87</td>\n",
" <td>317.68</td>\n",
" <td>318.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2014</th>\n",
" <td>397.93</td>\n",
" <td>398.10</td>\n",
" <td>399.47</td>\n",
" <td>401.33</td>\n",
" <td>401.88</td>\n",
" <td>401.31</td>\n",
" <td>399.07</td>\n",
" <td>397.21</td>\n",
" <td>395.40</td>\n",
" <td>395.65</td>\n",
" <td>397.23</td>\n",
" <td>398.79</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2015</th>\n",
" <td>399.85</td>\n",
" <td>400.31</td>\n",
" <td>401.51</td>\n",
" <td>403.45</td>\n",
" <td>404.10</td>\n",
" <td>402.88</td>\n",
" <td>401.61</td>\n",
" <td>399.00</td>\n",
" <td>397.50</td>\n",
" <td>398.28</td>\n",
" <td>400.24</td>\n",
" <td>401.89</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2016</th>\n",
" <td>402.65</td>\n",
" <td>404.16</td>\n",
" <td>404.85</td>\n",
" <td>407.57</td>\n",
" <td>407.66</td>\n",
" <td>407.00</td>\n",
" <td>404.50</td>\n",
" <td>402.24</td>\n",
" <td>401.01</td>\n",
" <td>401.50</td>\n",
" <td>403.64</td>\n",
" <td>404.55</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2017</th>\n",
" <td>406.07</td>\n",
" <td>406.64</td>\n",
" <td>407.06</td>\n",
" <td>408.95</td>\n",
" <td>409.91</td>\n",
" <td>409.12</td>\n",
" <td>407.20</td>\n",
" <td>405.24</td>\n",
" <td>403.27</td>\n",
" <td>403.64</td>\n",
" <td>405.17</td>\n",
" <td>406.75</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2018</th>\n",
" <td>408.05</td>\n",
" <td>408.34</td>\n",
" <td>409.25</td>\n",
" <td>410.30</td>\n",
" <td>411.30</td>\n",
" <td>410.88</td>\n",
" <td>408.90</td>\n",
" <td>407.10</td>\n",
" <td>405.59</td>\n",
" <td>405.99</td>\n",
" <td>408.12</td>\n",
" <td>409.23</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2019</th>\n",
" <td>410.92</td>\n",
" <td>411.66</td>\n",
" <td>412.00</td>\n",
" <td>413.52</td>\n",
" <td>414.83</td>\n",
" <td>413.96</td>\n",
" <td>411.85</td>\n",
" <td>410.08</td>\n",
" <td>408.55</td>\n",
" <td>408.43</td>\n",
" <td>410.28</td>\n",
" <td>411.85</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2020</th>\n",
" <td>413.37</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>63 rows × 12 columns</p>\n",
"</div>"
],
"text/plain": [
"Mn 1 2 3 4 5 6 7 8 9 \\\n",
"Yr \n",
"1958 NaN NaN 315.70 317.45 317.51 NaN 315.86 314.93 313.21 \n",
"1959 315.58 316.49 316.65 317.72 318.29 318.15 316.54 314.80 313.84 \n",
"1960 316.43 316.98 317.58 319.03 320.04 319.59 318.18 315.90 314.17 \n",
"1961 316.89 317.70 318.54 319.48 320.58 319.77 318.57 316.79 314.99 \n",
"1962 317.94 318.55 319.68 320.57 321.02 320.62 319.61 317.40 316.25 \n",
"1963 318.74 319.07 319.86 321.38 322.25 321.48 319.74 317.77 316.21 \n",
"1964 319.57 NaN NaN NaN 322.26 321.89 320.44 318.69 316.71 \n",
"... ... ... ... ... ... ... ... ... ... \n",
"2014 397.93 398.10 399.47 401.33 401.88 401.31 399.07 397.21 395.40 \n",
"2015 399.85 400.31 401.51 403.45 404.10 402.88 401.61 399.00 397.50 \n",
"2016 402.65 404.16 404.85 407.57 407.66 407.00 404.50 402.24 401.01 \n",
"2017 406.07 406.64 407.06 408.95 409.91 409.12 407.20 405.24 403.27 \n",
"2018 408.05 408.34 409.25 410.30 411.30 410.88 408.90 407.10 405.59 \n",
"2019 410.92 411.66 412.00 413.52 414.83 413.96 411.85 410.08 408.55 \n",
"2020 413.37 NaN NaN NaN NaN NaN NaN NaN NaN \n",
"\n",
"Mn 10 11 12 \n",
"Yr \n",
"1958 NaN 313.33 314.67 \n",
"1959 313.33 314.81 315.58 \n",
"1960 313.83 315.00 316.19 \n",
"1961 315.31 316.10 317.01 \n",
"1962 315.42 316.69 317.70 \n",
"1963 315.99 317.07 318.35 \n",
"1964 316.87 317.68 318.71 \n",
"... ... ... ... \n",
"2014 395.65 397.23 398.79 \n",
"2015 398.28 400.24 401.89 \n",
"2016 401.50 403.64 404.55 \n",
"2017 403.64 405.17 406.75 \n",
"2018 405.99 408.12 409.23 \n",
"2019 408.43 410.28 411.85 \n",
"2020 NaN NaN NaN \n",
"\n",
"[63 rows x 12 columns]"
]
},
"execution_count": 92,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df1 = pd.crosstab(df['Yr'],df['Mn'],values=df['s1'],aggfunc=pd.Series)\n",
"df1"
]
},
{
"cell_type": "code",
"execution_count": 95,
"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>Mn</th>\n",
" <th>1</th>\n",
" <th>2</th>\n",
" <th>3</th>\n",
" <th>4</th>\n",
" <th>5</th>\n",
" <th>6</th>\n",
" <th>7</th>\n",
" <th>8</th>\n",
" <th>9</th>\n",
" <th>10</th>\n",
" <th>11</th>\n",
" <th>12</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>count</th>\n",
" <td>62.000000</td>\n",
" <td>60.000000</td>\n",
" <td>61.000000</td>\n",
" <td>61.000000</td>\n",
" <td>62.000000</td>\n",
" <td>61.000000</td>\n",
" <td>62.000000</td>\n",
" <td>62.000000</td>\n",
" <td>62.000000</td>\n",
" <td>61.000000</td>\n",
" <td>62.000000</td>\n",
" <td>62.000000</td>\n",
" </tr>\n",
" <tr>\n",
" <th>mean</th>\n",
" <td>355.559355</td>\n",
" <td>355.965000</td>\n",
" <td>356.105902</td>\n",
" <td>357.471475</td>\n",
" <td>357.509355</td>\n",
" <td>357.575410</td>\n",
" <td>355.428871</td>\n",
" <td>353.441129</td>\n",
" <td>351.778871</td>\n",
" <td>352.420656</td>\n",
" <td>353.175645</td>\n",
" <td>354.473387</td>\n",
" </tr>\n",
" <tr>\n",
" <th>std</th>\n",
" <td>28.938877</td>\n",
" <td>28.066464</td>\n",
" <td>28.301402</td>\n",
" <td>28.452242</td>\n",
" <td>28.657777</td>\n",
" <td>28.333317</td>\n",
" <td>28.395757</td>\n",
" <td>28.315031</td>\n",
" <td>28.340139</td>\n",
" <td>28.284758</td>\n",
" <td>28.700183</td>\n",
" <td>28.819885</td>\n",
" </tr>\n",
" <tr>\n",
" <th>min</th>\n",
" <td>315.580000</td>\n",
" <td>316.490000</td>\n",
" <td>315.700000</td>\n",
" <td>317.450000</td>\n",
" <td>317.510000</td>\n",
" <td>318.150000</td>\n",
" <td>315.860000</td>\n",
" <td>314.800000</td>\n",
" <td>313.210000</td>\n",
" <td>313.330000</td>\n",
" <td>313.330000</td>\n",
" <td>314.670000</td>\n",
" </tr>\n",
" <tr>\n",
" <th>25%</th>\n",
" <td>329.612500</td>\n",
" <td>331.227500</td>\n",
" <td>331.480000</td>\n",
" <td>332.650000</td>\n",
" <td>332.632500</td>\n",
" <td>332.250000</td>\n",
" <td>330.947500</td>\n",
" <td>329.330000</td>\n",
" <td>327.450000</td>\n",
" <td>327.370000</td>\n",
" <td>328.235000</td>\n",
" <td>328.872500</td>\n",
" </tr>\n",
" <tr>\n",
" <th>50%</th>\n",
" <td>353.210000</td>\n",
" <td>353.885000</td>\n",
" <td>353.680000</td>\n",
" <td>355.420000</td>\n",
" <td>354.945000</td>\n",
" <td>355.120000</td>\n",
" <td>353.140000</td>\n",
" <td>351.050000</td>\n",
" <td>349.270000</td>\n",
" <td>349.990000</td>\n",
" <td>350.685000</td>\n",
" <td>351.930000</td>\n",
" </tr>\n",
" <tr>\n",
" <th>75%</th>\n",
" <td>377.952500</td>\n",
" <td>377.930000</td>\n",
" <td>378.390000</td>\n",
" <td>380.500000</td>\n",
" <td>380.052500</td>\n",
" <td>379.550000</td>\n",
" <td>377.470000</td>\n",
" <td>375.492500</td>\n",
" <td>373.782500</td>\n",
" <td>374.220000</td>\n",
" <td>375.467500</td>\n",
" <td>377.002500</td>\n",
" </tr>\n",
" <tr>\n",
" <th>max</th>\n",
" <td>413.370000</td>\n",
" <td>411.660000</td>\n",
" <td>412.000000</td>\n",
" <td>413.520000</td>\n",
" <td>414.830000</td>\n",
" <td>413.960000</td>\n",
" <td>411.850000</td>\n",
" <td>410.080000</td>\n",
" <td>408.550000</td>\n",
" <td>408.430000</td>\n",
" <td>410.280000</td>\n",
" <td>411.850000</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
"Mn 1 2 3 4 5 6 \\\n",
"count 62.000000 60.000000 61.000000 61.000000 62.000000 61.000000 \n",
"mean 355.559355 355.965000 356.105902 357.471475 357.509355 357.575410 \n",
"std 28.938877 28.066464 28.301402 28.452242 28.657777 28.333317 \n",
"min 315.580000 316.490000 315.700000 317.450000 317.510000 318.150000 \n",
"25% 329.612500 331.227500 331.480000 332.650000 332.632500 332.250000 \n",
"50% 353.210000 353.885000 353.680000 355.420000 354.945000 355.120000 \n",
"75% 377.952500 377.930000 378.390000 380.500000 380.052500 379.550000 \n",
"max 413.370000 411.660000 412.000000 413.520000 414.830000 413.960000 \n",
"\n",
"Mn 7 8 9 10 11 12 \n",
"count 62.000000 62.000000 62.000000 61.000000 62.000000 62.000000 \n",
"mean 355.428871 353.441129 351.778871 352.420656 353.175645 354.473387 \n",
"std 28.395757 28.315031 28.340139 28.284758 28.700183 28.819885 \n",
"min 315.860000 314.800000 313.210000 313.330000 313.330000 314.670000 \n",
"25% 330.947500 329.330000 327.450000 327.370000 328.235000 328.872500 \n",
"50% 353.140000 351.050000 349.270000 349.990000 350.685000 351.930000 \n",
"75% 377.470000 375.492500 373.782500 374.220000 375.467500 377.002500 \n",
"max 411.850000 410.080000 408.550000 408.430000 410.280000 411.850000 "
]
},
"execution_count": 95,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df1.describe()"
]
},
{
"cell_type": "code",
"execution_count": 133,
"metadata": {
"scrolled": true
},
"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>Mn</th>\n",
" <th>1</th>\n",
" <th>2</th>\n",
" <th>3</th>\n",
" <th>4</th>\n",
" <th>5</th>\n",
" <th>6</th>\n",
" <th>7</th>\n",
" <th>8</th>\n",
" <th>9</th>\n",
" <th>10</th>\n",
" <th>11</th>\n",
" <th>12</th>\n",
" </tr>\n",
" <tr>\n",
" <th>Yr</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1958</th>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>0.2027</td>\n",
" <td>0.2877</td>\n",
" <td>0.3699</td>\n",
" <td>NaN</td>\n",
" <td>0.5370</td>\n",
" <td>0.6219</td>\n",
" <td>0.7068</td>\n",
" <td>NaN</td>\n",
" <td>0.8740</td>\n",
" <td>0.9562</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1959</th>\n",
" <td>0.0411</td>\n",
" <td>0.1260</td>\n",
" <td>0.2027</td>\n",
" <td>0.2877</td>\n",
" <td>0.3699</td>\n",
" <td>0.4548</td>\n",
" <td>0.5370</td>\n",
" <td>0.6219</td>\n",
" <td>0.7068</td>\n",
" <td>0.7890</td>\n",
" <td>0.8740</td>\n",
" <td>0.9562</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1960</th>\n",
" <td>0.0410</td>\n",
" <td>0.1257</td>\n",
" <td>0.2049</td>\n",
" <td>0.2896</td>\n",
" <td>0.3716</td>\n",
" <td>0.4563</td>\n",
" <td>0.5383</td>\n",
" <td>0.6230</td>\n",
" <td>0.7077</td>\n",
" <td>0.7896</td>\n",
" <td>0.8743</td>\n",
" <td>0.9563</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1961</th>\n",
" <td>0.0411</td>\n",
" <td>0.1260</td>\n",
" <td>0.2027</td>\n",
" <td>0.2877</td>\n",
" <td>0.3699</td>\n",
" <td>0.4548</td>\n",
" <td>0.5370</td>\n",
" <td>0.6219</td>\n",
" <td>0.7068</td>\n",
" <td>0.7890</td>\n",
" <td>0.8740</td>\n",
" <td>0.9562</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1962</th>\n",
" <td>0.0411</td>\n",
" <td>0.1260</td>\n",
" <td>0.2027</td>\n",
" <td>0.2877</td>\n",
" <td>0.3699</td>\n",
" <td>0.4548</td>\n",
" <td>0.5370</td>\n",
" <td>0.6219</td>\n",
" <td>0.7068</td>\n",
" <td>0.7890</td>\n",
" <td>0.8740</td>\n",
" <td>0.9562</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
"Mn 1 2 3 4 5 6 7 8 9 \\\n",
"Yr \n",
"1958 NaN NaN 0.2027 0.2877 0.3699 NaN 0.5370 0.6219 0.7068 \n",
"1959 0.0411 0.1260 0.2027 0.2877 0.3699 0.4548 0.5370 0.6219 0.7068 \n",
"1960 0.0410 0.1257 0.2049 0.2896 0.3716 0.4563 0.5383 0.6230 0.7077 \n",
"1961 0.0411 0.1260 0.2027 0.2877 0.3699 0.4548 0.5370 0.6219 0.7068 \n",
"1962 0.0411 0.1260 0.2027 0.2877 0.3699 0.4548 0.5370 0.6219 0.7068 \n",
"\n",
"Mn 10 11 12 \n",
"Yr \n",
"1958 NaN 0.8740 0.9562 \n",
"1959 0.7890 0.8740 0.9562 \n",
"1960 0.7896 0.8743 0.9563 \n",
"1961 0.7890 0.8740 0.9562 \n",
"1962 0.7890 0.8740 0.9562 "
]
},
"execution_count": 133,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"m = df['Date 2'] - df['Yr']\n",
"df2 = pd.crosstab(df['Yr'],df['Mn'],values= m,aggfunc=pd.Series)\n",
"df2.head(5)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"On observe pour certaine année un pas différent à partir du mois de février cela est du aux années bissextiles. Dans l'intervalle d'extrapolation, nous avons l'année 2020 et 2024 en année bissextile."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nous créons une nouvelle DataFrame pour l'interpolation, elle ne contient que les colonnes nécessaire à l'extrapolation."
]
},
{
"cell_type": "code",
"execution_count": 126,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"period\n",
"1959-01-01 0.0411\n",
"1959-02-01 0.1260\n",
"1959-03-01 0.2027\n",
"1959-04-01 0.2877\n",
"1959-05-01 0.3699\n",
"1959-06-01 0.4548\n",
"1959-07-01 0.5370\n",
"1959-08-01 0.6219\n",
"1959-09-01 0.7068\n",
"1959-10-01 0.7890\n",
"1959-11-01 0.8740\n",
"1959-12-01 0.9562\n",
"dtype: float64\n",
"period\n",
"1960-01-01 0.0410\n",
"1960-02-01 0.1257\n",
"1960-03-01 0.2049\n",
"1960-04-01 0.2896\n",
"1960-05-01 0.3716\n",
"1960-06-01 0.4563\n",
"1960-07-01 0.5383\n",
"1960-08-01 0.6230\n",
"1960-09-01 0.7077\n",
"1960-10-01 0.7896\n",
"1960-11-01 0.8743\n",
"1960-12-01 0.9563\n",
"dtype: float64\n"
]
}
],
"source": []
},
{
"cell_type": "code",
"execution_count": 101,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.1246"
]
},
"execution_count": 101,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"0.0416+0.083"
]
},
{
"cell_type": "code",
"execution_count": 103,
"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>index</th>\n",
" <th>Yr</th>\n",
" <th>Mn</th>\n",
" <th>Date 1</th>\n",
" <th>Date 2</th>\n",
" <th>s1</th>\n",
" <th>s2</th>\n",
" <th>s3</th>\n",
" <th>s4</th>\n",
" <th>s5</th>\n",
" <th>s6</th>\n",
" </tr>\n",
" <tr>\n",
" <th>period</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>2020-01-01</th>\n",
" <td>744</td>\n",
" <td>2020</td>\n",
" <td>1</td>\n",
" <td>43845</td>\n",
" <td>2020.041</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" <td>413.11</td>\n",
" <td>413.05</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" index Yr Mn Date 1 Date 2 s1 s2 s3 s4 \\\n",
"period \n",
"2020-01-01 744 2020 1 43845 2020.041 413.37 413.33 413.11 413.05 \n",
"\n",
" s5 s6 \n",
"period \n",
"2020-01-01 413.37 413.33 "
]
},
"execution_count": 103,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df[-1:]"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment