<!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>Invasion Detection Model using Two-Stage Criterion of Detection of Network Anomalies</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Borys Grinchenko Kyiv University</institution>
          ,
          <country country="UA">Ukraine</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Kharkiv National University of Radio Electronics</institution>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <fpage>0000</fpage>
      <lpage>0002</lpage>
      <abstract>
        <p>The article considers the methods that affect the operation of the intrusion detection system. Using the “disconnection task,” a two-stage criterion for detecting anomalies in computer networks has been formed, which provides an analysis of network infrastructure characteristics and their identification with specific computer attacks and provides the ability to respond to possible attacks in real-time.</p>
      </abstract>
      <kwd-group>
        <kwd>Intrusion Detection Systems</kwd>
        <kwd>Anomalies</kwd>
        <kwd>Signature Analysis Method</kwd>
        <kwd>Statistical Analysis Method</kwd>
        <kwd>Neural Networks</kwd>
        <kwd>Invasion Detection Model</kwd>
        <kwd>Attack</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Based on the annual rate of increase in the number of incidents in cyberspace, we can
conclude that it is necessary to independently include in the integrated information
security system critical infrastructure of automated means of detecting computer
attacks and other dangerous divisions, including threats of man-made and natural
(random) nature.</p>
      <p>
        Modern Intrusion Detection Systems (IDS) to achieve information security goals
constantly monitor the functioning of hardware and software platforms of the network
infrastructure and record several quantitative and qualitative indicators of their work:
{X1(t), …, XN(t)}
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
where t is the time at which the measurement of the indicator Xk(t).
      </p>
      <p>Such indicators, in particular in the NIDES system [1], are:
 Ability to use the CPU separately by the system and the user.
 Time to complete the process.
 Total amount of memory used during the execution of the process and its
maximum size during execution.
 Number of open files at runtime; the number of page failures.
 Amount of information read from the disk.
 Number of I/O characters during application execution.
 Username changed during the execution of the application.
 Start time of the application.
 Number of signals received during the execution of the application.
 Application was running on the remote station and the name of this station
 Name of the application that was used on the remote station.
 Application was running on the local station and the name of this station, the name
of the application used on the local station, etc.
2</p>
    </sec>
    <sec id="sec-2">
      <title>Problem Formulation</title>
      <p>With this in mind, this article aims to build a model of intrusion detection system that
would detect network anomalies in real-time and with the least error, as well as find
answers to the following questions [2]:
 What happened online?
 What was attacked and how dangerous is the attack?
 When and where did the attack start?
 Who is the malefactor?
 How and as a result of what there was the invasion?
Accordingly, the complexity of the attacking party’s tasks to create cryptographic
software tools for attacks is significantly simplified, and the effectiveness of solving
problems of recovery (decryption) of cryptosystem data, which are encrypted with
software tools for attacks, has an objective tendency to constant decline from year to
year. At the same time, there are such groups of threats related to reliability [3–5]:
 Cryptographic primitive programming errors.
 Errors in the use of cryptographic primitives in crypto providers.
 Errors in transferring parameters to the crypto provider and returning processing
results.
 Errors in the failure of the hardware platform.
 Accidental and intentional violations of the integrity of programs and data of
crypto providers.</p>
      <p>These factors may be prerequisites for the formation of hidden channels of leakage of
information about the operation of cryptographic software tools for attacks. Consider
two options for attacks on the implementation of software tools for attacks that do not
involve the presence of critical information about their work.
3</p>
    </sec>
    <sec id="sec-3">
      <title>Crypto Attack Mechanics</title>
      <p>Premise 1. There is a certain cryptographic software tool for attacks (cryptosystem)
with secret keys based on a stable block symmetric cryptographic algorithm, which
provides encryption in output feedback (OFB) mode or otherwise tampering with
feedback [6]. It is necessary: to build an attack on the implementation of the specified
cryptosystem to ensure the possibility of partial or complete recovery of information
(decryption).</p>
      <p>This mode of operation of the crypto algorithm is widely used to build fully
