Artifacting
What You Will Learn
Heart rate variability recordings are vulnerable to artifacts that can dramatically distort your measurements and compromise clinical decisions. In this unit, you will learn to identify artifacts in raw signals and interbeat interval (IBI) data, understand why prevention far surpasses correction, and master two primary correction approaches for handling artifacts when they occur. Whether you are conducting publication-quality research or performing routine clinical assessments, this knowledge will ensure your data accurately reflects your client's physiology.
This unit covers artifact recognition across four categories: ectopic beats, movement, environmental interference, and hardware limitations. A dedicated section explains why a real premature beat becomes an artifact only inside a normal-to-normal interval series, how to tell true ectopy from a detection error, and when correction turns into invention. You will develop practical skills for each correction approach—including the use of Kubios noise detection, beat correction, and deviation threshold settings—and learn when to apply them. By the end, you will be able to make informed decisions about data quality that support valid clinical conclusions.
HRV recordings are vulnerable to diverse artifacts that can produce interbeat interval (IBI) values that are too long or too short. The IBI, also called the NN (normal-to-normal) interval when artifacts have been removed, represents the time between consecutive R-spikes in the ECG waveform. Artifacts are false values produced by the client's body (ectopic beats), the client's actions (movement), the environment (line current), and hardware limitations (light leakage).

The clinical stakes are high: a single artifactual IBI value in a 2-minute epoch can markedly distort time-domain indices like SDNN, the standard deviation of normal-to-normal intervals (Berntson et al., 1997). This sensitivity to outliers (i.e., extreme values) makes artifact detection and correction essential for valid HRV assessment. Consider that SDNN reflects overall autonomic function—a distorted value could lead you to overestimate or underestimate your client's regulatory capacity.

Client assessment and training are only as good as your data. As the familiar saying goes, "garbage in, garbage out." This principle applies whether you are working with a veteran managing PTSD symptoms, an athlete optimizing performance, or a cardiac rehabilitation patient rebuilding autonomic resilience.

We recommend the "close the barn door before the horse escapes" strategy to ensure data integrity. Prevent artifacts before you record data so that you don't have to remove them later. This proactive approach—proper sensor placement, controlled recording conditions, and client preparation—saves time and produces more reliable results than even the most sophisticated correction algorithms. Please review our BVP Hardware and ECG Hardware units for a review of artifacts and strategies to mitigate them.
Manually inspect the raw signal and IBI values to ensure that artifacts do not contaminate your measurements. You can do this within programs like Kubios Scientific. Both provide visualization tools that display the ECG or PPG waveforms and computed IBI values, allowing you to spot discrepancies between what the sensor recorded and what the software calculated.

The relationship between artifacts in the raw signal and false IBI values can be challenging for beginners to grasp. For ECG recording, the software analyzes the raw signal to detect R-spikes—the tall, sharp peaks that mark ventricular depolarization—and measures the time intervals between them. When artifacts distort this process, the resulting IBI values no longer reflect true cardiac timing.
IBIs are stored in a time-stamped table of consecutive values, similar to an Excel spreadsheet column. The ECG waveform circled above at 29.171 milliseconds corresponds to an IBI value of 818.359 milliseconds at the bottom of the table below. Understanding this correspondence is essential for tracing artifacts from their source in the raw signal to their impact on computed metrics.

Artifacts can cause the software to "see" extra beats, resulting in abnormally short IBIs. They can also cause it to miss beats, resulting in abnormally long IBIs. These outliers stand out when you review a table of IBI values. For example, we used the color red to identify a short IBI of 274.9 milliseconds produced by an extra beat—far below the typical resting range of 600-1200 milliseconds.

Conversely, software missed a beat before 00:00:31, resulting in a long IBI of 1334.4 milliseconds. When software fails to detect an R-spike, it measures the interval from one detected beat to the next, effectively doubling (or more) the true interval.

Short and long values stand out compared to typical adult values. Developing an eye for these outliers is a fundamental skill for any clinician working with HRV data.
BCIA Blueprint Coverage
This unit addresses III. HRV Instrumentation: E. Artifacting Strategies from the BCIA certification blueprint. Whether you are preparing for certification or expanding your clinical competencies, this material provides the foundation for defensible data quality practices.
Topics include Artifacting Overview, Plan A (Segment Removal and Detrending) and Plan B (Manual Correction and Detrending). Each approach suits different clinical scenarios, and knowing when to apply each one will improve both your efficiency and your data integrity.
🎧 Listen to the Full Chapter Lecture
Artifacting Overview: Understanding the Process
This section introduces the fundamental principles of artifacting and explains why certain correction approaches are inappropriate. You will learn what artifacting is, when data should be discarded, and why the "splicing" method fails to preserve data integrity.
Artifacting is the process of removing or replacing false values found within a data record. Think of it as quality control for your HRV data: you are ensuring that every value in your dataset represents a genuine cardiac event rather than noise or error. This process directly parallels the artifact rejection procedures familiar to clinicians who practice neurofeedback.

Researchers may discard the first and last minutes of data due to client adjustment to the monitoring environment and performance degradation after a long recording session. The initial adjustment period reflects the client settling into the recording situation—electrode gel stabilizing, the client becoming accustomed to the sensors, and initial anxiety subsiding. Toward the end of extended sessions, fatigue and reduced attention can affect data quality, particularly if the client shifts position or loses focus.

The Splicing Method: What Not to Do
Never remove a part of an epoch (a defined data collection period, typically 5 minutes for standard HRV assessment) and rejoin the remaining sections. This approach destroys the temporal relationship between heartbeats, which contains meaningful physiological information (Combatalade, 2010). When you cut out a section and splice the remaining parts together, you create artificial transitions that never occurred in the client's actual cardiac rhythm.
The consequences of splicing are particularly severe for frequency-domain analysis. HRV analysis algorithms expect a continuous time series with natural beat-to-beat relationships. When you splice data together, they interpret the artificial junction as a sudden physiological change, potentially creating spurious power in specific frequency bands that misrepresent your client's autonomic function.
Ectopic Beats: When a Real Heartbeat Becomes an Artifact
This section gathers everything you need to know about ectopic beats in one place, before you choose between Plan A and Plan B. You will learn how a real heartbeat becomes an analytic artifact, how to recognize premature ventricular and atrial beats on the strip, how to separate true ectopy from detector error, how much damage a single ectopic beat does to your metrics, and where careful repair turns into invention.
An ECG shows the heart's electrical activity, but HR and HRV throw most of that picture away and keep only the timing between beats (Kligfield et al., 2007; Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology [Task Force], 1996). That simplification works well until the heart adds an early beat. Every recording produces two lists of numbers, and confusing them causes most HRV errors. The RR interval list holds the time between every detected R spike, whatever produced it, while the NN interval list holds only the time between normal sinus beats (Task Force, 1996). Software often labels both lists as RR, so an RR file can quietly contain ectopic beats, missed beats, extra markers, and mistimed markers, even though every standard HRV metric assumes you handed it the second list (Berntson et al., 1997; Lipponen & Tarvainen, 2019; Weinstein et al., 2026).
One recording therefore asks you to do two jobs at once. The clinical job asks what the heart did, and the measurement job asks how much the sinus node varied. Keep three layers of the same recording and both jobs stay honest: the raw ECG stays untouched and answers clinical questions, an annotation layer records what you or the software called each beat, and a separate NN copy holds the estimated sinus timing you send to your HRV software (Laborde et al., 2017; Peltola, 2012). That habit protects the patient as well as the statistics, because ectopic burden, the count or percentage of ectopic beats in a record, carries clinical meaning of its own (Guichard et al., 2022; Marcus, 2020). Never let cleanup erase the rhythm history of the person in front of you.
Recognizing Ectopic Beats on the Strip
Premature ventricular contractions (PVCs) originate from an ectopic focus in the ventricles rather than the sinoatrial node. The heart's normal rhythm is briefly disrupted: an early ventricular depolarization produces a short interval, followed by a compensatory pause while the cardiac electrical system resets to its normal rhythm. This produces a characteristic short-long pattern that is recognizable once you know what to look for.

