我正在使用 R 包glmnet在具有大约 60'000 个二进制预测变量的二进制响应上拟合套索 logit 模型(alpha=1)。
根据http://www.stat.cmu.edu/~ryantibs/papers/lassounique.pdf中的结果,该解决方案不应该是唯一的,但是,通过glmnet和cv.glmnet的估计结果都很好,没有问题。
所以我想知道glmnet如何处理非唯一解决方案。它是否计算无限线性预测变量之一(因为即使使用非唯一解也是唯一的)并将其作为解返回?
我正在使用 R 包glmnet在具有大约 60'000 个二进制预测变量的二进制响应上拟合套索 logit 模型(alpha=1)。
根据http://www.stat.cmu.edu/~ryantibs/papers/lassounique.pdf中的结果,该解决方案不应该是唯一的,但是,通过glmnet和cv.glmnet的估计结果都很好,没有问题。
所以我想知道glmnet如何处理非唯一解决方案。它是否计算无限线性预测变量之一(因为即使使用非唯一解也是唯一的)并将其作为解返回?