Fix bug #15

RSS need to be calculated for each clustering

bin/reports/scenic_report.ipynb

@@ -341,11 +341,6 @@
     "### Show all dimensionality reduction plots with all clusterings clustering side-by-side"
    ]
   },
-  {
-   "cell_type": "markdown",
-   "metadata": {},
-   "source": []
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -391,26 +386,6 @@
     "import seaborn as sns"
    ]
   },
-  {
-   "cell_type": "markdown",
-   "metadata": {},
-   "source": [
-    "### Calculate RSS"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "for clustering in meta[\"clusterings\"]:\n",
-    "    printmd(f'#### {clustering[\"name\"]}')\n",
-    "    clusters = pd.Series(ca_clusterings[str(clustering[\"id\"])], index=ca_cell_id)\n",
-    "    rss = regulon_specificity_scores( auc_mtx, clusters )\n",
-    "    rss"
-   ]
-  },
   {
    "cell_type": "markdown",
    "metadata": {},
@@ -432,7 +407,10 @@
     "    fig_nrows = ceil(len(cats)/fig_ncols)\n",
     "    fig_asp = 3.75 # final figure aspect ratio\n",
     "\n",
-    "\n",
+    "    ## ############\n",
+    "    # Calculate RSS\n",
+    "    rss = regulon_specificity_scores( auc_mtx, clusters )\n",
+    "    \n",
     "    fig = plt.figure(figsize=(15,15/fig_asp*fig_nrows) )\n",
     "    for c,num in zip(cats, range(1,len(cats)+1)):\n",
     "        x=rss.T[c]\n",