《指数平滑法》のソースを表示

'''【しすうへいかつほう (exponential smoothing) 】'''

1. 指数平滑法の基本式

　時間の変化に従って与えられるデータ群, すなわち, 時系列データ (例えば需要量系列など) を用いた予測方式の共通点は, 先行するデータ群をつぎに続くデータ群に関連づけて推定を行うことである. この方法の一つで[[指数平滑法]]の基礎となる方式が移動平均法である. いま, {\it t}時点の移動平均値を$m_t$, 用いられるデータ群の項数を{\it l}, ({\it t}－{\it n})時点のデータを$d_{t-n}$とし, $d_{t-n}$の係数を$a_{n+1}$とすれば, {\it t}時点の移動平均値$m_t$は次式で表される. 


\begin{equation}
m_{t}=a_{1}d_{t}+a_{2}d_{t-1}+{\cdots}+a_{l}d_{t-l+1}
\label{B-H-02+siki1}
\end{equation}


\begin{displaymath}
\mbox{ただし, }\sum_{n=1}^{l}a_{n}=1
\end{displaymath}


　さて, 指数平滑法([1], [2])はこの(1)式の係数\begin{math}a_{1} ,a_{2},\cdots ,a_{l}\end{math}に対して, 現在時点に近い程ウエイトを大きくし, 過去にさかのぼる程ウエイトを小さくしていく指数型の考え方を導入したもので, 指数型加重移動平均法とも呼ばれている. すなわち, この方式による{\it t}時点の推定値$m_t$は次式のように表わされる. (ただし,

\begin{math}
0 \leq 1-a \leq 1
\end{math}
)
\begin{equation}
m_{t}=a\{d_{t}+(1-a)d_{t-1}+{\cdots}+(1-a)^{n}d_{t-n}+{\cdots}\}
\label{B-H-02+siki2}
\end{equation}

　この(2)式の重み係数\begin{math}a, a(1-a), a(1-a)^{2}, \cdots \cdots\end{math}の総計は1となる. なお, この場合の$m_t$は({\it t}+1)時点の予測値として用いられる. ここで, (\ref{B-H-02+siki2})式と同様な考え方で$m_{t-1}$を算出すると次式が得られる. 

\begin{equation}
m_{t-1}= a \{d_{t-1}+ (1-a) d_{t-2}+{\cdots}+ (1-a) ^{n-1} d_{t-n}+{\cdots}\}
  \label{B-H-02+siki3}
\end{equation}

　(3)式の両辺に(1－{\it a})を掛けて(2)式と対比すると, 次式のような指数平滑法の基本式(定数型モデル)が導出される. 

\begin{equation}
m_{t}= a d_{t}+ (1-a) m_{t-1}
  \label{B-H-02+siki4}
\end{equation}

　この(\ref{B-H-02+siki4})式の係数{\it a}を平滑化定数と呼んでいる. 

2. 平滑化定数の値

　上記の平滑化定数{\it a}は, 原則として, 0と1の間の値をとる. この中, {\it a}=1のときは\begin{math}m_{t}=d_{t}\end{math}となり推定値は同時点のデータと等しくなる. 一方, {\it a}=0のときは\begin{math}m_{t}=m_{t-1}\end{math}となり推定値は一時点前の先行する推定値と等しくなる. {\it a}を中間の値(0＜{\it a}＜1)にとった場合にはある程度ランダムな変動の影響を受けることになる. {\it a}の値を最適に決めることは難しい問題であるが, 0.5より若干小さい値をとる場合が比較的多い. 

3. 傾向を考慮した場合

　時系列データに傾向がない定数型モデルの場合には, (4)式の$m_t$が({\it t}+1)時点の有効な予測値となるが, もし, 上昇, あるいは下降の傾向がある場合は, この値は不満足なものになる. この問題を解決するために傾向を考慮した指数平滑法(直線型傾向モデル)が提案されている. このモデルでは2番目の変数として$t$時点の傾向の推定値$r_t$を導入している. この直線型傾向モデルはR. G. Brown[3], [4] により提示されたが, いま, ({\it t}+1)時点の予測値を$y_t$とすれば(\ref{B-H-02+siki5})～(\ref{B-H-02+siki7})式のように表わされる.


\begin{equation}
y_{t}= m_{t}+{ \frac{1-a}{a}} r_{t}
  \label{B-H-02+siki5}
\end{equation}
  \begin{equation}
m_{t}= a d_{t}+ (1-a) m_{t-1}
  \label{B-H-02+siki6}
\end{equation}
\begin{equation}
r_{t}= a (y_{t}- y_{t-1})+ (1-a) r_{t-1}
  \label{B-H-02+siki7}
\end{equation}

　このモデルの{\it k}時点先の予測値$y_{t+k}$は次式のように示される. 
　
   \begin{equation}
y_{t+k}= y_{t}+ k r_{t}
  \label{B-H-02+siki8}
\end{equation}

　なお, (4)式や(5)式などで用いられている$m_t$や$r_t$の初期値は, それまでのデータにより推定される. それらの値はとくに重要な値ではないので比較的単純な近似法を用いればよい. 

4. 季節変動を考慮した場合

 　時系列データに季節変動がある場合には, 季節型モデルが利用される. このモデルの代表的なものはP. R. Winters [5] によって提示されているが, いま, {\itt}時点の季節変動指数値$S_t$, 3つの平滑化定数をそれぞれ{\it a}, {\it b}, {\it c}とすればこのモデルは(9)～(11)式のように表わされる. 

\begin{equation}
y_{t}= a \frac{d_{t}}{S_{t-L}}+(1-a) (y_{t-1}+ r_{t-1})
  \label{B-H-02+siki9}
\end{equation}
\begin{equation}
S_{t}= b \frac{d_{t}}{y_{t}}+(1-b) S_{t-L}
  \label{B-H-02+siki10}
\end{equation}
\begin{equation}
r_{t}= c (y_{t}- y_{t-1})+(1-c) r_{t-1}
  \label{B-H-02+siki11}
\end{equation}

　ただし, $t-L$は{\it t}より1年前の時点. 

このモデルの{\it k}時点先の予測値$y_{t+k}$は次式のように示される. 

 \begin{equation}
y_{t+k}= (y_{t}+k r_{t})S_{t+k-L}
  \label{B-H-02+siki12}
\end{equation}

5. その他の指数平滑法

　上記以外のモデルとしては, $y_{t-2}$までを用いた2次のモデルや定数型モデルの推定値をデータとして同じモデルを繰り返し用いる2重や3重のモデルも提案されている. 



----

'''参考文献'''

[1] I. C. I Monograph, No.2, ''Short-Term Forecasting'', Imperial Chemical Industries Limited, 1964.

[2] C. C. Holt, ''Forecasting Seasonals and Trends by Exponentially Weighted Moving Averages'', Carnegie Institute of Technology, Pittsburgh, Pennsylvania, 1957.

[3] R. G. Brown, ''Statistical Forecasting for Inventory Control'', McGraw-Hill, 1959.

[4] R. G. Brown and R. F. Meyer, "The Fundamental Theorem of Exponential Smoothing," ''Operations Research'', '''19''' (1961), 673-685.

[5] P. R. Winters, "Forecasting Sales by Exponentially Weighted Moving Averages," ''Management Science'', '''6''' (1960), 324-342.