傾向推定のソースを表示

{{出典の明記|date=2023年3月}}
'''傾向推定'''（けいこうすいてい、{{lang-en-short|trend estimation}}）とは、ある過程（プロセス）を測定したものを[[時系列]]として扱い、そのデータの傾向を推定する[[統計学|統計的]]手法である。完全には解明されていない物理的系に対しては、何らかのモデルを構築して測定結果を説明しようと試みる。特に測定結果が増加傾向や減少傾向にあるかを知ることでランダムな振る舞いではないことを判断しようとする。例えば、ある地点での毎日の気温を測ることで季節による変化の傾向や長期的な気象変化の傾向を読み取る。

特に、[[等質性]]の問題は重要である（その時系列は全測定区間で等しく信頼できるか？）。以下では、単純化のためそのような観点をあえて避ける。

== 傾向への適合: 最小二乗法 ==
データ群が与えられ、そのデータから何らかのモデル（この場合、データに適合する関数を意味する）を構築したい場合、選択可能な関数は様々である。しかしそのデータについて何らかの事前の解釈が存在しない場合、最も単純な直線的関数を適合させるのが基本である。

直線に適合させると決めた場合にも様々な手法が存在する。しかし圧倒的に多く使われるのは[[最小二乗法]]である。データの地点 <math>x_i</math> とそのデータ値 <math>y_i</math> について <math>a</math> と <math>b</math> を選択することで次の式を最小化する。

:<math>\sum \{[(ax_i + b) - y_i]^2\}</math> 

解法については[[最小二乗法]]の項目を参照されたし。

以下では、最小二乗法で求めた「傾向」について述べる。問題は、その傾向の有意性であり、「[[有意]]」とはどういうことか、である。

== 無作為データにおける傾向 ==
実データにおける傾向を考える前に、無作為データにおける傾向を理解する必要がある。

[[Image:Random-trends-histogram.png|thumb|right|赤い部分は上位1%、青は5%、緑は10% を示す。この場合本文で述べられている95%の信頼度のV値は 0.2 である。]]

無作為であることが分かっているデータ列（例えばサイコロを振った結果やコンピュータが生成したランダムな数列）があるとき、その傾向を求めるとゼロ傾向となることはほとんどない。しかし、その傾向が極めて小さいことは予測される。ある決まった程度のノイズを含む決まったサイズ（例えば100個）のデータ列があり、それを多数生成する（例えば10万組）と、その10万組のデータ列から傾向を計算することができ、傾向の分布があることを経験的に知ることになる（右図参照）。その分布は（完全にランダムなら）ゼロを中心とする[[正規分布]]となるだろう（[[中心極限定理]]）。以上の手順からある程度の統計的確かさ <math>S</math> を設定することができる（95%が典型的だが、より正確には99%、より大まかなら 90%）。そして、<math>S</math>% の傾向が含まれる範囲を指定する傾向値 <math>V</math> を求めることができる。細かいことを言えば、分布は正と負の両方に広がっており、両方を対象と考える場合もあるし、一方だけを対象と考える場合もある。

以上のように多数回の試行によって経験的に経験の分布を計算することを示した。単純な場合（正規分布の無作為なノイズ）、傾向の分布は正確に求められる。

ここで、それまでのランダムデータ列とおおよそ同じ分散特性の新たなデータ列を考える。そのデータ列が実際に傾向を持つかどうかは分からないので、傾向 <math>T</math> を計算し、それが <math>V</math> より小さいと判明したとする。そこで、確からしさ <math>S</math> の範囲でこのデータの傾向はランダムノイズと区別できないと言える。

しかし、<math>S</math> を選んだとき、残りの <math>1-S</math> の部分がある傾向を持っていると（誤って）宣言する可能性があることに注意されたい。逆に本当に傾向を持つデータ列の残り部分は、傾向を持たないと宣言される可能性がある。

== 傾向＋ノイズとしてのデータ ==
時系列データを解析するため、データ列は傾向要素とノイズ要素から成ると仮定する。

:<math>x_i = at_i + b + e_i\,</math>

