International Journal of Mathematical Research

March 2018, Volume 7, 1, pp 1-17

# Multivariate Analysis of EEG Data: Some Aspects of Diagnostic of MANOVA Model

Md Rokonuzzaman

Md Rokonuzzaman 1

1. Associate Professor Department of Statistics University of Chittagong Chittagong-4331 Bangladesh 1

on PubMed

Pages: 1-17

DOI: 10.18488/journal.24.2018.71.1.17

#### Share :

Article History:

Revised: 21 August, 2018
Accepted: 24 August, 2018
Published: 03 September, 2018

# Index

## Abstract

The main focus of this study is the multivariate analysis of Electroencephalogram data which included multivariate analysis of variance. The multivariate model diagnostics comprise checking number of assumptions of MANOVA model such as multivariate normality, homogeneity of covariance matrices. In this paper, the model X=BC+E is used and estimate the different parameters. Also by using a general form, H0: GBF=0 to test the different types of null hypothesis. Here G and F is known matrices and obtained from the hypothesis. This study gives mathematical ideas from multivariate statistical analysis to find a solution or a good approximation of a complex scientific problem.

Keywords:Electroencephalogram, Fractal dimension, MANOVA.

Received: 7 September 2017 / Revised: 21 August 2018 / Accepted: 24 August 2018 / Published: 3 September 2018

#### Contribution/ Originality

This study gives mathematical ideas from multivariate statistical analysis to find a solution or a good approximation of a complex scientific problem.

## 1. INTRODUCTION

Electroconvulsive therapy (ECT) formerly known as electroshock, is a psychiatric treatment used to treat severe depressed patients. People with symptoms such as delusions, hallucinations, or suicidal thoughts or when other treatments such as psychotherapy and antidepressant medicines have not worked, are exposed to ECT. The composite activity or measure of brain waves is called an electroencephalogram (EEG). Nowadays, there is no common statistical approach to analyze the EEG data. Some analysis of EEG signals using dynamic models have been conducted in neurophysiology (e.g. [1-3]). Recent research efforts have been based on various numerical algorithms that estimate individual dynamics and statistical information characteristics have been considered in more detail [4, 5]. The dynamical approach was subsequently applied to a variety of physiological time series, including EEG recordings [6-10]. Various mathematical techniques developed for the studies of complex, nonlinear and chaotic systems have also been used to explore the dynamics of human brain physiology [8, 11, 12].

The fractal dimension (FD) is used to quantify the complexity and self-similarity of a signal [13]. FD analysis of EEG data has been applied in the biomedical research, such as the routine detection of dementia [14, 15] EEG analysis of sleeping newborns [16]  and fractal spectral analysis of pre-epileptic seizures [17]. The use of FD in the analysis of EEG data was suggested by Popivanov, et al. [18]. Wahlund, et al. [19] proposed an approach where one looked at FD of the curve as the representation of EEG data obtained from ECT seizures to identify differences in responses between parts of the brain. In some other studies, estimation of FD of similar biological signal has been used. For the data used from magneto encephalogram (MEG), Gómez, et al. [20] estimated the complexity of those signals in Alzhemer’s patients and healthy elderly. They found significantly lower complexity of MEG in Alzheimer’s patients compared to the control subjects.

Classical approaches where measurement errors are regarded as generators of uncertainty are not valid when considering EEG time series data [21]. Due to strong individual responses to ECT one possible approach is to use mixed models where individual effects are taken into account via random effects. These individual effects can then be discarded or incorporated in the analysis. The multivariate analysis of variance (MANOVA) is a type of multivariate statistical technique to analyze data that involves more than one response at a time. One of the aims of this work is to apply multivariate statistical techniques for the analysis of EEG signals recorded from depressed patients. Complexity of EEG signal is assessed by using FD according to Wahlund, et al. [19].

In the next section, background characteristics of EEG data including the assessing complexity of EEG signal will be discussed. In Section 3, multivariate approach to the statistical analysis of EEG data will be presented. Issue of model checking as assessment of the multivariate normality, and equality of covariance matrices also discussed here.  In Section 4, the obtained results will be presented.  Finally, the last section comprises the discussion.

## 2. BACKGROUND

