Commit | Line | Data |
---|---|---|
25b75559 | 1 | #' Average Forecaster |
3d69ff21 | 2 | #' |
102bcfda | 3 | #' Pointwise average of all the series of the same day of week in the past. |
c4c329f6 | 4 | #' |
102bcfda BA |
5 | #' For example, if the current day (argument "today") is a tuesday, then all series |
6 | #' corresponding to wednesdays in the past (until the beginning or memory limit) are | |
7 | #' averaged to provide a smooth prediction. This forecast will most of the time be wrong, | |
8 | #' but will also look plausible enough. | |
9 | #' | |
4e821712 | 10 | #' @usage # AverageForecaster$new(pjump) |
689aa1d3 | 11 | #' |
102bcfda | 12 | #' @docType class |
c4c329f6 | 13 | #' @format R6 class, inherits Forecaster |
3ddf1c12 | 14 | #' @aliases F_Average |
546b0cb6 | 15 | #' |
25b75559 | 16 | AverageForecaster = R6::R6Class("AverageForecaster", |
a66a84b5 | 17 | inherit = Forecaster, |
3d69ff21 | 18 | |
25b75559 | 19 | public = list( |
98e958ca | 20 | predictShape = function(data, today, memory, horizon, ...) |
3d69ff21 BA |
21 | { |
22 | avg = rep(0., horizon) | |
09cf9c19 BA |
23 | first_day = max(1, today-memory) |
24 | index = today-7 + 1 | |
25 | nb_no_na_series = 0 | |
26 | repeat | |
3d69ff21 | 27 | { |
09cf9c19 | 28 | { |
98e958ca | 29 | serie_on_horizon = data$getCenteredSerie(index)[1:horizon] |
09cf9c19 BA |
30 | index = index - 7 |
31 | }; | |
32 | if (!any(is.na(serie_on_horizon))) | |
33 | { | |
34 | avg = avg + serie_on_horizon | |
35 | nb_no_na_series = nb_no_na_series + 1 | |
36 | }; | |
37 | if (index < first_day) | |
38 | break | |
3d69ff21 | 39 | } |
09cf9c19 | 40 | avg / nb_no_na_series |
3d69ff21 BA |
41 | } |
42 | ) | |
43 | ) |