Dear nest community,
I have a network of 'iaf_cond_alpha' neurons and want to use the 'dc_generator' to apply two different current amplitudes to the network. (see code below)
- A subset of the neurons should receive a sub-threshold amplitude so that they won't spike. - A subset of the neurons should receive a suprathreshold amplitude, making them spike.
I created a random number of indices of neurons that should receive a suprathreshold current, the rest of the neurons should receive the subthreshold current.
However, there doesn't seem any way to set the current amplitude to a subset of neurons. I have looked at all the examples using the dc_generator but couldn't find anything about it.
How can this be done? Is it possible? If it's not possible could you point me to the source code so I could potentially write a something for it?
Thank you.
Kind regards, Daphne
*=== code ===*
np.random.seed(0) # for reproducibility
# amplitudes Asub = 300.6 Asupra = 367.4
# list of all neurons neuron_idx_lst = np.arange(0,N_total)
# create random array, these are the neurons that should be stimulated with supra Amplitude supra_indices = np.random.randint(1, N_total, size=int(N_total/neuron_groups-1))
# the rest should receive sub threshold current sub_indices = np.setdiff1d(neuron_idx_lst,supra_indices)
# here is where it goes wrong nest.SetStatus(dc_gen[supra_indices], params={'amplitude':Asupra, 'start':stim_start, 'stop':stim_end})
nest.SetStatus(dc_gen[sub_indices], params={'amplitude':Asub, 'start':stim_start, 'stop':stim_end})
TypeError: tuple indices must be integers or slices, not tuple