我有一个数据集,我使用自动。Arima在预测包中做预测。我想把预测结果保存到一个矩阵中,以便进一步分析。在使用以下代码行进行分析之前,我创建了一个空数据帧:
predictions <- data.frame(point = numeric(), Lo80= numeric(), High80= numeric(), Lo95= numeric(), High95= numeric())
然后我想添加(例如)的预测。
model <- auto.arima(x)
pred <- forecast(model,h=7,l=c(50,80))
for (kk in 1:7){
predictions <- rbind(predictions, data.frame(point = pred[kk,1], Lo50= pred[kk,2], High50= pred[kk,3], Lo80= pred[kk,4], High80= pred[kk,5])}
然而,我得到一个错误,原因是pred值不是矩阵格式,没有为它们定义维度。无论如何,我可以存储我想要的结果。"pred"的示例如下:
Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
53.14286 19.84681 6.384444 33.30918 -0.742103
53.42857 24.46304 10.016633 38.90945 2.369167
library(forecast)
model <- auto.arima(x)
pred <- as.data.frame(forecast(model, h=7, level=c(50,80)))
summary(pred)
会给你一个你想要的data.frame
,你可以用as.matrix
把它转换成一个矩阵。
library(forecast)
model <- auto.arima(WWWusage)
pred <- forecast(model, h=7, level=c(50,80))
predDF <- summary(pred)
predMatrix <- as.matrix(predDF)