p = 100
n = 100
beta1 = numeric(p)
beta1[c(1,5,10)] = c(1,1,1)
x = matrix(0,ncol=p,nrow=n)
xtest = matrix(0,ncol=p,nrow=n)
for(j in 1:p){
x[,j] = rnorm(n,0,1)
xtest[,j] = rnorm(n,0,1)
}
y = x %*% beta1 + rnorm(n)
p1 = regRSM(y,x)
predict(p1,xtest)Run the code above in your browser using DataLab