The power spectrum and bispectrum of the Sunyaev-Zel’dovich effect have the potential to constrain modified gravity theories. In this work, we do a combined Fisher forecast usind these to observables, to check how well we can constrain f(R) gravity models in an ideal, cosmic variance limited experiment. We find that the constraining power of the bispectrum, in particular, is significant and merits further research, not only for checks of modified gravity, but also constraints on ICM physics.