R 中 auto.arima 模型选择的解释和再现
Interpretation and reproduction of auto.arima model selection in R
在RStudio中,我运行auto.arima上内置了航空旅客数据(AirPassengers)。数据似乎有乘法趋势。预测结果好像还算准确
# Load the corresponging library.
library(forecast)
# Save the default parameters.
defaultPar <- par(no.readonly = TRUE)
# Prepare a graph of four (2x2) subgraphs
par(mfrow = c(3, 1))
# Fit the model. Use the built in AirPassengers data.
modelAA <- auto.arima(AirPassengers)
plot(AirPassengers,
main = 'Air passengers',
col = 'purple',
ylab = 'Passengers number',
xlab = 'Year')
plot(forecast(modelAA, 24),
main = 'Air passengers + auto.arima forecasting',
col = 'red',
ylab = 'Passengers number',
xlab = 'Year')
qqnorm(modelAA$residuals, col = 'red')
qqline(modelAA$residuals, col = 'green')
# Restore the default parameters.
par(defaultPar)
如以下代码示例所示,自动 ARIMA 预测已选择以下模型:(2, 1, 1)(0, 1, 0)[12]。
print(modelAA)
Series: AirPassengers
ARIMA(2,1,1)(0,1,0)[12]
Coefficients:
ar1 ar2 ma1
0.5960 0.2143 -0.9819
s.e. 0.0888 0.0880 0.0292
sigma^2 estimated as 132.3: log likelihood=-504.92
AIC=1017.85 AICc=1018.17 BIC=1029.35
如何解释自动选择 (2, 1, 1)(0, 1, 0)[12] 并通过调用 arima(x, order = c()...) 重现它?谢谢
它是 季节性华宇模型。您可以从 enter link description here.
获取详细信息
调用 arima
function (x, order = c(0L, 0L, 0L), seasonal = list(order = c(0L, 0L, 0L), period = NA), xreg = NULL, include.mean = TRUE, transform.pars = TRUE, fixed = NULL, init = NULL, method = c("CSS-ML", "ML", "CSS"), n.cond, SSinit = c("Gardner1980", "Rossignol2011"), optim.method = "BFGS", optim.control = list(), kappa = 1e+06)
时应设置参数 order/sensonal/period
其中 period = 12
。
在RStudio中,我运行auto.arima上内置了航空旅客数据(AirPassengers)。数据似乎有乘法趋势。预测结果好像还算准确
# Load the corresponging library.
library(forecast)
# Save the default parameters.
defaultPar <- par(no.readonly = TRUE)
# Prepare a graph of four (2x2) subgraphs
par(mfrow = c(3, 1))
# Fit the model. Use the built in AirPassengers data.
modelAA <- auto.arima(AirPassengers)
plot(AirPassengers,
main = 'Air passengers',
col = 'purple',
ylab = 'Passengers number',
xlab = 'Year')
plot(forecast(modelAA, 24),
main = 'Air passengers + auto.arima forecasting',
col = 'red',
ylab = 'Passengers number',
xlab = 'Year')
qqnorm(modelAA$residuals, col = 'red')
qqline(modelAA$residuals, col = 'green')
# Restore the default parameters.
par(defaultPar)
如以下代码示例所示,自动 ARIMA 预测已选择以下模型:(2, 1, 1)(0, 1, 0)[12]。
print(modelAA)
Series: AirPassengers
ARIMA(2,1,1)(0,1,0)[12]
Coefficients:
ar1 ar2 ma1
0.5960 0.2143 -0.9819
s.e. 0.0888 0.0880 0.0292
sigma^2 estimated as 132.3: log likelihood=-504.92
AIC=1017.85 AICc=1018.17 BIC=1029.35
如何解释自动选择 (2, 1, 1)(0, 1, 0)[12] 并通过调用 arima(x, order = c()...) 重现它?谢谢
它是 季节性华宇模型。您可以从 enter link description here.
获取详细信息调用 arima
function (x, order = c(0L, 0L, 0L), seasonal = list(order = c(0L, 0L, 0L), period = NA), xreg = NULL, include.mean = TRUE, transform.pars = TRUE, fixed = NULL, init = NULL, method = c("CSS-ML", "ML", "CSS"), n.cond, SSinit = c("Gardner1980", "Rossignol2011"), optim.method = "BFGS", optim.control = list(), kappa = 1e+06)
时应设置参数 order/sensonal/period
其中 period = 12
。