This paper aims to investigate the stationary probability density functions of the Duffing oscillator with time delay subjected to combined harmonic and white noise excitation by the method of stochastic averaging and equivalent linearization. By the transformation based on the fundamental matrix of the degenerate Duffing system, the paper shows that the displacement and the velocity with time delay in the Duffing oscillator can be computed approximately in non-time delay terms. Hence, the stochastic system with time delay is transformed into the corresponding stochastic non-time delay equation in Ito sense. The approximate stationary probability density function of the original system can be found by combining the stochastic averaging method, the equivalent linearization method, and the technique of auxiliary function. The response of Duffing oscillator is investigated. The analytical results are verified by numerical simulation results.