sample_c_two {ddtlcm}

Sample divergence function parameter c for a(t) = c / (1-t)^2 through Gibbs sampler

Description

Sample divergence function parameter c for a(t) = c / (1-t)^2 through Gibbs sampler

Usage

sample_c_two(shape0, rate0, tree_structure)

Arguments

shape0

shape of the Inverse-Gamma prior

rate0

rate of the Inverse-Gamma prior

tree_structure

a data.frame containing the divergence times and number of data points to the left and right branches of internal nodes on the tree

Value

a numeric value of the newly sampled c