A Skew-t space-varying regression model for the spectral analysis of resting state brain activity

It is known that in many neurological disorders such as Down syndrome, main brain rhythms shift their frequencies slightly, and characterizing the spatial distribution of these shifts is of interest. This article reports on the development of a Skew-t mixed model for the spatial analysis of resting state brain activity in healthy controls and individuals with Down syndrome. Time series of oscillatory brain activity are recorded using magnetoencephalography, and spectral summaries are examined at multiple sensor locations across the scalp. We focus on the mean frequency of the power spectral density, and use space-varying regression to examine associations with age, gender and Down syndrome across several scalp regions. Spatial smoothing priors are incorporated based on a multivariate Markov random field, and the markedly non-Gaussian nature of the spectral response variable is accommodated by the use of a Skew-t distribution. A range of models representing different assumptions on the association structure and response distribution are examined, and we conduct model selection using the deviance information criterion. 1 Our analysis suggests region-specific differences between healthy controls and individuals with Down syndrome, particularly in the left and right temporal regions, and produces smoothed maps indicating the scalp topography of the estimated differences.

[1]  Neal Alexander,et al.  Bayesian Disease Mapping: Hierarchical Modeling in Spatial Epidemiology , 2011 .

[2]  A. Azzalini,et al.  Statistical applications of the multivariate skew normal distribution , 2009, 0911.2093.

[3]  Hao Zhang,et al.  Kriging and cross‐validation for massive spatial data , 2009 .

[4]  Hao Zhang,et al.  On spatial skew‐Gaussian processes and applications , 2009 .

[5]  Ernesto San Martín,et al.  Linear mixed models with skew-elliptical distributions: A Bayesian approach , 2008, Comput. Stat. Data Anal..

[6]  Richard J. Haier,et al.  Neuroimaging of individuals with Down’s syndrome at-risk for dementia: Evidence for possible compensatory events , 2008, NeuroImage.

[7]  Mark F. J. Steel,et al.  Non-Gaussian Bayesian Geostatistical Modeling , 2006 .

[8]  B. Mallick,et al.  A Bayesian prediction using the skew Gaussian distribution , 2004 .

[9]  Renato M. Assunção,et al.  Space varying coefficient models for small area data , 2003 .

[10]  A. Azzalini,et al.  Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t‐distribution , 2003, 0911.2342.

[11]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[12]  Y. MacNab,et al.  Autoregressive Spatial Smoothing and Temporal Spline Smoothing for Mapping Rates , 2001, Biometrics.

[13]  Dani Gamerman,et al.  Space-varying regression models: specifications and simulation , 2001, Comput. Stat. Data Anal..

[14]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[15]  Andrew Thomas,et al.  WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility , 2000, Stat. Comput..

[16]  M. Steel,et al.  On Bayesian Modelling of Fat Tails and Skewness , 1998 .

[17]  H. Wiśniewski,et al.  Alzheimer's disease in Down's syndrome , 1985, Neurology.

[18]  Q. A.,et al.  Linear effects mixed models with skew-elliptical distributions : A Bayesian approach , 2008 .

[19]  P. Embrechts Skew-Elliptical Distributions and Their Applications: A Journey Beyond Normality , 2005 .

[20]  A. Gelfand,et al.  Proper multivariate conditional autoregressive models for spatial data analysis. , 2003, Biostatistics.

[21]  Bradley P. Carlin,et al.  Hierarchical Multivarite CAR Models for Spatio-Temporally Correlated Survival Data , 2002 .