https://github.com/cran/ftsa
Tip revision: 4682eed2913f88477018153fc0004a471e3b7f86 authored by Han Lin Shang on 15 September 2018, 07:00:03 UTC
version 5.0
version 5.0
Tip revision: 4682eed
mdase.R
mdase <- function(forecast, outsampletrue, insampletrue)
{
if(length(forecast) != length(outsampletrue))
stop("MdASE: the lengths of input vectors must be the same.")
n = length(insampletrue)
insamplerr = vector(, (n - 1))
for(i in 1:(n - 1))
{
insamplerr[i] = abs(insampletrue[i+1] - insampletrue[i])
}
qt = (outsampletrue - forecast)/(sum(insamplerr)/(n - 1))
scalederror = median(abs(qt))
return(round(scalederror, 6))
}