Abstract : The variability of a set of biomedical signals is analysed following some criterions. Firstly, when the signal source is unique, the variability may be, for one sensor in function of the time occurrence of the signals, for several sensors at the same time, a spatial variability. The case of several sources corresponds to the inter subject variability. Then the distinction between amplitude variability, time variability and shape variability is introduced. Functional Data Analysis is viewed as modelling variability through time warping functions. These functions are either composed with the signals, leading to a lot of curve registration techniques, or composed with the normalised integral functions of the signals when assumed to be positive. The later case is well suited to shape analysis and shape variability estimation around an averaged shape. The Integral Shape Averaging (ISA) and Corrected Integral Shape Averaging (CISA) algorithms are recalled. Two applications are then presented, including shape clustering. The first one concerns the correlation between the occurrence of an obstructive sleep apnoea and a shape change of the ECG P-waves. The second one is related to the Ensemble Spontaneous Activity (ESA) recorded near the cochlea. It is shown, by simulation, that the shape of the histogram of the ESA amplitude is very sensitive to a localised correlated firing of the fibres, possibly associated to Tinnitus, on the contrary of the PSD amplitude at 1 kHz.