Three features identify a PVC without any measurement at all. The beat arrives early, its QRS complex, the waveform produced by ventricular depolarization, is wide at roughly 120 ms or more, and no P wave precedes it (Marcus, 2020; Zeppenfeld et al., 2022). Because the ectopic impulse usually leaves the sinoatrial node ticking undisturbed, the next sinus impulse arrives while the ventricles are still refractory and produces nothing visible, which gives you a compensatory pause in which the short interval plus the pause add up to about two normal cycles (Marcus, 2020). Two exceptions break that arithmetic and are worth checking before you correct anything. An interpolated PVC squeezes between two sinus beats with no pause at all, and a fusion beat blends sinus and ventricular activation into a complex that looks intermediate (Marcus, 2020). Confirm that a full compensatory pause is actually present before you assume that a repair can simply redistribute two cycles worth of elapsed time.
Premature atrial complexes (PACs) are initiated by an ectopic pacemaker in the atria rather than the ventricles. Because a PAC usually resets the sinus node and produces no full compensatory pause, there is no simple time-preserving way to remove it through manual artifacting (Clinical EEG Interpretation, 2018; Mond & Haqqani, 2019). PACs are therefore more challenging to correct than the symmetrical short-long pattern that makes an isolated PVC relatively straightforward to address.

A PAC starts in atrial tissue outside the SA node, so its P wave, the deflection representing atrial depolarization, differs in shape or direction from its neighbors, while the QRS that follows usually looks narrow and ordinary (Mond & Haqqani, 2019). A conducted PAC usually resets the sinus node, which pushes the next sinus beat later than scheduled and produces a noncompensatory pause spanning less than two normal cycles (Mond & Haqqani, 2019).
Two variations cause trouble. When the ventricular conduction system has not yet recovered, aberrant conduction gives the PAC a wide QRS that mimics a PVC, and when the atrioventricular node has not recovered, no QRS follows at all, leaving a blocked PAC whose early P wave hides inside the previous T wave (Mond & Haqqani, 2019). The blocked PAC is the one that fools software, because a detector watching only R spikes reads that hidden P wave as a dropped beat or a sinus pause (Mond & Haqqani, 2019; Peltola, 2012).
Telling True Ectopy From Detector Error
A tachogram will show you a short-long pattern immediately, but it will not tell you what caused it (Peltola, 2012; Weinstein et al., 2026). Ectopy, a missed R spike, a double-counted T wave, and a client shifting in the chair can all produce similar spikes, so look at the waveform underneath every candidate before you decide (Laborde et al., 2017; Peltola, 2012). Detector errors leave three arithmetic signatures, and each one calls for a different fix.
A missed beat leaves one interval roughly twice as long as its neighbors, a false extra marker leaves two short intervals that add up to one normal interval, and a marker placed slightly off inside the QRS keeps the beat count correct but the timing wrong (Lipponen & Tarvainen, 2019; Weinstein et al., 2026).
Movement and poor electrode contact look different again, because they usually distort several beats in a row, shift the baseline, and produce shapes that no cardiac event would make (Friesen et al., 1990; Kligfield et al., 2007).
Repeating Ectopy
Ectopic beats often arrive in patterns, and the pattern changes what you can do about them. Bigeminy alternates sinus and ectopic beats.


A couplet is a pair of consecutive ectopic beats with no sinus beat between them, three in a row make a triplet, and longer runs belong in a rhythm report rather than an HRV analysis (Sattar & Hashmi, 2025; Zeppenfeld et al., 2022).

Couplets are more disruptive than isolated ectopy for a structural reason: an isolated PVC corrupts the ectopic interval and the compensatory pause that follows it, whereas a couplet adds the interval between the two ectopic beats and yields at least three consecutive corrupted intervals (Peltola, 2012). Interpolation is appropriate for occasional ectopic beats, but long stretches of consecutive corrupted intervals cannot reasonably be rebuilt that way, and excluding the affected segment becomes the more defensible choice (Peltola, 2012).
How Much Damage a Single Ectopic Beat Does
HRV metrics exist to magnify differences between intervals, which is exactly why a few bad intervals can swamp them (Choi & Shin, 2018; Stapelberg et al., 2018). SDNN widens because the short beat and the long pause stretch the distribution, while the RMSSD, the root mean square of successive differences, and pNN50, the percentage of successive intervals differing by more than 50 ms, react even more strongly because they measure the jump from each beat to the next (Bourdillon et al., 2022; Choi & Shin, 2018; Peltola, 2012).
The published magnitudes are striking. Bourdillon et al. (2022) added a single artifact and watched RMSSD rise by 413% supine and 269% standing while frequency-domain values barely moved, and Choi and Shin (2018) found significant distortion in every analysis domain at ectopic rates below 1% when no interpolation was applied.
Your cleanup setting matters as much as the ectopy itself, because Bassi-Dibai et al. (2026) showed that the strictest correction filter in a widely used analysis package changed time-domain, frequency-domain, and nonlinear results in the same recordings, which makes your filter selection part of your method rather than a buried software preference.
Repair Rather Than Delete
Deleting an interval does not restore sinus rhythm. It shortens the record, shifts the time axis, and stitches together two beats that were never neighbors in the heart, which damages frequency-domain estimates most of all because spectral analysis depends on continuous, correctly spaced timing (Clifford & Tarassenko, 2005; Peltola, 2012; Salo et al., 2001).
Interpolation, the estimation of rejected interval values from the beats around them, is usually the better choice, although it is not neutral, because linear, nearest-neighbor, spline, and model-based methods each make different assumptions about what the sinus node would have done (Choi & Shin, 2018; Lippman et al., 1994; Rincon Soler et al., 2018).
An isolated PVC with a full compensatory pause gives you an unusual advantage, since its short interval and pause together span about two sinus cycles, so a time-preserving replacement can redistribute that elapsed time into two estimated sinus intervals and keep the clock honest (Mateo & Laguna, 2003; Peltola, 2012). A PAC that resets the sinus node offers no such anchor, so its replacement values remain an estimate rather than a recovered measurement.
One post-ectopic effect is genuine physiology rather than error: heart rate turbulence (HRT), a brief speeding of the sinus rhythm after a PVC followed by a slower return, reflects the baroreflex response to a mechanically weak ectopic contraction, so decide in advance whether you will exclude it, model it, or study it rather than smoothing it away (Bauer et al., 2008).
Knowing When to Stop Correcting
Interpolation works when the damage is isolated and the surrounding rhythm is stable enough to model, and it stops working as ectopy becomes frequent, clustered, or varied in shape (Peltola, 2012; Tarkiainen et al., 2007). Short recordings need the strictest scrutiny, because editing even a small proportion of intervals, under 5%, can shift HF and LF power as well as pNN50 and RMSSD, and most investigators require at least 80% normal intervals before trusting a frequency-domain analysis (Peltola, 2012; Salo et al., 2001).
Do not expect a universal cutoff, since published thresholds differ with the metric, record length, ectopic pattern, correction method, and population, so set a rule in advance and justify it for your own design (Bassi-Dibai et al., 2026; Choi & Shin, 2018; Stapelberg et al., 2018). When the algorithm must invent more sinus rhythm than the ECG actually recorded, stop correcting, exclude the segment, and report the ectopic burden as a finding in its own right (Peltola, 2012; Task Force, 1996).
Referral follows the patient rather than the preprocessing problem, so syncope, exertional symptoms, complex or sustained ventricular ectopy, and evidence of structural disease deserve cardiac evaluation whether or not anyone plans an HRV analysis (Guichard et al., 2022; Marcus, 2020).
Check Your Understanding: Ectopic Beats
- What distinguishes an RR interval series from an NN interval series, and why does the difference matter for HRV analysis?
- Which three layers of a recording should you preserve, and what question does each one answer?
- How does the pause after a PVC differ from the pause after a PAC, and how does that difference change your repair options?
- Why is a blocked PAC especially likely to fool a detector that watches only R spikes?
- Why does a couplet resist interpolation more stubbornly than an isolated PVC?
- What conditions should prompt you to stop correcting a segment and exclude it instead?
Plan A: Segment Removal and Detrending
This section covers Plan A: when you have sufficient clean data, analyze only the artifact-free segment. You will learn how Kubios Scientific cleans data in two stages (noise detection and beat correction), how deviation thresholds and the automatic correction algorithm work, how cubic spline interpolation repairs identified intervals, and how detrending mathematically restores stationarity using the smoothness priors method.
When an epoch is at least 6 minutes long and the artifact occurs in the first or last 3 minutes, you can analyze the clean 3-minute segment (R. Gevirtz, personal communication, 2018). This approach preserves the temporal integrity of your data while avoiding the contaminated portion entirely. The clean segment contains a set of IBI values measured over time—what we call a time series. Visualize this as an Excel column where each cell contains a millisecond value (e.g., 818.3 ms) representing the time between consecutive heartbeats.

