init.compute_betas_from_shifts function

Initialization for the computation of the optimal values