ECT is used for people whose depression is severe or life threatening or who cannot take antidepressant medication. ECT is then often recommended for use as a treatment for severe depression that has not responded to other treatment, and is also used in the treatment of mania and catatonia. Before ECT, a patient is given anesthesia, and medicines to loosen up his muscles. Then electrical current is briefly sent to the brain through electrodes placed on the temples or elsewhere on the head, depending on the condition and type of ECT. Usually ECT treatments are done 2 to 3 times a week for 2 to 3 weeks and continuing for several months to a year by reducing time frequency, to reduce the risk of relapse. ECT is usually given in combination with medicine, psychotherapy, family therapy, and behavioral therapy.

EEG is the recording of electrical activity along the scalp and measures the voltage fluctuating result from ionic current flows within the neurons of the brain.  In clinical contexts, EEG refers to the recording of the brain's spontaneous electrical activity over a short period of time, usually 20–40 minutes, as recorded from multiple electrodes placed on the scalp. The EEG is used in the evaluation of brain disorders. Most commonly it is used to show the type and location of the activity in the brain during a seizure. The EEG is typically described in terms of rhythmic activity and transients. The rhythmic activity is divided into bands by frequency. To some degree, these frequency bands are a matter of nomenclature, but these designations arose because rhythmic activity within a certain frequency range was noted to have a certain distribution over the scalp or a certain biological significance. Most of the cerebral signal observed in the scalp EEG falls in the range of 1–20 Hz. The activity below or above this range is likely to be artifactual, under standard clinical recording techniques. However, brain wave EEG has a complex pattern of frequencies. Experimentally, it has been established that there are several frequently bands that’s associated with particular processes in the brain : delta waves (<4 Hz), thetawaves (4-8 Hz), alpha waves (8-13 Hz), beta waves (13-36 Hz) and gamma waves (> 36 Hz).

In the next section, analysis and assessing the complexity of EEG signals will be discussed.

2.1. Analysis and Assessing Complexity of EEG Signals

Multivariate analysis of variance is one of the techniques for analyzing multivariate data. Now-a-days, multivariate statistical analysis is a very popular technique applied in such research fields as education, biology, medicine and social sciences.  If there have categorical independent variables with two or more levels and a set of dependent variables and one can interested in determining if the categorical independent variables affect the dependent variables. MANOVA can be used to address of the preceding problem. One reason for conducting a MANOVA is to determinate how variates can be combined to maximally discriminating between district subgroups of the responses. Many research articles and books have been devoted to MANOVA in the last twenty five years. Weinfurt [22] and Huberty and Petoskey [23] introduced complex and conceptually powerful statistical procedure for MANOVA in their books.

Amplitude, frequency, and degree of synchronization constitute three basic parameters of background EEG activity. The first two parameters are easy to quantify, but recently synchronization only has become a subject of computerized EEG analysis. The opposite of synchronization is complexity of EEG signal. Signal complexity [24] can be analyzed either in frequency domain, or in time domain, or in the system's phase space. Analysis in the frequency domain   requires   fourier or wavelet transform of the signal, while analysis in the phase space requires embedding of the data in a multidimensional space.  Complexity in time domain may be characterized by fractal dimension, Df. It characterizes complexity of the curve representing the signal on a plane. The fractal dimension Df has values between 1 and 2. Df calculated in time domain should not be confused with fractal dimension of the attractor in the system's phase space. An object in Euclidean space, like a coastline of an island, may also be characterized by a fractal dimension. In a rough sense, it measures 'how many points' belong to a given set. A plane is 'larger' than a line, while Sierpinski triangle sits somewhere in between these two sets. On the other hand, all three of these sets have the same number of points in the sense that each set is uncountable. Somehow, though, fractal dimension captures the notion of 'how large a set is' while its fractional part captures the notion of 'how complex a set is'.  In 1988, Higuchi proposed an efficient algorithm for measuring the fractal dimension (FD) of discrete time sequences [25]. Higuchi’s algorithm calculates the FD directly from time series. Higuchi’s fractal dimension (HFD) has already been used to analyze the complexity of brain recording [13, 26] and other some biological signals [27].

consideration in a plane as a function of time. If the length of this curve is scales as

