Skip to content

Commit 20db4a4

Browse files
committed
Ignore half-called GTs
1 parent 58842c3 commit 20db4a4

File tree

1 file changed

+4
-0
lines changed

1 file changed

+4
-0
lines changed

SequenceAnalysis/src/org/labkey/sequenceanalysis/run/variant/PlinkPcaStep.java

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -227,6 +227,10 @@ private void runBatch(File inputVCF, File outputDirectory, VariantProcessingStep
227227
args.add(String.valueOf(maxRam));
228228
}
229229

230+
// NOTE: tools like sawfish can report half-called genotypes, like 0/.. For now, be most conservative in PCA:
231+
args.add("--vcf-half-call");
232+
args.add("missing");
233+
230234
args.addAll(getClientCommandArgs());
231235

232236
getWrapper().execute(args);

0 commit comments

Comments
 (0)