ABSTRACTWe use Bayesian analysis in order to constrain the equation of state for nuclear matter from astrophysical data related to the recent measurements from the NICER mission, LIGO/Virgo collaboration, and probability distributions of mass and radius from other 12 sources, including thermonuclear busters, and quiescent low-mass X-ray binaries. For this purpose, we base our study on a relativistic hadronic mean field model including an ω − ρ interaction. Our results indicate optimal ranges for some bulk parameters at the saturation density, namely, effective mass, incompressibility, and symmetry energy slope (L0). For instance, we find $L_0 = 50.79^{+15.16}_{-9.24}$ MeV (Case 1) and $L_0 = 75.06^{+8.43}_{-4.43}$ MeV (Case 2) in a 68 per cent confidence interval for the two cases analysed (different input ranges for L0 related to the PREX-II data). The respective parametrizations are in agreement with important nuclear matter constraints, as well as observational neutron star data, such as the dimensionless tidal deformability of the GW170817 event. From the mass–radius curves obtained from these best parametrizations, we also find the ranges of 11.97 km ≤ R1.4 ≤ 12.73 km (Case 1) and 12.34 km ≤ R1.4 ≤ 13.06 km (Case 2) for the radius of the $1.4\, \mathrm{M}_\odot$ neutron star.