generateTree {HMPTrees}

Generate Test Trees

Description

This function will take several initial trees and will randomly populate new trees based on the originals.

Usage

generateTree(data, numReadsPerSamp, theta = NULL, level = "genus", split = ".")

Arguments

data

A data frame in which each column contains the rdp read counts for every taxa given in the row names.

numReadsPerSamp

A vector specifying the number of reads or sequence depth for each sample.

theta

When theta is not NULL the base tree is generated by using the dirmult function.

level

The depth the tree will go down to (kingdom, phylum, class, order, family, genus, species, subspecies).

split

This is the character that separates the taxa levels in the row names.

Value

A data frame containing the generated tree(s).

Author(s)

Patricio S. La Rosa, Elena Deych, Berkley Shands, William D. Shannon

Examples

	data(saliva)
	
	### Generate a the number of reads per sample
	### The first number is the number of reads and the second is the number of subjects
	nrs <- rep(10000, 2)
	
	gendata <- generateTree(saliva, nrs)