ggplot2绘制 arima诊断图

library(ggfortify)
autoplot(acf(gold[,2], plot = FALSE))
ggtsdiag(auto.arima(gold[,2]))

将数据改为时间格式

gold <- as.xts(gold[, 2], order.by = gold[, 1])

设置时间格式 绘制时间趋势图

gdp <- ts((0.001 * (gdp[, 2])), frequency = 4, start = c(1992, 1))
plot(gdp)

每年每月图

monthplot(gdp, xlab="", ylab="")

每年每季度图

seasonplot(gdp, s = 4, year.labels = T, col = rainbow(12))