Reduced χ² for the fit of the redshift space halo power spectrum multipoles, computed with the TNS and the EFTofLSS prescriptions, to the measured multipoles from the DGP and f(R) ELEPHANT simulations, as a function of k_max and at redshift z = 1. From top to bottom, TNS fit to DGP, EFTofLSS fit to DGP, TNS fit to f(R) with |f_R0| = 10⁻⁵, TNS fit to f(R) with |f_R0| = 10⁻⁶, and EFTofLSS fit to f(R). Crosses, circles and pentagons indicate the exact computation, the EdS and USA approximation and the ΛCDM-scr approximation, respectively. The ΛCDM cases are shown as black (TNS) and grey (EFTofLSS) squares. The error bars are the 1σ errors on the χ² statistic with N_dof = 2N_k − N_x degrees of freedom, where N_k are the number of wave modes used in calculating χ². We use N_x = 4 for the TNS model using Eulerian bias and N_x = 5 for EFTofLSS and the TNS model with Q-bias. All EFTofLSS exact and f(R) calculations are performed using PyBird, while the DGP EdS and USA calculations are performed using PBJ. The ΛCDM predictions are made using the exact kernel calculations.