Agresti_compute_pi function

Computes the matrix pi of model-based proportions