minrss.asv

这个小程序可以简单的选取回归分析中box-cox变换的lambda值

function f=minRSS(x,y)
minRSS=M;%初始RSS值为一极大数
bestlambda=0;%初始lambda
for lambda=0:0.05:1;
transdat=boxcox(lambda,y);
J=1;
for i=1:length(y)  %这道题有53个数据，这个循环是为了将53个数据连乘而构造的
    J=J*(y(i)^(lambda-1));
end
    z=transdat./J^(1/53);
    E=eye(53);
    RSS=z'*(E-x*inv(x'*x)*x')*z;
    if RSS<minRSS