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

CC = matrix(c(1, -1, rep(0, 10), 1, -1, rep(0, 10), 1, -1, rep(0, 4)), nrow = 3, byrow = TRUE); 
dd = c(0, 1, -1); 
EE = matrix(c(rep(0, 6), 1, -1, rep(0, 10), 1, -1), nrow = 2, byrow = TRUE); 
ff = c(1, 2);

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 for quantreg.genlasso.fit.gurobi()
######################################################################

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

fit1 = genlasso.fit.gurobi(xx, yy, lambda = 2, 
	 		Cmat = CC, dvec = dd, Emat = EE, fvec = ff);

######################################################################
# example for quantreg.genlasso.fit.cvxr()
######################################################################

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

fit2 = genlasso.fit.cvxr(xx, yy, lambda = 2,  
			Cmat = CC, dvec = dd, Emat = EE, fvec = ff);

fit3 = genlasso.fit.cvxr(xx, yy, lambda = 2, 
			Cmat = CC, dvec = dd, Emat = EE, fvec = ff, solver = "GUROBI");

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