{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "0b85a598",
"metadata": {},
"outputs": [],
"source": [
"from brian2 import *\n",
"\n",
"from adex_sine import *\n",
"\n",
"defaultclock.dt = 10 * us\n",
"\n",
"\n",
"class Model:\n",
" # Model type flags\n",
" SINE = 0\n",
" EXP2SYN = 1\n",
"\n",
" # Noise flags\n",
" HIGH = 0\n",
" LOW = 1\n",
" OFF = -1\n",
"\n",
" # Noise parameters\n",
" EXCITATORY_NOISE_VARIANCE = {HIGH: 0.5 * nS, LOW: 0.25 * nS, OFF: 0 * nS}\n",
" INHIBITORY_NOISE_VARIANCE = {HIGH: 1.25 * nS, LOW: 0.625 * nS, OFF: 0 * nS}\n",
"\n",
" # Noise mean conductance\n",
" EXCITATORY_CONDUCTANCE = 1 * nS\n",
" INHIBITORY_CONDUCTANCE = 4 * nS\n",
"\n",
" DEFAULT_PARAMETERS = {\n",
" \"sigma_flux\" : 6.75*pA, \n",
" \"c\": 85 * pF,\n",
" \"tau_w\": 18 * ms,\n",
" \"b\": 0.25 * nA,\n",
" \"a\": 1.3 * nS,\n",
" \"v_T\": -45 * mV,\n",
" \"v_thresh\": 0 * mV,\n",
" \"DeltaT\": 0.2 * mV,\n",
" # EQUILIBRIUM POTENTIAL\n",
" \"e_l\": -65 * mV,\n",
" \"e_ex\": 0 * mV,\n",
" \"e_in\": -70 * mV,\n",
" # CONDUCTANCES\n",
" \"g_l\": 3 * nS,\n",
" \"mu_ex\": 0 * nS,\n",
" \"mu_in\": 0 * nS,\n",
" # EXCITATORY NOISE\n",
" \"sigma_ex\": 0 * nS,\n",
" \"tau_noise_ex\": 3 * ms,\n",
" # INHIBITORY NOISE\n",
" \"sigma_in\": 0 * nS,\n",
" \"tau_noise_in\": 10 * ms,\n",
" # SINE INPUT\n",
" \"f\": 100 * Hz,\n",
" \"A\": 0 * pA,\n",
" \"i_injected\": 0 * pA,\n",
" \"v_reset\": -70 * mV,\n",
" # m current\n",
" \"g_adapt\": 10 * nS,\n",
" \"e_k\": -90*mV,\n",
" \"beta_z\": -35*mV,\n",
" \"gamma_z\": 4*mV, #5\n",
" \"tau_z\": 100*ms,\n",
" }\n",
"\n",
" def __init__(\n",
" self, n, *, stim=None, noise=None, resistance=None, additional_vars=()\n",
" ):\n",
" if resistance is None:\n",
" raise ValueError(\"Resistance must be specified\")\n",
"\n",
" if noise is None:\n",
" raise ValueError(\"Noise must be specified\")\n",
"\n",
" self.stim_type = stim\n",
" self._input_resistance = None\n",
" self._noise_level = None\n",
" self._duration = 0\n",
" self.recorded_vars = (\"v\",) + additional_vars\n",
"\n",
" self.neurons = self.set_default(n_neuron=n)\n",
" self.set_resistance(resistance)\n",
" self.set_noise(noise)\n",
"\n",
" self.spikes = None\n",
" self.spiker = None\n",
" self.synapses = None\n",
" self.inhib_synapses = None\n",
" self.smon = None\n",
" self.network = None\n",
" self.build_network()\n",
"\n",
" def create_model(self):\n",
" return ADEX_MODEL, self.DEFAULT_PARAMETERS\n",
"\n",
" def set_default(self, n_neuron):\n",
" model, parameters = self.create_model()\n",
"\n",
" neurons = NeuronGroup(\n",
" n_neuron,\n",
" model=model,\n",
" method=\"Euler\",\n",
" name=\"neurons\",\n",
" threshold=\"v > v_thresh\",\n",
" reset=\"v = v_reset; w += b\",\n",
" )\n",
"\n",
" for parameter, value in parameters.items():\n",
" neurons.__setattr__(parameter, value)\n",
"\n",
" neurons.v = neurons.e_l # remove most of transient\n",
"\n",
" return neurons\n",
"\n",
" def set_resistance(self, level):\n",
" if level == self.LOW:\n",
" exc_conductance = self.EXCITATORY_CONDUCTANCE\n",
" inhib_conductance = self.INHIBITORY_CONDUCTANCE\n",
"\n",
" else:\n",
" exc_conductance = inhib_conductance = 0\n",
"\n",
" self._input_resistance = level\n",
" self._set_variable(\"mu_ex\", exc_conductance)\n",
" self._set_variable(\"mu_in\", inhib_conductance)\n",
"\n",
" def set_noise(self, level):\n",
" if level == self.HIGH or level == self.LOW:\n",
" exc_noise = self.EXCITATORY_NOISE_VARIANCE[level]\n",
" inhib_noise = self.INHIBITORY_NOISE_VARIANCE[level]\n",
"\n",
" else:\n",
" exc_noise = inhib_noise = 0\n",
"\n",
" self._noise_level = level\n",
" self._set_variable(\"sigma_ex\", exc_noise)\n",
" self._set_variable(\"sigma_in\", inhib_noise)\n",
"\n",
" def set_injected_current(self, amplitude):\n",
" self._set_variable(\"i_injected\", amplitude)\n",
" self._set_variable(\"A\", 0 * pA)\n",
"\n",
" def set_stimulus_current(self, amplitude):\n",
" self._set_variable(\"A\", amplitude)\n",
" self._set_variable(\"i_injected\", 0 * pA)\n",
"\n",
" @property\n",
" def f(self):\n",
" return self.neurons.f\n",
"\n",
" @f.setter\n",
" def f(self, new_f):\n",
" self._set_variable(\"f\", new_f) # this will reset smon\n",
" if self.stim_type == self.EXP2SYN:\n",
" self.spiker.T = 1 / new_f\n",
"\n",
" def run(self, duration, report=\"stdout\"):\n",
" self._duration = duration\n",
" self.network.run(duration, report=report)\n",
"\n",
" def build_network(self):\n",
" self.smon = StateMonitor(\n",
" self.neurons, self.recorded_vars, record=True, name=\"smon\"\n",
" )\n",
" self.spikes = SpikeMonitor(self.neurons, name=\"spikes\")\n",
"\n",
" self.network = Network(self.neurons, self.smon, self.spikes)\n",
"\n",
" def _set_variable(self, name, value):\n",
" self.neurons.__setattr__(name, value)\n",
" self.reset_recording()\n",
"\n",
" def reset_recording(self):\n",
" try:\n",
" self.network\n",
" except AttributeError:\n",
" return # network not yet initialized\n",
"\n",
" self.network.remove(self.smon, self.spikes)\n",
"\n",
" self.smon = StateMonitor(\n",
" self.neurons, self.recorded_vars, record=True, name=\"smon\"\n",
" )\n",
" self.spikes = SpikeMonitor(self.neurons, name=\"spikes\")\n",
"\n",
" self.network.add(self.smon, self.spikes)\n",
"\n",
" @property\n",
" def spike_train(self):\n",
" return self.spikes.spike_trains()\n",
"\n",
" @property\n",
" def firing_rate(self):\n",
" return self.spikes.count / self.duration\n",
"\n",
" @property\n",
" def duration(self):\n",
" return self._duration\n",
"\n",
" @property\n",
" def input_resistance(self):\n",
" if self._input_resistance == self.HIGH:\n",
" return \"HIGH\"\n",
" else:\n",
" return \"LOW\"\n",
"\n",
" @property\n",
" def noise_level(self):\n",
" if self._noise_level == self.HIGH:\n",
" return \"HIGH\"\n",
" elif self._noise_level == self.LOW:\n",
" return \"LOW\"\n",
" else:\n",
" return \"NO\"\n",
"\n",
" def __repr__(self):\n",
" return f\"{self.neurons.N} Neurons with {self.input_resistance} input resistance and {self.noise_level} noise\"\n",
"\n",
" def __str__(self):\n",
" return self.__repr__()\n",
"\n",
" def store(self, name):\n",
" self.network.store(name)\n",
"\n",
" def restore(self, name):\n",
" self.network.restore(name)\n",
"\n",
" @property\n",
" def v(self):\n",
" return self.smon.v\n",
"\n",
" @property\n",
" def t(self):\n",
" return self.smon.t\n",
"\n",
" @property\n",
" def injected_current(self):\n",
" return self.neurons.i_injected\n",
"\n",
" @property\n",
" def stimulus_amplitude(self):\n",
" return self.neurons.A\n",
"\n",
"\n",
"class CurrentModel(Model):\n",
" def __init__(self, **kwargs):\n",
" super().__init__(stim=self.SINE, **kwargs)\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += CURRENT_INPUT\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SineModel(CurrentModel):\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += SINE_INPUT\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SawModel(CurrentModel):\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += SAW_INPUT\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SynapticModel(Model):\n",
" def __init__(self, **kwargs):\n",
" super().__init__(stim=self.EXP2SYN, **kwargs)\n",
"\n",
" SYNAPTIC_PARAMETERS = {\n",
" \"tau_input_1\": 0.4 * ms,\n",
" \"tau_input_2\": 4 * ms,\n",
" \"offset_A\": 1.48793507e-11,\n",
" \"offset_B\": -2.66359562e-08,\n",
" \"offset_C\": 1.77538800e-05,\n",
" \"offset_D\": -8.05925810e-04,\n",
" \"offset_E\": -3.51463644e-02,\n",
" \"offset_switch\": 0,\n",
" }\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += EXP2SYN_WAVEFORM + SUMMATION_OFFSET\n",
" parameters = {**parameters, **self.SYNAPTIC_PARAMETERS}\n",
"\n",
" return model, parameters\n",
"\n",
" def build_network(self):\n",
" super().build_network()\n",
" self.spiker = NeuronGroup(\n",
" self.neurons.N,\n",
" \"\"\"T : second (constant)\n",
" lastspike : second\"\"\",\n",
" threshold=\"timestep(t-lastspike, dt)>=timestep(T, dt)\",\n",
" reset=\"lastspike=t\",\n",
" )\n",
" self.spiker.T = 1 / self.neurons.f\n",
" self.synapses = Synapses(\n",
" self.spiker, self.neurons, on_pre=\"input_aux += 1\"\n",
" ) # connect input to neurons\n",
" self.synapses.connect(\"i==j\") # one synapse goes to one neuron\n",
"\n",
" self.network.add(self.spiker, self.synapses)\n",
"\n",
"\n",
"class SynapticCurrentModel(SynapticModel):\n",
" def __init__(self, offset=True, **kwargs):\n",
" self.offset = 1 if offset else 0\n",
" super().__init__(**kwargs)\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += CURRENT_INPUT + SYNAPTIC_INPUT_CURRENT\n",
" parameters = {**parameters, **{\"offset_switch\": self.offset}}\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SynapticConductanceModel(SynapticModel):\n",
" FLAT = 0\n",
" ACTIVE = 1\n",
"\n",
" CONDUCTANCE_PARAMETERS = {\n",
" \"A\": 0 * nS, # overwrite A to be conductance\n",
" \"g_i\": 1 * nS,\n",
" }\n",
"\n",
" INHIBITION_PARAMETERS = {\n",
" \"tau_inhibition_1\": 1 * ms,\n",
" \"tau_inhibition_2\": 10 * ms,\n",
" }\n",
"\n",
" def __init__(self, offset=ACTIVE, **kwargs):\n",
" self.offset = offset\n",
" super().__init__(**kwargs)\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" if self.offset == self.FLAT:\n",
" model += CONDUCTANCE_INPUT + SYNAPTIC_CONDUCTANCE_FLAT\n",
" parameters = {\n",
" **parameters,\n",
" **self.SYNAPTIC_PARAMETERS,\n",
" **self.CONDUCTANCE_PARAMETERS,\n",
" **{\"offset_switch\": 1},\n",
" }\n",
"\n",
" elif self.offset == self.ACTIVE:\n",
" model += CONDUCTANCE_INPUT + SYNAPTIC_CONDUCTANCE_STIM\n",
" parameters = {\n",
" **parameters,\n",
" **self.SYNAPTIC_PARAMETERS,\n",
" **self.CONDUCTANCE_PARAMETERS,\n",
" **self.INHIBITION_PARAMETERS,\n",
" }\n",
"\n",
" return model, parameters\n",
"\n",
" def build_network(self):\n",
" super().build_network()\n",
" if self.offset != self.ACTIVE:\n",
" return\n",
"\n",
" self.inhib_synapses = Synapses(\n",
" self.spiker, self.neurons, on_pre=\"input_inhib_aux += 1\", delay=2 * ms\n",
" ) # connect input to neurons\n",
" self.inhib_synapses.connect(\"i==j\") # one synapse goes to one neuron\n",
"\n",
" self.network.add(self.inhib_synapses)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "c82769f2",
"metadata": {
"scrolled": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING Cannot use Cython, a test compilation failed: Microsoft Visual C++ 14.0 or greater is required. Get it with \"Microsoft C++ Build Tools\": https://visualstudio.microsoft.com/visual-cpp-build-tools/ (DistutilsPlatformError) [brian2.codegen.runtime.cython_rt.cython_rt.failed_compile_test]\n",
"INFO Cannot use compiled code, falling back to the numpy code generation target. Note that this will likely be slower than using compiled code. Set the code generation to numpy manually to avoid this message:\n",
"prefs.codegen.target = \"numpy\" [brian2.devices.device.codegen_fallback]\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Starting simulation at t=0. s for a duration of 12. s\n",
"0.56071 s (4%) simulated in 10s, estimated 3m 24s remaining.\n",
"1.0924 s (9%) simulated in 20s, estimated 3m 20s remaining.\n",
"1.60078 s (13%) simulated in 30s, estimated 3m 15s remaining.\n",
"2.12902 s (17%) simulated in 40s, estimated 3m 5s remaining.\n",
"2.6601 s (22%) simulated in 50s, estimated 2m 56s remaining.\n",
"3.19165 s (26%) simulated in 1m 0s, estimated 2m 46s remaining.\n",
"3.72407 s (31%) simulated in 1m 10s, estimated 2m 36s remaining.\n",
"4.25818 s (35%) simulated in 1m 20s, estimated 2m 25s remaining.\n",
"4.78971 s (39%) simulated in 1m 30s, estimated 2m 16s remaining.\n",
"5.31406 s (44%) simulated in 1m 40s, estimated 2m 6s remaining.\n",
"5.84306 s (48%) simulated in 1m 50s, estimated 1m 56s remaining.\n",
"6.37003 s (53%) simulated in 2m 0s, estimated 1m 46s remaining.\n",
"6.90143 s (57%) simulated in 2m 10s, estimated 1m 36s remaining.\n",
"7.43018 s (61%) simulated in 2m 20s, estimated 1m 26s remaining.\n",
"7.95861 s (66%) simulated in 2m 30s, estimated 1m 16s remaining.\n",
"8.48854 s (70%) simulated in 2m 40s, estimated 1m 6s remaining.\n",
"9.01487 s (75%) simulated in 2m 50s, estimated 56s remaining.\n",
"9.54432 s (79%) simulated in 3m 0s, estimated 46s remaining.\n",
"10.07578 s (83%) simulated in 3m 10s, estimated 36s remaining.\n",
"10.59778 s (88%) simulated in 3m 20s, estimated 26s remaining.\n",
"11.11216 s (92%) simulated in 3m 30s, estimated 17s remaining.\n",
"11.63647 s (96%) simulated in 3m 40s, estimated 7s remaining.\n",
"12. s (100%) simulated in 3m 47s\n"
]
},
{
"data": {
"text/plain": [
"[<matplotlib.collections.EventCollection at 0x19dd5321490>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5331a00>,\n",
" <matplotlib.collections.EventCollection at 0x19dd533d1f0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53459d0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53511c0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5359970>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5361f40>,\n",
" <matplotlib.collections.EventCollection at 0x19dd536d4f0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5375af0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5381370>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5389cd0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5394640>,\n",
" <matplotlib.collections.EventCollection at 0x19dd539eca0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53aa3a0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53b1c70>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53bd4f0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53c7eb0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53d3760>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53dbe80>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53e7640>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53eeca0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd53fa280>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5402b20>,\n",
" <matplotlib.collections.EventCollection at 0x19dd540e4f0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54179d0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54231c0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd542a8b0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5432f40>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5441670>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5446d30>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54536d0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd545cf10>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54687f0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54741c0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd547c850>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5488250>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54919a0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd549d070>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54a46d0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54adfd0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54b9580>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54c0d30>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54cc550>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54d4b80>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54e0280>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54e98e0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54f5100>,\n",
" <matplotlib.collections.EventCollection at 0x19dd54fc7f0>,\n",
" <matplotlib.collections.EventCollection at 0x19dd5505e80>,\n",
" <matplotlib.collections.EventCollection at 0x19dd55116a0>]"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD4CAYAAAAaT9YAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAijUlEQVR4nO2dbaxmVXXH/6sMiEUpM2UgUzAdTAgp0VTsrXFqY6iIWEuENrVgqp2mtGPSNtHaRIb6oeMn0TbGNCYVojbT1heoL2VCbJGgpDEhyCW+8DoOviF16lylEbUfKrr64Z6Dmz1777XWOc+9c/fT/y+Z3PPss/faa+1znp15/nudfURVQQghpE9+5kQ7QAghZDqcxAkhpGM4iRNCSMdwEieEkI7hJE4IIR2zbTM7O/PMM3X37t2b2SUhhHTPvffe+x1V3Vk6t6mT+O7du7G6urqZXRJCSPeIyDdq5yinEEJIx3ASJ4SQjuEkTgghHcNJnBBCOoaTOCGEdIwrO0VEvg7g+wB+DOBJVV0RkR0AbgKwG8DXAfyeqv73xrhJCCGkROR/4r+hqi9Q1ZXh834Ad6jq+QDuGD4TQgjZRObIKVcAODgcHwRw5WxvCCGEhPBO4grgUyJyr4jsG8rOVtWjADD8PavUUET2iciqiKyura3N95gQQshTeJ/YfImqfktEzgJwu4g87O1AVW8EcCMArKys8A0UhBCyQFyTuKp+a/h7TEQ+AeBFAL4tIrtU9aiI7AJwbAP9fIrnH7gNAHDfgcuOO3fVDXcBAG56w57i57ysdDxSKstJ7dZ8yetafY51Hjz6BC7cdbpZrxZfSmqr5uNYx0PLt5qPrfOpvfH6jr6U7NTiKeGJMb8XWn7X4iv5Yt1/pbJxPAA8zVfP/Vqi5V96vhaXB+veqfleux+t61G6/tb1y8es5XPrPin1VYohb7ORmHKKiJwmIs8ejwG8AsD9AA4B2DtU2wvglo1ykhBCSBnP/8TPBvAJERnrf0hV/11E7gFws4hcA+BRAK/ZODcJIYSUMCdxVf0qgF8ulH8XwCUb4RQhhBAfsplvu19ZWdGpW9G29MCRlm7o0XFb5DpcquXW/Ci18fidM0Ujj2qhFpbuXuqrpUt69cVUH54Se23sWnqs1bd3vCPx5/Hk56zr713fye3n93BrzEr2Sn7ntmvrH60Yc1ulfku+l+pHNPtSPFNjzW3P0cdF5N7kGZ2nwcfuCSGkYziJE0JIx3ASJ4SQjulCE2/lPqfHnlzZXCOr6Wqlfr35rqU+W8c53jxgS+eM5gbXbJbseuNI2wJojrvHZq1tKUZvvZrPqS+evGLPNapp717Nv9RvflzCimOqXlv7jkS+kyU7Jb88ueN5u+j3yDNGkXt/UVATJ4SQJYWTOCGEdAwncUII6ZhuNXHA1p0srWuqbpa38eba1nT8SO5yar917IkxknPtaVfyu5U/HdHd8xgiOcU1W54845Zvnj08cqbmk9d8bvleWovwjIe1JtHyqzYWLU3f0sVrfeZxWrY8Y2et1aRE17eYJ04IIeQ4OIkTQkjHcBInhJCO6U4Tr+l2U/cAjuhwebuS3liipY17+vDElZbPyTG3sDTrlm7c0jFbevdcPd9jx8pvXkQ+eal9at+TL56zqPWQkn+R8QPq+wsBqOrlrf5qcZY+l+qWYsrbep9N8NbJ+1oU1MQJIWRJ4SROCCEdw0mcEEI6pgtNHJiXQzxFC7XyWluapWe/DMufiE1LZ4/mm7ditXKJp+yvEs3nLsVYitMbe+p3DY/G2dLUa317cp9bPlg6fG1/lVZedw3ruQErR9yjZ+cxeb+PHqw1nFr9NFbPOc/5KNTECSFkSeEkTgghHdONnALEfu5FUq5q5z2P3OZtLf8j26O25Jfc3hSJp/Qz3JMi13oVlScOz6PI3kega6mekZ/hXvnLIz/V+rdkBc/18abh1bDSIC1ppFQ/778l2VjfzZaM4r1nvL7n/kUkoUhaYhrfHCinEELIksJJnBBCOoaTOCGEdEwXmrhH+41qfa2Uxaj2FnlcOHK+Fbcn3dKzBtDSs6e85s2je7bOt2znsVkpeYBPq7Y+13R14PgtFzzbDVv3RMqU9Zk81potz9pGdFxafrTWhCKpkq3YrPPRdSKv7h/pdwrUxAkhZEnhJE4IIR3DSZwQQjqmG018xPuobVQj8+aBWrqtVzP0vAorkrOat7H0vWhud8m+VS+SY7+oMZiynWvJL2+OdqSN596YM1Y1Wusmeb8trbzmq7XelJflfkXWu7xrWhG9vzROgL2FdEq+RuLdbtoLNXFCCFlSOIkTQkjHuCdxETlJRD4vIrcOn3eIyO0icmT4u33j3CSEEFLCrYmLyJsBrAA4XVUvF5F3AnhcVa8Xkf0AtqvqtS0bc7eiBabnbXo1wZpdq60n79irg3v03ryPaHxpnDUbpXoRWq+vy8dlrFfyv7VtqreshlfHnXLfRJ8j8OT3T8k9j4xB/rn13fLaacXr3R+n1o+ls7fwaPH5ce5DKbZaH3OYrYmLyLkAfgvA+5LiKwAcHI4PArhyho+EEEIm4JVT3g3gLQB+kpSdrapHAWD4e1apoYjsE5FVEVldW1ub4yshhJAMcxIXkcsBHFPVe6d0oKo3quqKqq7s3LlziglCCCEVTE1cRN4O4PUAngRwKoDTAXwcwK8CuFhVj4rILgB3quoFLVuLyhMHUNQF07pe/bqmo5byPtN+87peDbTkU+qbp32r71bdkY1YN7B02dyPKdeu5m/Nv1q/rfOprWjucWSfDG/+d65JA+W1A09Oeu2+Tfux1mFKbUv+W+s3U/dCiTyPkNvyjllOa/0mcj/PYZYmrqrXqeq5qrobwNUAPq2qrwNwCMDeodpeALfM9pQQQkiIOXni1wO4VESOALh0+EwIIWQT2RaprKp3ArhzOP4ugEsW7xIhhBAvXe2d0srfTetOyaf2aKC5P0B9/+ha/Vrfrc81Svqi57jW1tNX6TiPs5an7dlL2rufRopXX67h3eci0o+Ve5z2682N9mq20WtuXa+p1z4vt77Hc+4p6/mC2vMINVr3ba2s5Wd+Lgr3TiGEkCWFkzghhHQMJ3FCCOmYLjRxAHj+gdsAAPcduOypsinaqXe/hJTouztrWJq5tzyNBbA1Pk9eb4kp+c65/WjMJZut9QvLd8/+HN5c4dZ18t4j3lzq0udandazAq3708o7z+tb42GNQ+qTZ3/w0t47U55HaPVrtZ2T376ofVMAauKEELK0cBInhJCO4SROCCEd040mXtOzpujJpXqAnVfryWWdu3/E1Hxby24ea0pE/8zjsvLyS3a95a3cX8u3Gp57p+WDpUV7dW1LL209s1BbD4mMp0e79vhSsxfRk6M52V59uhaDFV/p2KrfisnyzQM1cUIIWVI4iRNCSMd0IadMeXy5Vu55fNnji/WYb6t96yd/bavL3Mep6Xb5Z+/j7p5HtRf1cz3H84o0r6wW+fkbkTFar6IrEUkLzMu8MdVs5edSf/L4S3Y9qaOlNrXPXqmwNRYlv2uySy1ey88cT0y5/alQTiGEkCWFkzghhHQMJ3FCCOmYLjRxwJfSFU1DTM+NtDQ6TyqblfZnpfOleLYI8Or+LY3XkwbpSYscfQZg9mFp/6Vx9d4DrfbWOe+4zrmOre0SWmMTWadIfSuNVWsMcnuRFDzv92mOz5YWDvi20Y2uK0XWnqKvLbSgJk4IIUsKJ3FCCOkYTuKEENIxXWjikbxgy4al5Vm5qlFtPmXKK9Ks179Z+dF5X3lZzc6cnFfPGoBXw2z54tnmNqLnRuL23pPW9a3FEr1+rbhL9nIiefRTt4vI8a5RROq26tV88qwVpW29OfWLzBWnJk4IIUsKJ3FCCOkYTuKEENIxXWni0fzQtO14DsDTNOb0c8tuibk6dIkp+cqe/N7aZ++WqICd79zyo6Ttp7RstdYFWjGlfrS06VYOfN5HNCc9qrlbmq6V22/1UfMxmt+d+5XH2HrGAWg/A1HCep7BQ3S9IaLL5z6lME+cEEJIFU7ihBDSMZzECSGkY7rSxIH2niVzclmtfUpyP0p20jq1fFOvVtp6BZa3bYp3r5e0rjcvvjQ2qf3onhOAvU5Ru4ZjW0vjt/a5yPuMvPYt7adVz6uVztHRveslNZ+9+263bNQoXWtLi7b2tR9t5nZLMbVy+0ux5X14XzU3RwsfoSZOCCFLCidxQgjpGHMSF5FTReRzIvJFEXlARN42lO8QkdtF5Mjwd/vGu0sIISTF1MRFRACcpqo/EJGTAXwWwBsB/A6Ax1X1ehHZD2C7ql7bsjV3P3HAzukt1c2Pa3j21/CeL9kez7c061IOdqlexHdPvrj3czRfuNaupGW23jnZype39t+28swt3bfUrharReteafnuHfvo/Tkn992TW537lJZZ6z7e5w5qeNY6Ini+Y7mvi2KWJq7r/GD4ePLwTwFcAeDgUH4QwJXzXSWEEBLBpYmLyEki8gUAxwDcrqp3AzhbVY8CwPD3rErbfSKyKiKra2trC3KbEEII4JzEVfXHqvoCAOcCeJGIPM/bgareqKorqrqyc+fOiW4SQggpEc4TF5G/BvBDAH8C4GJVPSoiuwDcqaoXtNouQhMfsfKwAXv/7WiOaOu8Nz+8RSnP2dp3YqwzZa/yEp7c3FLZ1HzlUl/A8bneEc3Zmwseye2O6NGtdt5nDqx85JbPHo2/ZKtkM7LPuTev29L9vffAnL5bZTle/dur6U9hliYuIjtF5Izh+JkAXg7gYQCHAOwdqu0FcMtkDwkhhExim6POLgAHReQkrE/6N6vqrSJyF4CbReQaAI8CeM0G+kkIIaSAOYmr6pcAXFQo/y6ASzbCKUIIIT662DsF8L03Ma034t0/wqtd5uc9fnv23bBytUu2vPtK1/ppaYHW/s2WBu/RBK09KqzxyMfEk/tdiqXmQ14/j7UVW8uOpfVH9zUv2fdoxdb6jXd9pHTO2uenFNvU9Zu0PeDfi8W6Hz1zjvf6ee6zFtw7hRBClhRO4oQQ0jFdyCneny9eWSHH+rk05eet14dF/PSNplJGfza2fE3x/mwttan14ZU+rDRFz3Vs+RCRIDzXySvReFMzS3Wmyo7RtL607tStl2vxlmKq3e/erZprtNpHHvPPfV4ElFMIIWRJ4SROCCEdw0mcEEI6pgtNfCSq0dba5Oci6YSt9DoPLe1xjl2grZ9H9NOSP96UvbztnLTClv3xszc11FovsGjFW8ObXpq3Scsj6ZItP3Kbpc+WfUszr6UYRtabanjXJmpxe9YbSufSfkasbQTyMVoE1MQJIWRJ4SROCCEdw0mcEEI6phtNPPKYc0Tz9Gq7QOzVaLXtNGuxtPRGTy6y99HovL98HDy53pF8Zk+sEb3Ssl+qF9VeR/LrWKpT8z8t9+S+5/62fPKc8+apR7XlFpF7Z+4zCaV+a3j09Np33coZ38jc8BRq4oQQsqRwEieEkI7hJE4IIR3ThSbu2RcjJbL3xWgzsm2pBytPNqdVz7v9Z8v3FM+eGjUtuaaZztkzo+a7tbVoGn9Ly83redZX5uQYe/HmXU/xr2Y/7yNvVzrX8qPmi/d+nrLfiXc9pBafte9Nrd+Sv5E1nzl6OTVxQghZUjiJE0JIx3ASJ4SQjulCEx+JannpuZpGndqN5MG28sJTvK9qmqqxenJ8Pftot8an5J9XW07btnJ/Pesa0fM1P4Bpr7YrlZf6i9ynU++5ks+lNi1b3n1IpvQb1datOK3vUWStJfenVM/yy3NdWnajUBMnhJAlhZM4IYR0DCdxQgjpmC408Wjuael8qX1q39qHxJMbPPd8xEYrx9Vj3zOWpXLg+P1gvLphKd+9tcdHLfa03KNxttY8Wte+1IcVZ+6nZ78d7znPPj6pjTzO2vnI2oNVv+aHZ/2l1re1XlPzycox9+6nU6L1XtscauKEEEKqcBInhJCO4SROCCEd060mnhPRI9Nzlr452my9P9Czn0IeS8mXGh69sNRX2l9k7Dy50953KLZ89miqNVqab0RzzdtHtOKc2rMDQHk/8qhe29LXS0TWeqwc8Zrvpfpz3ivqvT9KYxV99qDmQ6t86v07VxenJk4IIUuKOYmLyHNE5DMi8pCIPCAibxzKd4jI7SJyZPi7fePdJYQQkuL5n/iTAP5SVX8JwIsB/JmIXAhgP4A7VPV8AHcMnwkhhGwiYU1cRG4B8J7h38WqelREdgG4U1UvaLVddJ546Xz6eWTK3glT8ofzc1Yu6ZT3PuZYml2pX6sv73hPaZv7Xst9jux1M/VzyqL3Zc/rlOx6sa6p1b9nLAH/GlLed8mf3AcPnucBrPhqfXrz/b1rZd41lEXkiAML1MRFZDeAiwDcDeBsVT0KAMPfs2b6SQghJIh7EheRZwH4GIA3qeoTVv2k3T4RWRWR1bW1tSk+EkIIqeCSU0TkZAC3ArhNVd81lB3GJsopkUdkrVQsb9qd52dsydf8XFSeiPrVSrPKfYukb5V8y+u0pAnr57H3ceopKZxz6pbSN61XenleOVbC+sk+RRprXZuSz610Ve/4RSQ27/bMKVFZKiLJ1GLxfkdaPkTuhRaz5BQREQDvB/DQOIEPHAKwdzjeC+CWWV4SQggJs81R5yUAXg/gPhH5wlD2VwCuB3CziFwD4FEAr9kQDwkhhFQxJ3FV/SwAqZy+ZLHuEEIIidDVY/cjXk2vhOeVTKX+I69yyut4U5XmPuZttfWkk0XTEqPXoKWptraS9dqNppaV7LR0+Ij+Xoo5b+dJh2ylvUW3abX89r5CzWO/Rm38Pdv21u7P6Pa2LVt5/VaMkfWAOfCxe0IIWVI4iRNCSMdwEieEkI7pQhMH2ttVpudHpmqh43Gp31o/NR/z86X2Ht+8vub9TNmi0+rDk2vuzdMv1fU+Hm35WOo3978Vq7UmYeUAT7kXouNZi6/UV+18JL/bamfl0lv9lGxZ9VKfSn7N1cg96witeLZEnjghhJCtCydxQgjpGE7ihBDSMd1o4iNTdbGSDUvnjORJp31F9TSP7x6NtdZ+av2Sf3k7KyaPNp/XBcqvM0vbtWxarzZrradYOnfrdWc1v6P3Y3TNx8qdzvFq3JadUh1vrrjXv7wsv44AXK+fG22UaK231HwsXe+W73OhJk4IIUsKJ3FCCOkYTuKEENIxXWjilsaclqV49862dOhauVe7a533aOgp1l4all+leEp497+wdMmRyJ4Y3v5L9b15vNG1j5rfpbLWWoj39V9WTN5XrtX8Ldnw7J3v/Z55rkPah6VhW/n9JfsRnd+jh7fWbmqvGSzZngI1cUIIWVI4iRNCSMdwEieEkI7pQhMH6vmvU46Bukab9zc3t7qk93l10VbO8xiDpcVaWqul/eVM0U0tvdubW2ytSbR0bs/6gSePuNZ/VItvxehdI5nSl2d/E6s9gKqNHK/v1nMF1jpTzaZnrcEah1LdiO49Vw8HqIkTQsjSwkmcEEI6hpM4IYR0TBea+JxcYm8OsDfn1KttArF9oFOm5M169/iO6vel+HNfc1tTc7S9ey9HtU+rXl7X8/7RsU1rT5bU7tTc8lIcJb+imm7NTimumv+tc7Wc6qnj3vr+Wdr/lDUG79pB3tYap6lQEyeEkCWFkzghhHQMJ3FCCOmYrjTxGl69MbUVzS21/Jj6vkWPvtfyt+S3FUs0Xzqaa9vyyetvqX7Lf8C3Z7hXF43GWuozJ5pPX+srMq7edRPP+lEkT752b5X6i+rM1jqE1b7kb6k8+q6AKes7XqiJE0LIksJJnBBCOqYrOSVNLwLqj53n7UZaP/GsbScjj9jWbOWxTIm1FZvH3xqen32en9E1e96f2Sme17q1UhZb/Xnlqda1zP0p2Sidj6YPpvWislg0HW5q2qP3Ho/IWGN5JN1vSjpmbVuMFM92A4v4jtWgnEIIIUsKJ3FCCOkYcxIXkQ+IyDERuT8p2yEit4vIkeHv9o11kxBCSAlTExeRlwL4AYB/VNXnDWXvBPC4ql4vIvsBbFfVa63OFpli2EqdKmloNV1tbA+0t6etaa9ebbt0XIrPm6ZW89NTLzqepX5q4xdJO2sR0dBb17sWd/7qOEtfjT6G3oqj5Hspjrx9HltUi/Wkw3nWKzxrFED5+zT2N563tPGIZt3SzVtjW4trjh3vfe5lliauqv8B4PGs+AoAB4fjgwCunOMgIYSQaUzVxM9W1aMAMPw9q1ZRRPaJyKqIrK6trU3sjhBCSIkNX9hU1RtVdUVVV3bu3LnR3RFCyP8rXHniIrIbwK2JJn4YwMWqelREdgG4U1UvsOwsKk+8dJzWG4nkJI9EXlU2t6/Uft5HJD/Zm/NsHZdstLTAWp81/djKn5+yva83F9lad4iMfSQvPNfeSzY96yueXOmUmi5t6dWtOFoatzc279h5NOja/ePR2Wt+pv5Y+fR5v94xmcJG5IkfArB3ON4L4JaJdgghhMzAk2L4YQB3AbhARB4TkWsAXA/gUhE5AuDS4TMhhJBNZptVQVVfWzl1yYJ9IYQQEqS7Jza9GvWDR5/AVTfchZvesOepNrkuleqV4/F9By47rq/xOP9cOq75WztO4xj7uOkNe4r641U33IXnH7gNV91w19PajWV5bGl8qY20Xc3/1L/cn7T/UtvaGJa0zdTWfQcuO26cczujjbReaqcUUz42lu5e01VLumjt2qf30/gvj2X0N79Opb7y8+P1Ha9P7tN9By4rxnnhrtPNe7zUfzrupTjHv2n71EYrJ7u1rjH6mtps+Vmzk9rI/UxjH2Np2c3Pp/dLPo+UvoOLprtJnBBCyE/hJE4IIR3DSZwQQjqmq/3Ec1q5zZ72Vv1IuZXj7NkLwpvH6n0lVd5uLPf2WWqfx1Mag5rNVk65Jw/YauPdf7zlUysn2LJZIppz3soDr+UmR/ZTb/VZO+fJofbk+EfiL+HN5a7dS61xaH0nrXz1ku9W7n8U7idOCCFLCidxQgjpGE7ihBDSMV1o4oCtJdd0Q6D+fspoWWlPDO97Pi3tOXrOo/NH9N6SrjlnbwmgvFd26bx3/Lxjaemknjqtvmp6sFc/9dzHI16fanVqTF1HqmngES07t9WyW2pTsl3rP3rfWmscaV/ee2QRUBMnhJAlhZM4IYR0DCdxQgjpmG418bS8hZUzWqof0WEtfS31HcDk/PC0zWgHaOvKnvzlOXm7pXb5ONQ0Tq+2bxFda4ho3dZxyX5e5vG/pjFP0Xtb51u2I88DRPybs6+2Z1ysZw9q96LnPbwtv6xYFq2PUxMnhJAlhZM4IYR0DCdxQgjpmK40cWBa7nJUl7L094jOHs1tbxHJD7bGIbXpyc/N61tlNd9Tf2tj02rj6cObB1zTTCP6uFfz92jPkTz4yDtPS/Hl9a172lrr8H4vx7ZAeV93S6OOfo9zmzmt98G2fPRcx1psU6AmTgghSwoncUII6Zhu5RTPT2tPSiFgp8GV/Ci1j/ykTO1bP+ktPKmMY1mtb+tx5qgEMcV26ps3vc4rW+TXxtoSoOaHNZ4tP6P3hEeyqKWr1vpupbe25KQpsqAVS0Tys2JqjXV0G1lr7Dw+58yRUgDKKYQQsrRwEieEkI7hJE4IIR3ThSY+V1Ob06ZkY/Sp1dZKnxvrAO3tbD1+RbfvrNmekpZp6fyl+K0taD0aeUSHjejRU1LHLLwplGkcnjoef2tbO7TsW754bUTXPqyUzch3wbN+UIsrouXXYq2VTYWaOCGELCmcxAkhpGM4iRNCSMd0o4kDdv6xt/7UnM7II8hT8thbW9NaObFztMSWzufJpy5hbQc65dp584dbtiL50ZZ/U9ZUUjzrJpZ9b+7zlNz8uWtOLb8819Va1wBir/jL7efnS2PjfW7D+7rGqVATJ4SQJYWTOCGEdMysSVxEXikih0XkERHZvyinCCGE+JisiYvISQC+DOBSAI8BuAfAa1X1wVqbRe2dYu3rUMKrd0VzpVv9R/Zs8ObCe/R3Sz+s+VaLz2O71EceX/o5gjfuHG9es1cTrdny5MrX4mr52eqjVMfCo4NHfYxo7C27Vt8efd66VjnRdRZrDcuKYQ4bpYm/CMAjqvpVVf1fAB8BcMUMe4QQQoLMmcTPAfDN5PNjQ9nTEJF9IrIqIqtra2szuiOEEJIzZxKXQtlx2oyq3qiqK6q6snPnzhndEUIIyZmjie8BcEBVLxs+XwcAqvr2Wps5mjghhPx/ZaM08XsAnC8i54nIKQCuBnBohj1CCCFBtk1tqKpPisifA7gNwEkAPqCqDyzMM0IIISaTJ3EAUNVPAvjkgnwhhBAShE9sEkJIx3ASJ4SQjuEkTgghHcNJnBBCOoaTOCGEdAwncUII6RhO4oQQ0jGcxAkhpGM29R2bIrIG4BsTm58J4DsLdOdEsiyxLEscAGPZqjCWdX5RVYs7CG7qJD4HEVmtbQDTG8sSy7LEATCWrQpjsaGcQgghHcNJnBBCOqanSfzGE+3AAlmWWJYlDoCxbFUYi0E3mjghhJDj6el/4oQQQjI4iRNCSMds+UlcRF4pIodF5BER2X+i/SkhIs8Rkc+IyEMi8oCIvHEo3yEit4vIkeHv9qTNdUNMh0XksqT8V0TkvuHc34lI6YXUGx3PSSLyeRG5tfM4zhCRj4rIw8O12dNxLH8x3Fv3i8iHReTUXmIRkQ+IyDERuT8pW5jvIvIMEblpKL9bRHZvcix/M9xjXxKRT4jIGZsai6pu2X9Yf+3bVwA8F8ApAL4I4MIT7VfBz10AXjgcPxvAlwFcCOCdAPYP5fsBvGM4vnCI5RkAzhtiPGk49zkAewAIgH8D8JsnIJ43A/gQgFuHz73GcRDAHw/HpwA4o8dYAJwD4GsAnjl8vhnAH/YSC4CXAnghgPuTsoX5DuBPAbx3OL4awE2bHMsrAGwbjt+x2bFs6pdqwoDtAXBb8vk6ANedaL8cft8C4FIAhwHsGsp2AThcigPr7yndM9R5OCl/LYAbNtn3cwHcAeBl+Okk3mMcp2N94pOsvMdYzgHwTQA7sP5KxVuHiaObWADszia+hfk+1hmOt2H9qUjZrFiyc78N4IObGctWl1PGm3fksaFsyzL8/LkIwN0AzlbVowAw/D1rqFaL65zhOC/fTN4N4C0AfpKU9RjHcwGsAfiHQRp6n4ichg5jUdX/BPC3AB4FcBTA91T1U+gwloRF+v5UG1V9EsD3APz8hnne5o+w/j/rp/k1sCGxbPVJvKTXbdmcSBF5FoCPAXiTqj7Rqloo00b5piAilwM4pqr3epsUyk54HAPbsP6z9+9V9SIAP8T6z/YaWzaWQS++Aus/yX8BwGki8rpWk0LZlojFwRTft0RcIvJWAE8C+OBYVKi28Fi2+iT+GIDnJJ/PBfCtE+RLExE5GesT+AdV9eND8bdFZNdwfheAY0N5La7HhuO8fLN4CYBXi8jXAXwEwMtE5J/RXxwYfHhMVe8ePn8U65N6j7G8HMDXVHVNVX8E4OMAfg19xjKySN+faiMi2wD8HIDHN8zzAiKyF8DlAH5fBy0EmxTLVp/E7wFwvoicJyKnYF3oP3SCfTqOYWX5/QAeUtV3JacOAdg7HO/FulY+ll89rESfB+B8AJ8bflZ+X0RePNj8g6TNhqOq16nquaq6G+tj/WlVfV1vcQyx/BeAb4rIBUPRJQAeRIexYF1GebGI/OzgwyUAHkKfsYws0vfU1u9i/b7dzF+wrwRwLYBXq+r/JKc2J5bNWNSYuYjwKqxne3wFwFtPtD8VH38d6z95vgTgC8O/V2Fdy7oDwJHh746kzVuHmA4jyRAAsALg/uHce7CBCzRGTBfjpwubXcYB4AUAVofr8q8Atnccy9sAPDz48U9Yz3joIhYAH8a6lv8jrP9P85pF+g7gVAD/AuARrGd9PHeTY3kE6zr2+N1/72bGwsfuCSGkY7a6nEIIIaQBJ3FCCOkYTuKEENIxnMQJIaRjOIkTQkjHcBInhJCO4SROCCEd838NIB+f7OmAywAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"#from models import Model, SynapticConductanceModel\n",
"from brian2 import *\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"\n",
"model = SynapticConductanceModel(resistance=Model.LOW, # Model.LOW, Model.HIGH\n",
" noise=Model.HIGH, n=50, # Model.OFF, Model.LOW, Model,HIGH\n",
" offset=SynapticConductanceModel.ACTIVE)\n",
"\n",
"\n",
"# model = SineModel(resistance=Model.LOW, # Model.LOW, Model.HIGH\n",
"# noise=Model.OFF, n=1, # Model.OFF, Model.LOW, Model,HIGH\n",
"# )\n",
"\n",
"\n",
"model.f = 400 * Hz\n",
"model.set_stimulus_current(400 * nS) # current should be scaled by 100x for Active Offset so (500nS is actually 5nS)\n",
"model._set_variable(\"i_injected\", 65 * pA)\n",
"\n",
"\n",
"model.run(12*second)\n",
"\n",
"spike_times = [s/ms for s in model.spike_train.values()]\n",
"\n",
"plt.eventplot(spike_times)\n"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "7f22e21e",
"metadata": {},
"outputs": [],
"source": [
"#print(spike_times)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "20e3191c",
"metadata": {},
"outputs": [],
"source": [
"import scipy.io\n",
"import numpy as np\n",
"\n",
"\n",
"file_path = 'data0.mat'\n",
"scipy.io.savemat(file_path, {'data0': spike_times[0]})\n",
"\n",
"file_path = 'data1.mat'\n",
"scipy.io.savemat(file_path, {'data1': spike_times[1]})\n",
"\n",
"file_path = 'data2.mat'\n",
"scipy.io.savemat(file_path, {'data2': spike_times[2]})\n",
"\n",
"file_path = 'data3.mat'\n",
"scipy.io.savemat(file_path, {'data3': spike_times[3]})\n",
"\n",
"file_path = 'data4.mat'\n",
"scipy.io.savemat(file_path, {'data4': spike_times[4]})\n",
"\n",
"file_path = 'data5.mat'\n",
"scipy.io.savemat(file_path, {'data5': spike_times[5]})\n",
"\n",
"file_path = 'data6.mat'\n",
"scipy.io.savemat(file_path, {'data6': spike_times[6]})\n",
"\n",
"file_path = 'data7.mat'\n",
"scipy.io.savemat(file_path, {'data7': spike_times[7]})\n",
"\n",
"file_path = 'data8.mat'\n",
"scipy.io.savemat(file_path, {'data8': spike_times[8]})\n",
"\n",
"file_path = 'data9.mat'\n",
"scipy.io.savemat(file_path, {'data9': spike_times[9]})\n",
"\n",
"file_path = 'data10.mat'\n",
"scipy.io.savemat(file_path, {'data10': spike_times[10]})\n",
"\n",
"file_path = 'data11.mat'\n",
"scipy.io.savemat(file_path, {'data11': spike_times[11]})\n",
"\n",
"file_path = 'data12.mat'\n",
"scipy.io.savemat(file_path, {'data12': spike_times[12]})\n",
"\n",
"file_path = 'data13.mat'\n",
"scipy.io.savemat(file_path, {'data13': spike_times[13]})\n",
"\n",
"file_path = 'data14.mat'\n",
"scipy.io.savemat(file_path, {'data14': spike_times[14]})\n",
"\n",
"file_path = 'data15.mat'\n",
"scipy.io.savemat(file_path, {'data15': spike_times[15]})\n",
"\n",
"file_path = 'data16.mat'\n",
"scipy.io.savemat(file_path, {'data16': spike_times[16]})\n",
"\n",
"file_path = 'data17.mat'\n",
"scipy.io.savemat(file_path, {'data17': spike_times[17]})\n",
"\n",
"file_path = 'data18.mat'\n",
"scipy.io.savemat(file_path, {'data18': spike_times[18]})\n",
"\n",
"file_path = 'data19.mat'\n",
"scipy.io.savemat(file_path, {'data19': spike_times[19]})\n",
"\n",
"file_path = 'data20.mat'\n",
"scipy.io.savemat(file_path, {'data20': spike_times[20]})\n",
"\n",
"file_path = 'data21.mat'\n",
"scipy.io.savemat(file_path, {'data21': spike_times[21]})\n",
"\n",
"file_path = 'data22.mat'\n",
"scipy.io.savemat(file_path, {'data22': spike_times[22]})\n",
"\n",
"file_path = 'data23.mat'\n",
"scipy.io.savemat(file_path, {'data23': spike_times[23]})\n",
"\n",
"file_path = 'data24.mat'\n",
"scipy.io.savemat(file_path, {'data24': spike_times[24]})\n",
"\n",
"file_path = 'data25.mat'\n",
"scipy.io.savemat(file_path, {'data25': spike_times[25]})\n",
"\n",
"file_path = 'data26.mat'\n",
"scipy.io.savemat(file_path, {'data26': spike_times[26]})\n",
"\n",
"file_path = 'data27.mat'\n",
"scipy.io.savemat(file_path, {'data27': spike_times[27]})\n",
"\n",
"file_path = 'data28.mat'\n",
"scipy.io.savemat(file_path, {'data28': spike_times[28]})\n",
"\n",
"file_path = 'data29.mat'\n",
"scipy.io.savemat(file_path, {'data29': spike_times[29]})\n",
"\n",
"file_path = 'data30.mat'\n",
"scipy.io.savemat(file_path, {'data30': spike_times[30]})\n",
"\n",
"file_path = 'data31.mat'\n",
"scipy.io.savemat(file_path, {'data31': spike_times[31]})\n",
"\n",
"file_path = 'data32.mat'\n",
"scipy.io.savemat(file_path, {'data32': spike_times[32]})\n",
"\n",
"file_path = 'data33.mat'\n",
"scipy.io.savemat(file_path, {'data33': spike_times[33]})\n",
"\n",
"file_path = 'data34.mat'\n",
"scipy.io.savemat(file_path, {'data34': spike_times[34]})\n",
"\n",
"file_path = 'data35.mat'\n",
"scipy.io.savemat(file_path, {'data35': spike_times[35]})\n",
"\n",
"file_path = 'data36.mat'\n",
"scipy.io.savemat(file_path, {'data36': spike_times[36]})\n",
"\n",
"file_path = 'data37.mat'\n",
"scipy.io.savemat(file_path, {'data37': spike_times[37]})\n",
"\n",
"file_path = 'data38.mat'\n",
"scipy.io.savemat(file_path, {'data38': spike_times[38]})\n",
"\n",
"file_path = 'data39.mat'\n",
"scipy.io.savemat(file_path, {'data39': spike_times[39]})\n",
"\n",
"file_path = 'data40.mat'\n",
"scipy.io.savemat(file_path, {'data40': spike_times[40]})\n",
"\n",
"file_path = 'data41.mat'\n",
"scipy.io.savemat(file_path, {'data41': spike_times[41]})\n",
"\n",
"file_path = 'data42.mat'\n",
"scipy.io.savemat(file_path, {'data42': spike_times[42]})\n",
"\n",
"file_path = 'data43.mat'\n",
"scipy.io.savemat(file_path, {'data43': spike_times[43]})\n",
"\n",
"file_path = 'data44.mat'\n",
"scipy.io.savemat(file_path, {'data44': spike_times[44]})\n",
"\n",
"file_path = 'data45.mat'\n",
"scipy.io.savemat(file_path, {'data45': spike_times[45]})\n",
"\n",
"file_path = 'data46.mat'\n",
"scipy.io.savemat(file_path, {'data46': spike_times[46]})\n",
"\n",
"file_path = 'data47.mat'\n",
"scipy.io.savemat(file_path, {'data47': spike_times[47]})\n",
"\n",
"file_path = 'data48.mat'\n",
"scipy.io.savemat(file_path, {'data48': spike_times[48]})\n",
"\n",
"file_path = 'data49.mat'\n",
"scipy.io.savemat(file_path, {'data49': spike_times[49]})\n",
"\n",
"\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d353eff4",
"metadata": {},
"outputs": [],
"source": [
"\n"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.8"
}
},
"nbformat": 4,
"nbformat_minor": 5
}