When applying the probit transformation on the Dirichlet samples, this function is used as part of the stick-breaking algorithm. It computes the length of the remaining stick when the current element is broken off.
.computeLengthOfRemainingStick(theta_mat, k, hyp_direction)
matrix with samples from truncated Dirichlet density
current parameter index
specifies whether the imposed inequality constrained imposes an increasing (i.e., 'smaller') or decreasing (i.e., 'larger') trend
stick length