Function_ComputePairWiseModelDistances_BM: Function_ComputePairWiseModelDistances_BM: function to...

In radamsRHA/PRDATR: Computing probabilistic distances between phylogenetic models of continuous trait evolution

Description

This function returns a list of pairwise model distances

Usage

Function_ComputePairWiseModelDistances_BM(
  phylogeny.TreeSet,
  matrix.ModelParams
)

Arguments

phylogeny.TreeSet	Set of phylogenetic trees
matrix.ModelParams	Matrix of model parameters (BM) for each tree in the given set

Value

List List of pairwise model distances

Examples

################
# Load depends #
################
library(ape)
library(geiger)
library(gaussDiff)

#########################
# Simulate random trees #
#########################
handle.SimulatedRandomTrees <- rmtree(N = 10, n = 10)
vector.SimulatedDataset <- fastBM(tree = handle.SimulatedRandomTrees[[1]], a = 0, mu = 1)

#######################
# Fit models to trees #
#######################
handle.FittedPhylogeneticModels <- Function_FitCont_TreeSet_BM(phylogeny.TreeSet = handle.SimulatedRandomTrees, vector.InputData = vector.SimulatedDataset)

##############################
# Compute pairwise distances #
##############################
handle.Results <- Function_ComputePairWiseModelDistances_BM(phylogeny.TreeSet = handle.SimulatedRandomTrees, matrix.ModelParams = handle.FittedPhylogeneticModels)

radamsRHA/PRDATR documentation built on Dec. 25, 2021, 11:44 a.m.