ICgenlasso: Calculates the Bayesian information criterion or Akaike's...

In liudoubletian/phyloMDA: An R package for phylogeny-aware microbiome data analysis

Calculates the Bayesian information criterion or Akaike's information criterion for a genlasso object.

Description

This function extracts BIC or AIC from a genlasso object.

Usage

ICgenlasso(fit, IC = "BIC")

Arguments

fit	a genlasso object.
IC	IC = BIC for Bayesian information criterion and IC = AIC for Akaike's information criterion.

Value

Returns the minimal AIC or BIC and the corresponding lambda.

Examples

library(phyloseq)
otu.tab <- t(otu_table(combo.phyloseq.obj))
X <- eBay_comps(otu.tab, prior = "Dirichlet")
tree <- phy_tree(combo.phyloseq.obj)
metadata <- sample_data(combo.phyloseq.obj)
y <- metadata$bmi
fit <- TASSO(y, X, tree)
IC <- ICgenlasso(fit, IC="AIC")
IC <- ICgenlasso(fit, IC="BIC")

liudoubletian/phyloMDA documentation built on April 28, 2022, 6:51 p.m.