With recent advances in astronomical observations, major progress has been made in determining the pressure of neutron star matter at high density. This pressure is constrained by the neutron star deformability, determined from gravitational waves emitted in a neutron-star merger, and the mass-radii relation of two neutron stars, determined from a new X-ray observatory on the International Space Station. Previous studies have relied on nuclear theory calculations to constrain the equation of state at low density. Here we use a combination of constraints composed of three astronomical observations and twelve nuclear experimental constraints that extend over a wide range of densities. A Bayesian inference framework is then used to obtain a comprehensive nuclear equation of state. This data-centric result provides benchmarks for theoretical calculations and modeling of nuclear matter and neutron stars. Furthermore, it provides insights into the microscopic degrees of freedom of the nuclear matter equation of state and on the composition of neutron stars and their cooling via neutrino radiation.