Last data update: 2014.03.03

 R: multivariate Bernoulli LASSO model fitting
mvblpsR Documentation

multivariate Bernoulli LASSO model fitting

Description

fit multivariate Bernoulli LASSO model accelerated block-coordinate relaxation algorithm.

Usage

mvblps(x, y, maxOrder = 2, lambda = NULL, nlambda = 100,
       lambda.min.ratio = ifelse(nobs<nvars, .01, .0001),
       output = 0, printIter = 100, search = c('nm', 'grid'),
       tune = c("AIC", "BIC", "GACV", "BGACV"))

Arguments

x

input design matrix.

y

output binary matrix with number of columns equal to the number of outcomes per observation.

maxOrder

maximum order of interactions to be considered in outcomes.

lambda

a user specified tuning sequece. Typical usage is to have the program compute its own lambda.

nlambda

the number of lambda values, default is 100.

lambda.min.ratio

Smallest value for lambda, as a fraction of lambda.max. The default depends on the sample size nobs relative to the number of variables.

output

with values 0 or 1, indicating whether the fitting process is muted or not.

printIter

Number of iterations to be printed if output is true.

search

Tuning search approach, nm for Nelder Mead and grid for grid search.

tune

tuning approach, available methods including AIC, BIC, GACV, BGACV.

Details

The mvblps utilize the class structure of the underlying C++ code and fitted the model with accelerated block-coordinate relaxation algorithm.

Value

An object of classes mvbfit and lps, for which some methods are available.

See Also

mvbfit, unifit, stepfit, mvb.simu

Examples

# fit a simple MVB log-linear model
n <- 1000
p <- 5
kk <- 2
tt <- NULL
alter <- 1
for (i in 1:kk) {
  vec <- rep(0, p)
  vec[i] <- alter
  alter <- alter * (-1)
  tt <- cbind(tt, vec)
}
tt <- 1.5 * tt
tt <- cbind(tt, c(rep(0, p - 1), 1))

x <- matrix(rnorm(n * p, 0, 4), n, p)
res <- mvb.simu(tt, x, K = kk, rep(.5, 2))
fitMVB <- mvblps(x, res$response, output = 1)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(MVB)
Loading required package: Rcpp
Loading required package: RcppArmadillo
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/MVB/mvblps.Rd_%03d_medium.png", width=480, height=480)
> ### Name: mvblps
> ### Title: multivariate Bernoulli LASSO model fitting
> ### Aliases: mvblps
> 
> ### ** Examples
> 
> # fit a simple MVB log-linear model
> n <- 1000
> p <- 5
> kk <- 2
> tt <- NULL
> alter <- 1
> for (i in 1:kk) {
+   vec <- rep(0, p)
+   vec[i] <- alter
+   alter <- alter * (-1)
+   tt <- cbind(tt, vec)
+ }
> tt <- 1.5 * tt
> tt <- cbind(tt, c(rep(0, p - 1), 1))
> 
> x <- matrix(rnorm(n * p, 0, 4), n, p)
> res <- mvb.simu(tt, x, K = kk, rep(.5, 2))
> fitMVB <- mvblps(x, res$response, output = 1)

lambda =    1.6095   1.6095   1.6095

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.3177   1.3177   1.3177

iter 1 gpnorm = 0.0776186 nonzeros = 2(11.1111%)   alpha = 1
iter 8 gpnorm = 1.16415e-09 nonzeros = 4(22.2222%)   alpha = 51.2
Converged **

lambda =    1.0789   1.0789   1.0789

iter 1 gpnorm = 0.337801 nonzeros = 4(22.2222%)   alpha = 1
iter 7 gpnorm = 4.2075e-10 nonzeros = 4(22.2222%)   alpha = 25.6
Converged **

lambda =    0.8833   0.8833   0.8833

iter 1 gpnorm = 0.276568 nonzeros = 4(22.2222%)   alpha = 1
iter 7 gpnorm = 9.48843e-10 nonzeros = 4(22.2222%)   alpha = 64
Converged **

lambda =    0.7232   0.7232   0.7232

iter 1 gpnorm = 0.226435 nonzeros = 4(22.2222%)   alpha = 1
iter 13 gpnorm = 2.07548e-13 nonzeros = 5(27.7778%)   alpha = 4096
Converged **

lambda =    0.5921   0.5921   0.5921

iter 1 gpnorm = 0.227055 nonzeros = 5(27.7778%)   alpha = 1
iter 6 gpnorm = 4.3659e-11 nonzeros = 5(27.7778%)   alpha = 32
Converged **

