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

 Download zip file 
Help downloading and running models
Generates typical EEG data of sleeping Humans for sleep stages N2/N3 as well as wakefulness
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;
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];
Search NeuronDB for information about:  Neocortex L2/3 pyramidal GLU cell; AMPA; Gaba; I_K,Na; Na/K pump;
function Data_Time_Series(mode)
% This function creates the model data depicted in Figure 5-6 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

% To ensure availability of the simulation routine it should be called from Create_Data.m

% Parameter settings
case 1:
	Param_N3 = [ 6.5;	% sigma_e
		     2;		% g_KNa
		     0E-3];	% dphi

case 2:
	Param_N2 = [ 4.6;	% sigma_e
		     1.33;	% g_KNa
		     0E-3];	% dphi

% No stimulation so set everything to zero
var_stim      = zeros(8,1);
% Duration of the stimulation
T     	      = 15;

% Run the two simulations
[Ve_N2{i}, ~] = Cortex_mex(T, Param_N2, var_stim);
[Ve_N3{i}, ~] = Cortex_mex(T, Param_N3, var_stim);

save('Data/EEG_Data_N2.mat', 'Ve_N2');
save('Data/EEG_Data_N3.mat', 'Ve_N3');

Loading data, please wait...