Enhanced Data Representation by Kernel Metric Learning for Dementia Diagnosis

Alzheimer's disease (AD) is the kind of dementia that affects the most people around the world. Therefore, an early identification supporting effective treatments is required to increase the life quality of a wide number of patients. Recently, computer-aided diagnosis tools for dementia using Magnetic Resonance Imaging scans have been successfully proposed to discriminate between patients with AD, mild cognitive impairment, and healthy controls. Most of the attention has been given to the clinical data, provided by initiatives as the ADNI, supporting reliable researches on intervention, prevention, and treatments of AD. Therefore, there is a need for improving the performance of classification machines. In this paper, we propose a kernel framework for learning metrics that enhances conventional machines and supports the diagnosis of dementia. Our framework aims at building discriminative spaces through the maximization of center kernel alignment function, aiming at improving the discrimination of the three considered neurological classes. The proposed metric learning performance is evaluated on the widely-known ADNI database using three supervised classification machines (k-nn, SVM and NNs) for multi-class and bi-class scenarios from structural MRIs. Specifically, from ADNI collection 286 AD patients, 379 MCI patients and 231 healthy controls are used for development and validation of our proposed metric learning framework. For the experimental validation, we split the data into two subsets: 30% of subjects used like a blindfolded assessment and 70% employed for parameter tuning. Then, in the preprocessing stage, each structural MRI scan a total of 310 morphological measurements are automatically extracted from by FreeSurfer software package and concatenated to build an input feature matrix. Obtained test performance results, show that including a supervised metric learning improves the compared baseline classifiers in both scenarios. In the multi-class scenario, we achieve the best performance (accuracy 60.1%) for pretrained 1-layered NN, and we obtain measures over 90% in the average for HC vs. AD task. From the machine learning point of view, our proposal enhances the classifier performance by building spaces with a better class separability. From the clinical application, our enhancement results in a more balanced performance in each class than the compared approaches from the CADDementia challenge by increasing the sensitivity of pathological groups and the specificity of healthy controls.

1. Introduction

Alzheimer's Disease (AD) corresponds to a progressive cognitive impairment and loss of memory functions, becoming the kind of dementia with the largest prevalence in elderly subjects with nearly 44 million patients worldwide. Besides, Mild Cognitive Impairment (MCI), a previous stage of AD, affects 10–20% of people aged 65 or older (Alzheimer's Association, 2016). Hence, dementia diagnosis and treatment demand reliable biomarkers providing an objective and early characterization of the different AD stages (Shi et al., 2015). Among these biomarkers, structural magnetic resonance imaging (MRI) data became frequently used to develop computer-aided diagnosis (CAD) tools due to its wide availability and non-invasiveness (Jack et al., 2013). CAD tools learn to discriminate dementia through MRI features, benefiting from large amounts of neuroimaging data.

Particularly, the Alzheimer's Disease Neuroimaging Initiative (ADNI) focuses its researches on discriminating pathologies with a variety of classification tools from neuroimaging data, genetic information, and other biomarkes. However, insufficient attention has been given to build appropriate metrics from the training data that could maximize the performance of several classifiers (Shi et al., 2015).

In this regard, learning distance metrics tuned for classification tasks from given prior information, known as Metric Learning (ML), implicitly transforms input features into discriminative ones. As an example, linear dimensionality reduction changes the metric in the original space to maximize an objective function (Fukumizu et al., 2004). Therefore, introducing a metric learning stage into the discrimination process can significantly improve the performance of distance-based machines as k-nn, k-means, and even SVMs (Xu et al., 2012). Depending on the transformation to be sought, ML can be divided into linear and nonlinear models.

Linear models aim at estimating an optimal affine transformation for the input space, with straightforward interpretation, optimization simplicity, and reduced overfitting (Bellet et al., 2013). Among these models, principal component analysis (PCA) finds the subspace best preserving the variance of the data. Wachinger et al. (2016) stacks PCA matrices and logistic regressors in a multi-layer architecture. However, the generative properties of the resulting machine highlight over the discriminative ones. In the same way, Khedher et al. (2015) proposes a dimension reduction approach based on performing an independent component analysis (ICA) over high-dimensional voxel-wise features. Despite improving the accuracy of an SVM classifier, the maximization of statistical independencence of ICA may eliminate relationships between the dependent features that hold discriminative information, and restricts the resulting dimension to the number of involved classes. In addition, when handling data distributions with nonlinear structures, linear models show inherently limited performance and class separation capability (Shi et al., 2015).

