使用数据帧时从ARIMA预测特征值时出错

时间:2018-04-17 05:36:10

标签: r

我是一名ARIMA模特

> data=read.csv("R1.csv")
> resp=data$RESP
> respNA=resp[1:100]
> model=arima(respNA,order = c(0,0,1),include.mean = FALSE)
> p=predict(model$model,na.ahead=10)

但是对于最后一行,我得到以下错误:

Error in ets(object, lambda = lambda, biasadj = biasadj, allow.multiplicative.trend = allow.multiplicative.trend,  : 
  y should be a univariate time series

上面的代码出了什么问题?

当我调用print(respNA)时,我猜数据是coorect,输出如下:

  [1] 0.00000000 0.18333333 0.93333333 0.00000000 0.25000000 0.86666667 0.00000000 0.31666667 0.80000000 0.00000000
 [11] 0.38333333 0.73333333 0.00000000 0.45000000 0.66666667 0.00000000 0.51666667 0.60000000 0.00000000 0.58333333
 [21] 0.53333333 0.00000000 0.65000000 0.46666667 0.00000000 0.71666667 0.40000000 0.00000000 0.78333333 0.33333333
 [31] 0.00000000 0.85000000 0.26666667 0.00000000 0.91666667 0.20000000 0.00000000 0.98333333 0.13333333 0.05000000
 [41] 1.00000000 0.06666667 0.11666667 1.00000000 0.00000000 0.18333333 0.93333333 0.00000000 0.25000000 0.86666667
 [51] 0.00000000 0.31666667 0.80000000 0.00000000 0.38333333 0.73333333 0.00000000 0.45000000 0.66666667 0.00000000
 [61] 0.51666667 0.60000000 0.00000000 0.58333333 0.53333333 0.00000000 0.65000000 0.46666667 0.00000000 0.71666667
 [71] 0.40000000 0.00000000 0.78333333 0.33333333 0.00000000 0.85000000 0.26666667 0.00000000 0.91666667 0.20000000
 [81] 0.00000000 0.98333333 0.13333333 0.05000000 1.00000000 0.06666667 0.11666667 1.00000000 0.00000000 0.18333333
 [91] 0.93333333 0.00000000 0.25000000 0.86666667 0.00000000 0.31666667 0.80000000 0.00000000 0.38333333 0.73333333

1 个答案:

答案 0 :(得分:1)

predict的正确语法是predict(your_model, ...),即在您的情况下

# Fit ARIMA model
model <- arima(respNA,order = c(0,0,1),include.mean = FALSE);

# Predict next 10 steps
p <- predict(model, n.ahead = 10);
p;
$pred
Time Series:
Start = 101
End = 110
Frequency = 1
 [1] 0.2006601 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
 [8] 0.0000000 0.0000000 0.0000000

$se
Time Series:
Start = 101
End = 110
Frequency = 1
 [1] 0.4155428 0.5847784 0.5847784 0.5847784 0.5847784 0.5847784 0.5847784
 [8] 0.5847784 0.5847784 0.5847784

样本数据

respNA <- c(0.00000000, 0.18333333, 0.93333333, 0.00000000, 0.25000000, 0.86666667, 0.00000000, 0.31666667, 0.80000000, 0.00000000,
0.38333333, 0.73333333, 0.00000000, 0.45000000, 0.66666667, 0.00000000, 0.51666667, 0.60000000, 0.00000000, 0.58333333,
0.53333333, 0.00000000, 0.65000000, 0.46666667, 0.00000000, 0.71666667, 0.40000000, 0.00000000, 0.78333333, 0.33333333,
0.00000000, 0.85000000, 0.26666667, 0.00000000, 0.91666667, 0.20000000, 0.00000000, 0.98333333, 0.13333333, 0.05000000,
1.00000000, 0.06666667, 0.11666667, 1.00000000, 0.00000000, 0.18333333, 0.93333333, 0.00000000, 0.25000000, 0.86666667,
0.00000000, 0.31666667, 0.80000000, 0.00000000, 0.38333333, 0.73333333, 0.00000000, 0.45000000, 0.66666667, 0.00000000,
0.51666667, 0.60000000, 0.00000000, 0.58333333, 0.53333333, 0.00000000, 0.65000000, 0.46666667, 0.00000000, 0.71666667,
0.40000000, 0.00000000, 0.78333333, 0.33333333, 0.00000000, 0.85000000, 0.26666667, 0.00000000, 0.91666667, 0.20000000,
0.00000000, 0.98333333, 0.13333333, 0.05000000, 1.00000000, 0.06666667, 0.11666667, 1.00000000, 0.00000000, 0.18333333,
0.93333333, 0.00000000, 0.25000000, 0.86666667, 0.00000000, 0.31666667, 0.80000000, 0.00000000, 0.38333333, 0.73333333)