我正在使用以下r函数来拟合分位数套索。
library(quantreg)
p = 100
n = 50
beta0=c(1, 0, -1, 0, -2, rep(0,35), -2, rep(0,35), -1, rep(0,23))
X = matrix(rnorm(n * p, 0,1), ncol = p)
y=vector(mode="numeric",length=n)
ei_N=rnorm(n,0,1) # the model errors
y= X%*%beta0+ei_N
q.lasso = rq(y ~ X, method="lasso",lambda = 30)
q.scad = rq(y ~ X, method="scad",lambda = 30)
使用method = scad时，出现以下错误消息。
Error in rq.fit.fnb(x, y, tau = tau) :
Error info = 51 in stepy: singular design
我还在同一rq.fit.scad包上使用了r。但是我仍然收到相同的错误消息
q.scad = rq.fit.scad(X, y, tau, alpha = 3.7, lambda = seq(1, 100, length.out = 100), start = "rq", beta = 0.9995, eps = 1e-06)
这是什么原因？预先非常感谢。