def get_model_parts_sweep_from_number(sn,data_set,sweep_numbers,specimen_id):
sweep_numbers = sweep_numbers['Square - 2s Suprathreshold']
rheobase = -1
above_threshold_sn = []
currents = {}
sweep_data = data_set.get_sweep(sn)
stimulus = sweep_data['stimulus']
spike_times = data_set.get_spike_times(sn)
reponse = sweep_data['response']
sampling_rate = sweep_data['sampling_rate']
vmm = AnalogSignal([v*1000 for v in reponse],sampling_rate=sampling_rate*qt.Hz,units=qt.mV)
vmm = vmm[0:int(len(vmm)/2.1)]
vmm.sn = None
vmm.sn = sn
return vmm,stimulus,sn,spike_times
def get_model_parts_sweep_from_spk_cnt(spk_cnt, data_set, sweep_numbers,
specimen_id):
sweep_numbers = sweep_numbers['Square - 2s Suprathreshold']
rheobase = -1
above_threshold_sn = []
for sn in sweep_numbers:
spike_times = data_set.get_spike_times(sn)
if len(spike_times) >= spk_cnt:
# stimulus is a numpy array in amps
sweep_data = data_set.get_sweep(sn)
stimulus = sweep_data['stimulus']
reponse = sweep_data['response']
if len(spike_times):
thresh_ = len(np.where(reponse > 0))
sampling_rate = sweep_data['sampling_rate']
vmm = AnalogSignal([v * 1000 for v in reponse],
sampling_rate=sampling_rate * qt.Hz,
units=qt.mV)
vmm = vmm[0:int(len(vmm) / 2.1)]
vmm.sn = None
vmm.sn = sn
return vmm, stimulus, sn, spike_times
return None, None, None, None
if (len(spike_times) and np.max(reponse)<0) or len(spike_times) != thresh_:
reponse = [v+np.mean([0,np.abs(np.max(reponse))]) for v in reponse]
if np.max(reponse)<0:
reponse = [v+np.mean([0,np.abs(np.max(reponse))]) for v in reponse]
assert np.max(reponse)>0
assert np.min(reponse)<0
'''
<<<<<<< HEAD
=======
>>>>>>> 6072f03bf1e58526d60ee74e6513066a07f5a462
sampling_rate = sweep_data['sampling_rate']
vmm = AnalogSignal([v*1000 for v in reponse],sampling_rate=sampling_rate*qt.Hz,units=qt.mV)
vmm = vmm[0:int(len(vmm)/2.1)]
vmm.sn = None
vmm.sn = sn
<<<<<<< HEAD
#print(vmm)
=======
>>>>>>> 6072f03bf1e58526d60ee74e6513066a07f5a462
return vmm,stimulus,sn,spike_times
return None,None,None,None
#from elephant.signal_processing import threshold_detection
def get_model_parts_sweep_from_number(sn,data_set,sweep_numbers,specimen_id):
sweep_numbers = sweep_numbers['Square - 2s Suprathreshold']
rheobase = -1
above_threshold_sn = []
currents = {}
