Skip to content

Commit ea3d4f5

Browse files
committed
Ignore half-called GTs in KING
1 parent 20db4a4 commit ea3d4f5

File tree

1 file changed

+4
-0
lines changed

1 file changed

+4
-0
lines changed

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

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -115,6 +115,10 @@ public Output processVariants(File inputVCF, File outputDirectory, ReferenceGeno
115115
plinkArgs.add("--max-alleles");
116116
plinkArgs.add("2");
117117

118+
// NOTE: tools like sawfish can report half-called genotypes, like 0/.. For now, be most conservative in PCA:
119+
plinkArgs.add("--vcf-half-call");
120+
plinkArgs.add("missing");
121+
118122
Integer threads = SequencePipelineService.get().getMaxThreads(getPipelineCtx().getLogger());
119123
if (threads != null)
120124
{

0 commit comments

Comments
 (0)