connected communication networks, in which each subscriber must send a message to
any other subscriber of this network.</p>
      <p>
        This mode is characterized by the presence of an initialization vector IV is random
number, which provides at an acceptable level the probability of not overlapping the
cipher in the case of using one key K for some time [6]. The idea of the attack on this
cryptosystem is to implement a fictitious random data generator of the crypto
algorithm, which is included instead of the real pseudorandom number generator, provides
a hidden repetition of initialization vectors or keys. As a result, some messages will
be encrypted in the same way:
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
where  ̃1, … ,  ̃  are open messages;  ̃ ,  = 0̅̅.̅2̅̅̅̅−̅̅̅1̅ is a sequence of gamma
characters formed from a single initialization vector and an algorithm key;  ̃1, … ,  ̃  are
corresponding encrypted messages.
      </p>
      <p>In this case, depending on the number of equally encrypted messages   and
redundancy of the source texts of the messages, the corresponding encrypted messages
can be partially or completely decrypted [7]. To repeat the initialization vectors IV the
random bit generator should be used as follows:</p>
      <p>
        IV   1, 2,..., k   1,2,...,b
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
where  1,  2, … ,   is the sequence of random bits;  (… ) some unambiguous
function that expands the random sequence to a given size b;  1,  2, … ,   is coordinates
of the initialization vector.
      </p>
      <p>To solve the problem, we first calculate the allowable number of random bits to
provide the required number of repetitions of equally encrypted messages. Let in the
network M subscribers, each of whom sends on average daily  messages for each
connection direction. With T is public key validity period (days). Then the average
number of messages sent  is the value:
 S~1  T1  Г j</p>
      <p>~ ~

 ... ... ...
S~m j  T~m j  Г~j</p>
      <p> M </p>
      <p>N   2    T
During the specified period, the average number of repetitions  each of 2 the values
of vectors of random bits 〈 1,  2, … ,   〉 in the case of their even distribution will be:
2k</p>
      <p>.</p>
      <p>  M  
k  log 2  2     T  log 2 R </p>
      <p>
         
Based on equation (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) we obtain an estimate of the allowable number of random bits:
Using inequality (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) in each case of the attacked network, it is possible to calculate
the allowable number of random bits to provide an average static number of  equally
encrypted messages. Of course, random data used for cryptographic transformations
are subject to statistical testing, but the length of the vector is quite small. IV (block
length for most standard cryptographic algorithms is 64 or 128 bits) significantly
limits the ability to apply statistical criteria. To check the uniformity of the
distribution of sequences of this length, the criteria of frequencies of signs and bigram are
mainly used [8], which puts forward the appropriate requirements for the expansion
function. Therefore, it is clear that to ensure the equal occurrence of bigrams in the
binary vector 〈 1,  2, … ,   〉 appropriate value k choose a multiple of four (there are
four options for bigram 00, 01, 10, 00).
      </p>
      <p>For example, in the case of a special network that includes M = 10 subscribers with
an average intensity of sending messages during the validity of one network key
 ∙  = 10 to obtain the average number of repetitions of encryption  = 2 it is
necessary to have several random bits no more  ≤ 7.</p>
      <p>To hide statistical dependencies in vector IV and obtain the required number of bits
in the corresponding expression as a function , you should use a stable hash
function, such as MD5, which is part of the standard set of Windows algorithms, or
similar. Thus from 128 bits of the digest which we receive using a hash function, we will
use the necessary quantity.</p>
      <p>Premise 2. Suppose that there is a software tool for implementing attacks
(cryptosystem), built on a stable symmetric block cryptographic encryption algorithm in
OFB mode using public key distribution based on the El Gamal protocol [6]. The task
is to organize an attack on the implementation of the cryptosystem in such a way as to
ensure complete decryption.</p>
      <p>
        Initially, the operation of the cryptosystem in the El Gamal protocol using a cyclic
