update rho growth normalization

ArcInstitute · Dec 11, 2023 · de465e0 · de465e0
1 parent decc0c0
commit de465e0
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/screenpro/__init__.py b/screenpro/__init__.py
@@ -60,8 +60,9 @@ def calculateDrugScreen(self, t0, untreated, treated, db_untreated, db_treated,
             n_reps=self.n_reps,
             math=self.math, test=self.test, score_level=score_level
         )
+        # TO-DO: warning / error if db_untreated and db_treated are too close, i.e. growth_rate ~= 0.
         rho_name, rho = ps.runPhenoScore(
-            self.adata, cond1=untreated, cond2=treated, growth_rate=np.abs(db_untreated - (1 - db_treated)),
+            self.adata, cond1=untreated, cond2=treated, growth_rate=np.abs(db_untreated - db_treated),
             n_reps=self.n_reps,
             math=self.math, test=self.test, score_level=score_level
         )