dat4/mcembirds.R

devtools::install_github("franciscorichter/dmea")
library('dmea')
tree = read.nexus(file='BirdTree.tre')
tree = phylo2p(phylo)
time = proc.time()
EMbirds = EM_phylo(wt=tree$wt,init_par=c(0.2,0.02,34000),parallel = TRUE,n_trees = 10,n_it = 30)
print(proc.time()-time)
save(EMbirds, file = 'EMbirsd.RData')
franciscorichter/dmea documentation built on May 16, 2019, 1:54 p.m.