如何使脊回归八度

Question

我有八度音阶4.2.0并且知道它没有实现这个功能。我尝试制作一个代码，但它很有效。你能帮我解决一下吗？还是有现成的代码，我可以使用吗？

我的代码：

matrix=xlsread('path'); %predictors
Y=xlsread('path');      
[K,J]=size(matrix); 
mid=mean(matrix);      
standard=std(matrix);
Jn=J-1; 
for len= 1:K
for ll=1:J
Z(len,ll)=matrix(len,ll)-mid(ll))/standard(ll);
end
end

I=eye(J,J);
n=1;

for kfact=0:0.0001:0.5

teta(n,:)=inv(matrix'*matrix+kfakt*I)*matrix'*Y;
k(n)=kfakt;
n=n+1;
end

有人有个主意，因为我无法找到如何制作它？

我试着指出我的问题。在本书中是参数的值，它们以k ver small（例如0.02）的依赖性计算。我得到的值大约是20.我不明白为什么？

谢谢！