efischer/BiGpairSEQ
1
0
Fork 0
Code Issues Pull Requests Projects Releases 13 Wiki Activity

checking for possible divide by zero error.

Browse Source
This commit is contained in:
efischer
2022-02-23 08:54:14 -06:00
parent dd2164c250
commit e58f7b0a55
1 changed files with 9 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/main/java/Simulator.java
View File

@@ -13,6 +13,8 @@ import java.time.Duration;
import java.util.*; import java.util.*;
import java.util.stream.IntStream; import java.util.stream.IntStream;
import static java.lang.Float.*;
//NOTE: "sequence" in method and variable names refers to a peptide sequence from a simulated T cell //NOTE: "sequence" in method and variable names refers to a peptide sequence from a simulated T cell
public class Simulator { public class Simulator {
private static final int cdr3AlphaIndex = 0; private static final int cdr3AlphaIndex = 0;
@@ -246,7 +248,13 @@ public class Simulator {
BigDecimal attemptRateTrunc = new BigDecimal(attemptRate, mc); BigDecimal attemptRateTrunc = new BigDecimal(attemptRate, mc);
//rate of pairing error //rate of pairing error
double pairingErrorRate = (double) falseCount / (trueCount + falseCount); double pairingErrorRate = (double) falseCount / (trueCount + falseCount);
BigDecimal pairingErrorRateTrunc = new BigDecimal(pairingErrorRate, mc); BigDecimal pairingErrorRateTrunc;
if(pairingErrorRate == NaN || pairingErrorRate == POSITIVE_INFINITY || pairingErrorRate == NEGATIVE_INFINITY) {
pairingErrorRateTrunc = new BigDecimal(-1, mc);
}
else{
pairingErrorRateTrunc = new BigDecimal(pairingErrorRate, mc);
}
//get list of well populations //get list of well populations
Integer[] wellPopulations = data.getWellPopulations(); Integer[] wellPopulations = data.getWellPopulations();
//make string out of populations list //make string out of populations list