Neural mass model of the sleeping cortex (Weigenand et al 2014)

 Download zip file 
Help downloading and running models
Accession:226472
Generates typical EEG data of sleeping Humans for sleep stages N2/N3 as well as wakefulness
Reference:
1 . Weigenand A, Schellenberger Costa M, Ngo HV, Claussen JC, Martinetz T (2014) Characterization of K-complexes and slow wave activity in a neural mass model. PLoS Comput Biol 10:e1003923 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neural mass;
Brain Region(s)/Organism: Neocortex;
Cell Type(s): Neocortex L2/3 pyramidal GLU cell; Neocortex layer 2-3 interneuron;
Channel(s): I_K,Na; Na/K pump;
Gap Junctions:
Receptor(s): AMPA; Gaba;
Gene(s):
Transmitter(s):
Simulation Environment: C or C++ program (web link to model); MATLAB;
Model Concept(s): Bifurcation; Sleep; Activity Patterns;
Implementer(s): Schellenberger Costa, Michael [mschellenbergercosta at gmail.com];
Search NeuronDB for information about:  Neocortex L2/3 pyramidal GLU cell; AMPA; Gaba; I_K,Na; Na/K pump;
function Plot_Compare_N2()
% This function creates Figure 5 of
%
% Characterization of K-Complexes and Slow Wave Activity in a Neural Mass Model
% A Weigenand, M Schellenberger Costa, H-VV Ngo, JC Claussen, T Martinetz
% PLoS Computational Biology. 2014;10:e1003923
% Load the data

% Path to utility functions
if(isempty(strfind(path, [pwd, '/Tools'])))
    addpath([pwd, '/Tools']);
end

% Load the data
load('Data/EEG_Data_N2.mat')
load('Data/Simulation_N2.mat')

% Zscore the simulation data
Ve_N2 = zscore(Ve_N2); %#ok<NODEF>

% Set the time axis
timeaxis = linspace(0, 30, 3000);

% Create figure handle
figure(1)
clf,shg

% Plot the data
plot(timeaxis,Data_N2+4,'-k','LineWidth',0.5)
hold on
plot(timeaxis,Ve_N2-4,'-k','LineWidth',0.5);

% Remove the axes 
axis off

% Add the xlabel and the figure caption
text(-1.5,7,'(a)')
text(-1.5,-2,'(b)')

% Add a 5s scale
plot([2,7], [-9,-9],'-k','LineWidth',2)
text(4.5,-10,'5s')
hold off
end

Loading data, please wait...