raab/GP2023_chirp_detection
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

new introplot{

Browse Source
This commit is contained in:
weygoldt 2023-01-26 14:54:37 +01:00
parent 5cf6e3f9c0
commit 5ffaf2054b
2 changed files with 20 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

39
code/plot_kdes.py
View File

@ -249,15 +249,15 @@ def main(dataroot):
# time_before=time_before, # time_before=time_before,
# time_after=time_after)) # time_after=time_after))
loser_offsets_jackknife = jackknife( # loser_offsets_jackknife = jackknife(
loser, # loser,
nresamples=nbootstraps, # nresamples=nbootstraps,
subsetsize=0.5, # subsetsize=0.9,
kde_time=kde_time, # kde_time=kde_time,
kernel_width=kernel_width, # kernel_width=kernel_width,
event_times=offsets, # event_times=offsets,
time_before=time_before, # time_before=time_before,
time_after=time_after) # time_after=time_after)
if plot_all: if plot_all:
@ -289,17 +289,17 @@ def main(dataroot):
ax[i].axvline(0, color=ps.gray, linestyle='--') ax[i].axvline(0, color=ps.gray, linestyle='--')
ax[i].fill_between( # ax[i].fill_between(
kde_time, # kde_time,
np.percentile(loser_offsets_jackknife, 5, axis=0), # np.percentile(loser_offsets_jackknife, 5, axis=0),
np.percentile(loser_offsets_jackknife, 95, axis=0), # np.percentile(loser_offsets_jackknife, 95, axis=0),
color=ps.blue, # color=ps.blue,
alpha=0.5) # alpha=0.5)
ax[i].plot(kde_time, np.median(loser_offsets_jackknife, axis=0), # ax[i].plot(kde_time, np.median(loser_offsets_jackknife, axis=0),
color=ps.white, linewidth=2) # color=ps.white, linewidth=2)
ax[i].set_xlim(-60, 60) ax[i].set_xlim(-60, 60)
fig.supylabel('Chirp rate (Hz)', fontsize=14) fig.supylabel('Chirp rate (a.u.)', fontsize=14)
fig.supxlabel('Time (s)', fontsize=14) fig.supxlabel('Time (s)', fontsize=14)
# fig, ax = plt.subplots(2, 3, figsize=( # fig, ax = plt.subplots(2, 3, figsize=(
@ -519,7 +519,8 @@ def main(dataroot):
# loser_physicals_boot_quarts[2], # loser_physicals_boot_quarts[2],
# color=ps.gray, # color=ps.gray,
# alpha=0.5) # alpha=0.5)
plt.subplots_adjust(bottom=0.21)
plt.savefig('../poster/figs/kde.pdf')
plt.show() plt.show()

BIN
poster/figs/kde.pdf Normal file
View File

Binary file not shown.