1

我有可以通过以下方式模仿的数据：

set.seed(1234)

foo <- data.frame(month = rep(month.name, each = 24),
              hour = rep(seq(1:24), 12),
              value1 = rnorm(nrow(foo), 60, 1),
              value2 = rnorm(nrow(foo), 60, 1))

foo <- melt(foo, id = c('month', 'hour'))

我想使用 ggplot 创建一个全年的图，显示每个变量每月的 24 小时周期。

这是我到目前为止所尝试的：

t.plot <- ggplot(foo,
             aes(interaction(month,hour), value, group = interaction(variable,hour)))

t.plot <- t.plot + geom_line(aes(colour = variable))
print(t.plot)

我明白了，这会使数据错位。对于这么小的 SD，您会看到前 24 个值应该更接近 60，但它们到处都是。我不明白是什么导致了这种差异。

https://www.dropbox.com/s/rv6uxhe7wk7q35w/foo.png

当我绘制时：

plot(interaction(foo$month,foo$hour)[1:24], foo$value[1:24])

我得到了我期望的形状，但是 xaxis 非常奇怪，而不是我所期望的。

有什么帮助吗？

4

1 回答 1

4

解决方案是将您的日期设置为日期（不是因素的相互作用）

例如

library(lubridate)
library(reshape2)
Date <- as.Date(dmy('01-01-2000') + seq_len(24*365)*hours(1))
foo <- data.frame(Date = Date, 
  value1 = arima.sim(list(order = c(1,1,0), ar = 0.7), n = 24*365-1), 
   value2 = arima.sim(list(order = c(1,1,0), ar = 0.7), n = 24*365-1))
foo_melt <- melt(foo, id = 'Date')

# then you can use `scale_x_date` and `r` and ggplot2 will know they are dates
# load scales library to access date_format and date_breaks
library(scales) 
ggplot(foo_melt, aes(x=Date, y=value, colour = variable)) + 
 geom_line() +
 scale_x_date(breaks = date_breaks('month'), 
              labels = date_format('%b'), expand =c(0,0))

在此处输入图像描述

每月平均编辑 1 天

您可以使用 facet_wrap 按月分面

# using your created foo data set
levels(foo$month) <- sort(month.abb)
foo$month <- factor(foo$month, levels = month.abb)
ggplot(foo, aes(x = hour, y=value, colour = variable)) + 
 facet_wrap(~month) + geom_line() + 
 scale_x_continuous(expand = c(0,0)))

在此处输入图像描述

于 2012-10-04T04:51:31.833 回答