マローズのCp

提供: testwiki
ナビゲーションに移動 検索に移動

MallowsのCp[1] [2]は、最小二乗法によって推定された回帰モデルの適合度を評価するために用いられる指標である。名前はコリン・リングウッド・マローズにちなむ。モデル選択を行う際に用いられ、ある複数の変数から出力を予測することができるとき、その中から一部の変数を選んで最も良いモデルを見つけることが目的である。Cpの値が小さいほど、モデルが比較的正確であることを意味する。

マローズのCpは、ガウス線形回帰という特殊な場合において赤池情報量基準に相当することが示されている。[3]

定義と性質

マローズのCpは、過剰適合の問題に対する方法である。一般にモデルの変数が増えれば増えるほど、残差平方和などのモデル適合度の指標は常に小さくなる。したがって、残差平方和が最小となるモデルを選択する場合、常にすべての変数を含むモデルが選択されてしまう。代わりに、データのサンプルで計算されたC p統計は、 母集団ターゲットとして平均二乗予測誤差 (MSPE)を推定する。

Ej(Y^jE(YjXj))2σ2

ただし、Y^jj 番目のケースのフィット値、E (Yj | Xj) は j 番目ケースの期待値であり、σ2は誤差分散(全ケース共通の定数とみなされる)である。変数が追加されても、MSPEは自動的に小さくなることはない。この基準での最適なモデルは、サンプルサイズ、さまざまな予測変数の効果量、および変数間の共線性の程度によって決まる。

P個の変数がK>PであるようなK個の変数から選択された場合、Cpは次のように定義される。

Cp=SSEpS2N+2P,

ただし、

  • SSEp=i=1N(YiYpi)2は、P個の変数を持つモデルの残差平方和
  • Y piは、 P リグレッサからのYの i番目の観測の予測
  • S 2は、 K個すべての変数を用いて回帰分析を行った場合の残差平均平方(residual mean square)であり、平均二乗誤差(MSEによって推定される。
  • Nは標本サイズ

その他の定義

次のような線形モデルがあるとする。

Y=β0+β1X1++βpXp+ε

ただし、

  • β0,,βpは予測変数X1,,Xpの係数
  • εは誤差を表す

Cp以下のようにも定義される[4]

Cp=1n(RSS+2dσ^2)

ただし、

  • RSSは、教師データセットの残差平方和
  • テンプレート:Mvarは予測変数の数
  • σ^2は線形モデルの各応答に関連する分散の推定値を指す(すべての予測子を含むモデルで推定される)

この定義によるCpの値は、前掲の定義によるCpの値と等しくないが、いずれの定義においてもCpを最小にするようなモデルは同一である。

制約

Cp基準には主に2つの制約がある[5]

  1. Cp近似は大きなサンプルサイズに対してのみ有効である。
  2. Cp変数選択(または特徴選択)の問題のようなモデルの複雑な集合を扱うことができない[5]

実用

関連項目

参考文献

テンプレート:Reflist

参照