This just writes new values into P as if it were an output variable. Internal function used within estimate_m_rd.
estimate_m_rd
gibbsP(p, X, pri)