私はリッジ回帰に取り組んでおり、独自の関数を作成したいと考えています。以下を試してみました。k の個々の値では機能しますが、値のシーケンスの配列では機能しません。
dt<-longley
attach(dt)
library(MASS)
X<-cbind(X1,X2,X3,X4,X5,X6)
X<-as.matrix(X)
Y<-as.matrix(Y)
sx<-scale(X)/sqrt(nrow(X)-1)
sy<-scale(Y)/sqrt(nrow(Y)-1)
rxx<-cor(sx)
rxy<-cor(sx,sy)
for (k in 0:1){
res<-solve(rxx+k*diag(rxx))%*%rxy
k=k+0.01
}
最適化されたコードについても助けが必要です。