compute_reliabilities function

Compute reliabilities