2

我有两个问题：

  1. 希望是一个简单的问题，但我想在我的极坐标图中放置一条垂直线，表示数据的开始。我有一个数据集，可以收集一年中每小时的读数，并且我已经设法（在帮助下）获得了我想要的极坐标图。我的数据采用格式（x = yyyy-mm-dd HH:MM:SS 和 y = 数值）。我已经旋转了原点，因此起点不是垂直的（我需要保持这种旋转）。因此，为了表示一年中数据的开始和结束位置，我想添加一条垂直线（见下图）。

  2. 如果这是不可能的，是否可以在开始日期和结束日期之间放置一个空白空间，以便在数据中留下空白，作为开始和结束时间的可视化表示？

这是我当前的示例数据集和代码：

library(lubridate) #make an example dataset
NoOfHours <- as.numeric(ymd_hms("2019-6-1 17:00:00") - ymd_hms("2018-3-01 8:00:00"))*24 
data <- as.data.frame(ymd_hms("2018-3-01 8:00:00") + hours(0:NoOfHours))
colnames(data) <- 'date' 
set.seed(10)
data$level <- runif(nrow(data), min = 0, max = 250)

library(openair)
yeardata <- selectByDate(data, start = "2018-3-1", end = "2019-2-28", year = 2018:2019)

library(ggplot2)

plot <- ggplot(yeardata, aes(x=date, y=level, color = level)) +
  geom_hline(yintercept = seq(0, 300, by = 50), colour = "black", size = 0.75, alpha = 0.3)+ #make my own gridlines so that when on a white background, the gridlines wont cross the text.
  geom_vline(xintercept = as.POSIXct(data$date[1], origin = "1970-01-01 00:00:00 UTC"), colour = "black", size = 0.75, alpha = 0.3, )+
  scale_color_gradient(limits = c(0,200), low="green", high="red", oob = scales::squish)+ #need oob = scales::squish to get values over 200 to be red.
    geom_jitter(alpha = 0.2, size = 2) +# Use a slightly darker palette than normal
 theme(axis.title=element_text(size=16,face="bold"), axis.text.x = element_text(size = 16), axis.text.y = element_text(size = 12))+
   labs(x = NULL, y = bquote('Levels '~(m^2)), color = "Level")+ #bquote to allow superscripts
  scale_y_continuous(breaks = seq(0, 300, 50),
                     limits = c(-100,310))
plot
plot + coord_polar(start = ((2*60/365)*pi))+ #need to have the number of radians to get my start position. If march 1st is the start date, then 60 days have past since Jan 1.
  theme(legend.title = element_text(color = "black", size = 14, face = "bold"), panel.background = element_rect(fill = "white"), panel.grid  = element_blank())

我可以得到一个正常的散点图，我想要的垂直线（在数据集的开头）看起来像这样： 散点图

但是当我尝试把它放到极坐标中时，我得到了这个错误：

“as.POSIXct.numeric(value) 中的错误：必须提供‘原点’”

我试图以几种格式提供来源，最新的是在上面的代码中看到的。理想情况下，该图应如下所示： 理想情节

任何帮助，将不胜感激。

谢谢

4

2 回答 2

1

我无法发布更正的代码，因为 openair 库不会加载到我的系统中，但是您收到此错误的原因是：

“as.POSIXct.numeric(value) 中的错误：必须提供‘原点’”

是因为您在添加数据之前尝试在 vlines 上。如果没有数据已经​​在绘图上，系统就无法弄清楚如何在绘图上放置线条。在大多数情况下，这无关紧要，我正在努力寻找向我解释这一点的参考资料，但 ggplot2 对将对象添加到绘图的顺序有一定的依赖性，特别是对于使用 vline 和日期x 轴。

于 2020-05-10T10:45:19.833 回答
1

这是一种添加线条和文本标签的方法。添加第一个文本注释（在第一个日期之前带有 x）的意外结果似乎会产生您在 #2 中请求的间隙。但不利的一面是，它略微改变了休息时间的轮换，这可能需要对coord_polar(start =术语进行更多调整。

 ggplot(yeardata, aes(x=date, y=level, color = level)) +
   geom_hline(yintercept = seq(0, 300, by = 50), colour = "black", size = 0.75, alpha = 0.3)+ 
   # this annotation layer is drawn over the gridlines, but under the points
   annotate("segment",x = ymd_hms(data$date[1]), xend = ymd_hms(data$date[1]),
                      y = 0, yend = 310, colour = "black", size = 0.75, alpha = 0.3) +


  scale_color_gradient(limits = c(0,200), low="green", high="red", oob = scales::squish)+       
  geom_jitter(alpha = 0.2, size = 2) +# Use a slightly darker palette than normal

  # this text is prior to existing x data and creates a gap
  annotate("text", label = "my text", x = ymd_hms(data$date[1]) - ddays(5), 
            y = 300, angle = 0) +
  annotate("text", label = "my text2", x = ymd_hms(data$date[1]) + ddays(5), 
            y = 300, angle = 0) +
  theme(axis.title=element_text(size=16,face="bold"), axis.text.x = element_text(size = 16), axis.text.y = element_text(size = 12))+
  labs(x = NULL, y = bquote('Levels '~(m^2)), color = "Level")+ #bquote to allow superscripts
  scale_y_continuous(breaks = seq(0, 300, 50),
                     limits = c(-100,310))+
  coord_polar(start = ((2*60/365)*pi))+      
  theme(legend.title = element_text(color = "black", size = 14, face = "bold"),   
        panel.background = element_rect(fill = "white"), 
        panel.grid  = element_blank())

在此处输入图像描述

于 2019-09-17T19:05:51.550 回答