<!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>Optimizing Factors In uencing on Accuracy of Biometrical Cardiometry</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Marat R. Bogdanov</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Aleksander A. Dumchikov</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Vadim M. Kartak</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Aigul I. Fabarisova</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>M. Akmullah named after Bashkir State Pedagogical University</institution>
          ,
          <addr-line>Ufa</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Ufa State Aviation Technical University</institution>
          ,
          <addr-line>Ufa</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <fpage>61</fpage>
      <lpage>65</lpage>
      <abstract>
        <p>The paper is about some aspects concerning person biometric identi cation based on using of electrocardiograms. Signal preprocessing routing is considered in the paper. Classi cation was carried out with support vector machines algorithm. Tuning of hyper parameters of classi cation is considering. The problem of person biometric identi cation concerns classi cation problems. To solve it, we have to consider algorithms from some nite set and choose an algorithm that gives the least error of the forecast [3]. Let's introduce some notation. Let us suppose X is a space of objects.</p>
      </abstract>
      <kwd-group>
        <kwd>Biometric person identi cation</kwd>
        <kwd>Electrocardiogram chine learning</kwd>
        <kwd>Support vector machines</kwd>
        <kwd>Hyper parameter tuning</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>
        Various biometric methods of person identi cation are getting more popular.
Fingerprinting, face, voice and retina recognition are widely used in various
security systems. The vulnerabilities of traditional methods of biometric identi
cation were revealed over time. Researchers are increasingly turning their attention
to such person biometric features as electrocardiograms, electroencephalograms
and DNA [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. In this paper, we would like to discuss some practical aspects of
person identi cation using ECG.
      </p>
    </sec>
    <sec id="sec-2">
      <title>Motivation and Aim</title>
      <p>Copyright c by the paper's authors. Copying permitted for private and academic purposes.
In: S. Belim et al. (eds.): OPTA-SCL 2018, Omsk, Russia, published at http://ceur-ws.org
Xl = (xi; yi)li=1</p>
      <p>yi = y (xi);
At = fa : X ! Y g
t : (X</p>
      <p>Y )l ! At
are a models of algorithms, t 2 T , T is a number of algorithms under
consideration.
are learning methods. It is required to nd a method t with the best generalizing
power.</p>
      <p>When nding a method t, we often have to solve the following subtasks:
{ Choice of the best model At (model selection).
{ Choice of learning method t for a given model At (in particular,
optimization of hyperparameters).
{ Features selection:</p>
      <p>F = ffj : X ! Dj : j = 1; :::; ng
is a set of features. The method of learning j uses only features J
To assessment the quality of learning by precedents it's used:
L(a; x) is a cost function of algorithm a on the object x.</p>
      <p>F .</p>
      <p>
        Y is a set of answers.
is a training set, l is a sample size.
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
is a functional of accuracy a on X. In this case we consider an internal quality
criterion that is measured on the training set Xl:
and an external criterion evaluating the quality of learning on hold-out set Xk[
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]:
l
Q(a; Xl) = 1 X L(a; xi)
l
      </p>
      <p>i=1</p>
      <p>Q (Xl) = Q( (Xl); Xl)</p>
      <p>Q (Xl; Xk) = Q( (Xl); Xk):
In the paper presented we will consider such aspects of person biometric
identi cation as feature selection, model selection, choice of methods (tuning of
hyperparameters), assessment of the quality of learning.
3</p>
    </sec>
    <sec id="sec-3">
      <title>Feature Selection</title>
      <p>
        We used the MGH/MF Waveform Database hosted at physionet.org
resource [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. The Massachusetts General Hospital/Marquette Foundation
(MGH/MF) Waveform Database is a comprehensive collection of electronic
recordings of hemodynamic and electrocardiographic waveforms of stable and
unstable patients in critical care units, operating rooms, and cardiac
catheterization laboratories. It is the result of a collaboration between physicians,
biomedical engineers and nurses at the Massachusetts General Hospital. The database
consists of recordings from 250 patients and represents a broad spectrum of
physiologic and pathophysiologic states. Individual recordings vary in length from 12
to 86 minutes, and in most cases are about an hour long [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ].
      </p>
      <p>
        The typical recording includes three ECG leads, arterial pressure, pulmonary
arterial pressure, central venous pressure, respiratory impedance, and airway
CO2 waveforms. The raw sampling rate of 1440 samples per second per signal
was reduced by a factor of two to yield an e ective rate of 360 samples per
second per signal relative to real time [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ].
      </p>
      <p>
        When preprocessing stage we used a biopsy python library by John Reid. The
package enables the development of Pattern Recognition and Machine Learning
work ows for the analysis of biosignals including ECG [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. Using biopsy we
extracted rst lead from electrocardiogram and performed a low pass lter for
reducing of redundancy. After applying of low-pass lter R-peaks was extracted
from ECG-signal using a wfdb python library by Chen Xie and Julien Dubiel [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ].
The software allow extract peaks and QRS -cycles from electrocardiograms. We
choice amplitude and temporal features of Q,R and S -peaks ( Qx, Qy, Rx, Ry,
Sx, Sy). In total we had 6 features. Feature table together label class vector
were randomly splitted into training set and testing set in the ratio of 75:25 for
further cross validation. We learned a classi er on training set and performed
measuring of classifying accuracy on testing set.
4
      </p>
    </sec>
    <sec id="sec-4">
      <title>Model Selection</title>
      <p>
        We used a Support Vector Machines (SVM) algorithm for classi cation.
Support Vector Machines are based on the concept of decision planes that de ne
decision boundaries. A decision plane is one that separates between a set of
objects having di erent class memberships. SVM is primarily a classier method that
performs classi cation tasks by constructing hyperplanes in a multidimensional
space that separates cases of di erent class labels. SVM supports both
regression and classi cation tasks and can handle multiple continuous and categorical
variables [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ].
      </p>
      <p>To construct an optimal hyperplane, SVM employs an iterative training
algorithm, which is used to minimize an error function. According to the form of
the error function, SVM models can be classi ed into four distinct groups:
{ Classi cation SVM Type 1 (also known as C-SVM classi cation)
{ Classi cation SVM Type 2 (also known as nu-SVM classi cation)
{ Regression SVM Type 1 (also known as epsilon-SVM regression)
{ Regression SVM Type 2 (also known as nu-SVM regression)</p>
      <p>We used a Classi cation SVM Type 1 (also known as C-SVM classi cation)
model.</p>
    </sec>
    <sec id="sec-5">
      <title>Classi cation SVM</title>
    </sec>
    <sec id="sec-6">
      <title>Type 1</title>
      <p>For this type of SVM, training involves the minimization of the error function:</p>
      <p>N
1 wT w + C X
2
i=1
i
subject to the constraints:
yi(wT (xi) + b)
1
i
0; i = 1; :::; N;
where C is the capacity constant, w is the vector of coe cients, b is a constant,
and i represents parameters for handling nonseparable data (inputs). The index
i labels the N training cases. Note that y 2 1 represents the class labels and
xi represents the independent variables. The kernel is used to transform data
from the input (independent) to the feature space. It should be noted that the
larger the C, the more the error is penalized. Thus, C should be chosen with
care to avoid over tting.
(9)
(10)
(11)
6</p>
    </sec>
    <sec id="sec-7">
      <title>Kernel Functions</title>
      <p>K(Xi; Xj ) =
8
&gt;
&gt;
&lt;</p>
      <p>Xi Xj
( Xi Xj + C)d
&gt; exp( jXi Xj j2)
&gt;: tanh( Xi Xj + C)</p>
      <p>Linear &gt;9
P olynomial =&gt;</p>
      <p>SiRgBmFoid ;&gt;&gt;
;
where K(Xi; Xj ) = (Xi) (Xj ) that is, the kernel function, represents a
dot product of input data points mapped into the higher dimensional feature
space by transformation .
7</p>
    </sec>
    <sec id="sec-8">
      <title>Gamma is an Adjustable Parameter of Certain Kernel</title>
    </sec>
    <sec id="sec-9">
      <title>Functions</title>
      <p>
        The RBF is by far the most popular choice of kernel types used in Support
Vector Machines. This is mainly because of their localized and nite responses
across the entire range of the real x-axis [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ].
      </p>
      <p>
        Support vector machine classi er supported by sklearn python library uses as
default following hyper parameters: C=1.0, kernel='rbf', gamma='auto'. When
using of default parameters while performing of classi cation of
electrocardiograms we had accuracy score equal to 0.93. We tuned hyper parameters of
classi cation with Grid Search procedure varying C parameter in range of [
        <xref ref-type="bibr" rid="ref1">1, 10,
100, 1000</xref>
        ], kernel in range of ['linear, 'rbf"], and gamma in range of [1e-3, 1e-4].
After performing of tuning we had the following best parameters set: 'kernel':
'rbf', 'C': 10, 'gamma': 0.001. Using these parameters we had accuracy score
equal to 0.99.
      </p>
    </sec>
    <sec id="sec-10">
      <title>Results and Discussion</title>
      <p>During the preprocessing of electrocardiograms we extracted rst leads of
signal and performed low-pass lter for reducing redundancy. Then we extracted
cardiac cycles from the leads and extracted Q, R and S peaks from cardiac cycles.
Using amplitude and temporal features of peaks we composed a feature table
containing 6 features (Qx, Qy, Rx, Ry, Sx, Sy) and class labels vector y. Further
we randomly splitted a feature table and class labels vector into training set
and testing set on ration of 75:25 for further cross-validation. Training set was
used for learning a classi er and testing set was used for assessment of quality
of learning. SVM classi er supported by sklearn python library using default
options show accuracy score equal to 0.93. We found the best hyper parameters
are following: 'kernel': 'rbf', 'C': 10, 'gamma': 0.001. Using these parameters we
could improved accuracy score up to 0.99.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Abdulmonam</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ahlal</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fawzia</surname>
          </string-name>
          , E.:
          <article-title>Vulnerabilities of biometric authentication</article-title>
          .
          <source>Threats and countermeasures. IJICT</source>
          <volume>4</volume>
          (
          <issue>11</issue>
          ),
          <fpage>947</fpage>
          -
          <lpage>958</lpage>
          (
          <year>2014</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Goldberger</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Amaral</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Glass</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hausdor</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ivanov</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mietus</surname>
          </string-name>
          , J., Moody, G.,
          <string-name>
            <surname>Peng</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Stanley</surname>
          </string-name>
          , H.: 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>1</fpage>
          -
          <lpage>6</lpage>
          (
          <year>2000</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          <article-title>3. Machine Learning course by K</article-title>
          .V. Voroncov: http://www.MachineLearning.ru/ wiki [On-line; accessed 12-April-2011]
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <given-names>Support</given-names>
            <surname>Vector</surname>
          </string-name>
          <article-title>Machines (SVM) Introductory Overview</article-title>
          : http://www.statsoft. com/Textbook/Support-Vector-Machines [On-line; accessed 26-May-2012]
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>5. The biosppy Toolbox: http://biosppy.readthedocs.io/en/stable/ [On-line; accessed 24-March-2016]</mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>6. The WFDB Python Toolbox: https://pypi.python.org/pypi/wfdb [On-line; accessed 11-July-2015]</mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Ting-Fan</surname>
            ,
            <given-names>W.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chih-Jen</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Weng</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          , Singer.,
          <string-name>
            <surname>Y</surname>
          </string-name>
          . (eds.):
          <article-title>Probability estimates for multi-class classi cation by pairwise coupling</article-title>
          .
          <source>Journal of Machine Learning Research</source>
          <volume>5</volume>
          ,
          <fpage>975</fpage>
          -
          <lpage>1005</lpage>
          (
          <year>2004</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Welch</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ford</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Teplick</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rubsamen</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          :
          <article-title>The Massachusetts General Hospital-Marquette Foundation hemodynamic and electrocardiographic database comprehensive collection of critical care waveforms</article-title>
          .
          <source>J Clinical Monitoring</source>
          <volume>7</volume>
          (
          <issue>1</issue>
          ),
          <fpage>96</fpage>
          -
          <lpage>97</lpage>
          (
          <year>1991</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>