Nothing
#' Evolve two traits up a tree
#'
#' Given a phylogeny, will generate associated trait data for two traits following a
#' Brownian motion evolution model.
#'
#' @param tree A phylogeny.
#' @param sigma The Brownian motion rate parameter.
#'
#' @details Evolves two traits independently up phylogeny with Brownian motion evolution
#' process. Sigma is currently set to 0.1 as default.
#'
#' @return A list where the first object is a phylogeny with the desired number of species
#' and the second object is a matrix of trait values for those species.
#'
#' @export
#'
#' @importFrom ape rTraitCont
#'
#' @references Miller, E. T., D. R. Farine, and C. H. Trisos. 2016. Phylogenetic community
#' structure metrics and null models: a review with new methods and software.
#' Ecography DOI: 10.1111/ecog.02070
#'
#' @examples
#' tree <- geiger::sim.bdtree(b=0.1, d=0, stop="taxa", n=50)
#'
#' results <- evolveTraits(tree)
evolveTraits <- function(tree, sigma=0.1)
{
trait1 <- ape::rTraitCont(tree, model="BM", sigma=sigma)
trait2 <- ape::rTraitCont(tree, model="BM", sigma=sigma)
traits <- cbind(trait1, trait2)
output <- list(tree, traits)
return(output)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.