liks: Calculate the log-likelihoods of all possible models

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

Calculate the log-likelihoods of all possible models

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)

CenterForStatistics-UGent/RCM documentation built on April 24, 2023, 8:26 p.m.