Classification of Alzheimer’s Disease Stages: An Approach Using PCA-Based Algorithm

Abstract

Early diagnosis of Alzheimer’s disease (AD) allows individuals and their health managers to manage healthier medication. We proposed an approach for classification of AD stages, with respect to principal component analysis (PCA)-based algorithm. The PCA has been extensively applied as the most auspicious face-recognition algorithm. For the proposed algorithm, 100 images of 10 children were transformed for feature extraction and covariance matrix was constructed to obtain eigenvalues. The eigenvector provided a useful framework for face recognition. For the classification of AD stages, magnetic resonance imaging (MRI) and functional magnetic resonance imaging (fMRI) data were obtained from Alzheimer’s Disease Neuroimaging Initiative database. Hippocampus is one of the most affected regions by AD. Thus, we selected clusters of voxels from the “hippocampus” of AD screening stage (mild cognitive impairment), AD stage 1, AD stage 2, and AD stage 3. By using eigenvectors corresponding to maximum eigenvalues of fMRI data, the purposed algorithm classified the voxels of AD stages effectively.

Introduction

Faces are complex and multidimensional and have meaningful visual stimuli. Developing a system for face recognition is fairly problematic because of related shape of faces collective with several variations between images of the identical face. ¹ Nawaf et al developed a method for face recognition based on modular principal component analysis (PCA). ² Algorithms based on PCA are standard because of simplicity in implementation and have turned into the benchmark of algorithms. ³ Moon and Phillips established PCA-based algorithm comprises of normalization, PCA projection, and recognition modules. ³ The combination of optic flow and PCA can be used for the face recognition. ⁴ An algorithm was developed based on Hidden Markov model and optic flow for face recognition. ⁵ Cottrell and Metcalfe applied backpropagation neural network and PCA for recognition of facial expressions from still images. ⁶

Alzheimer’s disease (AD) is utmost numerous neurodegenerative disorder and emergent health problem. ⁷ Primary and accurate diagnosis/identification of AD is challenging but is decisive for the future treatments. ⁷ A method was purposed to classify AD (or mild cognitive impairment [MCI]) and healthy controls using kernel combination method by combining the magnetic resonance imaging (MRI), fluorodeoxyglucose positron emission tomography, and cerebrospinal fluid biomarkers. ⁸ The predetection of AD stages is perceived as significant because treatment might be more effectual if detected as timely as possible. ⁹ The MRI’s high-dimensional mode classification identifies patterns of brain structure that characterize MCI, which is usually the precursor to AD. ¹⁰ For the AD, a molecular test can lead for the better cure and rehabilitations, and 18 signaling proteins founded that could be used for the discrimination of AD and healthy controls blinded samples. ¹¹

Many MRI studies in AD stages have relied on measurements of volumes of particular brain regions, specifically the hippocampus and entorhinal cortex which shows histopathological changes at initial stages of AD. ⁹

Hippocampus is one of the primary regions of brain that is affected in AD, leading to loss of memory and confusion. ¹² The core function of hippocampus is primarily related with memory, and other major functions are spatial navigation, spatial memory, and behavioral inhibition. ¹² Impairment to the hippocampus can lead to loss of memory and difficulty in building new memories. This study used 2 data sets, first had 100 images of 10 children (10 × 10 = 100) with different facial expressions for face recognition while other had functional magnetic resonance imaging (fMRI) data of AD. The objective of this study was classification of faces (face recognition) and classification of AD stages with PCA-based algorithm. For classification of AD stages, brain’s region “hippocampus,” located in medial temporal, was analyzed.

Methodology

Principal component analysis is a well-known multivariate technique and well-recognized method for dimensionality reduction and feature extraction. ¹³ In PCA, we transform high-dimensional data into low-dimensional space. By reducing the degree of freedom and the complexity of space and time, PCA signifies data into a space that best expresses the variation. The PCA calculates the PCs that are base vectors in decreasing order of variability. ¹⁴

The method involves covariance matrix of input data, its eigenvectors and corresponding eigenvalues, and the projected data in new reductional space, defined by PCs. The first eigenvalue gives the direction of largest variance of data. Algebraic definition and process of PCs is given below.

Step 1: Given a sample of n observations on a vector of p variables. ¹⁴

$X=\left(X,X_2,X_3,\hspace{0.17em}\dots , X_p\right).$

1

Calculate the covariance matrix of data such as: ^{13 -15}

$\text{S}=\frac{1}{n-1}{\left(X-m\right)}^T\left(X-m\right),$

2

