標準偏差の閉じた形の不偏推定量はどの分布にありますか?


16

正規分布の場合、標準偏差の不偏推定量があります:

σ^不偏=Γn12)Γ(n2)12k=1n(xix¯)2

この結果があまり知られていない理由は、それが大部分の重要な輸入の問題ではなくむしろ骨cur品であるように思われます。証明はこのスレッドでカバーされています。正規分布の重要な特性を利用します。

1σ2k=1n(xix¯)2χn12

そこから、少しの作業で、、およびこの回答を倍数として識別することにより、の結果を推測でき。σ σ公平E(k=1n(xix¯)2)σσ^unbiased

これにより、他のどの分布が標準偏差の閉形式の不偏推定量を持っているのか興味があります。分散の不偏推定量とは異なり、これは明らかに分布固有です。さらに、他の分布の推定量を見つけるために証明を適応させるのは簡単ではありません。

スキュー正規分布には、2次形式の優れた分布特性がいくつかあります。これは、使用した正規分布特性が事実上特別なケースです(正規分布は特殊なタイプのスキュー正規分布なので)。このメソッドをそれらに拡張します。しかし、他のディストリビューションでは、まったく異なるアプローチが必要と思われます。

そのような推定量が知られている他の分布はありますか?


1
技術的な注意散漫を無視すると、答えの性質がより明確になります。通常の場合、あなたが書いたもののほとんどは、結論に本当に関連していません。重要なのは、この特定の推定量のバイアスの量がのみの関数であることです(データから推定する必要がある他の分布パラメーターに依存しません)。n
whuber

@whuberあなたがほのめかしている一般的な考えを見ることができると思います。そして明らかに「関数だけ」が必要です。しかし、私はそれが十分だとは思わない-私たちがいくつかの素晴らしい分布結果にアクセスできなければ、「閉じた形」の側面がどのように扱いやすいかを見ることができない。n
シルバーフィッシュ

5
「閉じた形式」の意味に依存します。たとえば、ある人にとってシータ関数は「閉じている」かもしれませんが、別の人にとってはそれは無限の積、べき級数、または複素積分です。考えてみると、それがまさにガンマ関数です:-)。
whuber

@whuber良い点!「この特定の推定量のバイアスの量」とは、のバイアス(質問にリストされている推定量がゼロバイアスであるのではなく)はnの関数であることを意味します(幸いなことにσでも)公平な推定量を見つけるために簡単に再配置できるような方法で)?snσ
シルバーフィッシュ

1
@whuber:ロケーションスケールファミリには同様の式が必要です。ただし、の関数は扱いにくい積分である可能性があることを指摘しました。n
西安

回答:


10

これは、直接質問に接続されていないが、あるピータービッケルとエーリッヒレーマン1968論文は、その分布の凸ファミリー、状態、機能の不偏推定量が存在するQ Fのサンプルサイズのために(N場合)十分な大きさとする場合にのみQ α F + 1 - α Gがの多項式であり、0 α 1Fq(F)nq(αF+(1α)G)0α1。この定理は、ガウス分布のコレクションが凸ではないため(ガウス分布の混合はガウス分布ではないため)、ここの問題には当てはまりません。

問題の結果の拡張は、任意の電源ことである標準偏差をunbiasedly推定することができ、十分な観察が提供されるα < 0。これは、結果 1σαα<0 ということσは、スケール(ユニークな)であるためのパラメータΣ N K = 1XI- ˉ X2

1σ2k=1n(xix¯)2χn12
σk=1n(xix¯)2

この通常の設定は、次に、任意の位置規模ファミリーに拡張することができる 有限の分散とσ 2。確かに、

X1,,Xniidτ1f(τ1{xμ})
σ2
  1. 分散のみの関数であり、τ
    varμ,τ(X)=Eμ,τ[(Xμ)2]=τ2E0,1[X2]
    τ
  2. 二乗和の形式の期待有するτ2ψnは
    Eμ,τ[k=1n(XiX¯)2]=τ2Eμ,τ[k=1nτ2(XiμX¯+μ)2]=τ2E0,1[k=1n(XiX¯)2]
    τ2ψ(n)
  3. 同様に、任意の電力のためのであり、期待値は有限です。
    Eμ,τ[{k=1n(XiX¯)2}α]=τ2αE0,1[{k=1n(XiX¯)2}α]

10

おそらくよく知られているケースですが、それでもケースです。
連続的な一様分布考えます。iidサンプルが与えられた場合、最大次数統計X n )には期待値がありますU(0,θ)X(n)

E(X(n))=nn+1θ

分布の標準偏差は

σ=θ23

だから、推定 σ = 1

σ^=123n+1nX(n)

に対して明らかに偏りがない。 σ

これは、範囲の不偏推定量を使用できるため、分布の下限も不明な場合に一般化され、標準偏差は再び範囲の線形関数になります(基本的に上記のとおりです)。

これは@whuberのコメントを例示しており、「バイアスの量はのみの関数」(および、おそらく既知の定数)であるため、決定論的に修正できます。そして、これはこの場合です。n


4
ここで難しいのは、世界で均一分布の標準偏差に関心があるときですか?(+1)
シャドウトーカー

1
@ssdecontrolそれは素晴らしい質問です!
ください...-アレコスパパドプロス

2
θ^

1
@Silverfishどういう意味で貧しい?いくつかのクイックシミュレーションは、これが通常の標準偏差よりも低いMSEを持っていることを示しています(驚いた)。
デイブ

@Dave面白い!私は最大次数の統計だけを見ていたので貧弱だという結論に飛びつきましたが、私も驚いています!シミュレーションを行うことの価値を示しています...
Silverfish
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.