Clearing the Static with Noise Detection
Kubios uses a feature called noise detection to act as a high-powered filter for your heart rate data. If a segment of your recording becomes too messy—perhaps because a sensor slipped or you moved suddenly—the software marks that entire block as unreadable. Instead of trying to guess what happened during those distorted moments, it simply sets that section aside so it won't ruin your final results. This tool is essential for handling external interference and ensuring that only the highest quality data makes it through to the final analysis. By removing these large, unreadable chunks, the software protects the accuracy of your heart rate variability scores.

Kubios allows you to select a color-coded clean segment for further analysis, skipping corrupted values.

Polishing the Details with Beat Correction
While noise detection handles the big, messy gaps, beat correction acts like a precision editor for the smaller details. It doesn't throw away whole sections of data but instead looks for tiny, individual hiccups like a single "extra" heartbeat. These brief stutters are often natural or caused by a minor sensor glitch, and the software is smart enough to recognize them. Rather than deleting these moments, beat correction surgically adjusts the timing of those specific intervals to keep the rhythm mathematically consistent. This allows you to keep as much of your heart’s story as possible without one weird beat throwing everything off.
A Two-Step Cleaning Process
These two features work together in a specific sequence to ensure your data is as clean as possible. First, the noise detection clears out the "trash" caused by movement or poor signal quality. Once those big gaps are handled, the beat correction tool goes through the remaining data to smooth out any minor, individual quirks. This two-stage process ensures that the software isn't trying to "fix" something that is actually just pure noise. The result is a polished, reliable picture of your heart's health that reflects your actual physiology rather than technical errors.
Decoding the Noise Level Graph
To see exactly how much work the software did, you can look at the "Noise Level" graph in your final Kubios report. This graph usually appears as a series of bars or a shaded area along the bottom of your heart rate timeline. When the bar spikes or changes color, it indicates a period where the signal quality dropped and the software had to intervene. If you see large sections marked as noise, it might mean your heart rate strap was too loose or your skin was too dry during the recording. Understanding this graph helps you decide if your data is "clean" enough to trust or if you need to try a more stable recording next time.
The Primary Stage of Error Detection
Before any mathematical transformations can occur, the researcher must act as a guardian of the data by identifying and removing outliers. Outliers are any data points that are so far outside the expected range that they are likely mistakes rather than real biological events. In a practical sense, this means scanning the collection of heartbeats for values that are physically impossible.
For a beginner, this involves setting a window of acceptable time between beats—usually between 400 and 1,300 milliseconds. If the software detects a value of 200 milliseconds, these data are likely "noise," which is electrical interference from the environment or the body that the computer mistakenly labeled as a heartbeat. Conversely, a value of 2,000 milliseconds represents data that are likely "missed beats," where the heart actually fired, but the sensor failed to record it. Removing these spikes first is essential because they represent infinite frequency changes that would confuse the software during the next stage.
How Kubios Determines When a Beat Is an Artifact
Kubios HRV Scientific identifies abnormal beats using a statistical method developed by Lipponen and Tarvainen that analyzes the RR-interval time series, which is the sequence of time intervals between successive heartbeats (Lipponen & Tarvainen, 2019). In electrocardiography (ECG), an RR interval is the time between two consecutive R waves, the large spikes that represent ventricular depolarization. In HRV analysis, these intervals form a tachogram, a time series that reflects how heart period fluctuates over time.
The Lipponen–Tarvainen algorithm does not rely on the ECG waveform shape to identify abnormal beats. Instead, it compares each interval to the local rhythm, which is estimated using a short moving window that calculates the expected interval duration based on neighboring beats. If the interval deviates too far from this expected value, Kubios assumes that something abnormal occurred, such as a premature beat, a missed detection, or motion artifact. The degree of deviation that is allowed before a beat is classified as abnormal is called the deviation threshold.
A deviation threshold defines the maximum difference between the observed RR interval and the expected interval that Kubios will accept as normal. When the difference exceeds this threshold, the software labels the beat as an artifact or ectopic event and replaces the abnormal interval using interpolation so that the HRV calculations are based on an estimated sinus rhythm (Lipponen & Tarvainen, 2019).
The Meaning of Deviation Thresholds in Kubios
Kubios provides several artifact-correction strength settings. Each setting corresponds to a specific deviation threshold expressed in milliseconds (ms). One millisecond equals one thousandth of a second. For reference, a typical resting RR interval in a healthy adult is around 600–1200 ms, corresponding to heart rates of about 50–100 beats per minute.
The threshold therefore represents how large a deviation Kubios will tolerate relative to the local rhythm before classifying a beat as abnormal. Larger thresholds allow greater variation and therefore detect fewer artifacts, while smaller thresholds are stricter and detect more artifacts.
No Correction
When the None option is selected, Kubios performs no automatic artifact detection or correction. Every interval detected in the RR series is accepted as valid, regardless of how abnormal it may appear. This setting is rarely appropriate unless the recording has already been manually inspected and edited. Because HRV metrics are highly sensitive to abnormal intervals, even a single premature beat can inflate measures such as RMSSD or SDNN if it is left uncorrected (Shaffer & Ginsberg, 2017).

