如何预测仅使用从R中的glmnet获得的权重？

Question

可以使用predict(fit,newx)方法进行预测。但是如何预测我是否没有拟合对象本身，只有预测变量的权重作为矢量文件和预测变量作为矩阵文件的新观察？预测因子和结果都是连续变量。

Answer 1

你可以使用矩阵乘法，这是glmnet的作用。在预测函数中，它是：as.matrix(cbind2(1, x) %*% coefs)

示例：

library(glmnet)

x=matrix(rnorm(100*20),100,20)
y=rnorm(100)
fit1=glmnet(x,y)
coefs <- coef(fit1,s=0.01) # extract coefficients at a single value of lambda

manaul_pred <- as.matrix(cbind2(1, x) %*% coefs)
pred <- predict(fit1,newx=x,s=0.01)

manual_pred - pred # there is a negligible difference due to numeric precision