# NOT RUN {
set.seed(33)
n = 200
p = 20
sigma = 1
x = matrix(rnorm(n*p),n,p)
beta = c(rep(3,10),rep(0,p-10))
y = x%*%beta + sigma*rnorm(n)
# run forward stepwise and predict functions
obj = fs(x,y)
fit = predict(obj,x,s=3)
# }
Run the code above in your browser using DataCamp Workspace