themetagenomics (version 1.0.2)

extract: Extract summary statistics

Description

Extract summary statistics

Usage

extract(object, ...)
# S3 method for effects
extract(object, ...)

Arguments

object

Object of class effects, fit via hmc.

...

Additional arguments for methods.

Value

A list containing

summary

Rstan summary of parameters from model.

flagged

Vector of parameter names with Rhat > 1.1.

Methods (by class)

• effects: Extract summary statistics from HMC effects object

  Extracts the summary information in a form conducive with vis methods, specifically in cases when return_summary was set to FALSE.

Examples

# NOT RUN {
formula <- ~DIAGNOSIS
refs <- 'Not IBD'

dat <- prepare_data(otu_table=GEVERS$OTU,rows_are_taxa=FALSE,tax_table=GEVERS$TAX,
                    metadata=GEVERS$META,formula=formula,refs=refs,
                    cn_normalize=TRUE,drop=TRUE)

# }
# NOT RUN {
topics <- find_topics(dat,K=15)

functions <- predict(topics,reference_path='/references/ko_13_5_precalculated.tab.gz')

function_effects <- est(functions,level=3,
                        iters=500,method='hmc',
                        prior=c('laplace','t','laplace'),
                        return_summary=FALSE)

function_effects_summary <- extract(function_effects)
# }
# NOT RUN {
# }