Question

我有过去几个月的测量数据：

变量

 x <- df$DatoTid
 y <- df$Partikler
 color <- df$Opgave

我试图根据时间戳绘制数据，以便在x轴上显示一天中的时间，而不是特定的POSIXct日期时间。我希望x轴的标签和刻度为fx“ 00:00”，“ 01:00”，...“ 24:00”。这样中午在x轴的中间。

到目前为止，我尝试将datetime值转换为字符。看起来还不太好（如您所见，轴刻度和标签都消失了。可能其他情况也有误）。

有人可以帮助我吗？并且请让我知道如何为您上传数据。我不知道如何添加巨大的.csv文件。...

# Rounding up to nearest 10 min: 
head(df)
df$Tid2 <- format(strptime("1970-01-01", "%Y-%m-%d", tz="CET") + 
round(as.numeric(df$DatoTid)/300)*300 + 3600, "%Y-%m-%d %H:%M:%S")
head(df)
df$Tid2 <- as.character(df$Tid2)
str(df)

x <- df$Tid2
y <- df$Partikler
color <- df$Opgave

plot2 <- ggplot(data = df, aes(x = x, y = y, color = color)) +
  geom_point(shape=16, alpha=0.6, size=1.8) +
  scale_y_continuous(labels=function(x) format(x, big.mark = ".", decimal.mark = ",", scientific = FALSE)) +
  scale_x_discrete(breaks=c("00:00:00", "06:00:00", "09:00:00", "12:00:00", "18:00:00", "21:00:00")) +
  scale_color_discrete(name = "Case") +
  xlab(" ") +
  ylab(expression(paste("Partikelkoncentration [pt/cc]"))) +
  myTheme + 
  theme(legend.text=element_text(size=8), legend.title=element_text(size=8))
plot2

Answer 1

我会通过制作一个新的时间戳来解决这一问题，该时间戳使用一天的时间，但要使用现有时间戳的小时/分钟/秒。

首先，这是数据的虚构版本，这里使用Partikler中的线性趋势：

library(tidyverse); library(lubridate)
df <- data_frame(Tid2 = seq.POSIXt(from = ymd_h(2019010100), 
                                   to = ymd_h(2019011500), by = 60*60),
                 Partikler = seq(from = 0, to = 2.5E5, along.with = Tid2),
                 Opgave = as.factor(floor_date(Tid2, "3 days")))

# Here's a plot that's structurally similar to yours:
ggplot(df, aes(Tid2, Partikler, col = Opgave)) + 
  geom_point() +
  scale_color_discrete(name = "Case")

现在，如果我们将时间戳记更改为同一天，则可以像在ggplot中一样正常地控制它们，但是它们会折叠成一天的时间。我们还可以更改x轴，使其不提及时间戳记的日期部分：

df2 <- df %>%
  mutate(Tid2_sameday = ymd_hms(paste(Sys.Date(), 
                                      hour(Tid2), minute(Tid2), second(Tid2))))


ggplot(df2, aes(Tid2_sameday, Partikler, col = Opgave)) + 
  geom_point() +
  scale_color_discrete(name = "Case")  +
  scale_x_datetime(date_labels = "%H:%M")

时间序列ggplot中x轴上的时间戳

1 个答案: