Skip to content
Snippets Groups Projects
Commit 44654f47 authored by Peter van 't Hof's avatar Peter van 't Hof
Browse files

Style changes

parent 21f352a2
No related branches found
No related tags found
No related merge requests found
...@@ -66,25 +66,26 @@ class GwasTest(val root: Configurable) extends QScript with BiopetQScript with R ...@@ -66,25 +66,26 @@ class GwasTest(val root: Configurable) extends QScript with BiopetQScript with R
val cvTotal = new CatVariants(this) val cvTotal = new CatVariants(this)
cvTotal.assumeSorted = true cvTotal.assumeSorted = true
cvTotal.outputFile = new File(outputDirGens, "merge.gens.vcf.gz") cvTotal.outputFile = new File(outputDirGens, "merge.gens.vcf.gz")
val chrGens = inputGens.groupBy(_.contig).map { case (contig, gens) => val chrGens = inputGens.groupBy(_.contig).map {
val cvChr = new CatVariants(this) case (contig, gens) =>
cvChr.assumeSorted = true val cvChr = new CatVariants(this)
//cvChr.isIntermediate = true cvChr.assumeSorted = true
cvChr.outputFile = new File(outputDirGens, s"${contig}.merge.gens.vcf.gz") //cvChr.isIntermediate = true
gens.zipWithIndex.foreach { gen => cvChr.outputFile = new File(outputDirGens, s"${contig}.merge.gens.vcf.gz")
val gensToVcf = new GensToVcf(this) gens.zipWithIndex.foreach { gen =>
gensToVcf.inputGens = gen._1.genotypes val gensToVcf = new GensToVcf(this)
gensToVcf.inputInfo = gen._1.info gensToVcf.inputGens = gen._1.genotypes
gensToVcf.contig = gen._1.contig gensToVcf.inputInfo = gen._1.info
gensToVcf.samplesFile = phenotypeFile gensToVcf.contig = gen._1.contig
gensToVcf.outputVcf = new File(outputDirGens, gen._1.genotypes.getName + s".${gen._2}.vcf.gz") gensToVcf.samplesFile = phenotypeFile
gensToVcf.isIntermediate = true gensToVcf.outputVcf = new File(outputDirGens, gen._1.genotypes.getName + s".${gen._2}.vcf.gz")
add(gensToVcf) gensToVcf.isIntermediate = true
cvChr.inputFiles :+= gensToVcf.outputVcf add(gensToVcf)
} cvChr.inputFiles :+= gensToVcf.outputVcf
add(cvChr) }
cvTotal.inputFiles :+= cvChr.outputFile add(cvChr)
contig -> cvChr.outputFile cvTotal.inputFiles :+= cvChr.outputFile
contig -> cvChr.outputFile
} }
add(cvTotal) add(cvTotal)
(cvTotal.outputFile, Map[String, File]()) (cvTotal.outputFile, Map[String, File]())
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment