{"id":37465,"date":"2021-03-30T10:24:19","date_gmt":"2021-03-30T10:24:19","guid":{"rendered":"https:\/\/biomedpharmajournal.org\/?p=37465"},"modified":"2021-04-09T06:43:01","modified_gmt":"2021-04-09T06:43:01","slug":"estimation-of-number-of-levels-of-scaling-the-principal-components-in-denoising-eeg-signals","status":"publish","type":"post","link":"https:\/\/biomedpharmajournal.org\/staging\/vol14no1\/estimation-of-number-of-levels-of-scaling-the-principal-components-in-denoising-eeg-signals\/","title":{"rendered":"Estimation of Number of Levels of Scaling the Principal Components in Denoising EEG Signals"},"content":{"rendered":"<p style=\"text-align: justify;\"><strong>Introduction<\/strong><\/p>\n<p style=\"text-align: justify;\">The study of Electroencephalogram is very much helpful in diagnosing different disorders of the nervous system. EEG is the electrical action recorded from the scalp surface, which is picked up by conductive media and electrodes [1-3]. EEG has been performing a vital role in investigating brain activities in clinical application and scientific research for several years\u00a0[4-6]. The EEG signals can be contaminated by various artifacts, of which the major noise source is ocular artifact, which includes Eye-movement and eye-blink\u2019s [7]. However, artifacts are the major enemies of high-class EEG signals.\u00a0 The mixing up of these ocular artifacts with the EEG signal at the time of recording causes the problems in the accurate estimation of EEG signal. These artifacts will plunge into either of the 2 categories namely, technical and physiological artifacts. Power line noise 50\/60Hz falls into technical artifact category while the artifacts that crop up\u00a0 because of ocular(EOG), heart(ECG) and muscular activity(EMG)\u00a0 falls into\u00a0 physiological artifacts category respectively [8].<\/p>\n<p style=\"text-align: justify;\">Regression in the time domain and frequency domain [9-11] methods were proposed in removing eye blinks artifacts. These methods require a reliable reference channel. This channel can be contaminated by EEG. So, EEG has to be removed from the reference channel by regression techniques. Hence, the regression methods are not the finest to remove EOG artifacts.<\/p>\n<p style=\"text-align: justify;\">Principal Component Analysis is one of the available techniques for extracting the information from the data and has found applications in a wide range of disciplines [12]. PCA\u00a0 was introduced by Pearson in 1901[13] and developed by Hotelling [14]in the year 1933.<\/p>\n<p style=\"text-align: justify;\">In this research paper, up to how many levels the principal components have to be scaled for obtaining the better denoised EEG signal is elaborated using MSPCA and WAVELETS [15] and later estimated the SNR.<\/p>\n<p style=\"text-align: justify;\"><strong>Methodology<\/strong><\/p>\n<p style=\"text-align: justify;\"><strong>Principal Component Analysis<\/strong><\/p>\n<p style=\"text-align: justify;\">Principal Component Analysis (PCA) can be applied to EEG data that contains a large number of measured variables to develop into a smaller number of artificial variables called principal components(PC). Obtaining a smaller number of variables from a large number of measured variables is to reduce the redundancy in the measured variables. Here, redundancy means some of the variables in the measured data are correlated with one another,because they are measuring the same construct.The main idea of PCA is to reduce the dimensionality of the data set, as the data set consists of a large number of interrelated variables, and trying to retain as much as\u00a0 possible variation present in the data set.<\/p>\n<p style=\"text-align: justify;\">These principal components are uncorrelated,orthogonal\u00a0 and ordered in such a way that the first few components retain most variation present in all\u00a0 of the original variables. PCA is performed by eigen value decomposition of data covariance matrix.This is usually done after mean centering the data for each attribute.<\/p>\n<p style=\"text-align: justify;\">If the variables in a data set are already uncorrelated, PCA is of no value. In addition to being uncorrelated, the principal components are orthogonal and are ordered in terms of the variability they represent. That is, the first principal component represents, for a single dimension (i.e., variable), the greatest amount of variability in the original data set. Each succeeding orthogonal component accounts for as much of the remaining variability as possible.<\/p>\n<p style=\"text-align: justify;\">In other words, Principal component analysis (PCA) is a multivariate data analysis procedure that transforms a set of \u2018n\u2019 correlated variables, X = (x<sub>1<\/sub>, x<sub>2<\/sub>&#8230; x<sub>n<\/sub>,), into a set of uncorrelated variables called principal components (p<sub>1<\/sub>, p<sub>2<\/sub><em>, <\/em>&#8230;, p<sub>n<\/sub>). The first principal component accounts for most of the variability in the data, while each of the succeeding components in turn account for the highest amount of the remaining variability. Each principal component is a linear combination of the variables, X. The i<sup>th<\/sup> principal component can thus be expressed as:<\/p>\n<p style=\"text-align: justify;\">y<sub>i <\/sub>= e<sub>i<\/sub><sup>T<\/sup>X\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0\u2026\u2026\u2026.\u00a0 \u00a0 \u00a0 \u00a0 (1)<\/p>\n<p style=\"text-align: justify;\">where,<\/p>\n<p style=\"text-align: justify;\">e<sub>i<\/sub> is the eigenvector of the covariance matrix (R) of X (e<sub>i<\/sub><sup>T <\/sup>is the transpose of e<sub>i<\/sub>).<\/p>\n<p style=\"text-align: justify;\">The variance of the ith principal component is given by<\/p>\n<p style=\"text-align: justify;\">Var(Y<sub>i<\/sub>) =e<sub>i<\/sub><sup>T<\/sup> R ;\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0(2)<\/p>\n<p style=\"text-align: justify;\">where<\/p>\n<p style=\"text-align: justify;\">e =\u03bb<sub>i<\/sub>;\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 <em>i <\/em>= 1, 2. . . n<\/p>\n<p style=\"text-align: justify;\">\u03bb<sub>i\u00a0 <\/sub>is the i<sup>th\u00a0\u00a0 <\/sup>eigenvalue.<\/p>\n<p style=\"text-align: justify;\">PCA makes one stringent but powerful assumption, linearity. This assumption simplifies the problem by restricting the number of variables from the measured data. Hence PCA is used to re-express the data, which is a linear combination of original basis and which is explained here in terms of linear algebra.<\/p>\n<p style=\"text-align: justify;\">Consider an original data set X, which is an mxn matrix, where &#8216;m&#8217; corresponds to the number of measurement types and &#8216;n&#8217; is the number of samples. The main goal is to find an orthonormal matrix &#8216;P&#8217; in Y=PX such that C<sub>Y<\/sub>=1\/n(YY<sup>T\u00a0 <\/sup>) is a diagonal matrix. The matrix &#8216;P&#8217; transforms &#8216;X&#8217; into &#8216;Y&#8217;. The rows of orthonormal matrix &#8216;P&#8217; represent principal components of original data set &#8216;X&#8217;.<\/p>\n<p style=\"text-align: justify;\"><strong>First Level\u00a0<\/strong><strong>PCA<\/strong><\/p>\n<p style=\"text-align: justify;\">Two data sets, namely, EEG data set1(X1) and EEG data set2(X2), each of size 1&#215;1000, were collected from physionet.org website [16]. These two data sets were down sampled by a factor of 2. This reduces the size of each data set to 1&#215;500. Each data set is normalized using the following formula:<\/p>\n<p style=\"text-align: justify;\">X=(X-mean(X))\/std(X)<\/p>\n<p style=\"text-align: justify;\">where, mean(X) is mean of X<\/p>\n<p style=\"text-align: justify;\">std(X) is standard deviation of X<\/p>\n<p style=\"text-align: justify;\">The mean of each data set is calculated using the following formula:<\/p>\n<p style=\"text-align: justify;\"><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_eq1.jpg\"><img decoding=\"async\" class=\"alignnone size-full wp-image-37468\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_eq1.jpg\" alt=\"Vol14No1_Estt_Kri_eq1\" width=\"154\" height=\"107\" \/><\/a><\/p>\n<p style=\"text-align: justify;\">Where, n corresponds to number of samples in the data set X.<\/p>\n<p style=\"text-align: justify;\">The variance of each data set is calculated using the following formula:<\/p>\n<p style=\"text-align: justify;\"><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_eq2.jpg\"><img decoding=\"async\" class=\"alignnone size-full wp-image-37469\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_eq2.jpg\" alt=\"Vol14No1_Estt_Kri_eq2\" width=\"240\" height=\"101\" \/><\/a><\/p>\n<p style=\"text-align: justify;\">After obtaining the normalized data sets, a noise signal (EOG signal collected from physionet.org website), whose variance is of 0.4 and of length 500, is added to the two data sets. This results in two noisy data sets. These two noisy data sets are of size 1&#215;500 are converted into a column vector.<\/p>\n<p style=\"text-align: justify;\">X=[X1; X2].<\/p>\n<p style=\"text-align: justify;\">The size of the column vector will be equal to 500&#215;2.This column vector is treated as noisy EEG signal. For each column of data matrix wavelet decomposition is done to a level of 6 using sym8 wavelet. The wavelet decomposition gives approximate and detailed coefficients of noisy EEG signal. At each level of wavelet decomposition, i.e., on approximations and as well as details in wavelet domain, corresponding covariance matrices are computed. For each covariance matrix, PCA is performed. After performing the PCA on each covariance matrix, at each level, the most significant Principal Components (PC\u2019S) are selected. Here the selection of Principal Components is done using the Kaiser\u2019s rule. This rule retains the Principal Components that are associated with the Eigen values greater the mean of all Eigen values [17].<\/p>\n<p style=\"text-align: justify;\">The Principal Components at each level of decomposition and corresponding Principal Component variances vectors are provided in the Table 1<\/p>\n<p><strong>Table 1: Principal Components at Each Level of Decomposition and Corresponding Principal Component Variances Vectors-[FIRST SCALE PCA].<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td style=\"text-align: center;\" width=\"103\"><strong>Level No<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"220\"><strong>Principal Components Vector At Each Level\u00a0\u00a0 (Eigen Vectors)<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"166\"><strong>PC Variances<\/strong><\/p>\n<p><strong>(Eigen Values)<\/strong><\/td>\n<\/tr>\n<tr>\n<td width=\"103\">\n<p style=\"text-align: center;\">1<\/p>\n<\/td>\n<td style=\"text-align: center;\" width=\"220\">-0.71426\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0 0.69987<\/p>\n<p>0.69987 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a00.71426<\/td>\n<td width=\"166\">\n<p style=\"text-align: center;\">0.18818<\/p>\n<p style=\"text-align: center;\">0.14141<\/p>\n<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"103\">2<\/td>\n<td style=\"text-align: center;\" width=\"220\">0.31328\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0 0.94965<\/p>\n<p>0.94965\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 &#8211; 0.31328<\/td>\n<td style=\"text-align: center;\" width=\"166\">0.23541<\/p>\n<p>0.17906<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"103\">3<\/td>\n<td style=\"text-align: center;\" width=\"220\">-0.03374\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 0.99943<\/p>\n<p>-0.99943\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.03374<\/td>\n<td style=\"text-align: center;\" width=\"166\">0.79591<\/p>\n<p>0.55813<\/td>\n<\/tr>\n<tr>\n<td width=\"103\">\n<p style=\"text-align: center;\">4<\/p>\n<\/td>\n<td style=\"text-align: center;\" width=\"220\">0.33820\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.94107<\/p>\n<p>0.94107\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 0.33820<\/td>\n<td width=\"166\">\n<p style=\"text-align: center;\">4.19395<\/p>\n<p style=\"text-align: center;\">1.614760<\/p>\n<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"103\">5<\/td>\n<td style=\"text-align: center;\" width=\"220\">-0.46860\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 0.88341<\/p>\n<p>-0.88341\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.46860<\/td>\n<td style=\"text-align: center;\" width=\"166\">8.69436<\/p>\n<p>6.47539<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"103\">6<\/td>\n<td style=\"text-align: center;\" width=\"220\">0.04964\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.99877<\/p>\n<p>0.99877\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 0.04964<\/td>\n<td style=\"text-align: center;\" width=\"166\">11.31132<\/p>\n<p>2.59938<\/td>\n<\/tr>\n<tr>\n<td width=\"103\">\n<p style=\"text-align: center;\">7<\/p>\n<\/td>\n<td style=\"text-align: center;\" width=\"220\">0.93849\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.34529<\/p>\n<p>0.34529\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 0.93849<\/td>\n<td width=\"166\">\n<p style=\"text-align: center;\">136.82706<\/p>\n<p style=\"text-align: center;\">11.83056<\/p>\n<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"103\">8<\/td>\n<td style=\"text-align: center;\" width=\"220\">-0.32884\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a00<\/p>\n<p>-0.94439\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0 0<\/td>\n<td style=\"text-align: center;\" width=\"166\">1.12188<\/p>\n<p>0.21859<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\">Principal Components (1.12188 and 0.21859) corresponding to level 8 are the number of retained principal components for final PCA after wavelet reconstruction.<\/p>\n<p style=\"text-align: justify;\">From level 8, it is observed that the original data in two dimensional spaces is reduced to one dimension and shown below for ready reference.<\/p>\n<p style=\"text-align: justify;\"><strong>\u00a0\u00a0 <\/strong>LEVEL 8<\/p>\n<p style=\"text-align: justify;\"><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel1.jpg\"><img decoding=\"async\" class=\"alignnone size-full wp-image-37470\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel1.jpg\" alt=\"Vol14No1_Estt_Kri_lel1\" width=\"417\" height=\"152\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel1-300x109.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel1.jpg 417w\" sizes=\"(max-width: 417px) 100vw, 417px\" \/><\/a><\/p>\n<p style=\"text-align: justify;\">Using these principal components one can reconstruct the denoised version of the input matrix X. The denoised versions of the input matrix X i.e., EEG data set1 and EEG data set2 are shown in the Fig. 1 (a) and Fig. 1(b) respectively.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1A.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-37477\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1A-150x150.jpg\" alt=\"Vol14No1_Estt_Kri_fig1(A)\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1A-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1A-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1A.jpg 867w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 1(a): Original and Artifact Removed EEG Signal of Data\u00a0<\/strong><strong>Set 1 using First Scale PCA.<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1A.jpg\" target=\"_blank\">Click here to view figure\u00a0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1B.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-37478\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1B-150x150.jpg\" alt=\"Vol14No1_Estt_Kri_fig1(B)\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1B-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1B-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1B.jpg 841w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 1(b): Original and Artifact Removed EEG Signal of Data\u00a0<\/strong><strong>Set 2 using First Scale PCA.<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig1B.jpg\" target=\"_blank\">Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\">The quality of column reconstructions is estimated by the relative mean square error are 22.1904933458124% and 73.9497428008459%, not closer to 100%.<\/p>\n<p style=\"text-align: justify;\">Since the quality of column reconstructions after first level PCA are not closer to 100%, hence the level of scaling principal components is taken to next level [18]. Hence retaining of principal components will be decided based on the quality of reconstruction of columns which is measured by relative mean square error.<\/p>\n<p style=\"text-align: justify;\"><strong>Second Level PCA<\/strong><\/p>\n<p style=\"text-align: justify;\">The simplified input matrix X, which was obtained from the first scale of PCA, is again decomposed to a level of 6 using sym8 wavelet. Now the wavelet coefficients obtained after the wavelet decomposition are thresholded using Heursure thresholding. For these wavelet coefficients PCA is performed and selected the significant principal components. Using these principal components one can reconstruct the much more denoised input matrix X. The quality of reconstruction of the columns estimated after the second time processing of the input matrix x are close to 100% and are 99.9981% and 99.9991%.<\/p>\n<p style=\"text-align: justify;\">The principal components and Principal Component variances vectors of the two data sets obtained after second scale PCA and wavelet denoising are shown in the Table.2.<\/p>\n<p style=\"text-align: justify;\"><strong>Table 2: Principal Components and Principal Component Variances Vectors Obtained After Second Scale PCA and Wavelet Denoising.<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td style=\"text-align: center;\" width=\"47\"><strong>Level No<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"217\"><strong>Principal Components Vector At Each Level (Eigen Vector)<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"227\"><strong>PC Variances (Eigen value)<\/strong><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">1<\/td>\n<td style=\"text-align: center;\" width=\"217\">-0.74792 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.66378<\/p>\n<p>-0.66378\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 0.74792<\/td>\n<td style=\"text-align: center;\" width=\"227\">0.17755<\/p>\n<p>0.14424<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">2<\/td>\n<td style=\"text-align: center;\" width=\"217\">-0.47075\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.88226<\/p>\n<p>0.88226\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0 -0.47075<\/td>\n<td style=\"text-align: center;\" width=\"227\">0.18846<\/p>\n<p>0.17321<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">3<\/td>\n<td style=\"text-align: center;\" width=\"217\">-0.03589 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.99936<\/p>\n<p>0.99936\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.03589<\/td>\n<td style=\"text-align: center;\" width=\"227\">0.73933<\/p>\n<p>0.43569<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">4<\/td>\n<td style=\"text-align: center;\" width=\"217\">0.45131\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0 0.89237<\/p>\n<p>0.89237\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.45131<\/td>\n<td style=\"text-align: center;\" width=\"227\">4.01686<\/p>\n<p>1.39216<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">5<\/td>\n<td style=\"text-align: center;\" width=\"217\">0.78691\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0 0.61707<\/p>\n<p>0.61707\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0 -0.78691<\/td>\n<td style=\"text-align: center;\" width=\"227\">9.15681<\/p>\n<p>5.49690<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">6<\/td>\n<td style=\"text-align: center;\" width=\"217\">0.04293\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0 -0.99907<\/p>\n<p>0.99907\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a00.04293<\/td>\n<td style=\"text-align: center;\" width=\"227\">11.28581<\/p>\n<p>2.15695<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">7<\/td>\n<td style=\"text-align: center;\" width=\"217\">0.93978\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0 0.34178<\/p>\n<p>0.34178\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 -0.93978<\/td>\n<td style=\"text-align: center;\" width=\"227\">175.29829<\/p>\n<p>6.19914<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"47\">8<\/td>\n<td style=\"text-align: center;\" width=\"217\">0.58069\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a00<\/p>\n<p>0.81412\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a00<\/td>\n<td style=\"text-align: center;\" width=\"227\">1.27604<\/p>\n<p>0.28103<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\">Principal Components (1.27604 and 0.28103) corresponding to level 8 are the number of retained principal components for final PCA after wavelet reconstruction.<\/p>\n<p style=\"text-align: justify;\">From level 8, it is observed that the original data in two dimensional spaces is reduced to one dimension and shown below for ready reference.<\/p>\n<p style=\"text-align: justify;\">Level 8<\/p>\n<p style=\"text-align: justify;\"><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel2.jpg\"><img decoding=\"async\" class=\"alignnone size-full wp-image-37471\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel2.jpg\" alt=\"Vol14No1_Estt_Kri_lel2\" width=\"378\" height=\"145\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel2-300x115.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_lel2.jpg 378w\" sizes=\"(max-width: 378px) 100vw, 378px\" \/><\/a><\/p>\n<p style=\"text-align: justify;\">On comparing Table.1 and Table2\u00a0 the final PC\u2019s in the level 8 are improved after second time PCA using wavelet denoising [19].Using the selected principal components one can reconstruct noise free EEG signals back. The waveforms of Denoised EEG signal constructed after second scale PCA are shown in Fig 2(a) and 2(b).<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2A.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-37479\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2A-150x150.jpg\" alt=\"Vol14No1_Estt_Kri_fig2(A)\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2A-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2A-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2A.jpg 740w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 2(a): Original and Artifact Removed EEG Signal of data\u00a0<\/strong><strong>set 1 with MSPCA<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2A.jpg\" target=\"_blank\">Click here to view figure\u00a0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2B.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-37480\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2B-150x150.jpg\" alt=\"Vol14No1_Estt_Kri_fig2(B)\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2B-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2B-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2B.jpg 727w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 2(b): Original and Artifact Removed EEG Signal of data\u00a0<\/strong><strong>Set 2 with MSPCA.<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/02\/Vol14No1_Estt_Kri_fig2B.jpg\" target=\"_blank\">Click here to view figure\u00a0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>Results<\/strong><\/p>\n<p style=\"text-align: justify;\">The results obtained after performing MSPCA on different data sets were tabulated in Table 3 and also compared with the previous results [18].<\/p>\n<p style=\"text-align: justify;\"><strong>Table 3: Comparison of SNR of Denoised EEG Signal Obtained using MSPCA and Different Wavelets.<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td width=\"174\"><\/td>\n<td style=\"text-align: center;\" width=\"277\"><strong>Wavelet family\/Threshold type<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"152\"><strong>SNR (dB)<\/strong><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" rowspan=\"6\" width=\"174\">1)B Krishna Kumar<\/p>\n<p>Analysis of EEG Signals Using Multi-Scale<\/p>\n<p>Principal Component<\/p>\n<p>Analysis-<\/p>\n<p>2) B. K. Kumar and K. V. S. V. R. Prasad, &#8220;Performance comparison of IST and multi scale principal component analysis in the EEG signal processing,&#8221;\u00a0<em>2017 International Conference on Computing Methodologies and Communication (ICCMC)<\/em><\/td>\n<td style=\"text-align: center;\" width=\"277\">Sym8 (soft)<\/td>\n<td style=\"text-align: center;\" width=\"152\">34.95<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"277\">Sym8 (hard)<\/td>\n<td style=\"text-align: center;\" width=\"152\">24.52<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"277\">Haar (soft)<\/td>\n<td style=\"text-align: center;\" width=\"152\">20.29<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"277\">Haar (hard)<\/td>\n<td style=\"text-align: center;\" width=\"152\">25.10<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"277\">dB10 (soft)<\/td>\n<td style=\"text-align: center;\" width=\"152\">28.89<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"277\">dB10 (hard)<\/td>\n<td style=\"text-align: center;\" width=\"152\">23.74<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"174\">Proposed method<\/td>\n<td style=\"text-align: center;\" width=\"277\">MSPCA ( for DATA SETS-1 AND 2)<\/td>\n<td style=\"text-align: center;\" width=\"152\">37.43<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"174\">Proposed method<\/td>\n<td style=\"text-align: center;\" width=\"277\">MSPCA( for DATA SETS-2 AND 3)<\/td>\n<td style=\"text-align: center;\" width=\"152\">35.36<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"174\">Proposed method<\/td>\n<td style=\"text-align: center;\" width=\"277\">MSPCA( for DATA SETS -4 AND 5 )<\/td>\n<td style=\"text-align: center;\" width=\"152\">35.16<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>Conclusion<\/strong><\/p>\n<p style=\"text-align: justify;\">The MSPCA is providing better SNR as the Relative Mean Square Error(RMSE) of the columns is closure to 100%. Hence, it is important to check the Relative Mean Square Error (RMSE) of the columns before reconstructing the Denoised EEG signal and estimating the SNR.<\/p>\n<p style=\"text-align: justify;\"><strong>Acknowledgement\u00a0\u00a0 <\/strong><\/p>\n<p style=\"text-align: justify;\">Not applicable<\/p>\n<p style=\"text-align: justify;\"><strong>Conflict of Interest\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 <\/strong><\/p>\n<p style=\"background: white; margin: 0in 0in 9.6pt; text-align: justify;\"><span lang=\"EN-IN\" style=\"font-size: 13.5pt;\">There is no conflict of interest.<\/span><\/p>\n<p style=\"text-align: justify;\"><strong>Funding Source\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 <\/strong><\/p>\n<p style=\"text-align: justify;\">Not applicable<\/p>\n<p style=\"text-align: justify;\"><strong>References<\/strong><\/p>\n<ol>\n<li style=\"text-align: justify;\">Krishna Kumar \u201c Denoising of EEG Signals using Wavelets and SIMULINK Techniques \u201cInternational Journal of Recent Technology and Engineering (IJRTE) , Volume-8 Issue-5, January 2020 pp.335-339<br \/>\n<a href=\"https:\/\/doi.org\/10.35940\/ijrte.C5113.018520\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">E.Niedermeyer and FH Silva. \u201cElectroencephalography: Basic principles, clinical applications and related fields\u201d, Lippincott, Williams &amp; Wilkins, (2004),<\/li>\n<li style=\"text-align: justify;\">MR Arab, AA Suratgar, VMM Hernandez, AR Ashtiani \u201cElectroencephalogram Signals Processing for the Diagnosis of Petit mal and Grand mal Epilepsies Using an Artficial Neural Network.Journal of Applied Research and Technology, 8 (2010), pp. 120-129.<br \/>\n<a href=\"https:\/\/doi.org\/10.22201\/icat.16656423.2010.8.01.483\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">GL Holmes, ad CT Lombroso. \u201cPrognostic value of background patterns in the neonatal EEG\u201d.J. Clin. Neurophysiol, 10 (1993), pp. 323-352<br \/>\n<a href=\"https:\/\/doi.org\/10.1097\/00004691-199307000-00008\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">S Almubarak, and PK Wong.\u201d Long-Term Clinical Outcome of Neonatal EEG Findings\u2019\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 J.Clin.Neurophysiol,28(2011),pp.185- 189http:\/\/dx.doi.org\/10.1097\/WNP.0b013e3182121731\u00a0|\u00a0Medline<br \/>\n<a href=\"https:\/\/doi.org\/10.1097\/WNP.0b013e3182121731\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">ASM Muthanantha Murugavel, S Ramakrishnan.\u201dTree Based Wavelet Transform and DAG SVM for Seizure Detection.Signal and Image Processing\u201d: An International Journal, 3 (2012), pp. 115-125<br \/>\n<a href=\"https:\/\/doi.org\/10.5121\/sipij.2012.3111\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">N.V. Thakor et al. (1993). \u201cMulti resolution Wavelet Analysis of Evoked Potentials\u201d, IEEE Transactions on Biomedical Engineering, Vol. 40, No 11, pp. 1085-1093,November.<br \/>\n<a href=\"https:\/\/doi.org\/10.1109\/10.245625\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">S. Ventakaramanan, P. Prabhat, S.R Choudhury, H.B Nemade, and J.S. Sahambi. (2000). \u201cBiomedical Instrumentation Based On Electrooculogram (EOG) Signal Processing And Application To A Hospital Alarm System\u201d, Indian Institute Of Technology (IIT) Gauhati, Proceedings of IEEE ICISEP, pp.535-539.<\/li>\n<li style=\"text-align: justify;\">Schlogl A, Keinrath C, Zimmermann D et al. 2007. A fully automated correction method of EOG artifacts in EEG recordings. Clinical Neurophysiology.118, 98\u2013104.<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/j.clinph.2006.09.003\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Jung T-P, Makeig S, Humphries C, Lee T-W, McKeown MJ,Iragui V and Sejnowski TJ. 2000a. Removing Electroencephalographic artifacts by blind source separation.Psychophysiology.37, 163 \u2013178.<br \/>\n<a href=\"https:\/\/doi.org\/10.1111\/1469-8986.3720163\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Gratton G, Coles MG and Donchin E. 1983. A new method for off-line removal of oucular artifact. Electroencephalography Clin.Neurophysiol. 55, 484\u2013 486.<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/0013-4694(83)90135-9\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">M.Aminghafari,N. cheze, J.M Poggi, \u201dMultivariate denoising using wavelets and principal component Analysis,\u201d computational statistics &amp; Data Analysis<em>, <\/em>vol. 50 pp. 2381-2398,2006.<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/j.csda.2004.12.010\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Pearson, K, &#8220;On Lines and Planes of Closest Fit to Systems of Points in Space&#8221;, Phil. Magazine, vol.6, issue. 2, pp.559\u2013572, 1901.<br \/>\n<a href=\"https:\/\/doi.org\/10.1080\/14786440109462720\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Hotelling. H,&#8221; Analysis of a Complex of Statistical Variables Into Principal Components&#8221;, J. Educ. Psychol., 24, pp.417\u2013441, pp.498\u2013520, 1933.<br \/>\n<a href=\"https:\/\/doi.org\/10.1037\/h0070888\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">B. K. Kumar and K. V. S. V. R. Prasad, &#8220;Performance comparison of IST and multi scale principal component analysis in the EEG signal processing,&#8221;\u00a0<em>2017 International Conference on Computing Methodologies and Communication (ICCMC)<\/em>, Erode, 2017, pp. 536-542, doi: 10.1109\/ICCMC.2017.8282523.<br \/>\n<a href=\"https:\/\/doi.org\/10.1109\/ICCMC.2017.8282523\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">www.physionet.org<\/li>\n<li style=\"text-align: justify;\">www.mathworks.com\/products\/matlab\/, Discrete Wavelet Transform-2010.<\/li>\n<li style=\"text-align: justify;\">B.KrishnaKumar ,Dr.K.V.S.V.R.Prasad, Dr.K.Kishan Rao\u00a0 and\u00a0 Narsimha Baddiri, &#8220;Analysis of EEG Signals Using Multi-Scale Principal Component Analysis&#8221;, IEEE Conference, Cape Institute of Technology, Levingapuram, Kanyakumari, pp.359-363, December15-18, 2011.<\/li>\n<li style=\"text-align: justify;\">W.M.Bukhari W, Daud and Rubita Sudirman, &#8220;Time Frequency Analysis of Electrooculography (EOG) Signal of Eye Movement Potentials Based on Wavelet Energy Distribution&#8221;,2011 Fifth Asia Modeling Symposium, pp.81-86,2011.<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Introduction The study of Electroencephalogram is very much helpful in  [&#8230;]<\/p>\n","protected":false},"author":14,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[88],"tags":[],"class_list":["post-37465","post","type-post","status-publish","format-standard","hentry","category-vol14no1"],"_links":{"self":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/37465","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/users\/14"}],"replies":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/comments?post=37465"}],"version-history":[{"count":5,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/37465\/revisions"}],"predecessor-version":[{"id":38494,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/37465\/revisions\/38494"}],"wp:attachment":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/media?parent=37465"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/categories?post=37465"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/tags?post=37465"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}