The 4-loop β-function in the 2D non-Abelian Thirring model, and comparison with its conjectured “exact” form