数理統計学における尤度比検定の定義
📂数理統計学数理統計学における尤度比検定の定義
定義
H0:H1:θ∈Θ0θ∈Θ0c
上記の仮説検定について、次の統計量 λ を 尤度比検定統計量likelihood Ratio test statisticと言う。
λ(x):=supΘL(θ∣x)supΘ0L(θ∣x)
与えられたc∈[0,1]に対して棄却域 {x:λ(x)≤c} を持つすべての仮説検定を 尤度比検定likelihood Ratio testと言い、通常、LRTと略す。
説明
λ の定義で、分子はsupΘ0で、分母はsupΘでスプレマムを探している。帰無仮説下の部分パラメータ空間Θ0は全パラメータ空間の部分集合Θ0⊆Θで、自然に0≤λ(x)≤1が成立する。この比が0に近いということは、その分、帰無仮説下でのパラメータがそんなにもっともらしくないという意味になる。
統計学を初めて学ぶ時に、基本的な確率分布論から始めて、t分布、F分布、カイ二乗分布などの検定統計量を別々に学ぶよりずっとわかりやすいと感じるだろう。もちろん、尤度比検定にはそれなりのモチベーションがあるが、前述の検定と違って、どんなビルドアップもなくても話が通じる。
例: 正規分布
実際にLRTを行う場合は、スプレマムsupが反映された状態でなければならない。分母はパラメータ空間Θ全体で最大になるので、最尤推定量を使用し、分子は帰無仮説下で最大になるように設定する。
H0:H1:θ=θ0θ=θ0
分散が既知の正規分布N(θ,σ2)のランダムサンプルX1,⋯,Xnに対して、上記の仮説検定を考えよう。ここで、分母は母平均θに対する最尤推定量である標本平均xˉを使用しなければならず、分子は帰無仮説のパラメータ空間が単集合Θ0={θ0}であるため、θ0をそのまま使用すればよい。数式的に次のように求められる。
λ(x)====supΘL(θ∣x)supΘ0L(θ∣x)L(xˉ∣x)L(θ0∣x)(2π)−n/2exp(−∑(xk−xˉ)2/2)(2π)−n/2exp(−∑(xk−θ0)2/2)exp(−n(xˉ−θ0)2/2)
このとき、λ(x)は正確にxˉ=θ0のとき1であり、差が大きくなるほど0に近づくことに注意しよう。もちろん、LRTの数式的な定義では、0から1までであることは分かっているが、計算された結果を見ると、実際に標本分散が帰無仮説下の母平均とどれほど類似しているかに応じて、直感的に仮説検定を行うことができる。