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 <- 2
tau <- 0.5
sol <- quantreglp(y,X,tau,lambda)
Run the code above in your browser using DataLab