The free-energy of discrete nonlinear Klein-Gordon (NKG) systems with anharmonic interparticle interactions is derived by means of the transfer integral operator method, with the first lattice corrections and kink-kink interactions taken into account. Two particular substrate potentials are considered: the φ − f our and the sine-Gordon (sG). We show that, in the general case where the system exhibits the kink soliton like excitations, the correction factors, due to the lattice discreteness, appearing in the free-energy and in the lattice corrected static kink soliton energy, depend on the temperature through a coupling of the interparticle anharmonicity strength to the temperature. Similarly, in the purely anharmonic NKG systems, characterized by the absence of the linear dispersion, where thermodynamic properties are sensitive to kink compactons, we find also that the correction factors are temperature dependent. In both cases, they decrease with increasing temperatures, although the correction factors verify different temperature laws.