Nonparametric collective spectral density estimation with an application to clustering the brain signals

Stat Med. 2018 Dec 30;37(30):4789-4806. doi: 10.1002/sim.7972. Epub 2018 Sep 26.

Abstract

In this paper, we develop a method for the simultaneous estimation of spectral density functions (SDFs) for a collection of stationary time series that share some common features. Due to the similarities among the SDFs, the log-SDF can be represented using a common set of basis functions. The basis shared by the collection of the log-SDFs is estimated as a low-dimensional manifold of a large space spanned by a prespecified rich basis. A collective estimation approach pools information and borrows strength across the SDFs to achieve better estimation efficiency. Moreover, each estimated spectral density has a concise representation using the coefficients of the basis expansion, and these coefficients can be used for visualization, clustering, and classification purposes. The Whittle pseudo-maximum likelihood approach is used to fit the model and an alternating blockwise Newton-type algorithm is developed for the computation. A web-based shiny App found at "https://ncsde.shinyapps.io/NCSDE" is developed for visualization, training, and learning the SDFs collectively using the proposed technique. Finally, we apply our method to cluster similar brain signals recorded by the for identifying synchronized brain regions according to their spectral densities.

Keywords: Whittle likelihood; roughness penalty; time series clustering.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms
  • Brain / physiology*
  • Cluster Analysis*
  • Data Interpretation, Statistical
  • Electroencephalography* / methods
  • Fourier Analysis
  • Humans
  • Likelihood Functions
  • Models, Statistical
  • Normal Distribution
  • Statistics, Nonparametric*
  • Time Factors