Very Low Correction (±450 ms)
The Very Low correction level allows intervals to deviate by approximately ±450 milliseconds from the expected interval before Kubios flags them as artifacts. In practical terms, this means that an RR interval could be nearly half a second longer or shorter than the surrounding rhythm and still be treated as normal.
This extremely permissive threshold is intended for data sets that contain very large natural variability or when the user wishes to minimize automatic corrections. Because the threshold is so wide, only the most obvious abnormalities, such as missed beats or extremely premature contractions, are detected.

Low Correction (±350 ms)
The Low setting reduces the allowable deviation to approximately ±350 milliseconds. This means that if the expected RR interval is 900 ms, Kubios will flag intervals shorter than roughly 550 ms or longer than about 1250 ms.
At this level, Kubios detects clear ectopic beats or missed detections but still allows relatively large physiological variability. This setting is commonly used for high-quality ECG recordings in which only occasional artifacts are expected.

Medium Correction (±250 ms)
The Medium correction level allows a deviation of approximately ±250 milliseconds from the expected interval. In the previous example of a 900 ms rhythm, intervals shorter than about 650 ms or longer than about 1150 ms would be classified as abnormal.
This level represents a balance between sensitivity and conservatism. It detects most premature beats and many motion artifacts while preserving legitimate physiological variability. For many research recordings, particularly those collected under resting conditions, the Medium setting is often recommended.
Strong Correction (±150 ms)
The Strong setting tightens the allowable deviation to approximately ±150 milliseconds. If the expected interval is 900 ms, Kubios will flag intervals shorter than about 750 ms or longer than about 1050 ms.
This stricter threshold captures smaller irregularities in the RR sequence and is often used when recordings contain significant noise or motion artifact. However, it can occasionally misclassify legitimate physiological changes as artifacts, particularly during large respiratory sinus arrhythmia, which is the natural fluctuation in heart rate associated with breathing.

Very Strong Correction (±50 ms)
The Very Strong setting allows only ±50 milliseconds of deviation from the expected interval before Kubios marks a beat as abnormal. In a rhythm of 900 ms, intervals shorter than 850 ms or longer than 950 ms would be corrected.
This extremely strict threshold aggressively removes irregular intervals and is typically used when analyzing noisy signals from wearable photoplethysmography (PPG) devices or recordings contaminated by substantial motion artifact. Because the threshold is so narrow, it can also remove genuine physiological variability if applied to clean data.

Why the Threshold Matters for HRV Analysis
Heart rate variability metrics quantify the natural fluctuations in cardiac timing produced by autonomic regulation of the sinoatrial node. When abnormal beats remain in the data, they artificially increase variability measures such as the RMSSD (root mean square of successive differences) and the SDNN (standard deviation of normal intervals). Conversely, overly aggressive artifact correction can remove legitimate variability and reduce these metrics.
The deviation threshold therefore determines how Kubios balances preserving physiological variability against removing abnormal intervals. Selecting the appropriate setting requires considering the recording quality, the measurement method (ECG versus PPG), and the expected magnitude of heart rate fluctuations.
Practical Perspective
In practical HRV analysis, researchers typically begin with Low or Medium correction for clean ECG recordings. Stronger settings are used when signals contain substantial artifact, particularly in wearable devices where motion interference is common. After correction, analysts should visually inspect the tachogram to ensure that genuine physiological rhythms have not been removed.
Understanding these deviation thresholds helps the user interpret how aggressively Kubios is editing the RR-interval series and how that editing may influence the final HRV metrics.
Automatic Beat Correction
The automatic system provides a standardized, repeatable way to ensure that every file in a study is treated exactly the same way, which is vital for the scientific validity of research.

How Kubios' Automatic Algorithm Thinks
The reason Kubios developers favor the automatic method over older, manual thresholds is its ability to adapt to the unique rhythm of the individual. Traditional methods often use a fixed threshold, such as a rule that flags any beat that is 20% different from its neighbor. However, a fixed rule like this can be problematic; it might be too strict for an athlete with a very dynamic heart or too loose for an elderly patient with a very steady rhythm. The automatic algorithm instead uses a time-varying approach, which means it constantly recalculates what a "normal" change looks like based on the heartbeats immediately surrounding the current one. It essentially learns the "language" of that specific heart in real-time, allowing it to distinguish between a natural, healthy jump in heart rate and a technical glitch.
Precision in Identifying Ectopic Beats
One of the most difficult challenges for a beginner is identifying ectopic beats, which are premature heartbeats that originate from a different part of the heart than the standard pacemaker. While these are real physical events, they act as statistical outliers that can ruin heart rate variability calculations. Kubios recommends the automatic setting because it is incredibly precise at identifying these specific patterns. Kubios does not examine waveform morphology (i.e., shape). Instead it examines interval lengths. The software looks for a very short interval followed immediately by a long, "compensatory" pause— a classic signature of an ectopic beat. By automatically detecting these, the software can remove the influence of the irregular beat while accurately estimating what the heart's rhythm would have been if the interruption hadn't occurred.
The Role of Cubic Spline Interpolation
Once the automatic system identifies an error, it doesn't just leave a hole in the data. Kubios uses a sophisticated technique called cubic spline interpolation to heal the gap. To understand this, imagine you have a series of dots on a page representing heartbeats, and one dot is missing. Interpolation is the process of drawing a smooth, curved line between the existing dots to estimate where the missing dot should have been. A cubic spline is just a specific mathematical way of drawing that curve so that it remains perfectly smooth and follows the natural acceleration and deceleration of the heart. Kubios recommends the automatic system because it performs this complex math instantly and accurately, ensuring that the final data set remains continuous and ready for analysis.
Practical Reliability for Modern Research
Ultimately, the recommendation for automatic correction comes down to its proven reliability across thousands of different types of recordings. In formal testing, this algorithm has shown nearly 100% accuracy in finding missed or extra beats, which are the most common errors in heart rate data.
For a beginner, this provides a safety net that allows you to focus on interpreting your results rather than worrying about whether you manually clicked the right button on a jagged graph. By trusting the automatic process, you are using a tool that has been refined through years of peer-reviewed research, ensuring that your data are as clean and accurate as possible before you draw any conclusions about the body's physiological state.
Although applications like Kubios have become increasingly refined, beginners and experienced professionals should always "trust, but verify" by examining the the resulting tachogram (beat-to-beat map) for artifacts.
Understanding Trends in HRV Data
Once the sharpest errors have been scrubbed away, the data are ready for the second stage: detrending. To understand detrending, one must first understand stationarity. In heart rate research, we prefer that our data be stationary, meaning their average value and their variance stay relatively constant over time.
During long measurement periods, there is often a gradual change in the mean or variance of IBI values, called a trend. Trends occur for various physiological and environmental reasons.
Human bodies are naturally non-stationary. If a person breathes deeply or shifts their weight, the heart rate may slowly drift upward or downward over several minutes. This slow drift is called a trend. If we left these trends in the data, they would appear as massive waves in our final analysis, drowning out the faster, more subtle heart rate changes we actually want to measure. Detrending is the technical process of identifying this slow, underlying drift and subtracting it from the data so that the remaining heartbeats sit on a level, horizontal line.