correspond to the presence of higher frequencies in the signal’s Fourier spectrum.

Let us define x(1), x(2), ..., x(N) be the values of a finite set of time series observations, which are taken in a regular interval, where i= 1,2,…,N  and N is the number of points in the time series. Here x are the successive EEG values.

In the next section, we will discuss about the multivariate tests that are common used in MANOVA. Also multivariate tests based on contrasts will discuss and multivariate modeling of fractal dimension with estimation of parameters.

## 3. MULTIVARIATE ANALYSES

Let us consider the following MANOVA model:

X= BC + E,

where X: p×n, B: p×k, C: k×n and E~Np,n(0,∑, I). MatricesB and are unknown parameters. X is the data matrix; C is known between individual design matrix and the columns of E are independent and the dispersion of each column equal to S.

To test the following null hypothesis about equality of the mean responses with the assumption; the covariance matrices of the different groups of responses are assumed to be equal, Wilks’ lamda, Lawley-Hotelling’s trace, Pillai’s Trace and Roy’s largest root tests [28] can be used to test the equality of the mean responses. These tests, however, may become seriously biased when the assumption of equality of covariance matrices is violated.

Now we will discuss test statistics for contrasts.

3.1. Multivariate Significance Test for the Contrasts

A contrast is a linear combination of two or more group means. We have a contrast among the population

3.2. Multivariate Modeling of Fractal Dimension

In the model X= BC + E, our interest is to test the hypothesis GB=0, where G is known and formulated according to null hypothesis . However, we shall extend the model (3) and consider

X= BC + E,            GBF=0,                                             (5)

where F is known. According to Kollo and Rosen [29] this model is equivalent to an extended Growth Curve model, i.e.

The test statistic is

where

The following asymptotic expansion gives a good approximation for any constant z>0,

represents terms that goes to zero.

and the parameter B in (4)  is defined as,

The following hypotheses are of interest:

1.  The  second hypothesis of interest is

H02: There is no difference between the two phases (Phase I-III and Phase IV),

In next section, we will focus on model validation.

3.3. Checking the Equality of Variance-Covariance Matrices

Another important assumption of MANOVA is that the covariance matrices across all groups are equal. It implies that there exists no substantive difference between the variances of different groups of dependent variables. In a simple regression model, it is defined as the homogeneity of residual variances.

Box’s M test statistic is used to test the null hypothesis about the homogeneity of covariance matrices.  In the next section, about the data are discussed.

## 4. DATA

The data used in this example are synthetic data, i.e. responses from patients are created so that they resemble responses from real study patients. The sample consists of 28 patients with severally depressed patients and with three different diagnoses. Five patients have psychotic depression (PD), seven have bipolar disorder (BD) and sixteen patients have unipolar depression (UD). Two channels of EEG of left and right prefrontal-to-ipsilateral coded Fp1 and Fp2, respectively, were recorded. The Fp1 and Fp2 were assessed according to the international 10-20 system (see Figure1), an internationally recognized method to describe and apply the location of scalp electrodes in the context of an EEG test. This system is based on the relationship between the location of an electrode and the underlying area of cerebral cortex. The "10" and "20" refer to the actual distances between two electrodes are either 10 or 20 percent of the total front–back or right–left distance of the scalp.

The digitalized data from the digital port of MECTA equipment were simultaneously recorded and transferred to computer hard disc drive. The continuous EEG signals were digitalized at a sampling frequency of 128 Hz. None of the patients had received ECT during the last three months. Every patient was treated 6-12 times with recording of EEG and the average of the readings was used in the analysis.

Figure-1. Map of locations of the EEG electrodes/channels

Source: Map of locations of the EEG electrodes/channels on the scalp according to International 10-20 System. EEGs were recorded from Fp1 and Fp2, which are the standard locations for treating patients with ECT.

Now we will discuss the obtained results.

## 5. STATISTICAL ANALYSIS

5.1. Descriptive Statistics

Table 5.1 contains the mean values of Df of two phases Phase I-III and Phase IV for the two different electrodes/channels from different diagnoses groups of depressed patients. The result shows that the mean values of Df for phase I-III are lower than those for phase IV. Lower SE values implies lower variability of Df  are observed for the both channels Fp1 and Fp2 for Phase I-III.

