minrss.txt

这个小程序可以简单的选取回归分析中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)