<!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>Anomaly detection to predict failures in server systems *</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Tyumen State University</institution>
          ,
          <addr-line>6, Volodarskogo ave., Tyumen, 625003, Russian Federation</addr-line>
        </aff>
      </contrib-group>
      <abstract>
        <p>The article discusses the problem of predicting failures to ensure the uninterrupted state of the server and the application of anomaly detection methodology to solve this problem. Some solutions in this area are briefly analyzed and the advantage of the method based on the Holt-Winters forecasting model is described. A mathematical formulation of the problem of detecting anomalies and a formalized description of the method for solving the research objectives are carried out. Several shortcomings of the standard method for detecting anomalies in the operation of server systems have been identified. Several additions are proposed that allow you to adapt the method for solving work goals and reduce the number of false positives. To improve the forecasting accuracy, point anomalies entering the model input are additionally smoothed by a weighted moving average. To eliminate redundant detections associated with the simultaneous appearance of anomalies, a comprehensive assessment of the server state is introduced. To exclude false alarms associated with noisy data, only those anomalous events are recognized, in which an abnormal state of the server remains for a period of time. Computational experiments were carried out to evaluate the resulting improved method. It is concluded that the proposed additions make it possible to improve the forecasting accuracy of the model and reduce the number of false positives of the method, and the method can be used for early detection of gradual failures in the operation of server systems.</p>
      </abstract>
      <kwd-group>
        <kwd>Anomaly detection</kwd>
        <kwd>Proactive monitoring</kwd>
        <kwd>Failure prediction</kwd>
        <kwd>HoltWinters model</kwd>
        <kwd>Brutlag method</kwd>
        <kwd>Server health</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        Server systems are often one of the central elements in the corporate information system
of enterprises. Server failures can lead to the loss of valuable information and
significant financial costs, which necessitates ensuring the smooth operation of such
systems. One of the relevant approaches for these purposes is the introduction of
proactive monitoring systems [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], which are aimed at predicting failures in the
operation of server systems, which allows you to eliminate possible problems at the
stage of their inception.
      </p>
      <p>
        One of the tools for implementing proactive monitoring systems are anomaly
