Passive concept drift handling via variations of learning vector quantization

Concept drift is a change of the underlying data distribution which occurs especially with streaming data. Besides other challenges in the field of streaming data classification, concept drift has to be addressed to obtain reliable predictions. Robust Soft Learning Vector Quantization as well as Generalized Learning Vector Quantization has already shown good performance in traditional settings and is modified in this work to handle streaming data. Further, momentum-based stochastic gradient descent techniques are applied to tackle concept drift passively due to increased learning capabilities. The proposed work is tested against common benchmark algorithms and streaming data in the field and achieved promising results.


Introduction
A key concept in machine learning is to separate the training step of the model and the evaluation phase. However, this is not applicable in the domain of stream classification. In this field, it is assumed that data arrive continuously, making the storage of data in memory unfeasible. Further, not all training data are available at training time, raising the need for constantly updating the model, e.g., online or incremental. Stream classification algorithms [1] tackle these requirements.
However, stream classifiers are prone to concept drift of streaming data, which is a change of underlying distribution and could lead to a collapse in prediction performance. There are various types of concept drift, i.e., incremental, abrupt, gradual and reoccurring, and as a consequence, a variety of approaches addressing these issues have been proposed [13]. In general, these strategies are separated into active and passive [22]. Active adaptation changes a model noticeable. The passive ones use no explicit detection strategy, but are continually updating the model, without awareness of concept drift.
The family of prototype-based classification algorithms, the Learning Vector Quantization (LVQ) [19], receives much attention as a potential stream classification algorithm due to its online learning capabilities [30]. In [30], it was also shown that the usage of weight decay has the potential to improve the generalization error of the LVQ in non-stationary environments. Other LVQ variants have not been considered as stream classifiers yet, except from our adaptive Robust Soft Learning Vector Quantization (RSLVQ) in [16]. Therefore, we propose adaptive versions of RSLVQ and Generalized Learning Vector Quantization (GLVQ) which maximize the objective function with momentum-based gradient descent/ascent in Sect. 4.1. This is applied to increase learning speed and rapid adaptation to occurring concept drifts. As we show in the experiments, the LVQ variants with constant learning rate perform considerably worse than with weight decay or momentumbased gradient descent.
In summary, we provide the following contributions in this article: 1. we apply momentum-based gradient techniques known from deep learning to prototype-based stream learning 2. we provide two adaptive RSLVQ versions using Adadelta and Adamax 3. we provide two adaptive GLVQ versions using Adadelta and Adamax 4. we provide experiments comparing our proposed methods against state-of-the-art classifiers without hyperparameter optimization to get a realistic comparison in non-stationary environments The paper is structured as follows: Sect. 2 presents related work on stream classification and adaptive gradient optimization. In Sect. 3, the fundamentals and particularities of streaming data and concept drift are shown. In Sect. 4, we modify RSLVQ and GLVQ by replacing the gradient optimization technique by adaptive optimization algorithms. The comparison of our proposed methods to baseline RSLVQ/GLVQ and standard algorithms on common data streams in the field is presented in Sect. 5. Section 6 summarizes up the results of this work.

Related work
In the field of streaming data, different kinds of algorithms successfully apply passive drift handling to evolving data streams, which is often done by using a fixed size window of recent samples. In [22], a modern approach of a self-adjusting memory version of the K-Nearest Neighbor (KNN) which is called SAM-KNN is proposed. SAM-KNN is developed to handle various types of concept drift, using biologically inspired memory models and their coordination. The basic idea is to store dedicated models for current and former concepts, used according to the demands of the given situation [22].
The Adaptive Windowing algorithm (ADWIN) [2] is a drift detector and works by keeping updated statistics of a variable sized window, such that changes can be detected. It performs cuts in its window to better adapt to the learning algorithms. Kolmogorv-Smirnov Windowing (KSWIN) [25] follows a similar concept, but uses a more sensitive statistical test.
For evolving data stream, classification tree-based algorithms, like the famous Hoeffding Tree (HT) [9], are common [6]. To address the problem of concept drift, an adaptive HT with ADWIN as drift detector was published in [3] and showed better prediction performance on evolving data streams as the classical HT.
Also, ensemble models are used to combine multiple classifiers [14] in the streaming domain. The OzaBagging (OB) [23] algorithm is an online ensemble model that uses different base classifiers. For the use of concept drift detection, this ensemble model can be again combined with the ADWIN algorithm to the OzaBagging ADWIN (OBA) algorithm. For a comprehensive data stream description, see [21].
Also, LVQ algorithms first introduced by [19] have received attention as potential stream classification algorithms [30]. The Robust Soft Learning Vector Quantization (RSLVQ) [29] is a promising probabilistic classifier, which assumes class distributions as Gaussian mixture models learned via Stochastic Gradient Ascent (SGA) and has only been evaluated as a stream classifier in a previous version of this article [16] so far. Also, Generalized Learning Vector Quantization has not been considered as a stream classifier yet.
In the advent of deep learning, variations in SGD algorithms receive more and more attention. A comparison is given in [27]. It has been shown that momentum-based gradient descent algorithms like Adadelta [34] and RMSprop [33] converge faster to better optimums as traditional approaches. Both of these algorithms are extensions of the Adagrad [10] algorithm. Note that momentumbased gradient descents have been applied to a prototypebased learner for stationary environments in [20], showing Adam gradient update technique as well working algorithm. Furthermore, there is a proposed extension of Adam which is called AdaMax [18], which has not been considered as update technique of a LVQ algorithm yet.
In this work, we reformulate Adadelta and AdaMax as SGA optimizer and apply it to RSLVQ and to GLVQ as a gradient descent optimizer.
3 Streaming data and concept drift

