R语言真是博大精深
方法一
Acf(gold[,2], type = "correlation",lag.max = 100)
Acf(gold[,2], type = "partial")方法二
library(ggfortify)
autoplot(acf(gold[,2], plot = FALSE))方法三
bacf <- acf(gold[, 2], plot = FALSE)
bacfdf <- with(bacf, data.frame(lag, acf))
ggplot(data = bacfdf, mapping = aes(x = lag, y = acf)) +
geom_segment(mapping = aes(xend = lag, yend = 0)) +
geom_hline(aes(yintercept = 0.05), linetype = 2, color = 'darkblue')
# pacf
bpacf <- pacf(gold[, 2], plot = FALSE)
bpacf <- with(bpacf, data.frame(lag, acf))
ggplot(data = bpacf, mapping = aes(x = lag, y = acf)) +
geom_segment(mapping = aes(xend = lag, yend = 0)) +
geom_hline(aes(yintercept = 0.05), linetype = 2, color = 'darkblue')方法四
library("forecast")
forecast::ggtsdisplay(gold[,2],lag.max=30)方法五
ggAcf(gold[,2])
ggPacf(gold[,2])内容来源于网络如有侵权请私信删除
- 还没有人评论，欢迎说说您的想法！
客服
