Изменение оси Y на плотность
По умолчанию ось Y показывает количество наблюдений, попавших в каждый интервал. Это удобно и легко интерпретируется. Но что, если мы хотим рассматривать график как настоящую кривую плотности — ту, которую он пытается оценить? То есть сумма произведений (ширина столбца) × (высота столбца) должна быть равна 1.
Для этого достаточно добавить y = stat(density) в маппинг эстетик — это пересчитает ось Y от количества наблюдений к эмпирической оценке вероятности. Обратите внимание: форма графика при этом не изменится, изменится лишь интерпретация оси Y.
Давайте применим это к переменной hour_of_day — часу суток, в который был остановлен нарушитель скоростного режима. Кроме того, немного снизьте непрозрачность столбцов, чтобы сквозь них просматривались линии сетки — это облегчит сравнение.
Это упражнение является частью курса
Лучшие практики визуализации в R
Инструкции к упражнению
- Задайте для эстетики x значение
hour_of_day. - Задайте для эстетики y значение
stat(density). - Измените значение параметра
alphaв функцииgeom_histogram()на 0.8.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
ggplot(md_speeding) +
geom_histogram(
# set x and y aesthetics to hour_of_day and stat(density) respectively.
___
# make points see-through by setting alpha to 0.8
___
)