Streaming data
In the context of supervised learning, a data stream is given as a sequence S ¼ fs 1 ; . . .; s t ; . . .g of tuples s i ¼ fx i ; y i g, with potentially infinite length. A tuple s i ¼ fx i ; y i g contains the data point x i 2 R d and the respective label y i ¼ f1; . . .; Cg whereby s t arrives at time t. A classifier predicts labelsŷ t of unseen data x t 2 R employing a prior model h tÀ1 , i.e.,ŷ t ¼ h tÀ1 ðx t Þ. The prior model and s t are subsequently included into the new model h t ¼ learnðh tÀ1 ; s t Þ. This behavior is called immediate or test-then-train evaluation. It also means that the learning algorithm trains for an infinite time and that at time t the classifier is (re-)trained on the just arrived tuple.

Concept drift
Concept drift is the change of joint distributions of a set of samples X and corresponding labels y between two points in time: 9X : pðX; yÞ t 6 ¼ pðX; yÞ tÀ1 ð1Þ The term virtual drift refers to a change in distribution pðXÞ for two points in time, without affecting pðyjXÞ. Note that we can rewrite Eq.(1) to 9X : pðXÞ t pðyjXÞ t 6 ¼ pðXÞ tÀ1 pðyjXÞ tÀ1 ; ð2Þ A virtual concept drift may appear in conjunction with a real concept drift. Figure 1 shows the four common drift types. For a comprehensive study of these concept drift types, see [13]. The stability-plasticity dilemma [13] defines the trade-off between incorporating new knowledge into models (plasticity) and preserves prior knowledge (stability). This prevents stable performance over time because on the edge of a drift, significant efforts go into learning and testing against new distributions.

Robust Soft Learning Vector Quantization
The Robust Soft Learning Vector Quantization (RSLVQ) [29] is a probabilistic prototype-based classification algorithm. Given a labeled dataset X ¼ fðx i ; y i Þ 2 R d Â f1; . . .; Cgg n i¼1 , with data points x i , labels y i , C as the number of classes, d the number of features and n as the number of samples, the RSLVQ algorithm trains a   [13] prototype model such that the error on the classification task is minimized. The RSLVQ model consists of a set of m prototypes H ¼ fðh j ; y j Þ 2 R d Â f1; . . .; Cgg m j¼1 . Each prototype represents a multi-variate Gaussian model, i.e., N ðh j ; rÞ, approximating an assumed class dependent Gaussian mixture of X. The goal of RSLVQ algorithms is to learn prototypes, representing the class dependent distribution, i.e., corresponding class samples x i should be mapped to the correct class or Gaussian mixture based on the highest probability. The RSLVQ algorithms maximize the maximum likelihood ratio as objective function: where pðx i ; y k jHÞ is the probability density function that x is generated by the mixture model of the same class and pðx i jHÞ is the overall probability density function of x given H. Equation (3) will be optimized with SGA, i.e., h l ðt þ 1Þ ¼ h l ðtÞ À a Ã ðP y ðljxÞ À PðljxÞÞðx À h l Þ; c l ¼ y; ÀPðljxÞðx À h l Þ; c l 6 ¼ y: where P y ðljxÞ is the assignment probability that x is assigned to the component l of the mixture model of same class prototype h y and PðljxÞ is the probability that x is assigned to the component l of the mixture using all classes. The parameter a Ã ¼ a r 2 is the learning rate. For a more comprehensive derivation, see [29].

