simStudy: Simulate OTU table for presence-absence or abundance-weighted...

In smvargas/PERMANOVA: Power estimation for microbiome studies using pairwise distances and PERMANOVA.

Description

Incorporates within-group distance simulation from OTU number and depth of subsampling (rarefying) with group-level differences modeled by segregating OTU membership between groups.

Usage

simStudy(group_size_vector = c(100, 100, 100), otu_number = 1000,
  sequence_depth = 1, rare_depth = 0.4, effect = 0.2)

Arguments

group_size_vector	numeric vector representing subjects per exposure/intervention group
otu_number	number of simulated OTUs
sequence_depth	number of sequence counts per OTU bin
rare_depth	proportion of sequence counts to retain after subsampling
effect	proportion of unique community membership in affected group of subjects

Details

For weighted analysis, specify a sequence depth greater than 1; for unweighted analyses, sequence depth is 1 (default).

Value

two-dimensional matrix OTU table, with row and column names to suit downstream analysis

See Also

calcUJstudy, calcWJstudy

Examples

simStudy(c(16,16,16),100,10,0.8,0.1)

smvargas/PERMANOVA documentation built on May 20, 2019, 12:37 p.m.