Ventilatory thresholds (VT1 and VT2) are useful in many fields of medicine and sports. Nevertheless, their measurement is cumbersome and needs trained personnel. This work proposes an alternative method to predict VT1, VT2 and maximum loads in incremental maximal tests based on heart rate variability (HRV) analysis. Twelve competitive male cyclists executed an incremental exhaustive test. During the test, RR time series and gas concentrations were recorded. After artifact correction the power spectrum was estimated in a sliding window, and central frequency (CF) and bandwidth that contains half the total power (BW) were computed. An automatic algorithm recognized the loads where CF and BW undergo a significant change. These loads were used as inputs in linear regression models to predict VT1, VT2 and maximum loads. The errors of the predictions are similar to the load resolution.