Context. Recent developments in multi-messenger astronomy bring important information on matter far beyond the one reachable in terrestrial laboratories, as matter existing inside neutron stars (NS). As a consequence, there is a large interest on the problem of extracting the equation of state of stellar matter as well as of determining the internal composition of NS from the knowledge of observational data generally using statistical methods. Aims. Determine possible signatures of the presence of hyperons inside neutron stars, using a set of EOS that satisfy laboratory, theoretical and observational constraints. Methods. Three sets of models of 15000 equations of state (EOS) based on a density dependent relativistic mean field description of hadronic matter are built. One of the sets is restricted to nucleonic matter, the second set includes Λ-hyperons and nucleons and the third set includes Λ and Ξ − -hyperons and nucleons. Electrical neutrality and β-equilibrium are imposed in the three sets. Nuclear matter properties, hypernuclei properties and observational information are used to constraint the models within a Bayesian inference approach.Results. General properties of neutron stars such as the maximum mass, radius, tidal deformability, proton fraction, hyperon fraction and speed of sound are discussed. It is shown that the two solar mass constraint imposes that neutron stars described by equations of state that include hyperons have in average a larger radius, 0.5km, and a larger tidal deformability, 150, than the stars determined from a nucleonic equation of state, while the speed of sound at the center of the star is more than 25% smaller. Conclusions. If a 1.4 M star with a radius 12.5 km is measured it is quite improbable that a massive star described by the same model contains hyperons. A similar conclusion is drawn if a two solar mass star with a radius 11.5 km or a neutron star with a mass above 2.2 M is observed: the possible hyperon content of these stars is ruled out or very reduced. The hyperon presence inside neutron stars is compatible with the present NICER mass-radius observations of the pulsars PSR J0030+0451 and PSR J0740+6620 and the gravitational wave detection GW170817.