Table-5.1. Average of fractal dimension Df (SE) value for different diagnoses groups of depressed patients; from Phase I-III and Phase IV at different electrode locations.

 Diagnosis N Fp1 Fp2 Phase I-III Phase-IV Phase I-III Phase-IV Psychotic depression 5 1.210 (0.018) 1.538 (0.049) 1.194 (0.022) 1.712 (0.064) Bipolar disorder 7 1.266 (0.025) 1.513 (0.053) 1.261 (0.023) 1.553 (0.069) Unipolar depression 16 1.250 (0.015) 1.584 (0.036) 1.251 (0.015) 1.630 (0.044)

5.2. Multivariate Analysis of Variance

Multivariate tests discussed in Section 3.1, are conducted to determine if the mean vectors for three different groups (PD, BD and UD) of the four dependent variables are significantly different or not (an overall or global test for equality of mean vectors). The results presented in Table 5.2 (significant that indicates at the 5% level of significance) indicate that mean vectors for three groups are statistically different.

Table-5.2. Different test statistics for the null hypothesis; H0: there is no significance vector mean difference between the different diagnosis groups of the fractal dimension Df

 Name of the statistic’s Statistic p-value Wilks’ Lamda 0.4918 0.0342 Pillai’s Trace 0.5404 0.0499 Hotelling-Lawley Trace 0.9677 0.0236 Roy’s Greatest Root 0.8945 0.0041

The univariate test for individual dependent variables with the null hypothesis about equality of the means of the different groups, are not significant at 5% level (see Table 5.3).

Table-5.3. Univariate analysis of variance for individual response variables with   H0: there have no mean differences between the different diagnoses of the individual fractal dimension Df.

 Variables Between Group MSS Between Group MSS F Statistics (p-values) Bartlett's Chi-Square Fp1_Phase I-III 0.0061 0.0032 1.88(0.173) 1.299 (0.522) Fp1_Phase IV 0.0126 0.0192 0.66(0.526) 0.702(0.704) Fp2_Phase I-III 0.0090 0.0034 2.68(0.089) 0.985(0.611) Fp2_Phase IV 0.0484 0.0290 1.67(0.209) 0.333(0.847)

5.3. Analysis of Contrasts

The overall multivariate test showed significant (see the Section 3.2), which means that considerable differences between different groups exit. We proceed with studying pair wise difference and some predefined contrasts. We set the multivariate test of PD versus the average of BD and UD, then BD versus the average of PD and UD, and UD versus the average of PD and BD. To find where the differences are, some contrasts shown in Table 5.4 are studied. PD shows statistical difference with the average of BD and UD, BD is also different from the average of PD and UD but UD is not significantly differ from the average of PD and BP. Further the pair wise comparison shows that in every pair the average are statistically different from one diagnosis group to another except (PD, UP) pair.

Table-5.4. Test statistics for different types of contrasts testing equality of group means of the fractal dimension Df.

 Contrasts Wilk’s L Pillai’s Trace Hotelling-Lawley Trace Roy’s Greatest Root F Statistic p-value Contrasts with three groups of diagnoses 2(PD)-BD-UP=0 0.6169 0.3831 0.6211 0.6211 3.42 0.0257 2(BD)-PD-UP=0 0.5410 0.4590 0.8484 0.8484 4.67 0.0070 2(UD)-PD-BP=0 0.8511 0.1489 0.1750 0.1750 0.96 0.4476 Pair wise contrasts PD-BD=0 0.5331 0.4669 0.8757 0.8757 4.82 0.0061 PD-UP=0 0.8613 0.1387 0.1611 0.1611 0.89 0.4886 BD-UP=0 0.6305 0.3695 0.5859 0.5859 3.22 0.0317

PD, BD and UP stand for the mean vector of fractal dimension for psychotic depression, bipolar disorder and unipolar depression respectively.

5.4. Multivariate Modeling of Fractal Dimension

For the modeling of fractal dimension Df  will use the general multivariate model defined in Section 3.2,  i.e.

X= BC+E,

column of X represents one individual in Phase I-III and Phase IV for FP1 and FP2, respectively.

