Carrying capacity depends on number of drivers

Author: melkebir

src/simulation/cell.cpp

@@ -37,7 +37,7 @@ Cell::Outcome Cell::performGeneration(const AnatomicalSiteFactorVector& anatomic
 
   double logistic_factor = map.find(_driverMutations)->second;
 
-  double birthRate = 0.65 * (1 + s_0) * (_driverMutations.size() + 1) * logistic_factor;
+  double birthRate = 0.5 * (1 + s_0) * (_driverMutations.size() + 1) * logistic_factor;
   if (birthRate < 0)
     birthRate = 0;
   else if (birthRate > 1)

src/simulation/celltree.cpp

@@ -297,7 +297,7 @@ void CellTree::updateAnatomicalSiteFactors()
       nrCells_s += nrCells_sX;
 
 //      _anatomicalSiteFactors[s][X] = double(_K - nrCells_sX) / _K;
-      _anatomicalSiteFactors[s][X] = 1 - double(nrCells_sX) / _K;
+      _anatomicalSiteFactors[s][X] = 1 - double(nrCells_sX) / (_K * (1 + X.size()));
 
       if (nrCells_sX > 0)
       {