Linear and nonlinear analysis of heart rate variability during propofol anesthesia for short-duration procedures in children

Using sensor-fusion and machine-learning algorithms to assess acute pain in non-verbal infants: a study protocol

INTRODUCTION

Objective pain assessment in non-verbal populations is clinically challenging due to their inability to express their pain via self-report. Repetitive exposures to acute or prolonged pain lead to clinical instability, with long-term behavioural and cognitive sequelae in newborn infants. Strong analgesics are also associated with medical complications, potential neurotoxicity and altered brain development. Pain scores performed by bedside nurses provide subjective, observer-dependent assessments rather than objective data for infant pain management; the required observations are labour intensive, difficult to perform by a nurse who is concurrently performing the procedure and increase the nursing workload. Multimodal pain assessment, using sensor-fusion and machine-learning algorithms, can provide a patient-centred, context-dependent, observer-independent and objective pain measure.

METHODS AND ANALYSIS

In newborns undergoing painful procedures, we use facial electromyography to record facial muscle activity-related infant pain, ECG to examine heart rate (HR) changes and HR variability, electrodermal activity (skin conductance) to measure catecholamine-induced palmar sweating, changes in oxygen saturations and skin perfusion, and electroencephalography using active electrodes to assess brain activity in real time. This multimodal approach has the potential to improve the accuracy of pain assessment in non-verbal infants and may even allow continuous pain monitoring at the bedside. The feasibility of this approach will be evaluated in an observational prospective study of clinically required painful procedures in 60 preterm and term newborns, and infants aged 6 months or less.

ETHICS AND DISSEMINATION

The Institutional Review Board of the Stanford University approved the protocol. Study findings will be published in peer-reviewed journals, presented at scientific meetings, taught via webinars, podcasts and video tutorials, and listed on academic/scientific websites. Future studies will validate and refine this approach using the minimum number of sensors required to assess neonatal/infant pain.

TRIAL REGISTRATION NUMBER

ClinicalTrials.gov Registry (NCT03330496).

Nociception monitoring tools using autonomic tone changes for intraoperative analgesic guidance in pediatric patients

Nociception monitoring devices using changes in autonomic nervous system activity have been developed in numerous ways. Although there have been few studies conducted on children, compared to the relatively higher number of studies on adults, most of the nociception monitors in children, as in adults, appear to be more useful than the standard clinical practice that uses hemodynamic parameters in the evaluation and treatment of intraoperative nociception (pain) during general anesthesia. Particularly, when monitoring the surgical pleth index (SPI) in anesthetized children, the application of a new target range of SPI values (≤ 40) to the SPI monitoring criteria seems to be necessary for providing a more proper intraoperative analgesia. The analgesia nociception index (ANI) shows promising results in anesthetized adults, and recently, positive results along with cardiorespiratory coherence have been reported in pediatric patients. Newborn infant parasympathetic evaluation (NIPE) could be useful for providing adequate analgesia in newborns, infants, and children under 2 years of age in anesthetized or awake states. In cases of skin conductance and pupillometry, further studies are needed. Understanding the pros, cons, and limitations of these nociception monitoring tools will provide more effective and safe intraoperative analgesia to pediatric patients undergoing general anesthesia, and it may also help to plan and conduct promising research on the use of perioperative nociception monitoring in pediatric patients in the future.

Multimodality monitoring: informatics, integration data display and analysis

The goal of multimodality neuromonitoring is to provide continuous, real-time assessment of brain physiology to prevent, detect, and attenuate secondary brain injury. Clinical informatics deals with biomedical data, information, and knowledge including their acquisition, storage, retrieval, and optimal use for clinical decision-making. An electronic literature search was conducted for English language articles describing the use of informatics in the intensive care unit setting from January 1990 to August 2013. A total of 64 studies were included in this review. Clinical informatics infrastructure should be adopted that enables a wide range of linear and nonlinear analytical methods be applied to patient data. Specific time epochs of clinical interest should be reviewable. Analysis strategies of monitor alarms may help address alarm fatigue. Ergonomic data display that present results from analyses with clinical information in a sensible uncomplicated manner improve clinical decision-making. Collecting and archiving the highest resolution physiologic and phenotypic data in a comprehensive open format data warehouse is a crucial first step toward information management and two-way translational research for multimodality monitoring. The infrastructure required is largely the same as that needed for telemedicine intensive care applications, which under the right circumstances improves care quality while reducing cost.

