<!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>
      <journal-title-group>
        <journal-title>Nanotechnology (ITNT-2015), CEUR Workshop Proceedings</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.18287/1613-0073-2015-1490-346-354</article-id>
      <title-group>
        <article-title>The Big Data mining to improve medical diagnostics quality</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Ilyasova N.Yu.</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Kupriyanov A.V.</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Samara State Aerospace University, Image Processing Systems Institute, Russian Academy of Sciences</institution>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2015</year>
      </pub-date>
      <volume>1490</volume>
      <fpage>346</fpage>
      <lpage>354</lpage>
      <abstract>
        <p>The paper offers a method of the big data mining to solve problems of identification of cause-and-effect relationships in changing diagnostic information on medical images with different kinds of diseases. As integrated indices of the fundus vessels and coronary heart blood vessels we have used a global set of geometric features which is supposed to be a rather complete characteristic of diagnostic images and allows to make a successful diagnosis of vascular malformations. To evaluate informativity of vascular diagnostic features based on a classification efficiency criterion and in order to form new features required to improve a diagnostics quality a method of discriminative analysis of sample data has been considered. A filtration method of invalid data is proposed using a clustering algorithm to improve a performance quality of the developed algorithm of the discriminative analysis of feature vectors.</p>
      </abstract>
      <kwd-group>
        <kwd>human vascular system</kwd>
        <kwd>image processing</kwd>
        <kwd>diagnostic features</kwd>
        <kwd>discriminative analysis</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        The very large data (‘big data’) mining is a key problem of modern information
technologies. In accordance with the Forecast of the Scientific and Technology
Development of the Russian Federation for the period until 2030 approved by the
Chairman of the Government of the Russian Federation Dmitry A. Medvedev
promising research areas include the “Technologies of data processing and analysis”
including methods and techniques of collection, processing, analysis and storage of
very large volumes of information. The purpose of the paper is to develop methods
and algorithms of the bid data mining to solve problems of identification of
causeand-effect relationships in changing diagnostic information on medical images with
different kinds of diseases. It is offered to use a single approach to the analysis of
different image classes based on evaluation of combined geometric parameters of
selected areas of interest which are considered to be a basic feature set for further
diagnostic analysis [
        <xref ref-type="bibr" rid="ref1 ref2">1-2</xref>
        ]. In order to identify images based on the big data mining
using methods of the discriminative analysis a technique of efficient feature space
formation has been developed [
        <xref ref-type="bibr" rid="ref3 ref4 ref5">3-5</xref>
        ]. As integrated indices of the fundus vessels and
coronary heart blood vessels it is proposed to use the global set of geometric features
which allow to make the successful diagnosis of vascular malformations [
        <xref ref-type="bibr" rid="ref6 ref7">6-7</xref>
        ]. Based
on the specified methods new distributed technologies and software are developed to
provide the remote image processing, analysis and understanding which are intended
for implementation in automated telemedicine systems. The methods being developed
are aimed at improvement of the medical diagnostics quality due to obtaining
objective numerical estimates of biomedicine image parameters using large volumes
of accessible information.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. Diagnostic Image Mining Information Technology</title>
      <p>The diagnostic image mining information technology includes a method of
formation of the efficient feature space to classify a predefined image set.</p>
      <p>
        A highlighting technique of diagnostically significant information on blood vessel
images is based on a new generalized mathematical model of blood vessels for two
classes of diagnostic images, i.e. the fundus blood-vessel system and coronary blood
vessels, characterized by a set of geometric parameters [
        <xref ref-type="bibr" rid="ref6 ref7 ref8">6-8</xref>
        ].
      </p>
      <p>A geometric approach to formation of diagnostic features, which, unlike traditional
abstract spectral-and-correlation features are well accustomed and understandable by
medical professionals, demonstrate a good visual effect and take into account an
object’s specific character, allows to finally increase a diagnostics efficiency.</p>
      <p>In order to select the most successful features we use their correlation with results
of expert evaluations, as well as a dispersion analysis of learning samples or the
diagnostic error analysis using particular characteristics. The effectiveness of different
features is evaluated for automated diagnosis problems and proper recommendations
are articulated on how to use different groups of features in medical practice.</p>
      <p>The diagnostic image mining information technology includes the following
