diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantstovcf/PlinkToVCF.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantstovcf/PlinkToVCF.java index 420c6068b..df971cced 100644 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantstovcf/PlinkToVCF.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/variantstovcf/PlinkToVCF.java @@ -194,6 +194,9 @@ public class PlinkToVCF extends RodWalker { } } + // TODO -- REWRITE THIS TO WORK WITH VARIANT CONTEXT + /****** + private String smartHardy(ReferenceContext ref, VCFRecord rec) { HashMap> genotypesByPopulation = new HashMap>(10); HashMap hardyWeinbergByPopulation = new HashMap(10); @@ -211,8 +214,8 @@ public class PlinkToVCF extends RodWalker { for ( String population : samplesToPopulation.values() ) { VCFVariationCall v = new VCFVariationCall(ref.getBase(),ref.getLocus(),VCFVariationCall.VARIANT_TYPE.SNP); - // TODO -- v.setGenotypeCalls(genotypesByPopulation.get(population)); - // TODO -- hardyWeinbergByPopulation.put(population,HWCalc.annotate(null,ref,null,v)); + v.setGenotypeCalls(genotypesByPopulation.get(population)); + hardyWeinbergByPopulation.put(population,HWCalc.annotate(null,ref,null,v)); } return smartHardyString(hardyWeinbergByPopulation); @@ -228,6 +231,8 @@ public class PlinkToVCF extends RodWalker { return String.format("%s",maxH); } + *********/ + private HashMap parsePopulationFile(File file) { HashMap samplesToPopulation = new HashMap(); try { @@ -237,7 +242,7 @@ public class PlinkToVCF extends RodWalker { String[] populationSamples = line.split("\t"); String population = populationSamples[0]; for ( int i = 1; i < populationSamples.length; i ++ ) { - samplesToPopulation.put(populationSamples[i],population); + samplesToPopulation.put(populationSamples[i], population); } } } catch ( IOException e) {