lambda =    0.4848   0.4848   0.4848

iter 1 gpnorm = 0.185897 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 1.4583e-15 nonzeros = 5(27.7778%)   alpha = 64
Converged **

lambda =    0.3969   0.3969   0.3969

iter 1 gpnorm = 0.152199 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 1.09163e-09 nonzeros = 5(27.7778%)   alpha = 25.6
Converged **

lambda =    0.3249   0.3249   0.3249

iter 1 gpnorm = 0.12461 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 1.52258e-09 nonzeros = 5(27.7778%)   alpha = 25.6
Converged **

lambda =    0.2660   0.2660   0.2660

iter 1 gpnorm = 0.102022 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 1.56909e-10 nonzeros = 5(27.7778%)   alpha = 25.6
Converged **

lambda =    0.2178   0.2178   0.2178

iter 1 gpnorm = 0.0835287 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 4.9478e-16 nonzeros = 5(27.7778%)   alpha = 25.6
Converged **

lambda =    0.1783   0.1783   0.1783

iter 1 gpnorm = 0.0683875 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 2.20088e-09 nonzeros = 5(27.7778%)   alpha = 64
Converged **

lambda =    0.1460   0.1460   0.1460

iter 1 gpnorm = 0.055991 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 2.57881e-09 nonzeros = 5(27.7778%)   alpha = 25.6
Converged **

lambda =    0.1195   0.1195   0.1195

iter 1 gpnorm = 0.0458415 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 1.06121e-15 nonzeros = 5(27.7778%)   alpha = 64
Converged **

lambda =    0.0979   0.0979   0.0979

iter 1 gpnorm = 0.0375319 nonzeros = 5(27.7778%)   alpha = 1
iter 7 gpnorm = 4.64562e-09 nonzeros = 5(27.7778%)   alpha = 64
Converged **

lambda =    0.0801   0.0801   0.0801

iter 1 gpnorm = 0.0307285 nonzeros = 5(27.7778%)   alpha = 1
iter 8 gpnorm = 2.58091e-09 nonzeros = 6(33.3333%)   alpha = 128
Converged **

lambda =    0.0656   0.0656   0.0656

iter 1 gpnorm = 0.0295436 nonzeros = 6(33.3333%)   alpha = 1
iter 8 gpnorm = 4.68853e-13 nonzeros = 6(33.3333%)   alpha = 128
Converged **

lambda =    0.0537   0.0537   0.0537

iter 1 gpnorm = 0.0237844 nonzeros = 6(33.3333%)   alpha = 1
iter 12 gpnorm = 6.08119e-14 nonzeros = 7(38.8889%)   alpha = 2048
Converged **

lambda =    0.0537   0.0656   0.0656

iter 1 gpnorm = 0.0237844 nonzeros = 7(38.8889%)   alpha = 1
iter 7 gpnorm = 1.40963e-16 nonzeros = 6(33.3333%)   alpha = 64
Converged **

lambda =    0.0656   0.0537   0.0656

iter 1 gpnorm = 0.0205979 nonzeros = 6(33.3333%)   alpha = 1
iter 11 gpnorm = 1.55825e-14 nonzeros = 7(38.8889%)   alpha = 1024
Converged **

lambda =    0.0656   0.0656   0.0537

iter 1 gpnorm = 0.0237844 nonzeros = 7(38.8889%)   alpha = 1
iter 6 gpnorm = 3.21931e-11 nonzeros = 6(33.3333%)   alpha = 32
Converged **

lambda =    0.0193   0.0193   0.0193

iter 1 gpnorm = 0.0872953 nonzeros = 6(33.3333%)   alpha = 1
iter 18 gpnorm = 5.08658e-13 nonzeros = 11(61.1111%)   alpha = 52428.8
Converged **

lambda =    0.3249   0.2940   0.3249

iter 1 gpnorm = 0.554761 nonzeros = 7(38.8889%)   alpha = 1
iter 10 gpnorm = 5.182e-10 nonzeros = 5(27.7778%)   alpha = 204.8
Converged **

lambda =    0.2940   0.3249   0.3249

iter 1 gpnorm = 0.509431 nonzeros = 6(33.3333%)   alpha = 1
iter 10 gpnorm = 8.97097e-10 nonzeros = 5(27.7778%)   alpha = 204.8
Converged **

lambda =    0.3249   0.3249   0.2940

iter 1 gpnorm = 0.509431 nonzeros = 6(33.3333%)   alpha = 1
iter 10 gpnorm = 9.57014e-16 nonzeros = 5(27.7778%)   alpha = 512
Converged **

