Mean absolute scaled error

In statistics, the mean absolute scaled error (MASE) is a measure of the accuracy of forecasts. It is the mean absolute error of the forecast values, divided by the mean absolute error of the in-sample one-step naive forecast. It was proposed in 2005 by statistician Rob J. Hyndman and Professor of Decision Sciences Anne B. Koehler, who described it as a "generally applicable measurement of forecast accuracy without the problems seen in the other measurements."^[1] The mean absolute scaled error has favorable properties when compared to other methods for calculating forecast errors, such as root-mean-square-deviation, and is therefore recommended for determining comparative accuracy of forecasts.^[2]

Rationale

The mean absolute scaled error has the following desirable properties:^[3]

1. Scale invariance: The mean absolute scaled error is independent of the scale of the data, so can be used to compare forecasts across data sets with different scales.
2. Predictable behavior as ${\displaystyle y_{t}\rightarrow 0}$ : Percentage forecast accuracy measures such as the Mean absolute percentage error (MAPE) rely on division of ${\displaystyle y_{t}}$, skewing the distribution of the MAPE for values of ${\displaystyle y_{t}}$ near or equal to 0. This is especially problematic for data sets whose scales do not have a meaningful 0, such as temperature in Celsius or Fahrenheit, and for intermittent demand data sets, where ${\displaystyle y_{t}=0}$ occurs frequently.
3. Symmetry: The mean absolute scaled error penalizes positive and negative forecast errors equally, and penalizes errors in large forecasts and small forecasts equally. In contrast, the MAPE and median absolute percentage error (MdAPE) fail both of these criteria, while the "symmetric" sMAPE and sMdAPE^[4] fail the second criterion.
4. Interpretability: The mean absolute scaled error can be easily interpreted, as values greater than one indicate that in-sample one-step forecasts from the naïve method perform better than the forecast values under consideration.
5. Asymptotic normality of the MASE: The Diebold-Mariano test for one-step forecasts is used to test the statistical significance of the difference between two sets of forecasts.^[5][6][7] To perform hypothesis testing with the Diebold-Mariano test statistic, it is desirable for ${\displaystyle DM\sim N(0,1)}$, where ${\displaystyle DM}$ is the value of the test statistic. The DM statistic for the MASE has been empirically shown to approximate this distribution, while the mean relative absolute error (MRAE), MAPE and sMAPE do not.^[2]

Non seasonal time series

For a non-seasonal time series,^[8] the mean absolute scaled error is estimated by

${\displaystyle \mathrm {MASE} =\mathrm {mean} \left({\frac {\left|e_{j}\right|}{{\frac {1}{T-1}}\sum _{t=2}^{T}\left|Y_{t}-Y_{t-1}\right|}}\right)={\frac {{\frac {1}{J}}\sum _{j}\left|e_{j}\right|}{{\frac {1}{T-1}}\sum _{t=2}^{T}\left|Y_{t}-Y_{t-1}\right|}}}$^[3]

where the numerator e_j is the forecast error for a given period (with J, the number of forecasts), defined as the actual value (Y_j) minus the forecast value (F_j) for that period: e_j = Y_j − F_j, and the denominator is the mean absolute error of the one-step "naive forecast method" on the training set (here defined as t = 1..T),^[8] which uses the actual value from the prior period as the forecast: F_t = Y_t−1^[9]

Seasonal time series

For a seasonal time series, the mean absolute scaled error is estimated in a manner similar to the method for non-seasonal time series:

${\displaystyle \mathrm {MASE} =\mathrm {mean} \left({\frac {\left|e_{j}\right|}{{\frac {1}{T-m}}\sum _{t=m+1}^{T}\left|Y_{t}-Y_{t-m}\right|}}\right)={\frac {{\frac {1}{J}}\sum _{j}\left|e_{j}\right|}{{\frac {1}{T-m}}\sum _{t=m+1}^{T}\left|Y_{t}-Y_{t-m}\right|}}}$^[8]

The main difference with the method for non-seasonal time series, is that the denominator is the mean absolute error of the one-step "seasonal naive forecast method" on the training set,^[8] which uses the actual value from the prior season as the forecast: F_t = Y_t−m,^[9] where m is the seasonal period.

This scale-free error metric "can be used to compare forecast methods on a single series and also to compare forecast accuracy between series. This metric is well suited to intermittent-demand series (a data set containing a large amount of zeros) because it never gives infinite or undefined values^[1] except in the irrelevant case where all historical data are equal.^[3]

When comparing forecasting methods, the method with the lowest MASE is the preferred method.

Non-time series data

For non-time series data, the mean of the data (${\displaystyle {\bar {Y}}}$) can be used as the "base" forecast.^[10]

${\displaystyle \mathrm {MASE} =\mathrm {mean} \left({\frac {\left|e_{j}\right|}{{\frac {1}{J}}\sum _{j=1}^{J}\left|Y_{j}-{\bar {Y}}\right|}}\right)={\frac {{\frac {1}{J}}\sum _{j}\left|e_{j}\right|}{{\frac {1}{J}}\sum _{j}\left|Y_{j}-{\bar {Y}}\right|}}}$

In this case the MASE is the Mean absolute error divided by the Mean Absolute Deviation.

See also

• Mean squared error
• Mean absolute error
• Mean absolute percentage error
• Root-mean-square deviation
• Test set
• Fraction of variance unexplained

References