Computational Depth of Anesthesia via Multiple Vital Signs Based on Artificial Neural Networks

This study evaluated the depth of anesthesia (DoA) index using artificial neural networks (ANN) which is performed as the modeling technique. Totally 63-patient data is addressed, for both modeling and testing of 17 and 46 patients, respectively. The empirical mode decomposition (EMD) is utilized to purify between the electroencephalography (EEG) signal and the noise. The filtered EEG signal is subsequently extracted to achieve a sample entropy index by every 5-second signal. Then, it is combined with other mean values of vital signs, that is, electromyography (EMG), heart rate (HR), pulse, systolic blood pressure (SBP), diastolic blood pressure (DBP), and signal quality index (SQI) to evaluate the DoA index as the input. The 5 doctor scores are averaged to obtain an output index. The mean absolute error (MAE) is utilized as the performance evaluation. 10-fold cross-validation is performed in order to generalize the model. The ANN model is compared with the bispectral index (BIS). The results show that the ANN is able to produce lower MAE than BIS. For the correlation coefficient, ANN also has higher value than BIS tested on the 46-patient testing data. Sensitivity analysis and cross-validation method are applied in advance. The results state that EMG has the most effecting parameter, significantly.

[Evaluation of heart rate variability for monitoring the depth of anaesthesia in dogs. Investigations based on total intravenous anaesthesia using propofol alone or in combination with dexmedetomidine or remifentanil]

OBJECTIVE

Evaluation of heart-rate variability (HRV) as an indicator for autonomous activity to monitor anaesthesia in dogs during three different total intravenous anaesthetic protocols and three anaesthetic depth levels as well as before and after electrical nociceptive stimulation.

MATERIAL AND METHODS

Seven beagle dogs (14.3±1.7 kg) were used in a randomised experimental trial with a complete cross-over design. Each dog went through all three anaesthetic protocols, which were propofol alone (group P) and propofol combined with dexmedetomidine (3 µg/kg/h, group PD) or remifentanil (18 µg/kg/h, group PR). Propofol was given using target-controlled infusion. Three anaesthetic depth levels (light, medium, deep) were defined by target concentrations for propofol in the blood and were adapted to the individual animal and treatment (mean of 7, 9 and 11 µg/ml, and in combination with dexmedetomidine or remifentanil, a mean of 3, 5 and 7 µg/ml). During each anaesthetic level, a standardised supramaximal nociceptive electric stimulus (50 Hz, 50 V, 10 ms) was applied medially to the right forearm. The bipolar-derived electrocardiogram (ECG) was recorded continuously. For each anaesthetic depth, the RR-intervals recorded 2 minutes before and after each stimulation were included in the statistical analysis. Using an HRV analytical program (Kubios HRV), the frequency domain HRV-parameters low (LF) and high (HF) frequency and the time-domain HRV-parameters RR-intervals, standard deviation of all RR-intervals (SDNN) and the square root of the mean of the sum of the squares of the differences between consecutive RR-intervals (RMSSD) were determined.

RESULTS

Neither the RR-intervals nor the currently available HRV-parameters which were derived from the RR-intervals were able to discriminate between the different anaesthetic depths levels. Nociception could only be represented by the RR-intervals.

CONCLUSION

Overall, the investigated standard HRV parameters offered no additional information for the monitoring of anaesthetic depths at the investigated, clinically used dose rates.

Identification of Criticality in Neuronal Avalanches: II. A Theoretical and Empirical Investigation of the Driven Case

The observation of apparent power laws in neuronal systems has led to the suggestion that the brain is at, or close to, a critical state and may be a self-organised critical system. Within the framework of self-organised criticality a separation of timescales is thought to be crucial for the observation of power-law dynamics and computational models are often constructed with this property. However, this is not necessarily a characteristic of physiological neural networks-external input does not only occur when the network is at rest/a steady state. In this paper we study a simple neuronal network model driven by a continuous external input (i.e. the model does not have an explicit separation of timescales from seeding the system only when in the quiescent state) and analytically tuned to operate in the region of a critical state (it reaches the critical regime exactly in the absence of input-the case studied in the companion paper to this article). The system displays avalanche dynamics in the form of cascades of neuronal firing separated by periods of silence. We observe partial scale-free behaviour in the distribution of avalanche size for low levels of external input. We analytically derive the distributions of waiting times and investigate their temporal behaviour in relation to different levels of external input, showing that the system's dynamics can exhibit partial long-range temporal correlations. We further show that as the system approaches the critical state by two alternative 'routes', different markers of criticality (partial scale-free behaviour and long-range temporal correlations) are displayed. This suggests that signatures of criticality exhibited by a particular system in close proximity to a critical state are dependent on the region in parameter space at which the system (currently) resides.

