The equation of state provided by effective models of strongly interacting matter should comply with the restrictions imposed by current astrophysical observations on compact stars. We explore the possibility, that there is quark matter in the center of neutron stars. For the equation of state of the neutron star matter, we use a relativistic mean field model for the hadronic one and the (axial‐)vector meson extended linear sigma model for the quark phase. We assume that the transition is a cross‐over between these phases. We connect the two branches in this cross‐over regime continuously by a polynomial. This equation of state should comply with the restrictions imposed by current astrophysical observations on compact stars. Our model has four free parameters that should be fitted to experimental data. By using a Bayesian analysis, we found that our model could explain all the experimental constraints; furthermore, we have obtained the most probable parameter set for our model.