2

我想绘制一个时间序列,在图中不包括中间的一段时间。如果我单独绘制系列,x 轴上只有一个索引,这就是我得到的。内部排除点集不会出现。

x <- rnorm(50)
Dates <- seq(as.Date("2008-1-1"), by = "day", length.out = length(x))
dummy <- c(rep(1, 25), rep(0, 10), rep(1, length(x) - 35))
plot(x[dummy == 1])

然而,一旦日期位于 x 轴上,R 就会尽职尽责地呈现准确的真实时间刻度,包括排除的日期。这会在绘图上产生一个空白区域。

plot(Dates[dummy == 1], x[dummy == 1])

如何在 x 轴上获取日期,但不显示排除日期的空白区域?

4

2 回答 2

1

三种选择:

1. ggplot2显然,ggplot2 不允许使用不连续的轴,但您可以使用它facet_wrap来获得类似的效果。

# get the data
  x = rnorm(50)
  df <- data.frame( x = x,
                  Dates = seq(as.Date("2008-1-1"), by = "day", length.out = length(x)) ,
                  dummy = c(rep(1, 25), rep(0, 10), rep(1, length(x) - 35)))

  df$f <- ifelse(df$Dates <= "2008-01-25",  c("A"), c("B")) 

# plot
  ggplot( subset(df, dummy==1)) +
    geom_point(aes(x= Dates, y=x)) +
    facet_wrap(~f , scales = "free_x")

在此处输入图像描述

2.基础R

plot(df$x ~ df$Dates, col= ifelse( df$f=="A", "blue", "red"), data=subset(df, dummy==1))

在此处输入图像描述

3. plotrix 另一种选择是使用gap.plot{plotrix}. 代码将类似于下面的代码。但是,我无法弄清楚如何在具有date值的轴上进行中断。也许这会和其他问题。

library(plotrix)

gap.plot(Dates[dummy == 1], x[dummy == 1], gap=c(24,35), gap.axis="x")
于 2016-05-27T21:19:39.853 回答
0

我想我想通了。axis按照我上面提出的思路,我不得不长时间摆弄这个命令才能让它把日期标签放在正确的位置。这是我使用的:

plot(x[dummy == 1], xaxt = "n", xlab = "")  # plot with no x-axis title or tick labels
Dates1_index <- seq(1,length(Dates1), by = 5)  # set the tick positions
axis(1, at = Dates1_index, labels = format(Dates1[Dates1_index], "%b %d"), las = 2)

成功后,我现在同意@alistaire 的观点,它看起来非常具有误导性。也许如果我在休息处放一条垂直虚线......

于 2016-05-28T05:32:47.557 回答