套索GLM

Question

後之後做了一些迴歸林嘗試做套索迴歸，但上鍵入我不斷收到找不到對象錯誤在這裏的各種命令是到目前爲止我的代碼

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) 

所以basical後勤公式LY什麼是我下一步要做LASSO（與Logit模型），並之後得到的係數和任何漂亮的圖形它

有幾個包在那裏被預先

混淆初學者和我一樣

謝謝

邁克爾

Answer 1

該glmnet包應該是有用的。作者有一個偉大的tutorial。但是，這裏是使用你的代碼的快速入門。

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] 

希望這有助於！