2次元カーネル密度推定の混乱R

Question

カーネル密度推定器については、下記のように、確率密度関数を推定するために使用されます。

• http://www.mvstat.net/tduong/research/seminars/seminar-2001-05/
• http://scikit-learn.org/stable/modules/density.html

私の混乱は正確にやってkde2dを何についてですか？下の例では、2つの確率変数 - f（a、b）の共分散確率密度関数を推定していますか？そして色は何を意味しますか？

ここでは、私が参照しているコード例を示します。事前に

b<-log10(rgamma(1000,6,3)) 
a<-log10((rweibull(1000,8,2))) 
density<-kde2d(a,b,n=100) 
filled.contour(density,color.palette=colorRampPalette(c('white','blue','yellow','red','darkred'))) 

おかげで、 林

Answer 1

カーネル密度推定量は何ですか？ 本質的に、それはデータのあらゆる点（その点である標準密度の中心）に少しの標準密度曲線をフィットさせ、次にすべての小さな通常密度をカーネル密度推定器に追加します。

説明のため、1次元カーネル密度推定器の画像をone of your linksから追加します。

これまでは、1Dカーネル密度推定値の概要を素早く説明しました。

2次元カーネル密度はどうですか？

# library(MASS) 
b <- log10(rgamma(1000,6,3)) 
a <- log10((rweibull(1000,8,2))) 
# a and b contain 1000 values each. 

density <- kde2d(a,b,n=100) 

機能はmin(a)からmax(b)にmax(a)とmin(b)にグリッドを作成します。 aまたはbのすべての値に小さな1Dの標準密度をフィッティングする代わりに、kde2dはグリッド内のすべての点で小さな2D標準密度に適合するようになりました。 1次元のカーネル密度と同様に、すべての密度値を加算します。

色は何を意味しますか？コメントの中で指摘されているように、推定確率は2つの変数に依存するので、現在は3つの軸があります（a,bおよびestimated probability）。 3軸を視覚化する1つの方法は、等確率の等高線を使用することです。これは空想的ですが、基本的に天気予報から知っている高/低圧画像と同じです。

あなたは推定確率の最も高い値のため、だから、ローからハイに、プロットはwhiteに色付けされます

filled.contour(density,color.palette = colorRampPalette(c('white', 'blue', 
                  'yellow','red', 
                  'darkred'))))) 

を使用してblue、yellow、red、最終的にdarkredされています。