绘制2个ACF ggplot2 - 显示最小值颜色

时间:2016-02-26 14:59:28

标签: r ggplot2 advanced-custom-fields

我正在尝试使用ggplot在R中绘制2个ACF。我使用geom_segment()并且它工作正常,但是一个令人烦恼的细节是,根据我声明geom_segment的顺序,如果我在最大值之前声明这些,则绘图将覆盖最小的ACF。如果一个ACF总是具有更大的值(在之前我会对其进行geom_segment的情况),这不会是问题。我知道这听起来令人困惑,所以这是一个有效的例子:

library(ggplot2)
library(grid)
library(plyr)

set.seed(123)
data1 <- arima.sim(n = 10000, model = list(ar = 0.4,ma=0.4))
data2 <- arima.sim(n = 10000, model = list(ma = c(0.9,0.3)))

acf1 <- acf(data1, plot = FALSE,lag.max=15)
acf1DF <- with(acf1, data.frame(lag, acf))
acf1DF$lag <- as.integer(acf1DF$lag)
acf2 <- acf(data2, plot = FALSE,lag.max=15)
acf2DF <- with(acf2, data.frame(lag, acf))
acf2DF$lag <- as.integer(acf2DF$lag)

q <- ggplot(data = acf1DF, mapping = aes(x = lag, y = acf)) +
  geom_hline(aes(yintercept = 0)) +
  geom_segment(mapping = aes(xend = lag, yend = 0,color="red")) +
  geom_segment(data=acf2DF,mapping = aes(x=lag,y=acf,xend = lag, yend = 0),
    inherit.aes=FALSE,show.legend=FALSE)
q

请注意,我看不到第一个滞后的最小值。

感谢。

1 个答案:

答案 0 :(得分:3)

是否可以为您创建一个条形图?分段重叠的问题使您的情节不易理解,可能指向使用不同的可视化作为解决方案。这是一个条形图:

首先,我们通常在ggplot中合并数据,当你多次调用geom时,你做得不对。

acf1DF$group <- 1
acf2DF$group <- 2

dat <- rbind(acf1DF,acf2DF)

然后我们可以制作一个条形图:

p1 <- ggplot(dat,aes(x=factor(lag), y=acf,fill=factor(group))) + 
  geom_bar(width=0.3,stat="identity", position=position_dodge())

enter image description here 或者,如果我们真的想要段而不是条形(例如因为我们想要一个连续的x轴,或者像线条的外观一样作为图例),我们可以在绘制之前在第2组中添加一个稍微偏移的延迟

dat$lag_offset <- dat$lag + as.numeric(dat$group==2) * 0.1

然后绘制细分:

p2 <- ggplot(dat, aes(x=lag_offset, xend=lag_offset, y=0,yend=acf, color=factor(group))) +
  geom_segment() +
  scale_color_manual(values=c(`1`="black",`2`="red"))

enter image description here