Using a relativistic effective Lagrangian at the hadronic level, near-threshold ω and φ meson productions in proton proton (pp) collisions, pp → ppω/φ, are studied within the distorted wave Born approximation. Both initial and final state pp interactions are included. In addition to total cross section data, both ω and φ angular distribution data are used to constrain further the model parameters. For the pp → ppω reaction we consider two different possibilities: with and without the inclusion of nucleon resonances. The nucleon resonances are included in a way to be consistent with the π − p → ωn reaction. It is shown that the inclusion of nucleon resonances can describe the data better overall than without their inclusion. However, the SATURNE data in the range of excess energies Q < 31 MeV are still underestimated by about a factor of two. As for the pp → ppφ reaction it is found that the presently limited available data from DISTO can be reproduced by four sets of values for the vector and tensor φN N coupling constants. Further measurements of the energy dependence of the total cross section near threshold energies should help to constrain better the φN N coupling constant. PACS number(s): 25.40. Ve, 24.10.Jv, 13.75.Cs