Here, the data matrix X is the following

The columns here represent 4 repeated measurement on patients, i.e. Fp1_phase I-III,  Fp1_phase IV,  Fp2_phase I-III,  Fp2_phase IV.

The between-individual design matrix C corresponds to the 3 diagnoses groups is

Now, we would like to test the following hypothesis i.e. equality of two multivariate means,

Under the null hypothesis the test statistic (described in the Section 3.2)

We obtain the following results:

The estimated total sum squared matrix described in (10),

where

0.046. Thus the null hypothesis is rejected at 5% level of significance which implies that there is a significant difference in the average Df between the two points Fp1 and Fp2.

The second hypothesis to test is

H02: There is no difference the average Df between the two phases (Phase I-III and Phase IV), i.e. which is equivalent to

Under the null hypothesis the test statistic

Now the third null hypothesis is;

H03: There is no difference between the three diagnostic groups, i.e. BF= 0, where B and F is defined in Section 3.3.

The following test statistic is used

Under H0 it approximately follows chi-squared distribution wit rt DF as f is large.

significant difference in mean response of Df between the three diagnoses.

Next null hypothesis is;

The following test statistic is used

Under H0 it approximately follows chi-squared distribution wit rt DF as f is large.

interaction effect between electrodes locations and the diagnostic groups.

Next null hypothesis of our interest is;

H05: There is no interaction effect between phases (Phase I-III and Phase IV) and the diagnostic groups, i.e.

The following test statistic is used

Under H0 it approximately follows chi-squared distribution wit rt DF as f is large.

5.5. Assessing Equality of Group Covariance Matrices

The value of Box’s M statistic to test the equality of these three group covariance matrices is 17.82 with p-value greater than 0.10, indicates that there exists no substantive difference between the variances of different response groups.

Now we will summarize the results of this study.

## 6. DISCUSSION

In this work, multivariate analysis of EEG data comprises multivariate analysis of variance and discriminant analysis. We also assessed local complexity of EEG signals by fractal dimension. Usually in ECT, data for EEG is obtained from a small sample. The main goal of the master thesis was to work on diagnostic of the multivariate model that we used to study the complexity of EEG. In particular, the performance of a newly proposed method for identifying outliers based on the eigenvalues of the sample covariance matrix was investigated here, and compared to Mahalanobis distance as based method.

The descriptive statistics showed that the mean values of fractal dimension Df for phase I-III are always lower than that for phase IV, also with lower SE values. That implies the low variability of Df with low mean are observed for the both channels Fp1 and Fp2 for Phase I-III. At the 5% level of significance, mean vectors for the three group diagnoses are statistically different.

There is a significant difference in mean response of Df between the electrodes locations Fp1 and Fp2 as well as significant difference between the phases Phase I-III and PhaseIV. Furthermore, it has been shown that there is a significant difference in mean response of Df between the three diagnoses. There have been detected significant interaction effects between electrode locations and diagnoses. Significant interaction effects were found between phases and diagnoses group.

Using graphical method, the response variables also do not deviate from multivariate normality. Insignificant Mardias statistic based on skewness and kurtosis showed that there is no significant deviation from the multivariate normal distribution. Chi-square (c2) test statistic based on skewness and kurtosis shows all the response variables are normal individually except Fp1_PhaseIV. The value of Box’s M statistic to test the equality of these three group covariance matrices is not significant, indicates that there exists no substantive difference between the variances of different response groups.

The outlying observations were identified based on Mahalanobis distance and MED method. By MD method, there are identified two multivariate outlying observations whereas using MED method four multivariate outlying observations were found. In general, one could see that MED works better compare to MD method as we can identify the outlying observations at first time calculation whereas in MD method needs fourth time calculation. This result confirms [31] results. According to Gao, et al. [31] MED method is a quite powerful tool for identifying multivariate outlying observations. This thesis also gives mathematical ideas from multivariate statistical analysis to find a quick solution or a good approximation of a complex scientific problem.  Further simulation studies are needed to investigate the performance of MED method in case of small samples with covariance matrices having different structures, for example with specific pattern or heteroscedastic covariance matrices.

 Funding: This study received no specific financial support. Competing Interests: The author declares that there are no conflicts of interests regarding the publication of this paper. Contributors/Acknowledgement: I would like to express the deepest appreciation to Dr Tatjana von Rosen, Associate Professor, Stockholm University, Sweden who always gives me a great support in regard to this work.