On the other hand, the most popular nonlinear models are built through kernel-based methods, that embed input features into higher dimensional spaces with an increased linear separability. For instance, Awate et al. (2017) develop a kernel-based framework for building cortex descriptors that detect population differences and regress clinical variables. Specifically for AD diagnosis, Zhang et al. (2011) combines three different biomarkers using a simple-while-effective multiple-kernel learning for improving the SVM-based classification of AD and MCI. However, optimization of kernel weighting is carried out by a grid search, which is very time consuming when the number of features and samples gets large (Liu et al., 2014). Another approach uses multiple-kernels, each of them representing anatomically meaningful brain regions, to improve AD discrimination (Young et al., 2015). The kernel combination weights, that are computed via Gaussian Process framework, can be clinically read. Despite enhancing the performance of linear methods for many highly nonlinear problems, these kind of kernel-based solutions are prone to over-fitting (Bellet et al., 2013), and their utilization is inherently limited by the sizes of the kernel matrices (He et al., 2013). Another kind of non-linear models combine multiple local metrics to learn a feature transformation based on either local neighborhoods or class memberships (Wang et al., 2012). Although these multi-metric strategies are usually more powerful in modeling nonlinear structures, generalizing these methods to train robust classifiers is not trivial. In addition, the estimation of high-level complexity distances on such metric manifolds also yields to high computational cost (Shi et al., 2015).

In this paper, we introduce a kernel-based metric learning framework for supporting the dementia diagnosis task. The proposed approach searches for projections into more discriminative spaces so that the resulting data distribution resembles as much as possible the label distribution. Hence, we incorporate kernel theory for assessing the affinity between projected data and available labels through the Center Kernel Alignment (CKA) criterion. Maximization of CKA produces discriminative features aimed at improving the discrimination neurological classes. Furthermore, the proposed metric learning is introduced into three commonly used supervised classification machines, namely, k-nearest neighbors (k-nn), support vector machines (SVM) and neural networks (NN). For the latter, we generalize the supervised linear projection as a multilayer architecture. Our approach is evaluated on two scenarios of dementia diagnosis from structural MRIs (multi-class and bi-class). To this end, we use morphological measurements (volume, area, and thickness) computed by the widely used FreeSurfer suite (Tustison et al., 2014). Attained test classification results of several performance measures show that introducing a supervised metric learning improves other the baseline classifiers in both considered scenarios. Therefore, the proposal builds spaces with a better class separability that increases both the sensitivity and specificity, with the additional benefit of a more balanced the performance for each class.

The agenda of this paper is organized as follows: Firstly, we describe the mathematical framework and our proposed approach in Section 2. Then, Section 3 illustrates the carried out evaluations on the well-know ADNI collection for tuning the parameters of the CKA and classifiers. Lastly, we describe all obtained results and discuss concluding remarks in Section 4 and 5, respectively.

2. Materials and Methods

2.1. Centered Kernel Alignment

Kernel functions are bivariate measures of similarity, which are based on the inner product between samples embedded in a Hilbert space. For an input feature space ${X}$, a kernel ${\kappa }_X:{X}\times {X}\to {ℝ}^{+}$ is a positive-definite function that defines an implicit mapping ${\phi }_X:{X}\to {{H}}_X$, aiming to embed a data point x ∈ ${X}$ into the element φ_X(x) ∈ ${H}$_X of some Reproducing Kernel Hilbert Space (RKHS) (noted as ${H}$_X). Within a supervised learning framework, a kernel ${\kappa }_L:{L}\times {L}\to {ℝ}^{+}$ is also introduced that acts over the target space ${L}$ to account for the attribute labeling information so that κ_L defines the implicit mapping φ_L(l):${L}$→${H}$_L. Due to each function (κ_X and κ_L) reflects a different notion of similarity extracted from a distinct sample set, the concept of alignment between mappings can be introduced to measure the degree of agreement between the input and target kernels. To unify both tasks into a coherent optimization problem, we employ the Centered Kernel Alignment (CKA) that assesses the kernel affinity through the expected value of their normalized inner product over all data points as follows:

$\begin{array}{cc}\rho \left({\kappa }_X,{\kappa }_L\right)=\frac{𝔼{}_{x{x}^{\prime }l{l}^{\prime }}\left\{{\stackrel{-}{\kappa }}_X\left(x,x^{\prime }\right){\stackrel{-}{\kappa }}_L\left(l,l^{\prime }\right)\right\}}{\sqrt{{𝔼}_{x{x}^{\prime }}\left\{{\stackrel{-}{\kappa }}_X^2\left(x,x^{\prime }\right)\right\}{𝔼}_{l{l}^{\prime }}\left\{{\stackrel{-}{\kappa }}_L^2\left(l,l^{\prime }\right)\right\}}},& (1)\end{array}$

where notation 𝔼_z{·} stands for the expected value of the random variable z, ${\stackrel{-}{\kappa }}_Z(z,z^{\prime })$ is the centered version of the kernel function ${\kappa }_Z(z,z^{\prime })={({\phi }_Z(z)-{\stackrel{-}{\phi }}_Z)}^{\top }({\phi }_Z(z^{\prime })-{\stackrel{-}{\phi }}_Z)$, being ${\stackrel{-}{\phi }}_Z\in {{H}}_Z$ the expected value of the data distribution on ${H}$_Z.