Trends can violate the stationarity assumption underlying many HRV analysis methods. Stationarity means that the statistical properties of the IBI values—such as the mean and standard deviation—do not change systematically over time. When this assumption is violated, frequency-domain metrics may not accurately reflect underlying autonomic function because the algorithm cannot distinguish between true oscillatory patterns and the drift superimposed upon them.
The graphics below illustrate acceptable versus weak stationarity (adapted from Kubios). In the acceptable case, statistical properties remain relatively stable throughout the recording. In the weak stationarity example, some drift is present but may still be acceptable depending on your clinical purposes and the specific metrics you need.
Detrending: Removing Systematic Changes
Detrending is the mathematical removal of a trend from a time series to ensure stationarity. Kubios provides a smoothness priors function for this purpose, allowing you to select the smoothness of the removed trend by entering a lambda value. Lambda is the smoothing factor; higher values remove smoother (lower-frequency) trends. Kubios displays the corresponding cutoff frequency beside the lambda value you enter, helping you understand exactly what you are removing.
Implementing Smoothness Priors in Practice
In Kubios, smoothness priors is a sophisticated digital filter that acts like a flexible ruler held against the data points. The software calculates a smooth, curving line that follows the general trend of the heart rate without being distracted by the quick, beat-to-beat changes.
This smooth line represents the noise of the body’s slow drifts. Once this line is calculated, the software performs a subtraction: it takes the original data and removes the values represented by that smooth line. What remains are detrended data where every heartbeat is measured against a stable baseline. This ensures that when we later look at the frequencies of the heart—such as those associated with breathing or stress—we are looking at the pure rhythm of the heart's internal clock rather than the slow physical movements of the person being tested.
Adjusting Filter Intensity
When setting up this process, a beginner will encounter a technical parameter known as lambda. Lambda is the setting that determines how stiff or flexible that digital ruler should be. If the lambda is set very high, the ruler is stiff and only removes the very slowest, most gradual drifts.
If the lambda is set low, the ruler becomes more flexible and starts to follow even the quicker changes in heart rate. For most standard research, a lambda of 500 is the gold standard because it is specifically tuned to ignore any movement that happens slower than about 0.035 Hertz, which is a technical way of saying it ignores anything that takes longer than 30 seconds to complete one cycle. By keeping this setting consistent, researchers ensure that they are only removing the non-biological drift while preserving the rapid fluctuations that reveal how the nervous system is responding to the world.
The following graphics demonstrate how different lambda values affect the FFT spectrum (Tarvainen & Niskanen, 2020). All lambda values markedly reduce VLF power (0-0.04 Hz) and may also lower LF power (0.04-0.15 Hz). Understanding this trade-off is essential: aggressive detrending can inadvertently remove genuine physiological variation along with the unwanted trend.




The upper section of the detrending graphic below shows the original time series with a smoothed trend (red line). The lower section shows the same RR series with the trend removed. Notice how the mean becomes stable after detrending—you can draw a straight horizontal line through the mean values, indicating that the systematic drift has been eliminated while preserving the beat-to-beat variability you need to analyze.

Executing the Workflow in Kubios Scientific
To begin this process in the Kubios Scientific interface, you first open your data file and look to the left-hand side of the screen where the analysis settings are housed. You will first click on the Artifact Correction section to address your outliers. In the dropdown menu, you can select a correction level, such as "Strong" or "Very Strong," which uses the software’s algorithms to compare each beat to its neighbors and automatically replace the outliers with estimated values.
Once these data are corrected, you move your mouse to the "Trend Removal" or "Detrending" section immediately below. Here, you will click the dropdown menu and select "Smoothness Priors." After making this selection, you will see a small box labeled "Lambda." You should ensure this value is set to 500 for a standard five-minute recording. Once these selections are made, the software will process the data, first purging the outliers and then smoothing the trends, leaving you with a clean, level signal ready for final interpretation.

Check Your Understanding: Plan A
- What is a time series in the context of HRV analysis?
- Why is stationarity important for HRV measurement?
- What happens if you select a lambda value that is too low during detrending?
- When is Plan A the appropriate choice for handling artifacts?
- What is the functional difference between noise detection and beat correction in Kubios Scientific?
- How does decreasing the deviation threshold affect how many beats Kubios flags as artifacts, and what is the risk of using a threshold that is too strict?
- Why is the automatic correction algorithm generally preferred over fixed-threshold deviation settings for most recordings?
A Practical Perspective
This section addresses the reality that clinical practice often differs from research protocols. You will learn how to balance thoroughness with efficiency, when simplified approaches are appropriate, and how to develop practical artifacting skills through mentorship.
Artifacting is as vital in HRV biofeedback as it is in neurofeedback. However, the intensive data analysis required for publication-quality research may exceed what routine clinical practice demands. The key distinction lies in your purpose: research requires demonstrable data integrity to withstand peer review, while clinical practice requires data quality sufficient to support valid treatment decisions. Both matter, but they involve different standards and workflows.
Clinicians can often limit their review to pre- and post-baseline recordings rather than scrutinizing every training session. You may not need to detrend 5-minute datasets if you can find 3 consecutive minutes of clean data—the Task Force (1996) established this as the minimum for valid short-term HRV assessment. If artifacts have corrupted all 5 minutes, examine automatic artifacting within your analysis software and reserve manual artifacting for cases where automated methods prove inadequate.
Artifacting skills can be readily mastered through guidance from an experienced mentor. Contact the Biofeedback Certification International Alliance (BCIA) or your equipment vendor for mentor recommendations. A skilled mentor can help you develop clinical judgment about when to invest additional time in careful artifacting versus when simpler approaches suffice.

Plan B: Manual Correction and Detrending
This section covers manual artifact correction, the gold standard when automated methods are insufficient. You will learn to use visualization tools for artifact identification and master specific correction strategies for different artifact types. Review the ectopic beats section above before you begin, because the repair you choose depends on whether the outlier came from a real ectopic beat or from a detection error.
Manual correction and detrending is the best approach when you cannot find a clean 5-minute or longer segment and fewer than 5% of all beats are corrupted (Berntson & Stowell, 1998). Pre-process your data using Excel or a dedicated HRV analysis program like Thought Technology's CardioPro®. This threshold exists because beyond 5% corruption, even careful manual correction may not produce data that accurately represents the client's true autonomic function.
Kubios' Time-Domain tab heart rate plot helps you visualize invalid beats (R. Gevirtz, personal communication, 2018). The plot displays heart rate over time, and artifacts typically appear as dramatic spikes or drops that deviate sharply from surrounding values. Because heart rate is the inverse of IBI (HR = 60,000/IBI in ms), a missed beat appears as a sudden heart rate drop, while an extra beat appears as a sudden spike.
Kubios' Non-Linear tab provides a Poincaré plot, a powerful tool for identifying false beats (R. Gevirtz, personal communication, 2018). In a Poincaré plot, each IBI is plotted against the subsequent IBI, creating a visual map of beat-to-beat variability (Sassi et al., 2015). Normal beats cluster together in an elliptical pattern oriented along a 45-degree line of identity. Artifacts appear as outliers—points that fall far from this main cluster, making them immediately visible even when they might be difficult to spot in a simple time-series display.
Correcting Extra Beats and Missed Beats
Understanding the two fundamental artifact patterns is essential for appropriate correction. Adjacent skeletal muscle activity (SEMG), a double-counted T wave, and similar noise sources produce false extra beats that appear as very short IBIs. True ectopic beats also shorten an interval, but they are genuine cardiac events and are handled as described in the ectopic beats section above. When the software incorrectly detects a non-cardiac event as a heartbeat, it generates two artificially short intervals where one normal interval should exist. Think of it as the software "seeing" a beat that isn't really there, splitting one true interval into two false ones.

