## Not run: ------------------------------------
# data(GLESsmall)
#
# ## extract data and center covariates for better interpretability
# Y <- GLESsmall$Y
# X <- scale(GLESsmall$X, scale = FALSE)
# Z1 <- scale(GLESsmall$Z1, scale = FALSE)
#
# ## vector of subtitles, containing the coding of the X covariates
# subs.X <- c("", "female (1); male (0)")
#
# ## vector of tuning parameters
# lambda <- exp(seq(log(61), log(1), length = 30)) - 1
#
#
# ## compute BTLLasso model
# m.gles <- BTLLasso(Y = Y, X = X, Z1 = Z1, lambda = lambda)
# print(m.gles)
#
# singlepaths(m.gles, x.axis = "loglambda", subs.X = subs.X)
# paths(m.gles, y.axis = "L2")
#
# ## Cross-validate BTLLasso model
# m.gles.cv <- cv.BTLLasso(Y = Y, X = X, Z1 = Z1, lambda = lambda)
# print(m.gles.cv)
#
# singlepaths(m.gles.cv, subs.X = subs.X)
## ---------------------------------------------
Run the code above in your browser using DataLab