liks: Calculate the log-likelihoods of all possible models

In RCM: Fit row-column association models with the negative binomial distribution for the microbiome

Description

Calculate the log-likelihoods of all possible models

Usage

liks(rcm, Sum = TRUE)

Arguments

rcm	an object of the RCM class
Sum	a boolean, should log-likelihoods be summed?

Value

If Sum is FALSE, a named array log-likelihoods of the independence model and all models with dimension 1 to k, including after filtering on confounders. Otherwise a table with log-likelihoods, deviance explained and cumulative deviance explained.

Examples

data(Zeller)
require(phyloseq)
tmpPhy = prune_taxa(taxa_names(Zeller)[1:100],
prune_samples(sample_names(Zeller)[1:50], Zeller))
zellerRCM = RCM(tmpPhy, round = TRUE)
liks(zellerRCM)

RCM documentation built on Nov. 8, 2020, 5:22 p.m.