In practice, the characterizing kernel matrices, $K_X\in {ℝ}^{N\times N}$ and $K_L\in {ℝ}^{N\times N}$, are extracted from a provided input dataset X ∈ ℝ^{N × P}, holding samples $x_n\in {ℝ}^P$, along with its corresponding target vector $l=\left\{l_n\subset \mathbb{Z}:n\in \left[1,N\right]\right\}\in {\mathbb{Z}}^N$. Hence, the empirical estimate for the CKA value can be computed as follows:

$\begin{array}{cc}\widehat{\rho }\left({\stackrel{-}{K}}_X,{\stackrel{-}{K}}_L\right)=\frac{{\langle {\stackrel{-}{K}}_X,{\stackrel{-}{K}}_L\rangle }_F}{\sqrt{{\langle {\stackrel{-}{K}}_X,{\stackrel{-}{K}}_X\rangle }_F{\langle {\stackrel{-}{K}}_L,{\stackrel{-}{K}}_L\rangle }_F}},& (2)\end{array}$

where notation 〈·, ·〉_F stands for the matrix-based Frobenius product, and $\stackrel{-}{K}={({\phi }_Z-{\stackrel{-}{\phi }}_Z)}^{\top }({\phi }_Z-{\stackrel{-}{\phi }}_Z)$ is the centered kernel matrix (associated with ${\stackrel{-}{\kappa }}_Z(,)$) computed as $\stackrel{-}{K}=\stackrel{-}{I}K\stackrel{-}{I}$, being 1 ∈ ℝ^{N × 1} the all-ones vector, I the identity matrix, and $\stackrel{-}{I}=\left[I-{\mathrm{11}}^{\top }/N\right]$.

The inner product of both kernel functions in Equation (2) estimates the agreement between the data jointly sampled from the spaces ${X}$ and ${L}$ through their statistical dependence ρ ∈ [0, 1]. Therefore, the larger the value of the CKA, the more similar the distributions of the input and target data. As a result, using ρ as the cost function to maximize the dependence between the characterizing kernels of the input data and the class label promotes the class discrimination from the given input samples (Gretton et al., 2005).

2.2. Supervised Metric Learning for Classification

The CKA dependence using the Mahalanobis metric learning is developed for three commonly used approaches to supervised classification: k-nn and SVM that are fed into a Gaussian kernel optimization, and NN that are optimized by a backpropagation strategy with cross-entropy as the cost function.

2.2.1. Gaussian Kernel Optimization for Classification

In general, the Gaussian kernel is preferred in pattern classification applications since it aims at finding an RKHS with universal approximating ability, not to mention its mathematical tractability. Nonetheless, to account for the variance of each space when measuring the pairwise distance between samples x_n and $x_{n^{\prime }}$, the Gaussian kernel relies on the generalized Euclidean metric that is parameterized by a linear projection matrix W in the form:

$\begin{array}{cc}{\kappa }_X\left(W,\sigma \right)=\exp \left(-\left(x_n-x_{n^{\prime }}\right)W{W}^{\top }{\left(x_n-x_{n^{\prime }}\right)}^{\top }/2{\sigma }^2\right)& (3)\end{array}$

where σ ∈ ℝ⁺ is the kernel bandwidth that rules the observation window within the similarity distance is assessed.

In terms of the projection matrix, the formulation of the CKA-based optimizing function in Equation (2) can be integrated into the following kernel-based learning problem:

$\begin{array}{cc}\begin{array}{l}\widehat{W}=\arg \underset{W}{\min }\left\{-\log \left(\widehat{\rho }\left({\stackrel{-}{K}}_X(W),{\stackrel{-}{K}}_L\right)\right)\right\},\\ =\arg \underset{W}{\min }\{\log \left(\text{tr}\left(K_X(W)\stackrel{-}{I}{K}_L\stackrel{-}{I}\right)\right)\end{array}& (4\text{a})\end{array}$

$\begin{array}{cc}-\frac{1}{2}\log \left(\text{tr}\left(K_X(W)\stackrel{-}{I}{K}_X(W)\stackrel{-}{I}\right)\right)\},& (4\text{b})\end{array}$

where the logarithm function is used for mathematical convenience. Therefore, the first term in Equation (4b) assesses the similarity between input and target kernels while the second one works as a regularization term minimizing the norm of the input kernel.

2.2.2. Network Pretraining for NN Classification

From Equation (3), it holds that the projection matrix defines a new space linearly dependent on the input features, that is, y = xW. Consequently, the generalized Euclidean metric also defines in the new space the following kernel function:

$\begin{array}{cc}{\kappa }_Y(y,y^{\prime })={\kappa }_X(xW,x^{\prime }W)& (5)\end{array}$

