compute_ph_y_mctq_all function

Compute all the MCTQ variables