Constraining the Neutron Star Equation of State using Observational Data