R/mdse.R

In ftsa: Functional Time Series Analysis

mdse <- function(forecast, true)
{
    if (length(forecast) != length(true))
        stop("MSE: the lengths of input vectors must be the same.")
    err = median((true - forecast)^2)
    return(round(err, 4))
}