lambda =    0.0728   0.0728   0.0982

iter 1 gpnorm = 0.404907 nonzeros = 5(27.7778%)   alpha = 1
iter 15 gpnorm = 1.44368e-15 nonzeros = 7(38.8889%)   alpha = 6553.6
Converged **

lambda =    0.0036   0.0049   0.0081

iter 1 gpnorm = 0.536081 nonzeros = 5(27.7778%)   alpha = 1
alpha = 1.15292e+08 alphaMax = 1e+08

lambda =    0.1058   0.1058   0.1292

iter 1 gpnorm = 0.351734 nonzeros = 5(27.7778%)   alpha = 1
iter 10 gpnorm = 1.79403e-14 nonzeros = 5(27.7778%)   alpha = 512
Converged **

lambda =    0.3422   0.3422   0.3976

iter 1 gpnorm = 0.106425 nonzeros = 5(27.7778%)   alpha = 1
iter 6 gpnorm = 1.29267e-11 nonzeros = 5(27.7778%)   alpha = 32
Converged **

lambda =    0.6879   0.7232   0.7232

iter 1 gpnorm = 0.68726 nonzeros = 5(27.7778%)   alpha = 1
iter 15 gpnorm = 2.47064e-13 nonzeros = 5(27.7778%)   alpha = 6553.6
Converged **

lambda =    0.7232   0.6879   0.7232

iter 1 gpnorm = 0.68726 nonzeros = 5(27.7778%)   alpha = 1
iter 15 gpnorm = 2.4992e-13 nonzeros = 5(27.7778%)   alpha = 6553.6
Converged **

lambda =    1.8733   1.8733   1.7819

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    0.5235   0.5235   0.5785

iter 1 gpnorm = 0.296413 nonzeros = 5(27.7778%)   alpha = 1
iter 6 gpnorm = 3.57928e-11 nonzeros = 5(27.7778%)   alpha = 32
Converged **

lambda =    1.0522   1.0789   1.0789

iter 1 gpnorm = 0.621075 nonzeros = 5(27.7778%)   alpha = 1
iter 13 gpnorm = 9.49843e-09 nonzeros = 4(22.2222%)   alpha = 4096
Converged **

lambda =    1.0789   1.0522   1.0789

iter 1 gpnorm = 0.621075 nonzeros = 5(27.7778%)   alpha = 1
iter 14 gpnorm = 4.04247e-15 nonzeros = 4(22.2222%)   alpha = 3276.8
Converged **

lambda =    0.7422   0.7422   0.8000

iter 1 gpnorm = 0.694103 nonzeros = 5(27.7778%)   alpha = 1
iter 14 gpnorm = 1.15986e-10 nonzeros = 5(27.7778%)   alpha = 8192
Converged **

lambda =    1.4642   1.5782   1.7884

iter 1 gpnorm = 0.0619347 nonzeros = 0(0%)   alpha = 1
iter 9 gpnorm = 2.00713e-13 nonzeros = 3(16.6667%)   alpha = 256
Converged **

lambda =    1.1645   1.1645   1.2242

iter 1 gpnorm = 0.141152 nonzeros = 4(22.2222%)   alpha = 1
iter 7 gpnorm = 1.67578e-09 nonzeros = 4(22.2222%)   alpha = 64
Converged **

lambda =    1.0929   1.0929   1.1347

iter 1 gpnorm = 0.598424 nonzeros = 5(27.7778%)   alpha = 1
iter 18 gpnorm = 2.00388e-15 nonzeros = 4(22.2222%)   alpha = 131072
Converged **

lambda =    1.3013   1.3177   1.3177

iter 1 gpnorm = 0.345138 nonzeros = 4(22.2222%)   alpha = 1
iter 12 gpnorm = 1.20381e-14 nonzeros = 4(22.2222%)   alpha = 2048
Converged **

lambda =    1.3177   1.3013   1.3177

iter 1 gpnorm = 0.345138 nonzeros = 4(22.2222%)   alpha = 1
iter 15 gpnorm = 3.68004e-11 nonzeros = 4(22.2222%)   alpha = 16384
Converged **

lambda =    6.9696   6.9696   6.8830

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.7368   1.7368   1.7807

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.4472   1.4563   1.4563

iter 1 gpnorm = 0.201214 nonzeros = 4(22.2222%)   alpha = 1
iter 9 gpnorm = 5.76972e-09 nonzeros = 3(16.6667%)   alpha = 102.4
Converged **

