我正在尝试使用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
请注意,我看不到第一个滞后的最小值。
感谢。
答案 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())
或者,如果我们真的想要段而不是条形(例如因为我们想要一个连续的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"))