advanced techniques and algorithms:
─ the technique and the algorithm of increasing a degree of the feature informativity
based on the discriminative analysis and formation of an optimal learning sample
to learn a disease diagnosis expert system;
─ an estimation method of a class separatability which is not influenced by
distribution of objects through classes and is independent from a used classifier;
─ the algorithm of decrease of feature space dimensions and formation of new
informative features that maximize a separatability criterion based on methods of
the discriminative analysis and allow to increase a diagnosis accuracy of a
pathology degree;
─ a technique of optimum learning sample formation to learn a diagnostic system
based on exemption of anomalous observations that will also enable to increase a
disease diagnosis accuracy.</p>
      <p>
        Problem-oriented and distributed software solutions of the analysis of medical and
diagnostic images are developed to detect pathological changes including software
tools for quantitative estimates of the pathology degree based on expert findings and
the proposed methods of classification. They are intended to ensure the user with an
opportunity to manage a decision-making and analysis process [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]. Automated
systems of the quantitative feature analysis make it possible to standardize a
diagnosing process, considerably reduce an examination time and decrease its cost.
The systems allow to carry out the analysis of subclinical morphological changes of
pathomorphological components, computerize diagnosis stages and make a
quantitative monitoring of pathological changes in diagnostic images. The peculiarity
is the use of expert system components, i.e. a database of diagnostic features, the
correlative, discriminative and cluster analysis of the feature space, and a prognosis of
the pathology degree based on expert estimates.
      </p>
      <p>
        A system of classification and diagnostic researches [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ] provides tools for the
correlative and discriminative analysis to form the informative feature space,
resources to form an optimal feature sample based on the criterion of separation
efficiency by pathology groups, and facilities for the cluster analysis to filtrate the
learning sample in order to eliminate invalid data and to obtain standard feature
values in accordance with pathology groups. The data mining system allows the user
to obtain a proper pathology degree, standard feature values for each degree of
disease pathology and a predicted disease development, and will provide proper
diagnostic decisions.
      </p>
    </sec>
    <sec id="sec-3">
      <title>3. The Discriminative Analysis to Form the Informative Feature Space</title>
      <p>
        Proper researches have been performed together with medical professionals from
the Medical and Stomatological University (Moscow), Ophthalmology Academic
Department, based on a digital image analysis of the fundus. Diagnostics techniques
for ophthalmic diseases have been developed based on evaluation of global vascular
characteristics (features). The paper considers geometric characteristics proposed in
[
        <xref ref-type="bibr" rid="ref10 ref6 ref7 ref8 ref9">6-10</xref>
        ]. These characteristics involve the following: the mean diameter, straightness,
beads-looking shape, amplitude of thickness variations, frequency of thickness
variations, thickness tortuosity, amplitude of path variations, frequency of path
variations and path tortuosity, which correspond to diagnostic features of the fundus
vessels. If two or more classes available (in our case there are 5 classes including a
normal feature and 4 degrees of diabetic retinopathy and diabetes), the objective of
feature selection is to select those features which are the most efficient in accordance
with the class separability [
        <xref ref-type="bibr" rid="ref5 ref8">5,8</xref>
        ]. In the discriminative analysis the class separability
criteria are formed using scatter matrixes inside classes and scatter matrixes between
classes [
        <xref ref-type="bibr" rid="ref11 ref12">11,12</xref>
        ].
      </p>
      <p>The scatter matrix inside classes demonstrates a variety of objects with respect to
mean vectors of classes: W  kg1(Xk  xk )(Xk  xk ) , where k – class data will
correspond to the mean vector xk  [x1k x21k
xpk ] , and g – is a total number of
classes.</p>
      <p>
        Elements of the scatter matrix between classes B is counted by the formula:
bij  kg1 nk  xik  xi  xjk  xj , i, j  1, p , xi  1 n k1 nk xik – is a mean feature
g
value of  - feature in all classes, nk – is a number of objects in k - class,
xik  1 nk nmk1 xikm – is the mean feature value in class k , and xikm – is a value of i
feature for m -object in k - class. The matrixes W and B contain all basic information
about interrelationships inside and between classes. In order to obtain the class
separability criterion some number is to be associated to these matrixes. This number
should be increased with the increase of scattering between classes or with the
decrease of scattering inside classes. For this purpose the following criterion is more
frequently used [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ]: J1  tr(T-1B) , where T  B  W . The greater the value of the
criterion - the more separability of classes. The following algorithm presented in
Figure 1 has been developed to form new features.
      </p>
      <sec id="sec-3-1">
        <title>Preparatory stage</title>
      </sec>
      <sec id="sec-3-2">
        <title>Receiving model data</title>
      </sec>
      <sec id="sec-3-3">
        <title>Objects classification</title>
      </sec>
      <sec id="sec-3-4">
        <title>Computation of features</title>
      </sec>
      <sec id="sec-3-5">
        <title>Formation of new features</title>
      </sec>
      <sec id="sec-3-6">
        <title>Identification of eigen values and vectors</title>
      </sec>
      <sec id="sec-3-7">
        <title>Identification of discriminative functions coefficients</title>
      </sec>
      <sec id="sec-3-8">
        <title>Formation of features based on values of discriminative functions</title>
      </sec>
      <sec id="sec-3-9">
        <title>Findings interpretation</title>
      </sec>
      <sec id="sec-3-10">
        <title>Analysis of basic features</title>
        <p>Identification of correlation
matrix
Identification of scatter matrix
Computation of separability
criteria
Feature selection to form new
features
Classification using selected
features
Computation of structural
coefficients</p>
      </sec>
      <sec id="sec-3-11">
        <title>Analysis of formed features</title>
        <p>Analysis of discriminative
functions</p>
        <p>Identification of correlation
matrix and scatter matrixes for
new features
Computation of criteria for new
features
Classification using new
features</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Experimental Studies</title>
      <p>A number of studies of vascular malformations have been undertaken for diabetic
retinopathy (DR) of 151 patients suffered with diabetes (D) based on the digital image
analysis of the fundus. After processing of images the sample amounted to 8175
measurements including 1490 arteriols of the 1st-order, 2345 arteriols of the 2nd-order,
1960 venules of the 1st-order, and 2380 venules of the 2nd-order. Medical
professionals separately consider venules and arteriolas of the 1st and 2nd order
(GROUP 1 – GROUP 4) since different tendencies of vessel changes can be observed
in these GROUPS at different pathology stages. Figure 2 shows some examples of
diagnostic images of the fundus for different stages of diabetes (1D, 4D) and a
measurement sketch of the vascular system that gives an order of the vessel
examination accepted by ophthalmologists.
4 D</p>
      <p>At the examination of features it was concluded that there are two high-correlated
characteristic groups. The first group includes those features which describe path
parameters (e.g. the path straightness and tortuosity), and the second group includes
the features characterizing a vessel radius function (e.g. the radius tortuosity and a
beads-factor). Figure 3 gives values of the separability criteria of single features for
each of the four GROUPS (venules and arteriols of the 1st and 2nd orders).</p>
      <p>Figure 3 shows that maximum feature values belong to the vessel features which
are not particularly significant for medical professionals in diagnosing. Besides, in
proper pathological conditions each group may have vessels which do not comply
with the given pathology (for instance, the norm feature). It may be thus concluded
that the sample may contain some noisy data.</p>
      <p>In order to eliminate noise in the sample it is proposed to filter down an original
sample using the clustering approach of k-means. Each GROUP was divided into 5
clusters in accordance with original clusters inside each of the GROUPS. Feature
vectors which were not included into a correct cluster were removed from the sample.
Values of individual separability criteria for features inside the GROUPS after
filtration are given in Figure 4.</p>
      <p>
        Analyzing Figure 4 we may conclude that in the sample obtained after filtration the
features with the greatest separability criterion are the characteristics which are
considered by medical professionals as specifically diagnostically significant in visual
diagnosing of pathology that corresponds to an information letter [
        <xref ref-type="bibr" rid="ref13">13</xref>
        ]. The
discriminative analysis algorithm has been applied for original and filtered samples.
In order to form new features we have completely enumerated original features to
search a combination of new features that has maximized the separability criterion.
The result was that we have obtained a set of four features for all GROUPS of either
sample.
      </p>
      <p>
        The classification error was evaluated for the GROUPS obtained thereby. The error
was estimated by means of the U-approach [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ]. Two samples were formed, i.e. the
learning and test samples. The classifier was tuned using the learning sample based on
the SVM approach (Support Vector Machines), and the test sample was thereby
classified. Objects of the learning sample which are not contained within the test
sample are only used for the classifier synthesis. There are many opportunities to
implement the U-approach. In order to evaluate probability of the classification error
in researches a one-object-elimination method was used. The findings are illustrated
in Table 1 where criterion values before the discriminative analysis are presented for
the combination of original sample quads with the best separability values inside the
GROUPS.
      </p>
      <p>Analyzing the obtained findings it may be concluded that filtration allows not only
to increase feature separability criteria and decrease the classification error, but also to
identify diagnostically significant features. Researches made on four GROUPS of
blood vessels have shown that it is important for each GROUP to have its own set of
diagnostic features that is proved by clinical research studies of medical practitioners.
For example, the mean diameter of the blood vessel acts differently for venules and
arteriols under pathological changes. The research results have showed that the
application of the proposed feature formation technique allowed not only to eliminate
invalid data, but also led to a reduction in errors of classification. This resulted to the
increase of the separability criterion by 39% for GROUP 1, by 42% - for GROUP 2,
by 24% - for GROUP 3 and by 39% - for GROUP 4. Besides, we have obtained some
additional information based on the used features such as their informativity and have
identified relationships between some features.</p>
    </sec>
    <sec id="sec-5">
      <title>5. Conclusion</title>
      <p>A new technology has been developed for the data mining to identify
cause-andeffect relationships of changing diagnostic information on medical images with
different kinds of diseases including the filtration of invalid data and the
discriminative analysis based on maximizing the separability criterion.</p>
      <p>The algorithm based on the selection of features with the biggest separability
criterion and on a complete enumeration followed by new features maximizing this
criterion has been developed. Separability criteria are formed using the scattering
matrixes between and inside classes. As a result of the discriminative analysis the
best features have been defined for each group of vessels based on the separability
criterion. It is shown that each of the four GROUPS of vessels has efficiently used its
own set of global geometric features that has been proved by clinical researches. The
classification error was calculated for each group of vessels before and after the
algorithm performance. It is shown that the proposed technology of the feature space
analysis by groups, including the algorithm of filtration of sampled data and the
algorithm of formation of the efficient feature space, allowed to increase the
classification efficiency of vessels by “normal-feature” classes and different degrees
of “pathology”. The classification error was hereby reduced to 2.3%-3.5% for
different pathology GROUPS.</p>
    </sec>
    <sec id="sec-6">
      <title>Acknowledgements</title>
      <p>This work was partially supported by the Ministry of education and science of the
Russian Federation in the framework of the implementation of the Program of
increasing the competitiveness of SSAU among the world’s leading scientific and
educational centers for 2013-2020 years; by the Russian Foundation for Basic
Research grants (# 14-01-00369-a, # 14-07-97040-p_ povolzh'e_a, # 15-29-03823, #
15-29-07077); by the ONIT RAS program # 6 “Bioinformatics, modern information
technologies and mathematical methods in medicine” 2015.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <given-names>Ilyasova</given-names>
            <surname>NYu</surname>
          </string-name>
          .
          <article-title>Diagnostic computer complex for vascular fundus image analysis</article-title>
          ,
          <source>Biotehnosfera</source>
          ,
          <year>2014</year>
          ;
          <volume>33</volume>
          (
          <issue>3</issue>
          ):
          <fpage>20</fpage>
          -
          <lpage>24</lpage>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <given-names>Ilyasova</given-names>
            <surname>NYu</surname>
          </string-name>
          .
          <article-title>Methods for digital analysis of human vascular system. Literature review</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2013</year>
          ;
          <volume>37</volume>
          (
          <issue>4</issue>
          ):
          <fpage>517</fpage>
          -
          <lpage>541</lpage>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Simcher</surname>
          </string-name>
          , VM.
          <article-title>Methods of multivariate statistical analysis</article-title>
          .
          <source>M.: Finance and statistics</source>
          ,
          <year>2008</year>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Mookiaha</surname>
            <given-names>MRK</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Acharyaa</surname>
            <given-names>UR</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lima</surname>
            <given-names>CM</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Petznickb</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jasjit</surname>
            <given-names>S</given-names>
          </string-name>
          .
          <article-title>Data mining technique for automated diagnosis of glaucoma using higher order spectra and wavelet energy features</article-title>
          .
          <source>Knowledge-Based Systems</source>
          ,
          <year>2012</year>
          ;
          <volume>33</volume>
          :
          <fpage>73</fpage>
          -
          <lpage>82</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Ilyasova</surname>
            <given-names>NYu</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kupriyanov</surname>
            <given-names>AV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Paringer</surname>
            <given-names>RA</given-names>
          </string-name>
          .
          <article-title>Formation features for improving the quality of medical diagnosis based on the discriminant analysis methods</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2014</year>
          ;
          <volume>38</volume>
          (
          <issue>4</issue>
          ):
          <fpage>851</fpage>
          -
          <lpage>855</lpage>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <given-names>Ilyasova</given-names>
            <surname>NYu</surname>
          </string-name>
          .
          <article-title>Estimating the geometric features of a 3D vascular structure</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2014</year>
          ;
          <volume>38</volume>
          (
          <issue>3</issue>
          ):
          <fpage>529</fpage>
          -
          <lpage>538</lpage>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <given-names>Ilyasova</given-names>
            <surname>NYu</surname>
          </string-name>
          .
          <article-title>Methods to Evaluate the Three-Dimensional Features of Blood Vessels</article-title>
          .
          <source>Optical Memory and Neural Networks (Information Optics)</source>
          ,
          <year>2015</year>
          ;
          <volume>24</volume>
          (
          <issue>1</issue>
          ):
          <fpage>36</fpage>
          -
          <lpage>41</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Ilyasova</surname>
            <given-names>NYu</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kupriyanov</surname>
            <given-names>AV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khramov</surname>
            <given-names>AG</given-names>
          </string-name>
          .
          <article-title>Information technologies of image analysis in medical diagnostics</article-title>
          .
          <source>M: Radio and Communication</source>
          ,
          <year>2012</year>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Ilyasova</surname>
            <given-names>N.</given-names>
          </string-name>
          <article-title>Computer Systems for Geometrical Analysis of Blood Vessels Diagnostic Images</article-title>
          .
          <source>Optical Memory and Neural Networks (Information Optics)</source>
          ,
          <year>2014</year>
          ;
          <volume>23</volume>
          (
          <issue>4</issue>
          ):
          <fpage>278</fpage>
          -
          <lpage>286</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Ilyasova</surname>
            <given-names>NYu</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kupriyanov</surname>
            <given-names>AV</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ananin</surname>
            <given-names>MA</given-names>
          </string-name>
          .
          <article-title>Measurement of the biomechanical vessels parameter for the diagnostics of the early stages of the retina vascular pathology</article-title>
          .
          <source>Computer Optics</source>
          ,
          <year>2005</year>
          ;
          <volume>27</volume>
          :
          <fpage>165</fpage>
          -
          <lpage>169</lpage>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Fukunaga</surname>
            <given-names>K.</given-names>
          </string-name>
          <article-title>Introduction to statistical pattern recognition</article-title>
          . New York and London: Academic Press,
          <year>1972</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Kim</surname>
            <given-names>JA</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Myuller</surname>
            <given-names>ChU</given-names>
          </string-name>
          , Klekka WR.
          <article-title>Factor, discriminant and cluster analysis</article-title>
          .
          <source>M.: Finance and Statistics</source>
          ,
          <year>1989</year>
          . [in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          13.
          <string-name>
            <surname>Moshetova</surname>
            <given-names>LK</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Yushchuk</surname>
            <given-names>ND</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tsyganov</surname>
            <given-names>DI</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sister</surname>
            <given-names>VG</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Branchevsky</surname>
            <given-names>SL</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ilyasova</surname>
            <given-names>NYu</given-names>
          </string-name>
          , Pavlova Y.
          <article-title>Method for digital image processing fundus</article-title>
          .
          <source>Information letter</source>
          ,
          <year>2004</year>
          ;
          <fpage>1</fpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>