We studied the thermal effects on the behavior of incommensurate solid 4He at low temperatures using the path integral Monte Carlo method. Below a certain temperature, depending on the density and the structure of the crystal, the vacancies delocalize and a finite condensate fraction appears. We calculated the vacancy formation energy as a function of the temperature and observed a behavior compatible with a two-step structure, with a gap of few K appearing at the onset temperature of off-diagonal long-range order. Estimation of the energy cost of creating two vacancies seems to indicate an effective attractive interaction among the vacancies but the large error inherent to its numerical estimation precludes a definitive statement.