0

少しロジスティック回帰を行った後、ラッソ回帰を実行しようとしていますが、さまざまなコマンドを入力すると、オブジェクトが見つからないというエラーが発生し続けますこれまでのコードは次のとおりです

data <- read.csv("ahw.csv", 
+                  colClasses = c("factor",
+                                 rep("numeric", 15)))

それは私に与えます

 $ PlacedN: Factor w/ 2 levels "0","1": 1 1 1 1 1 2 1 1 1 1 ...
 $ TrA    : num  9 13 11 9 -1 13 10 10 10 6 ...
  $ JoA    : num  0 2 0 10 6 15 8 9 9 6 ...
  $ aPr    : num  0 0 0 0 0 0 0 0 0 0 ...
  $ bPr    : num  0 1 0 0 1 0 0 0 0 0 ...
 $ mPr    : num  0 0 0 0 1 0 0 0 0 0 ...
 $ Vdw    : num  0 0 0 0 0 0 0 0 0 0 ...
 $ ALL    : num  32 46 15 16 48 50 32 9 28 4 ...
 $ COD    : num  -1 80 0 25 -1 50 47 -1 -1 -1 ...
 $ DIS.   : num  32 46 23 15 48 50 32 9 28 0 ...
 $ cIM    : num  -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 ...
 $ cFE    : num  -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 ...
  $ LAT    : num  38 38 9 20 61 58 37 12 25 5 ...
  $ CLA    : num  40 50 9 25 -1 75 61 -1 -1 -1 ...
  $ CLC    : num  -1 -1 0 25 -1 75 61 -1 -1 -1 ...
  $ LAC    : num  38 50 9 25 -1 60 61 -1 -1 -1 ...

明らかに、次を使用してロジスティック式を呼び出すことができます

est <- glm(PlacedN ~ ., data=data, family="binomial")
summary(est)

したがって、基本的にLASSO（ロジットを使用）を実行し、その後に係数と素敵なグラフを取得するための次のステップは何ですか

私のような初心者を混乱させるパッケージがいくつかあります

前もって感謝します

マイケル

4

1 に答える 1

2

パッケージはglmnet役に立つはずです。著者による素晴らしいチュートリアルがあります。ただし、ここでコードを使用して簡単に開始します。

require(glmnet)
est <- glmnet(as.matrix(data[,2:16]), data$PlacedN , family="binomial")
summary(est)
plot(est)
last <- dim(coef(est))[2]
coef(est)[last]

お役に立てれば！

于 2014-12-29T17:19:27.343 に答える