  We compare analytic predictions for real and Fourier space two-point
statistics for biased tracers from a variety of Lagrangian Perturbation Theory
approaches against those from state of the art N-body simulations in $f(R)$
Hu-Sawicki and the nDGP braneworld modified gravity theories.
  We show that the novel physics of gravitational collapse in scalar tensor
theories with the chameleon or the Vainshtein screening mechanism can be
effectively factored in with bias parameters analytically predicted using the
Peak-Background Split formalism when updated to include the environmental
sensitivity of modified gravity theories as well as changes to the halo mass
function.
  We demonstrate that Convolution Lagrangian Perturbation Theory (CLPT) and
Standard Perturbation Theory (SPT) approaches provide accurate analytic methods
to predict the correlation function and power spectra, respectively, for biased
tracers in modified gravity models and are able to characterize both the BAO,
power-law and small scale regimes needed for upcoming galaxy surveys such as
DESI, Euclid, LSST and WFIRST.