element g of some field and the secret key of the asymmetric algorithm x is formed
public key y:
y  g x mod p
(
        <xref ref-type="bibr" rid="ref7">7</xref>
        )
where p is some large prime number, 1024 bits long.
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref6">6</xref>
        )
To encrypt the session key of the symmetric algorithm, a random number k is
generated, the value is calculated:
and the secret key of the symmetric algorithm is encrypted K:
      </p>
      <p>
        y1  g k mod p .
  K  yk mod p .
(
        <xref ref-type="bibr" rid="ref8">8</xref>
        )
(
        <xref ref-type="bibr" rid="ref9">9</xref>
        )
Pair (δ,y1) together with the encrypted message is transmitted to the owner of the
secret key x, which is based on the obtained pair (δ,y1) calculates the secret key of a
symmetric algorithm K:
      </p>
      <p>
        K    y1x mod p .
(
        <xref ref-type="bibr" rid="ref10">10</xref>
        )
We organize an attack on the specified protocol by substituting a random number k
with a pseudo-random one so that the testing system does not reveal this fact. Let
k{t., t=1,2,…,}, at the same time power  the set of valid values of the
pseudorandom number k choose a fairly large, but less than the performance of a specialized
computer system for the search of keys for a reasonable period.
      </p>
      <p>
        Calculation of possible key variants K у in this case we carry out using
International Planned Parenthood Federation (IPPF) proceeding from the equation (
        <xref ref-type="bibr" rid="ref10">10</xref>
        ) based on
a set of admissible values of pseudorandom number {t} [9].
      </p>
      <p>As a result, it should be noted that access to information of the “true” random data
generator and methods of its substitution in each case depends on the hardware and
software platforms on which the IPPF, the peculiarities of its implementation, the
presence of an automated system of protection against unauthorized interference. her
work, etc. To create a fictitious random number generator, it is possible to use a
sequence of Mersenne numbers   = 2 − 1, if  is a simple number,   also simple.
4</p>
    </sec>
    <sec id="sec-4">
      <title>Methods of Analysis</title>
      <p>It is known that the efficiency of the IDS system significantly depends on the applied
method of analysis of the original data. Currently, the following main methods are
distinguished [1]:
1. Signature methods of analysis.
2. Statistical methods of analysis.</p>
      <p>3. Hybrid methods of analysis with the function of self-learning.
4.1</p>
      <sec id="sec-4-1">
        <title>Signature Method</title>
        <p>Signature methods of analysis are based on the fact that most attacks and their
scenarios are generally known. In this approach, the signatures of invasions determine the
characteristics and conditions of the objects, the occurrence of events that are signs of
attempted attacks (invasion), and their relationship. Of course, signature methods of
analysis use intrusion signature databases that are supported by the security system. In
this case, the order (sequence) of actions performed or initiated by the user or
information process (program) is compared with known signatures. A sign of an attempted
security breach may be partial compliance of the sequence of events with the
signature.</p>
        <p>Anti-virus scanners that work with a database of virus signatures are typical
representatives which implement this idea. Their advantage is a fairly high speed of
analysis. The effectiveness of signature methods of IDS can be increased through the use of
methods of artificial intelligence. The advantages of statistical methods of analysis in
IDS include:
 No need for a large amount of memory to store controlled variables.
 Ease of detecting deviations in the data characterizing the behavior of users and
processes.
 Ability to analyze quantitative and qualitative data of different nature of origin as a
parameter in the analysis.
 Disadvantages of statistical methods include difficulties in generating statistics of
the normal behavior of users and processes.
4.2</p>
      </sec>
      <sec id="sec-4-2">
        <title>Hybrid Method</title>
        <p>Hybrid methods of analysis with the function of self-learning include:
 Teaching methods for the classification of examples.
 Neural networks.
 Genetic algorithms.</p>
        <p>To achieve the goal set in the work, it is advisable to use statistical methods exactly.