Comparison of heart rate variability response in children undergoing elective endotracheal intubation with and without neuromuscular blockade: a randomized controlled trial

BACKGROUND

The routine use of neuromuscular blocking drugs (NMBD) for endotracheal intubation in children is the subject of much controversy. The analysis of heart rate variability (HRV) can reveal information about the functional state of the autonomic nervous system (ANS).

AIM

The purpose of this study was to determine if HRV elucidates differences in the sympathovagal balance of children undergoing elective endo-tracheal intubation with and without neuromuscular blockade (NMB).

METHODS

In this prospective study, 38 children (2-6 years) scheduled for adenotonsillectomy were randomized into two groups to receive fentanyl 2 μg·kg(-1) and propofol 4 mg·kg(-1) , with either mivacurium 0.25 mg·kg(-1) (NMB group) or saline solution (NoNMB group) for anesthesia induction. The same experienced, blinded anesthesiologist performed endotracheal intubation. Heart rate variability, RR intervals, ECG as well as an electroencephalogram were recorded with HRV and BIS XP monitors, respectively. Heart rate variability was analyzed in the frequency domain.

RESULTS

There was no significant difference in HRV changes immediately after mivacurium administration compared with an administration of saline. The groups were comparable for the bispectral index value (NMB 35 [33-41] vs NoNMB 34 [32-42]) during endotracheal intubation. Changes in both the low-frequency power and the low-/high-frequency ratio immediately after endotracheal intubation compared with the unstimulated state before laryngoscopy were significantly higher without NMB (P = 0.015 and P = 0.006, respectively), whereas there was no significant difference with respect to the high-frequency power.

CONCLUSIONS

The stress response during endotracheal intubation in pediatric patients represented by the frequency domain analysis of HRV was found to be higher without NMB. When mivacurium was added to a propofol-fentanyl induction regimen, the ANS alterations during endotracheal intubation decreased significantly.

Effects of isoflurane with and without dexmedetomidine or remifentanil on heart rate variability before and after nociceptive stimulation at different multiples of minimum alveolar concentration in dogs

OBJECTIVE

To evaluate the influence of 3 anesthetic protocols and multiples of minimum alveolar concentration (MAC) on heart rate variability (HRV) with and without nociceptive stimulation in dogs.

ANIMALS

6 healthy adult Beagles.

PROCEDURES

Each dog was anesthetized 3 times: with isoflurane alone, with isoflurane and a constant rate infusion of dexmedetomidine (IsoD; 3 μg/kg/h, IV), and with isoflurane and a constant rate infusion of remifentanil (IsoR; 18 μg/kg/h, IV). Individual MAC was determined via supramaximal electrical stimulation. Sinus rhythm-derived intervals between 2 adjacent R-R intervals were exported from ECG recordings. Selected HRV time and frequency domain variables were obtained (at 2-minute intervals) and analyzed offline with signed rank tests before and after stimulation at 0.75, 1.0, and 1. 5 MAC for each anesthetic session.

RESULTS

The isoflurane session had the overall lowest prestimulation SDNN (SD of all R-R intervals) values. Prestimulation SDNN values decreased significantly with increasing MAC in all sessions. For the IsoD session, SDNN (milliseconds) or high-frequency power (milliseconds(2)) was inversely correlated with MAC (Spearman rank correlation coefficient for both variables, -0.77). In the isoflurane and IsoR sessions, heart rate increased significantly after stimulation. In the IsoD session, poststimulation SDNN was increased significantly, compared with prestimulation values, at 0.75 and 1.0 MAC.

CONCLUSIONS AND CLINICAL RELEVANCE

On the basis of SDNN and high-frequency power values, anesthetic levels between 0.75 and 1.5 MAC within the same anesthetic protocol could be differentiated, but with a large overlap among protocols. Usefulness of standard HRV variables for assessment of anesthetic depth and nociception in dogs is questionable.

Inattentive and hyperactive preschool-age boys have lower sympathetic and higher parasympathetic activity

