<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Using Visual Analytics of Heart Rate Variation to Aid in Diagnostics</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Stephen McIntyre</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>J. Mikael Eklund</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Christopher Collins</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>University of Ontario Institute of Technology</institution>
        </aff>
      </contrib-group>
      <fpage>20</fpage>
      <lpage>27</lpage>
      <abstract>
        <p>We present an interactive visualization tool for exploring RR interval data (the time between consecutive heart beats) to support diagnostics. An RR interval sequence diagram allows us to reduce hours of data into a general overview opposed to using short term ECG strips. A simple moving average is applied to the sequence diagram to smooth out short-term variance and highlight long-term trends. The moving average is surrounded by standard deviation bands which allow us to see the lfuctuations in variance. After a brief training period using these tools coupled with RR interval and RR interval diference histograms, nonexpert participants (undergraduate students) were able to diferentiate between normal, atrial fibrillation, and congestive heart failure.</p>
      </abstract>
      <kwd-group>
        <kwd>HRV</kwd>
        <kwd>Time domain</kwd>
        <kwd>RR intervals</kwd>
        <kwd>Information visualization</kwd>
        <kwd>Clinical decision support system</kwd>
        <kwd>Atrial fibrillation</kwd>
        <kwd>Congestive heart failure</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>Our visualization is comprised of two diferent sections, the main RR interval
sequence diagram, and the RR interval histograms. Figure 1 shows the RR
sequence diagram and a seconds by hours scatter plot of each RR interval in
the signal. With 24 hours of data, this sequence diagram could contain over
200,000 RR intervals causing them to be clustered together making it virtually
impossible to trace between consecutive RR intervals. To address this, each RR
interval point has an opacity of 0.1 allowing us to focus on the dense regions of
RR intervals by reducing the distraction caused by outliers within the signals.
We also provide the ability to apply a simple central moving average (mAvg)
across the sequence diagram allowing us to smooth out short-term fluctuations
and highlight long term-trends. The window size of the mAvg can be defined to
suit the needs of the user. A larger window size will reduce the efect of short
term variance highlighting the long term trend, while a smaller window size will
highlight the short term variance.</p>
      <p>The mAvg is accompanied by upper and lower standard deviation bands
calculated by finding the standard deviation of the subset of values above the
mean and below the mean for each average within the mAvg. The upper standard
deviation is then added to the respective mean and plotted on the sequence
diagram while the lower one is subtracted from the respective mean. By plotting
these upper and lower deviations for each average within the mAvg, we obtain two
bands that surround the mAvg. The amount of variation that exists throughout
the sequence can be inferred from a quick glance at the width of the bands. A
larger band represents a greater amount of variance than a thinner band.</p>
      <p>
        Two histograms are included, an RR interval histogram (IH) and an RR