This allows, using the well-known “problem of disorder” [10]: to form a two-stage
criterion for the detection of anomalies in computer networks, to overcome the
connection of statistical methods to the model of usual (normal) users’ behavior.</p>
        <p>
          Step 1. Let’s suppose that two complex hypotheses are considered concerning the
random sequence  = {  ,   ∈  ,  = ̅1̅,̅̅̅} under analysis: H0: the sequence X is
stationary with a single probability distribution function, H1: the sequence X is a
concatenation (the result of “gluing”) of two stationary random sequences with different
distribution functions:
 =  1|| 2,
(
          <xref ref-type="bibr" rid="ref11">11</xref>
          )
where  1 = {  ,  = ̅1̅,̅̅̅̅∗},  2 = {  ,  = ̅̅∗̅̅+̅̅̅̅1̅,̅̅̅},  ∗ = [ ], 0 &lt;  &lt; 1
        </p>
        <p>
          It is necessary to estimate the point of “gluing” n*. It is believed that the sequences
X1 and X2 differ in one of the two-dimensional distribution functions, namely, the
probability distribution of the vector (  ,   +2):
 ( 0,  1) =  {  ≤  0,   +2 ≤  1}
(
          <xref ref-type="bibr" rid="ref12">12</xref>
          )
|| 1( )−  2( )|| ≥  &gt; 0,
where ||…|| is normal sup-norm.
        </p>
        <p>It is known that the distribution function of a finite-dimensional random vector can
be approximated uniformly with any accuracy to the probability distribution function
of a random vector with a finite number of values. It follows that if we give the set R
is a combination of a sufficiently large number of domains { ,  = ̅1̅,̅̅}, that do not
intersect   ∩   = ∅ for i ≠ j, then the vector (  ,   +2) can be approximated by a
distribution vector with a finite number of values [11].</p>
        <p>Therefore, if we are to enter new random sequences</p>
        <p>=  (  ∈   ,   +2 ∈   ), де 1 ≤  ≤  , 1 ≤  ≤  ,
where I(A) is an indicator of the set A, then at least in one of them there is a change in
mathematical expectation.</p>
        <p>Therefore, if we use an algorithm that detects a change in mathematical
expectation, the same algorithm will detect a change in the distribution function. This fact
allowed the work to be limited to the development of only one, a basic algorithm that
can detect changes in the mathematical expectation. To do this, to identify moments
of “disorders” a family of statistics of the form is proposed:
to the moment  1∗ =  ∗ − 2 including is equal to F1(u), and at  ≥  2∗ =  ∗ + 1 is
equal to F2(u), and
(13)
(14)
(15)
(16)
  ( ,  ) = [ (1 −


 
)] ∙ [1 ∑


 =1   −</p>
        <p>1
 −</p>
        <p>∑ = +1   ],
where 0 &lt;  ≤ 1, 1 ≤  ≤  − 1,  = {  ,  = ̅1̅,̅̅̅}} is the sequence under
study.</p>
        <p>The given family of statistics in the case of a fixed n is a generalized variant of
Kolmogorov-Smirnov’s statistics, which is used to test hypotheses of coincidence or
difference of distribution functions in two samples. The paper also proves that the
statistics in the case of  = 1 for</p>
        <p>→ ∞ and maintaining the ratio between the
volumes of “glued” implementations minimizes the maximum possible probability of
error estimating the moment of “disorder” (minimax in order):
 {</p>
        <p>max
1≤ ≤ −1</p>
        <p>
          √ |  ( , 1)| &gt;  (
          <xref ref-type="bibr" rid="ref1">1</xref>
          ) } →
2 ∑∞=1(−1) +1 exp (−2 2 (С(
          <xref ref-type="bibr" rid="ref1">1</xref>
          ) 2
        </p>
        <p>
          ) ) ≡  ( (
          <xref ref-type="bibr" rid="ref1">1</xref>
          )),
 ∗
