grpLassoQTL: Give the result for a group lasso

Description Usage Arguments Value Examples

Description

Give the result for a group lasso

Usage

1
2
grpLassoQTL(X, Y, marker, map, PFER = 1, B = 500, cutoff = 0.85,
  nb.cores, sep = "\\.", mrk2lg, type_group = "both")

Arguments

X

a matrix with the value of different allele for different Marker

Y

a vector or a matrix with the values of different Trait

type_group

a character indicating the way of creating the group : "both" means that the variable are selected for all a Marker for all the Trait. "marker" means that the variable are selected for a Marker but for one Trait. "trait" means that a variable are selected for all the Trait but not for all the marker but just for one allele

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.