"cells": [
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",
{
"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",
{
"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, \"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\")\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": false
+ },
+ "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(e1,e2,e3,e4,e5,e6), cols=c(1,2,colors()[258], 4,5,6))"
]
},
{
"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",
{
"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",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"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",
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",
{
"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",
{
"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",
{
"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",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"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",
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",
{
"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",
{
"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",
{
"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",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"cell_type": "code",
"execution_count": null,
- "metadata": {},
+ "metadata": {
+ "collapsed": false,
+ "deletable": true,
+ "editable": true
+ },
"outputs": [],
"source": [
"par(mfrow=c(1,2))\n",
{
"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",
},
{
"cell_type": "markdown",
- "metadata": {},
+ "metadata": {
+ "deletable": true,
+ "editable": true
+ },
"source": [
"\n",
"\n",