Generalized Learning Vector Quantization
GLVQ is a prototype-based classifier and has many things in common with RSLVQ. However, in GLVQ prototypes are not representing a Gaussian model and prototypes are updated using the winner-takes-all rule. Assume h þ represents the nearest prototype of the same class as x and h À represents the nearest prototype belonging to a different class than x. Consider the relative distance difference: where d þ is the distance between h þ and x and d À is the distance between h À and x. lðxÞ is in ½À1; þ1, and if l is negative, x is classified correct, else it is classified incorrect. To reduce error rates, lðxÞ should be minimized for all input vectors. Hence, the goal is to minimize a cost function S: where n is the number of input vectors for training and f ðlÞ is a monotonically increasing function [28]. To minimize the cost function, the prototypes h þ and h À are updated by gradient descent, using a learning rate a: Using squared euclidean distance d l ¼ kx À h l k 2 , we can obtain the following learning rule [28]: of ol is a kind of update weight depending on x. To decrease the error rate, it is effective to update the prototypes mainly by the input vectors near the class boundaries, so that decision boundaries are shifted toward Bayes limits. Hence, f ðlÞ should be a nonlinear monotonically increasing function and it is considered that the classification ability depends on the definition of f ðlÞ. We choose of ol ¼ f ðl; tÞð1 À f ðl; tÞÞ, where t is the time step and f ðl; tÞ is a sigmoid function with 1 1þe Àlt as proposed in [28].

Momentum-based optimization
One of the most common algorithms to optimize error functions is the gradient descent/ascent algorithm, and in particular the stochastic formulations SGD and SGA [27].
In the field of Deep Learning, the classic SGA has been further modified to reduce common problems, like sensitivity to steep imbalanced valleys, i.e., areas where the cost surface curves are much more steep in one dimension than in another [32], which are common around local optima. Momentum [24] is a method that helps accelerate SGA in the relevant direction and dampens oscillations. It does this by adding a fraction c of the update vector v of the previous time step to the current update vector: The momentum term c is usually set to 0.9 or a similar value [18,27,33] but can be seen as a hyperparameter which can be optimized, e.g., via grid search. A common range for optimizing c is [0.9, 0.99] [34], while lower decay rates like c ¼ 0:5 are appropriate when having larger gradients, e.g., at the first iterations of optimization [33]. Effective momentum-based techniques are Adagrad, Adadelta, Adam, AdaMax, and RMSprop [27]. While Adagrad was the first publication of these three algorithms, RMSprop and Adadelta are further developments of Adagrad, which both try to reduce its aggressive, monotonically decreasing learning rate [34]. These momentumbased algorithms diverge to a local optima much faster and sometimes reach better optima than SGA [27].
Due to the fact that momentum-based algorithms make larger steps per iteration, it should adapt faster to new concepts. Thus, we implemented this idea into the RSLVQ and GLVQ and replaced the gradient update rule by Adadelta and AdaMax.

Adadelta
Instead of accumulating all past squared gradients, like in the Adagrads momentum approach, Adadelta restricts the window of accumulated past gradients to some fixed size w.
While it could inefficiently store all w squared gradients, the sum of gradients is recursively defined as a decaying average of all past squared gradients instead. The running average E½g 2 t at time step t of the squared gradient g 2 then depends (as a fraction c similarly to the Momentum term) only on the previous average and the current gradient: The decay rate c should be set to a value of around 0.9 [10].
Hence, in [10] another exponentially decaying average is introduced, this time not of squared gradients but squared parameter updates: Thus, the update of the root mean squared (RMS) error of parameters is: Because of the reason that RMS½Dh t is not known, it is approximated by the RMS of parameter updates until the previous time step RMS½Dh tÀ1 . Finally, the learning rate g of the previous update rule is replaced with RMS½Dh tÀ1 , and we receive the following equation for updating Adadelta: Due to the fact that the learning rate g has been eliminated, it does not have to be optimized when using Adadelta, which keeps the number of hyperparameters small [27].