The presented study aimed to clarify the relationship between the autonomic nervous system and disruptive behaviors among preschoolers. Possible gender differences in autonomic activity were also examined. A total of 88 preschool-age children were enrolled in this study. Autonomic activities were measured by power spectrum analysis of heart rate variability (HRV). The Swanson, Nolan and Pelham parents and teachers rating scale (SNAP-IV) was applied to evaluate each subject's severity of disruptive behavior. The relationship between the HRV results and the SNAP-IV was evaluated by correlation analysis, which disclosed that the scores for inattention, hyperactivity/impulsivity and oppositional defiant disorder showed a negative association with LF % and LF/HF. The above scales, except for the hyperactivity subscale, also showed a positive association with HF. On separating the two genders, only boys showed the above correlations. Preschool-age boys who show more inattentive and hyperactive features have lower sympathetic and higher parasympathetic activity.

Heart rate dynamics during acute pain in newborns

Autonomic nervous system modulation of heart rate is significantly altered during painful procedures in newborns. Most studies investigating pain employed only linear-based analysis methods, thus ignoring the complex, non-linear nature of heart rate control mechanisms. The emergences of dynamic, nonlinear analysis methods enable us to uncover information embedded in the fluctuations of heart rate not otherwise noticeable. Our objective was to examine how cardiac dynamics change in newborns who undergo heel lancing by analyzing linear and nonlinear characteristics of heart rate fluctuations. We used dynamic nonlinear analyses methods to reveal heart rate variability and complexity alterations during painful stimulus in newborns. Poincaré plots were applied to examine the dynamics of the system, sample entropy to investigate the complexity of the system, and detrended fluctuation analysis, to reveal the fractal properties of the system. Heart rate significantly increased (165 vs.123 beats per minute, p < 0.001) while variability decreased. Sample entropy and the quantitative measures of the Poincaré plots (SD1 and SD2) significantly decreased during heel lancing (0.75 vs. 1.0, p < 0.01; 6.4 vs. 12.8, p < 0.001; and 30.4 vs. 50.5, p < 0.01, respectively). Detrended fluctuation analysis showed a significant decrease in the short-term scaling exponent α1 (1.06 vs. 1.3, p < 0.001), and an increase in the long-term scaling exponent α2 (1.5 vs. 1.1, p < 0.001). Our results indicate altered complexity of heart rate variability during painful stimulus in newborns and disruption of the mechanisms that regularly control it. Such alterations resemble certain pathological conditions and may represent stress reaction.

Heart rate variability: a diagnostic and prognostic tool in anesthesia and intensive care

The autonomic nervous system (ANS) plays an important role in the human response to various internal and external stimuli, which can modify homeostasis, and exerts a tight control on essential functions such as circulation, respiration, thermoregulation and hormonal secretion. ANS dysfunction may complicate the perioperative course in the surgical patient undergoing anesthesia, increasing morbidity and mortality, and, therefore, it should be considered as an additional risk factor during pre-operative evaluation. Furthermore, ANS dysfunction may complicate the clinical course of critically ill patients admitted to intensive care units, in the case of trauma, sepsis, neurologic disorders and cardiovascular diseases, and its occurrence adversely affects the outcome. In the care of these patients, the assessment of autonomic function may provide useful information concerning pathophysiology, risk stratification, early prognosis prediction and treatment strategies. Given the role of ANS in the maintenance of systemic homeostasis, anesthesiologists and intensivists should recognize as critical the evaluation of ANS function. Measurement of heart rate variability (HRV) is an easily accessible window into autonomic activity. It is a low-cost, non-invasive and simple to perform method reflecting the balance of the ANS regulation of the heart rate and offers the opportunity to detect the presence of autonomic neuropathy complicating several illnesses. The present review provides anesthesiologists and intensivists with a comprehensive summary of the possible clinical implications of HRV measurements, suggesting that autonomic dysfunction testing could potentially represent a diagnostic and prognostic tool in the care of patients both in the perioperative setting as well as in the critical care arena.

Autonomic Nervous System and Cardiovascular Disease

