Schuster_compute_pi function

Compute matrix of model-based proportions pi.