tests/FittingConstantTest.R

In dbellettini/RobustHoltWinters:

library(RobustHoltWinters)

test.deseasonalized <- function () {
    series <- ts(frequency = 7, rep(10, 50))
    model <- RobustHoltWinters(series)

    checkEqualsNumeric(10.0, model$coefficients['a'], "Deseasonalized series should be 10.0, equal to the constant value", 0.001)
    checkEqualsNumeric(0.0, model$coefficients['b'], "Trend data should be nearly 0", 0.001)

    for (i in 1:7) {
        checkEqualsNumeric(0.0, model$coefficients[paste('s', i, sep='')], "Seasonal data should be nearly 0", 0.001)
    }
}