compute_p_si function

Compute choice probabilities at posterior samples