《指数平滑法》

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

1. 指数平滑法の基本式

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

$m_{t}=a_{1}d_{t}+a_{2}d_{t-1}+{\cdots }+a_{l}d_{t-l+1}\,$　　　　　$(1)\,$

ただし, $\ \sum _{n=1}^{l}a_{n}=1\,$

　さて, 指数平滑法([1], [2])はこの(1)式の係数構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle a_{1} ,a_{2},\cdots ,a_{l}\, } に対して, 現在時点に近い程ウエイトを大きくし, 過去にさかのぼる程ウエイトを小さくしていく指数型の考え方を導入したもので, 指数型加重移動平均法とも呼ばれている. すなわち, この方式による${\it {t}}\,$時点の推定値構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_t\, } は次式のように表わされる. (ただし, $0\leq 1-a\leq 1\,$)

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_{t}=a\{d_{t}+(1-a)d_{t-1}+{\cdots}+(1-a)^{n}d_{t-n}+{\cdots}\} \, } 　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (2)\, }

　この(2)式の重み係数構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle a, a(1-a), a(1-a)^{2}, \cdots \cdots\, } の総計は1となる. なお, この場合の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_t\, } は構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ({\it t}+1)\, } 時点の予測値として用いられる. ここで, (2)式と同様な考え方で構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_{t-1}\, } を算出すると次式が得られる.

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_{t-1}= a \{d_{t-1}+ (1-a) d_{t-2}+{\cdots}+ (1-a) ^{n-1} d_{t-n}+{\cdots}\} \, } 　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (3)\, }

　(3)式の両辺に構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (1-{\it a})\, } を掛けて(2)式と対比すると, 次式のような指数平滑法の基本式(定数型モデル)が導出される.

$m_{t}=ad_{t}+(1-a)m_{t-1}\,$　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (4)\, }

　この(4)式の係数${\it {a}}\,$を平滑化定数と呼んでいる.

2. 平滑化定数の値

　上記の平滑化定数構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {\it a}\, } は, 原則として, 0と1の間の値をとる. この中, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {\it a}=1\, } のときは構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_{t}=d_{t}\, } となり推定値は同時点のデータと等しくなる. 一方,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {\it a}=0\, } のときは$m_{t}=m_{t-1}\,$となり推定値は一時点前の先行する推定値と等しくなる. 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {\it a}\, } を中間の値構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (0<{\it a}<1)\, } にとった場合にはある程度ランダムな変動の影響を受けることになる. 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {\it a}\, } の値を最適に決めることは難しい問題であるが, 0.5より若干小さい値をとる場合が比較的多い.

3. 傾向を考慮した場合

　時系列データに傾向がない定数型モデルの場合には, (4)式の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m_t\, } が構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ({\it t}+1)\, } 時点の有効な予測値となるが, もし, 上昇, あるいは下降の傾向がある場合は, この値は不満足なものになる. この問題を解決するために傾向を考慮した指数平滑法(直線型傾向モデル)が提案されている. このモデルでは2番目の変数として構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t\, } 時点の傾向の推定値r_tを導入している. この直線型傾向モデルはR. G. Brown[3], [4] により提示されたが, いま, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ({\it t}+1)\, } 時点の予測値を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_t\, } とすれば(5)～(7)式のように表わされる.

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_{t}= m_{t}+{ \frac{1-a}{a}} r_{t} \, } 　　　　　$(5)\,$

$m_{t}=ad_{t}+(1-a)m_{t-1}\,$　　　　　$(6)\,$

$r_{t}=a(y_{t}-y_{t-1})+(1-a)r_{t-1}\,$　　　　　$(7)\,$

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

$y_{t+k}=y_{t}+kr_{t}\,$　　　　　$(8)\,$

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

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

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

$y_{t}=a{\frac {d_{t}}{S_{t-L}}}+(1-a)(y_{t-1}+r_{t-1})\,$　　　　　$(9)\,$

$S_{t}=b{\frac {d_{t}}{y_{t}}}+(1-b)S_{t-L}\,$　　　　　$(10)\,$

$r_{t}=c(y_{t}-y_{t-1})+(1-c)r_{t-1}\,$　　　　　$(11)\,$

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

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

$y_{t+k}=(y_{t}+kr_{t})S_{t+k-L}\,$　　　　　$(12)\,$

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.