Adaptive basis functions for prototype-based classification of functional data

We present a framework for distance-based classification of functional data. We consider the analysis of labeled spectral data and time series by means of generalized matrix relevance learning vector quantization (GMLVQ) as an example. To take advantage of the functional nature, a functional expansion of the input data is considered. Instead of using a predefined set of basis functions for the expansion, a more flexible scheme of an adaptive functional basis is employed. GMLVQ is applied on the resulting functional parameters to solve the classification task. For comparison of the classification, a GMLVQ system is also applied to the raw input data, as well as on data expanded by a different predefined functional basis. Computer experiments show that the methods offer potential to improve classification performance significantly. Furthermore, the analysis of the adapted set of basis functions give further insights into the data structure and yields an option for a drastic reduction of dimensionality.


Introduction
In a variety of modern systems, high-dimensional sensor data are produced that is usually difficult to handle with traditional methods. To cope with the high number of input dimensions, the use of prior knowledge of the data generating systems can be used for simplification and dimensionality reduction. In many cases, high-dimensional data vectors can be assumed to represent discretized versions of continuous functions. Although time series might come to mind first as an example, similar properties are also present in other types of sequential data, such as spectral image data.
Such data are frequently recorded to serve as input for a classification task. Various machine learning algorithms can be applied, having specific advantages and disadvantages. The popularity of prototype-and distance-based classification systems results from their intuitive interpretation and straightforward implementation [2,7]. In this paper, an extension of the popular learning vector quantization (LVQ) [21] is used. LVQ systems comprise different prototypes which represent characteristic properties of their corresponding classes. Together with an appropriate distance measure, they constitute an efficient classification method.
The selection of a suitable distance measure is a key issue in the design of a prototype-based classifier [6]. While using predefined distances such as the Euclidean or other Minkowski-like measures often yields reasonable results, the application of more flexible relevance schemes has shown to be beneficial [8]. In generalized matrix relevance LVQ (GMLVQ) [34], a parameterized distance measure is employed whose parameters are determined in a data-driven training process [35]. This offers greater flexibility, since only the general structure of the distance measure has to be specified in advance. The analysis of the emerging distance measure can give further insights into the nature of the input data [34][35][36].
In this paper, we introduce a further extension of GMLVQ with respect to functional data. An expansion of functional data by means of suitable basis functions in combination with GMLVQ has proven to make superior classification performance and drastic reduction of input dimensions possible. In [26,27], the representation of spectral data and time series as weighted sums of basis functions facilitated an reduction of input dimensions. As a first example, Chebyshev polynomials of first kind were examined as possible basis functions. In [35], a wavelet transformation was employed for the dimensionality reduction, but was mainly motivated as a definition of discriminative features for sharply peaked mass spectra.
The use of functional data analysis for classification has been investigated in recent years in combination with different approaches including support vector machines [31], random forests [16] and self-organizing maps [30]. The increasing availability of deep learning frameworks and suitable hardware has further motivated the analysis of functional data by using convolutional neural networks [19] and auto-encoders [18]. While the latter studies follow a data-driven approach to incorporate the functional nature of the data, the interpretation of the obtained results remains tough, due to the high number of free parameters in deep learning models. Most of the other previous studies take only a single predefined generic functional expansion into account. The problem of an appropriate functional basis selection (with respect to the classification task) remains and is scarcely discussed in the literature.
Since the choice of the functional basis has direct influence on the resulting distance measure, it becomes a key problem in the design of functional data classifiers. Along the lines of the adaptive distance matrix in GMLVQ, we propose an adaptive functional basis that is optimized during the same data-driven training process. Please note that GMLVQ serves here as an example classification system and is favored based on its intuitive interpretation. The principle idea of an adaptive functional expansion can be transferred to other classification systems. This paper extends the approach of adaptive functional basis in combination with GMLVQ that was originally presented in the context of the WSOM? Conference 2017 [4].
In the following section, we outline the mathematical framework including the classification using GMLVQ and its extension toward an adaptive functional basis. The presented approach is then applied to several publicly available benchmark datasets of time series and spectral data. Results of computer experiments are discussed in Sect. 4. To validate the plausibility of the adapted functional basis, results are transferred to other machine learning algorithm in Sect. 5 before we conclude with a summary and outlook.

The mathematical framework
In a general supervised classification task, a set of P labeled d-dimensional feature vectors v l 2 R d and labels S l ¼ Sðv l Þ 2 f1; 2; . . .; Cg ð1Þ for l ¼ 1; 2; . . .; P is considered as a training set which is exploited for the adaptation of a classifier. As an example of prototype-and distance-based classification, we consider generalized matrix relevance LVQ in the following.

