Question: SnpSift case control command for multiple samples
0
gravatar for kumari.indu31
8 weeks ago by
kumari.indu310 wrote:

We are trying to write a script which can analyse multiple samples for SnpSift CaseControl. There is "+" symbol in the command to represent the no. of samples. How to define this option of the command in the script. Here is the script, we are trying to run, but it is not working.

!/bin/bash

for f1 in *vcf do java -jar SnpSift.jar caseControl "$f1+" $f1 > $f1_CC.vcf done

Please help with this.

ADD COMMENTlink written 8 weeks ago by kumari.indu310
1

Just a few remarks:

  • if you use for loop in one line, the syntax should be different and contain semicolons. Otherwise, you could use indentations.
  • if I am not mistaken, the usual input for caseControl is a vcf file with multiple samples in it (meaning multiple sample columns in the vcf file). In this case, the symbols +, -, 0 specify whether the case is cancer, control, or neutral.
ADD REPLYlink written 8 weeks ago by e.rempel980

My question is regarding how can I define no. of samples in the command in script so that there is no need to add "+" symbol manually.

ADD REPLYlink written 8 weeks ago by kumari.indu310

You can't run caseControl without a multi sample vcf.

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by Arup Ghosh2.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1656 users visited in the last hour
_