library(phangorn) library(devtools) library(NELSI) library(TreeSim) library(phytools) library(EnvStats) library(phylobase) library(ggplot2) library(Rmisc) n_tips <- 8 lambda <- 1 mu <- 0 n_reps <- 1 ########## tip_vector <- seq(1, n_tips, 1) node_numbers <- seq(length(tip_vector)+1, (length(tip_vector)+2) + (length(tip_vector)-3), 1) ########## entire_tree <- sim.bd.taxa(n_tips, n_reps, lambda, mu, frac = 1, complete = TRUE, stochsampling = FALSE) entire_tree <- entire_tree[[1]] ########## source("simulate_sequences.R")