We present measurements of the spatial clustering statistics in redshift space of various scalar field modified gravity simulations. We utilise the two-point and the three-point correlation functions to quantify the spatial distribution of dark matter halos within these simulations and thus discern between the models. We compare $Lambda$CDM simulations to various modified gravity scenarios and find consistency with previous work in terms of 2-point statistics in real and redshift-space. However using higher order statistics such as the three-point correlation function in redshift space we find significant deviations from $Lambda$CDM hinting that higher order statistics may prove to be a useful tool in the hunt for deviations from General Relativity.