detection methods. An anomaly is understood as a piece of data in which the behavior
of an object significantly differs from the standard or expected behavior. In this context,
the task is reduced to the detection of contextual and collective new anomalies in real
time [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. This article discusses only hardware and software anomalies that are not
caused by deliberate actions of intruders. Prediction of gradual failures in this case is
based on the assumption that if an anomaly occurs on the server, after a while, its
functioning may be disrupted.
      </p>
      <p>
        Analysis of the works shows that there is no universal method for detecting
anomalies for solving any problem [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. To improve the accuracy of anomaly
recognition, it is necessary to take into account the peculiarities of the problem and the
data under investigation. Anomaly detection methods based on tagged data [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] are less
suitable in this area, since the change in the functional load on the server necessitates
constant retraining of the classifier. One Class SVM is used to detect anomalies [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ], but
for its application it is necessary to have an idea of the percentage of anomalies in
datasets, which is not always possible. Choosing this option unjustifiably can lead to
gaps or false detections. Methods based on artificial neural networks are quite
widespread, from simple perceptrons to complex recurrent networks, for example,
LSTM [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. The disadvantages of ANNs in this area are high computational complexity
and black box operation. Because of the latter, it becomes difficult to interpret the
results, which is why the localization of the problem in the server operation requires a
separate analysis. A number of methods based on forecasting models are also
widespread [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. The disadvantage of models from the ARIMA family [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] is that the
models work with stationary data, and converting the server parameter values to a
stationary form according to the Box-Jenkins methodology can be quite laborious.
      </p>
      <p>
        The paper considers an anomaly detection method based on the Holt-Winters
forecasting model [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ]. To use the method, no tagged data is required, which expands
the possible options for its use. The method does not require preliminary processing of
the input data, which is an advantage for analyzing many server parameters. The
method has a relatively low computational complexity, which makes it possible to
process information about a large number of server parameters in short periods and
thereby obtain up-to-date information about the state of the object. The results of the
method are clear and amenable to interpretation. More details on the application of the
Holt-Winters model for predicting the state of servers are described in another article
by the authors [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ], this work is based on the results of the previous study.
      </p>
      <p>The purpose of this work is to study and supplement the anomaly detection method
based on the Holt-Winters model for predicting failures in server systems. Within the
framework of this goal, the article first provides a mathematical formulation of the
problem of detecting anomalies and a formalized description of the application of the
problem to the server's operation. Then, a number of additions to the standard anomaly
detection method are proposed, which allow it to be adapted to achieve the set goal and
to reduce the number of redundant detections. Taking into account the additions, a
modified method for detecting anomalies is described. After that, computational
experiments are carried out to evaluate the proposed additions. At the last stage, the
results were discussed and tasks for further research were proposed.</p>
    </sec>
    <sec id="sec-2">
      <title>Materials and methods</title>
      <p>Each corporate server S can be described with a list of parameters P characterizing
its state, P = {x1, x2 ,.., xn} . Then the state of the server Zt at a moment in time t is
described by a vector of values of its parameters.</p>
      <p>= (⋮</p>
      <p>)
 1,
  ,</p>
      <p>The standard formulation of the problem of detecting anomalies consists in
constructing a functional dependence f : xi,t →{−1,1} such that:
x
i,t = 
normal, f (xi,t ) = 1
anomaly, f (xi,t ) = −1
,</p>
      <sec id="sec-2-1">
        <title>Where xi,t is the value of the i parameter at the moment t.</title>
        <p>To analyze the state of the server, it is necessary to collect statistical information
about its work. The values of each parameter xi  P known at discrete points in time
t = 1, T represent a time series Xi = {xi,1, xi,2 ,.., xi,T} .</p>
        <p>To analyze the state of the server, it is necessary to collect statistical information
about its work. The values of each parameter known at discrete points in time represent
a time series.</p>
        <p>
          The basic technique for detecting anomalies in time series, based on the forecasting
model, consists in generating a forecast for each parameter value and calculating the
forecast error. A significant deviation of the forecast from the actual value indicates the
presence of an anomaly [
          <xref ref-type="bibr" rid="ref10">10</xref>
          ].
        </p>
        <p>
          To simplify further notation, information about an arbitrary parameter x О P is
described without specifying an index. The Holt-Winters multiplicative model was
used as a forecasting model. The model is described by the following system of
equations [
          <xref ref-type="bibr" rid="ref8">8</xref>
          ]:
 
  −


  =  ∗
        </p>
        <p>+ (1 −  )∗ (  −1 +   −1)


 =  ∗   + (1 −  )∗   −
 =  ∗ (  −   −1)+ (1 −  )∗   −1 ,
{  + = (  +  ∗   )∗   − +</p>
        <p>
          Where Rt is an exponentially smoothed series; Tt - trend value; St - the seasonal
component of the series; xt+ m - forecast for m steps; a,b , g - coefficients of
smoothing of the series, trend and seasonality; L - seasonality period. Training of the
model (3) consists in the selection of the optimal values a,b , g . The selection is carried
(1)
(2)
(3)
out using an enumeration of possible combinations of parameters in order to minimize
the loss function for cross-validation [
          <xref ref-type="bibr" rid="ref11">11</xref>
          ].
        </p>
        <p>
          To detect anomalous values, the Brutlag method is used [
          <xref ref-type="bibr" rid="ref12">12</xref>
          ]. At the first stage, the
measure of deviation is determined dt : dt =  * | xt − xt | +(1 −  ) * dt −L . Values xt , xt
and parameters g, L are contained in model (3). Next, the upper and lower boundaries
of the predicted value are calculated:
 xmax t = xt + m * dt −L ,
 xmin t = xt − m * dt −L
        </p>
        <p>1, xt [ xmin t , xmax t ]
f ( xt ) = 
−1, xt [ xmin t , xmax t ]
.</p>
        <p>Where m is the coefficient of the interval width.</p>
        <p>As a result, the value xt is considered abnormal if it is outside the boundaries of the
predicted value (4), and normal in the opposite case:</p>
        <p>Using the Brutlag method to detect anomalies has a number of advantages: the
administrator does not need to manually create thresholds, the threshold values are not
static and are adjusted to the server state, and this method allows you to detect many
anomalies. But when applying the method to solve research problems, a number of
shortcomings were identified. Below is a number of author's additions that allow us to
eliminate the disadvantages indicated below and to adapt the use of the Holt-Winters
model and the Brutlag method for detecting anomalies in the operation of server
systems and some other types of complex information systems.</p>
        <p>
          Appendix 1. It has been investigated that the Holt-Winters model (3) shows less
accurate forecast results on the series with an increased noise level [
          <xref ref-type="bibr" rid="ref9">9</xref>
          ]. Reducing the
accuracy of predictions can lead to gaps or false positives in the anomaly detection
method. To eliminate this disadvantage, the following solution is proposed.
        </p>
        <p>
          Parameter values that are recognized by the Brutlag method as anomalous, and do
not retain their anomalous value in the future, are considered single random outliers in
the data. When forming the next forecast, such emissions are additionally smoothed by
replacing the actual value with a weighted moving average [
          <xref ref-type="bibr" rid="ref13">13</xref>
          ] based on K the
previous normal values:
        </p>
        <p>K −1
 (K − i) * xt −i
xt* = i=0</p>
        <p>,
K
 i
i=1</p>
        <p>Where xt* is the value of the point outlier supplied to the input of the Holt-Winters
model. This addition presumably leads to a decrease in the influence of noise on the
formation of a forecast and an increase in its accuracy.
(4)
(5)
Appendix 2. The server is a multi-parameter system; the occurrence of anomalies in its
operation can be caused by the exit of several parameters at once outside the permissible
intervals. The standard approach identifies each deviation as a new anomaly.</p>
        <p>To exclude redundant detections, not the value of an individual parameter is
analyzed, but the state of the object as a whole. The state of the server Zt is considered
normal Z N if all values of the server parameters are in the normal range, and abnormal</p>
      </sec>
      <sec id="sec-2-2">
        <title>Z A if at least one parameter is in the abnormal range.</title>
        <p>(6)
(7)
Zt = ZN , x  P, xt [ xmin t , xmax t ] .</p>
        <p>
          ZA , x  P, xt [ xmin t , xmax t ]
Appendix 3. Time series describing server states are not stationary and can be quite
noisy [
          <xref ref-type="bibr" rid="ref9">9</xref>
          ], which leads to false positives of the anomaly detection method. Such
redundant discoveries are not of interest for analyzing server health.
        </p>
        <p>To exclude such alarms, only those anomalous events A are recognized, in which the
abnormal state of the server remains at the time interval p:</p>
        <p>A = Z1,.., Z p  , Zt , Zt = Z A , t = 1, p .</p>
        <p>Taking into account the additions, the modified anomaly detection method consists
of the following steps:
─ Collect statistical data on the values of parameters that characterize the state of the
server;
─ Determine the value of seasonality L and smoothing parameters a, b, g for model
(3), the coefficient of the width of the interval for detecting anomalies;
─ Collect new actual values of the server parameters that require analysis;
─ For each parameter, generate a forecast based on the history of previous values,
taking into account the condition: if there are point anomalies in the history, replace
each outlier value with a weighted moving average based on K of previous normal
values;
─ Determine the measure of deviation of the predicted values from the actual ones;
─ For each forecast, determine the boundaries of the area of normal values (4);
─ For each parameter value, determine the belonging to the area of normal or abnormal
data (5);
─ Record the anomalous state of the server if there is at least one parameter that is in
the area of anomalous values (6);
─ Record the occurrence of an anomalous event, while maintaining an anomalous state
of the server during p time steps (7);
─ Repeat from step 3.</p>
        <p>To assess the correctness of the proposed improved method for detecting anomalies,
computational experiments were carried out, the results of which are presented in the
next section of the article.</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>Results</title>
      <p>To prepare an experimental study, the values of 4 parameters were recorded from 2
computer servers within 3 weeks with a discreteness of 1 minute. The sufficiency of
these parameters for assessing the state of servers is not considered in the article. The
statistics for the first two weeks served as a training sample for further experiments.
The selection of the optimal smoothing coefficients a, b, c for the Holt-Winters model
was made on the training sample. Seasonality parameter L corresponds to a weekly
period and is equal to 1440 measurements. The parameter of the interval width m for
the Brutlag method is empirically selected as 6 units. A smaller value leads to a larger
number of detections that slightly deviate from the normal range.</p>
      <p>
        On the control sample, for each parameter value, a short-term forecast one step ahead
was formed based on the previous history. The accuracy of the obtained model results
was assessed using the MAPE metric [
        <xref ref-type="bibr" rid="ref14">14</xref>
        ]:

=

1

∑
 =1
|  −   |
 
∗ 100,
      </p>
      <p>Where xt is the actual value of the series, xt is the predicted value of the series, N
is the number of measurements.</p>
      <p>Two options for the formation of forecasts are considered. In the first version, the
average accuracy of which is indicated in the table as MAPE1 , the classical application
of the Holt-Winters model was carried out. The implementation of the second version
of the forecasts ( MAPE2 ) was carried out taking into account the first addition
described in the previous section. The previously recognized point anomalous values,
when fed into the model (3), are replaced by a weighted moving average of the three
previous normal values. Information about server parameters, smoothing coefficients
and experiment results is presented in Table 1.</p>
      <p>According to the results of the experiment, it can be seen that, in general, the
HoltWinters model shows a fairly high accuracy in predicting the values of server
parameters. For parameters with a high random component, the forecast accuracy is
slightly lower. For such parameters, an increase in accuracy is observed when applying
the principle described in the first supplement of the previous section.</p>
      <p>For the first server:
For the second server:</p>
      <p>2.65  MAPE1  9.5, MAPE1 = 4.88.
2.65  MAPE2  9.08, MAPE2 = 4.77.</p>
      <p>2.4  MAPE1  12.7, MAPE1 = 5.5.
2.4  MAPE2  11.29, MAPE</p>
      <p>(8)
(9)
(10)</p>
      <p>Using the described additions, the method recognized only those anomalous values
that were simulated by software. The Brutlag method additionally recognized point
random outliers that are caused by the partially stochastic behavior of the server and
are not of interest for diagnosing the state of the object.</p>
      <p>
        The occurrence of an anomalous event on the Memory Usage graph of the first server
is shown in Figure 1. The blue line on the graph denotes actual values, orange
predicted values. The green and brown lines indicate the upper and lower boundaries
of the normal range. The occurrence of an anomalous event is marked in the figure with
red dots.
The development of anomaly detection methods based on the Holt-Winters forecasting
model is of practical and research interest, which is confirmed by research in this
direction [
        <xref ref-type="bibr" rid="ref16 ref17">16-17</xref>
        ]. The authors of the work made the following contribution in this area.
      </p>
      <p>
        Based on the results of the experiment, it can be argued that the Holt-Winters model
shows good results when predicting the values of server parameters. A study on the
analysis of the SARIMA model for predicting server load demonstrates less accurate
forecast results [
        <xref ref-type="bibr" rid="ref18">18</xref>
        ]. But in the presence of noise and emissions in the structure of time
series, the accuracy of the forecasts formed by the model (3) decreases [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]. To reduce
the impact of the random component, the following was undertaken. Recognition of
point anomalous values in datasets is performed. At the next generation of the forecast,
the detected outliers are additionally smoothed by replacing the actual value with a
weighted moving average. The results of the experiment demonstrated that such an
addition allows one to obtain more accurate forecasts for time series with an increased
level of noise. The improvement in the accuracy of the model for a given loss function
for individual parameters is up to 12.5%, the average forecast accuracy increases to
7.2%.
      </p>
      <p>
        When using the Brutlag method, it is recommended [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ] to choose the coefficient of
the width of the interval m (4) in the range from 2 to 3 units. Application of the method
in practice has shown that setting the coefficient in such a range leads to a large number
of method triggers, slightly deviating from the range of normal values. In order to
reduce the excess operation during the experiments, the coefficient was set equal to 6
units.
      </p>
      <p>It was revealed that the Brutlag method generates redundant triggers in the case when
several server parameters are simultaneously abnormal. To exclude redundant
detections, not the value of an individual parameter is analyzed, but the state of the
object as a whole. In addition, the method detects point random outliers that are caused
by the partially stochastic behavior of the server and are not of interest for diagnosing
its state. To exclude such alarms, only those anomalous events are recognized, in which
the anomalous state of the server persists for a time interval. The proposed additions to
the method make it possible to reduce the number of false positives, which is
demonstrated using a computational experiment.</p>
      <p>The application of the modified method makes it possible to recognize contextual
collective anomalies, which may indicate the occurrence of gradual failures in the
functioning of the server, which makes it possible to eliminate possible problems even
at the stage of their inception.
5</p>
    </sec>
    <sec id="sec-4">
      <title>Conclusion</title>
      <p>
        The work is based on the results of a previous study by the authors [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ], which
considered the application of the Holt-Winters model to predict the state of server
systems. The aim of the current work is to investigate and supplement an anomaly
detection method based on the Holt-Winters model for predicting failures in server
systems. At the first stage of the work, a mathematical formulation of the problem of
detecting anomalies and a formalized description of the application of the problem to
the work of the server are carried out. The existing method of detecting anomalies has
a number of shortcomings, on the basis of which a number of additions are proposed
that allow it to be adapted to achieve the set goal and reduce the number of redundant
detections. To improve forecasting accuracy, single outliers entering the model input
are smoothed with a weighted moving average. To eliminate redundant detections
associated with the simultaneous appearance of anomalies, a comprehensive
assessment of the server state is introduced. To exclude false alarms associated with
noisy data, only those anomalous events are recognized, in which an abnormal state of
the server remains for a period of time. Further, computational experiments are carried
out to evaluate the resulting improved method. Based on the results of the experiments,
it was concluded that the proposed additions make it possible to improve the forecasting
accuracy of the model and reduce the number of false positives of the method, and the
method can be used for early detection of gradual failures in the operation of server
systems. Further prospects for the development of the research area include the
classification of recognized anomalies and support for the administrator's
decisionmaking when such events are detected.
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Kothamasu</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Huang</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          , VerDuin W.:
          <article-title>System health monitoring and prognostics - a review of current paradigms and practices</article-title>
          .
          <source>The International Journal of Advanced Manufacturing Technology</source>
          ,
          <volume>28</volume>
          (
          <fpage>9</fpage>
          -
          <lpage>10</lpage>
          ),
          <fpage>1012</fpage>
          -
          <lpage>1024</lpage>
          (
          <year>2006</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Chandola</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Banerjee</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kumar</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          :
          <article-title>Anomaly detection: A survey</article-title>
          .
          <source>ACM computing surveys (CSUR)</source>
          ,
          <volume>41</volume>
          (
          <issue>3</issue>
          ),
          <fpage>1</fpage>
          -
          <lpage>58</lpage>
          (
          <year>2009</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Omar</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ngadi</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jebur</surname>
          </string-name>
          , H.:
          <article-title>Machine learning techniques for anomaly detection: an overview</article-title>
          .
          <source>International Journal of Computer Applications</source>
          ,
          <volume>79</volume>
          (
          <issue>2</issue>
          ) (
          <year>2013</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Li</surname>
            ,
            <given-names>K.L.</given-names>
          </string-name>
          et al.:
          <article-title>Improving one-class SVM for anomaly detection</article-title>
          .
          <source>In: Proceedings of the 2003 International Conference on Machine Learning and Cybernetics</source>
          ,
          <volume>5</volume>
          ,
          <fpage>3077</fpage>
          -
          <lpage>3081</lpage>
          , IEEE,
          <source>Xi'an (</source>
          <year>2003</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Ergen</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kozat</surname>
            <given-names>S.</given-names>
          </string-name>
          :
          <article-title>Unsupervised Anomaly Detection With LSTM Neural Networks</article-title>
          .
          <source>In: IEEE Transactions on Neural Networks and Learning Systems</source>
          ,
          <volume>31</volume>
          ,
          <issue>8</issue>
          ,
          <fpage>3127</fpage>
          -
          <lpage>3141</lpage>
          , IEEE (
          <year>2019</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Gupta</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          et al.:
          <article-title>Outlier detection for temporal data: A survey</article-title>
          .
          <source>In: IEEE Transactions on Knowledge and data Engineering</source>
          ,
          <volume>26</volume>
          ,
          <issue>9</issue>
          ,
          <fpage>2250</fpage>
          -
          <lpage>2267</lpage>
          , IEEE (
          <year>2013</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Moayedi</surname>
            ,
            <given-names>H.Z.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Masnadi-Shirazi</surname>
            ,
            <given-names>M.A.</given-names>
          </string-name>
          :
          <article-title>Arima model for network traffic prediction and anomaly detection</article-title>
          .
          <source>International Symposium on Information Technology</source>
          ,
          <volume>4</volume>
          ,
          <fpage>1</fpage>
          -
          <lpage>6</lpage>
          , IEEE, Kuala
          <string-name>
            <surname>Lumpur</surname>
          </string-name>
          (
          <year>2008</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Kalekar</surname>
            ,
            <given-names>P.S.:</given-names>
          </string-name>
          <article-title>Time series forecasting using holt-winters exponential smoothing</article-title>
          .
          <source>Kanwal Rekhi School of Information Technology</source>
          ,
          <volume>4329008</volume>
          (
          <issue>13</issue>
          ) (
          <year>2004</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Dubrovin</surname>
            ,
            <given-names>M.G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gluhih</surname>
            ,
            <given-names>I.N.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Karyakin</surname>
          </string-name>
          , I.Y.:
          <article-title>Forecasting the server status using the triple exponential smoothing model</article-title>
          .
          <source>J. Phys.: Conf. Ser.</source>
          ,
          <volume>1661</volume>
          ,
          <issue>012031</issue>
          (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Jinka</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          :
          <article-title>Anomaly detection for monitoring: A statistical approach to time series anomaly detection</article-title>
          .
          <source>O'Reilly Media, USA</source>
          ,
          <fpage>15</fpage>
          -
          <lpage>31</lpage>
          (
          <year>2015</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Hyndman</surname>
          </string-name>
          , R. J.:
          <article-title>Measuring forecast accuracy. Business forecasting: Practical problems</article-title>
          and solutions,
          <volume>177</volume>
          -
          <fpage>183</fpage>
          (
          <year>2014</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Brutlag</surname>
            ,
            <given-names>J.D.</given-names>
          </string-name>
          <article-title>Aberrant Behavior Detection in Time Series for Network Monitoring</article-title>
          . LISA,
          <volume>14</volume>
          (
          <year>2000</year>
          ),
          <fpage>139</fpage>
          -
          <lpage>146</lpage>
          (
          <year>2000</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          13.
          <string-name>
            <surname>Zhuang</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          et al.:
          <article-title>A weighted moving average-based approach for cleaning sensor data</article-title>
          .
          <source>In: 27th International Conference on Distributed Computing Systems (ICDCS'07)</source>
          ,
          <fpage>38</fpage>
          -
          <lpage>38</lpage>
          , IEEE, Toronto (
          <year>2007</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          14.
          <string-name>
            <surname>Shcherbakov</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          et al:
          <article-title>A survey of forecast error measures</article-title>
          .
          <source>World Applied Sciences Journal</source>
          ,
          <volume>24</volume>
          ,
          <fpage>171</fpage>
          -
          <lpage>176</lpage>
          (
          <year>2013</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          15.
          <string-name>
            <surname>Elmrabit</surname>
          </string-name>
          , N. et al.:
          <article-title>Evaluation of machine learning algorithms for anomaly detection</article-title>
          .
          <source>In: 2020 International Conference on Cyber Security and Protection of Digital Services (Cyber Security)</source>
          ,
          <fpage>1</fpage>
          -
          <lpage>8</lpage>
          , IEEE, Dublin, Ireland (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          16.
          <string-name>
            <surname>Szmit</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          et al.:
          <article-title>Implementation of Brutlag's algorithm in Anomaly Detection 3.0</article-title>
          .
          <source>In: 2012 Federated Conference on Computer Science and Information Systems (FedCSIS)</source>
          ,
          <fpage>685</fpage>
          -
          <lpage>691</lpage>
          , IEEE,
          <string-name>
            <surname>Wroclaw</surname>
          </string-name>
          (
          <year>2012</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          17.
          <string-name>
            <surname>Ekberg</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ylinen</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Loula</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          :
          <article-title>Network behaviour anomaly detection using Holt-Winters algorithm</article-title>
          .
          <source>International Conference for Internet Technology and Secured Transactions</source>
          ,
          <fpage>627</fpage>
          -
          <lpage>631</lpage>
          , IEEE, Abu
          <string-name>
            <surname>Dhabi</surname>
          </string-name>
          (
          <year>2011</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          18.
          <string-name>
            <surname>Debusschere</surname>
          </string-name>
          , V.:
          <article-title>Hourly server workload forecasting up to 168 hours ahead using seasonal ARIMA model</article-title>
          .
          <source>In: 2012 IEEE international conference on industrial technology</source>
          ,
          <fpage>1127</fpage>
          -
          <lpage>1131</lpage>
          , IEEE, Athens (
          <year>2012</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>