tau_star function

Calculates parameter $\tau^*_i$ of $q^*(\gamma_i)$ for $i=1,...,K$ clusters to optimize the evidence based lower bound (ELBO) in survregVB.frailty.fit.