将时间序列绘制为热图

时间:2013-10-03 20:48:10

标签: r ggplot2 time-series heatmap

我有一个带有时间戳和值的数据框(一个月的值大约每5分钟一次)。 现在我想要绘制x轴上的日期,y轴上的时间和颜色的值。如果我可以控制y轴上的分箱(5分钟,10分钟,15分钟......值,作为每个分档的平均值),那将是最好的。

我正在尝试ggplot2,但我无法得到合理的结果。

p_heat <- ggplot(data = data,
                 aes(x = days(timestamp),
                     y = minutes(timestamp) + hours(timestamp) * 60, fill = value)) +  
            geom_tile()

可以按如下方式生成测试数据:

data <- data.frame(timestamp = seq(from = as.POSIXct("2013-09-01 00:00:00"), 
                                   to = as.POSIXct("2013-10-01 00:00:00"), by = "5 mins"),
                   value = runif(8641, 0, 1))

这是迄今为止我能得到的最好的:) 我也尝试过使用scale_x / y_date()

我会很高兴看到一些促使我走向正确方向的提示。谢谢!

2 个答案:

答案 0 :(得分:7)

你正在寻找这样的东西吗?

library(ggplot2)   

# create date variable for the x-axis
df$date <- as.Date(df$timestamp, format = "%Y-%m-%d")

# get H:M components
df$hm <- format(df$timestamp, "%H:%M")    

# create y-axis breaks and labels
lab <- with(df, paste(format(df$timestamp, "%H"), "00", sep = ":"))

gg <- ggplot(data = df, aes(x = date, y = hm, fill = value)) +
  geom_tile() +
  scale_y_discrete(breaks = lab)

gg

enter image description here

答案 1 :(得分:0)

也许试试这个:

p_heat <- ggplot(data=data, aes(x=timestamp$mday ,y=timestamp$min+timestamp$hour*60)) 
    + stat_density2d(geom="tile", aes(fill = value), contour = FALSE)

以下是相关文档: http://docs.ggplot2.org/current/stat_density2d.html