The gravitational wave provides a new method to examine General Relativity and its alternatives in the high speed, strong field regime. Alternative theories of gravity generally predict more polarizations than General Relativity, so it is important to study the polarization contents of theories of gravity to reveal the nature of gravity. In this talk, we analyze the polarization contents of Horndeski theory and f(R) gravity. We find out that in addition to the familiar plus and cross polarizations, a massless Horndeski theory predicts an extra transverse polarization, and there is a mix of pure longitudinal and transverse breathing polarizations in the massive Horndeski theory and f(R) gravity. It is possible to use pulsar timing arrays to detect the extra polarizations in these theories. We also point out that the classification of polarizations using Newman–Penrose variables cannot be applied to massive modes. It cannot be used to classify polarizations in Einstein-æther theory or generalized Tensor-Vector-Scalar (TeVeS) theory, either.