Robust Soft Learning Vector Quantization
To exchange the SGA learning of the RSLVQ with the Adedelta algorithm, the learning rule of the RSLVQ is replaced by the update rule provided in Eq. (13). The prior and posterior probabilities are calculated the same way as in the RSLVQ SGA . Based on the posterior and prior, the gradient is calculated over the objective function of the RSLVQ: In the next step, the running average of past squared gradients E½g 2 is updated by Eq. (10). Now the gradient update at time step t can be calculated: Afterward, the prototype update Dh t has to be applied as gradient ascent: Finally, the squared parameter updates E½Dh 2 are stored by Eq. (11).

Generalized Learning Vector Quantization
Exchanging the gradient descent technique of GLVQ is done by replacing the update of Eq. (8) by a gradient calculation: The other steps are the same as for introducing Adadelta into RSLVQ. Only Eq. (14) is replaced by Eq. (17) and the prototype update is a gradient descent instead of ascent:

AdaMax
AdaMax is an extension of Adam algorithm [18]; thus, we will introduce Adam first. Adam stores an exponentially decaying average of past squared gradients v t like Adadelta. Additionally, it keeps an exponentially decaying average of past gradients m t . It also has two separate decay factors: b 1 decays m t and b 2 decays v t . The decaying average of past gradients is computed as follows: The corresponding decaying average of past squared gradients: As m t and v t are initialized as vectors of zeros, Adam is biased to 0, especially during the initial time steps and especially when the decay rates are small (i.e. b 1 and b 2 are close to 1). Thus, the authors of [18] introduce a bias correction for the first and second moment estimates: Finally, these parameters are used to update h: Here g is a learning rate like in vanilla gradient descent.
Recommended default values for Adam, as well as Ada-Max, are b 1 ¼ 0:9, b 2 ¼ 0:999 and ¼ 10 À8 . In Adam, v t scales the gradient inversely proportional to the L 2 norm of the past gradients and current gradients kg t k 2 : This update can be generalized to the L p norm: Norms for large p values become numerically unstable; thus, L 1 and L2 norm are most common in practice. However, L 1 also exhibits stable behavior [27]. Hence, in [18] AdaMax is proposed and shown that v t with L 1 converges to the following more stable value. We denote the infinity constrained v t of Adam as u t : Now Eq. (22) is transformed into: Due to the fact that u t relies on a max operation, it is not suggestible to bias toward zero as in Adam. Hence, a bias correction for u t is not necessary. In this paper, we use a learning rate g of 0.001.

Robust Soft Learning Vector Quantization
To transform the update rule of RSLVQ to AdaMax, we start by calculating g t by Eq. (14). Then, m t is calculated via Eq. (19). We then do the bias correction of m t , following Eq. (21) to obtainm t . Afterward, we obtain u t by using Eq. (25) and calculating the prototype update by an ascending version of Eq. (26): Note that u 0 and m 0 are vectors of zeros.

Generalized Learning Vector Quantization
To change the gradient descent technique of GLVQ to AdaMax, we also start by calculating g t via Eq. (14). Furthermore, we can again calculate m t by Eq. (19) and its bias corrected versionm t via Eq. (21). Afterward, we obtain u t via Eq. (25) and perform the final prototype update by Eq. (26). As we can see, it is the same procedure for exchanging prototypes for both LVQ versions expect for the gradient calculation and final update calculation. Note that we can now address the stability plasticity dilemma [13] of all presented algorithms by our decay factors c; b 1 ; b 2 .

