sigma_squared_star function

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