The major aim of our study is to demonstrate that a concerted combination of time-variant, frequency-selective, linear and nonlinear analysis approaches can be beneficially used for the analysis of heart rate variability (HRV) in epileptic patients to reveal premonitory information regarding an imminent seizure and to provide more information on the mechanisms leading to changes of the autonomic nervous system. The quest is to demonstrate that the combined approach gains new insights into specific short-term patterns in HRV during preictal, ictal, and postictal periods in epileptic children. The continuous Morlet-wavelet transform was used to explore the time-frequency characteristics of the HRV using spectrogram, phase-locking, band-power and quadratic phase coupling analyses. These results are completed by time-variant characteristics derived from a signal-adaptive approach. Advanced empirical mode decomposition was utilized to separate out certain HRV components, in particular blood-pressure-related Mayer waves (≈0.1 Hz) and respiratory sinus arrhythmia (≈0.3 Hz). Their time-variant nonlinear predictability was analyzed using local estimations of the largest Lyapunov exponent (point prediction error). Approximately 80-100 s before the seizure onset timing and coordination of both HRV components can be observed. A higher degree of synchronization is found and with it a higher predictability of the HRV. All investigated linear and nonlinear analyses contribute with a specific importance to these results.