Generalized matrix relevance LVQ
In learning vector quantization (LVQ), the classifier is composed of a set of M labeled prototypes in combination with an appropriate distance measure dðw V ; vÞ [6]. An unlabeled feature vector v is assigned to the class of the closest prototype w Ã V with respect to the employed distance measure with dðw Ã V ; vÞ ¼ min m fdðw m V ; vÞg M m¼1 . The prototypes are adapted during a training process that is either of heuristic nature like the popular LVQ1 algorithm [21], or guided by a suitable cost function, e.g., [32]. Please note that the subscript V of the prototypes indicates that they are defined in the same space as the feature vectors v l . The assumption that they serve as representatives for their corresponding class facilitates an intuitive interpretation of the classifier with respect to the data.
The selection of an appropriate distance measure is a key issue in the design of an LVQ system [6]. Instead of using a fixed measure, like the Euclidean or other Minkowski-like distances, the use of flexible relevance schemes has shown to improve classification performance [17]. In generalized matrix relevance LVQ (GMLVQ) [34], a parameterized quadratic distance measure of the form is employed. The parametrization K V ¼ X > V X V ensures positive semi-definiteness of the measure [34]. To avoid numerical problems, we impose a normalization TrðK V Þ ¼ P i;j ½X V 2 i;j ¼ 1, as suggested in [34]. For the sake of simplicity, we restrict the definition here to a single global square distance measure X V 2 R dÂd . Extensions of the system to multiple local distance measures or rectangular matrices have been studied in [10].
Prototypes w V and the distance matrix X V are adapted during a single data-driven training process. This optimization can be implemented as a stochastic gradient descent, as in [34]. Here, we employ a batch gradient descent algorithm together with an automatic step size adaption scheme described in [28]. The optimization is guided by the cost function introduced in [32] For a given feature vector v l , w J V denotes the closest prototype of the correct class Sðw J V Þ ¼ Sðv l Þ, while w K V denotes the closest prototype of any other class Sðw K V Þ 6 ¼ Sðv l Þ. Often the function Uð. . .Þ in (4) is chosen to be sigmoidal. To avoid the introduction of further parameters, we resort to the simplest setting with UðxÞ ¼ x.
All results presented here are achieved using a Python implementation of GMLVQ, based on the code available at [5].

Adaptive basis for functional data
In contrast to a generic classification task, where feature vectors simply represent the concatenation of more or less independent quantities, functional input data are assumed to be a discretized representation of an unknown continuous function At this point, we restrict the discussion to regularly sampled functions with equidistant x j . An example for prototype-based classification of functional data that incorporate irregular sampling can be found in [26]. Given a suitable set of basis functions g k ðxÞ it is possible to find a parameterized representation of f l as resulting in an alternative feature vector c l . Limiting the number of coefficients to n makes c l 2 R n an approximative representation of f l ðxÞ. The parameters c l;k can be determined by a suitable optimization criterion such as the quadratic or maximum deviation. In [26], a matrix representation for the functional expansion was introduced as where the rows of W 2 R nÂd represent the employed basis functions evaluated at the sampling points x j .
Under consideration of the functional expansion of the input data, the distance measure introduced in (3) becomes Note that the relevance matrix K c 2 R nÂn is defined with respect to the space of functional parameters. Given this relation, the functional approximation matrix W can be optimized in a similar way as the relevance matrix K in the training. In contrast to a predefined fixed functional basis, this facilitates a more flexible framework. Not only it does eliminate the need of a prior basis selection, also the interpretation of the emerging functional basis can give insights into the input data structure with respect to the classification task.
To ensure a functional characteristic for the rows of W the optimization has to be constrained. We consider a matrix W to have a functional characteristic if its rows represent a set of basis functions. A desirable functional basis is a set of continuous, linear independent or even orthogonal functions; In our case, a discretized representation, respectively.
To achieve these properties, we propose two methods. As a first approach, we alter the update rule of W to incorporate the desired properties. We make use of a Gram-Schmidt orthonormalization [33] coupled with a row-wise moving average filter, resulting in a set of smooth orthogonal basis functions. During training the modified update is repeatedly applied as with a; b 2 ½0; 1 and a þ b 1 ð9Þ whereW represents the row-wise smoothened version of W and W Ã its orthonormalization. Using the parameters a and b the behavior of functional constraint is controlled, where a ¼ 0; b ¼ 0 voids the functional constraint completely, while a ¼ 1 only incorporates a smoothing of the basis functions and b ¼ 1 forces a full orthonormalization of the matrix. Previous studies [4] have shown that the results of proposed algorithm are robust over a wide range of values for a and b.
A second approach extends the GMLVQ cost function (4) to favor models with the desired functional properties. Two additional terms representing the discontinuity and non-orthogonality of W are added to the cost function as As a measure for the discontinuity dðWÞ represents the mean local (discrete) gradient of each basis function w i The orthogonality score oðWÞ is computed as Similar to the first approach, the parameters c and g control the behavior of the functional constraint. Both approaches are evaluated in the following.

