Function_ConductSimExperimentOU_ModelTestBM_OU: Function_ConductSimExperimentOU_ModelTestBM_OU: function to...

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

Description

This function returns a matrix.SimulationResults for the simulation results

Usage

Function_ConductSimExperimentOU_ModelTestBM_OU(
  list.Model_Simulation_OU,
  numeric.NumberOfReps
)

Arguments

list.Model_Simulation_OU

List containing components of simulation model

Value

vector.Distances Vector containing the distances computed between the two focal tree models

Examples

################
# Load depends #
################
#library(ape)
#library(phytools)
#library(geiger)
#library(vioplot)
#library(gaussDiff)
#library(MASS)

##########################
# Build simulation model #
##########################
phylogeny.RandomTree <- pbtree(b = 1, d = 0, n = 50)
vector.Model_Simulation_Theta <- c(1, 10)
names(vector.Model_Simulation_Theta) <- c("Sig2", "alpha")

list.Model_Simulation_OU <- list(handle.Phylogeny = phylogeny.RandomTree,
                                 string.Model = "OU",
                                 vector.Z = rep(0, length(phylogeny.RandomTree$tip.label)),
                                 vector.Theta = vector.Model_Simulation_Theta)


handle.Sim_Results <- Function_ConductSimExperimentOU_ModelTestBM_OU(list.Model_Simulation_OU = list.Model_Simulation_OU, numeric.NumberOfReps = 5)

###################
# Compare this to #
###################
##########################
# Build simulation model #
##########################
phylogeny.RandomTree <- pbtree(b = 1, d = 0, n = 100)
vector.Model_Simulation_Theta <- c(1, 0.0001)
names(vector.Model_Simulation_Theta) <- c("Sig2", "alpha")

list.Model_Simulation_OU <- list(handle.Phylogeny = phylogeny.RandomTree,
                                 string.Model = "OU",
                                 vector.Z = rep(0, length(phylogeny.RandomTree$tip.label)),
                                 vector.Theta = vector.Model_Simulation_Theta)


handle.Sim_Results <- Function_ConductSimExperimentOU_ModelTestBM_OU(list.Model_Simulation_OU = list.Model_Simulation_OU, numeric.NumberOfReps = 5)

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