Как создать графики плотности ядра в r (с примерами)

График плотности ядра — это тип графика, который отображает распределение значений в наборе данных с помощью непрерывной кривой.

График плотности ядра похож на гистограмму , но он даже лучше отображает форму распределения, поскольку на него не влияет количество ячеек, используемых в гистограмме.

Мы можем использовать следующие методы для создания графика плотности ядра в R:

Способ 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. Создайте одноядерный график плотности.

Следующий код показывает, как создать график плотности ядра для набора данных в 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