https://github.com/cran/ftsa
Tip revision: e756b408445e4c9988ad29d10bea68f8390481ed authored by Han Lin Shang on 06 January 2012, 06:14:11 UTC
version 2.8
version 2.8
Tip revision: e756b40
mdase.R
mdase = function(forecast, true)
{
if (length(forecast) != length(true))
stop("MdASE: the lengths of input vectors must be the same.")
n = length(true)
truerror = vector(, n - 1)
for(i in 1:(n - 1)){
truerror[i] = abs(true[i + 1] - true[i])
}
qt = (true - forecast) / (sum(truerror) / (n - 1))
scalederror = median(abs(qt))
return(round(scalederror, 6))
}