Computer experiments
To validate the presented approaches, they are applied to four publicly available benchmark datasets comprising different kinds of functional data. The Tecator dataset [37] is a prominent benchmark dataset for spectral data classification. It contains 215 reflectance spectra of meat probes in the range from 850 to 1050 nm wavelength. Each of the spectra has 100 sampling points and is labeled according to the fat content (high/ low).
As a second dataset, we use the a subset of the Sugar dataset [23,24]. The Sugar dataset is a collection of spectral information from various sugar and sugar-related chemical compounds, acquired with different sensors. In particular, we use the data acquired with the FieldSpec sensor which offers the widest spectral bandwidth ranging from 350 to 2500 nm wavelength. It contains 450 spectra labeled to represent one of the nine compounds. The spectra are sampled at a resolution of 1 nm resulting in 2151 sampling points per spectrum.
In addition to these two spectral datasets, two time series datasets are selected from the UCR Time Series Repository [11], namely the Plane and Symbols data sets. Apart from the key properties of the datasets given in Table 1, no further information on the interpretation of the data or associated classification tasks is available at the repository.
For each of the datasets, we consider four alternative scenarios for the classification.
Scenario A To obtain a natural baseline for the classification performance in a first experiment, the functional nature of the data is ignored and it is used directly for the training of a GMLVQ system. Scenario B In a second set of experiments, the input data are expanded by fixed basis functions. To investigate and illustrate the influence of a proper choice Chebyshev polynomials of first kind (which yielded reasonable results in previous studies [25,27]) or Gaussian functions are employed. In both setups, the number of basis functions is varied as n ¼ 4; 8; 16; 32. For the Gaussian functions, the parameters are chosen as a k ¼ 1:0, b k ¼ À1 þ 2k=n and c k ¼ 0:25, what yields an even distribution of the functions with respect to the functional input data space. Scenario C As a third set of experiments, the proposed extension of GMLVQ introducing an adaptive functional basis W is considered using the altered update rule introduced in (9). It is applied every 64 optimization steps. In analogy to Scenario B the size of W and therefore the number of employed basis functions is varied n ¼ 4; 8; 16; 32. To investigate the influence of different initializations of W it is initialized as in Scenario B either by using Chebyshev polynomials of first kind or by distributing n Gaussian functions g k ðxÞ evenly over the input space. The control parameters for the update rule are kept constant a ¼ b ¼ 0:25, since previous studies [4] have shown minor influence on the results over a wide range of values.
Scenario D The final set of experiment resembles all experiments performed in Scenario C including all parameter variations but employs the altered cost function (10) instead of the modified update rule. The control parameters for the cost function are chosen as c ¼ g ¼ 0:25. Note that the parameters c and g serve as a weighting within the cost function. Assigning the same value to both parameters implies equal importance to orthogonality and smoothness of the emerging functional basis, which is a reasonable choice regarding the experiment.
Besides the different model setups, the same settings and parameters were used in the different scenarios. All GMLVQ systems comprised a single prototype per class. The distance matrix X was initialized proportional to the identity matrix and prototypes were initialized as the corresponding class conditional mean prior to the training.
Model parameters were optimized using a Python implementation of the Adam [20] optimizer included in the PyTorch library [1,29]. The step size for the optimization was initialized to 10 À3 . In each experiment 2048 optimization steps were performed, to ensure final convergence of the model parameters To avoid numerical problems, a row-wise normalization P d l¼1 ½W 2 k;l ¼ 1 was applied. Normalization of X was kept unchanged as described in Sect. 2.1.
In all experiments, the classification performance of the resulting GMLVQ systems is evaluated by means of the overall correct classification rate on a validation set. For the Sugar, Plane and Symbol datasets dedicated validation datasets are provided along with the data. Results for the experiments on the Tecator dataset are obtained by averaging over 16 random splits of the dataset into 60% training and 40% validation data. The achieved classification accuracies are shown in Table 2.
To highlight different benefits of the proposed method, Figs. 1, 2 and 3 depict additional details of the resulting GMLVQ systems for the Tecator dataset.

