INNER CODE UNIT · JavaScript
main
lh3/psmc · utils/split-time.js:174
function main(args) {
let opt = { round:20, d0:1e-5, scale:100 };
for (const o of getopt(args, "N:d:s:", [])) {
if (o.opt == "-N") opt.round = parseInt(o.arg);
else if (o.opt == "-d") opt.d0 = parseFloat(o.arg);
else if (o.opt == "-s") opt.scale = parseInt(o.arg);
}
if (args.length == 0) {
print("Usage: split-time.js [options] <in1.psmc> [...]");
print("Options:");
print(` -N INT pick the INT-th interation [${opt.round}]`);
print(` -d FLOAT sequence divergence cutoff [${opt.d0}]`);
print(` -s INT PSMC scaling used for psmcfa generation [${opt.scale}]`);
print(`Output:`);
print(` 1. index of the max time interval that is smaller than {-d}`);
print(` 2. lower bound of the interval`);
print(` 3. fraction of genome coalesced <{-d} (modeled)`);
print(` 4. fraction of genome coalesced <{-d} (posterior)`);