回帰係数の逆数の分布


9

我々は、線形モデルがあるとしすべての標準回帰(ガウス-マルコフ)前提条件を満たしています。我々は、に興味があるθ = 1 / β 1。yi=β0+β1xi+ϵiθ=1/β1

質問1:どのような仮定を配信するために必要であるθを明確に定義されるべき?β 1 ≠ 0重要であろう---他のもの?θ^β1≠0

質問2:エラーが正規分布に従うという仮定を追加します。場合我々は、それを知っているβ 1は、 MLEであり、Gは、(⋅ )単調関数であり、次に、G (β 1 )のためのMLEであるG (β 1)。単調性は近傍にのみ必要であるβ 1?言い換えれば、あるθ = 1 / β MLE?連続マッピング定理は、少なくともこのパラメーターが一貫していることを示しています。β^1g(⋅)g(β^1)g(β1)β1θ^=1/β^

質問3:デルタ方法と、ブートストラップの両方が配信見つけるための両方の適切な手段であるθを?θ^

質問4:どのようにパラメータのためにこれらの答えの変更を行う?γ=β0/β1

余談:問題を整理してx iを与えることを検討するかもしれません 直接パラメータを推定します。ガウス・マルコフの仮定がここではもはや意味をなさないので、これは私にはうまくいかないようです。たとえば、E[ϵ∣y]について話すことはできません。この解釈は正しいですか?

バツ私=β0β1+1β1y私+1β1ε私=γ+θy私+1β1ε私
E[ε|y]

「標準」の仮定には、正規性が含まれていますか?ε私
— whuber

いい視点ね; MLEの部分にその仮定を追加しました。ただし、他の人にとっては必要ありません。
— チャーリー

1
標本分布のどこから、正常でありθは、通常の逆数です。これは、二峰性の平均発散(無限大)と、の意味に関係なくβ 1であり、そして0ザ・デルタ方式で平坦無限になることがあり、従って恐ろしいとなり、通常の漸近MLE近似が悪く、さらにブートストラップう疑わしいかもしれません。β1θβ1
— whuber

@whuber、それについて詳しく教えてください。私の直感は、法線の逆数がどのように二峰性であるかを理解していません。私の推測では、すべての質量は、通常の平均の逆数であろうことになります(ここでは、)。質量が0に近いため、無限平均の可能性が心配でした。ブートストラップと漸近的な結果は、推定されるモーメントの存在を必要とするため、これが最終的にこの問題にかかっています。1/β^1
— チャーリー

1
相互の通常のPDFは、。0では、すべての導関数は0になります。その対数の臨界点を見つけることは、正と負のモードを識別します(σとμ/σに関して簡単に計算されます); |の積分 x| |の積分のように発散する x| /x2=1/| x| 。無限一次モーメントの問題は、すべての法線を含む、0で正の確率密度を持つ確率変数の逆数に関係します。exp⁡(−(1/バツ−μ)2/(2σ2))/(2πバツ2σ)dバツσμ/σ|バツ||バツ|/バツ2=1/|バツ|
— whuber

回答:


3

Q1。場合β 1のMLEであるβ 1、次にθのMLEでθ及びβ 1 ≠ 0、この推定は、よく定義するための十分条件です。β^1β1θ^θβ1≠0

Q2。θ = 1 / βのMLEであるθ MLEの不変性による。さらに、その逆を取得する必要がない場合は、gの単調性は必要ありません。各点でgを明確に定義する必要があるだけです。これは、Nitis Mukhopadhyayによる「Probability and Statistical Inference」の定理7.2.1 pp。350で確認できます。θ^=1/β^θgg

Q3。はい、両方の方法を使用できます。また、プロファイル尤度もチェックします。θ

(θ、γ)γγ^=β^0/β^1

最後に言及したアプローチは正しくありません。実際には、文献で確認できる「キャリブレーションモデル」を検討しています。必要なのは、関心のあるパラメーターに関してパラメーターを再設定することだけです。

これがお役に立てば幸いです。

敬具。


3
ご回答ありがとうございます。あなたが引用している本はありませんが、これらの特性は推定されるモーメントの存在を必要とします。通常の逆数に必要な瞬間があるかどうかはわかりません。私の質問では、この点をもっと明確にすべきだった。
— チャーリー
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.