unifrac<-function(data, phylo) {
util.pd <- beta.pd.utils(data,phylo)
betadiv <- betajac.pd(util.pd)
results=NULL
results <- dist.mat(data,betadiv[,1])
names(results)<- colnames(betadiv)[1]
results
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.