ARIMA预测一直出现错误，提示“数据”必须是向量类型，而不是“NULL”。

Question

在将我的ARIMA拟合到数据时，我一直遇到一个错误，即“data”必须是向量类型，却是“NULL”。

library(forecast)

foo <- read.csv("https://nofile.io/g/0qrJl41nhf3bQQFjBmM6JurzGJFQSioCTGEzZhWVl9zA1kXnAJsCsSsxN1ZN7F4D/data.csv/")

data <- data.frame(year, Car)
data <- ts(data[,2],start = c(1990,1),frequency = 1)

plot(data)
plot(diff(data),ylab='Differenced Car Usage')
plot(log10(data),ylab='Log (Car Usage)')
plot(diff(log10(data)),ylab='Differenced Log (Tractor Sales)')
par(mfrow = c(1,2))
acf(ts(diff(log10(data))),main='ACF Tractor Sales')
pacf(ts(diff(log10(data))),main='PACF Tractor Sales')

require(forecast)
ARIMAfit <- auto.arima(log10(data), approximation=FALSE,trace=FALSE)
summary(ARIMAfit)

par(mfrow = c(1,1))
pred <- predict(ARIMAfit, n.ahead = 3)

array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), : 'data'必须是向量类型，而不是'NULL'

我简单地不明白我做错了什么，如果有人看到问题，我会感激任何帮助。 谢谢 -MF

- farmerm3

---

你在哪里定义了 Car？我指的是 data <- data.frame(year, Car) 这一行。 - MKR

你的代码中没有使用过 foo。 - MKR

Answer 1

library(forecast)
foo <- read.table(file="data.csv", header=T, sep=",")
data <- ts(foo$Car,start = c(1990,1),frequency = 1)

# Use 'forecast' to get predition from the model estimated by 'auto.arima'
ARIMAfit1 <- auto.arima(log10(data), approximation=T, trace=FALSE, allowdrift=F)
summary(ARIMAfit1)
forecast(ARIMAfit1, h = 3)

#      Point Forecast    Lo 80    Hi 80     Lo 95    Hi 95
# 2017       1.415713 1.165870 1.665556 1.0336109 1.797815
# 2018       1.415713 1.128307 1.703119 0.9761635 1.855262
# 2019       1.415713 1.095115 1.736310 0.9254014 1.906024


# The same model estimated using 'arima'
# Here you can use 'predict'
ARIMAfit2 <- arima(log10(data), order=c(0,1,1))
summary(ARIMAfit2)
predict(ARIMAfit2, n.ahead=3)

# $pred
# Time Series:
# Start = 2017 
# End = 2019 
# Frequency = 1 
# [1] 1.415713 1.415713 1.415713
# $se
# Time Series:
# Start = 2017 
# End = 2019 
# Frequency = 1 
# [1] 0.1911677 0.2199090 0.2453055

- Marco Sandri

---

请问您能否详细解释一下您所设置的参数，错误背后的原因以及如何使用这些参数来解决它？ - InVinci

已经过去一年多了，但我仍然没有得到答案。 - InVinci

嗨@InVinci。错误的原因是使用predict来处理auto.arima返回的对象。正确的命令是forecast。predict命令应该用于由arima估计的模型。我答案中使用的数据集只是一个示例，与解决方案无关。 - Marco Sandri

即使使用了预测，我仍然面临这个问题。我从未使用过预测。 - InVinci

@InVinci 请在 Stack Overflow 上开一个新的问题，并描述您的问题。 - Marco Sandri