extract: Extract summary statistics

Description Usage Arguments Value Methods (by class) Examples

View source: R/extract_effects.R

Description

Extract summary statistics

Usage

1
2
3
4
extract(object, ...)

## S3 method for class '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)

Examples

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

## End(Not run)

EESI/themetagenomics documentation built on May 10, 2020, 1:40 a.m.