24小时热门版块排行榜    

查看: 2208  |  回复: 0

maitian9875

金虫 (小有名气)

[求助] 求计算AIC的程序

各位大侠,本人手中有几组模型的预测数据,每一组中有实验值和预测值。现想通过AIC的计算来看看这几组模型的预测效果。但是小弟才疏学浅,不会用Matlab编AIC计算程序。希望哪位大侠帮我编一个Matlab计算AIC的程序。之前也编过,但是发现这个AIC随n值增大而增大,不收敛,估计是残差平方和(RSS)计算错误。求指点啊!!!【附:之前编的Matlab程序】
function AIC=SEP1
RSS=0;
int m;
int n;
K=5;
A=xlsread('114.xls');
[m,n]=size(A);
if m>=200;
    for  i=1:m;
        b(i)=(A(i,1)-A(i,2))/(A(i,2));
        RSS=RSS+b(i)^2;
    end
    AIC=2*K+m*log(RSS/m);
else
    for  i=1:m;
        b(i)=A(i,1)-A(i,2);
        RSS=RSS+b(i)^2;
    end
    AIC=2*K+m*log(RSS/m)+2*K*(K+1)/(m-K-1);
end

$BMH`VO0~O]%)GYHWZHH7R7.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 maitian9875 的主题更新
信息提示
请填处理意见