  Bayesian methods are used to constrain the density dependence of the QCD
Equation of State (EoS) for dense nuclear matter using the data of mean
transverse kinetic energy and elliptic flow of protons from heavy ion
collisions (HIC), in the beam energy range $\sqrt{s_{\mathrm{NN}}}=2-10 GeV$.
The analysis yields tight constraints on the density dependent EoS up to 4
times the nuclear saturation density. The extracted EoS yields good agreement
with other observables measured in HIC experiments and constraints from
astrophysical observations both of which were not used in the inference. The
sensitivity of inference to the choice of observables is also discussed.
