From f50befaff8a186b7002709b78d108e38307ea358 Mon Sep 17 00:00:00 2001
From: Yuka Suzuki <yuka.suzuki@pasteur.fr>
Date: Mon, 9 Oct 2023 14:14:43 +0200
Subject: [PATCH] initial push

---
 .../scripts/Fig_SUP19_pi_neff_dependence.py     | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)
 create mode 100644 Figures_manuscript/scripts/Fig_SUP19_pi_neff_dependence.py

diff --git a/Figures_manuscript/scripts/Fig_SUP19_pi_neff_dependence.py b/Figures_manuscript/scripts/Fig_SUP19_pi_neff_dependence.py
new file mode 100644
index 0000000..bfa2669
--- /dev/null
+++ b/Figures_manuscript/scripts/Fig_SUP19_pi_neff_dependence.py
@@ -0,0 +1,17 @@
+import pandas as pd
+import matplotlib.pyplot as plt
+import seaborn as sns
+import numpy as np
+
+df = pd.read_csv("../inputs/72trait_data_2023-07-07.csv",delimiter=',')
+df['log10_pi'] = np.log10(df['polygenicity (unadjusted)'])
+
+plt.rcParams.update({'font.family':'arial','font.size':18})
+plt.figure(figsize=[6,4])
+sns.regplot(x='Neff',y='log10_pi',data=df,ci=95)
+plt.ylim([-6,-0.9])
+plt.ylabel("$log_{10}$ polygenicity (unadjusted)")
+plt.xlabel("effective sample size")
+
+plt.tight_layout()
+plt.savefig("../outputs/Figure_SUP_pi_neff_dependence.png",dpi=300)
-- 
GitLab