We constrain Einstein-æther gravity – a Lorentz-violating extension of General Relativity in which a dynamical, unit timelike vector field selects a preferred frame – using updated high-precision pulsar timing observations of PSR J1738+0333 from EPTA second Data Release and the NANOGrav 9year release, in combination with ToAs from Arecibo, Green Bank, Nançay, Parkes, and Westerbork. Our method accounts for both conservative and dissipative first post-Newtonian corrections arising from Lorentz violation; here we apply it to PSR J1738+0333 using the Bayesian timing pipeline Vela to process the full ToA dataset. We sample the joint posterior over binary component masses, post-Keplerian parameters and center-of-mass velocity components, and then apply a resampling scheme to propagate posteriors into robust constraints on the fundamental theory parameters, obtaining the most stringent strong-field bounds on the Einstein-æther coupling constants from a single binary pulsar system to date.