matlab很小的数,matlab中aicbic确定阶数的太小

想用aicbic确定ARMA的阶数,但用这个确定的阶数时会出现如下错误

Error using arima/validateModel (line 1314)

The non-seasonal moving average polynomial is non-invertible.

Error in arima/setLagOp (line 391)

Mdl = validateModel(Mdl);

Error in arima/estimate (line 1183)

Mdl = setLagOp(Mdl, 'MA' , LagOp([1  coefficients(iMA)' ], 'Lags', [0 LagsMA ]));

解决方法:

论坛有人说把p,q的阶数设置小点,我有300个数据呢,最后阶数只能设置为pMax=5; qMax=2; pMax设置为6都出错,如果这样的话,就感觉太小了,而且和eviews自相关函数和偏自相关函数的结果相差好大,能有人帮看下原因吗,谢谢

代码如下:

LogL = zeros(pMax,qMax);%pMax,qMax来自于自相关、偏自相关序列图显著阶数最(较)大值

SumPQ = LogL;

for p = 1:pMax

for q = 1:qMax

Mdl = arima(p,0,q);

[fit,~,LogL(p,q)] = estimate(Mdl,X,'print',false);%X是待检验序列

SumPQ(p,q) = p+q;

end

end

logL = reshape(LogL,pMax*qMax,1);

numParams = reshape(SumPQ,pMax*qMax,1) + 1;

[~,bic] = aicbic(logL,numParams,500);

BIC = reshape(bic,pMax,qMax);

minBIC = min(bic)

[BICbestP,BICbestQ] = find(BIC == minBIC)

Mdl = arima(BICbestP,0,BICbestQ);

%[fit,~,~] = estimate(Mdl,lnR,'print',false)

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值