Setup
We compared our classifiers against other state-of-the-art stream classifiers. Evaluation is done using the test-thentrain method as described in Sect. 3. Since accuracy can be misleading on datasets with class imbalances, we also report Kappa statistics j. Kappa is a statistic for imbalanced classes, which compares the classifiers performance with those of a chance classifier. If the classifier is always correct, then j ¼ 1. If its predictions coincide with the correct ones as often as those of a chance classifier, then j ¼ 0 [4].
To test if there are significant differences between the performance of the algorithms, the Friedman [11] test with a 95 % significance level is performed followed by the Bonferroni-Dunn [8] post hoc test.
Nine synthetic and three real data streams are used in the experiments. The synthetic data streams include abrupt, gradual, incremental drifts and one stationary data stream. The real-world data streams have been thoroughly used in the literature [15,22] to evaluate the classification performance of data stream classifiers and exhibit multi-class, temporal dependencies and imbalanced data streams with different drift characteristics. Note that the gradual and abrupt drifts are generated by a concept drift generator which switches the class data generator functions as described below.
LED The LED data set simulates both abrupt and gradual drifts based on the used generator. The generator was first introduced in [7]. This data set yields instances with 24 Boolean features, 17 of which are irrelevant. The remaining 7 features correspond to each segment of a sevensegment LED display. The goal is to predict the digit displayed on the LED display, where each feature has a 10 % chance of being inverted. To simulate drifts in this data set, the relevant features are swapped with irrelevant features. The gradual drift, as well as the abrupt drift, happens at the 250,000th instance of the stream. The first drift swaps three features and is replaced with a drift that swaps seven features. LED g simulates one gradual drift, while LED a simulates three abrupt drifts.
SEA The SEA generator is an implementation of the data stream with abrupt concept drift, first described by Street and Kim in [31]. It produces data streams with three continuous attributes ðf 1 ; f 2 ; f 3 Þ. The range of values that each attribute can assume lies between 0 and 10. Only the first two attributes ðf 1 ; f 2 Þ are relevant, i.e., f 3 does not influence the class value determination. New instances are obtained through randomly setting a point in a two-dimensional space, such that these dimensions correspond to f 1 and f 2 . This two-dimensional space is split into four blocks, each of which corresponds to one of the four different functions. In each block, a point belongs to class 1 if f 1 þ f 2 h and to class 0 otherwise. The threshold h is used to split instances between class 0 and 1, assumes values 8 (block 1), 9 (block 2), 7 (block 3), and 9.5 (block 4). Two important features are the possibility to balance classes, which means the class distribution will tend to a uniform one, and the possibility to add noise, which will, according to some probability, change the chosen label for an instance. In this experiment, the SEA generator is used with 10 % noise in the data stream. SEA g simulates one gradual drift, while SEA a simulates an abrupt drift. Both drifts happen at the 250,000th instance. Figure 2 shows the sliding mean per class of one run of the SEA a generator.
RBF This generator produces data sets utilizing the Radial Basis Function (RBF). This generator creates several centroids, having a random central position and associates them with a standard deviation value, a weight, and a class label. To create new instances, one centroid is selected at random, where centroids with higher weights have more chances to be selected. The new instance input values are set according to a random direction chosen to offset the centroid. The extent of the displacement is randomly drawn from a Gaussian distribution according to the standard deviation associated with the given centroid. Incremental drift is introduced by moving centroids at a continuous rate, effectively causing new instances that ought to belong to one centroid to another with (maybe) a different class. Both RBF m and RBF f were parametrized with 50 centroids, and all of them drift. RBF m simulates a moderate incremental drift (speed of change set to 0.0001) while RBF f simulates a faster incremental drift (speed of change set to 0.001).
Sine The Sine generator [12] creates 4 numerical attributes that vary from 0 to 1, where only 2 of them are relevant to the classification task. A classification function is chosen among four possible ones: The abrupt drift is generated by changing the classification function, thus changing the threshold. In our experiments, we switch the classification function between SINE1 and SINE2.
HYPER The HYPER data set simulates an incremental drift and it was generated based on the hyperplane generator [17]. A hyperplane is a flat, n À 1-dimensional subset of that space that divides it into two disconnected parts. It is possible to change a hyperplane orientation and position by slightly changing its relative size of the weights w i . This generator can be used to simulate time-changing concepts, by varying the values of its weights as the stream progresses [4]. HYPER was parametrized with ten attributes and a magnitude of change of 0.001. Also, 10 % noise was added.
GMSC The Give Me Some Credit (GMSC) data set 1 is a credit scoring data set where the objective is to decide whether a loan should be allowed or not. This decision is important for banks since erroneous loans lead to the risk of default and unnecessary expenses on future lawsuits. The data set contains historical data on 150,000 borrowers, each described by ten attributes.
Electricity The Electricity data set 2 was collected from the Australian New South Wales Electricity Market, where prices are not fixed. These prices are affected by the demand and supply of the market itself and set every five minutes. The Electricity data set contains 45, 312 instances, where class labels identify the changes in the price (2 possible classes: up or down) relative to a moving average of the last 24 hours. An important aspect of this data set is that it exhibits temporal dependencies [12].
Poker Hand The Poker Hand data set consists of 1,000,000 instances and eleven attributes. Each record of the Poker Hand data set is an example of a hand consisting of five playing cards drawn from a standard deck of 52. Each card is described using two attributes (suit and rank), for a total of ten predictive attributes. There is one class attribute that describes the ''Poker Hand'' .
This data set has no drift in its original form since the poker hand definitions do not change and the instances are randomly generated. Thus, the version presented in [5] is used, in which virtual drift is introduced via sorting the instances by rank and suit. Duplicate hands were also removed.
For comparison purposes, the configuration of synthetic stream generators is taken from [15]. In most real-world stream scenarios, we only have limited information about our data at the beginning. Thus, hyperparameter optimization cannot be done before creating the model. Hence, we do not optimize hyperparameters in our experiments. Note that in [16] hyperparameters of all classifiers were tuned, which lead to different results.
As a rule of thumb, we use 2 prototypes per class for RSLVQ classifiers; we use the default decay values for Adadelta [34] and Adamax [18]. For HAT, OBA, and SAMKNN, we used the default parameters as provided by the scikit-multiflow framework. Table 1 presents the data streams which are used in our experiment, as well as their configuration by the drift types. For a detailed description of the used data streams, see [15] and [21]. The tests of RSLVQ and GLVQ are performed with two prototypes per class.

