In this work we parameterize the equation of state of dense neutron star (NS) matter with four pressure parameters of {p 1 ,p 2 ,p 3 ,p 4 } and then set the combined constraints with the data of GW 170817 and the data of six low-mass X-ray binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range ofp 1 , the gravitational-wave data plays the leading role in boundingp 2 , and the LMXB data as well as the lower bound on the maximal gravitational mass of non-rotating NSs govern the constraints onp 3 andp 4 . Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability (Λ), moment of inertia (I), and binding energy (BE) of NSs. For a canonical mass of 1.4M , we have I 1.4 = 1.43 +0.30 −0.13 × 10 38 kg · m 2 , Λ 1.4 = 390 +280 −210 , R 1.4 = 11.8 +1.2 −0.7 km and BE 1.4 = 0.16 +0.01 −0.02 M if the constraints from the nuclear data and the gravitational-wave data have been jointly applied. For the joint analysis of gravitational-wave data and the LMXB data, we have I 1.4 = 1.28 +0.15 −0.08 × 10 38 kg · m 2 , Λ 1.4 = 220 +90 −90 , R 1.4 = 11.1 +0.7 −0.6 km, and BE 1.4 = 0.18 +0.01 −0.01 M . These results suggest that the current constraints on Λ and R still suffer from significant systematic uncertainties, while I 1.4 and BE 1.4 are better constrained. arXiv:1909.06944v2 [astro-ph.HE]