5 "execution_count": null,
16 "execution_count": null,
22 "data = getData(ts_data=\"../data/pm10_mesures_H_loc.csv\", exo_data=\"../data/meteo_extra_noNAs.csv\",\n",
23 " input_tz = \"Europe/Paris\", working_tz=\"Europe/Paris\", predict_at=7)"
27 "cell_type": "markdown",
30 "## Pollution par chauffage"
35 "execution_count": null,
41 "p_ch_nn = getForecast(data, seq(as.Date(\"2015-01-18\"),as.Date(\"2015-01-24\"),\"days\"), Inf, 17,\n",
42 " \"Neighbors\", \"Neighbors\", simtype=\"mix\", same_season=FALSE, mix_strategy=\"mult\")\n",
43 "p_ch_pz = getForecast(data, seq(as.Date(\"2015-01-18\"),as.Date(\"2015-01-24\"),\"days\"), Inf, 17,\n",
44 " \"Persistence\", \"Zero\")\n",
45 "p_ch_az = getForecast(data, seq(as.Date(\"2015-01-18\"),as.Date(\"2015-01-24\"),\"days\"), Inf, 17,\n",
46 " \"Average\", \"Zero\")"
51 "execution_count": null,
57 "e_ch_nn = getError(data, p_ch_nn, 17)\n",
58 "e_ch_pz = getError(data, p_ch_pz, 17)\n",
59 "e_ch_az = getError(data, p_ch_az, 17)\n",
60 "options(repr.plot.width=9, repr.plot.height=6)\n",
61 "plotError(list(e_ch_nn, e_ch_pz, e_ch_az), cols=c(1,2,colors()[258]))"
66 "execution_count": null,
72 "par(mfrow=c(1,2))\n",
73 "options(repr.plot.width=9, repr.plot.height=4)\n",
74 "plotPredReal(data, p_ch_nn, 3)\n",
75 "plotPredReal(data, p_ch_nn, 4)"
80 "execution_count": null,
86 "par(mfrow=c(1,2))\n",
87 "plotFilaments(data, p_ch_nn$getIndexInData(3))\n",
88 "plotFilaments(data, p_ch_nn$getIndexInData(4))"
93 "execution_count": null,
99 "par(mfrow=c(1,3))\n",
100 "plotSimils(p_ch_nn, 3)\n",
101 "plotSimils(p_ch_nn, 4)\n",
102 "plotSimils(p_ch_nn, 5)"
106 "cell_type": "markdown",
109 "## Pollution par épandage"
114 "execution_count": null,
120 "p_ep_nn = getForecast(data, seq(as.Date(\"2015-03-15\"),as.Date(\"2015-03-21\"),\"days\"), Inf, 17,\n",
121 " \"Neighbors\", \"Neighbors\", simtype=\"mix\", same_season=FALSE, mix_strategy=\"mult\")\n",
122 "p_ep_pz = getForecast(data, seq(as.Date(\"2015-03-15\"),as.Date(\"2015-03-21\"),\"days\"), Inf, 17,\n",
123 " \"Persistence\", \"Zero\")\n",
124 "p_ep_az = getForecast(data, seq(as.Date(\"2015-03-15\"),as.Date(\"2015-03-21\"),\"days\"), Inf, 17,\n",
125 " \"Average\", \"Zero\")"
130 "execution_count": null,
136 "e_ep_nn = getError(data, p_ep_nn, 17)\n",
137 "e_ep_pz = getError(data, p_ep_pz, 17)\n",
138 "e_ep_az = getError(data, p_ep_az, 17)\n",
139 "options(repr.plot.width=9, repr.plot.height=6)\n",
140 "plotError(list(e_ep_nn, e_ep_pz, e_ep_az), cols=c(1,2,colors()[258]))"
145 "execution_count": null,
151 "par(mfrow=c(1,2))\n",
152 "options(repr.plot.width=9, repr.plot.height=4)\n",
153 "plotPredReal(data, p_ep_nn, 3)\n",
154 "plotPredReal(data, p_ep_nn, 4)"
159 "execution_count": null,
165 "par(mfrow=c(1,2))\n",
166 "plotFilaments(data, p_ep_nn$getIndexInData(3))\n",
167 "plotFilaments(data, p_ep_nn$getIndexInData(4))"
172 "execution_count": null,
178 "par(mfrow=c(1,3))\n",
179 "plotSimils(p_ep_nn, 3)\n",
180 "plotSimils(p_ep_nn, 4)\n",
181 "plotSimils(p_ep_nn, 5)"
185 "cell_type": "markdown",
188 "## Semaine non polluée"
193 "execution_count": null,
199 "p_np_nn = getForecast(data, seq(as.Date(\"2015-04-26\"),as.Date(\"2015-05-02\"),\"days\"), Inf, 17,\n",
200 " \"Neighbors\", \"Neighbors\", simtype=\"mix\", same_season=FALSE, mix_strategy=\"mult\")\n",
201 "p_np_pz = getForecast(data, seq(as.Date(\"2015-04-26\"),as.Date(\"2015-05-02\"),\"days\"), Inf, 17,\n",
202 " \"Persistence\", \"Zero\")\n",
203 "p_np_az = getForecast(data, seq(as.Date(\"2015-04-26\"),as.Date(\"2015-05-02\"),\"days\"), Inf, 17,\n",
204 " \"Average\", \"Zero\")"
209 "execution_count": null,
215 "e_np_nn = getError(data, p_np_nn, 17)\n",
216 "e_np_pz = getError(data, p_np_pz, 17)\n",
217 "e_np_az = getError(data, p_np_az, 17)\n",
218 "options(repr.plot.width=9, repr.plot.height=6)\n",
219 "plotError(list(e_np_nn, e_np_pz, e_np_az), cols=c(1,2,colors()[258]))"
224 "execution_count": null,
230 "par(mfrow=c(1,2))\n",
231 "options(repr.plot.width=9, repr.plot.height=4)\n",
232 "plotPredReal(data, p_np_nn, 3)\n",
233 "plotPredReal(data, p_np_nn, 4)"
238 "execution_count": null,
244 "par(mfrow=c(1,2))\n",
245 "plotFilaments(data, p_np_nn$getIndexInData(3))\n",
246 "plotFilaments(data, p_np_nn$getIndexInData(4))"
251 "execution_count": null,
257 "par(mfrow=c(1,3))\n",
258 "plotSimils(p_np_nn, 3)\n",
259 "plotSimils(p_np_nn, 4)\n",
260 "plotSimils(p_np_nn, 5)"
271 "codemirror_mode": "r",
272 "file_extension": ".r",
273 "mimetype": "text/x-r-source",
275 "pygments_lexer": "r",