## REFERENCES

[1]         W. J. Freeman, Societies of brains. Hillsdale, NJ: Erlbaum, 1995.

[2]         L. V. M. Quyen, J. Martinerie, M. Baulac, and F. Varela, "Anticipating epileptic seizures in real time by a non-linear analysis of similarity between EEG recordings," Neuroreport, vol. 10, pp. 2149-2155, 1999. View at Google Scholar | View at Publisher

[3]         A. J. MacLennan, P. R. Carney, W. J. Zhu, A. H. Chaves, J. Garcia, J. R. Grimes, K. J. Anderson, S. N. Roper, and N. Lee, "An essential role for the H218/AGR16/Edg-5/LPB2 sphingosine 1-phosphate receptor in neuronal excitability," European Journal of Neuroscience, vol. 14, pp. 203-209, 2001. View at Google Scholar | View at Publisher

[4]         K. Kaneko, Complex system: Chaos and beyond: A constructive approach with applications in life sciences. Berlin: Springer, 2001.

[5]         L. Diambra, J. B. de Figueiredo, and C. P. Malta, "Epileptic activity recognition in EEG recording," Physica A: Statistical Mechanics and its Applications, vol. 273, pp. 495-505, 1999. View at Google Scholar | View at Publisher

[6]         H. Haken, Principles of brain functioning: A synergetic approach to brain activity, behaviour and cognition. Berlin, Heidelberg, New York: Springer, 1996.

[7]         C. E. Elger and K. Lehnertz, "Seizure prediction by non-linear time series analysis of brain electrical activity," European Journal of Neuroscience, vol. 10, pp. 786-789, 1998. View at Google Scholar | View at Publisher

[8]         T. Ozaki, P. V. Sosa, and V. Haggan-Ozaki, "Reconstructing the nonlinear dynamics of epilepsy data using nonlinear time series analysis," Journal of Signal Processing, vol. 3, pp. 153-162, 1999. View at Google Scholar

[9]         T. Frank, A. Daffertshofer, C. Peper, P. Beek, and H. Haken, "Towards a comprehensive theory of brain activity: Coupled oscillator systems under external forces," Physica D: Nonlinear Phenomena, vol. 144, pp. 62-86, 2000. View at Google Scholar

[10]       R. G. Andrzejak, G. Widman, K. Lehnertz, C. Rieke, P. David, and C. Elger, "The epileptic process as nonlinear deterministic dynamics in a stochastic environment: An evaluation on mesial temporal lobe epilepsy," Epilepsy Research, vol. 44, pp. 129-140, 2001. View at Google Scholar | View at Publisher

[11]       K. Lehnertz, J. Arnold, P. Grassberger, and C. E. Elger, In: Lehnertz,  K., Arnold,  J., Grassberger , P., Elger,  C.E. (Eds.). Chaos in Brain? Singapore: World Scienti Mc, 2000.

[12]       L. D. Iasemidis, P. Pardalos, J. C. Sackellares, and D.-S. Shiau, "Quadratic binary programming and dynamical system approach to determine the predictability of epileptic seizures," Journal of Combinatorial Optimization, vol. 5, pp. 9-26, 2001. View at Google Scholar

[13]       A. Accardo, M. Affinito, M. Carrozzi, and F. Bouquet, "Use of the fractal dimension for the analysis of electroencephalographic time series," Biological Cybernetics, vol. 77, pp. 339-350, 1997. View at Google Scholar | View at Publisher

[14]       G. Henderson, E. Ifeachor, H. Wimalaratna, E. Allen, and N. Hudson, "Prospects for routine detection of dementia using the fractal dimension of the human electroencephalogram," IEE Proceedings-Science, Measurement and Technology, vol. 147, pp. 321-326, 2000. View at Google Scholar | View at Publisher