Because anesthesia affects the integrity of the autonomic nervous system, anesthesiologists use vital signs to maintain respiratory and circulatory homeostasis. However, patients with genetic predispositions or with autonomic dysfunctions are at risk of severe complications from anesthesia. For these patients, the monitoring of vital signs may not give sufficient warning to avoid complications. The development of methods to measure autonomic tone could be of interest to anesthesiologists because they could warn of changes in autonomic tone before vital signs are affected. New noninvasive methods are being developed to obtain measurements of parasympathetic and sympathetic output allowing for the monitoring of perioperative autonomic tone. These measurements are based on analysis of heart rate and blood pressure variability. In this report, the principals of the analysis of heart rate and blood pressure variability will be explained and the usefulness of these methods to anesthesiologists will be discussed.

Propofol versus alfentanil to prevent movement responses during uterine curettage

BACKGROUND

Propofol anaesthesia is frequently associated with movement responses in non-paralysed patients. Opioids decrease the probability of movement during noxious stimulation. Heart rate variability and frontal electromyography (EMG), which are related to subcortical functions, may be more closely related than surface electroencephalography (EEG) to movement responses to noxious stimulation.

METHODS

Eighty-two patients scheduled for uterine dilatation and curettage were randomized to receive at the first intra-operative movement either a supplemental alfentanil bolus, 0.5 mg intravenously, or a supplemental propofol bolus, 0.7 mg/kg intravenously. The incidences of recurring movement during the procedure were compared between the two groups. The associations of a measure of heart rate variability (Anemon index), heart rate, EMG and two EEG variables with movement responses were evaluated.

RESULTS

The incidences of recurring movement were 73% and 38% in the alfentanil and propofol groups, respectively [difference, 35%; 95% confidence interval, 9-56%; P= 0.014 between the groups). The Anemon index, heart rate, EMG and surface EEG variables displayed mainly reactive associations with movement responses.

CONCLUSION

During uterine curettage under propofol-alfentanil-nitrous oxide anaesthesia, a propofol bolus of 0.7 mg/kg was more effective in preventing the recurrence of movement responses than an alfentanil bolus of 0.5 mg. Several physiological variables may be used to track significant arousal reactions, but not to predict them.

Heart rate variability does not discriminate between different levels of haemodynamic responsiveness during surgical anaesthesia †

BACKGROUND

Hypnotic depth but not haemodynamic responsiveness is measured with EEG-based monitors. In this study we compared heart rate variability (HRV) in unstimulated patients and stimulation-induced HRV at different levels of anaesthesia.

METHODS

A total of 95 ASA I or II patients were randomly assigned to five groups (Group 1: BIS 45(5), remifentanil 1 ng ml(-1); Group 2: BIS 45(5), remifentanil 2 ng ml(-1); Group 3: BIS 45(5), remifentanil 4 ng ml(-1); Group 4: BIS 30(5), remifentanil 2 ng ml(-1); Group 5: BIS 60(5), remifentanil 2 ng ml(-1)). A time- and frequency-domain analysis of the RR interval (RRI) from the electrocardiogram was performed. HRV before induction, before and after a 5 s tetanic stimulus of the ulnar nerve, and before and after tracheal intubation was compared between groups, between stimuli, and between responders to intubation [systolic arterial pressure (SAP) increase >20 mm Hg, a maximal heart rate (HR) after intubation >90 min(-1) or both] and non-responders (anova).

RESULTS

Induction of anaesthesia significantly lowered HR and HRV. Mean RRI before stimulation was higher in G3 than in G1, G2, and G4 (P < 0.001), whereas the other HRV parameters were similar. Intubation induced a greater HRV response than tetanic stimulation. The mean RRI after intubation was lower in G3 compared with the other groups and the sd of the RRI after tetanic stimulation was lower in G3 compared with G5. Otherwise, unstimulated HRV and stimulation-induced HRV were similar in responders and non-responders.

CONCLUSION

HRV parameters discriminate between awake and general anaesthesia, are different after tracheal intubation and a 5 s ulnar nerve stimulation, but do not discriminate between different levels of haemodynamic responsiveness during surgical anaesthesia.

Loss of complexity characterizes the heart rate response to experimental hemorrhagic shock in swine

OBJECTIVE

To improve our ability to identify physiologic deterioration caused by critical illness, we applied nonlinear and frequency-domain analytical methods to R-to-R interval (RRI) and systolic arterial pressure (SAP) time series during hemorrhagic shock.

DESIGN

Prospective, randomized, controlled trial.

SETTING

Animal laboratory of a government research institute.

SUBJECTS

Twenty swine (weight 36.4+/-0.11 kg).

INTERVENTIONS

