[talweg.git] / R / S_Persistence.R

#' @include ShapeForecaster.R
#'
#' @title Persistence Shape Forecaster
#'
#' @description Return the last centered last (similar) day curve.
#'   Inherits \code{\link{ShapeForecaster}}
PersistenceShapeForecaster = setRefClass(
	Class = "PersistenceShapeForecaster",
	contains = "ShapeForecaster",

	methods = list(
		initialize = function(...)
		{
			callSuper(...)
		},
		predict = function(today, memory, horizon, ...)
		{
			#return centered last (similar) day curve, avoiding NAs until memory is run
			first_day = max(1, today-memory)
			index = today-7 + 1
			repeat
			{
				{
					last_similar_serie = data$getCenteredSerie(index)[1:horizon]
					index = index - 7
				};
				if (!any(is.na(last_similar_serie)))
					return (last_similar_serie);
				if (index < first_day)
					return (NA)
			}
		}
	)
)
Commit	Line	Data
3d69ff21 BA	1	#' @include ShapeForecaster.R
	2	#'
	3	#' @title Persistence Shape Forecaster
	4	#'
09cf9c19 BA	5	#' @description Return the last centered last (similar) day curve.
09cf9c19 BA	6	#' Inherits \code{\link{ShapeForecaster}}
3d69ff21 BA	7	PersistenceShapeForecaster = setRefClass(
	8	Class = "PersistenceShapeForecaster",
	9	contains = "ShapeForecaster",
	10
	11	methods = list(
	12	initialize = function(...)
	13	{
	14	callSuper(...)
	15	},
	16	predict = function(today, memory, horizon, ...)
	17	{
09cf9c19 BA	18	#return centered last (similar) day curve, avoiding NAs until memory is run
	19	first_day = max(1, today-memory)
	20	index = today-7 + 1
	21	repeat
	22	{
	23	{
	24	last_similar_serie = data$getCenteredSerie(index)[1:horizon]
	25	index = index - 7
	26	};
	27	if (!any(is.na(last_similar_serie)))
	28	return (last_similar_serie);
	29	if (index < first_day)
	30	return (NA)
	31	}
3d69ff21 BA	32	}
	33	)
	34	)