We compare Wiener chaos and stochastic collocation methods for linear advection-reaction-diffusion equations with multiplicative white noise. Both methods are constructed based on a recursive multi-stage algorithm for long-time integration. We derive error estimates for both methods and compare their numerical performance. Numerical results confirm that the recursive multi-stage stochastic collocation method is of order $Δ$ (time step size) in the second-order moments while the recursive multi-stage Wiener chaos method is of order $Δ^{\mathsf{N}}+Δ^2$ ($\mathsf{N}$ is the order of Wiener chaos) for advection-diffusion-reaction equations with commutative noises, in agreement with the theoretical error estimates. However, for non-commutative noises, both methods are of order one in the second-order moments.