Four-loop―order β function for two-dimensional nonlinear sigma models