r - glmnet 错误 (nulldev == 0) 停止 ("y is constant; gaussian glmnet fails at standardization step")

Question

我在 R 中使用 glmnet 运行以下(截断的)代码

# do a lot of things to create the design matrix called x.design

> glmnet(x.design, y, thresh=1e-11)

其中 x.design 是 n x p 设计矩阵，其中 n > p 和 y 是使用核密度估计获得的 n x 1 响应向量。 x.design 和 y 都包含真实的条目。运行代码时收到以下错误消息:

Error in if (nulldev == 0) stop("y is constant; gaussian glmnet fails at 
standardization step") : missing value where TRUE/FALSE needed 

我访问过并阅读过

Running glmnet package in R, getting error "missing value where TRUE/FALSE needed", maybe due to missing values?

但是我想不出解决我的问题的方法。

有人可以提出解决方案吗？

Answer 1

看来您的响应向量 y是恒定的。 GLMNET 尝试对其进行标准化(可能减去均值，然后除以当前的标准差)，但不能，因为标准差为 0。打印 y 及其方差以确保。
您还应该检查您的内核估计程序。