data(primates_ebd_log)
posterior <- read.RevBayes(primates_ebd_log, max_t = 65, n_samples = 10)
samples <- sample.congruence.class.posterior(posterior,
num.samples = 5,
rate.type = "extinction",
rate0.median = 0.1,
model = "MRF",
max.rate = 1.0)
print(samples)
Run the code above in your browser using DataLab