######################################################################
# examples
######################################################################

xx = matrix(rnorm(400 * 10), nrow = 400, ncol = 10); 
bb = c(1, -1, 0.5, 2, rep(0, 6));
yy = xx %*% bb + rnorm(400, 0.3); 

######################################################################
# example using glmnet
######################################################################

library(glmnet)
fit.glmnet = glmnet(xx, yy, family = "gaussian", alpha = 1, 
                    standardize = FALSE, intercept = FALSE);
b.glmnet = as.matrix(coef(fit.glmnet, s = 1.2 / nrow(xx), exact = TRUE, x = xx, y = yy))[-1];

######################################################################
# example for lasso.fit.cvxr()
######################################################################

source("http://www.auburn.edu/~zengpen/papers/Rcodes/lasso_fit_cvxr.R");

fit1 = lasso.fit.cvxr(xx, yy, lambda = 1.2); 
fit2 = lasso.fit.cvxr(xx, yy, lambda = 1.2, solver = "GUROBI");

######################################################################
# example for lasso.fit.gurobi()
######################################################################

source("http://www.auburn.edu/~zengpen/papers/Rcodes/lasso_fit_gurobi.R");

fit3 = lasso.fit.gurobi(xx, yy, lambda = 1.2); 

######################################################################
# THE END
######################################################################