joshuacoles/compb-dla-data-analysis
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

sp-fd graph with low p sampling

Browse Source
This commit is contained in:
Joshua Coles 2023-03-19 22:56:04 +00:00
parent cf8fb79331
commit ec2fd7d98f
3 changed files with 671 additions and 374 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
figures/sp-fd.png
View File

Binary file not shown.
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 30 KiB

After

Width:  |  Height:  |  Size: 31 KiB

1036
figures/sp-fd.svg
View File

File diff suppressed because it is too large Load Diff
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 33 KiB

After

Width:  |  Height:  |  Size: 40 KiB

9
graphs/sp-fd.py
View File

@ -1,10 +1,12 @@
import pandas as pd
from matplotlib import pyplot as plt
from lib.lib import read_sp, aggregate_sp_fd
data_2d_sp = read_sp("../data/rust-stick-probability")
data_2d_sp = read_sp("../data/rust-sticking-probability")
data_2d_sp_low_p = read_sp("../data/rust-sp-low-p")
sp_fd_data = aggregate_sp_fd(data_2d_sp)
full_data = pd.concat([data_2d_sp_low_p, data_2d_sp])
sp_fd_data = aggregate_sp_fd(full_data)
# %%
plt.fill_between(sp_fd_data.index, sp_fd_data.fd - sp_fd_data.fd_std, sp_fd_data.fd + sp_fd_data.fd_std, alpha=0.2, label=f"Standard error band")
@ -13,6 +15,9 @@ plt.plot(sp_fd_data.index, sp_fd_data.fd, color='tab:blue', label='NF, fd mean,
plt.xlabel("$p_{stick}$")
plt.ylabel("$fd$")
plt.legend()
# plt.show()
# plt.loglog(sp_fd_data.index, sp_fd_data.fd, color='tab:blue', label='NF, fd mean, seeds = 100')
# plt.show()
# %%
plt.savefig('../figures/sp-fd.svg')