Working memory (WM) is a mechanism that temporarily stores and manipulates information in service of behavioral goals and is a highly dynamic process. Previous studies have considered decoding WM load using EEG but have not investigated the contribution of sequential information contained in the temporal patterns of the EEG data that can differentiate different WM loads. In our study, we develop a novel method of investigating the role of sequential information in the manipulation and storage of verbal information at various time scales and localize topographically the sources of the sequential information based decodability. High density EEG (128-channel) were recorded from twenty subjects performing a Sternberg verbal WM task with varying memory loads. Long Short-Term Memory Recurrent Neural Networks (LSTM-RNN) were trained to decode memory load during encoding, retention, activity-silent, and retrieval periods. Decoding accuracy was compared between ordered data and a temporally shuffled version that retains pattern based information of the data but not temporal relation to assess the contribution of sequential information to decoding memory load. The results show that (1) decoding accuracy increases with increase in the length of the EEG time series given to the LSTM for both ordered and temporally shuffled cases, with the increase being faster for ordered than temporally shuffled time series, and (2) according to the decoding weight maps, the frontal, temporal and some parietal areas are an important source of sequential information based decodability. This study, to our knowledge, is the first study applying a LSTM-RNN approach to investigate temporal dynamics in human EEG data in encoding WM load information.
High-density EEG data was recorded while the subject was performing the Sternberg working meory task. The recording took place in an electrically and acoustically shielded room with a 128-channel BioSemi Active Two System at a 1 KHz sampling rate. As shown in (A) paradigm for the verbal working memory task. (B) Time period of interest for the working memory stages.
Time course of decoding accuracy for load 2 vs load 6 conditions with SVM during the working memory task. Decoding accuracy is used to assign time periods of interest for the four memory stages in which time series are constructed for LSTM decoding. Shaded region shows the standard error across subjects.
LSTM decoding. (A) LSTM decoding accuracy for all memory stages across different input time series lengths. The decoding accuracy of the ordered data (thicker line) is compared to the temporally shuffled data (thinner line). The horizontal lines at the bottom represent time lengths associated with decoding accuracies that are significantly different via a two-sample t-test (p < 0.05) in the ordered and shuffled scenario. (B) The difference between the decoding accuracy of the ordered and shuffled data across variable time lengths.
For each memory stage, the topographical plot of weights for the temporally shuffled scenario is subtracted from the topographical plot of weights for the ordered scenario. A z-score of the topographical plot is computed.