我正在尝试使用acf
制作ggplot2
图。我的代码如下:
library(ggplot2)
x = lh
conf.level = 0.95
ciline = qnorm((1 - conf.level)/2)/sqrt(length(x))
bacf = acf(x, plot = FALSE)
bacfdf = with(bacf, data.frame(lag, acf))
ggplot(data=bacfdf, mapping=aes(x=lag, y=acf)) +
geom_bar(stat="identity", position = "identity") +
ggtitle("Orders")
通过这个,我能够生成ggplot2
acf
自相关图,其滞后范围从零到20。
如何编辑此项以使ggplot2
的图示滞后范围从-10到10?
我基本上使用了此来源的代码来提出上面发布的代码:http://ask.programmershare.com/387_17805747/
答案 0 :(得分:1)
Autocorrelation function为calculated为
mx <- mean(x)
sum((x[1:(N-1)] - mx)*(x[2:N] - mx)) / sum((x-mx)^2)
对于滞后= -1,唯一会改变的是x[1:(N-1)]
和x[2:N]
改变他们的位置。它是对称,因此x[t]
与x[t-h]
的比较与x[t+h]
的{{1}}相同,因为我们的想法是比较 t -th值与另一个值相距 h 步骤的值。