X-Git-Url: https://git.auder.net/assets/img/logo_Westcastle.png?a=blobdiff_plain;f=reports%2Freport.ipynb;h=f0a06d007b84dd633d9399e2c5fcb11ac654e663;hb=ee8b1b4e3c13f8dcf13a2c8da6a3bef1520c8252;hp=74d6880f31f696aabe484b08c93a80d49d30d76f;hpb=55639673dd1510a02671c4646813ae346cdca4d6;p=talweg.git
diff --git a/reports/report.ipynb b/reports/report.ipynb
index 74d6880..f0a06d0 100644
--- a/reports/report.ipynb
+++ b/reports/report.ipynb
@@ -2,7 +2,10 @@
"cells": [
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",
@@ -29,82 +32,155 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"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\",\n",
- "\tpredict_at=7) #predict from P+1 to P+H included\n",
+ "# Predict from P+1 to P+H included\n",
+ "H = 17\n",
+ "data = getData(ts_data, exo_data, input_tz = \"GMT\", working_tz=\"GMT\", 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\")\n"
+ "indices_np = seq(as.Date(\"2015-04-26\"),as.Date(\"2015-05-02\"),\"days\")"
]
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
- "\n",
- "\n",
"
Pollution par chauffage
"
]
},
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
- "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)"
+ "reload(\"../pkg\")\n",
+ "#p1 = computeForecast(data, indices_ch, \"Neighbors\", \"Zero\", horizon=H, simtype=\"exo\")\n",
+ "#p2 = computeForecast(data, indices_ch, \"Neighbors\", \"Zero\", horizon=H, simtype=\"endo\")\n",
+ "p3 = computeForecast(data, indices_ch, \"Neighbors\", \"Zero\", horizon=H, simtype=\"mix\")\n",
+ "p4 = computeForecast(data, indices_ch, \"Neighbors\", \"Neighbors\", horizon=H, simtype=\"mix\")\n",
+ "#p4 = computeForecast(data, indices_ch, \"Neighbors2\", \"Zero\", horizon=H, simtype=\"exo\")\n",
+ "#p5 = computeForecast(data, indices_ch, \"Neighbors2\", \"Zero\", horizon=H, simtype=\"endo\")\n",
+ "#p6 = computeForecast(data, indices_ch, \"Neighbors2\", \"Zero\", horizon=H, simtype=\"mix\")"
]
},
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
- "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",
- "# Noir: neighbors_mix, bleu: neighbors_exo, vert: moyenne, rouge: persistence\n",
- "\n",
- "i_np = which.min(e_nn_exo$abs$indices)\n",
- "i_p = which.max(e_nn_exo$abs$indices)"
+ "getSimilarDaysIndices(1000,10,TRUE,data)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
+ "outputs": [],
+ "source": [
+ "as.POSIXlt(data$getTime(1000)[1])"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
+ "outputs": [],
+ "source": [
+ "#e1 = computeError(data, p1, H)\n",
+ "#e2 = computeError(data, p2, H)\n",
+ "e3 = computeError(data, p3, H)\n",
+ "e4 = computeError(data, p4, H)\n",
+ "#e5 = computeError(data, p5, H)\n",
+ "#e6 = computeError(data, p6, H)\n",
+ "plotError(list(e3,e4), cols=c(1,2))"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
+ "outputs": [],
+ "source": [
+ "\tfirst_day = 1\n",
+ "params=p3$getParams(3)\n",
+ "\tfilter = (params$indices >= first_day)\n",
+ "\tindices = params$indices[filter]\n",
+ "\tweights = params$weights[filter]\n",
+ "\n",
+ "\n",
+ "\tgaps = sapply(indices, function(i) {\n",
+ "\t\tdata$getSerie(i+1)[1] - tail(data$getSerie(i), 1)\n",
+ "\t})\n",
+ "\tscal_product = weights * gaps\n",
+ "\tnorm_fact = sum( weights[!is.na(scal_product)] )\n",
+ "\tsum(scal_product, na.rm=TRUE) / norm_fact\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [],
+ "source": [
+ "hist(weights)"
]
},
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"options(repr.plot.width=9, repr.plot.height=4)\n",
"par(mfrow=c(1,2))\n",
"\n",
- "plotPredReal(data, p_nn_exo, i_np); title(paste(\"PredReal nn exo day\",i_np))\n",
- "plotPredReal(data, p_nn_exo, i_p); title(paste(\"PredReal nn exo day\",i_p))\n",
- "\n",
- "plotPredReal(data, p_nn_mix, i_np); title(paste(\"PredReal nn mix day\",i_np))\n",
- "plotPredReal(data, p_nn_mix, i_p); title(paste(\"PredReal nn mix day\",i_p))\n",
+ "plotPredReal(data, p3, 3); title(paste(\"PredReal nn exo day\",3))\n",
+ "plotPredReal(data, p3, 5); title(paste(\"PredReal nn exo day\",5))\n",
"\n",
- "plotPredReal(data, p_az, i_np); title(paste(\"PredReal az day\",i_np))\n",
- "plotPredReal(data, p_az, i_p); title(paste(\"PredReal az day\",i_p))\n",
+ "plotPredReal(data, p4, 3); title(paste(\"PredReal nn mix day\",3))\n",
+ "plotPredReal(data, p4, 5); title(paste(\"PredReal nn mix day\",5))\n",
"\n",
"# Bleu: prévue, noir: réalisée"
]
@@ -112,7 +188,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -126,7 +206,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -140,7 +224,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -156,7 +244,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -172,7 +264,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"# Fenêtres sélectionnées dans ]0,10] / endo à gauche, exo à droite\n",
@@ -185,7 +281,10 @@
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",
@@ -195,7 +294,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"p_nn_exo = computeForecast(data, indices_ep, \"Neighbors\", \"Neighbors\",\n",
@@ -211,7 +314,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"e_nn_exo = computeError(data, p_nn_exo, 3)\n",
@@ -230,7 +337,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"options(repr.plot.width=9, repr.plot.height=4)\n",
@@ -251,7 +362,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -265,7 +380,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -279,7 +398,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -295,7 +418,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -311,7 +438,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"# Fenêtres sélectionnées dans ]0,10] / endo à gauche, exo à droite\n",
@@ -324,7 +455,10 @@
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",
@@ -334,7 +468,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"p_nn_exo = computeForecast(data, indices_np, \"Neighbors\", \"Neighbors\",\n",
@@ -350,7 +488,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"e_nn_exo = computeError(data, p_nn_exo, 3)\n",
@@ -369,7 +511,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"options(repr.plot.width=9, repr.plot.height=4)\n",
@@ -390,7 +536,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -404,7 +554,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -418,7 +572,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -434,7 +592,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
@@ -450,7 +612,11 @@
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"# Fenêtres sélectionnées dans ]0,10] / endo à gauche, exo à droite\n",
@@ -463,7 +629,10 @@
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",