Discussion
The observed classification accuracy as given in Table 2 for scenarios B, C and D shows that a functional representation of the data has the potential to increase the classification performance significantly with respect to the naive approach which ignores the functional nature of the data (Scenario A). For all examined datasets, the classification accuracy achieved using a functional expansion of the data exceeds the one obtained by training in the original input feature space. For the datasets Tecator, Sugar and Symbols, the best observed accuracy was achieved using the adaptive functional expansion (Scenario C and D). The best accuracy for the Plane dataset was observed in Scenario B, where a fixed functional basis was employed. Nevertheless, the use of an adaptive functional basis achieved classification performances exceeding the one achieved in Scenario A. As already discussed in [27], the positive effect of the functional representation of the input data is not only due to its noise reducing character. Therefore, the results confirm that a functional encoding of the input data yields a representation more suitable for prototype-based classification.
Furthermore, Table 2 shows that the resulting classification performance achieved in Scenario C is sensitive to the initialization of W. The moving average filter introduced in the context of the functional constraint in (9) should void the influence of the initial basis choice. However, in the examined experimental setup, the functional constraint is applied every 64 optimization steps, so over the complete training of 2048 steps only 32 times. Therefore, the initialization of the functional expansion matrix W is not eroded completely over the training period. Scenario D overcomes this problem by directly incorporating the desired functional characteristics of W into an altered cost function, cf. (10). Since the smoothness constraint is taken into account in every optimization step, the influence of the initial values are voided comparatively fast. Therefore, the achieved classification performance is more stable with respect to the initialization of W.
Besides the classification accuracy, the interpretability of machine learning models is a key issue in the design of reliable classification systems. Prototype-based models are  often straightforward to interpret, since prototypes, which serve as characteristic exemplars of their class, are adapted in the original feature space. To validate the plausibility of the obtained models incorporating a functional transformation, we analyze the prototypes and relevance profiles of the resulting GMLVQ systems. In Fig. 1, the prototypes and the corresponding relevance profiles resulting from training on the Tecator dataset in the different scenarios are depicted. Please note that the relevance profile is obtained as R ¼ diagðK V Þ for Scenario A and R ¼ diagðW > K c WÞ for scenarios B, C and D, respectively. All obtained prototypes have a very similar shape.
The obtained relevance profiles, as illustrated in Fig. 1b, show more diversity with respect to the different scenarios. Note that the presented results were obtained using a set of only 4 basis functions. The limited capabilities of a small fixed functional basis to represent complex curves become clearly visible in the obtained relevance profiles of Scenario B. In the following scenarios the basis functions are adaptive, thus capable to represent more complex but data specific patterns even though only 4 basis functions are employed. All scenarios incorporating an adaptive functional basis create a relevance profile that is similarly interpretable as the one obtained in the original feature space. In particular, a concentration of relevance around feature index 40 is clearly visible in all obtained profiles. This observed maximum has been associated with a specific peak of spectral absorbance associated to fat alone in previous studies [30]. The relevance profiles resulting from Scenario D resemble the original relevance profile best. Not only does this emphasize the conversation of model interpretability, it also backs the reduced sensitivity with respect to the initialization of W discussed before.
Please note that the analysis of model interpretability here assumes prototypes and relevance profile of Scenario A to be correct. A detailed discussion that the prototypes and relevance profiles reflect the underlying biochemistry most faithfully is beyond the scope of this contribution.
Since the functional basis (in Scenarios C and D) is adapted during training, it incorporates specific features of the data. In other words, we expect the functional basis not only to yield a good representation of the data, but also to give further insights into it. To follow this idea, it is essential to understand that the presented approach can be separated into three different aspects. During the training phase, representative examples for each class are determined (prototypes) while a distance measure is optimized. The distance measure itself is separated into the adaptation of a functional basis and a weighting of the different basis functions realized by the relevance matrix K.
As an illustrative example, the relation between the adapted basis functions and their particular relevance is given in Fig. 2. As a first interesting point, we observe the small relevance assigned to the first basis function which is represented by a solid line in Fig. 2a. This particular function is approximately constant and thus can be understood as an offset in the input data. A small relevance implies that this offset present in the data plays a minor role for the classification task, which confirms the results of previous studies [25]. The basis function with the highest relevance according to Fig. 2b is the second basis functions, which is depicted as a dashed line in Fig. 2a. This basis function clearly shows a linear trend superimposed by a local maximum centered at feature index 40. In accordance with the previous discussion, this emphasizes the importance of the fat-related spectral peak once more.
As stated before, the resulting classification performance is influenced by the initialization of W. This motivates the question if the interpretation of the most relevant basis function for the Tecator dataset by means of the spectral characteristics of fat is based on a lucky coincidence. Figure 3 depicts the two most relevant basis functions obtained from the different combinations of W initialization and adaptive learning rule. Although the specific shape of each set of basis functions is rather unique, the most relevant basis function displays the characteristic fat peak at feature index 40 in all experiments. The presence of this most prominent feature (with respect to the classification task) together with the high assigned relevance suggests that all of the obtained basis function sets are plausible.
Besides the superior classification accuracy and enhanced possibilities of model interpretation, the presented approach also reduces the computational effort for the model generation. The introduction of a rectangular matrix W reduces the number of free parameters for the classification model drastically. While there are d 2 free parameters for a simple quadratic X in GMLVQ, with the presented approach the number reduces to n Á d þ n 2 . Note that this is only beneficial for n\dð ffiffi ffi 5 p À 1Þ=2 % 0:62d. Hence, for a reasonable number of basis functions the advantage of fewer free parameters persists. This does not only reduce the computational effort of the optimization process but also reduces the risk of overfitting and convergence problems.

