randomseqsig: Compositional Significance Testing

In escamero/mirlyn: Multiple Iterations of Rarefying for Library Normalization

Compositional Significance Testing

Description

The randomseqsig() function will identify whether a taxonomic group of interest is significantly over-represented or underrepresented in the community.

Usage

randomseqsig(x, taxlevel = "Phylum", group = "Cyanobacteria", nshuff = 1000)

Arguments

x	is a phyloseq object.
taxlevel	is the taxonomic level of the group of interest.
group	is the name of the taxonomic group of interest.
nshuff	is the number of times data shuffling is performed.

Value

A list object which contatins: 1) objshuff: The shuffled sequence variant counts. 2) stats: The mean, standard deviation, real counts and p-values of the shuffled data.

Examples

library(mirlyn)
data(example)
significance_example <- randomseqsig(example, nshuff = 10)

escamero/mirlyn documentation built on Aug. 6, 2023, 9:30 p.m.