X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=reports%2Freport.ipynb;h=74d6880f31f696aabe484b08c93a80d49d30d76f;hb=d09b09b0c811305405caf871d467f1a5ab34770f;hp=05f51de7d0b28c376070dc2e999ab043b7486b59;hpb=ff5df8e310b73883565761ab4b1aa5a0672e9f27;p=talweg.git

diff --git a/reports/report.ipynb b/reports/report.ipynb
index 05f51de..74d6880 100644
--- a/reports/report.ipynb
+++ b/reports/report.ipynb
@@ -29,20 +29,19 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "library(talweg)\n",
     "\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\",\n",
+    "\tpredict_at=7) #predict from P+1 to P+H included\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"
    ]
   },
   {
@@ -57,29 +56,29 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
-    "p_nn_exo = computeForecast(data, indices_ch, \"Neighbors\", \"Neighbors\", simtype=\"exo\", horizon=H)\n",
-    "p_nn_mix = computeForecast(data, indices_ch, \"Neighbors\", \"Neighbors\", simtype=\"mix\", horizon=H)\n",
-    "p_az = computeForecast(data, indices_ch, \"Average\", \"Zero\", horizon=H) #, memory=183)\n",
-    "p_pz = computeForecast(data, indices_ch, \"Persistence\", \"Zero\", horizon=H, same_day=TRUE)"
+    "p_nn_exo = computeForecast(data, indices_ch, \"Neighbors\", \"Neighbors\",\n",
+    "\thorizon=3, simtype=\"exo\")\n",
+    "p_nn_mix = computeForecast(data, indices_ch, \"Neighbors\", \"Neighbors\",\n",
+    "\thorizon=3, simtype=\"mix\")\n",
+    "p_az = computeForecast(data, indices_ch, \"Average\", \"Zero\",\n",
+    "\thorizon=3)\n",
+    "p_pz = computeForecast(data, indices_ch, \"Persistence\", \"Zero\",\n",
+    "\thorizon=3, same_day=TRUE)"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "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, 3)\n",
+    "e_nn_mix = computeError(data, p_nn_mix, 3)\n",
+    "e_az = computeError(data, p_az, 3)\n",
+    "e_pz = computeError(data, p_pz, 3)\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",
@@ -92,9 +91,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "options(repr.plot.width=9, repr.plot.height=4)\n",
@@ -115,9 +112,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -131,9 +126,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -147,9 +140,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -165,9 +156,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -175,7 +164,7 @@
     "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"
    ]
@@ -183,9 +172,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "# Fenêtres sélectionnées dans ]0,10] / endo à gauche, exo à droite\n",
@@ -208,29 +195,29 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
-    "p_nn_exo = computeForecast(data, indices_ep, \"Neighbors\", \"Neighbors\", simtype=\"exo\", horizon=H)\n",
-    "p_nn_mix = computeForecast(data, indices_ep, \"Neighbors\", \"Neighbors\", simtype=\"mix\", horizon=H)\n",
-    "p_az = computeForecast(data, indices_ep, \"Average\", \"Zero\", horizon=H) #, memory=183)\n",
-    "p_pz = computeForecast(data, indices_ep, \"Persistence\", \"Zero\", horizon=H, same_day=TRUE)"
+    "p_nn_exo = computeForecast(data, indices_ep, \"Neighbors\", \"Neighbors\",\n",
+    "\thorizon=3, simtype=\"exo\")\n",
+    "p_nn_mix = computeForecast(data, indices_ep, \"Neighbors\", \"Neighbors\",\n",
+    "\thorizon=3, simtype=\"mix\")\n",
+    "p_az = computeForecast(data, indices_ep, \"Average\", \"Zero\",\n",
+    "\thorizon=3)\n",
+    "p_pz = computeForecast(data, indices_ep, \"Persistence\", \"Zero\",\n",
+    "\thorizon=3, same_day=TRUE)"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "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, 3)\n",
+    "e_nn_mix = computeError(data, p_nn_mix, 3)\n",
+    "e_az = computeError(data, p_az, 3)\n",
+    "e_pz = computeError(data, p_pz, 3)\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",
@@ -243,9 +230,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "options(repr.plot.width=9, repr.plot.height=4)\n",
@@ -266,9 +251,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -282,9 +265,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -298,9 +279,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -316,9 +295,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -326,7 +303,7 @@
     "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"
    ]
@@ -334,9 +311,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "# Fenêtres sélectionnées dans ]0,10] / endo à gauche, exo à droite\n",
@@ -359,29 +334,29 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
-    "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_nn_exo = computeForecast(data, indices_np, \"Neighbors\", \"Neighbors\",\n",
+    "\thorizon=3, simtype=\"exo\")\n",
+    "p_nn_mix = computeForecast(data, indices_np, \"Neighbors\", \"Neighbors\",\n",
+    "\thorizon=3, simtype=\"mix\")\n",
+    "p_az = computeForecast(data, indices_np, \"Average\", \"Zero\",\n",
+    "\thorizon=3)\n",
+    "p_pz = computeForecast(data, indices_np, \"Persistence\", \"Zero\",\n",
+    "\thorizon=3, same_day=FALSE)"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "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, 3)\n",
+    "e_nn_mix = computeError(data, p_nn_mix, 3)\n",
+    "e_az = computeError(data, p_az, 3)\n",
+    "e_pz = computeError(data, p_pz, 3)\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",
@@ -394,9 +369,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "options(repr.plot.width=9, repr.plot.height=4)\n",
@@ -417,9 +390,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -433,9 +404,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -449,9 +418,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -467,9 +434,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "par(mfrow=c(1,2))\n",
@@ -477,7 +442,7 @@
     "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"
    ]
@@ -485,9 +450,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": false
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "# Fenêtres sélectionnées dans ]0,10] / endo à gauche, exo à droite\n",