interval diferences histogram (IDH) (see Fig. 1). Significant amounts of data
from a signal can be encoded into these normalized histograms, allowing easy
comparison across time scales. For example a large amount of data (24 hour
signal) can be compared to a small amount (4 hours) because we care more
about the shape (distribution) of the histograms than we do the total values.
It has been shown that when comparing or classifying signals, using both the
IH and IDH combined provides more reliable qualitative information than the
IH alone [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. The IH displays the spread of the heart rates around the mean
value while the IDH shows the smoothness of the rate changes. Our visualization
expands on this concept by incorporating the mAvg into the creation of each
histogram. During exploration, the RR sequence diagram can be brushed to
select a particular period of interest within the signal. The histograms will then
update to represent the selected subset of data. Depending on whether or not
a mAvg has been applied to the signal, the histograms will represent either the
RR average intervals or just the RR intervals themselves. Without a mAvg the
IDH will show the diferences between each consecutive RR interval allowing
us to see the immediate short term variance. With a mAvg it will display the
diferences between the RR intervals and the respective mean within the mAvg,
generalizing the IDH which lowers the amount short term variance can afect the
ifnal distribution.
3
      </p>
    </sec>
    <sec id="sec-2">
      <title>Use Case Scenario</title>
      <p>
        A decreased amount of short term variability has been found to be one of the
indicators of congestive heart failure [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. By looking at the RR interval diferences
histogram in Fig. 2 we can see that it has a thin shape. The RR interval histogram
is wider because it shows us the total distribution throughout the whole signal —
the long term variation. The IDH shows the diference of each consecutive RR
interval which isn’t afected by the overall baseline shifts throughout the signal
and allows us to see the short term variance. Therefore, considering low short
term variation is an indicator of CHF, the IDH may be useful in diagnosing CHF
as the short term variation can be directly inferred from the chart.
      </p>
      <p>By adding a moving average and standard deviation bands to the CHF signal
(Fig. 3) we obtain a figure that we would expect to see from a signal with low
short term variation. The short term variation causes the bands to be incredibly
tight around the moving average. Considering the RR interval points are
semitransparent, these tight bands allow us to confirm that there are no clusters of
RR intervals surrounding the moving average as they would increase the width
of the bands. Such clusters would increase the amount of short term variance
lowering the probability of CHF. It should also be noted that when adding a
moving average the IDH is calculated as diferences from the moving average,
which is why the IDH in Fig. 3 is wider than the IDH in Fig. 2.</p>
      <p>Therefore by using both techniques together we can infer a low short term
variability from the IDH when no moving average has been applied and that no
hidden clusters surround the moving average once the bands have been applied.
Both of these are evidence that CHF could be the cause of such cardiac patterns.</p>
      <p>
        Atrial fibrillation characteristically appears in the RR interval sequence as
a curtain of points (essentially a large amount of short term variation). Atrial
ifbrillation can appear for short or long periods of time and it can be almost
impossible to determine a baseline average at a glance during a period of atrial
ifbrillation. Without the added moving average it has been shown that atrial
ifbrillation produces wide triangular histogram patterns within the IDH [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. This
can be seen in Fig. 4: notice how a specific region within the signal that has been
speculated to be an instance of AF has been brushed to analyze the histograms
for that section.
      </p>
      <p>
        The moving average after being applied to the AF signal allows us to see
the baseline throughout the period of AF as shown in Fig. 5. Notice how the
standard deviation bands around the moving average during the period of atrial
ifbrillation are wide, showing that there is a high amount of short term variation.
The IDH also produces a poly-modal distribution during the period of AF. Not
all instances of AF have been shown to produce these poly-modal distributions
but (after analyzing all the signals in the nsrdb [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]) we did not find a case where
a normal sinus rhythm has produced a distribution such as this. So it may be
inferred that if the IDH has a poly-modal distribution is after adding the moving
average, the signal is not from the normal sinus rhythm database.
We performed a counterbalanced within subjects study to observe if 20 non-expert
participant students would be able to accurately diagnose signals as normal or
abnormal.
Training and testing datasets were drawn from several data bases within the
PhysioNet catalogue [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. Each dataset was a patient ECG record lasting a
minimum of 4 hours (10 hours on average). Records were collected for congestive
heart failure (CHF) [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], atrial fibrillation (AF) [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ], unknown cause abnormal heart
rate (Abnormal) [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ], and normal sinus rhythm (Normal) [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] were included. 10
testing datasets were randomly selected from each dataset, as well as 4 training
datasets for each of CHF, AF, and Normal.
Participants were not experts in heart rate analysis or experienced with interactive
visual analytics. Participants were all undergraduate students studying in a range
of programs, from health science to business and information technology.
      </p>
      <p>The main factor varied in the experiment was interface style (2 levels): Basic
style, and Statistical style. Participants were trained on how to use each style
and how to find characteristics of atrial fibrillation (AF), congestive heart failure
(CHF), and normal sinus rhythm using that style. After hands-on training,
participants were able to explore 6 training datasets provided for each style to
practice making a diagnosis. During this stage participants were aware of the
proper diagnosis for each signal provided and the exploration was intended for
them to identify patterns within the signals they could recognize in the trial
datasets.</p>
      <p>
        After training they were tasked with diagnosing datasets as normal or
abnormal. If they believed it to be abnormal, they also stated whether any instances
of AF or CHF exist. After selecting their diagnosis they were asked to rate the
confidence of their diagnosis on a scale of 1 to 7, 1 being not confident at all and
7 being completely confident. The Basic style was created to replicate the time
domain analysis tools provided by a state-of-the-art HRV analysis toolkit [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
The Basic style tasked participants with diagnosing signals while only having
access to the RR interval sequence diagram and the IH. The Statistical style
included everything that was available in the Basic style as well as the moving
average, standard deviation bands, and the IDH.
One of the greatest strengths of visualizing RR interval data is that we can get
an overview of an extended period. To bridge the connection between a long-term
general overview and short-term ECG strips, we plan to incorporate the ECG
strip into the visualization. Furthermore, discovering where to investigate the
short-term data is a challenge. We plan a classifier to determine which areas of
the ECG signal are of interest. These areas of interest would be highlighted to
aid in the exploration of the signal, allowing faster drill down to the ECG data.
6
      </p>
    </sec>
    <sec id="sec-3">
      <title>Conclusion</title>
      <p>
        Given the success of the participants in diferentiating between normal and
abnormal signals after a brief training period, it can be said that patterns can be
easily recognized when searching for AF or CHF. Assuming these conditions are
not the only ones with easily recognizable characteristics, our visualization will
be able to support clinical decisions as patterns of other conditions are discovered.
A promising condition that could be investigated would be diabetic autonomic
neuropathy due to the low variation of heart rate caused by this condition [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ].
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Baim</surname>
            ,
            <given-names>D.S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Colucci</surname>
            ,
            <given-names>W.S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Monrad</surname>
            ,
            <given-names>E.S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Smith</surname>
            ,
            <given-names>H.S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Wright</surname>
            ,
            <given-names>R.F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lanoue</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          , Gauthier,
          <string-name>
            <given-names>D.F.</given-names>
            ,
            <surname>Ransil</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.J.</given-names>
            ,
            <surname>Grossman</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            ,
            <surname>Braunwald</surname>
          </string-name>
          , E.:
          <article-title>Survival of patients with severe congestive heart failure treated with oral milrinone</article-title>
          .
          <source>J American College of Cardiology (3)</source>
          ,
          <fpage>661</fpage>
          -
          <lpage>670</lpage>
          (
          <year>Mar 1986</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Cashman</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          :
          <article-title>The use of RR interval and diference histograms in classifying disorders of sinus rhythm: Original articles</article-title>
          .
          <source>Journal of Medical Engineering &amp; Technology</source>
          <volume>1</volume>
          (
          <issue>1</issue>
          ),
          <fpage>20</fpage>
          -
          <lpage>28</lpage>
          (
          <year>1977</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Goldberger</surname>
            ,
            <given-names>A.L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Amaral</surname>
            ,
            <given-names>L.A.N.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Glass</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hausdorf</surname>
            ,
            <given-names>J.M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ivanov</surname>
            ,
            <given-names>P.C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mark</surname>
            ,
            <given-names>R.G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mietus</surname>
            ,
            <given-names>J.E.</given-names>
          </string-name>
          , Moody, G.B.,
          <string-name>
            <surname>Peng</surname>
            ,
            <given-names>C.K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Stanley</surname>
          </string-name>
          , H.E.: PhysioBank, PhysioToolkit, and
          <article-title>PhysioNet: Components of a new research resource for complex physiologic signals</article-title>
          .
          <source>Circulation</source>
          <volume>101</volume>
          (
          <issue>23</issue>
          ),
          <fpage>e215</fpage>
          -
          <lpage>e220</lpage>
          (
          <year>2000</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Greenwald</surname>
          </string-name>
          , S.D.:
          <article-title>Development and analysis of a ventricular fibrillation detector</article-title>
          .
          <source>Master's thesis</source>
          , Massachusetts Institute of Technology (
          <year>1986</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5. Moody, G.B.,
          <string-name>
            <surname>Mark</surname>
          </string-name>
          , R.G.:
          <article-title>A new method for detecting atrial fibrillation using r-r intervals</article-title>
          .
          <source>Computers in Cardiology 10</source>
          ,
          <fpage>227</fpage>
          -
          <lpage>230</lpage>
          (
          <year>1983</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>van</surname>
            Ravenswaaij-Arts,
            <given-names>C.M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kollee</surname>
            ,
            <given-names>L.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hopman</surname>
            ,
            <given-names>J.C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Stoelinga</surname>
            , G.B., van Geijn,
            <given-names>H.P.</given-names>
          </string-name>
          :
          <article-title>Heart rate variability</article-title>
          .
          <source>Annals of internal medicine 118(6)</source>
          ,
          <fpage>436</fpage>
          -
          <lpage>447</lpage>
          (
          <year>1993</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Singh</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bharti</surname>
          </string-name>
          , N.:
          <article-title>Software tools for heart rate variability analysis</article-title>
          .
          <source>International Journal of Recent Scientific Research</source>
          <volume>6</volume>
          (
          <issue>4</issue>
          ),
          <fpage>3501</fpage>
          -
          <lpage>3506</lpage>
          (
          <year>2015</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Tarvainen</surname>
            ,
            <given-names>M.P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Niskanen</surname>
            ,
            <given-names>J.P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lipponen</surname>
            ,
            <given-names>J.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ranta-Aho</surname>
            ,
            <given-names>P.O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Karjalainen</surname>
            ,
            <given-names>P.A.</given-names>
          </string-name>
          :
          <article-title>Kubios hrv-heart rate variability analysis software</article-title>
          .
          <source>Computer Methods and Programs in Biomedicine</source>
          <volume>113</volume>
          (
          <issue>1</issue>
          ),
          <fpage>210</fpage>
          -
          <lpage>220</lpage>
          (
          <year>2014</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>