In this paper we develop and compare different real-time methods to calculate spectral functions. These include classical-statistical simulations, the Gaussian state approximation (GSA), functional renormalization group (FRG) formulated on Keldysh closed-time path. Our test-bed system is quartic anharmonic oscillator, a single self-interacting bosonic degree of freedom, coupled an external heat...