The fully relativistic multiconfiguration Dirac-Hartree-Fock method is used to compute excitation energies and lifetimes for the 143 lowest states of the 3s^2^3p^3^, 3s3p^4^, 3s^2^3p^2^3d, 3s3p^3^3d, 3p^5^, 3s^2^3p3d^2^ configurations in P-like ions from Cr X to Zn XVI. Multipole (E1, M1, E2, M2) transition rates, line strengths, oscillator strengths, and branching fractions among these states are also given. Valence-valence and core-valence electron correlation effects are systematically accounted for using large basis function expansions. Computed excitation energies are compared with the NIST ASD and CHIANTI compiled values and previous calculations. The mean average absolute difference, removing obvious outliers, between computed and observed energies for the 41 lowest identified levels in Fe XII, is only 0.057%, implying that the computed energies are accurate enough to aid identification of new emission lines from the Sun and other astrophysical sources. The amount of energy and transition data of high accuracy are significantly increased for several P-like ions of astrophysics interest, where experimental data are still very scarce.