Within the framework of thermo field dynamics, we derive the partition functions for a harmonic oscillator obeying parafermi and parabose statistics. The thermodynamic functions are studied for different limits of the order of the statistics and are shown to have the expected behaviour.