Generalization of the learned basis
The proposed method to find a data-and problem-specific functional basis is tightly coupled to the optimization of a GMLVQ system. The observed classification performance benefits from the compact representation of the functional data by the means of abstract approximation coefficients. Nevertheless, the question, if this beneficial effect is caused by the problem-specific functional basis or results from the combination with GMLVQ, remains.
To investigate this issue, we perform additional experiments in which the adapted functional basis is used in combination with other machine learning approaches. The original input data from the dataset are transformed to a new set of feature vectors according to (7) that is then used as input to the training of a classifier. We restrict the discussion to the Sugar dataset as an illustrative example for the generalization of the adapted functional basis, since it has shown the most impressive reduction of input dimensionality while preserving or even enhancing the achieved classification performance in the previous experiments.
A wide variety of classifiers is available each having different advantages and disadvantages. To cover a broad range, we selected four basic algorithms, namely k-nearestneighbor (kNN) [14], multilayer perceptron (MLP) [22], linear support vector machine (SVM) [12] and decision tree (Tree) [9] as illustrative examples. Please note that an in-depth description of the mathematical background for all of the employed classifiers goes beyond the scope of this document. For details on the training and classification algorithms consult the according references.
For each of the classification algorithms, an additional set of computer experiments is performed. The first experiment serves, similar to the previous setup, as baseline for the classification performance and trains the classifier using the original input data without further preprocessing. In a second set of experiments, the input data are expanded by a predefined set of basis functions, namely Chebyshev polynomials of first kind and Gaussian functions, equivalent to Scenario B in the previous experiments. The number of basis functions is also varied as n ¼ 4; 8; 16; 32. In a final set of experiments, the adapted functional basis from the experiments in Scenario D in Sect. 3 is employed to transform the input data into problem-specific approximation coefficients that are then used to optimize the classification system. The problem-specific expansion coefficients are obtained by transformation of the original feature vectors according (7). The employed transformation matrices W result from the corresponding experiments in Scenario D of Sect. 3 with Chebyshev polynomial initialization. As in the previous experiments the size of the functional basis is varied as n ¼ 4; 8; 16; 32. Please note that the transformation matrix W serves in these experiments as a fixed functional basis and is not altered by any means during the training of the classifier.
Apart from the different representations of the input data, all experiments use the same set of parameters. For the kNN classifier k ¼ 1 together with the squared Euclidean distance serves as the simplest classifier possible. The MLP model employs a single hidden layer containing 16 neurons with sigmoidal activation. Since SVMs are only suited for binary classification problems, an ensemble of linear kernel SVMs is trained using error-correcting output codes [13] to map the nine different classes to multiple binary classification tasks. A standard classification and regression tree (CART) algorithm introduced in [9] employing Gini's diversity index [15] is used to construct the decision tree classifier.
All other parameters for the classification are kept at the default values specified in the corresponding MAT-LAB TM toolboxes. The used implementation for kNN, SVM and decision tree classification is provided by the Statistics and Machine Learning Toolbox (v10.2), while the MLP is provided by the Neural Network Toolbox (v9.0).
Equivalent to the previous experiments the classification performance of the different algorithms is measured by the accuracy with respect to the validation set provided with the input data. All experiments are repeated 64 times to account for random model initializations and the mean accuracy over all experiments is evaluated. Table 3 shows the observed accuracies. The table results indicate that a superior classification performance is achieved for all examined classifiers when a functional representation of the data is used. The use of a generic functional basis can be beneficial, but the performance is highly dependent on the selected basis function as well as the classification algorithm. While the use of Chebyshev polynomials of first kind can be beneficial for the kNN, Tree and ANN classifier, the achieved accuracy is significantly harmed for the SVM classification system. Even worse, all of the evaluated classification algorithms showed poor performance using simple Gaussian basis functions. This observation confirms the hypothesis that Gaussian basis functions are not suitable for a classification focused approximation of the sugar dataset.
Using the learned problem-specific functional basis outperforms the classification in the original input space in all experiments. The constant superior performance with respect to the classification in the original input space proves that the functional basis provides a meaningful, low-dimensional and problem-specific feature space which is not limited to use of GMLVQ classifiers. The drastic reduction in input dimensions (here d ¼ 2151 ! n ¼ 8) reduces the complexity of the resulting classifiers significantly.

