%%AIC准则定阶
maxLags = 4;
AICSet = zeros(maxLags,maxLags);
parfor i = 1:maxLags%使用并行工具箱进行计算
for j=1:maxLags
mdl=arima('ARLags',[1:i],'MALags',[1,j]);
[EstMdl, EstParamCov, logL, info]=estimate(mdl, pct,'display','off');
AICSet(i,j)=aicbic(logL,length(info.X));
end
end
在进行完定阶时报错
Nonseasonal moving average lags must be unique, positive integers.
请帮忙