We use the path integral approach to a two-dimensional noncommutative harmonic oscillator to derive the
partition function of the system at finite temperature. It is shown that the result based on the Lagrangian formulation
of the problem, coincides with the Hamiltonian derivation of the partition function.