"\n",
"ts_data = read.csv(system.file(\"extdata\",\"pm10_mesures_H_loc_report.csv\",package=\"talweg\"))\n",
"exo_data = read.csv(system.file(\"extdata\",\"meteo_extra_noNAs.csv\",package=\"talweg\"))\n",
- "data = getData(ts_data, exo_data, input_tz = \"Europe/Paris\", working_tz=\"Europe/Paris\", predict_at=13)\n",
+ "data = getData(ts_data, exo_data, input_tz = \"Europe/Paris\", working_tz=\"Europe/Paris\", predict_at=7)\n",
"\n",
"indices_ch = seq(as.Date(\"2015-01-18\"),as.Date(\"2015-01-24\"),\"days\")\n",
"indices_ep = seq(as.Date(\"2015-03-15\"),as.Date(\"2015-03-21\"),\"days\")\n",
- "indices_np = seq(as.Date(\"2015-04-26\"),as.Date(\"2015-05-02\"),\"days\")"
+ "indices_np = seq(as.Date(\"2015-04-26\"),as.Date(\"2015-05-02\"),\"days\")\n",
+ "\n",
+ "H = 3 #predict from 2pm to 4pm"
]
},
{
},
"outputs": [],
"source": [
- "e_nn_exo = computeError(data, p_nn_exo)\n",
- "e_nn_mix = computeError(data, p_nn_mix)\n",
- "e_az = computeError(data, p_az)\n",
- "e_pz = computeError(data, p_pz)\n",
+ "e_nn_exo = computeError(data, p_nn_exo, H)\n",
+ "e_nn_mix = computeError(data, p_nn_mix, H)\n",
+ "e_az = computeError(data, p_az, H)\n",
+ "e_pz = computeError(data, p_pz, H)\n",
+ "\n",
"options(repr.plot.width=9, repr.plot.height=7)\n",
"plotError(list(e_nn_mix, e_pz, e_az, e_nn_exo), cols=c(1,2,colors()[258], 4))\n",
"\n",
"plotSimils(p_nn_exo, i_p); title(paste(\"Weights nn exo day\",i_p))\n",
"\n",
"plotSimils(p_nn_mix, i_np); title(paste(\"Weights nn mix day\",i_np))\n",
- "plotSimils(p_nn_mix, i_p); title(paste(\"Weights nn mix day\",i_p)\n",
+ "plotSimils(p_nn_mix, i_p); title(paste(\"Weights nn mix day\",i_p))\n",
"\n",
"# - pollué à gauche, + pollué à droite"
]
},
"outputs": [],
"source": [
- "e_nn_exo = computeError(data, p_nn_exo)\n",
- "e_nn_mix = computeError(data, p_nn_mix)\n",
- "e_az = computeError(data, p_az)\n",
- "e_pz = computeError(data, p_pz)\n",
+ "e_nn_exo = computeError(data, p_nn_exo, H)\n",
+ "e_nn_mix = computeError(data, p_nn_mix, H)\n",
+ "e_az = computeError(data, p_az, H)\n",
+ "e_pz = computeError(data, p_pz, H)\n",
"options(repr.plot.width=9, repr.plot.height=7)\n",
"plotError(list(e_nn_mix, e_pz, e_az, e_nn_exo), cols=c(1,2,colors()[258], 4))\n",
"\n",
"plotSimils(p_nn_exo, i_p); title(paste(\"Weights nn exo day\",i_p))\n",
"\n",
"plotSimils(p_nn_mix, i_np); title(paste(\"Weights nn mix day\",i_np))\n",
- "plotSimils(p_nn_mix, i_p); title(paste(\"Weights nn mix day\",i_p)\n",
+ "plotSimils(p_nn_mix, i_p); title(paste(\"Weights nn mix day\",i_p))\n",
"\n",
"# - pollué à gauche, + pollué à droite"
]
"p_nn_exo = computeForecast(data, indices_np, \"Neighbors\", \"Neighbors\", simtype=\"exo\", horizon=H)\n",
"p_nn_mix = computeForecast(data, indices_np, \"Neighbors\", \"Neighbors\", simtype=\"mix\", horizon=H)\n",
"p_az = computeForecast(data, indices_np, \"Average\", \"Zero\", horizon=H) #, memory=183)\n",
- "p_pz = computeForecast(data, indices_np, \"Persistence\", \"Zero\", horizon=H, same_day=TRUE)"
+ "p_pz = computeForecast(data, indices_np, \"Persistence\", \"Zero\", horizon=H, same_day=FALSE)"
]
},
{
},
"outputs": [],
"source": [
- "e_nn_exo = computeError(data, p_nn_exo)\n",
- "e_nn_mix = computeError(data, p_nn_mix)\n",
- "e_az = computeError(data, p_az)\n",
- "e_pz = computeError(data, p_pz)\n",
+ "e_nn_exo = computeError(data, p_nn_exo, H)\n",
+ "e_nn_mix = computeError(data, p_nn_mix, H)\n",
+ "e_az = computeError(data, p_az, H)\n",
+ "e_pz = computeError(data, p_pz, H)\n",
"options(repr.plot.width=9, repr.plot.height=7)\n",
"plotError(list(e_nn_mix, e_pz, e_az, e_nn_exo), cols=c(1,2,colors()[258], 4))\n",
"\n",
"plotSimils(p_nn_exo, i_p); title(paste(\"Weights nn exo day\",i_p))\n",
"\n",
"plotSimils(p_nn_mix, i_np); title(paste(\"Weights nn mix day\",i_np))\n",
- "plotSimils(p_nn_mix, i_p); title(paste(\"Weights nn mix day\",i_p)\n",
+ "plotSimils(p_nn_mix, i_p); title(paste(\"Weights nn mix day\",i_p))\n",
"\n",
"# - pollué à gauche, + pollué à droite"
]