# \donttest{
# run bGMYC
bgmyc_res <- bGMYC::bgmyc.singlephy(ape::as.phylo(geophagus_beast),
mcmc = 11000,
burnin = 1000,
thinning = 100,
t1 = 2,
t2 = ape::Ntip(geophagus_beast),
start = c(1, 0.5, 50)
)
# create a tibble
bgmyc_df <- bgmyc_tbl(bgmyc_res, ppcutoff = 0.05)
# check
bgmyc_df
# }
Run the code above in your browser using DataLab