# load data set
data(GLESsmall)
# define response and covariate matrix
X <- scale(GLESsmall[, 11:14])
Y <- as.matrix(GLESsmall[, 1:10])
# vector of subtitles, containing the coding of the single covariates
subs <- c("(in years)","female (1); male (0)",
"East Germany (1); West Germany (0)","(very) good (1); else (0)")
# vector of tuning parameters
lambda <- exp(seq(log(31),log(1),length=50))-1
# compute BTLLasso model, increase accuracy of the algorithm
m <- BTLLasso(Y = Y, X = X, lambda = lambda, control = BTLLasso.ctrl(epsilon = 1e-05))Run the code above in your browser using DataLab