公開日: 11/25/2021

Image shown here収束スコア検定

反復計算の収束に失敗すると、「収束スコア検定」というレポートが表示されます。この検定では、「最尤推定値は、現在の反復計算で得られた最終ステップでのパラメータ値と等しい」という帰無仮説に対するスコア検定が表示されます。この検定は、収束基準として用いている相対的な傾き(gradient)を、検定統計量として用いています。なお、このスコア検定は、本当の最尤推定値がどのような値を取るかを知る必要がありません。

Image shown hereスコア検定

最初に、推定するパラメータθが1個の場合を見てみましょう。θにおける対数尤度関数をl、得られたデータをxとします。θにおける対数尤度関数の導関数のベクトルを、「スコア」と呼びます。

Equation shown here

以下のように、対数尤度関数を2回微分したことによって得られる行列を、「(観察された)情報量行列」と呼びます。

Equation shown here

H0: θ = θ0 に対するスコア検定の検定統計量は次のとおりです。

Equation shown here

この検定統計量は、帰無仮説のもとで、自由度が1の漸近的なカイ2乗分布に従います。

スコア検定は、複数のパラメータに対して一般化できます。複数のパラメータを含むベクトルをθとすると、H0: θ = θ0 に対するスコア検定の検定統計量は次のとおりです。

Equation shown here

ここで

Equation shown here

および

Equation shown here

また、Uは、行列Uの転置を表します。

この検定統計量は、自由度kのカイ2乗分布に漸近的に従います。ここで、kは、制約されていない自由パラメータの個数です。

Image shown here相対的な勾配

[混合モデル]手法での収束基準は相対的な勾配gH1gに基づいています。ここで、g(θ) = U(θ)は対数尤度関数の勾配、H(θ) = I(θ)はヘッセ行列です。

θ0 を、反復計算が終了したときのθの値とします。θ0 で評価される相対的な勾配が、スコア検定の検定統計量となります。p値は、自由度kのカイ2乗分布に基づいて計算されます。このp値によって、「最尤推定値は、θ0と等しい」という帰無仮説が検定されます。なお、kは、「変量効果の共分散パラメータ推定値」レポートにリストされている、制約されていないパラメータの個数です。

より詳細な情報が必要な場合や、質問があるときは、JMPユーザーコミュニティで答えを見つけましょう (community.jmp.com).