where the parameter  ∗ is the standard deviation, С(
          <xref ref-type="bibr" rid="ref1">1</xref>
          ) is the limit of the criterion, the
excess of which will be perceived as the occurrence of “disorder,” the value of  ∗, for
which it occurred, is the desired moment of “disorder.”
        </p>
        <p>
          Step 2. Having fixed the level of probability  “false alarm” about the disorder,
during the statistical processing of real data we determine the level of the threshold of
the first level C(
          <xref ref-type="bibr" rid="ref1">1</xref>
          ):
where  ̅∗ is the estimate of the parameter  ∗ (standard deviation), and N is the sample
volume of the sequence under study [12].
        </p>
        <p>In the case of hypothesis Н1 on the representation of the original sequence of
measurements in the form of the concatenation of several stationary random
sequences with different probability distribution functions, we apply the criterion of the
allowable number of triggers (“disorders”) Z. That is, we assume that the following
equation takes place:</p>
        <p>1 +  2 + … +   =  , де 2 ≤  &lt;  − 2.</p>
        <p>For a small number of sequences, the Chebyshev inequality can be used to determine
the boundary of the criterion if the random variable Z has a mathematical expectation
 and a standard deviation  for a given  &gt; 0:
 {| −  | ≥  } ≤
 {| −  | ≥  } ≤
 2
 2
4
9 2
If Z is a random variable with a single-mode probability distribution with a
mathematical expectation  and a standard deviation of 0 &lt;  &lt; ∞, then for any  &gt; √8⁄3 ≈
1.63299, there is a Vysochansky-Petunin inequality, which improves the estimation
of the deviation probability:
In particular, for a typical deviation of 3σ (three sigmas), the significance level of the
criterion is the probability of “false alarm” is calculated as [13]:
(17)
(18)
(19)
(20)
(21)
(22)
(23)
(24)
4
 = 9 2 ≈ 0.0494.</p>
        <p>=  1 +  2 + … +   .</p>
        <p>
          A more accurate result can be obtained for the case of a random variable Z, which is
the sum of a sufficiently large number of independent random variables (
→ ∞):
According to the integral limit theorem, the probability of deviation can be
approximated by the normal probability distribution function  (
          <xref ref-type="bibr" rid="ref1">0,1</xref>
          ):
        </p>
        <p>{|
 − | ≥  1− } =
1
√2
∫−∞1− exp(−  2) .</p>
        <p>
          2
Based on (
          <xref ref-type="bibr" rid="ref12">12</xref>
          ), we calculate the boundary of the second level criterion:
 (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ) =  +  1−  .
        </p>
        <p>Thus, the algorithm for calculating the criterion includes the following steps (Fig. 1).</p>
        <p>It is possible to see that the complexity of the algorithm that implements the proposed
criterion is estimated by the value of O(N). Thus, its implementation on the principle
of the “sliding window” will not cause a significant load on the computer system.</p>
        <p>The most difficult issue for the implementation of the method of detecting
anomalies, as in other statistical methods, is the formation of regular behavior of the network
infrastructure and its users.
5</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>Conclusions</title>
      <p>In this regard, of particular interest are modern technologies that increase the
efficiency of classical decryption methods [6], which using vulnerabilities in a particular
cryptographic algorithm and/or protocol provide recovery of the secret key used to
create encrypted messages, and/or disclosing the original value of the encrypted
information. Thus recovery of the initial information can be full or partial therefore
speak about full or partial decryption.</p>
      <p>A two-level statistical criterion for detecting anomalies based on measurements of
