grp_lasso_st: Description of the function

Description Usage Arguments Value Examples

Description

Description of the function

Usage

1
grp_lasso_st(X, Y, group, nb.cores = 7, B = 500, PFER = 1)

Arguments

group

a vector with two levels. (The group of the ANCOVA)

response

a vector response variable

regressors

a quantitative matrix of regressor

a

the parameters that indicate how much the coefficients will be fused

lambda

if the user wants to use it owns values of lambdas

Value

The coefficients of the fused lasso ANCOVA for the different value of lambda

Examples

1
2
3
4
5
6
7
8
9
B <- c(1, -1, 1.5, 1.5, rep(0, 6), 2, 0, 2, 0)
group <- c(rep('M1', 10), rep('M2', 10))
regressors <- matrix(rnorm(6*20), ncol = 6)
X  <- model.matrix(~group + group:regressors - 1)
y <- X%*%B + rnorm(20)
y <- scale(y)
mod <- fl2(y, regressors, group)
colors <- c(rep("grey",2), rep('green',2),rep('black', 6), rep(c("orange","blue"), 2), 'darkgreen', rep('yellow',3), rep('purple',2))
matplot(mod$lambda ,t(mod$beta),type='l',col=colors)

Marie-PerrotDockes/ModQTL documentation built on May 16, 2019, 7:26 a.m.