カーネル密度推定（R, MASS パッケージを使用）（書きかけ）

前準備

R システムのインストール

• Windows での Rシステム のインストールは，別ページ »で説明
• Ubuntu での Blender のインストールは，別ページ »で説明

R のプログラミングについては，別のページにまとめている．

MASS パッケージのインストール

1. R の起動
   
2. install.packages("MASS") の実行

   if (!require("MASS")) install.packages("MASS")
   

   初回実行時には，インストールが始まる．

   1. Windows では，初めて何らかのパッケージをインストールするとき，管理者権限がない場合，ユーザ領域にインストールするという表示が出る．「はい」をクリック．
      

      次の表示が出た場合には，「はい」をクリック．

      
   2. ミラーサイトは，日本から選ぶ
      
3. library(MASS) の実行

   library(MASS)
   

4. サンプルの表示

   example(MASS)
   

5. ヘルプの表示

   help(MASS)
   

カーネル密度推定の例

ここでは，２次元空間内の２つの正規分布の重ね合わせを作り，カーネル密度推定を行う．

1. データの準備

   mu1 = c(1, 1)
   mu2 = c(5, 3)
   a1 = data.frame(x = rnorm(N, mu1[1]), y = rnorm(N, mu1[2]))
   a2 = data.frame(x = rnorm(N, mu2[1]), y = rnorm(N, mu2[2]))
   a = rbind(a1, a2)
   plot(a)
   

2. カーネル密度推定

   library(MASS)
   den = kde3d(x = a$x, y = a$y, h=c(bandwidth.nrd(a$x), bandwidth.nrd(a$y)), n=80)
   plot(den)
   image(den)