McCullagh_compute_pi_from_beta function

Computes matrix of p-values pi based on x and current value of beta.