|
|
[求助]
求计算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 |
» 猜你喜欢
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有9人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有8人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有9人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
OpenAI又攻克数学难题了,离AI做有机合成还有多远?
已经有6人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
|