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))