The temperature dependence of the methane 12CH4 partition function was investigated up to T=3000 K. The ro-vibrational energy levels were calculated using high-order contact transformation (CT) method that permitted building accurate effective Hamiltonian models for highly excited vibration–rotation states from the molecular potential energy surface. Comparisons with variational calculations using large basis sets as well as convergence and extrapolation issues are discussed. The accuracy of the present calculations for the methane partition function was estimated as 1–2% at 2000 K. This estimated error is considerably smaller than the difference with other calculations available in the literature.