R ggplot2 geom_boxplot grouper 对于自定义分组器失败

Question

我正在尝试使用 geom_boxplot 中的 group= 选项，它适用于一个分组功能，但不适用于另一个。第一个地块运行，第二个和第三个地块（实际上相同，称为不同）都未能产生 2017 年之前的 2 个月箱线图和 2017 年的一个月箱线图，正如石斑鱼打算的那样。对于 grouper 函数 ggplot 声明 警告消息：position_dodge 需要不重叠的 x 间隔“ 但图形之间的 X 值相同。显然与我的 groupdates 函数相关，但组似乎构造正确。欢迎提出建议。与谢谢。

library(tidyverse)
library(lubridate)
# I want two month groups before 2017, and one-month groups in 2017

groupdates <- function(date) {
  month_candidate <-case_when(
    year(date) < 2017 ~ paste0(year(date), "-", (floor(((0:11)/12)*6)*2)+1),
    TRUE ~ paste0(year(date), "_", month(date))
  )
  month_candidate2 <-case_when(
    (str_length(month_candidate)==6) ~ paste0(str_sub(month_candidate,1,5), "0", str_sub(month_candidate,6)),
    TRUE ~ month_candidate
  )
  return(month_candidate2)
}

generate_fake_date_time <- function(N, st="2015/01/02", et="2017/02/28") {
       st <- as.POSIXct(as.Date(st))
       et <- as.POSIXct(as.Date(et))
       dt <- as.numeric(difftime(et,st,unit="sec"))
       ev <- sort(runif(N, 0, dt))
       rt <- st + ev
}

n=5000
set.seed(250)
test <-as.data.frame(generate_fake_date_time(n))
colnames(test) <- "posixctdate"
test$ranvalue <- month(test$posixctdate)+runif(length(test), 0,1)
test$grouped_time <-groupdates(test$posixctdate)
table(test$grouped_time)

ggplot(test)+geom_boxplot(aes(x=posixctdate, y=ranvalue, group=paste0(year(posixctdate), "_", month(posixctdate))))
#ggplot(test)+geom_violin(aes(x=posixctdate, y=ranvalue, group=junk))
ggplot(test)+geom_boxplot(aes(x=posixctdate, y=ranvalue, group=grouped_time))
ggplot(test)+geom_boxplot(aes(x=posixctdate, y=ranvalue, group=groupdates(posixctdate)))

    sessionInfo()

Answer 1

如果我正确理解了您的问题，您应该考虑修改您的 groupdates 函数。

我只修改了第三行：

• ceiling 而不是 floor
• month(date) 而不是 0:11

导致：

groupdates <- function(date) {
    month_candidate <-case_when(
        year(date) < 2017 ~ paste0(year(date), "-", (ceiling(((month(date))/12)*6)*2)+1),
        TRUE ~ paste0(year(date), "_", month(date))
    )
    month_candidate2 <-case_when(
        (str_length(month_candidate)==6) ~ paste0(str_sub(month_candidate,1,5), "0", str_sub(month_candidate,6)),
        TRUE ~ month_candidate
    )
    return(month_candidate2)
}

我还修改了ranvalue 的计算以获得更好的分布，我敢打赌你想使用nrow 而不是length：

test$ranvalue <- month(test$posixctdate) + runif(nrow(test), 0, 1)
test$grouped_time <-groupdates(test$posixctdate)
table(test$grouped_time)

输出（没有变化）：

ggplot(test)+geom_boxplot(aes(x=posixctdate, y=ranvalue, group=grouped_time))