R でカーネル密度プロットを作成する方法 (例付き)

カーネル密度プロットは、連続曲線を使用してデータセット内の値の分布を表示するプロットの一種です。

カーネル密度プロットはヒストグラムに似ていますが、ヒストグラムで使用されるビンの数の影響を受けないため、分布の形状を表示するのにさらに適しています。

次の方法を使用して、R でカーネル密度プロットを作成できます。

方法 1: 1 カーネル密度プロットを作成する

 #define kernel density
kd <- density(data)

#create kernel density plot
plot(kd)

方法 2: 充填カーネル密度プロットを作成する

 #define kernel density
kd <- density(data)

#create kernel density plot
plot(kd)

#fill in kernel density plot with specific color
polygon(kd, col=' blue ', border=' black ')

方法 3: 複数のカーネル密度プロットを作成する

 #plot first kernel density plot
kd1 <- density(data1)
plot(kd1, col=' blue ')

#plot second kernel density plot
kd2 <- density(data2)
lines(kd2, col=' red ')

#plot third kernel density plot
kd3 <- density(data3)
lines(kd3, col=' purple ')

...

次の例は、各メソッドを実際に使用する方法を示しています。

方法 1: 1 カーネル密度プロットを作成する

次のコードは、R のデータセットのカーネル密度プロットを作成する方法を示しています。

 #create data
data <- c(3, 3, 4, 4, 5, 6, 7, 7, 7, 8, 12, 13, 14, 17, 19, 19)

#define kernel density
kd <- density(data)

#create kernel density plot
plot(kd, main=' Kernel Density Plot of Data ') 

X 軸はデータセット内の値を示し、Y 軸は各値の相対頻度を示します。グラフ上の最高点は、値が最も頻繁に表示される場所を示します。

方法 2: 充填カーネル密度プロットを作成する

次のコードは、特定の境界線の色と塗りつぶしの色を使用してカーネル密度プロットを作成する方法を示しています。

 #create data
data <- c(3, 3, 4, 4, 5, 6, 7, 7, 7, 8, 12, 13, 14, 17, 19, 19)

#define kernel density
kd <- density(data)

#create kernel density plot
plot(kd)

#add color
polygon(kd, col=' steelblue ', border=' black ')

方法 3: 複数のカーネル密度プロットを作成する

次のコードは、R の単一のプロットで複数のカーネル密度プロットを作成する方法を示しています。

 #create datasets
data1 <- c(3, 3, 4, 4, 5, 6, 7, 7, 7, 8, 12, 13, 14, 17, 19, 19)
data2 <- c(12, 3, 14, 14, 4, 5, 6, 10, 14, 7, 7, 8, 10, 12, 17, 20)

#plot first kernel density plot
kd1 <- density(data1)
plot(kd1, col=' blue ', lwd= 2 )

#plot second kernel density plot
kd2 <- density(data2)
lines(kd2, col=' red ', lwd= 2 ) 

同様の構文を使用して、単一のグラフに必要な数のカーネル密度プロットを作成できることに注意してください。

追加リソース

次のチュートリアルでは、R で他の一般的なプロットを作成する方法を説明します。

R で複数の箱ひげ図をプロットする方法
R で複数のヒストグラムをプロットする方法
Rで複数の線を描く方法