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

Add read depth simulation options to CLI

Browse Source
This commit is contained in:
eugenefischer
2022-09-27 14:35:55 -05:00
parent 7c3c95ab4b
commit 423c9d5c93
1 changed files with 18 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
src/main/java/CommandLineInterface.java
View File

@@ -146,7 +146,17 @@ public class CommandLineInterface {
CellSample cells = getCells(cellFilename); CellSample cells = getCells(cellFilename);
//get plate //get plate
Plate plate = getPlate(plateFilename); Plate plate = getPlate(plateFilename);
GraphWithMapData graph = Simulator.makeGraph(cells, plate, 1, 0.0, 0.0, false); GraphWithMapData graph;
if(line.hasOption("rd")){
String[] readDepthArgs = line.getOptionValues("rd");
Integer readDepth = Integer.parseInt(readDepthArgs[0]);
Double readErrorRate = Double.parseDouble(readDepthArgs[1]);
Double errorCollisionRate = Double.parseDouble(readDepthArgs[2]);
graph = Simulator.makeGraph(cells, plate, readDepth, readErrorRate, errorCollisionRate, false);
}
else{
graph = Simulator.makeGraph(cells, plate, 1, 0.0, 0.0, false);
}
if (!line.hasOption("no-binary")) { //output binary file unless told not to if (!line.hasOption("no-binary")) { //output binary file unless told not to
GraphDataObjectWriter writer = new GraphDataObjectWriter(outputFilename, graph, false); GraphDataObjectWriter writer = new GraphDataObjectWriter(outputFilename, graph, false);
writer.writeDataToFile(); writer.writeDataToFile();
@@ -388,11 +398,18 @@ public class CommandLineInterface {
.longOpt("no-binary") .longOpt("no-binary")
.desc("(Optional) Don't output serialized binary file") .desc("(Optional) Don't output serialized binary file")
.build(); .build();
Option simulateReadDepth = Option.builder("rd")
.longOpt("read-depth")
.desc("(Optional) Simulate read depth and read errors. Has three arguments: the read depth (integer >= 1), the read error probability (0.0 - 1.0), and the error collision probability (0.0 - 1.0)")
.hasArgs()
.numberOfArgs(3)
.build();
graphOptions.addOption(cellFilename); graphOptions.addOption(cellFilename);
graphOptions.addOption(plateFilename); graphOptions.addOption(plateFilename);
graphOptions.addOption(outputFileOption()); graphOptions.addOption(outputFileOption());
graphOptions.addOption(outputGraphML); graphOptions.addOption(outputGraphML);
graphOptions.addOption(outputSerializedBinary); graphOptions.addOption(outputSerializedBinary);
graphOptions.addOption(simulateReadDepth);
return graphOptions; return graphOptions;
} }