simulateSignatures: Simulate a list of signatures based on a universe of taxa...

View source: R/montecarlo.R

simulateSignaturesR Documentation

Simulate a list of signatures based on a universe of taxa equal in size and length to some smaller set of signatures

Description

Simulate a list of signatures based on a universe of taxa equal in size and length to some smaller set of signatures

Usage

simulateSignatures(relevant.sigs, siglengths)

Arguments

relevant.sigs

A list of signatures providing a relevant background for drawing taxa

siglengths

An integer vector, the length of which provides the number of signatures to be simulated, and the integers of which provide the number of taxa to be simulated in each signature

Value

A list of signatures of the same number and individual lengths as found in my.dat

Examples

full.dat <- bugsigdbr::importBugSigDB()
my.dat <- full.dat[full.dat$Curator == "Mst Afroza Parvin", ]
relevant.dat <- full.dat[full.dat$`Body site` %in% my.dat$`Body site`, ]
relevant.sigs <- bugsigdbr::getSignatures(relevant.dat, tax.level = "genus")
siglengths <- sapply(my.dat, length)
simulateSignatures(relevant.sigs, siglengths)

waldronlab/bugSigSimple documentation built on Feb. 1, 2024, 12:22 a.m.