dev_tests.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 50,
   "id": "d03d274e-6792-4bbf-93bf-b8c7259c1d7f",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The autoreload extension is already loaded. To reload it, use:\n",
      "  %reload_ext autoreload\n"
     ]
    }
   ],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2\n",
    "import pandas as pd\n",
    "from src.dysregnet.dysregnet import run"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "b6408bef-4768-42c3-88e3-21e250102240",
   "metadata": {},
   "outputs": [],
   "source": [
    "meta = pd.read_csv(\"test_data/tpm_meta.csv\")\n",
    "expr = pd.read_csv(\"test_data/tpm.csv\")\n",
    "grn=pd.read_csv(\"test_data/HTRIdb_data.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "0b876315-e883-450e-95a4-88b7caf76085",
   "metadata": {},
   "outputs": [],
   "source": [
    "meta[\"birth_days_to\"] = meta[\"birth_days_to\"].fillna(meta[\"birth_days_to\"].mean())\n",
    "meta[\"race\"] = meta[\"race\"].fillna(\"not reported\")\n",
    "meta[\"race\"] = meta[\"race\"].replace({\"[Unknown]\": \"not reported\", \"[Not Evaluated]\":\"not reported\"})\n",
    "\n",
    "expr = expr.set_index(expr.columns[0])\n",
    "expr = expr.T\n",
    "expr.insert(0, \"sample\", expr.index)\n",
    "assert all(expr.iloc[:, 0].values == meta.iloc[:,0].values)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "208046d3-5d5f-4ea4-aae7-c575094bc88e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>sample</th>\n",
       "      <th>_PATIENT</th>\n",
       "      <th>cancer type abbreviation</th>\n",
       "      <th>age_at_initial_pathologic_diagnosis</th>\n",
       "      <th>gender</th>\n",
       "      <th>race</th>\n",
       "      <th>ajcc_pathologic_tumor_stage</th>\n",
       "      <th>clinical_stage</th>\n",
       "      <th>histological_type</th>\n",
       "      <th>histological_grade</th>\n",
       "      <th>...</th>\n",
       "      <th>DSS.time</th>\n",
       "      <th>DFI</th>\n",
       "      <th>DFI.time</th>\n",
       "      <th>PFI</th>\n",
       "      <th>PFI.time</th>\n",
       "      <th>Redaction</th>\n",
       "      <th>sample_type_id</th>\n",
       "      <th>sample_type</th>\n",
       "      <th>_primary_disease</th>\n",
       "      <th>condition</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>TCGA-55-7995-01</td>\n",
       "      <td>TCGA-55-7995</td>\n",
       "      <td>LUAD</td>\n",
       "      <td>73.0</td>\n",
       "      <td>FEMALE</td>\n",
       "      <td>WHITE</td>\n",
       "      <td>Stage IA</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Lung Adenocarcinoma</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>889.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>468.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>468.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Primary Tumor</td>\n",
       "      <td>lung adenocarcinoma</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>TCGA-38-4625-11</td>\n",
       "      <td>TCGA-38-4625</td>\n",
       "      <td>LUAD</td>\n",
       "      <td>66.0</td>\n",
       "      <td>FEMALE</td>\n",
       "      <td>WHITE</td>\n",
       "      <td>Stage IB</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Lung Adenocarcinoma</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>2973.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2973.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2973.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>11.0</td>\n",
       "      <td>Solid Tissue Normal</td>\n",
       "      <td>lung adenocarcinoma</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>TCGA-69-7761-01</td>\n",
       "      <td>TCGA-69-7761</td>\n",
       "      <td>LUAD</td>\n",
       "      <td>84.0</td>\n",
       "      <td>MALE</td>\n",
       "      <td>WHITE</td>\n",
       "      <td>Stage IB</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Lung Adenocarcinoma</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>186.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.0</td>\n",
       "      <td>186.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Primary Tumor</td>\n",
       "      <td>lung adenocarcinoma</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>TCGA-67-6216-01</td>\n",
       "      <td>TCGA-67-6216</td>\n",
       "      <td>LUAD</td>\n",
       "      <td>57.0</td>\n",
       "      <td>FEMALE</td>\n",
       "      <td>WHITE</td>\n",
       "      <td>Stage IA</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Lung Adenocarcinoma</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>141.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>141.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>141.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Primary Tumor</td>\n",
       "      <td>lung adenocarcinoma</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>TCGA-44-6148-01</td>\n",
       "      <td>TCGA-44-6148</td>\n",
       "      <td>LUAD</td>\n",
       "      <td>60.0</td>\n",
       "      <td>MALE</td>\n",
       "      <td>WHITE</td>\n",
       "      <td>Stage IA</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Lung Adenocarcinoma</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>704.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>704.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>704.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Primary Tumor</td>\n",
       "      <td>lung adenocarcinoma</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 38 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            sample      _PATIENT cancer type abbreviation  \\\n",
       "0  TCGA-55-7995-01  TCGA-55-7995                     LUAD   \n",
       "1  TCGA-38-4625-11  TCGA-38-4625                     LUAD   \n",
       "2  TCGA-69-7761-01  TCGA-69-7761                     LUAD   \n",
       "3  TCGA-67-6216-01  TCGA-67-6216                     LUAD   \n",
       "4  TCGA-44-6148-01  TCGA-44-6148                     LUAD   \n",
       "\n",
       "   age_at_initial_pathologic_diagnosis  gender   race  \\\n",
       "0                                 73.0  FEMALE  WHITE   \n",
       "1                                 66.0  FEMALE  WHITE   \n",
       "2                                 84.0    MALE  WHITE   \n",
       "3                                 57.0  FEMALE  WHITE   \n",
       "4                                 60.0    MALE  WHITE   \n",
       "\n",
       "  ajcc_pathologic_tumor_stage  clinical_stage    histological_type  \\\n",
       "0                    Stage IA             NaN  Lung Adenocarcinoma   \n",
       "1                    Stage IB             NaN  Lung Adenocarcinoma   \n",
       "2                    Stage IB             NaN  Lung Adenocarcinoma   \n",
       "3                    Stage IA             NaN  Lung Adenocarcinoma   \n",
       "4                    Stage IA             NaN  Lung Adenocarcinoma   \n",
       "\n",
       "   histological_grade  ...  DSS.time  DFI  DFI.time  PFI PFI.time  Redaction  \\\n",
       "0                 NaN  ...     889.0  1.0     468.0  1.0    468.0        NaN   \n",
       "1                 NaN  ...    2973.0  0.0    2973.0  0.0   2973.0        NaN   \n",
       "2                 NaN  ...     186.0  NaN       NaN  0.0    186.0        NaN   \n",
       "3                 NaN  ...     141.0  0.0     141.0  0.0    141.0        NaN   \n",
       "4                 NaN  ...     704.0  0.0     704.0  0.0    704.0        NaN   \n",
       "\n",
       "   sample_type_id          sample_type     _primary_disease  condition  \n",
       "0             1.0        Primary Tumor  lung adenocarcinoma          1  \n",
       "1            11.0  Solid Tissue Normal  lung adenocarcinoma          0  \n",
       "2             1.0        Primary Tumor  lung adenocarcinoma          1  \n",
       "3             1.0        Primary Tumor  lung adenocarcinoma          1  \n",
       "4             1.0        Primary Tumor  lung adenocarcinoma          1  \n",
       "\n",
       "[5 rows x 38 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "meta.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "35febb56-1441-4049-a452-450c25afebb1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th>sample</th>\n",
       "      <th>sample</th>\n",
       "      <th>RP11-34P13.7</th>\n",
       "      <th>RP11-34P13.8</th>\n",
       "      <th>CICP27</th>\n",
       "      <th>RP11-34P13.15</th>\n",
       "      <th>RP11-34P13.16</th>\n",
       "      <th>RP11-34P13.13</th>\n",
       "      <th>FO538757.2</th>\n",
       "      <th>AP006222.2</th>\n",
       "      <th>RP4-669L17.10</th>\n",
       "      <th>...</th>\n",
       "      <th>MT-CO2</th>\n",
       "      <th>MT-ATP8</th>\n",
       "      <th>MT-ATP6</th>\n",
       "      <th>MT-CO3</th>\n",
       "      <th>MT-ND3</th>\n",
       "      <th>MT-ND4L</th>\n",
       "      <th>MT-ND4</th>\n",
       "      <th>MT-ND5</th>\n",
       "      <th>MT-ND6</th>\n",
       "      <th>MT-CYB</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>TCGA-55-7995-01</th>\n",
       "      <td>TCGA-55-7995-01</td>\n",
       "      <td>0.496253</td>\n",
       "      <td>0.562163</td>\n",
       "      <td>1.819967</td>\n",
       "      <td>1.945791</td>\n",
       "      <td>1.013340</td>\n",
       "      <td>1.489258</td>\n",
       "      <td>2.809733</td>\n",
       "      <td>2.720060</td>\n",
       "      <td>1.448061</td>\n",
       "      <td>...</td>\n",
       "      <td>0.620082</td>\n",
       "      <td>-0.891616</td>\n",
       "      <td>0.031685</td>\n",
       "      <td>0.590567</td>\n",
       "      <td>-0.847172</td>\n",
       "      <td>0.537382</td>\n",
       "      <td>0.866038</td>\n",
       "      <td>0.486255</td>\n",
       "      <td>-0.036586</td>\n",
       "      <td>-0.059644</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TCGA-38-4625-11</th>\n",
       "      <td>TCGA-38-4625-11</td>\n",
       "      <td>-0.110067</td>\n",
       "      <td>0.014913</td>\n",
       "      <td>-1.962869</td>\n",
       "      <td>-0.813442</td>\n",
       "      <td>-0.416048</td>\n",
       "      <td>-1.713702</td>\n",
       "      <td>-0.940562</td>\n",
       "      <td>-1.039181</td>\n",
       "      <td>-1.737634</td>\n",
       "      <td>...</td>\n",
       "      <td>0.151399</td>\n",
       "      <td>-0.685961</td>\n",
       "      <td>-0.306375</td>\n",
       "      <td>0.289113</td>\n",
       "      <td>0.895153</td>\n",
       "      <td>-0.455624</td>\n",
       "      <td>-0.069085</td>\n",
       "      <td>-0.627751</td>\n",
       "      <td>-0.722210</td>\n",
       "      <td>-0.002905</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TCGA-69-7761-01</th>\n",
       "      <td>TCGA-69-7761-01</td>\n",
       "      <td>-0.110067</td>\n",
       "      <td>0.578951</td>\n",
       "      <td>0.341960</td>\n",
       "      <td>-0.030543</td>\n",
       "      <td>0.101305</td>\n",
       "      <td>0.602601</td>\n",
       "      <td>-1.764752</td>\n",
       "      <td>-1.116434</td>\n",
       "      <td>0.600953</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.695560</td>\n",
       "      <td>0.842728</td>\n",
       "      <td>0.008627</td>\n",
       "      <td>-0.768548</td>\n",
       "      <td>-2.694695</td>\n",
       "      <td>-0.043072</td>\n",
       "      <td>-0.480471</td>\n",
       "      <td>-0.274686</td>\n",
       "      <td>-0.205755</td>\n",
       "      <td>-1.255776</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TCGA-67-6216-01</th>\n",
       "      <td>TCGA-67-6216-01</td>\n",
       "      <td>-0.414540</td>\n",
       "      <td>-2.486586</td>\n",
       "      <td>-0.101231</td>\n",
       "      <td>1.162652</td>\n",
       "      <td>0.742126</td>\n",
       "      <td>0.936864</td>\n",
       "      <td>-1.340492</td>\n",
       "      <td>-1.116434</td>\n",
       "      <td>0.747482</td>\n",
       "      <td>...</td>\n",
       "      <td>0.152220</td>\n",
       "      <td>-0.215530</td>\n",
       "      <td>0.056840</td>\n",
       "      <td>0.926111</td>\n",
       "      <td>-2.544535</td>\n",
       "      <td>0.081287</td>\n",
       "      <td>-0.070977</td>\n",
       "      <td>-0.378933</td>\n",
       "      <td>-1.056571</td>\n",
       "      <td>-0.535812</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TCGA-44-6148-01</th>\n",
       "      <td>TCGA-44-6148-01</td>\n",
       "      <td>1.098818</td>\n",
       "      <td>0.595261</td>\n",
       "      <td>0.728725</td>\n",
       "      <td>0.670473</td>\n",
       "      <td>0.439598</td>\n",
       "      <td>0.729149</td>\n",
       "      <td>0.512986</td>\n",
       "      <td>1.445461</td>\n",
       "      <td>0.868339</td>\n",
       "      <td>...</td>\n",
       "      <td>0.323585</td>\n",
       "      <td>-0.281548</td>\n",
       "      <td>0.094381</td>\n",
       "      <td>0.249448</td>\n",
       "      <td>0.895798</td>\n",
       "      <td>0.406246</td>\n",
       "      <td>0.483969</td>\n",
       "      <td>-0.844608</td>\n",
       "      <td>-1.144351</td>\n",
       "      <td>-0.227255</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 22579 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "sample                    sample  RP11-34P13.7  RP11-34P13.8    CICP27  \\\n",
       "TCGA-55-7995-01  TCGA-55-7995-01      0.496253      0.562163  1.819967   \n",
       "TCGA-38-4625-11  TCGA-38-4625-11     -0.110067      0.014913 -1.962869   \n",
       "TCGA-69-7761-01  TCGA-69-7761-01     -0.110067      0.578951  0.341960   \n",
       "TCGA-67-6216-01  TCGA-67-6216-01     -0.414540     -2.486586 -0.101231   \n",
       "TCGA-44-6148-01  TCGA-44-6148-01      1.098818      0.595261  0.728725   \n",
       "\n",
       "sample           RP11-34P13.15  RP11-34P13.16  RP11-34P13.13  FO538757.2  \\\n",
       "TCGA-55-7995-01       1.945791       1.013340       1.489258    2.809733   \n",
       "TCGA-38-4625-11      -0.813442      -0.416048      -1.713702   -0.940562   \n",
       "TCGA-69-7761-01      -0.030543       0.101305       0.602601   -1.764752   \n",
       "TCGA-67-6216-01       1.162652       0.742126       0.936864   -1.340492   \n",
       "TCGA-44-6148-01       0.670473       0.439598       0.729149    0.512986   \n",
       "\n",
       "sample           AP006222.2  RP4-669L17.10  ...    MT-CO2   MT-ATP8   MT-ATP6  \\\n",
       "TCGA-55-7995-01    2.720060       1.448061  ...  0.620082 -0.891616  0.031685   \n",
       "TCGA-38-4625-11   -1.039181      -1.737634  ...  0.151399 -0.685961 -0.306375   \n",
       "TCGA-69-7761-01   -1.116434       0.600953  ... -0.695560  0.842728  0.008627   \n",
       "TCGA-67-6216-01   -1.116434       0.747482  ...  0.152220 -0.215530  0.056840   \n",
       "TCGA-44-6148-01    1.445461       0.868339  ...  0.323585 -0.281548  0.094381   \n",
       "\n",
       "sample             MT-CO3    MT-ND3   MT-ND4L    MT-ND4    MT-ND5    MT-ND6  \\\n",
       "TCGA-55-7995-01  0.590567 -0.847172  0.537382  0.866038  0.486255 -0.036586   \n",
       "TCGA-38-4625-11  0.289113  0.895153 -0.455624 -0.069085 -0.627751 -0.722210   \n",
       "TCGA-69-7761-01 -0.768548 -2.694695 -0.043072 -0.480471 -0.274686 -0.205755   \n",
       "TCGA-67-6216-01  0.926111 -2.544535  0.081287 -0.070977 -0.378933 -1.056571   \n",
       "TCGA-44-6148-01  0.249448  0.895798  0.406246  0.483969 -0.844608 -1.144351   \n",
       "\n",
       "sample             MT-CYB  \n",
       "TCGA-55-7995-01 -0.059644  \n",
       "TCGA-38-4625-11 -0.002905  \n",
       "TCGA-69-7761-01 -1.255776  \n",
       "TCGA-67-6216-01 -0.535812  \n",
       "TCGA-44-6148-01 -0.227255  \n",
       "\n",
       "[5 rows x 22579 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "expr.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "9dd72783-5b5d-488d-937e-f1d34535b29a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>SYMBOL_TF</th>\n",
       "      <th>SYMBOL_TG</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>PARP1</td>\n",
       "      <td>BRCA2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>AHR</td>\n",
       "      <td>CYP1A1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>AHR</td>\n",
       "      <td>CYP1A2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>AHR</td>\n",
       "      <td>CYP1B1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>AHR</td>\n",
       "      <td>FOS</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  SYMBOL_TF SYMBOL_TG\n",
       "0     PARP1     BRCA2\n",
       "1       AHR    CYP1A1\n",
       "2       AHR    CYP1A2\n",
       "3       AHR    CYP1B1\n",
       "4       AHR       FOS"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "grn.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "abda2742-f17d-440c-b5b2-e021d2c26f9c",
   "metadata": {},
   "outputs": [],
   "source": [
    "CatCov=['gender']  \n",
    "ConCov=['birth_days_to',]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "id": "62b921d3-27e0-4f54-a77b-d3b35d5eedfb",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "14979it [00:48, 307.28it/s]\n"
     ]
    }
   ],
   "source": [
    "data=run(expression_data=expr,\n",
    "                   meta=meta, \n",
    "                   CatCov=CatCov,\n",
    "                   ConCov=ConCov,\n",
    "                   GRN=grn,\n",
    "                   conCol='condition',\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "c3aede6a-7044-456e-affd-3804f54eb9d7",
   "metadata": {},
   "outputs": [],
   "source": [
    "data_old = data.get_results()\n",
    "data_binary_old = data.get_results_binary()\n",
    "model_stats_old = data.get_model_stats()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "id": "f93aab62-29c1-483f-b08c-8b7505917834",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n",
      "0\n",
      "128425\n"
     ]
    }
   ],
   "source": [
    "print((data.get_results() != data_old).sum().sum())\n",
    "print((data.get_results_binary() != data_binary_old).sum().sum())\n",
    "print((data.get_model_stats() != model_stats_old).sum().sum())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "id": "a77d69fa-9401-466a-a20a-29d738649016",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1.6502188025642124\n"
     ]
    }
   ],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}