`smdape` <- function(forecast, true)
{
if (length(forecast) != length(true))
stop("SMDAPE: the lengths of input vectors must be the same.")
err = median(200 * (true - forecast) / (true + forecast))
return(round(err, 6))
}