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. 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).

Artifacts in HRV data

Artifacts represent false values that do not reflect actual cardiac activity and can significantly distort HRV measurements.

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.

Effect of multiple artifacts in a recording segment

Multiple artifacts in a recording segment can produce dramatically distorted HRV metrics, emphasizing the importance of data quality.

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.

Dilbert cartoon about data quality

The importance of data quality cannot be overstated. © Scott Adams.

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.

Prevention is always preferable to correction. A closed barn door represents proper sensor placement and recording conditions, while the escaping horse represents artifacts that contaminate your data after the fact.

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.

Kubios HRV review screen

Kubios Scientific provides comprehensive tools for reviewing raw signals and identifying artifacts. The review screen allows you to examine both the ECG waveform and IBI values simultaneously.

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.

ECG signal with R-spike detection

Kubios Scientific software detects valid R-spikes in the raw ECG signal, designated by a red plus sign, and calculates the time interval between consecutive peaks.

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.

CardioPro IBI table

IBI values are stored in a time-stamped table, allowing you to identify specific beats and their corresponding intervals.

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.

Short IBI caused by extra beat

An extra beat after 00:00:32 produces an abnormally short IBI value (274.9 ms) that stands out from typical values.

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.

Long IBI caused by missed beat

A missed beat before 00:00:31 results in an abnormally long IBI value (1334.4 ms) because the software measures from one detected beat to the next detected beat, skipping the undetected beat. This exceeds the normal adult resting range (600-1200 ms) and the previous IBI by more than 20%.

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.

Typical IBI values

Typical resting IBI values fall within a normal adult range of 600-1200 ms and don't change more than 20% across successive beats. This makes artifactual values easy to identify.

Artifacts are false IBI values produced by ectopic beats (initiated outside of the S-A node), movement, environmental interference, or hardware limitations. A single artifact can dramatically distort HRV metrics like SDNN because these measures are sensitive to outliers. Prevention through proper recording technique is always preferable to correction. Understanding the connection between raw signal distortions and computed IBI values is essential for identifying artifacts effectively.

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.

Removing artifacts concept

Artifacting involves identifying and correcting false values to ensure data integrity.

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.

Person adjusting to monitoring

Clients typically need time to adjust to the monitoring environment, and the first minute of data often contains more artifacts. Graphic © Vitaly Titov/Shutterstock.com.

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.

don't splice

The splicing approach may seem intuitive, but it destroys the temporal integrity of your data.

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.

