X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=pkg%2FR%2FcomputeError.R;h=0b1771f2e1769b16603c9837291f5b916facf881;hb=d2ab47a744d8fb29c03a76a7ca2368dae53f9a57;hp=7da103288805fa64d84feecdad8689b60e55de66;hpb=2057c793ad9929ed5bef8663ea28b896c84df0fc;p=talweg.git

diff --git a/pkg/R/computeError.R b/pkg/R/computeError.R
index 7da1032..0b1771f 100644
--- a/pkg/R/computeError.R
+++ b/pkg/R/computeError.R
@@ -1,18 +1,21 @@
 #' Compute error
 #'
-#' Obtain the errors between forecast and data
+#' Compute the errors between forecasted and measured series.
 #'
-#' @param data Dataset, object of class \code{Data} output of \code{getData}
-#' @param pred Forecast object, class \code{Forecast} output of \code{computeForecast}
+#' @param data Object of class \code{Data} output of \code{getData}
+#' @param pred Object of class \code{Forecast} output of \code{computeForecast}
 #' @param horizon Horizon where to compute the error
 #'   (<= horizon used in \code{computeForecast})
 #'
-#' @return A list (abs,MAPE) of lists (day,indices)
+#' @return A list (abs,MAPE) of lists (day,indices). The "indices" slots contain series
+#'   of size L where L is the number of predicted days; i-th value is the averaged error
+#'   (absolute or MAPE) on day i. The "day" slots contain curves of errors, for each time
+#'   step, averaged on the L forecasting days.
 #'
 #' @export
-computeError = function(data, pred, horizon=data$getStdHorizon())
+computeError = function(data, pred, predict_from, horizon=length(data$getSerie(1)))
 {
-	L = forecast$getSize()
+	L = pred$getSize()
 	mape_day = rep(0, horizon)
 	abs_day = rep(0, horizon)
 	mape_indices = rep(NA, L)
@@ -21,9 +24,9 @@ computeError = function(data, pred, horizon=data$getStdHorizon())
 	nb_no_NA_data = 0
 	for (i in seq_len(L))
 	{
-		index = forecast$getIndexInData(i)
-		serie = data$getSerie(index+1)[1:horizon]
-		forecast = pred$getForecast(i)[1:horizon]
+		index = pred$getIndexInData(i)
+		serie = data$getSerie(index)[predict_from:horizon]
+		forecast = pred$getForecast(i)[predict_from:horizon]
 		if (!any(is.na(serie)) && !any(is.na(forecast)))
 		{
 			nb_no_NA_data = nb_no_NA_data + 1