R/seasadj.R

In pli2016/forecast: Forecasting Functions for Time Series and Linear Models

## Generic seasadj functions

seasadj <- function(object,...) UseMethod("seasadj")

seasadj.stl <- function(object, ...)
{
   return(trendcycle(object) + remainder(object))
}

seasadj.decomposed.ts <- function(object, ...)
{
  if(object$type=="additive")
    return(object$x - object$seasonal)
  else
    return(object$x / object$seasonal)
}

seasadj.tbats <- function(object, ...)
{
  return(object$y - seasonal(object))
  # comp <- tbats.components(object)
  # scols <- grep("season",colnames(comp))
  # sa <- comp[,"observed"] - rowSums(comp[,scols,drop=FALSE])
  # # Back transform if necessary
  # if (!is.null(object$lambda))
  #   sa <- InvBoxCox(sa, object$lambda)
  # return(sa)
}

seasadj.seas <- function(object, ...)
{
  return(seasextract_w_na_action(object, "final"))
}