lambda =    1.4563   1.4472   1.4563

iter 1 gpnorm = 0.201214 nonzeros = 4(22.2222%)   alpha = 1
iter 10 gpnorm = 1.89358e-08 nonzeros = 3(16.6667%)   alpha = 204.8
Converged **

lambda =    1.3263   1.3263   1.3514

iter 1 gpnorm = 0.33002 nonzeros = 4(22.2222%)   alpha = 1
iter 13 gpnorm = 8.51161e-09 nonzeros = 4(22.2222%)   alpha = 1638.4
Converged **

lambda =    6.5532   6.6772   6.8892

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.1211   2.1211   2.1477

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.4610   1.4610   1.4748

iter 1 gpnorm = 0.190562 nonzeros = 4(22.2222%)   alpha = 1
iter 10 gpnorm = 9.67154e-09 nonzeros = 3(16.6667%)   alpha = 512
Converged **

lambda =    1.5262   1.5310   1.5310

iter 1 gpnorm = 0.0789629 nonzeros = 3(16.6667%)   alpha = 1
iter 5 gpnorm = 5.6824e-16 nonzeros = 2(11.1111%)   alpha = 16
Converged **

lambda =    1.5310   1.5262   1.5310

iter 1 gpnorm = 0.0746662 nonzeros = 3(16.6667%)   alpha = 1
iter 5 gpnorm = 4.61591e-16 nonzeros = 2(11.1111%)   alpha = 16
Converged **

lambda =    8.4126   8.4918   8.6256

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.3440   2.3440   2.3587

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.5334   1.5334   1.5407

iter 1 gpnorm = 0.072424 nonzeros = 3(16.6667%)   alpha = 1
iter 5 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 16
Converged **

lambda =    1.5673   1.5697   1.5697

iter 1 gpnorm = 5.6824e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.92726e-16 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    1.5697   1.5673   1.5697

iter 1 gpnorm = 4.61591e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 7.2073e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    10.2242   10.2242   10.2083

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.4641   2.4641   2.4718

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.5882   1.5895   1.5895

iter 1 gpnorm = 4.92726e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.9394e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.5895   1.5882   1.5895

iter 1 gpnorm = 7.2073e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.94778e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.5710   1.5710   1.5747

iter 1 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    10.5080   10.5080   10.4997

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.5265   2.5265   2.5304

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.5988   1.5994   1.5994

iter 1 gpnorm = 4.9394e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.94778e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.5994   1.5988   1.5994

iter 1 gpnorm = 4.94778e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 3.91034e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.5901   1.5901   1.5920

iter 1 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    10.6528   10.6528   10.6486

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.5582   2.5582   2.5602

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.6041   1.6044   1.6044

iter 1 gpnorm = 4.94778e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.23315e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6044   1.6041   1.6044

iter 1 gpnorm = 3.91034e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 3.70234e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.5998   1.5998   1.6007

iter 1 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    10.7259   10.7259   10.7238

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.5742   2.5742   2.5752

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.6068   1.6070   1.6070

iter 1 gpnorm = 4.23315e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 5.56405e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6070   1.6068   1.6070

iter 1 gpnorm = 3.70234e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 3.70234e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6046   1.6046   1.6051

iter 1 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    10.7627   10.7627   10.7616

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.5823   2.5823   2.5828

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.6081   1.6082   1.6082

iter 1 gpnorm = 5.56405e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.30242e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6082   1.6081   1.6082

iter 1 gpnorm = 3.70234e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 4.9394e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6070   1.6070   1.6073

iter 1 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    10.7811   10.7811   10.7806

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    2.5863   2.5863   2.5866

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **

lambda =    1.6088   1.6088   1.6088

iter 1 gpnorm = 4.30242e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 3.98522e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6088   1.6088   1.6088

iter 1 gpnorm = 4.9394e-16 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 7.64893e-16 nonzeros = 2(11.1111%)   alpha = 0.8
Converged **

lambda =    1.6082   1.6082   1.6084

iter 1 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 1
iter 2 gpnorm = 1.80728e-10 nonzeros = 2(11.1111%)   alpha = 2
Converged **

lambda =    1.6095   1.6095   1.6095

iter 1 gpnorm = 0.0471699 nonzeros = 0(0%)   alpha = 1
iter 6 gpnorm = 2.17987e-11 nonzeros = 2(11.1111%)   alpha = 32
Converged **
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>


Created & Maintained by Osamu Ogasawara (osamu.ogasawara@gmail.com) and IMSBIO