From ae6e78f9a749481a2ba56035c0bfe17b900549c4 Mon Sep 17 00:00:00 2001 From: Ramona Date: Mon, 19 Nov 2018 17:12:09 +0100 Subject: [PATCH] still problems with access --- code/spikes_analysis.py | 14 +++++++++----- 1 file changed, 9 insertions(+), 5 deletions(-) diff --git a/code/spikes_analysis.py b/code/spikes_analysis.py index e059c38..2d72447 100644 --- a/code/spikes_analysis.py +++ b/code/spikes_analysis.py @@ -31,21 +31,25 @@ for deltaf in df_map.keys(): for idx in range(len(phase_vec)-1): if phase[1] > phase_vec[idx] and phase[1] < phase_vec[idx+1]: if phase_vec[idx] in phase_df.keys(): - phase_df[phase_vec[idx]].append(phase[1]) + phase_df[phase_vec[idx]].append(phase) else: - phase_df[phase_vec[idx]] = [phase[1]] + phase_df[phase_vec[idx]] = [phase] #spikes_one_chirp = spikes[rep][phase] phase_mat_df[deltaf] = phase_df -embed() -exit() -plot_spikes = spikes[(0, '-50Hz', '20%', '100Hz')][(0, 0.789)] +#df_spikes = spikes[df_map['-50Hz'][0]] +#phase_spikes = df_spikes[phase_mat_df['-50Hz'][0.0][0]] +trial_spikes = np.asarray(spikes[(0, '-50Hz', '20%', '100Hz')][(0, 0.789)]) +plot_spikes = trial_spikes[(trial_spikes < 50.0) & (trial_spikes > -50.0)] + +hist_data = plt.hist(plot_spikes, color='w') fig, ax = plt.subplots() ax.scatter(plot_spikes, np.ones(len(plot_spikes))*10, marker='|', color='k') +ax.plot(hist_data[1][:-1], hist_data[0]) plt.show() #mu = 1