diff --git a/code/base_chirps.py b/code/base_chirps.py new file mode 100644 index 0000000..aa8f838 --- /dev/null +++ b/code/base_chirps.py @@ -0,0 +1,43 @@ +from read_chirp_data import * +#import nix_helpers as nh +import matplotlib.pyplot as plt +import numpy as np +from IPython import embed + + +data_dir = "../data" +dataset = "2018-11-09-ad-invivo-1" +#data = ("2018-11-09-ad-invivo-1", "2018-11-09-ae-invivo-1", "2018-11-09-ag-invivo-1", "2018-11-13-aa-invivo-1", "2018-11-13-ac-invivo-1", "2018-11-13-ad-invivo-1", "2018-11-13-ah-invivo-1", "2018-11-13-ai-invivo-1", "2018-11-13-aj-invivo-1", "2018-11-13-ak-invivo-1", "2018-11-13-al-invivo-1", "2018-11-14-aa-invivo-1", "2018-11-14-ac-invivo-1", "2018-11-14-ad-invivo-1", "2018-11-14-af-invivo-1", "2018-11-14-ag-invivo-1", "2018-11-14-ah-invivo-1", "2018-11-14-ai-invivo-1", "2018-11-14-ak-invivo-1", "2018-11-14-al-invivo-1", "2018-11-14-am-invivo-1", "2018-11-14-an-invivo-1") +spikes = read_chirp_spikes(os.path.join(data_dir, dataset)) +eod = read_chirp_eod(os.path.join(data_dir, dataset)) +times = read_chirp_times(os.path.join(data_dir, dataset)) + + +df_map = {} +for k in spikes.keys(): + df = k[1] + if df in df_map.keys(): + df_map[df].append(k) + else: + df_map[df] = [k] +print(df_map.keys()) + +e1 = eod[0, '-50Hz', '20%', '100Hz'] +#plt.plot(e1[1]) +#plt.show() +plt.title('EOD chirps') +plt.xlabel('Frequency') +plt.ylabel('Amplitude') +plt.plot(e1[0],e1[1]) +plt.show() + +ct = times[0, '-50Hz', '20%', '100Hz'] +#plt.scatter(ct*1000, np.ones(len(ct))) +plt.scatter(np.asarray(ct)*1000, np.ones(len(ct))) +plt.show() + +#plt.scatter(spikes[0, '-50Hz', '20%', '100Hz'][0, 0.614]) + +#print(len(spikes)) +#print(len(eod)) +#print(len(times)) diff --git a/code/liste.py b/code/liste.py index cba7d05..7a51eda 100644 --- a/code/liste.py +++ b/code/liste.py @@ -1,43 +1,43 @@ # 9.11.18 -aa: quality: poor, depth: -1341 -ab: quality: poor, depth: -1341 -ac: quality: good, depth: -1341 -ad: quality: good, depth: -200 -ae: quality: good, depth: -200 +aa: quality: poor, depth: -1341, base +ab: quality: poor, depth: -1341, base +ac: quality: good, depth: -1341, base +ad: quality: good, depth: -200, base, chirps +ae: quality: good, depth: -200, chirps af: quality: good, depth: -200 -ag: no info.dat +ag: no info.dat, chirps # 13.11.18 -aa: good, -30 µm, maybe no reaction -ab: good, -309 µm -ac: poor, -309 µm -ad: fair, -360 µm +aa: good, -30 µm, maybe no reaction, base, chirps +ab: good, -309 µm, base +ac: poor, -309 µm, chirps +ad: fair, -360 µm, base, chirps ae: fair, -350 µm -af: good, -440 µm, bursting -ag: fair, -174 µm -ah: good, -209 µm -ai: good, -66.9 µm -aj: good, -132 µm -ak: good, -284 µm -al: good, -286 µm +af: good, -440 µm, bursting, base +ag: fair, -174 µm, base +ah: good, -209 µm, base, chirps, FI, SAM +ai: good, -66.9 µm, base, chirps, SAM +aj: good, -132 µm, base, chirps +ak: good, -284 µm, base, chirps +al: good, -286 µm, base, chirps, SAM # 14.11.18 -aa: good, -184 µm -ab: fair, -279 µm, no reaction -ac: fair, -60 µm -ad: good, -357 µm -ae: fair, -357 µm -af: fair, -527 µm, only baseline -ag: fair, -533 µm -ah: poor, -505 µm -ai: good, -500 µm, still same cell 3x -aj: poor, -314 µm,no modulation -ak: good, -140 µm -al: good, -280 µm -am: good, -320 µm -an: good, -434 µm +aa: good, -184 µm, base, chirps, FI, SAM, noise +ab: fair, -279 µm, no reaction, base +ac: fair, -60 µm, base, chirps +ad: good, -357 µm, base, chirps +ae: fair, -357 µm, base +af: fair, -527 µm, base, (chirps) +ag: fair, -533 µm, base, chirps +ah: poor, -505 µm, chirps +ai: good, -500 µm, still same cell 3x, chirps, FI, noise +aj: poor, -314 µm, no modulation, base +ak: good, -140 µm, base, chirps, FI, SAM, noise +al: good, -280 µm, base, chirps, SAM +am: good, -320 µm, base, chirps, FI, SAM, noise +an: good, -434 µm, base, chirps, FI, SAM, noise