set.seed(12345)
n <- 20
p <- 20
X <- matrix(rnorm(n*p),n,p)
b0 <- double(p)
k <- 4
b0[sample(1:p,k,replace=FALSE)] <- 10*rnorm(k)
y <- X%*%b0 + 0.1*rnorm(n)
lambda <- 0
tau <- c(0.05,0.25,0.5,0.75,0.95)
sol <- quantreg(y~X,tau,lambda)
coef(sol)
Run the code above in your browser using DataLab