2017-05-15 19 views
0

リッジ回帰のラムダ行列を抽出したいと思います。私はそれをプロットすることができますが、私は値を抽出することはできません。必要な行列には、ラムダ(行内)と名前を含む変数(列)が含まれている必要があります。glmnetのラムダ行列を抽出する

library(glmnet) 

x1<-rnorm(200,5,4) 
x2<-rnorm(200,45,3) 
x3<-rnorm(200,-10,4) 
y<-80+x2*3+rnorm(200,0,1) 
X<-cbind(x1,x2,x3) 

l.vec<-10^seq(10,-2,length=100) 

my.ridge<-glmnet(x=X,y=y,alpha=0, 
       lambda=l.vec, 
       standardize=TRUE) 

plot(my.ridge,xvar="lambda") 

どのようにしてプロット情報をマトリックスにリダイレクトできますか?

Thxを&種類を使用すると、係数行列を意味するか

答えて

1

に関しては?

coef(my.ridge) 
+0

あなたが正しいです。これを転置する必要があるのは:-( cbind(t(coef(my.ridge))、lambda = my.ridge $ lambda)[1:4、] –

関連する問題