Abstract: Neutron star physics has wrestled with the longstanding challenge of the hyperon puzzle, attempting to reconcile lowered theoretical predictions of maximum masses due to hyperons with astrophysical observations based on the measured masses of the heaviest pulsars. In this study, we conduct a comprehensive statistical analysis of equations of state (EoS) for neutron stars with hyperons, including both laboratory data and astrophysical observations. Results from the statistical analysis reveal the important role of the correlations between the scalar and vector channels of hyperon-nucleon interactions deduced from available \varLambda-separation-energy data of single \varLambda hypernuclei. The analysis preliminarily quantifies uncertainties in hyperon star properties due to the uncertain hyperon-nucleon interaction in dense matter, and the maximum mass of hyperon star is found to reach ~2.2\,M_\odot, challenging the existence of the hyperon puzzle. As part of a broader initiative connecting nuclear physics and astronomy to quantitatively determine neutron star EoS, the study provides valuable insights into the hyperon puzzle and its implications for our understanding of neutron star interiors. Moreover, the investigation addresses the lack of precise knowledge regarding hyperonic interactions, emphasizing the need for additional hypernuclear data through a combined effort involving theory, experiments, and observations.