Conversely, when the pulse wave (BVP) peak or R-spike (ECG) amplitude is too low for software recognition, the result is missed beats that appear as prolonged IBIs. The software fails to detect a genuine heartbeat, measuring from one detected beat to the next, effectively skipping the missed one. This commonly occurs with low-amplitude R-spikes, poor electrode contact, or excessive movement.
Manual Correction Strategies
Effective manual correction requires matching the correction strategy to the artifact type. The three primary operations—Add, Split, and Average—each address different patterns of data corruption.
In CardioPro's data editing window (shown below), the Add, Split, and Average options appear at the bottom left. Highlight the abnormal beats (circled in blue) and select the appropriate correction. The Add function combines two consecutive short IBIs into one normal IBI—use this for extra beats. The Split function divides one long IBI into two normal IBIs—use this for missed beats. The Average function replaces abnormal values with an average of surrounding normal values—use this for PVCs and other patterns where the total time should remain constant.

At the current state of the art, expert manual editing can be superior to automated editing. While automated algorithms continue to improve, experienced clinicians can make nuanced judgments about artifact correction that software cannot replicate. Manual editing allows you to consider context—was the client moving, coughing, or talking during that segment? Is this pattern consistent with PVCs versus a missed beat? These contextual factors inform optimal correction choices.
Check Your Understanding: Plan B
- What percentage of corrupted beats is the threshold for using Plan B rather than discarding the data?
- How does a Poincaré plot help identify artifacts?
- Describe the IBI pattern produced by a PVC with a compensatory pause.
- Why are PACs more difficult to correct than PVCs?
- When would you use the Split function versus the Add function in manual correction?
Check Your Understanding: A Practical Perspective
- How do the artifacting requirements for clinical practice differ from those for research publication?
- What is the minimum amount of clean data typically needed for HRV biofeedback assessment?
- When should a clinician resort to manual artifacting in routine practice?
Cutting Edge Topics
Machine Learning Approaches to Artifact Detection
Researchers are increasingly applying machine learning algorithms to automate and improve artifact detection in HRV data. These approaches can learn from large datasets of expert-labeled recordings to identify artifacts with greater sensitivity and specificity than traditional threshold-based methods (Cherif et al., 2016). Deep learning models show particular promise for distinguishing between genuine ectopic beats and recording artifacts—a clinically important distinction because ectopic beats may warrant medical referral while recording artifacts simply need correction.
Real-Time Artifact Correction in Biofeedback Applications
Emerging biofeedback systems incorporate real-time artifact detection and correction, allowing clients to receive immediate feedback even when minor recording issues occur. These systems must balance rapid processing against the accuracy requirements of clinical assessment. Some developers are exploring hybrid approaches that apply simple corrections in real-time while flagging segments for more careful post-session review—a practical compromise that maintains session flow without sacrificing data quality.
Wearable Device Artifacts: New Challenges
The proliferation of consumer-grade HRV monitoring devices has introduced new artifacting challenges (Altini, 2016). Wrist-based optical sensors using photoplethysmography (PPG) are particularly susceptible to motion artifacts during daily activities—the very situations where ambulatory HRV monitoring would be most valuable (Elgendi, 2012). Researchers are developing specialized algorithms to address these unique artifact patterns, enabling more accurate HRV assessment outside controlled laboratory environments. Understanding these limitations is essential for clinicians who incorporate wearable device data into their assessments.
Assignment
Now that you have completed this module, review the missed beat and extra beat graphics and see whether you can identify abnormally long and short IBIs. Based on this module, how might you improve your artifacting practice? Consider the following questions as you reflect on your learning.
First, think about your current recording setup. What steps could you take to prevent artifacts before they occur? Consider factors like sensor placement, environmental interference from nearby electronics, and client preparation including skin cleaning and minimizing movement.
Second, evaluate your current artifacting workflow. Are you using the most appropriate plan (A or B) for your typical recording scenarios? What conditions would prompt you to switch between plans?
Third, consider how you would explain the importance of data quality to a client who is impatient with the careful setup procedures required for clean recordings. How can you help them understand that taking an extra minute for proper preparation ultimately produces more valid and useful results?
Acknowledgment
This article draws heavily on graphics published in Didier Combatalade's Basics of Heart Rate Variability Applied to Psychophysiology (Thought Technology Ltd.) and Tarvainen and Niskanen's Kubios HRV User's Guide. A special thanks to Dr. Richard Gevirtz for sharing Alliant University's artifacting protocol. Dr. Gevirtz is a gifted clinician, educator, and researcher whose contributions have shaped our understanding of HRV biofeedback.

