Multiscale analysis of human heartbeat dynamics has been proved effective in characterizeing cardiovascular control physiology in health and disease. However, estimation of multiscale properties can be affected by the interpolation procedure used to preprocess the unevenly sampled R-R intervals derived from the ECG. To this extent, in this study we propose the estimation of wavelet coefficients and wavelet leaders on the output of inhomogeneous point process models of heartbeat dynamics. The RR interval series is modeled using probability density functions (pdfs) characterizing and predicting the time until the next heartbeat event occurs, as a linear function of the past history. Multiscale analysis is then applied to the pdfs' instantaneo...