如何设置Y表示gmlnet几个预测

Question

我怎么才能找到使用岭回归模型中的截距项。

“使用岭回归具有α= 0到和λ= 0，离婚作为响应和所有其他变量作为预测”。

我知道我应该我的数据转换为矩阵模式，然后改造它以适应glmnet功能。我已经转换我对矩阵模式响应，但我不知道如何我所有的预测转化为矩阵模式，太。

set.seed(100)

require(faraway)
require(leaps)
require(glmnet)

mydata = divusa

mymodel = lm(divorce ~ year + unemployed + femlab + marriage + birth + 
military, data=mydata)
summary(mymodel)
.
.
.

y = model.matrix(divorce~.,mydata)

任何人都可以与我的X变量的代码帮助吗？我很新的R和发现它很难理解它。

Answer 1

你y = model.matrix(divorce~.,mydata)实际创建您的预测矩阵（通常称为X）。尝试

 X = model.matrix(divorce~.,mydata)
 y = mydata$divorce
 glmnet(X,y)
 glmnet(X,y,alpha=0,lambda=0)

我认为，如果你设置你实际上在做lambda=0普通的回归（即，你罚设置为零，所以岭 - > OLS）。