Кривая плотности для сгруппированных наблюдений

Question

У меня есть набор данных, где я имею почасовые наблюдения за различными видами.Для каждого часа все наблюдения объединяются для каждого вида, и это было сделано в течение 24 часов.Ниже приводится часть данных:

Hr  Sp1 Sp2 Sp3 Sp4 Sp5
1   8   25  1   7   7
2   3   14  0   50  4
3   5   25  1   4   10
4   7   26  1   4   11
5   5   19  0   7   11

Я хочу построить кривую плотности для наблюдений каждого вида в зависимости от времени.Следовательно, у меня есть 24 непрерывные группы (0100 часов, 0200 часов ... 2400 часов), начиная с 0000 часов как ось X, ось Y - это подсчет или наблюдения для каждого часа в течение 24 часов для каждого вида.Пожалуйста, предложите, как я могу сделать это в R?Я также хочу, чтобы все кривые плотности видов на одном участке.Как мне это сделать?

Большое спасибо!

Answer 1

library(reshape2)
library(ggplot2)

dt <- data.frame(Hr = 1:5,
                 Sp1 = c(8, 3, 5, 7, 5),
                 Sp2 = c(25, 14, 25, 26, 19),
                 Sp3 = c(1, 0, 1, 1, 0),
                 Sp4 = c(7, 50, 4, 4, 7),
                 Sp5 = c(7, 4, 10, 11, 11))

data <- melt(dt, id = c("Hr"))
data <- data.frame(Hr = rep(data$Hr, data$value),
                   variable = rep(data$variable, data$value))

ggplot(data, aes(Hr, colour = variable)) + geom_density()

http://ggplot2.tidyverse.org/reference/geom_density.html