# NOT RUN {
set.seed(12)
X=matrix(data=sample(1:3,1200*5, replace = TRUE) ,nrow=1200,ncol=5)
b=sample(1:3,5, replace = TRUE)
e= rnorm(1200, 0, 1)
Y=X%*%b+e
alpha=0.05
K=10
nk=1200/K
Opt1(X,Y,alpha,K,nk)
# }
Run the code above in your browser using DataLab