{
"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.57797 s (4%) simulated in 10s, estimated 3m 18s remaining.\n",
"1.15458 s (9%) simulated in 20s, estimated 3m 8s remaining.\n",
"1.7348 s (14%) simulated in 30s, estimated 2m 58s remaining.\n",
"2.30077 s (19%) simulated in 40s, estimated 2m 49s remaining.\n",
"2.8383 s (23%) simulated in 50s, estimated 2m 41s remaining.\n",
"3.37323 s (28%) simulated in 1m 0s, estimated 2m 33s remaining.\n",
"3.80617 s (31%) simulated in 1m 10s, estimated 2m 31s remaining.\n",
"4.3266 s (36%) simulated in 1m 20s, estimated 2m 22s remaining.\n",
"4.81076 s (40%) simulated in 1m 30s, estimated 2m 15s remaining.\n",
"5.34558 s (44%) simulated in 1m 40s, estimated 2m 5s remaining.\n",
"5.88416 s (49%) simulated in 1m 50s, estimated 1m 54s remaining.\n",
"6.42113 s (53%) simulated in 2m 0s, estimated 1m 44s remaining.\n",
"6.95157 s (57%) simulated in 2m 10s, estimated 1m 34s remaining.\n",
"7.47731 s (62%) simulated in 2m 20s, estimated 1m 25s remaining.\n",
"8.01432 s (66%) simulated in 2m 30s, estimated 1m 15s remaining.\n",
"8.55137 s (71%) simulated in 2m 40s, estimated 1m 5s remaining.\n",
"9.08927 s (75%) simulated in 2m 50s, estimated 54s remaining.\n",
"9.60162 s (80%) simulated in 3m 0s, estimated 45s remaining.\n",
"10.10329 s (84%) simulated in 3m 10s, estimated 36s remaining.\n",
"10.61838 s (88%) simulated in 3m 20s, estimated 26s remaining.\n",
"11.1533 s (92%) simulated in 3m 30s, estimated 16s remaining.\n",
"11.65781 s (97%) simulated in 3m 40s, estimated 6s remaining.\n",
"12. s (100%) simulated in 3m 46s\n"
]
},
{
"data": {
"text/plain": [
"[<matplotlib.collections.EventCollection at 0x231816a6af0>,\n",
" <matplotlib.collections.EventCollection at 0x231816b8df0>,\n",
" <matplotlib.collections.EventCollection at 0x231816c31c0>,\n",
" <matplotlib.collections.EventCollection at 0x231816ca520>,\n",
" <matplotlib.collections.EventCollection at 0x231816d27c0>,\n",
" <matplotlib.collections.EventCollection at 0x231816d99d0>,\n",
" <matplotlib.collections.EventCollection at 0x231816e1c10>,\n",
" <matplotlib.collections.EventCollection at 0x231816e7fd0>,\n",
" <matplotlib.collections.EventCollection at 0x231816f3370>,\n",
" <matplotlib.collections.EventCollection at 0x231816fc580>,\n",
" <matplotlib.collections.EventCollection at 0x231817037f0>,\n",
" <matplotlib.collections.EventCollection at 0x231817099a0>,\n",
" <matplotlib.collections.EventCollection at 0x23181711d60>,\n",
" <matplotlib.collections.EventCollection at 0x2318171e190>,\n",
" <matplotlib.collections.EventCollection at 0x231817254f0>,\n",
" <matplotlib.collections.EventCollection at 0x2318172c4f0>,\n",
" <matplotlib.collections.EventCollection at 0x23181732940>,\n",
" <matplotlib.collections.EventCollection at 0x2318173cbb0>,\n",
" <matplotlib.collections.EventCollection at 0x23181743e50>,\n",
" <matplotlib.collections.EventCollection at 0x2318174d220>,\n",
" <matplotlib.collections.EventCollection at 0x23181755640>,\n",
" <matplotlib.collections.EventCollection at 0x2318175e9a0>,\n",
" <matplotlib.collections.EventCollection at 0x23181765be0>,\n",
" <matplotlib.collections.EventCollection at 0x2318176ce50>,\n",
" <matplotlib.collections.EventCollection at 0x23181778100>,\n",
" <matplotlib.collections.EventCollection at 0x231817803d0>,\n",
" <matplotlib.collections.EventCollection at 0x23181787580>,\n",
" <matplotlib.collections.EventCollection at 0x2318178e880>,\n",
" <matplotlib.collections.EventCollection at 0x23181794a00>,\n",
" <matplotlib.collections.EventCollection at 0x2318179dcd0>,\n",
" <matplotlib.collections.EventCollection at 0x231817a4fa0>,\n",
" <matplotlib.collections.EventCollection at 0x231817b13a0>,\n",
" <matplotlib.collections.EventCollection at 0x231817b6790>,\n",
" <matplotlib.collections.EventCollection at 0x231817bfa90>,\n",
" <matplotlib.collections.EventCollection at 0x231817c7dc0>,\n",
" <matplotlib.collections.EventCollection at 0x231817d2160>,\n",
" <matplotlib.collections.EventCollection at 0x231817d93d0>,\n",
" <matplotlib.collections.EventCollection at 0x231817e26d0>,\n",
" <matplotlib.collections.EventCollection at 0x231817e9a90>,\n",
" <matplotlib.collections.EventCollection at 0x231817f0b80>,\n",
" <matplotlib.collections.EventCollection at 0x231817f7e50>,\n",
" <matplotlib.collections.EventCollection at 0x231817ffe50>,\n",
" <matplotlib.collections.EventCollection at 0x2318180b1f0>,\n",
" <matplotlib.collections.EventCollection at 0x23181812460>,\n",
" <matplotlib.collections.EventCollection at 0x2318181a610>,\n",
" <matplotlib.collections.EventCollection at 0x23181820a00>,\n",
" <matplotlib.collections.EventCollection at 0x23181827df0>,\n",
" <matplotlib.collections.EventCollection at 0x23181834220>,\n",
" <matplotlib.collections.EventCollection at 0x2318183b670>,\n",
" <matplotlib.collections.EventCollection at 0x23181843a00>]"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD4CAYAAAAaT9YAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAgLElEQVR4nO2dbaxmV1XH/8uZlmKhdsbemYwtcUrSNDYYKV4JFUMqQwGxodWILRG8xuqYqAmICZ3KB+ETBQ0hhkQ6AcyovLTy4kwmaJkUGkMyKb0T3vo2THkrI1fmQg0F/CCF5YfnnHK6u1/W2uc8z9z9+P8lk3ue/brWOedZmee/19lHVBWEEELa5KfOtgGEEELqYRAnhJCGYRAnhJCGYRAnhJCGYRAnhJCG2b7IyS666CLdu3fvIqckhJDmOXHixLdVdSVWt9AgvnfvXqyvry9ySkIIaR4R+XqqjnIKIYQ0DIM4IYQ0DIM4IYQ0DIM4IYQ0DIM4IYQ0jCk7RUS+BuB7AH4E4HFVXRWRnQBuB7AXwNcA/K6q/vd8zCSEEBLD8z/xX1fV56rqavf5AIC7VPUyAHd1nwkhhCyQMXLKdQAOdceHAFw/2hpCCCEurEFcAXxCRE6IyP6ubLeqbgBA93dXrKOI7BeRdRFZ39zcHG8xIYSQJ7A+sflCVf2miOwCcExEHrJOoKoHARwEgNXVVb6BghBCJsQUxFX1m93fMyLyMQDPB/AtEdmjqhsisgfAmXkZecNtx584fmDjMVyx54InHd/+J1c9pe2wLDZWrH44z7BNqk+pfDhOauzYWKXj3LgpUr7k7Eq19fibmyc2Xun6WefNlaVs8J7rUp9hu5jdpWti/Wwpy41xw23Hn/S9yvUZfudS/VL2xMjZHY5prc/NVYoZ1vu05IOnfixFOUVEzheRZ/bHAF4K4D4ARwCsdc3WAByei4WEEEKSWP4nvhvAx0Skb/8BVf13EbkXwB0ichOARwC8an5mEkIIiVEM4qr6FQC/FCn/DoB98zCKEEKIDVnk2+5XV1e1Zivaku47RnezaMO5+pKeFvYNNcQUFr07Z0PKr9w8pT5WHT/WJjXGsF1Jly3p96HeaVn/sGieMT9Cau9Bq0Zbe47H6Mc9sfMaauNWG63XLCQ2nufcxfyKjZvCs+ZQo6uXEJETg2d0ngQfuyeEkIZhECeEkIZhECeEkIZpQhMH4popMNPrAERzVft+fVlqDGv9cLxUuVULHM7pmSdGKQe4ZNeYvPacfd78Zo9+Hyvz6o1T6PhhfW3OfG7smnvIo/mn9OWYj57rVcKTp1/qM6wPbY/VxfB+J6x1U0BNnBBClhQGcUIIaRgGcUIIaZgmNHGP5uvRpi1zheUpLDnhNXp3zCbrXg3D8Wu1zVJee8omT171GE13yn0patYhwrapfPfQVqu2n1rTseZsp2yJ+VDy26pTx2wflpXWpcLjWH1p7th8OTtLawJDLHs25casgZo4IYQsKQzihBDSMAzihBDSME1p4j1WbdSTI13Sdr37QuR8KOmaMZ+H/XosbUu63ZCY7urJUc718+RKe3KcU/6nxo71s641lPp48+0tzxx48qitud6pOoue7tHRrd8Xz32YapOzwzK35VkBy/UKP0+1ZkNNnBBClhQGcUIIaRgGcUIIaZgmNfEQi3bo3bcipymn5g3bWfT5cN6YrTX6XWkO4KkatvWdpV4NcEwOsdVXj3Zc6p+zxarnpsaaUnuN9bf6F/avWTuxjJ/qH5vLos977rGSjxbN37OGYxmnFmrihBCypDCIE0JIwzQlp3h/OqdkAm9/b9tU/9zcKb/HPLrr/YnrTZ8bY4tXXqpN1fKkxKXsqvF9zL0Ss7V2jJQcENbn7K6RMHKU5JqUbbl5PXaHY3m3lgBsKbypshoopxBCyJLCIE4IIQ3DIE4IIQ3TnCbuSVfyPN6e0iQtW3ladV6PxuxJdfLo29Z+ubEsjFmrSFGj7Vv101h/IP0oezhOzp5Y25oUPUt6oLW+tN5S0qkturPnvvGuf9To9VabLenFPVYbx0JNnBBClhQGcUIIaRgGcUIIaZjmNPHh5x6P9pnSY3PzeWwc2mQZM/VIccq3KfTHnM21efPW8XPlw37e/HiLLuo9jzltvJRfnLumIZb87tJ5qLn/PHZ51lZy/S32eNZIPHN71mFK97rl2QFq4oQQQoowiBNCSMOYg7iIbBORz4rI0e7zThE5JiKnur875mcmIYSQGGZNXETeAGAVwAWqeq2IvB3Ao6p6q4gcALBDVW/OjTFlnrhF86rNea7Rv3L9a9rW7lthzTeeUvsu+TYc1+KbZUyPvh7W19wDuXqrdpp7RiE3pnXuGtstePLqc7ZY10hCvHaW7Cj1z62HeNYoxux7FDJaExeRSwD8JoD3DIqvA3CoOz4E4PoRNhJCCKnAKqe8E8AbAfx4ULZbVTcAoPu7K9ZRRPaLyLqIrG9ubo6xlRBCSEAxiIvItQDOqOqJmglU9aCqrqrq6srKSs0QhBBCEhQ1cRF5K4DXAngcwHkALgDwUQC/AuBqVd0QkT0A7lbVy3NjzVMTt+71EKuvqRtiec1Tyh6vvm7Vwi125+wC8n5NbXfOjli9dY5YPndJx7X6UaPte55nsI5j0f29awepdjkfrJp/OHbtukWMsWssuTpvTrrHjhKjNHFVvUVVL1HVvQBuBPBJVX0NgCMA1rpmawAOT2ItIYQQM2PyxG8FcI2InAJwTfeZEELIAtnuaayqdwO4uzv+DoB905tECCHESlN7p6TI7QeR0zJz/S26aW2Ot2Uer4ZcoiZHPGePp01Js875acnrBWzrEmNzvlNzW+bIzVNqn2Je6yk1Y1jmj9X1eO6JWB9LP48+7V1T8trkhXunEELIksIgTgghDcMgTgghDdOUJu7VNHPjpfY1qJmrx7NXhTVfOIVVv7Vo3qH9tdqi97yG9nvef2rR2YH4OzKtY3n0YIvGG1KjQeeuzQMbjwFAdr9sy5pIypawzLt+4Fl7iuX3e7Hmp+fywL1rLVPmvA+hJk4IIUsKgzghhDQMgzghhDRME5o4YNvDoien643R+UoaoEdPH34e2p3TAq1552FZbJy+PrdXslVnD8cdtrNq0rWacm3evNUW67tYrTndNfnhNfdurD5nf8kWTxvPMwkpP2rPe836UGnfb8/9N2RL7SdOCCFka8IgTgghDdOEnFLzU6unpo1FUhn2Lf38Dan9eTXVI7xe6WVMiqc1nTFGbXpX6Wd2yoYxskKtpGAhdh6mtKnULywf2mGV96zXwzKepZ/13AzHSNWl5vGmZ46BcgohhCwpDOKEENIwDOKEENIwTWniPSXtMvb4dulVYx79sCYVzzp/Cuujwql1A8tcNalp3muTsrFGdw/nHKOXp+YbtvGsfeRSK606tOVaW8cN+1ipSYkczmX5rpTWaGL3ybA8hUcP99wTNem4Y7VxauKEELKkMIgTQkjDMIgTQkjDuN6x2QqxHPJQUwu1qpTuOPx8w23Hq7Qtb55uWH7DbcefZEdOT7famNIZLXaGx97c8PA4Nndu/HCMVJ/YOevPT2ydJOZ7TnfNrWP0fYbXI2ZDXx47HzF/a9dqwvs+NVfqexH6PCb3PHZtUuMPy2N9xuSDx8iNk7sXwmu7SPg/cUIIaRgGcUIIaRgGcUIIaZim8sQt+yZY930I8ei5Kdtq6mNzWnO5c3OVcuNLc1m3qLWeb0sud8yP0mfv/KnxSnn8qTGtudw5/2rs9vhnmdfbzrN2Yek/HKevt2zJbOkXq0vZG7YJ7a0992NhnjghhCwpDOKEENIwDOKEENIwzWniqXKPbh72B56qAVt1MYuWbt2rITZvycecfj4c24J1P4hc33BOz3nx2muxZ8x5C+3x5Plb23vvVU+9xa+Shu+9rh4/Yp9jc5fWGzzxIUfu2QnLNfSsy3ihJk4IIUsKgzghhDRMMYiLyHki8hkR+byI3C8ib+nKd4rIMRE51f3dMX9zCSGEDClq4iIiAM5X1e+LyDkAPg3gdQB+G8CjqnqriBwAsENVb86NNdV+4kA+jznVZ4hFK7XkSOd0tHCu1LtBvbntVo3RogmGdaW1gZpzlvLVk2Ndo8/X5v/W3kdj7LaMAaTft1ka13tvWbVra/50zrZhe+/6l9Wf3BhT5fp7NXsPozRxnfH97uM53T8FcB2AQ135IQDXV1tICCGkCpMmLiLbRORzAM4AOKaq9wDYraobAND93ZXou19E1kVkfXNzcyKzCSGEAMYgrqo/UtXnArgEwPNF5DnWCVT1oKququrqyspKpZmEEEJiuPPEReSvAfwAwB8DuFpVN0RkD4C7VfXyXN8p37HpzW+N9U3tqZCac1g+Rgu12O7V/L06e03ea8oH67nLjW/RYD05xWGdZS+YEmPy80s6fa5P7fgxanXb0n1vzSf33MOe+S1MqfPn+ub2fKlllCYuIisicmF3/HQALwHwEIAjANa6ZmsADldbSAghpArLm332ADgkItswC/p3qOpRETkO4A4RuQnAIwBeNUc7CSGERCgGcVX9AoArI+XfAbBvHkYRQgix0dzeKTkt1tKnL4u1t+ZED+uGY3pybz0aojc326OlDv239rHm1Zby+0tY8/VzfpQ03Nix5T5IjTUGq4/esSw58lPYbF27qVmDiNXV2BO2iX3OtUudR++6kBfunUIIIUsKgzghhDRMc3JKrMz6ky720zL1E6nGjpRdPZ5Hmj1pft5UKI+8YPXdmqIW2lc7ztgUyVJ7y3W1Xpuxc3vtqbknYzbl2uT8sMgQ3u9pbo6UxDqP76vHhpCzlmJICCFk68IgTgghDcMgTgghDdOkJh57tLWUZlgax5JGCMQfMS+VTZ1u5F0DGB57NMOxaXzheBZfcufSSikVMCyz3AO1Gn7K35wNqXapuXNtrPbkNN7YHLF5Utcy1teLJ1XRssaTGt+61jLVeowVauKEELKkMIgTQkjDMIgTQkjDNKWJ91g0u1rdPFdX0sxr9Exr7mnKxtS4ublr+6TKcj7nGKMtp8Yazl+bf+6xreS3RVsu2WKxz/KawHDu3HWv3R7BUjacz/P9mwKLjg7kc72ta2vDPmOhJk4IIUsKgzghhDQMgzghhDRMU5r42PxP69g9Y/KGQ2q0SosWPayzjhuWl7Bop6kxPTm3Jbvmobtbrl1I7XMHsT5We633aazMun4xpQbtPa8WTd9zDw37lNayYm1zc85jb5QS1MQJIWRJYRAnhJCGYRAnhJCGaUoTB2yv9xqjoafmzY1R0ihrXt2Uy1n15mjXaJ3evOqwLJVrbNHvrc8FePR4734Y1r0xvDq61a4Qz7MHqTY5m0rUPleQG2dYltLES/vmWM9bbQ58OIZ1/Sk29xioiRNCyJLCIE4IIQ3DIE4IIQ3TlCbuzcOu2evCakfsc9jW8l7B3DwlLPtNWLVdT1lYnssX9+6NMWWufk0Oem0fi/0htfnq83hGIcZYPbpE7dqVN/89VR+2y/X3fjem3EsdoCZOCCFLC4M4IYQ0DIM4IYQ0THOaeOyzpY93HoveZc37tmqnHhti1OS9pmzJnfOaPSu8e2mEY6Xs9+wZksKSuzxsF9rlfZ7AslYT4sldT42dupbevcitWrb3e1qrSVttjLVL2ZKyNZVj7r3eXqiJE0LIklIM4iLyLBH5lIg8KCL3i8jruvKdInJMRE51f3fM31xCCCFDLP8TfxzAX6rqLwB4AYA/E5ErABwAcJeqXgbgru4zIYSQBeLWxEXkMIB3df+uVtUNEdkD4G5VvTzXd6p3bPZY9o8I25bmKOUzT5GL69WZS1pnrl1prwhPTnIJjzZb0mJz56RWw43VDe2aum+MKXRiyzMIVp3fMt8UWP3Oac6pvdprNG9PPn84XulcT6mF90ymiYvIXgBXArgHwG5V3QCA7u+ukXYSQghxYg7iIvIMAB8B8HpVfczRb7+IrIvI+ubmZo2NhBBCEpjkFBE5B8BRAHeq6ju6spPYInJKqp3lJ2tuzHn85Cq1y2FNMauZa56yRG4sy09eq70leSbV37s9Qumn/AMbs//jpMb0yGjetLtUnSfNtVTvSTEsyTc5m4flscfYpzo3KV/HlFn88jBKThERAfBeAA/2AbzjCIC17ngNwOFqCwkhhFSx3dDmhQBeC+CLIvK5ruyvANwK4A4RuQnAIwBeNRcLCSGEJCkGcVX9NABJVO+b1hxCCCEeLP8TbxpL+lTpkdlekxv2T+nAtZp5rCz8nNLUehtvuO14Mc2w/2zRRoflfZ+p0jlDe/p+lrkt86U09ZIe6tUvY/rnULsNz9eYezCcJ3Y8vFdy90Rsrtg1js0Z8yV1nNvKNhwzVx5b67B8B/vj0ivZcus6/by/+OY7s6mtqXM+ZaphDD52TwghDcMgTgghDcMgTgghDdPUVrQ9Vo2plHNszSWtyWv1lNdSytEdYs2XLY0fUpvnXcI7Tm49IZfba12TqMmBz7WJzenJ1U/ZmMuntubzW/LySzn1Xr9SfWrmCefMPQcQO2fDPqGNsXMR86c2ZqXgVrSEELKkMIgTQkjDMIgTQkjDNKGJA/F9GKyaq1VfG9b11OwnYZ0fsGlvVs0+Z5NVF47NY9m+tPa1VaW87Vi7MesVObwaeOn1bUP7Y/MM6y3abqpPCouGHtpnvWc8dsTsqsXzrMCwfep7av2+16wpedvloCZOCCFLCoM4IYQ0DIM4IYQ0TLN7p6T2augJ9y4Y7n0Q9huj6+X04JidYZuYHm7ViVPlOQ2uVi+P7QkRaz8kdY1Se0yk5g9tseq7MVtK16i0R85wrtS9FI4d3ndj9PvhufO28a4reNqFeL5bJZ17OFaNlt/X585vqs6aH1675jQF/J84IYQ0DIM4IYQ0DIM4IYQ0TBN54rV7UaT6W7Q+r8Zqzcet2ROi5FdNLnRPrXY3Jtc+1q7m/ZgePT/G1FqlNe++bzM2B9mT8+zJrS/1L41Xs9ZUk/Nv+U6MzStP4bnvprjPmCdOCCFLCoM4IYQ0DIM4IYQ0TJOauEXvsu4/EZunpJOn7LDqnCksenA4n9Wukg/WOSz6nkd/7duUNPGSDp6aN2Ss/bV2lI7H7guUwutvqqxWE6/xxXs9Pbp7OF9p/3XLOkFuvql0cWrihBCypDCIE0JIwzCIE0JIwzSliVs1xx6LJhkjtxdKP3/svXylOVO2x8awzjVsV6rzaJ6xdtYc3KEN1lxo6/Uonc9c2xKe9YiSj5b1iBQe7bgmV3qMVmv1y7KfkOfeK9nT49HIU5Ty62NzWdezaqEmTgghSwqDOCGENEyTcsqwrCeVqlXqb5U8rI8/ex779Y6dw/Oz0iOHWOeK2TsmVW/MNYt9TpWlfAqpOZ+e81crAUwhKZbOdwzr9yw1lrVfDM/3Lfc9y6UYliTVmF+1aaIWKKcQQsiSwiBOCCENUwziIvI+ETkjIvcNynaKyDEROdX93TFfMwkhhMQoauIi8iIA3wfwj6r6nK7s7QAeVdVbReQAgB2qenNpslpNHPBryzWPy3r02dw4sTkt/XKPY9fO4XnsOjZmzfmwpk/G7KxJ88u1tcxlve7e6x3rH9Ngw/ZhXemcTJGGmrM7VZbrm5tjSg3c03dMm5QdIWPSN3OM0sRV9T8APBoUXwfgUHd8CMD1YwwkhBBSR60mvltVNwCg+7sr1VBE9ovIuoisb25uVk5HCCEkxtwXNlX1oKququrqysrKvKcjhJD/V5jyxEVkL4CjA038JICrVXVDRPYAuFtVLy+NM8/H7nu8j7969O6h7mh5rDic3/p4sSfX1+JLrN776HZurByeHPzStbRuGzo2p9+a5+t9LNyyvlHzaLfnusbsirXxrBGU7psxGnvYNpbDncrV9q7ReL9X1vsnNpeXeeSJHwGw1h2vAThcOQ4hhJARWFIMPwjgOIDLReS0iNwE4FYA14jIKQDXdJ8JIYQsmO2lBqr66kTVvoltIYQQ4qSpvVN6LHtNePdAGM6T26+jtD+CVycb2mttb9kqNzdHqr03596T716yx3utrNSuJZTurxJjfAjPxdj9QVJtS7ZatHbLMwYlPX2sZm4ZL9Y2tNmraXvWV8bCvVMIIWRJYRAnhJCGYRAnhJCGaUoTT2lPubzRYb9QZ8zllHr2orDYUtLxS/PVaoRe+y2fa8ZK2Zc79551jCGWva5LGm3t68U82rF17SJlj9eOcPyp8rdDnyzfAe89VXNerfr7PJgyRxygJk4IIUsLgzghhDQMgzghhDRM05q4NQe7J9TQcjq0dx8Fjw/WOm+O8zxzgsP+NdfBk+fs0YtzttasJVjeu1iTOx+r92ru4eeUvTF/rOs5Jds9efY5/1LnwJqjX7N+VZqzdJ1Kvk6R8x6DmjghhCwpDOKEENIwDOKEENIwTWjiQF0+aaouHMOS4zymb6mdVTOz5raXtPBcHnRszSCVB2/NBU/h0RCn1MhL+9+kGJNHbxkz5XfuWlmekbDME/PHSu48lPb3ttyrKT/G2G5pb7l/c+2mzEunJk4IIUsKgzghhDQMgzghhDRME5q4V5uy7iVhyeWszYEu5b568khzue0W3dhrozVfvSb3PdYmZmfK55o9x626qmWcmjWBKXPYY34MseSEl2z05HaP9bGmvnZ9ymrfFPflVFp4DzVxQghZUhjECSGkYZqTU1KP2PZM8fM9VmZJlbLMk7LJIuOM3WI1Nb73MefYcWpMy/WwPnZt+VkMpCWnGkmtNiXRKjOENpbGGJOeZ6VGPivN6fHFe16t57Nkp1Uus8qWU0I5hRBClhQGcUIIaRgGcUIIaZjmNPEYY/QroPxqtHAsazpVKRUqNU8OTzpiDouOaNGmc6+8S42T0m0ten7YPpdymRojV1dzr3nWF7wpaGNT1Sz3bGrr3dhxbo5h/5jtnvu+tCZksc/TJqRGK/fM64WaOCGELCkM4oQQ0jAM4oQQ0jBNaeK53O2wbY9Ff/Xk7IZtavKyY+OV/PFo7qnxw765/GdrHnWIt583v9szr8We3HXI2eHN8x5+Hs7h1V69WrrFvqE9sfapdQjP99FyLceuj6TaW9Y+LN+L3LkJ8a59lKAmTgghSwqDOCGENMyoIC4iLxeRkyLysIgcmMooQgghNqo1cRHZBuBLAK4BcBrAvQBeraoPpPpMrYn3Zbn2sfIUFi0cSO+jUju/Je/Voy2W8Oa75/rn6odjlvp6NFDLuJbc7Rq9v6SXAuXX2uVs9q6LxGws+TqmT82zD2Pvt9oxSmtgpbFSPpbGSDFGF5+XJv58AA+r6ldU9X8BfAjAdSPGI4QQ4mRMEL8YwDcGn093ZU9CRPaLyLqIrG9ubo6YjhBCSMiYIC6RsqdoM6p6UFVXVXV1ZWVlxHSEEEJCxmjiVwF4s6q+rPt8CwCo6ltTfWo1cUII+f/MvDTxewFcJiKXisi5AG4EcGTEeIQQQpxsr+2oqo+LyJ8DuBPANgDvU9X7J7OMEEJIkeogDgCq+nEAH5/IFkIIIU74xCYhhDQMgzghhDQMgzghhDQMgzghhDQMgzghhDQMgzghhDQMgzghhDQMgzghhDTMQt+xKSKbAL5e2f0iAN+e0JyzybL4six+APRlq0JfZvy8qkZ3EFxoEB+DiKynNoBpjWXxZVn8AOjLVoW+lKGcQgghDcMgTgghDdNSED94tg2YkGXxZVn8AOjLVoW+FGhGEyeEEPJUWvqfOCGEkAAGcUIIaZgtH8RF5OUiclJEHhaRA2fbnhgi8iwR+ZSIPCgi94vI67rynSJyTEROdX93DPrc0vl0UkReNij/ZRH5Ylf3dyISeyH1vP3ZJiKfFZGjjftxoYh8WEQe6q7NVQ378hfdvXWfiHxQRM5rxRcReZ+InBGR+wZlk9kuIk8Tkdu78ntEZO+Cffmb7h77goh8TEQuXKgvqrpl/2H22rcvA3g2gHMBfB7AFWfbroidewA8rzt+JoAvAbgCwNsBHOjKDwB4W3d8RefL0wBc2vm4rav7DICrAAiAfwPwG2fBnzcA+ACAo93nVv04BOCPuuNzAVzYoi8ALgbwVQBP7z7fAeAPWvEFwIsAPA/AfYOyyWwH8KcA3t0d3wjg9gX78lIA27vjty3al4V+qSpO2FUA7hx8vgXALWfbLoPdhwFcA+AkgD1d2R4AJ2N+YPae0qu6Ng8Nyl8N4LYF234JgLsAvBg/CeIt+nEBZoFPgvIWfbkYwDcA7MTslYpHu8DRjC8A9gaBbzLb+zbd8XbMnoqURfkS1P0WgPcv0petLqf0N2/P6a5sy9L9/LkSwD0AdqvqBgB0f3d1zVJ+Xdwdh+WL5J0A3gjgx4OyFv14NoBNAP/QSUPvEZHz0aAvqvqfAP4WwCMANgB8V1U/gQZ9GTCl7U/0UdXHAXwXwM/OzfI8f4jZ/6yfZFfHXHzZ6kE8ptdt2ZxIEXkGgI8AeL2qPpZrGinTTPlCEJFrAZxR1RPWLpGys+5Hx3bMfvb+vapeCeAHmP1sT7Flfen04usw+0n+cwDOF5HX5LpEyraELwZqbN8SfonImwA8DuD9fVGk2eS+bPUgfhrAswafLwHwzbNkSxYROQezAP5+Vf1oV/wtEdnT1e8BcKYrT/l1ujsOyxfFCwG8UkS+BuBDAF4sIv+M9vxAZ8NpVb2n+/xhzIJ6i768BMBXVXVTVX8I4KMAfhVt+tIzpe1P9BGR7QB+BsCjc7M8goisAbgWwO9pp4VgQb5s9SB+L4DLRORSETkXM6H/yFm26Sl0K8vvBfCgqr5jUHUEwFp3vIaZVt6X39itRF8K4DIAn+l+Vn5PRF7Qjfn7gz5zR1VvUdVLVHUvZuf6k6r6mtb86Hz5LwDfEJHLu6J9AB5Ag75gJqO8QER+urNhH4AH0aYvPVPaPhzrdzC7bxf5C/blAG4G8EpV/Z9B1WJ8WcSixshFhFdglu3xZQBvOtv2JGz8Ncx+8nwBwOe6f6/ATMu6C8Cp7u/OQZ83dT6dxCBDAMAqgPu6undhjgs0BZ+uxk8WNpv0A8BzAax31+VfAexo2Je3AHios+OfMMt4aMIXAB/ETMv/IWb/07xpStsBnAfgXwA8jFnWx7MX7MvDmOnY/Xf/3Yv0hY/dE0JIw2x1OYUQQkgGBnFCCGkYBnFCCGkYBnFCCGkYBnFCCGkYBnFCCGkYBnFCCGmY/wOdJ9N3Tr7UUgAAAABJRU5ErkJggg==\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.LOW, 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 = 100 * 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": []
}
],
"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
}