The numerical codes for simulating a thalamic spiking network and mean-field as used in the paper 'A multi-scale study of thalamic responsiveness' are represented here. Two populations of thalamocortical relay neurons (TC) and thalamic reticular neurons (RE) are simulated. Two paramater sets are initiated corresponding to the 'awake' and 'sleep' states defined via the concentration of the neuromodulator Acetycholine (ACh) in the thalamus (see paper for more detail).
The files for simulating spiking network and mean-field respectively are: TNetwork.py, MF_sctipt.py.
The file for fitting the transfer function F on (simulated) single cell data is: fitTF.ipynb, with ExpTF.ipynb the single cell scan simulation for the fit data.