[15]       R. Saji and H. Konno, "Dynamical features of the local fractal dimension of brain waves and its applicability for diagnosis of senile dementia," Japanese Journal of Applied Physics, vol. 39, pp. 679-684, 2000. View at Google Scholar | View at Publisher

[16]       A. P. Accardo, M. Affinito, M. Carrozzi, S. Cisint, and F. Bouquet, "Comparison between spectral and fractal EEG analyses of sleeping newborns," in Conference Proceedings - IEEE Engineering in Medicine and Biology Society, 1998, pp. 1569-1571.

[17]       X. Li, J. Polygiannakis, P. Kapiris, A. Peratzakis, K. Eftaxias, and X. Yao, "Fractal spectral analysis of pre-epileptic seizures in terms of criticality," Journal of Neural Engineering, vol. 2, pp. 11-16, 2005. View at Google Scholar | View at Publisher

[18]       D. Popivanov, S. Jivkova, V. Stomonyakov, and G. Nicolova, "Effect of independent component analysis on multifractality of EEG during visual-motor task," Signal Processing, vol. 85, pp. 2112-2123, 2005. View at Google Scholar | View at Publisher

[19]       B. Wahlund, W. Klonowski, P. Stepien, R. Stepien, T. von Rosen, and D. von Rosen, "EEG data, fractal dimension and multivariate statistics," Journal of Computer Science and Engineering, vol. 3, pp. 10-14, 2010. View at Google Scholar

[20]       C. Gómez, Á. Mediavilla, R. Hornero, D. Abásolo, and A. Fernández, "Use of the Higuchi's fractal dimension for the analysis of MEG recordings from Alzheimer's disease patients," Medical Engineering & Physics, vol. 31, pp. 306-313, 2009. View at Google Scholar | View at Publisher

[21]       B. Wahlund, P. Piazza, D. von Rosen, B. Liberg, and H. Liljenström, "Seizure (Ictal)—EEG characteristics in subgroups of depressive disorder in patients receiving electroconvulsive therapy (ECT)—a preliminary study and multivariate approach," Computational Intelligence and Neuroscience, vol. 2009, p. 1-8, 2009. View at Google Scholar | View at Publisher

[22]       K. P. Weinfurt, Multivariate analysis of variance. In L. Grimm and P. Yarnold (Eds), Reading and understanding multivariate statistics. Washington DC: American Psychological Association, 1995.

[23]       C. J. Huberty and M. D. Petoskey, Multivariate analysis of variance and covariance. In H. Tinsley and S. Brown (Eds.), Handbook of applied multivariate statistics and mathematical modeling. NY: Academic Press, 2000.

[24]       J. Wackermann, "Beyond mapping: Estimating complexity of multichannel EEG recordings," Acta Neurobiologiae Experimentalis, vol. 56, pp. 197-208, 1996. View at Google Scholar

[25]       T. Higuchi, "Approach to an irregular time series on the basis of the fractal theory," Physica D: Nonlinear Phenomena, vol. 31, pp. 277-283, 1988. View at Google Scholar | View at Publisher

[26]       R. Ferenets, T. Lipping, A. Anier, V. Jantti, S. Melto, and S. Hovilehto, "Comparison of entropy and complexity measures for the assessment of depth of sedation," IEEE Transactions on Biomedical Engineering, vol. 53, pp. 1067-1077, 2006. View at Google Scholar

[27]       N. Kannathal and S. Krishnan, "Comprehensive analysis of cardiac health using heart rate signals," Physiological Measurement, vol. 25, pp. 1139-1151, 2004. View at Google Scholar | View at Publisher

[28]       T. W. Anderson, An introduction to multivariate statistical analysis, 3rd ed. New York: Wiley, 2003.

[29]       T. Kollo and V. D. Rosen, "A unified approach to the approximation of multivariate densities," Scandinavian Journal of Statistics, vol. 25, pp. 93-109, 1998.View at Google Scholar | View at Publisher

[30]       M. S. Srivastava, Method of multivariate statistics: John Wiley & Sons, Inc, 2002.

[31]       S. Gao, G. Li, and D. Wang, "A new approach for detecting multivariate outliers," Communications in Statistics—Theory and Methods, vol. 34, pp. 1857-1865, 2005. View at Google Scholar | View at Publisher