[求助]写了段R code,老有错，不知道怎么回事

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

> aic<-matrix(rep(0,36),6,6)
> for(i in 1:6) for(j in 1:6)
+ {fit.arima<-arima(WR,order=c(i,0,j))
+ aic[i,j]<-fit.arima$aic
+ }
Warning messages:
1: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
2: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
3: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
4: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
5: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
6: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
7: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1
8: In arima(WR, order = c(i, 0, j)) :
  possible convergence problem: optim gave code=1

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：code ODE COD 不知道 ARIMA code

你可能没用acf & pacf初判

ARIMA model可能的p q order.

不收敛

Explain of error infor:

ARMA(p,q), given p,q, then it generally uses m.l.e. to estimate the model parameters(phi_hat,theta_hat, and white noise variance),  probably(you should check maybe), I think they use the R base built-in optim() to get the nonlinear optimization solution. In optim(), if flag returns 1, which means that the algorithm does Not converge, in another words, it's not the m.l.e. estimate you want.

Hope this helps.

这个意思是36个模型中有8个不收敛？