write correlation to correlation.txt

chanzuckerberg · Aug 29, 2023 · 93e3039 · 93e3039
1 parent 140461d
commit 93e3039
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/workflows/benchmark/notebooks/short-read-mngs-benchmarks.ipynb b/workflows/benchmark/notebooks/short-read-mngs-benchmarks.ipynb
@@ -240,8 +240,12 @@
     "if ref_data: \n",
     "    correlations = joined[joined[\"max_rPM\"] >= min_rPM].fillna(0)\n",
     "    res = scipy.stats.spearmanr(correlations[\"max_rPM\"], correlations[\"max_rPM.REF\"])\n",
+    "    correlation = round(res.correlation, 2)\n",
+    "    with open(\"correlation.txt\", \"w\") as f:\n",
+    "        f.write(str(correlation))\n",
+    "\n",
     "    fig, ax = plt.subplots(1, 2, figsize=(12,6))\n",
-    "    fig.suptitle(f\"Correlation plots: Spearman's R {round(res.correlation, 2)}\")\n",
+    "    fig.suptitle(f\"Correlation plots: Spearman's R {correlation}\")\n",
     "    sns.scatterplot(ax=ax[0], data=correlations, x=\"max_rPM\", y=\"max_rPM.REF\")\n",
     "    ax[0].set_title(\"Raw scatterplot\")\n",
     "\n",