Summary and outlook
This contribution proposes an extension of GMLVQ for functional data that takes the special characteristics of the input data into account by using an adaptive functional expansion. We have shown that this approach has the potential to improve overall classification accuracy. For all of the examined benchmark datasets, models using an adaptive functional basis outperformed the classification accuracy achieved on the raw input data.
Furthermore, the interpretation of the adapted functional basis can give additional insights into the structure of the input. For the Tecator dataset, one of the learned basis functions can be interpreted as a data offset, while other basis functions can be associated with the spectrum of fatty tissue. The experiments have shown that the presented approach has the potential to decompose signal into different compounds with respect to the classification task. These findings may motivate a new approach for the recovery of spectral end-members that overcomes the assumption that end-members have to be present in the data, like in [3].
In combination with the relevance scheme used in GMLVQ, the importance of different data characteristics with respect to the classification task is weighted. For the Tecator example, the data offset is almost ignored with respect to the classification task, while the basis function associated with the fat spectrum is assigned the highest relevance.
In addition to the separate interpretation of functional basis and assigned relevances, it is possible to map the model parameters to the original feature space, which facilitates a direct analysis of relevances in the space of input data. Due to the limited number of basis functions employed, a perfect reconstruction is not possible and the back projected relevances (and relevance matrices) can only serve as an estimation for the input space relevances. Nevertheless, our experiments have shown that the general shape of the relevance profiles is preserved. Therefore, the presented approach does not only allow for the interpretation of model parameters in the space of functional coefficients, but preserves the interpretability of prototypes and the relevance matrix in the original input space. Additionally, the introduction of the rectangular matrix W drastically decreases the number of free model parameters. This leads to a reduced computational effort for the model generation and a lower risk to encounter overfitting effects and convergence problems.
In additional experiments, the use of the adapted functional basis in combination with different classification algorithms, namely k-nearest-neighbor, support vector machines, decision trees and artificial neural networks, was investigated. All classifiers have shown superior performance to the brute-force approach when trained on the data projected by the adaptive basis. In contrast, the classification performance in combination with generic basis functions could not significantly exceed the one on the original input data. This proves that the adapted set of basis All accuracies that exceed the accuracy achieved on the original input data are highlighted in bold face. The highest observed accuracy for each classification algorithm is additionally emphasized in bold italic functions provides a meaningful, data-and problemspecific domain with respect to the classification task. Still there are various questions remaining to be addressed in future research. For the introduction of the adaptive functional basis, the number of basis functions has to be specified in advance. The performed computer experiments have shown that the approach is sensitive to the size of the functional basis and an increase in the number of basis functions can lead to a decrease in classification performance. A consistent procedure to estimate the number of necessary basis functions is desirable. Furthermore, the experiments have shown that different initial values for W as well as different optimization criteria lead to significant differences in the adapted functional basis. The question if there is an unique optimal set of basis functions, given a set of training data and a classification task, remains. For the experiments performed in the scope of this work, the optimization was canceled after a defined number of optimization steps. It is possible that the convergence of W toward an unique optimal solution was therefore aborted too early.
Finally, the presented results are obtained on a comparatively low number of example datasets and classification algorithms. With the increasing availability of annotated functional datasets, the proposed approach should to be further validated.