Fixed-volume hemorrhage followed by resuscitation; off-line analysis of RRI and SAP data.

MEASUREMENTS AND MAIN RESULTS

Anesthetized swine (shock group, n=12) underwent withdrawal of 30 mL/kg blood in 10 mL/kg decrements. A control group (n=8) received maintenance fluids only. Electrocardiogram and arterial pressure waveforms were acquired at 500 Hz. Eight hundred-beat data sets were analyzed at six time points: at baseline, after each blood withdrawal, after lactated Ringer's resuscitation, and after infusion of shed blood. Nonlinear methods were used to estimate the complexity (approximate entropy, sample entropy, Lempel-Ziv entropy, normalized entropy of symbol dynamics), RRI bits per word, and fractal dimension by curve lengths and by dispersion analysis of the RRI and SAP time series. Fast Fourier transformation was used to measure the high-frequency and low-frequency powers of RRI and SAP. Baroreflex sensitivity was assessed in the time domain with the sequence method. Hemorrhagic shock caused decreases in RRI complexity as quantified by approximate entropy, sample entropy, and symbol dynamics; these changes were reversed by resuscitation. Similar but statistically insignificant changes in fractal dimension by curve lengths were seen. RRI high-frequency power decreased with hemorrhagic shock-indicating withdrawal of vagal cardiac input-and was restored by resuscitation. Similar changes in baroreflex sensitivity were seen. Hemorrhagic shock did not affect SAP complexity.

CONCLUSIONS

Hemorrhagic shock caused a reversible decrease in RRI complexity; these changes may be mediated by changes in vagal cardiac control. Assessment of RRI complexity may permit identification of casualties with hemorrhagic shock.

Smoothed Periodogram of Oxyhemoglobin Saturation by Pulse Oximetry in Sleep Apnea Syndrome

BACKGROUND

Variability of oxyhemoglobin saturation (Spo(2)) during sleep has been utilized as a diagnostic index for sleep apnea. Spectral analysis with its graphical presentation, the periodogram, is an approach for measuring such variability. This work examined the parameters on a smoothed periodogram created from series data for Spo(2) obtained by pulse oximetry during a sleep study.

DESIGN AND RESULTS

Spo(2) was recorded during polysomnography study of 273 subjects. Clinical data of subjects were collected retrospectively. A novel automated algorithm was created to measure the low-frequency (< 0.1 Hz) peak and the slope of spectral density vs frequency in the frequency region of 0.1 to 0.5 Hz (slope(0.1-0.5)). Two successive modified Daniell smoothers with span lengths of 3 to 121 in odd numbers were applied to determine the effect of smoothing on these parameters. slope(0.1-0.5) was least affected by smoothing and had a sensitivity of 78% and a specificity of 80% in diagnosing sleep apnea defined by a value of apnea-hypopnea index >/= 5. Combining slope(0.1-0.5) with parameters of the low-frequency peak enlarged the area under the receiver operating characteristic curve. A composite indicator comprised of slope(0.1-0.5) and ratio of the area under the curve of the low-frequency peak to that of whole periodogram (AUCratio) had a positive likelihood ratio of 15.25 in identifying patients with moderate-to-severe obstructive sleep apnea. The algorithm was validated in another 206 patients undergoing polysomnographic studies.

CONCLUSIONS

These analytical results demonstrate that the smoothed periodogram of Spo(2) is a useful tool for screening subjects with sleep apnea.

Complex systems and the technology of variability analysis

Characteristic patterns of variation over time, namely rhythms, represent a defining feature of complex systems, one that is synonymous with life. Despite the intrinsic dynamic, interdependent and nonlinear relationships of their parts, complex biological systems exhibit robust systemic stability. Applied to critical care, it is the systemic properties of the host response to a physiological insult that manifest as health or illness and determine outcome in our patients. Variability analysis provides a novel technology with which to evaluate the overall properties of a complex system. This review highlights the means by which we scientifically measure variation, including analyses of overall variation (time domain analysis, frequency distribution, spectral power), frequency contribution (spectral analysis), scale invariant (fractal) behaviour (detrended fluctuation and power law analysis) and regularity (approximate and multiscale entropy). Each technique is presented with a definition, interpretation, clinical application, advantages, limitations and summary of its calculation. The ubiquitous association between altered variability and illness is highlighted, followed by an analysis of how variability analysis may significantly improve prognostication of severity of illness and guide therapeutic intervention in critically ill patients.