Summary: The Duffing oscillator under external non-Gaussian excitations is investigated by means of statistical linearization. The input process is modelled as a polynomial of a Gaussian process or as a renewal-driven impulse process. Four criteria of statistical linearization are considered. The interarrival times of the renewal process are distributed according to a Pearson type III law. Predictions of the stationary variance are compared with Monte Carlo simulations.