「カーネル密度推定」は、何の畳み込みですか?


25

カーネル密度推定の理解を深めようとしています。

ウィキペディアの定義を使用:https : //en.wikipedia.org/wiki/Kernel_density_estimation#Definition

fh^(x)=1ni=1nKh(xxi)=1nhi=1nK(xxih)

レッツテイク矩形与える関数であるた場合間にあるととさもなければ、及び 1であると(ウィンドウサイズ)。K()1x0.50.50h

密度は2つの関数の畳み込みであることは理解していますが、これら2つの関数を定義する方法がわかりません。それらの1つは(おそらく)データの関数であり、Rのすべてのポイントに対して、その場所にあるデータポイントの数(ほとんど)を示します。そして、他の関数は、おそらくウィンドウサイズと組み合わされたカーネル関数の何らかの修正であるはずです。しかし、それをどのように定義するのか分かりません。0

助言がありますか?

Bellowは、(私が疑う)上記で定義した設定を(2つのガウスと混合で)複製するRコードの例であり、その上で、畳み込まれる関数が疑わしいという「証明」を見たい。n=100

# example code:
set.seed(2346639)
x <- c(rnorm(50), rnorm(50,2))
plot(density(x, kernel='rectangular', width=1, n = 10**4))
rug(x)

ここに画像の説明を入力してください


3
底の敷物は大まかな直観を与えます。からまでの各値が、関連する重み持つスパイクであると想像してください。ここで、カーネルの形状と幅を使用して各スパイクを塗り付け、スパイクが同じ形状と幅をとるように変換され、高さが下の領域がます。結果を追加すると、カーネル密度の推定値が得られます。 i = 1 n 1 / n 1 / nxii=1n1/n1/n
ニックコックス

こんにちはニック、コメントありがとうございます。これは、これまで私は、すでに得た直感で、それは私が見て興味があった畳み込みの形に正式にそれを回す:)(!私は今Whuberの答えを通過することを熱望してる)である
タルGalili

回答:


27

データのバッチに対応するはその「経験密度関数」ですX=(x1,x2,,xn)

fX(x)=1ni=1nδ(xxi).

ここで、は「一般化された関数」です。その名前にもかかわらず、それは関数ではありません。それは積分内でのみ使用できる新しい数学的なオブジェクトです。その定義プロパティは、近傍で連続しているコンパクトサポートの関数、δg0

Rδ(x)g(x)dx=g(0).

(名前には、「アトミック」または「ポイント」メジャーと「ディラックデルタ関数」が含まれます。次の計算では、この概念は、片側からのみ連続する関数を含むように拡張されます。)δグラムg

この特徴付けを正当化することは、fX

xfX(y)dy=x1ni=1nδ(yxi)dy=1ni=1nxδ(yxi)dy=1ni=1nRI(yx)δ(yxi)dy=1ni=1nI(xix)=FX(x)

ここで、は通常の経験的CDFで、は通常の特性関数です(引数がtrueの場合は、それ以外の場合は)。(私は上で定義された関数にコンパクトサポートの機能を移動するために必要な基本制限引数スキップ ;ので唯一の範囲内の値のために定義される必要があるコンパクトであり、これは問題ありません。)FXI10RIX

と他の関数の畳み込みは、定義により次のように与えられます。fX(x)k

(fXk)(x)=RfX(xy)k(y)dy=R1ni=1nδ(xyxi)k(y)dy=1ni=1nRδ(xyxi)k(y)dy=1ni=1nk(xix).

せるk(x)=Kh(x)(同じであるKh(x)ウィキペディア式は畳み込みである: -対称カーネルおよび最もカーネルが対称である)我々は主張結果を取得します。


1
2次元の状況は(より口語的に)説明され、GISサイトのgis.stackexchange.com/questions/14374/…で説明されています。
whuber

1
親愛なるフーバー、私はちょうどあなたの答えを喜んで読んだ!説明と詳細に感謝します。あなたの答え(これと他の一般的なもの)は本当に刺激的です。あなた、タル
タルガリリ

1
δg,xig(xi).

1
@whuberありがとう。文一般化された関数δは関数ではありません。それは積分内でのみ使用できる新しい数学的なオブジェクトです。より明確にしました。いつも通りに。;)
Jan Vainer

1
@Janご協力ありがとうございます。このアイデアをこの回答に取り入れました。
whuber
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.