Proposal Presentation

Methus Detchusananart and Munchuporn Techachaianun
12 November 2014

MIAD Forecasting Model

Multi-Intervention Aggregate-Disaggregate Forecasting Model

Methus Detchusananart and Munchuporn Techachaianun

Advisor: Asst.Prof.Seeronk Prichanont, Ph.D.
Co-Advisor: Asst. Prof.Krung Sinapiromsaran, Ph.D.

Information and Communication Engineering, 2014
Faculty of Engineering, Chulalongkorn University

Agenda

  • Introduction and Background
  • Literature Review
  • Objective and Scope
  • Methodology
  • Plan and Schedule
  • Expected Outcome

Automatic Teller Machine

Problems

Motivation

To optimize the amount of cash notes in the machine

Background

  • Automatic Teller Machine (ATM)
  • ATM replenishment
  • Dataset from NN5
  • Time series
    • Exponential smoothing
    • ARMA
    • ARIMA

History of ATM

ATM in Thailand

ATM Replenishment

  • Cassettes are replaced by prepared full cassettes and are returned to a financial institution

Dataset from NN5

plot of chunk unnamed-chunk-1

Time series

  • Time series is a collection of data recorded over a period of time—weekly, monthly, quarterly, or yearly.
  • Time series forecasting is the use of a model to predict future values based on previously observed values.
  • Examples: Dow Jones Industrial Average, Daily Temperature, Gold price

Time series

Patterns of Time Series

  • Trend: long term movements in the mean
  • Seasonality: cyclical fluctuations according to calendar
  • Long-run cycle: cyclical fluctuations such as a business cycles
  • Outliers: irregular patterns

Autoregressive model

To predict the value at the present time using the value at the previous time.

EQUATION AR(p)

\[ X_t = \Phi_1 X_{t-1} + \Phi_2 X_{t-2} + ... + \Phi_p X_{t-p} + \varepsilon_t \]

where \( \varepsilon_t \) is independent identically distributed (iid) with Normal Distribution of mean 0 and constant variance \( \sigma_\varepsilon^2 \)

\( \varepsilon_t \) ~ \( N(0,\sigma_\varepsilon^2) \)

Moving Average model

“Past Error”

EQUATION MA(q)

\[ X_t = \varepsilon_t + \theta_1 \varepsilon_{t-1} + ... + \theta_q \varepsilon_{t-q} \]

where \( \varepsilon_t \) ~ \( N(0,\sigma_\varepsilon^2) \)

Autoregressive Moving-average (ARMA)

EQUATION ARMA(p,q)

\[ X_t = \Phi_1 X_{t-1} + ... + \Phi_p X_{t-p} + \varepsilon_t + \theta_1 \varepsilon_{t-1} + ... + \theta_q \varepsilon_{t-q} \]

For example ARMA(1,1)

\[ X_t = \Phi_1 X_{t-1} + \varepsilon_t + \theta_1\varepsilon_{t-1} \]

Autoregressive Integrated Moving-average (ARIMA)

From ARMA(p,q)

\[ X_t = \Phi_1 X_{t-1} + ... + \Phi_p X_{t-p} + \varepsilon_t + \theta_1 \varepsilon_{t-1} + ... + \theta_q \varepsilon_{t-q} \]

\[ X_t - \Phi_1 X_{t-1} - ... - \Phi_p X_{t-p} = \varepsilon_t + \theta_1 \varepsilon_{t-1} + ... + \theta_q \varepsilon_{t-q} \]

  • B is the backward operator

\[ (1-\Phi_1 B - ... - \Phi_p B^p) X_t = (1 + \theta_1 B + ... + \theta_q B^q)\varepsilon_t \]

\[ \Phi(B)X_t = \theta(B)\varepsilon_t \]

\[ ARIMA(p,d,q): \Phi(B)(1-B)^d X_t = \theta(B)\varepsilon_t \]

Exponential Smoothing

EQUATION Holt-Winters(\( \alpha \))

\[ F_{t+1} = \alpha X_t + (1-\alpha)F_t \]

where \( \alpha \) is a smoothing constant and \( 0 \leq \alpha \leq 1 \)

Limitations

ARMA

  • Can not handle seasonality or trend

ARIMA

  • Not accurated long term forecast

Exponential smoothing

  • Hard to find the best alpha (smoothing constant)
  • Can not handle seasonality or trend

Literature review

  • Multi-intervention forecasting model for ATM, Ms. Wanicha Shithamarach, Ms. Apichaya Phannongwah, (2013).

Literature review (Cont.)

plot of chunk unnamed-chunk-3

Literature review (Cont.)

Literature review (Cont.)

Literature review (Cont.)

Subsequence Suitable model MAPE Rank Forecast value at t =
1 Additive model 0.1032391 6 33,34,35,36,…
2 ARIMA(0,1,0) model 0.0960504 5 33,35,37,39,…
3 Exponential Smoothing 0.2024615 10 34,36,38,40,…
4 Exponential Smoothing 0.1315447 7 34,37,40,…
5 Exponential Smoothing 0.1458527 8 35,38,41,…
6 ARIMA(0,1,0) model 0.08148158 2 33,36,39,…
7 ARIMA(0,1,0) model 0.08949418 4 33,37,41,…
8 ARIMA(0,1,0) model 0.08702252 3 34,38,42,…
9 ARIMA(0,1,0) model 0.0771577 1 35,39,43,…
10 ARIMA(0,1,0) model 0.17808997 9 36,40,44,…

Objective

  • To create new forecasting model to optimize the amount of cash note inventory and apply the new model to data from NN5
  • To compare with the existing model such as exponential smoothing, ARMA, ARIMA, and Multi-intervention model by using the mean absolute percentage error (MAPE) method.

Scope

  • Use R language on Rstudio software to formulate the model and forecast data
  • The data used in this project is from NN5
  • Dataset is portioned into 2 parts
    • 80% in-sample
    • 20% out-sample
  • This model will be created from using Exponential smoothing, ARMA, and ARIMA model
  • Model is tested by the mean absolute percentage error (MAPE) \[ M = \frac{1}{n}\sum_{t=1}^{n} \left|\frac{A_t-F_t}{A_t}\right| \]

Methodology

plot of chunk unnamed-chunk-4

Methodology

Methodology - Disaggregate data

  • Use the same distribution of data as the previous period

    EX. For week 9, use the distribution from week 1 and 5

Plan and schedule

Expected outcome

  • The accuracy of outcome from new forecasting model is expected to be better comparing with other existing models.
  • The opportunity to publish this project