Consider a clinician who notices a 30-second period of movement artifact in the middle of a 5-minute recording. If they simply cut out those 30 seconds and rejoin the remaining data, the resulting time series would show artificial transitions where the splice occurred. Frequency-domain analysis would be particularly affected because the algorithm expects a continuous time series with natural beat-to-beat relationships.
Artifacting is the systematic process of identifying and correcting false IBI values to ensure data integrity. The first and last minutes of recordings often require special attention due to client adjustment and fatigue. The splicing method of removing contaminated sections and rejoining remaining data is never appropriate because it destroys the temporal relationships essential for valid HRV analysis, particularly frequency-domain measures.

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 (Gevirtz, 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.

IBI column data

A time series of IBI values appears as a column of consecutive millisecond values.

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.

Editing these noise markings is designed to be quick and intuitive using a few simple mouse shortcuts. If the software has flagged a section that you know is actually clean, you can simply right-click that marking to remove it instantly. To mark a new section as noise—or to unmark an existing one—you can use a combination of your keyboard and mouse. Windows users simply hold the Control key while left-clicking the data area, while Mac users do the same using the Command key. These tools allow you to surgically clean your recording without needing to navigate complex menus or difficult settings.

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

There are two segments, colored green and blue, separated by corrupted IBIs. Kubios allows users to choose either clean segment.

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.

Trend in IBI data

A trend represents a gradual, systematic change in IBI values over time that can violate stationarity assumptions.

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.

Acceptable stationarity example

Acceptable stationarity: The statistical properties of the IBI series remain relatively stable throughout the recording period.

Weak stationarity example

Weak stationarity: Some drift is present in the data, but it may still be acceptable for analysis in clinical or performance applications.

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.

No smoothing applied

No Smoothing: The original data without any detrending applied.

Lambda 250 detrending

Lambda 250 with Cutoff = 0.051 Hz: A moderate amount of detrending is applied.

Lambda 500 detrending

Lambda 500 with Cutoff = 0.035 Hz: More aggressive detrending that may begin to affect LF power.

Lambda 1000 detrending

Lambda 1000 with Cutoff = 0.029 Hz: Very aggressive detrending that will significantly reduce VLF and may affect LF power.

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.

Kubios detrending visualization

Kubios detrending visualization shows the original series with trend (top) and the detrended series (bottom).

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.

You reach Kubios' Preferences menu from its File menu.

Plan A is appropriate when you have at least 6 minutes of data and artifacts are confined to the first or last 3 minutes. Within Kubios Scientific, artifact cleaning proceeds in two stages: noise detection removes large contaminated segments, and beat correction addresses individual outlier intervals. Kubios identifies abnormal beats using deviation thresholds—the maximum allowable difference between an observed and expected RR interval—with settings ranging from Very Low (±450 ms) to Very Strong (±50 ms). The automatic correction algorithm is recommended for most recordings because it adapts to the individual's cardiac rhythm using a time-varying approach and applies cubic spline interpolation to replace abnormal intervals. A time series is a consecutive sequence of IBI values. Trends are gradual changes in statistical properties that can violate the stationarity assumption required by many HRV analysis methods. Detrending mathematically removes these trends using the smoothness priors method; the lambda parameter controls filter intensity, with higher values removing smoother (lower-frequency) trends. A standard lambda of 500 is recommended for most 5-minute recordings. Exercise caution to avoid removing meaningful physiological variation along with the unwanted drift.

Check Your Understanding: Plan A

  1. What is a time series in the context of HRV analysis?
  2. Why is stationarity important for HRV measurement?
  3. What happens if you select a lambda value that is too low during detrending?
  4. When is Plan A the appropriate choice for handling artifacts?
  5. What is the functional difference between noise detection and beat correction in Kubios Scientific?
  6. 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?
  7. 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 of the European Society of Cardiology (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.

Mentorship concept

Working with an experienced mentor can accelerate your mastery of artifacting skills and help you develop clinical judgment about when and how to apply different correction strategies.
Rene, a biofeedback clinician in private practice, was initially overwhelmed by the complexity of research-level artifacting procedures. Her mentor helped her develop a practical workflow: first, scan the recording visually for obvious problems; second, check if at least 3 clean minutes exist; third, use automatic correction if needed and verify the results make physiological sense. This streamlined approach allowed Rene to maintain data quality without spending excessive time on every session recording.

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, understand the physiology of ectopic beats, and master specific correction strategies for different artifact types.

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.

Flagging abnormal IBIs

The first step in manual correction is identifying and flagging abnormal IBI values for review and correction.

Kubios' Time-Domain tab heart rate plot helps you visualize invalid beats (Gevirtz, 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 Time-Domain tab

The Kubios Time-Domain tab provides a heart rate plot that makes artifacts visually apparent as sharp deviations from the baseline pattern.

Kubios' Non-Linear tab provides a Poincaré plot, a powerful tool for identifying false beats (Gevirtz, 2018). In a Poincaré plot, each IBI is plotted against the subsequent IBI, creating a visual map of beat-to-beat variability. 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.

Kubios Non-Linear tab with Poincare plot

The Poincaré plot in Kubios' Non-Linear tab displays each IBI against the subsequent IBI. Artifacts appear as points that fall far from the main elliptical cluster of normal beats.

Correcting Extra Beats and Missed Beats

Understanding the two fundamental artifact patterns is essential for appropriate correction. Adjacent skeletal muscle activity (SEMG) or ectopic beats produce extra beats that appear as very short IBIs. 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.

Extra beats shortening IBIs

Extra beats shorten IBIs. The software detects an additional "beat" that is not a true cardiac event, resulting in two abnormally short intervals. Graphic © Kubios Oy.

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.

Missed beats lengthening IBIs

Missed beats lengthen IBIs. When the software fails to detect a beat, the resulting interval spans two or more actual cardiac cycles. Graphic © Kubios Oy.

Premature Ventricular Contractions and Their Consequences

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.

Premature Ventricular Contractions (PVCs). In the lower ECG panel, the ectopic beats occur prematurely (earlier than the next expected sinus complex) and are characterized by a wide, abnormal (“bizarre”) QRS complex consistent with ventricular origin (typically ≥ 120 ms), rather than normal conduction through the His–Purkinje system. The PVCs also show no clearly preceding sinus P wave and are followed by secondary ST–T changes that are discordant (the T-wave deflects opposite the main QRS polarity), a common repolarization pattern after ventricular ectopy. Finally, the rhythm demonstrates a post-PVC compensatory pause before the next sinus beat, reflecting interruption of the expected sinus timing and reinforcing the diagnosis of a classic PVC

Premature atrial complexes (PACs) are initiated by an ectopic pacemaker in the atria rather than the ventricles. Since PACs involve an extra beat without a compensatory pause, there is no easy way to remove them through manual artifacting (Clinical EEG Interpretation, 2018). PACs are more challenging to correct because they don't produce the symmetrical long-short pattern that makes PVCs relatively straightforward to address.

The lower panel illustrates Premature Atrial Contractions (PACs), characterized by a regular underlying sinus rhythm interrupted by ectopic beats originating from the atria. Some beats arrive early relative to the baseline rhythm (the R–R interval suddenly shortens), and each early beat is preceded by an ectopic P wave that looks different from the surrounding sinus P waves and may partially distort the preceding T wave, indicating the impulse started in the atria but outside the sinoatrial node. The QRS complexes that follow these early atrial impulses remain narrow and “normal-looking” (ventricular activation still uses the His–Purkinje system), which helps distinguish PACs from premature ventricular beats that typically widen the QRS. After each premature atrial beat, the tracing shows a noncompensatory pause (the pacemaker timing is “reset,” so the pause is not a full two-cycle gap), producing the characteristic irregular spacing that, when frequent, raises the apparent rate (here labeled 115 bpm). .

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.

Manual correction involves choosing the appropriate strategy based on the type of artifact present.

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.

CardioPro's data editing window provides Add, Split, and Average options for manual artifact correction. Highlight abnormal beats and select the appropriate correction strategy.

View and detrend workflow

Manually edit corrupted beats, view your data to verify the corrections, and then apply detrending if necessary.

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.

Imagine you are reviewing Sarah's HRV recording and notice a short IBI of 320 ms followed by a long IBI of 1180 ms. This classic pattern suggests a PVC with a compensatory pause. The total time for both intervals is approximately 1500 ms, which would represent about two normal cardiac cycles at her typical heart rate of 80 bpm (750 ms per beat). Using the Average function, you would replace these two values with two intervals of approximately 750 ms each, preserving the temporal structure while eliminating the artifact.
Plan B is appropriate when artifacts are scattered throughout the recording but affect fewer than 5% of beats. Visualization tools like heart rate plots and Poincaré plots help identify artifacts efficiently. Extra beats cause abnormally short IBIs; missed beats cause abnormally long IBIs. PVCs produce a characteristic short-long pattern with a compensatory pause, while PACs lack this pause and are more difficult to correct. Manual correction strategies include Add (combine two short IBIs), Split (divide one long IBI), and Average (replace with surrounding values). Manual editing remains the gold standard for most clinical applications.

Check Your Understanding: Plan B

  1. What percentage of corrupted beats is the threshold for using Plan B rather than discarding the data?
  2. How does a Poincaré plot help identify artifacts?
  3. Describe the IBI pattern produced by a PVC with a compensatory pause.
  4. Why are PACs more difficult to correct than PVCs?
  5. When would you use the Split function versus the Add function in manual correction?
Clinical practice typically requires less intensive artifacting than research publication. The minimum standard for valid short-term HRV assessment is 3 consecutive minutes of clean data. For routine clinical work, use Kubios' automatic correction as a first-line approach and always verify results by visually inspecting the tachogram. Reserve manual artifacting for cases where automated methods prove inadequate. Seek mentorship to develop practical artifacting skills and clinical judgment. The goal is data quality sufficient for valid clinical decisions, not perfection.

Check Your Understanding: A Practical Perspective

  1. How do the artifacting requirements for clinical practice differ from those for research publication?
  2. What is the minimum amount of clean data typically needed for HRV biofeedback assessment?
  3. 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.

Dr. Richard Gevirtz has made fundamental contributions to the science and practice of HRV biofeedback, including the development of practical artifacting protocols.

Glossary

artifacts: measurement errors in calculating the interbeat interval, produced by ectopic beats, movement, environmental interference, or hardware limitations.

artifacting: the process of removing or replacing false values found within a data record to ensure data integrity for HRV analysis.

cubic spline interpolation: a mathematical technique used to estimate corrected values in an RR-interval time series by generating a smooth curve between surrounding data points, preserving the natural continuity of cardiac rhythm; Kubios applies this method when replacing intervals identified as artifacts during automatic beat correction.

detrending: the mathematical removal of a trend from a time series to ensure stationarity, commonly performed using smoothness priors methods.

deviation threshold: the maximum allowable difference between an observed RR interval and the expected interval based on surrounding beats; in Kubios HRV Scientific, this parameter determines how aggressively the software flags and corrects abnormal beats, with settings ranging from Very Low (±450 ms) to Very Strong (±50 ms).

ectopic beats: beats initiated outside of the S-A node that produce consecutive IBI values that are abnormally long and short, including PVCs and PACs.

epoch: a data collection period, typically 5 minutes for standard HRV assessment.

extra beats: ECG artifacts that shorten the IBI when signal distortion causes the software to detect nonexistent beats.

heart rate variability (HRV): the beat-to-beat changes in heart rate that include changes in the RR intervals between consecutive heartbeats.

interbeat interval (IBI): the time interval between the peaks of successive R-spikes (initial upward deflections in the QRS complex). The IBI is also called the NN (normal-to-normal) interval.

lambda: the parameter that specifies the smoothness of the removed trend in the smoothness priors method; higher values remove smoother trends.

missed beats: an ECG artifact that lengthens the IBI when signal distortion causes the software to overlook a beat and use the next good beat.

NN (normal-to-normal) interval: the time interval between successive normal sinus beats after artifacts have been removed; used interchangeably with interbeat interval (IBI) in clean recordings. See interbeat interval (IBI).

outliers: IBI values that fall far outside the expected range of normal cardiac intervals; in artifacting, extreme outliers that are physiologically implausible (e.g., below 400 ms or above 1,300 ms) are the primary targets of the initial stage of data cleaning.

Poincaré plot: a nonlinear analysis method that plots each IBI against the subsequent IBI, revealing patterns of variability and helping identify artifacts as outliers from the main cluster.

premature atrial complexes (PACs): ectopic beats originating within a pacemaker in the atria, producing an extra beat without a compensatory pause.

premature ventricular contractions (PVCs): ectopic beats originating from the ventricles, typically producing a short IBI followed by a compensatory pause.

R-spikes: the tall, sharp peaks in the ECG waveform representing ventricular depolarization; HRV software detects R-spikes to measure the time intervals between them, producing the interbeat interval (IBI) time series.

SDNN: the standard deviation of normal-to-normal (NN) intervals; a time-domain HRV metric reflecting overall autonomic regulation that is highly sensitive to outliers, making thorough artifact removal essential before computation.

smoothness priors method: a mathematical method for removing trends of a specific smoothness from a time series, controlled by the lambda parameter.

stationarity: the assumption that the statistical properties of a series of values (e.g., mean and standard deviation) do not change over time; a requirement for many HRV analysis methods.

tachogram: a graphical or tabular representation of the RR-interval time series showing how the time between successive heartbeats fluctuates over a recording period; used in Kubios and other HRV software to visually inspect data quality and verify the results of artifact correction.

time series: a set of values (IBIs) measured over time, forming the basis for HRV analysis.

time-varying approach: an artifact detection method that continuously recalculates the expected RR interval based on the heartbeats immediately surrounding each data point, adapting to the individual's unique cardiac rhythm rather than applying a fixed threshold; Kubios uses this approach in its automatic correction algorithm.

trend: a gradual change in a statistical property (mean or standard deviation) of a time series that can violate stationarity assumptions.

waveform morphology: the shape characteristics of individual ECG waveform components; although some artifact-detection methods analyze waveform shape to identify ectopic beats, Kubios HRV Scientific identifies abnormal beats by comparing interval lengths rather than examining waveform morphology.

References

Altini, M. (2016). Issues in heart rate variability (HRV) analysis: Motion artifacts & ectopic beats. HRV4Training Blog. https://www.hrv4training.com/blog/issues-in-heart-rate-variability-hrv-analysis-motion-artifacts-ectopic-beats

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, 623-648.

Berntson, G. G., & Stowell, J. R. (1998). ECG artifacts and heart period variability: Don't miss a beat! Psychophysiology, 35, 127-132. PMID: 9499713

Cherif, S., Pastor, D., Nguyen, Q.-T., & L'Her, E. (2016). Detection of artifacts on photoplethysmography signals using random distortion testing. 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). https://doi.org/10.1109/EMBC.2016.7592148

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, 14-25. https://doi.org/10.2174/157340312801215782

Gevirtz, R. (2018). Personal communication about artifacting.

Kuusela, K. (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. CRC Press.

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, 577-584.

Lehrer, P. (2018a). Heart rate, interbeat interval, and the electrocardiogram (shared presentation).

Lehrer, P. (2018b). Personal communication about artifacting.

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

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), 16. https://doi.org/10.3390/jcdd3020016

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

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.

Stern, R. M., Ray, W. J., & Quigley, K. S. (2001). Psychophysiological recording (2nd ed.). Oxford University Press.

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, 1043-1065. PMID: 8598068

Return to Top