Results
In the following, we present the experimental results. Please note that in all tables, the values are represented in percent, which means a Kappa score of 1 equals 100 % in the table. At the beginning, we only compare different LVQ versions against each other, to see which performs best. Table 2 shows the performance of the LVQ variants based on their accuracy. Based on the overall average, RSLVQ Adamax performs best. Also, the Adadelta version of RSLVQ performs significantly better than the vanilla version. The improvements are statistically significant based on a Friedman and Bonferroni-Dunn test. The better performance does not seem to be related to a specific drift type and can be seen on nearly every stream. However, the vanilla RSLVQ performs better than the adaptive versions on ELEC and POKR. Note that this was not the case when hyperparameters were tuned as in [16]. The results for GLVQ are very similar-the Adamax version is superior to Adadelta and vanilla GLVQ. However, the improvements are not significant here. Table 3 shows the results of the LVQ variants according to Kappa statistics. Again, RSLVQ Adamax performs best overall and on the synthetic streams. On the real-world streams, the vanilla RSLVQ has superior results to every other LVQ version and also to all other tested classifiers (Table 5). This is due to its good performance on ELEC and POKR, while the performance on the very imbalanced GMSC dataset is near 0. In summary, RSLVQ Adamax performed statistically significant better than the vanilla RSLVQ. Once again, the GLVQ behaves very similar, which means overall the adaptive versions achieve a better Kappa score, while being worse than its vanilla version on ELEC and POKR. Also, the improvements between GLVQ Adamax and GLVQ are not significant here.
Both tables show that the adaptive versions lack accuracy and Kappa on ELEC and POKR streams. This is probably caused by the fact that the decay does not make sense on a dataset which is not time dependent and thus not change over time like the cards of POKR. On all other datasets, the adaptive versions are superior; in particular, Adamax leads to remarkable improvements compared to the both vanilla LVQ variants.
In the next experiment, we compare our LVQ versions to other state-of-the-art streaming classifiers. We only use GLVQ and RSLVQ based on Adamax, due to its performance in the previous experiment. As non-LVQ classifiers, we use OzaBaggingAdwin (OBA), Hoeffding Adaptive Tree (HAT), and Self-Adjusting Memory K-Nearest Neighbor (SAMKNN), which are common stream classifiers [13,14]. Table 4 shows the results comparing our LVQ variants against state-of-the-art classifiers. SAMKNN is the best performing classifier, followed by OBA and GLVQ Adamax . On the synthetic data, RSLVQ Adamax is close above HAT but still behind OBA and SAMKNN. However, overall HAT and RSLVQ Adamax are very similar according to their accuracy. Also, GLVQ Adamax is only close behind its RSLVQ implementation. Note that while the results tend to  Table 5 shows the performance of the classifiers w.r.t. Kappa statistics. The results are as we would expect them from the accuracy. SAMKNN is the leading classifier with an overall Kappa of 0.76, followed by OBA and

Conclusion
In summary, the integration of Adadelta and Adamax into RSLVQ and GLVQ leads to improvements in prediction performance over their vanilla versions. In particular, RSLVQ Adamax performed best of all LVQ variants and is statistically significant better than RSLVQ. RSLVQ Adamax achieved a better performance than HAT, which is a very widely used stream classifier. Furthermore, the complexity of the RSLVQ variants is much lower than those of HAT. Hence, our adaptive RSLVQ variant is a suitable stream classifier which can be easily interpreted by its prototypebased characteristic. In further work, the performance is verified by using an ensemble of RSLVQ Adamax and a concept drift detector like done by other famous stream classifiers (e.g., Hoeffding Adaptive Tree). Additionally, momentum-based gradient descent should handle different drift types more explicitly, which should be addressed in subsequent work.