tests/TesttuneLasso.R

#source("charge.R")
library("LINselect")

set.seed(15) # pour pouvoir comparer des exécutions sucessives
ex <- simulData(p=100,n=50,r=0.95,rSN=10)

res.tuneLasso <- tuneLasso(ex$Y,ex$X,method="Glasso",verbose=TRUE)
res.tuneLasso$Glasso$CV$support
res.tuneLasso$Glasso$Ls$support

res.tuneLasso <- tuneLasso(ex$Y,ex$X,method="lasso",verbose=TRUE)
res.tuneLasso$lasso$CV$support
res.tuneLasso$lasso$Ls$support

## Time consuming
# set.seed(5) # pour pouvoir comparer des exécutions sucessives
# ex <- simulData(p=100,n=100,r=0.95,rSN=10)
# res.tuneLasso <- tuneLasso(ex$Y,ex$X, dmax=20)
# 
# 
# set.seed(67) # pour pouvoir comparer des exécutions sucessives
# ex <- simulData(p=50,n=1000,r=0.95,rSN=10)
# res.tuneLasso <- tuneLasso(ex$Y,ex$X, dmax=20)
# res.tuneLasso$lasso$CV$coeff
# res.tuneLasso$lasso$Ls$coeff
# res.tuneLasso$Glasso$CV$support
# res.tuneLasso$Glasso$Ls$support

Try the LINselect package in your browser

Any scripts or data that you put into this service are public.

LINselect documentation built on Aug. 30, 2023, 9:10 a.m.