套索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] 

希望这有助于！