where X is the data matrix and m represents the mean vector. S is symmetric matrix.

Step 2: Define the first PC by the linear combination, ¹⁴

$z_1=a_1^{T}X$

3

where the vector $a_1=\left(a_{11}{a}_{21}{a}_{31},\dots ,a_{p1}\right)$ is chosen such that $var\left[z_1\right]$ is maximum. Likewise, define the kth PC of sample:

$z_k=a_k^T,k=1,\dots , p.$

Here, the vector $a_k=(a_{1k},a_{2k},a_{3k},\dots ,a_{pk}$ is chosen such that $var\left[z_k\right]$ is maximum. Subject to $cov\left[z_k,z_1\right]=0$ and $a_k^T{a}_1=1$ .

Step 3: To find $\hspace{0.17em}{a}_1$ , maximize $var\left[z_1\right]$ subject to $\hspace{0.17em}{a}_1^T{a}_1=1$ , where $a_1$ is an eigenvector of S corresponding to eigenvalue λ = λ₁. It is maximized as, ^{13 -15}

$var\left[z_1\right]=a_1^{T}S{a}_1=a_1^T{\text{λ}}_1{a}_1={\text{λ}}_1.$

4

So, λ₁ is the largest eigenvalue of S. The first PC, $z_1$ , retains the greatest amount of variation in the sample. Similarly, we find the second largest eigenvalue, λ₂ of covariance matrix S and second PC, $z_2$ , which contains the variation of sample less than first PC, $z_1$ . In general, ^{13 -15}

$var\left[z_k\right]=a_k^{T}S{a}_k=a_k^T{\text{λ}}_k{a}_k={\text{λ}}_k.$

5

The kth largest eigenvalue of S is the variance of kth PC. The kth PC $z_k$ retains the kth greater fraction of the variation in the sample. ^{13 -15}

Step 4: Arrange the eigenvectors $z_k$ in descending order with respect to eigenvalues in λ_k. ¹⁴

$\text{Projected data}={\left[z_k^T{\left(X-m\right)}^T\right]}^T,$

6

where each row of $z_k$ is an eigenvector. The PCA algorithm was applied to both data sets (face images and fMRI). The eigenvectors were selected on basis of largest eigenvalues of covariance matrix.

Analysis and Results

The PCA has been frequently used for problems related to faces recognition. ¹⁵ Ease and rapidity of PCA algorithm in face recognition is one advantage over other face-recognition algorithms. ¹⁵ In computation of PCA, there were 100 images of 10 children (10 images of each child). The face database is shown in Figure 1.

Figure 1.

Face database of children images.

Each grayscale image is converted into $n\times 1$ vector of pixels. By this, the data matrix for applying PCA consists 100 columns, each column is an image. Except one image for testing, calculate the mean image of faces, as shown in Figure 2.

Figure 2.

Mean face image of all images, except the test image.

Next, subtract the mean image from the data matrix and calculate the covariance/correlation matrix to apply PCA algorithm. As PCA used the eigenvalues and eigenvectors for recoginition, calculate the eigenvalues and eigenvectors/eigenfaces of covarinace matrix. Based on 10 largest eigenvalues, pick the eigenvector and calculate signature for each image. For the recognition of test image, subtract the mean image from test image and compare it with feature data matrix, by using PCA. The proposed algorithm will exactly match the test image, as shown in Figure 3.

Figure 3.

Test face (query face) and matched face from the database and number of PCAs used for best classification. PCA indicates principal component analysis.

We analyzed real fMRI data that were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (http://adni.loni.usc.edu). Investigators within ADNI contribute to the design and implementation of ADNI and provide the data but did not participate in analysis or write-up of this article. A complete list of ADNI investigators can be found at http://adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.

For preprocessing and analysis of fMRI AD data, MATLAB (version 2015, Statistics Toolbox: MathWorks, Massachusetts), Statistical Paramteric Mapping (SPM), version 2016 and MRIcron, version 2016, were used. The SPM involves 2 steps, preprocessing and model specification. During preprocessing of functional data, a low-pass filter option was applied to remove respiratory and cardiac noise effects, and generalized linear model was applied to analyze signal of subject. After preprocessing, to check the activation of fMRI data in SPM, 100 smoothed, realigned, and normalized images in model specification were used. Model specification involved some specifications or conditions about fMRI experiment to estimate the βs of fMRI brain data. In the next step, that is, results, because we used the resting-state fMRI data and there was no specific task-related condition in fMRI experiment, so for checking the activation of brain, we assumed dummy contrast, that is, [1 0]. Activations in brain regions during resting-state experiment are shown in Figure 4.

Figure 4.

Activation of brain’s region “hippocampus” in multislice view at MCI stage, AD stage 1, AD stage 2, and AD stage 3. AD indicates Alzheimer’s disease; MCI, mild cognitive impairment.

On the basis of activated voxels of the brain, region of interest, that is, hippocampus, was selected for the classsifaction of AD stages, as shown in Figure 5. While classifying the stages of AD, we selected the voxel clusters of hippocampus and used that voxels as vectors to apply the PCA algorithm. For classification of AD stages, we applied purposed alogorithm on that voxels. The objective of applying PCA was to develop an approach for classification of AD stages. The PCA-based algorithm classified the voxels of 4 stages of AD similar to classifcation of face images. This method involves the same steps as mentioned above for face recognition, that is, selecting a voxel for testing and calculation of covariance matrix of data, calcluation of eigenvalues and eigenvectors corresponding to largest eigenvalues of covarinace matrix, and the testing voxel matched with feature data matrix (which involves the voxels of 4 stages of AD). The PCA effectively matched the testing voxel with its actual stage of AD, as shown in Figure 6. Figure 7 showed the location of voxels recognized by proposed PCA-based algorithm, where red-, green-, and blue-colored voxels belongs to MCI, AD stage 1, and AD stage 2, respectively. The Scree plot of PCA shows the number of PCAs (approximately 10 PCAs) used for best classification, as shown in Figure 8.

Figure 5.

Red-colored area showing the region of interest (hippocampus) in coronal, axial, and sagittal view for classification of Alzheimer’s disease stages.

Figure 6.

Location of voxels in Montreal Neurological Institute coordinate system with detected stage of Alzheimer’s disease.

Figure 7.

Cutout brain showing the location of voxels recognized by proposed PCA-based algorithm where red-, green-, and blue-colored voxels belong to MCI, AD stage 1, and AD stage 2, respectively. AD indicates Alzheimer’s disease; MCI, mild cognitive impairment; PCA, principal component analysis.

Figure 8.

Scree plot showing that almost 10 PCAs used for best classification of Alzheimer’s disease stages. PCA indicates principal component analysis.

Discussion

The PCA-based algorithm has been widely used in problems related to face recognition. As PCA performs well in classification of faces, we made an approach by using PCA-based algorithm for classification of images and AD stages. Previously, we identified and classified the voxels of brain, used for rewardless-related decision-making using artificial neural network technique. ¹⁶ In our present approach, we used images as vectors of covariance matrix, and eigenvectors were obtained corresponding to largest eigenvalues for the classification and recognition of images. The AD has become common cause of dementia in all age groups. ¹⁷ The nature of AD is progressive, ¹⁸ so early diagnosis of AD is important and challenging for future treatments. ⁷ As we have mentioned earlier that hippocampus is one of the first regions of brain that is affected by AD, we analyzed the fMRI data using SPM and extracted the voxels of hippocampus region of 4 stages (MCI, AD stage 1, AD stage 2, and AD stage 3) of AD and obtained 95% classification rate. In our approach, we used that voxels as independent variables (vectors of covariance matrix). Eigenvectors based on largest eigenvalues were extracted and used for classification of AD stages. Our PCA-based algorithm for face-recognition approach classified the AD stages effectively, as shown Figure 7. This approach is simple and helpful for early diagnosis of AD. Early diagnosis of AD can help doctors healthier manage an individual’s conditions and sidestep medications that may deteriorate cognition or function. ¹⁸ Premature recognition of AD enables the best medical care and health outcomes for people with the disease, it also helps to doctors to identify and treat revocable conditions that mimic cognitive impairment and dementia, such as depression or insufficient vitamin. ^{18 -20}

Conclusion

In this article, an efficient human face-recognition technique based on PCA is proposed. On the basis of eigenvectors corresponding to largest eigenvalues, PCA provided a structure for face recognition. In addition, the technique presented here promises to distinguish the stages of AD. Rather than MRI, with the usage of fMRI data, we were able to identify and classify the stages of AD using activated voxels of hippocampus as input vectors of covariance matrix and through eigenvalues and eigenvectors which covers maximum variability of data, PCA provided an effective framework for the classification of AD. The PCA also provided best classification rate as compared to other methods. As hippocampus is one of the first regions of brain that is affected by AD, we purposed that, using fMRI data of a patient based on the voxels from hippocampus, PCA-based algorithm can identify and accurately classify the early diagnosis of AD in simple way using significant number of PCAs.