Glossary
aberrant conduction: supraventricular activation that reaches part of the ventricular conduction system while it is still refractory, producing a widened or altered QRS complex.
artifact: an unwanted signal component introduced by movement, environmental interference, hardware, processing, or biological activity outside the signal of interest.
artifacting: the process of identifying and excluding or correcting artifact-contaminated portions of a physiological record.
bigeminy: a repeating rhythm in which every other beat is ectopic.
blocked PAC: a premature atrial depolarization that fails to conduct through the atrioventricular node and therefore produces no following QRS complex.
compensatory pause: a post-ectopic pause in which the premature interval plus the following interval approximate two baseline sinus cycles, usually because the sinus node was not reset.
couplet: two consecutive ectopic beats with no intervening sinus beat.
cubic-spline interpolation: estimation of missing or corrected data with piecewise cubic polynomials joined smoothly at selected knots.
detrending: removal of a modeled slow trend from a time series; detrending can improve but does not guarantee stationarity.
deviation threshold: a software parameter specifying how far an observed beat interval may deviate from a local estimate before being flagged; values are algorithm-specific.
ectopic beat: a cardiac beat initiated outside the sinoatrial node, such as a premature atrial or ventricular beat; it is excluded from normal-to-normal interval series.
ectopic burden: the count or percentage of ectopic beats within a recording, reported separately from HRV metrics.
epoch: a defined time segment of a physiological signal selected for recording, averaging, artifact review, or analysis.
fusion beat: a ventricular complex produced by near-simultaneous activation from a sinus impulse and a ventricular ectopic impulse.
heart rate turbulence (HRT): the characteristic sinus acceleration and deceleration following an isolated PVC that reflects baroreflex activity.
heart rate variability (HRV): variation in successive cardiac intervals, usually quantified from normal sinus beats and shaped by autonomic, respiratory, circadian, behavioral, and other influences.
interbeat interval (IBI): the elapsed time between successive heartbeats, usually measured between adjacent R peaks in ECG and expressed in ms.
interpolated PVC: a PVC inserted between two conducted sinus beats without the usual compensatory pause.
interpolation: estimation of missing or rejected interval values from neighboring observations or from a timing model.
lambda wave: a positive occipital transient occurring during visual scanning and time-locked to saccadic eye movements.
missed-beat artifact: failure of a detection algorithm to identify a true heartbeat, producing an interbeat interval approximately equal to the sum of adjacent intervals.
noncompensatory pause: a post-PAC pause for which the interval spanning the premature beat is shorter than two expected sinus cycles, usually because the PAC reset the sinus node.
normal-to-normal interval (NN interval): the interval between successive normal sinus beats after ectopic beats, artifacts, and non-sinus intervals have been excluded or corrected.
outlier: an observation unusually distant from the main data pattern under a specified model or distribution.
P wave: the ECG deflection representing atrial depolarization.
percentage of successive NN interval differences exceeding 50 ms (pNN50): the percentage of consecutive accepted NN interval pairs that differ by more than 50 ms.
Poincaré plot: a scatterplot of each NN interval against the next, used to visualize beat-to-beat variability, patterns, and outliers.
premature atrial contraction (PAC): an early atrial depolarization arising outside the sinoatrial node, usually producing an abnormal P wave and a noncompensatory pause.
premature ventricular contraction (PVC): an early ventricular depolarization arising outside the normal conduction system, producing a wide QRS complex and often a compensatory pause.
QRS complex: the ECG waveform produced primarily by ventricular depolarization.
R wave: the first positive deflection of the QRS complex; some QRS complexes contain no R wave.
root mean square of successive differences (RMSSD): the square root of the mean squared differences between consecutive NN intervals, a short-term time-domain HRV measure.
RR interval: the time between consecutive detected R spikes, whether the associated beats are normal, ectopic, or falsely detected.
sinoatrial (SA) node: the usual cardiac pacemaker that initiates sinus rhythm.
sinus rhythm: a cardiac rhythm initiated by the SA node with expected atrial and ventricular conduction.
smoothness-priors method: a detrending method that estimates a smooth trend by penalized least squares, with the regularization parameter controlling the permitted curvature of the trend.
standard deviation of NN intervals (SDNN): the standard deviation of all NN intervals within a specified recording period, expressed in ms; interpretation depends strongly on recording duration.
stationarity: the property that specified statistical characteristics of a time series remain constant over time.
tachogram: a graph of successive RR or NN intervals, or their corresponding heart rates, plotted over time.
time series: an ordered sequence of observations indexed by time.
time-varying artifact detection: an adaptive method that updates the expected RR interval from neighboring beats as heart rate changes over time.
trend: a systematic change in the level or another statistical property of a time series across time.
trigeminy: a repeating rhythm in which an ectopic beat follows every two sinus beats.
triplet: three consecutive ectopic beats; longer runs are variously classified as runs of PVCs or ventricular tachycardia, so state the convention you apply.
waveform morphology: the structural features of a waveform, including shape, duration, amplitude, slope, sharpness, and polarity.
References
Altini, M. (2016). Issues in heart rate variability (HRV) analysis: Motion artifacts & ectopic beats. HRV4Training. https://www.hrv4training.com/blog/issues-in-heart-rate-variability-hrv-analysis-motion-artifacts-ectopic-beats
Bassi-Dibai, D., Santos-de-Araújo, A. D., Rocha, D. S., de Almeida, L. V., Figueiredo, J. K., Gondim, L. A. R., Diniz, M. R. de F., Frutuoso, V. P., Maia, M. C., Santos, P. M., & Borghi-Silva, A. (2026). Threshold-based artefact correction methods influence heart rate variability measurements in individuals with type 2 diabetes mellitus. Scientific Reports, 16, Article 11341. https://doi.org/10.1038/s41598-026-42255-y
Bauer, A., Malik, M., Schmidt, G., Barthel, P., Bonnemeier, H., Cygankiewicz, I., Guzik, P., Lombardi, F., Müller, A., Oto, A., Schneider, R., Watanabe, M., Wichterle, D., & Zareba, W. (2008). Heart rate turbulence: Standards of measurement, physiological interpretation, and clinical use. Journal of the American College of Cardiology, 52(17), 1353–1365. https://doi.org/10.1016/j.jacc.2008.07.041
Berntson, G. G., Bigger, J. T., Jr., Eckberg, D. L., Grossman, P., Kaufmann, P. G., Malik, M., Nagaraja, H. N., Porges, S. W., Saul, J. P., Stone, P. H., & van der Molen, M. W. (1997). Heart rate variability: Origins, methods, and interpretive caveats. Psychophysiology, 34(6), 623–648. https://doi.org/10.1111/j.1469-8986.1997.tb02140.x
Berntson, G. G., & Stowell, J. R. (1998). ECG artifacts and heart period variability: Don't miss a beat! Psychophysiology, 35, 127–132.
Bourdillon, N., Yazdani, S., Vesin, J.-M., Schmitt, L., & Millet, G. P. (2022). RMSSD is more sensitive to artifacts than frequency-domain parameters: Implication in athletes' monitoring. Journal of Sports Science and Medicine, 21(2), 260–266. https://doi.org/10.52082/jssm.2022.260
Cherif, S., Pastor, D., Nguyen, Q.-T., & L'Her, E. (2016). Detection of artifacts on photoplethysmography signals using random distortion testing. In 2016 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) (pp. 6214–6217). IEEE. https://doi.org/10.1109/EMBC.2016.7592148
Choi, A., & Shin, H. (2018). Quantitative analysis of the effect of an ectopic beat on the heart rate variability in the resting condition. Frontiers in Physiology, 9, Article 922. https://doi.org/10.3389/fphys.2018.00922
Clifford, G. D., & Tarassenko, L. (2005). Quantifying errors in spectral estimates of HRV due to beat replacement and resampling. IEEE Transactions on Biomedical Engineering, 52(4), 630–638. https://doi.org/10.1109/TBME.2005.844028
Clinical EEG Interpretation. (2018). Premature ventricular contractions (premature ventricular complex, premature ventricular beats). https://ecgwaves.com/premature-ventricular-contractions-complex-beats-ecg/
Combatalade, D. (2010). Basics of heart rate variability applied to psychophysiology. Thought Technology Ltd.
Combatalade, D. (2013). CardioPro Infiniti: HRV analysis module for BioGraph Infiniti. Thought Technology Ltd.
Elgendi, M. (2012). On the analysis of fingertip photoplethysmogram signals. Current Cardiology Reviews, 8(1), 14–25. https://doi.org/10.2174/157340312801215782
Friesen, G. M., Jannett, T. C., Jadallah, M. A., Yates, S. L., Quint, S. R., & Nagle, H. T. (1990). A comparison of the noise sensitivity of nine QRS detection algorithms. IEEE Transactions on Biomedical Engineering, 37(1), 85–98. https://doi.org/10.1109/10.43620
Guichard, J.-B., Guasch, E., Roche, F., Da Costa, A., & Mont, L. (2022). Premature atrial contractions: A predictor of atrial fibrillation and a relevant marker of atrial cardiomyopathy. Frontiers in Physiology, 13, Article 971691. https://doi.org/10.3389/fphys.2022.971691
Kligfield, P., Gettes, L. S., Bailey, J. J., Childers, R., Deal, B. J., Hancock, E. W., van Herpen, G., Kors, J. A., Macfarlane, P., Mirvis, D. M., Pahlm, O., Rautaharju, P., & Wagner, G. S. (2007). Recommendations for the standardization and interpretation of the electrocardiogram: Part I: The electrocardiogram and its technology. Circulation, 115(10), 1306–1324. https://doi.org/10.1161/CIRCULATIONAHA.106.180200
Kuusela, T. (2013). Methodological aspects of heart rate variability analysis. In M. V. Kamath, M. A. Watanabe, & A. R. M. Upton (Eds.), Heart rate variability (HRV) signal analysis (pp. 9–42). CRC Press.
Laborde, S., Mosley, E., & Thayer, J. F. (2017). Heart rate variability and cardiac vagal tone in psychophysiological research: Recommendations for experiment planning, data analysis, and data reporting. Frontiers in Psychology, 8, Article 213. https://doi.org/10.3389/fpsyg.2017.00213
Langner, P. H., Jr., & Geselowitz, D. B. (1960). Characteristics of the frequency spectrum in the normal electrocardiogram and in subjects following myocardial infarction. Circulation Research, 8(3), 577–584. https://doi.org/10.1161/01.RES.8.3.577
Lin, I.-M., & Peper, E. (2009). Keep cell phones and PDAs away from EMG sensors and the human body to prevent electromagnetic interference artifacts and cancer. Biofeedback, 37(3), 114–116. https://doi.org/10.5298/1081-5937-37.3.114
Lippman, N., Stein, K. M., & Lerman, B. B. (1994). Comparison of methods for removal of ectopy in measurement of heart rate variability. American Journal of Physiology-Heart and Circulatory Physiology, 267(1), H411–H418. https://doi.org/10.1152/ajpheart.1994.267.1.H411
Lipponen, J. A., & Tarvainen, M. P. (2019). A robust algorithm for heart rate variability time series artefact correction using novel beat classification. Journal of Medical Engineering & Technology, 43(3), 173–181. https://doi.org/10.1080/03091902.2019.1640306
Marcus, G. M. (2020). Evaluation and management of premature ventricular complexes. Circulation, 141(17), 1404–1418. https://doi.org/10.1161/CIRCULATIONAHA.119.042434
Mateo, J., & Laguna, P. (2003). Analysis of heart rate variability in the presence of ectopic beats using the heart timing signal. IEEE Transactions on Biomedical Engineering, 50(3), 334–343. https://doi.org/10.1109/TBME.2003.808831
Mond, H. G., & Haqqani, H. M. (2019). The electrocardiographic footprints of atrial ectopy. Heart, Lung and Circulation, 28(10), 1463–1471. https://doi.org/10.1016/j.hlc.2019.03.005
Nederend, I., Jongbloed, M. R. M., de Geus, E. J. C., Blom, N. A., & ten Harkel, A. D. J. (2016). Postnatal cardiac autonomic nervous control in pediatric congenital heart disease. Journal of Cardiovascular Development and Disease, 3(2), Article 16. https://doi.org/10.3390/jcdd3020016
Peltola, M. A. (2012). Role of editing of R-R intervals in the analysis of heart rate variability. Frontiers in Physiology, 3, Article 148. https://doi.org/10.3389/fphys.2012.00148
Peper, E., Shaffer, F., & Lin, I.-M. (2010). Garbage in; garbage out: Identify blood volume pulse (BVP) artifacts before analyzing and interpreting BVP, blood volume pulse amplitude, and heart rate/respiratory sinus arrhythmia data. Biofeedback, 38(1), 19–23. https://doi.org/10.5298/1081-5937-38.1.19
Rincon Soler, A. I., Silva, L. E. V., Fazan, R., Jr., & Murta, L. O., Jr. (2018). The impact of artifact correction methods of RR series on heart rate variability parameters. Journal of Applied Physiology, 124(3), 646–652. https://doi.org/10.1152/japplphysiol.00927.2016
Salo, M. A., Huikuri, H. V., & Seppänen, T. (2001). Ectopic beats in heart rate variability analysis: Effects of editing on time and frequency domain measures. Annals of Noninvasive Electrocardiology, 6(1), 5–17. https://doi.org/10.1111/j.1542-474X.2001.tb00080.x
Sassi, R., Cerutti, S., Lombardi, F., Malik, M., Huikuri, H. V., Peng, C.-K., Schmidt, G., & Yamamoto, Y. (2015). Advances in heart rate variability signal analysis: Joint position statement by the e-Cardiology ESC Working Group and the European Heart Rhythm Association co-endorsed by the Asia Pacific Heart Rhythm Society. Europace, 17(9), 1341–1353. https://doi.org/10.1093/europace/euv015
Sattar, Y., & Hashmi, M. F. (2025). Premature ventricular complex. In StatPearls. StatPearls Publishing. https://www.ncbi.nlm.nih.gov/books/NBK547713/
Shaffer, F., & Combatalade, D. C. (2013). Don't add or miss a beat: A guide to cleaner heart rate variability recordings. Biofeedback, 41(3), 121–130. https://doi.org/10.5298/1081-5937-41.3.04
Shaffer, F., & Ginsberg, J. P. (2017). An overview of heart rate variability metrics and norms. Frontiers in Public Health, 5, Article 258. https://doi.org/10.3389/fpubh.2017.00258
Stapelberg, N. J. C., Neumann, D. L., Shum, D. H. K., McConnell, H., & Hamilton-Craig, I. (2018). The sensitivity of 38 heart rate variability measures to the addition of artifact in human and artificial 24-hour cardiac recordings. Annals of Noninvasive Electrocardiology, 23(1), Article e12483. https://doi.org/10.1111/anec.12483
Stern, R. M., Ray, W. J., & Quigley, K. S. (2001). Psychophysiological recording (2nd ed.). Oxford University Press.
Tarkiainen, T. H., Kuusela, T. A., Tahvanainen, K. U. O., Hartikainen, J. E. K., Tiittanen, P., Timonen, K. L., & Vanninen, E. J. (2007). Comparison of methods for editing of ectopic beats in measurements of short-term non-linear heart rate dynamics. Clinical Physiology and Functional Imaging, 27(2), 126–133. https://doi.org/10.1111/j.1475-097X.2007.00726.x
Tarvainen, M. P., & Niskanen, J.-P. (2020). Kubios HRV version 3.2 user's guide. University of Eastern Finland.
Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. (1996). Heart rate variability: Standards of measurement, physiological interpretation, and clinical use. Circulation, 93(5), 1043–1065. https://doi.org/10.1161/01.CIR.93.5.1043
Weinstein, A., Rodiño, J., & Otero, M. (2026). Effects of electrocardiogram QRS detection algorithms in heart rate variability metrics. Scientific Reports, 16, Article 18589. https://doi.org/10.1038/s41598-026-49215-6
Zeppenfeld, K., Tfelt-Hansen, J., de Riva, M., Winkel, B. G., Behr, E. R., Blom, N. A., Charron, P., Corrado, D., Dagres, N., de Chillou, C., Eckardt, L., Friede, T., Haugaa, K. H., Hocini, M., Lambiase, P. D., Marijon, E., Merino, J. L., Peichl, P., Priori, S. G., . . . ESC Scientific Document Group. (2022). 2022 ESC guidelines for the management of patients with ventricular arrhythmias and the prevention of sudden cardiac death. European Heart Journal, 43(40), 3997–4126. https://doi.org/10.1093/eurheartj/ehac262
Zhao, L., Li, P., Li, J., & Liu, C. (2021). Influence of ectopic beats on heart rate variability analysis. Entropy, 23(6), Article 648. https://doi.org/10.3390/e23060648
Return to Top