In this work, we present a preliminary study about constraints on saturation properties of nuclear matter using the Bayesian analysis method, where it was considered an equation of state (EoS) that allows $\omega-\rho$ interactions. Regarding the priors, we have used a uniform distribution on the values of saturation density, binding energy, effective mass, incompressibility, symmetry energy, and its slope at nuclear saturation density, whose ranges agree with the current literature. The posterior has been computed using the procedure of Bayesian inference, where observational data of neutron star masses and radii from gravitational wave event (GW170817) and the most recent published NICER events are employed. The results indicate lower values of $L_0$ for the EoS with $\omega-\rho$ interactions if compared to $\sigma^3-\sigma^4$ ones, $65.4^{+18.1}_{-11.3}$ MeV and $78.8^{+8.0}_{-1.3}$ MeV respectively.