一组时间序列数据,通过了adf检验,纯随机性检验,但是自相关图是这样的
根据提供的自相关图可以看出,时间序列数据存在明显的季节性,因此可以考虑对其进行季节性差分。然后再进行ADF检验和白噪声检验,直到满足条件为止。具体代码如下:
# 载入数据
data <- read.csv("data.csv", header = TRUE, sep = ",")
ts_data <- ts(data[,2], start = c(2016, 1), frequency = 12)
# 季节性差分
diff_data <- diff(ts_data, differences = 1, 12)
# ADF检验
library(tseries)
adf.test(diff_data)
# 白噪声检验
library(lmtest)
Box.test(diff_data, type = "Ljung-Box")
# 自动定阶
library(forecast)
auto.arima(diff_data)
根据运行结果,可以得到最终的定阶结果为ARIMA(0,1,1)(0,1,1)[12],即进行一阶季节性差分,使用一阶差分和一阶季节性差分的组合,同时考虑季节性影响。