Commit 48243ca3 authored by Hanna  JULIENNE's avatar Hanna JULIENNE
Browse files

fix p-value qqplot

parent 439d21d9
This diff is collapsed.
This diff is collapsed.
......@@ -270,8 +270,8 @@ def create_qq_plot(work_file_path: str, qq_plot_path: str):
x = -np.log10(np.arange(1, pvalue.shape[0] + 1) / (pvalue.shape[0] + 2))
y = pvalue.sort_values()
plt.scatter(x[::-1], y, s=5)
pval_median = np.nanmedian(y)
lambda_value = pval_median / np.median(x)
pval_median = df.JASS_PVAL.median()
lambda_value = pval_median / 0.5
x_1 = np.linspace(0, 6)
y_1 = lambda_value * x_1
x_2 = np.linspace(0, 6)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment