各位大牛们,希望能帮个忙啊,最近在看chi2gof函数这块,可是不知道为什么返回的检验值P总为NaN,希望大牛们能帮我解决下,谢谢!
以下是我的程序:
clc,clear
a=[17.98 16.92 16.78 16.85 16.85 16.97 16.92 17.54 16.49 17.08 17.26 17.16 16.90];
a=nonzeros(a')'; %按照原来数据的顺序去掉零元素
da=diff(a);%消除趋势性的差分运算
daa=diff(da);%消除趋势性的差分运算
w=length(daa); %计算差分后的数据个数
n=20;%预报的数据个数
for i=0:3
for j=0:3
spec= garchset('R',i,'M',j,'Display','off'); %指定模型的结构
[coeffX,errorsX,LLFX] = garchfit(spec,daa); %拟合参数
num=garchcount(coeffX); %计算拟合参数的个数
%compute Akaike and Bayesian Information Criteria
[aic,bic]=aicbic(LLFX,num,w);
fprintf('R=%d,M=%d,AIC=%f,BIC=%f\n',i,j,aic,bic); %显示计算结果
end
end
r=input('输入阶数R=