minrss.m

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

function f=minRSS(x,y)
minRSS=999999;%初始RSS值为一极大数
bestlambda=0;%初始lambda
for lambda=-2:0.005:2;%比较[-2，2]之间的lambda值
transdat=boxcox(lambda,y);
J=1;
n=length(y);
for i=1:n  %这道题有53个数据，这个循环是为了将53个数据连乘
    J=J*(y(i)^(lambda-1));%变换的Jacobi行列式
end
    z=transdat./J^(1/n);
    E=eye(n);  %单位阵
    RSS=z'*(E-x*inv(x'*x)*x')*z;
    if (RSS<minRSS)
        minRSS=RSS;
        bestlambda=lambda;
    end
end
fprintf('the best lambda=%g',bestlambda)