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 = {}