Оценка плотности ядерным методом

Лучшие практики визуализации в R

Nick Strayer

Instructor

Когда гистограммы не справляются

Данные с несколькими выраженными пиками
Малый объём данных

Графики плотности ядерного метода

Помещаем «ядро» на каждую точку данных
Суммируем высоты всех перекрывающихся ядер

Построение KDE в ggplot

Достаточно заменить geom_histogram() на geom_density()

sample_n(md_speeding, 100) %>%
  ggplot(aes(x = percentage_over_limit)) +
  # Swap out geom_histogram() 
    geom_density(
      # Fill in curve with color
      fill = 'steelblue',
      # Standard deviation of kernel
      bw = 8 
    )

Новый параметр ширины

Необходимо настроить стандартное отклонение ядра для каждой точки

Отображение всех данных

Используйте geom_rug(), чтобы отобразить все данные под KDE в виде черточек

p <-sample_n(md_speeding, 100) %>%
  ggplot(aes(x = percentage_over_limit)) +
  geom_density(
    fill = 'steelblue', # fill in curve with color
    bw = 8 # standard deviation of kernel
  ) 

p + geom_rug(alpha = 0.4)

Давайте потренируемся!

Лучшие практики визуализации в R

Preparing Video For Download...