div_time: Sample divergence time on an edge uv previously traversed by...
In ddtlcm: Latent Class Analysis with Dirichlet Diffusion Tree Process Prior

div_time

R Documentation

Sample divergence time on an edge uv previously traversed by m(v) data points

Sample divergence time on an edge uv previously traversed by m(v) data points

div_time(t_u, m_v, c, c_order = 1, alpha = 0, theta = 0)

`t_u`	a number in the interval (0, 1) indicating the divergence time at node u
`m_v`	an integer for the number of data points traversed through node v
`c`	a positive number for the divergence hyperparameter. A larger value implies earlier divergence on the tree
`c_order`	equals 1 if using divergence function `a(t) = c / (1-t)`, or 2 if a(t) = c / (1-t)^2. Default is 1
`alpha`, `theta`	hyparameter of branching probability a(t) Gamma(m-alpha) / Gamma(m+1+theta) For DDT, alpha = theta = 0. For general multifurcating tree from a Pitman-Yor process, specify positive values to alpha and theta. It is, however, recommended using alpha = theta = 0 in inference because multifurcating trees have not been tested rigorously.