compute_ph_p_pds_all function

Compute all the ph_p_pds summary scores