A weak values approach for testing simultaneous Einstein–Podolsky–Rosen elements of reality for non-commuting observables