network infrastructure characteristics is proposed, which provides the possibility of
their further analysis to identify them with certain computer attacks. The application
of the criterion does not require significant computing resources and provides an
opportunity to respond to possible attacks in almost real-time.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Court</surname>
            ,
            <given-names>S. S.: Intruder</given-names>
          </string-name>
          <string-name>
            <surname>Detection Methods</surname>
          </string-name>
          (
          <year>2004</year>
          ). [Publication in Ukrainian]
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Kostrov</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          :
          <article-title>Intrusion detection systems</article-title>
          .
          <source>Byte</source>
          <volume>8</volume>
          (
          <issue>49</issue>
          ):
          <fpage>14</fpage>
          -
          <lpage>21</lpage>
          (
          <year>2002</year>
          ). [Publication in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Zubov</surname>
          </string-name>
          , A. Y.:
          <article-title>Cryptographic Methods of Information Protection</article-title>
          . Perfect
          <string-name>
            <surname>Ciphers</surname>
          </string-name>
          (
          <year>2005</year>
          ). [Publication in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Duhin</surname>
            ,
            <given-names>A. A.</given-names>
          </string-name>
          :
          <article-title>Information theory (</article-title>
          <year>2007</year>
          ). [Publication in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Gorbenko</surname>
            ,
            <given-names>I. D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gorbenko</surname>
            <given-names>Y. I.</given-names>
          </string-name>
          : Applied Cryptology: Theory. Practice.
          <string-name>
            <surname>Application</surname>
          </string-name>
          (
          <year>2012</year>
          )
          <string-name>
            <surname>Buriachok</surname>
            ,
            <given-names>V. L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gulak</surname>
            ,
            <given-names>G. M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Horoshko</surname>
            ,
            <given-names>V. O.</given-names>
          </string-name>
          :
          <article-title>On the issue of organizing and conducting intelligence in cyberspace</article-title>
          .
          <source>Sci. Def</source>
          .
          <volume>2</volume>
          :
          <fpage>19</fpage>
          -
          <lpage>23</lpage>
          (
          <year>2011</year>
          ). [Publication in Ukrainian]
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Olіynikov</surname>
          </string-name>
          , Р., et al.:
          <article-title>Principles of construction and main properties of the new national standard of block encryption of Ukraine</article-title>
          . Inf. Prot.
          <volume>17</volume>
          (
          <issue>2</issue>
          ):
          <fpage>142</fpage>
          -
          <lpage>157</lpage>
          (
          <year>2015</year>
          ). [Publication in Ukrainian]
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Kelsey</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Schneier</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Wagner</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hall</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          :
          <article-title>Side channel cryptanalysis of product ciphers</article-title>
          .
          <source>5th European Symposium on Research in Computer Security</source>
          :
          <fpage>97</fpage>
          -
          <lpage>111</lpage>
          (
          <year>1998</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8. Safety and Security Consultancy, https://www.ippf.org/about-us/jobs-and
          <article-title>-opportunities/ consultancy/safety-and-security-consultancy</article-title>
          ,
          <source>last accessed</source>
          <year>2020</year>
          /05/11.
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Lloyd</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          , et al.:
          <source>Reference Book on Applied Statistics</source>
          . Vol.
          <volume>2</volume>
          (
          <year>1990</year>
          ). [Publication in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Brodsky</surname>
            ,
            <given-names>B. E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Darkhovsky</surname>
            ,
            <given-names>B. S.</given-names>
          </string-name>
          :
          <article-title>Non-Parametric Statistical Diagnosis: Problems and Methods (</article-title>
          <year>2000</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Shiryaev</surname>
            ,
            <given-names>A. N.</given-names>
          </string-name>
          :
          <source>Probability</source>
          . Vol.
          <volume>1</volume>
          (
          <year>2007</year>
          ). [Publication in Russian]
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Vysochansky</surname>
            ,
            <given-names>D. F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Petunin</surname>
            ,
            <given-names>Y. I.</given-names>
          </string-name>
          :
          <article-title>Justification of the 3-sigma rule for unimodal distributions</article-title>
          .
          <source>Probab. Theor. Math. Stat</source>
          .
          <volume>21</volume>
          :
          <fpage>23</fpage>
          -
          <lpage>35</lpage>
          (
          <year>1979</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>