From dbd365efb08f23370063733d2bde080859165709 Mon Sep 17 00:00:00 2001 From: bbimber Date: Thu, 14 Dec 2023 12:39:02 -0800 Subject: [PATCH] Update logic for VcfToLuceneIndexer and fractionHet to include just samples with the allele --- .../java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/main/java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java b/src/main/java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java index 01db6282..4fbd5146 100644 --- a/src/main/java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java +++ b/src/main/java/com/github/discvrseq/walkers/VcfToLuceneIndexer.java @@ -264,7 +264,7 @@ else if (line.getCountType() == VCFHeaderLineCount.INTEGER || line.getCountType( doc.add(new IntPoint("nCalled", (int)nCalled)); doc.add(new StoredField("nCalled", (int)nCalled)); - float fractionHet = (float) nHet / (float) nCalled; + float fractionHet = (float) nHet / (float) (nHet + nHomVar); doc.add(new FloatPoint("fractionHet", fractionHet)); doc.add(new StoredField("fractionHet", fractionHet));