metaFunction: A statistical tool in profiling functions in a microbial community

The overall workflow of metaFunction is shown as below:

 

metaFunction package can be downloaded here, and the associated manual file.

The simulated sequence data (.fasta) can be found here.


R code for analyzing a built-in example dataset:

library(metaFunction)

require(data.table)

require(mefa)

require(lattice)

 

data(ExampleData)

data(refseq)

data(subs)

 

pre=preprocess(ExampleData)

 

epsilon=0.01

est=Twosteps(pre, epsilon)

 

bootsize=100   

boot.data=bootstrap(pre, bootsize, epsilon, method="B")

plot.abundance(boot.data)