The above metric learning can be read as a particular case of a multi-layer architecture, having a single layer and a linear transfer function. This type of architectures sequentially unravels input data through the stacked linear/non-linear mappings aiming to estimate the output data. Therefore, an M-layered Artificial Neural Network (ANN) can be employed to predict the label of an input sample through a battery of feed-forward deterministic transformations, which are implemented by the hidden layers h^m, mapping the input sample x to the network output h^M as below (Bengio, 2009):

$\begin{array}{cc}\{\begin{array}{l}{h}^m=\vartheta (b^m+h^{m-1}{W}^m),\forall m=1,\dots ,M-1\hfill \\ h^0=x\hfill \end{array}& (6)\end{array}$

where $b^m\in {ℝ}^{Q_{m+1}}$ is the m-th offset vector, $W^m\in {ℝ}^{Q_{m+1}\times Q_m}$ is the m-th linear projection, $Q_m\in {\mathbb{Z}}^{+}$ is the size of the m-th layer, and ϑ(·) ∈ ℝ is the transfer function that non-linearly maps, saturating each scalar input element. Besides, the first and last layers of the network are always fixed to the input and output dimensions of the learning task, that is, h⁰ ∈ ℝ^P and h^M ∈ ${L}$⊂{0, 1}^C with C mutually exclusive classes.

It is worth noting that the architectures in Equation (6) attempt to encode into each hidden space the most discriminative attributes of the available dataset, managing the layers. In this line of analysis, we propose a pre-training procedure to provide a suitable initial set of projections {W^m}. To this end, the matrix $H^m\in {ℝ}^{N\times Q_m}$ is defined to hold every m-th non-linear projection of the provided input data X, assuming H⁰ = X. As a result, all pairwise similarities between samples, at m-th layer, are gathered into the kernel matrix $K_m\in {ℝ}^{N\times N}$ with elements computed as ${\kappa }_m(d(h_n^m,h_{n^{\prime }}^m))$, being $d:{ℝ}^{Q_m}\times {ℝ}^{Q_m}\to {ℝ}^{+}$ a distance operator defined as follows:

$\begin{array}{cc}d\left(h_{n^{\prime }}^m,h_{n^{\prime }}^m\right)=\left(h_n^m-h_{n^{\prime }}^m\right){\left(h_n^m-h_{n^{\prime }}^m\right)}^{\top }& (7)\end{array}$

where every $h_n^m\in {ℝ}^{Q_m}$ holds each linear projection that we choose of a saturated class, that is, $h_n^m=\vartheta (h_n^{m-1}{W}^m)$.

With the aim of making our proposal suitable for supervised learning tasks, we introduce the pre-training of W^m at each layer along with the provided prior knowledge into the CKA optimization framework of Equation (4b). Hence, each matrix W^m is learned by maximizing the similarity between K_m and K_L through the real-valued function that evaluates the alignment between both kernels in Equation (2). Thus, the pre-trained set $\left\{{\hat{W}}^m\right\}$ that initializes the network layers provides an assembly of discriminative linear projections, matching the most the distribution of the projected data H^m and target information l.

For both the metric learning and network pre-training, the parameters of the kernel-based learning problem defined in this section are optimized in order solve the problem in Equation (4b). Since, the kernel-based objective is a nonlinear function of the parameters, we use the iterative stochastic gradient method for optimization (Brockmeier et al., 2014). In particular, the CKA optimization procedure is initialized from the principal component mapping of the data, which meets orthogonality conditions. Nonetheless, the span of the resulting projected space is not necessarily orthogonal since there are no restrictions on optimization. Resulting from the optimization, the proposed transformation applies any linear mapping of the data so that the classes are easier to identify. In the case of the metric learning, the projection matrix can be seen as an affine transformation where classes are spatially separated but the nonlinear properties are held. In the case of the pre-training, a nonlinear mapping unravels the data distribution to better classify the samples.

3. Experimental Set-up

This work introduces a metric learning framework, based on kernel alignment, to support classification tasks of computer-aided diagnosis of dementia. Specifically, validation is accomplished by grouping structural Magnetic Resonance Images (MRI) into one of the following three neurological classes: Normal Control (NC), Mild Cognitive Impairment (MCI), and Alzheimer's Disease (AD). Figure 1 outlines the training pipeline, where the highlighted box contains the proposed data representation stage based on CKA that can be introduced into three considered classification machines.

FIGURE 1

Figure 1. Schematic training pipeline: Input MRI is automatically segmented and a set of features are further extracted. CKA learns a metric or pre-trains network weigths. Tuned classifiers (k-nn, SVM, and ANN) are used for computer-aided diagnosis task.

3.1. MRI Database and Preprocessing

The proposed metric learning is evaluated on a set holding 3,304 MRIs obtained from 896 adult patients selected from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database¹. This data collection aims at measuring the progression of dementia from serial MRIs, positron emission tomography, other biosignals, and neuropsychological assessments. Table 1 briefly describes the demographic information of the tested cohort of ADNI MRIs, which is split into two subsets. The first one holds 30% of the subjects and is devoted to a blindfolded assessment of the performance framework. The remaining 70% of subjects is employed for framework parameter tuning, which is carried by a 5-fold cross-validation scheme to guarantee that all images of the same subject are assigned to a single group of data analysis (i.e., a validation fold or the test subset).

TABLE 1

Table 1. Demographic details of the selected ADNI dataset cohort, distributed over the three classes.

The tested set of structural MRIs if automatically pre-processed via the widely used FreeSurfer software package² that computes the needed morphological measurements with suitable test-retest reliability across scanner manufacturers and field strengths. Firstly, the bias field of the structural images is corrected as well as its intensity is normalized. Further, the gray and white matter are segmented, and the brain cortex is parcellated. The shapes of resulting structures are further tessellated. Then, structure-wise thickness, area, and volume measurements are computed. The volume features are normalized to the Total Intracranial Volume (Buckner et al., 2004). As a result, an input feature matrix X with size N = 3304 and P = 310 is built using the features from each MRI as summarized in Table 2.

TABLE 2

Table 2. Morphological features extracted by FreeSurfer.

3.2. Parameter Setting and Network Topology

As regards the CKA implementation, the kernel function is computed as below:

$\begin{array}{cc}\begin{array}{l}{\stackrel{-}{\kappa }}_Z\left(z,z^{\prime }\right)={\kappa }_Z\left(z,z^{\prime }\right)-{𝔼}_{z^{\prime }}\left\{{\kappa }_Z\left(z,z^{\prime }\right)\right\}-{𝔼}_z\left\{{\kappa }_Z\left(z,z^{\prime }\right)\right\}\\ +{𝔼}_{z{z}^{\prime }}\left\{{\kappa }_Z\left(z,z^{\prime }\right)\right\},\forall z,z^{\prime }\subset Z\end{array}& (8)\end{array}$

Therefore, provided a finite sample {(x_n, l_n):n ∈ N}, elements ${\stackrel{-}{k}}_{n{n}^{\prime }}={\stackrel{-}{\kappa }}_Z(z_n,z_{n^{\prime }})$ in Equation (8) are estimated by the following centered kernel function:

$\begin{array}{cc}{\stackrel{-}{k}}_{n{n}^{\prime }}=k_{n{n}^{\prime }}-{𝔼}_n\left\{k_{n{n}^{\prime }}\right\}-{𝔼}_{n^{\prime }}\left\{k_{n{n}^{\prime }}\right\}+{𝔼}_{n,n^{\prime }}\left\{k_{n{n}^{\prime }}\right\}.& (9)\end{array}$

In the validation stage, we consider two scenarios of diagnosis, namely, multi-class task (i.e., MCI, HC, and AD) and three bi-class tasks: HC vs. MCI, HC vs. AD, and MCI vs. AD. In all cases, each compared classifier is adjusted to the following parameter setting:

– k Nearest Neighbors (k-nn): The number of neighbors is tuned by heuristic searches according to the best performance reached within the specific range [1, 3, 5, 7, 9], resulting in k ={7, 9} as the optimal parameter for the baseline and metric learning approaches, respectively.

– Support Vector Machine (SVM): The regularization parameter of the soft margin is searched within the range C∈[0.1, 1, 10, 100, 1000]. The SVM-based classifier is trained by the sequential minimal optimization solver for a Gaussian kernel with bandwidth tuned following the algorithm proposed by Álvarez-Meza et al. (2014) that achieves the highest accuracy at C = 10 for the baseline and C=0.1 for the metric learning approach.

– Neural Networks (NN): Architectures with one and two hidden layers (M = {1, 2}) with m₁ = m₂ = m are considered for all experiments. Activation functions of the hidden states are the standard sigmoid, ϑ(z) = (tanh(z) + 1)/2, while at the output layer the following function is selected:

$h_c^M=\exp (b_c^m+w_c^M{h}^{M-1})/{\displaystyle {\sum }_j\{\exp (}{b}_j^m+w_j^M{h}^{M-1})\}$

where $b_c^M$ is the c-th element of b^M, $w_c^M$ the c-th row of W^M, and $\sum _c{h}_c^M=1$. Also, we tune the hidden layer size using the exhaustive search strategy for the highest reached accuracy. To this end, the NN classifier is iteratively trained through a range of log-spaced layer sizes bounded by m ∈ [C, P−1]. The exhaustive search for the layer size attained m₁ = 11 for 1-layered network, either randomly initialized or with supervised pre-training. The 2-layered architecture obtains the highest accuracy at m₁ = m₂ = 11 when randomly initialized and m₁ = m₂ = 5 when pre-trained.

The diagnosis performance is assessed by the classification accuracy a and true positive fraction for the corresponding class (τ_c), defined as follows:

$\begin{array}{cc}a={𝔼}_c\left\{t_c^p+t_c^n\right\}& (10\text{a})\end{array}$

$\begin{array}{cc}{\tau }_c=t_c/N_c,c=\{MCI,HC,AD\}& (10\text{b})\end{array}$

being $t_c^p$, $t_c^n$, and N_c the number of true positives, true negatives, and total samples for c-th class, respectively. Note that the true positive fractions provide information about the sensitivity and the specificity. Thus, for the multi-class case, τ_AD and τ_MCI can be interpreted as the two-class sensitivity, while τ_NC agrees with the two-class specificity (Bron et al., 2015).

4. Results and Discussion

For enhancing the computer-aided diagnosis of dementia, we explore a metric learning framework based on the centered kernel alignment function, aiming to estimate more discriminative spaces. Evaluation prove the benefits of the proposed CKA-based meatric learning framework to a certain degree. From the obtained experimental results the following considerations are to be taken into account:

The use of a kernel learning method, termed Centered Kernel Alignment, allows quantifying the similarity between the input sample space and corresponding class label set. However, the core problem of implementing CKA is the way of choosing an appropriate kernel and setting its parameters, tending to guarantee good generalization performance of the learning machine. Besides its generalization ability, the Gaussian kernel is employed due to the possibility to incorporate a more generalized Euclidean metric that accounts for the variance of the feature and label spaces through a projection matrix. Therefore, the applied supervised metric learning projects the feature set so that the CKA value becomes maximal. However, the alignment optimization procedure is strongly dependent on the metric model. That is, simple linear mappings are easier to tune but cannot encode higher order relations. On the other hand, multi-layer architectures can discover data non-linearities, but are more expensive to tune.

The CKA dependence using the Mahalanobis metric learning is developed for three commonly used approaches to supervised classification: k-nn and SVM that are fed into a gradient-based optimization, and NN that are optimized using a backpropagation strategy with cross-entropy as the cost function. For the NN classifier, we generalize the supervised linear projection as a multilayer neural network, having one layer provided with a linear transfer function. Thus, a supervised CKA-based network pre-training is developed by applying the CKA procedure sequentially from the input to the output layer.

A bootstrapping procedure on the test set with 1,000 resamples estimates the average and the 95% confidence interval (CI) of each performance measure, as suggested by Bron et al. (2015). Besides, a two-sample t-test assesses the significative difference between a classifier and its CKA-enhanced version at a significance level of (1%). From the obtained results, it follows that the either case of the CKA metric learning gives rise to the classification performance for the multi-class scenario. Table 3 displays the validation performance of the multi-class scenario presented by each classifier. Due to the assumed influence of the topology, we also consider two versions of NN: 1-layered and 2-layered. In the first case of classification (top four rows), when feeding just the input feature set, the 2-layered NN classifier produces the best performance (average 57.4%), while the baseline k-NN certainly yields the worst results (45.5%). The bottom four rows hold the results of using the suggested learned metric, showing that the performance improves for most of classifiers. The last column computes the standard deviation of the true-positive rates (std_τ) a bias measure. As seen, the SVM is the most balanced classifier and benefits from the CKA by decreasing its standard deviation (std_τ = 4.2%). The gray-shaded cells denote the best performance of each classifier accuracy for all cases, noting that the pre-trained 1-layered NN is the best configuration. Another aspect to remark is that the CKA-based metric learning improves every one of the performance measures regardless of the used classifier as shown in the row Performance after CKA.

TABLE 3

Table 3. Testing performance of the metric learning for the multi-class task.

Likewise, the use of CKA also enhances the classification of each considered bi-class task as shown in Table 4. The distinction between HC and AD classes presents the highest values of performance (above 92% in average), meaning that this task is likely to have the easiest discriminating representation. However, the more challenging the bi-class task, the lower the increasing of classification performance.

TABLE 4

Table 4. Testing performance of the metric learning for the bi-class task.

A detailed analysis of each performance measure shows that the developed metric increases the accuracy obtained by all classification machines for each bi-class task. Also, the proposed metric increases the percentage of true positive fractions achieved by most of classifiers for each neurological class, improving the sensitivity of pathological groups and the specificity of healthy controls, especially for the HC vs. AD and HC vs. MCI task. Specifically, for the HC vs. MCI bi-class task, the proposed metric learning enhances significantly the performance balance of the k-nn and SVM classifier (2.3 and 8.3%, respectively) in comparison with their performance without ML (10.5 and 30.3%, respectively). Therefore, CKA forces each hidden space to increase the class separability by sequentially increasing the match between input and target data distributions. Nonetheless, most of the classifiers underperform for the MCI class due to its heterogeneous distribution. This situation arises because MCI is an intermediate disease between HC and AD, so that some patients may eventually convert to any of the other two classes.

For the sake of comparison, the proposed methodology of training is contrasted against the classification results for CADDementia challenge of AD diagnosis in Bron et al. (2015). To this end, the accuracy, true-positive rate, and area under the ROC (Receiver operating curve) are recomputed following the same evaluation scheme above described. Tables 4, 5 report the obtained results of CKA-enhanced SVM and NN machines along with the best performing approach in the challenge for each measure. It is worth noting that we can contrast results because class distributions and data acquisition protocols are similar for both works. However, we make clear that the comparison is given only in terms of the average and confidence intervals since the hypothesis test can not be performed.

TABLE 5

Table 5. Accuracy performance measures for the multi-class scenario following the validation scheme in Bron et al. (2015).

Results in Table 5 show that no approach achieves the best performance in all measures. In the case of the accuracy, Soerensen-equal outperforms with a = 63.0%, followed by the CKA-enhanced SVM and NN machines. Regarding the true-positive rates, Soerensen-equal achieves the highest τ_NC = 96.9 at the cost of decreasing τ_MCI and τ_AD (e.g., τ_AD = 61.2). Similarly, Abdulkadir and Routier-train algorithms bias toward AD (τ_AD = 80.6) and MCI (τ_MCI = 65.6), so that they perform the worst in MCI (τ_MCI = 21.3) and HC (τ_HC = 45.7), respectively. Therefore, the within approach true-positive rates, that vary largely, prove that some algorithms incline toward one class but reducing the discrimination of the other classes. In contrast, the proposed CKA enhancement provides more similar τ and a more balanced performance. Such a fact is proved by the standard deviation of the true-positives (std_τ), where the enhanced SVM and pre-trained 1-layered NN vary the least in comparison with results reported in the CADDementia challenge.

From Table 6 it follows that the per-class AUCs also depends on the selected approach. For example, Ledig-All and Wachinger-step1 algorithms provide the best AUC for HC and AD class, respectively. However, the former achieves 59.7% for MCI and the latter becomes the worst-performing. On the contrary, the proposed approach, particularly the pre-trained NN, achieve the best AUC for the MCI with no compromise of the other two classes. These results prove that the CKA criterion aims at discriminating all classes simultaneously obtaining balanced results and avoiding class biasing.

TABLE 6

Table 6. Area under the ROC performance measures for the multi-class scenario following the validation scheme in Bron et al. (2015).

5. Conclusion

This work introduces a supervised metric learning to support MRI classification. The proposed learning decodes discriminant information based on the maximization of the similarity between the input distribution and the corresponding target (diagnosis classes), aiming at enhancing the class separability. Thus, the conventional centered kernel alignment (CKA) is introduced as the cost function to infer the metric paramters that feed two classifiers (k-nn and SVM) and initialize the training of an NN. Evaluation of the proposed metric learning framework is carried out on the well-known ADNI dataset, where several morphological measurements are extracted using FreeSurfer to represent each MRI scan. In order to avoid over-fitting and prevent a biased classification performance, we split the dataset into two groups: A training subset for parameter tuning, and a test subset for performance assessment. Experiments are conducted according two diagnosis scenarios. For the multi-class one, our proposed CKA improves the performance for all classifiers in terms of the classification accuracy and the true positive fraction of each neurological class. In particular, the 1-layered NN classifier achieves the best performance (average 57.3%), and k-nn reaches competitive results (46.4%). In the bi-class scenario, CKA also enhances the classification, attaining measures over 90% in average for HC vs. AD task. Therefore, introducing CKA in the supervised metric learning determines spaces with a better class separability that increases both the sensitivity and specificity. As an additional benefit, the proposal better balances the performance in each class.

Attained results in the work yield the following research directions: Since the span of the resulting space is not guaranteed to be orthogonal, CKA criterion may produce redundant information. Then we propose to further introduce orthogonality constraints between the bases to solve this issue. Also, because the considered input feature set corresponds to a standard biomarker extraction, it is interesting to analyze other kinds of image representation strategies aiming at finding their relevance for class discrimination assessed by the CKA criterion. Besides, we devise two improvements of the sequential pre-training, namely, an NN training using CKA as the cost function and the evaluation of other multi-layer machines as Generative Stochastic Networks. Finally, we note that the class-wise performance can be parameterized by the introduced kernel function in the target space, so that a larger similarity of a particular class should increase its true positive rate.

Author Contributions

All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This research is supported by COLCIENCIAS projects: 111974455497 and 111956933522.

Footnotes

1. ^www.adni-info.org

2. ^freesurfer.nmr.mgh.harvard.edu

References

Álvarez-Meza, A. M., Cárdenas-Peña, D., and Castellanos-Dominguez, G. (2014). “Unsupervised kernel function building using maximization of information potential variability,” in Iberoamerican Congress on Pattern Recognition (Cham: Springer), 335–342.

Alzheimer's Association (2016). 2017 Alzheimer's Disease Facts and Figures. Technical Report.

Awate, S. P., Leahy, R. M., and Joshi, A. A. (2017). “Kernel methods for riemannian analysis of robust descriptors of the cerebral cortex,” in International Conference on Information Processing in Medical Imaging (Cham: Springer), 28–40.

Google Scholar

Bellet, A., Habrard, A., and Sebban, M. (2013). A survey on metric learning for feature vectors and structured data. arXiv preprint. arXiv:1306.6709.

Google Scholar

Bengio, Y. (2009). Learning deep architectures for AI. Found. Trends Mach. Learn. 2, 1–127. doi: 10.1561/2200000006

CrossRef Full Text | Google Scholar

Brockmeier, A. J., Choi, J. S., Kriminger, E. G., Francis, J. T., and Principe, J. C. (2014). Neural decoding with kernel-based metric learning. Neural Comput. 26, 1080–1107. doi: 10.1162/NECO_a_00591

PubMed Abstract | CrossRef Full Text | Google Scholar

Bron, E. E., Smits, M., Van Der Flier, W. M., Vrenken, H., Barkhof, F., Scheltens, P., et al. (2015). Standardized evaluation of algorithms for computer-aided diagnosis of dementia based on structural MRI: the caddementia challenge. NeuroImage 111, 562–579. doi: 10.1016/j.neuroimage.2015.01.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L., Head, D., Parker, J., Fotenos, A. F., Marcus, D., Morris, J. C., et al. (2004). A unified approach for morphometric and functional data analysis in young, old, and demented adults using automated atlas-based head size normalization: reliability and validation against manual measurement of total intracranial volume. NeuroImage 23, 724–738. doi: 10.1016/j.neuroimage.2004.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukumizu, K., Bach, F. R., and Jordan, M. I. (2004). Dimensionality reduction for supervised learning with reproducing kernel hilbert spaces. J. Mach. Learn. Res. 5, 73–99. Available online at: http://www.jmlr.org/papers/volume5/fukumizu04a/fukumizu04a.pdf

Google Scholar

Gretton, A., Bousquet, O., Smola, A., and Schölkopf, B. (2005). “Measuring statistical dependence with hilbert-schmidt norms,” in Algorithmic Learning Theory (Berlin; Heidelberg: Springer), 63–77.

Google Scholar

He, Y., Chen, W., Chen, Y., and Mao, Y. (2013). “Kernel density metric learning,” in Data Mining (ICDM), 2013 IEEE 13th International Conference on (Danvers, MA: IEEE), 271–280.

Google Scholar

Jack, C. R., Knopman, D. S., Jagust, W. J., Petersen, R. C., Weiner, M. W., Aisen, P. S., et al. (2013). Tracking pathophysiological processes in alzheimer's disease: an updated hypothetical model of dynamic biomarkers. Lancet Neurol. 12, 207–216. doi: 10.1016/S1474-4422(12)70291-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Khedher, L., Ramírez, J., Górriz, J. M., Brahim, A., and Illán, I. (2015). “Independent component analysis-based classification of alzheimers disease from segmented MRI data,” in International Work-Conference on the Interplay between Natural and Artificial Computation (Cham: Springer), 78–87.

Liu, F., Zhou, L., Shen, C., and Yin, J. (2014). Multiple kernel learning in the primal for multimodal alzheimers disease classification. IEEE J. Biomed. Health Inform. 18, 984–990. doi: 10.1109/JBHI.2013.2285378

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, B., Chen, Y., Hobbs, K., Smith, C. D., and Liu, J. (2015). “Nonlinear metric learning for alzheimer's disease diagnosis with integration of longitudinal neuroimaging features,” in BMVC, eds X. Xie, M. W. Jones, and G. K. L. Tam (BMVA Press), 138.1–138.13.

Google Scholar

Tustison, N. J., Cook, P. A., Klein, A., Song, G., Das, S. R., Duda, J. T., et al. (2014). Large-scale evaluation of ants and freesurfer cortical thickness measurements. Neuroimage 99, 166–179. doi: 10.1016/j.neuroimage.2014.05.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Wachinger, C., Reuter, M., Initiative, A. D. N., et al. (2016). Domain adaptation for alzheimer's disease diagnostics. Neuroimage 139, 470–479. doi: 10.1016/j.neuroimage.2016.05.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Kalousis, A., and Woznica, A. (2012). “Parametric local metric learning for nearest neighbor classification,” in Advances in Neural Information Processing Systems (La Jolla, CA: Curran Associates, Inc.), 1601–1609.

Google Scholar

Xu, Z., Weinberger, K. Q., and Chapelle, O. (2012). Distance metric learning for kernel machines. arXiv preprint. arXiv:1208.3422.

Google Scholar

Young, J., Mendelson, A., Cardoso, M. J., Modat, M., Ashburner, J., and Ourselin, S. (2015). “Improving mri brain image classification with anatomical regional kernels,” in Medical Learning Meets Medical Imaging (Cham: Springer), 45–53.

Google Scholar

Zhang, D., Wang, Y., Zhou, L., Yuan, H., Shen, D., Initiative, A. D. N., et al. (2011). Multimodal classification of alzheimer's disease and mild cognitive impairment. Neuroimage 55, 856–867. doi: 10.1016/j.neuroimage.2011.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar