?? minrss.txt
字號:
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)
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -