<!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>October</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>A STUDY ON PERFORMANCE ASSESSMENT OF ESSENTIAL CLUSTERING ALGORITHMS FOR THE INTERACTIVE VISUAL ANALYSIS TOOLKIT INVEX</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>M.A. Titov</string-name>
          <email>mikhail.titov@cern.ch</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff4">4</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>M.A. Grigorieva</string-name>
          <email>maria.grigorieva@cern.ch</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff4">4</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>A.A. Alekseev</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff4">4</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>N.A. Belov</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>T.P. Galkin</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>D.V. Grin</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>T.A. Korchuganova</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>S.A. Zhumatiy</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Lomonosov Moscow State University</institution>
          ,
          <addr-line>Leninskie Gory, 1, Moscow, 119991</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Mikhail Titov</institution>
          ,
          <addr-line>Maria Grigorieva, Aleksandr Alekseev, Nikita Belov, Timofei Galkin, Dmitry Grin, Tatiana Korchuganova, Sergey Zhumatiy</addr-line>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>National Research Center “Kurchatov Institute”</institution>
          ,
          <addr-line>Akademika Kurchatova pl., 1, Moscow, 123182</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff3">
          <label>3</label>
          <institution>National Research Nuclear University “MEPhI”</institution>
          ,
          <addr-line>Kashirskoe shosse, 31, Moscow, 115409</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff4">
          <label>4</label>
          <institution>Plekhanov Russian University of Economics</institution>
          ,
          <addr-line>Stremyanny lane, 36, Moscow, 117997</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2019</year>
      </pub-date>
      <volume>4</volume>
      <issue>2019</issue>
      <fpage>434</fpage>
      <lpage>438</lpage>
      <abstract>
        <p>Interactive visual analysis tools bring the ability of the real-time discovery of knowledge in large and complex datasets using visual analytics. It involves multiple iterations of data processing using various data handling approaches and the efficiency of the whole chain of the analysis process depends on the performance of chosen techniques and related implementations, as well as the quality of applied methods. Stages, where data processing includes intellectual handling (i.e., data mining and machine learning), which are the most resource-intensive, require a distinct attention for evaluation of different approaches. Clustering is one such machine learning technique that is commonly used to discover groups of data objects for further analysis. This work is focused on evaluation of clustering algorithms within the interactive visual analysis toolkit InVEx (Interactive Visual Explorer). InVEx represents a visual analytics approach aimed at cluster analysis and in-depth study of implicit correlations between multidimensional data objects. It is originally designed to enhance the analysis of computing metadata of the ATLAS experiment at the LHC for operational needs, but it also provides the same capabilities for other domains to analyze large amounts of multidimensional data. The experiments and evaluation processes are carried out using operational data from the supercomputer at the Lomonosov Moscow State University. These processes include benchmark tests to assess the relative performance between chosen clustering algorithms and corresponding metrics to assess the quality of produced clusters. Obtained results will be used as guidelines in assisting users in a process of visual analysis using InVEx.</p>
      </abstract>
      <kwd-group>
        <kwd>visual analytics</kwd>
        <kwd>clustering</kwd>
        <kwd>benchmarks</kwd>
        <kwd>InVEx</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>The interactivity as an integral part of the visual analysis, besides its essential objective of the
knowledge discovery in real-time (with a focus on large and complex datasets) while doing the
analysis, also brings challenges of keeping up performance and efficiency. Every stage of the process
of interactivity should be evaluated to estimate the overall performance metrics. Thus, this work is
emphasized to the one of crucial processes of the analysis - intellectual data processing, e.g.,
application of machine learning (ML) algorithms for clustering.</p>
      <p>Clustering algorithms, used in the developed visual analysis toolkit InVEx, require a
significant effort to select the most relevant object attributes for different chosen algorithms [as for
user/analyst], and to adjust primary parameters of provided algorithms [as for developers]. In addition,
the outcome of different clustering algorithms should be compared between each other and being
evaluated by the quality of clustering based on expert reviews.</p>
      <p>This paper brings the assessment of the clustering algorithms usage in the InVEx toolkit to
enhance the user experience and to improve the quality of the analysis process.</p>
      <sec id="sec-1-1">
        <title>1.1 InVEx overview</title>
        <p>
          InVEx stands for the Interactive Visual Explorer toolkit [
          <xref ref-type="bibr" rid="ref1 ref2">1,2</xref>
          ]. It provides advanced interactive
data visualization tools, which are used for the analysis of large volumes of multidimensional data
with its core process as an interactive visual clustering. Its development has been started for the
ATLAS Distributed Computing project (the ATLAS experiment [
          <xref ref-type="bibr" rid="ref3">3</xref>
          ] at the Large Hadron Collider) to
enhance the analysis of computing metadata [
          <xref ref-type="bibr" rid="ref2">2</xref>
          ]. Large amount of ATLAS ProdSys2/PanDA metadata
provides means to test and prove the efficiency of applied technologies and methods. The ATLAS
Production System (ProdSys2) [
          <xref ref-type="bibr" rid="ref4">4</xref>
          ], in conjunction with the workload management system namely the
Production and Distributed Analysis system (PanDA) [
          <xref ref-type="bibr" rid="ref5">5</xref>
          ], represents a complex set of computing
components that are responsible for organizing, planning, starting and executing distributed computing
tasks and jobs. Initial integration of InVEx with PanDA includes the direct access to the information
about computing jobs from PanDA’s monitoring system.
        </p>
        <p>The stack of technologies for InVEx includes: Python-based Django web framework;
crossbrowser JavaScript library Three.js to create and display animated 3D computer graphics in a web
browser (uses WebGL); Python libraries for data handling and analysis such as Pandas (data
manipulation and analysis), SciPy (scientific computing and technical computing, as well includes
clustering algorithms: hierarchical clustering, vector quantization, K-means), Scikit-learn (ML library,
it features various classification, regression and clustering algorithms), Kmodes (clustering for
categorical data, implementations of k-modes and k-prototypes clustering algorithms), Intel Data
Analytics Acceleration Library / daal4py (optimized algorithmic building blocks for data analysis
stages), Prince (factor analysis that aims to find independent latent variables).</p>
      </sec>
      <sec id="sec-1-2">
        <title>1.2 Lomonosov-2 supercomputer overview</title>
        <p>
          The current study uses log data about computing jobs gathered from the Lomonosov-2
supercomputer. This supercomputer is designed by the T-Platforms company and installed at the
Lomonosov Moscow State University (MSU) [
          <xref ref-type="bibr" rid="ref6 ref7">6,7</xref>
          ] (its rank is #93 in the TOP500 list1). It is
characterized by the Intel Xeon/FDR InfiniBand cluster, accelerated with NVidia Tesla K40s and
Tesla P100 GPUs, and with overall 1696 nodes (Intel Haswell-EP E5-2697v3, 2.6GHz, 14 cores and
Intel Xeon Gold 6126 2.6GHz, 12 cores) with 64/96 GB of memory per node. Theoretical peak
performance is 4.946 petaFLOPS. (For more details please follow the references [
          <xref ref-type="bibr" rid="ref6 ref7">6,7</xref>
          ].)
1 Position of the Lomonosov-2 supercomputer within the TOP500 ranking of supercomputers. Available at:
https://www.top500.org/system/178444 (accessed on 20.11.2019)
        </p>
      </sec>
    </sec>
    <sec id="sec-2">
      <title>2. Methods and techniques</title>
      <sec id="sec-2-1">
        <title>2.1 Clustering overview</title>
        <p>
          Clustering is a ML technique, that is aimed at grouping similar objects into unlabeled groups
called clusters (unsupervised learning). The process of clustering in InVEx is implemented in two
stages: i) Level-of-Detail (LoD) generator [
          <xref ref-type="bibr" rid="ref2">2</xref>
          ] that brings the initial (optional) grouping to reduce the
amount of data objects provided to the user (algorithms that are used: scikit-learn/MiniBatchKMeans,
daal4py/KMeans, kmodes/KPrototypes); ii) cluster analysis, which is a core process to analyze data
objects similarities (algorithms that are used: same as for LoD, as well as scikit-learn/KMeans,
scikitlearn/DBSCAN, scipy/Hierarchical).
        </p>
      </sec>
      <sec id="sec-2-2">
        <title>2.2 Clustering validation measures</title>
        <p>Validation measures (i.e., quality metrics) are classified as internal and external. Internal
measures reflect compactness, connectedness and separation of the cluster partitions. The following
metrics were chosen: Silhouette coefficient (ranges from 1 to 1, where a high value indicates that the
object is well matched to its own cluster); Calinski-Harabaz Index (the maximum value for index
indicates a suitable partition for the data set); Davies-Bouldin Index (closer to 0 is better, it computes
the ratio between the within cluster distances and the between cluster distances). External measures,
which provide comparison of the identified clusters to external preset labels, are represented in this
paper by the following metric - Adjusted Rand Index - a function that measures the similarity of the
two assignments (given the knowledge of the ground truth class assignments and clustering algorithm
assignments). (Further extension of using external metrics considers Fowlkes-Mallows score that is
defined as the geometric mean of the pairwise precision and recall.)</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>3. Experiments</title>
      <sec id="sec-3-1">
        <title>3.1 Data pre-processing</title>
        <p>Gathered data for experiments required initial transformations because of the format of data
objects attributes. Log data from the Lomonosov-2 supercomputer was collected from the period of
300 days (from June 2018 to March 2019): 245K records with 12 attributes (user ID, execution time
duration, number of allocated nodes, CPU load during the job execution per user, GPU load during the
job execution per user, number of executed instructions per second, etc.). Almost all of the attributes
were of the categorical type (nominal and ordinal data) and most of them are with such nominal values
as: “none”, “low”, “average”, “high”. Thus, the technique of dimensionality reduction was applied,
that is used to map the data record to a lower-dimensionality space.</p>
        <p>The process of the data transformation was the following: i) apply multiple correspondence
analysis (MCA) to the dataset with all 11 categorical attributes to represent data objects in a
multidimensional Euclidean space with 5 dimensions; ii) apply principal component analysis (PCA) to
the dataset with 5 attributes from the previous step and 1 non-categorical attribute from the original
dataset. The outcome of this transformation is a dataset with 5 attributes (per record), which will be
used for clustering in the benchmarks.</p>
      </sec>
      <sec id="sec-3-2">
        <title>3.2 Benchmarks and quality metrics</title>
        <p>There are several essential algorithms (partitioning-based and density-based methods) that
were chosen for performance and quality evaluation. Figure 1 presents benchmarks for KMeans
algorithms of different implementations. Scikit-learn/KMeans was significantly inferior to other
implementations, and especially for input data over 150K records that took up to hundreds of seconds
to be executed, thus its benchmark was not included into the figure. Figure 1 shows that
daal4py/KMeans outperforms scikit-learn/MiniBatchKMeans. Quality metrics (internal measures) for
KMeans implementations gives better results compared to MiniBatchKMeans (Table 1). Table 2
shows that scikit-learn/KMeans and daal4py/KMeans give very close to each other labeling of
clusters, which makes daal4py implementation preferred to be used in the next versions of InVEx for
performance improvements and without loss of quality. Benchmarks for two density-based algorithms
are presented in Figure 2, which shows greater performance for HDBSCAN. These algorithms
produce very close results in terms of quality metrics (Tables 1,2), thus HDBSCAN is more suitable
for use (which is also more robust to parameter selection in comparison to DBSCAN).</p>
        <p>(a) (b)
Figure 1. Performance comparison for partitioning-based clustering algorithms with different number
of records as input data and different number of outcome clusters: a) scikit-learn/MiniBatchKMeans;
b) daal4py/KMeans</p>
        <p>(a) (b)
Figure 2. Performance comparison for density-based clustering algorithms with different number of
records as input data: a) HDBSCAN; b) OPTICS</p>
        <p>Silhouette Score
HDBSCAN
0.61
0.58
0.68
0.58
92853.7
88390.1
32503.6
211.5
0.86
0.94
1.08
1.37
HDBSCAN vs. OPTICS
0.55
0.95
0.84</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Conclusion</title>
      <p>The interactive visual analysis toolkit InVEx represents a visual analytics approach aimed at
cluster analysis and in-depth study of implicit correlations between multidimensional data objects and
object parameters interdependencies. Its capabilities were applied in analysis of log data from the
Lomonosov-2 supercomputer, which also were used to conduct experiments on performance
estimation for InVEx clustering algorithms.</p>
      <p>Experiments outcome is the process of evaluation essential clustering algorithms. Benchmark
tests assess the relative performance between chosen algorithms and corresponding metrics, and the
quality of produced clusters. Obtained results and the approach itself will be integrated into InVEx and
will be used as guidelines in assisting users in a process of visual analysis (as well as a self-adjustment
mechanism to configure initial parameters for clustering algorithms).</p>
    </sec>
    <sec id="sec-5">
      <title>5. Acknowledgement</title>
      <p>Many thanks to all members of the InVEx team and colleagues from the Research Computing
Center (RCC) of MSU for providing experimental data and for the continued support. This work was
financially supported by the Russian Science Foundation (grant No.18-71-10003).</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1] InVEx project, “InVEx” [software],
          <year>2019</year>
          . Available at: https://github.com/PanDAWMS/InVEx (accessed
          <source>on 20.11</source>
          .
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Grigorieva</surname>
            <given-names>M.A.</given-names>
          </string-name>
          et al.
          <article-title>Evaluation of the Level-of-Detail Generator for Visual Analysis of the ATLAS Computing Metadata //</article-title>
          <source>Lobachevskii Journal of Mathematics</source>
          , vol.
          <volume>40</volume>
          , no.
          <issue>11</issue>
          , pp.
          <fpage>1788</fpage>
          --
          <lpage>1798</lpage>
          (
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>ATLAS</given-names>
            <surname>Collaboration</surname>
          </string-name>
          .
          <source>The ATLAS Experiment at the CERN Large Hadron Collider // JINST</source>
          , vol.
          <volume>3</volume>
          ,
          <issue>S08003</issue>
          (
          <year>2008</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>Barreiro F.H</surname>
          </string-name>
          . et al.
          <article-title>The ATLAS Production System Evolution: New Data Processing and Analysis Paradigm for the LHC Run2 and High-</article-title>
          <source>Luminosity // J. Phys.: Conf. Ser.</source>
          , vol.
          <volume>898</volume>
          , no.
          <issue>5</issue>
          ,
          <issue>052016</issue>
          (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Barreiro F.H</surname>
          </string-name>
          . et al.
          <article-title>PanDA for ATLAS distributed computing in the next decade // J</article-title>
          . Phys.: Conf. Ser., vol.
          <volume>898</volume>
          , no.
          <issue>5</issue>
          ,
          <issue>052002</issue>
          (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Voevodin</surname>
            <given-names>V.V.</given-names>
          </string-name>
          et al. Supercomputer Lomonosov-
          <volume>2</volume>
          :
          <string-name>
            <given-names>Large</given-names>
            <surname>Scale</surname>
          </string-name>
          ,
          <article-title>Deep Monitoring and Fine Analytics for the User Community // Supercomputing Frontiers and Innovations</article-title>
          , vol.
          <volume>6</volume>
          , no.
          <issue>2</issue>
          , pp.
          <fpage>4</fpage>
          --
          <lpage>11</lpage>
          (
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <surname>Leonenkov</surname>
            <given-names>S.</given-names>
          </string-name>
          , Zhumatiy S. Supercomputer Efficiency: Complex Approach Inspired by Lomonosov-2 History Evaluation // RuSCDays 2018:
          <article-title>Supercomputing</article-title>
          .
          <source>Communications in Computer and Information Science</source>
          , vol.
          <volume>965</volume>
          . Springer, Cham (
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>