Description Usage Arguments Details Value See Also Examples

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

1 2 3 4 |

`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 |

`nlambda` |
the number of |

`lambda.min.ratio` |
Smallest value for |

`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, |

`tune` |
tuning approach, available methods including AIC, BIC, GACV, BGACV. |

The `mvblps`

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

An object of classes `mvbfit`

and `lps`

, for which some methods are
available.

`mvbfit`

, `unifit`

, `stepfit`

, `mvb.simu`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ```
# 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)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.