<math>a</math> と <math>b</math> は（通常、未知の）定数であり、<math>e</math> は無作為な誤差である。<math>e</math> が何らかの特殊な性質を持つと判明するまでは、[[正規分布]]であると仮定する。<math>e</math> が常に同じ分布であると仮定するのが最も単純だが、そうでない場合（いくつかのデータの分散が非常に大きいなど）、最小二乗法においてそれらのデータの分散の逆で重み付けすることで考慮することができる。

1つの時系列を分析するとき、傾向推定によって <math>e</math> の分散を推定することができる。つまり、傾向推定で求めた <math>at+b</math> に従って残差として <math>e</math> を取り出し、そこから分散を求める。多くの場合、これが <math>e</math> の分散を求める唯一の方法である。

特殊な例として気温の時系列がある。気温データは時間に対して均質でないことが分かっている。一般に気象観測データは最近になるに従って増えており、従って気温の推定に関わる誤差は時と共に減少している。このため気象データの傾向推定を行うにはこれを考慮する。

データ列のノイズが明らかになると、傾向 <math>a</math> が 0 とほとんど差異がないという[[帰無仮説]]によって傾向を検定することができる。上述の無作為データ列の傾向の分散の話から、無作為な（本来傾向のない）データからも傾向が得られることがあることが分かる。もし計算された傾向 <math>a</math> が <math>V</math> より大きければ、その傾向は <math>S</math> の水準においてゼロと有意な差があると言える。

== ノイズの多い時系列 ==
ノイズの多い時系列から傾向を抽出することは難しい。例えば、本来の時系列が 0, 1, 2, 3 という値で、それとは独立した正規分布ノイズ <math>e</math> の[[標準偏差]]を <math>E</math> とする。長さ50の時系列データがあるとき、<math>E=0.1</math> なら傾向は明らかだろう。<math>E=100</math> では傾向はおそらく分かるだろう。しかし、<math>E=10000</math> では傾向はノイズに埋もれてしまうだろう。

具体例として、[[気候変動に関する政府間パネル|IPCC]]が示した過去140年間の気温の記録[http://www.grida.no/climate/ipcc_tar/wg1/figspm-1.htm]を見てみよう。年間気温の分散は約 0.2°C で、傾向は約 0.6°C、95% 信頼度は 0.2°C である（年間の分散と同じ値となっているのは偶然である）。従ってこの傾向は統計的に 0 とは有意な差がある。もっとも、気温の変動の具体的原因はこのデータからは分からない。

== 傾向推定と残差 ==
[[Image:Random-data-plus-trend-r2.png|thumb|right|フィルタリングによって残差の二乗は変化するが、傾向は変わらない]]

最小二乗法による傾向推定では、残差の二乗を推定後に捨てる。それはつまり推定された傾向のラインで説明されるデータの分散の部分がどれだけかということでもある。それは傾向の有意性には関係しない（右図参照）。ノイズの多い系列では残差の二乗は非常に小さいこともあるが、推定の有意性が非常に大きいこともある。フィルタリングを行うと残差の二乗は増大する傾向があるが、推定される傾向そのものやその有意性にはあまり違いが生じない。

== 自己相関的な実データ ==
これまで、データ列は傾向とノイズから構成されるとしてきた。また、ノイズは各データで「独立」であった（[[マルコフ性]]、[[正規分布]]ノイズ）。ノイズが定常的な[[ガウス・マルコフ過程]]に従うという前提は情報最小の原理から生じた。これは統計の容易さという点で大きな意味がある。気象データのような実データはこの前提を満たさないかもしれない。

[[自己相関]]的時系列は[[自己回帰移動平均モデル]]を使ってモデル化される。

==関連項目==
*[[外挿]]
*[[予測]]
*[[予測区間]]
*[[回帰分析]]
*[[天気予報]]

==参考文献==
*Chatfield, C. (1993) "Calculating Interval Forecasts," ''Journal of Business and Economic Statistics,'' '''11(2)''' 121-135.

{{統計学}}
{{DEFAULTSORT:けいこうすいてい}}
[[Category:変化検出]]
[[Category:信号処理]]
[[Category:情報理論]]
[[Category:推計統計学]]
[[Category:計量経済学]]
[[Category:ノイズ]]
[[Category:時系列分析]]
[[Category:数学に関する記事]]