error in cuffdiff ---(mutiple sample comparision)
2
0
Entering edit mode
9.0 years ago

Hello Everyone,

I am using cuffdiff to compare my 2samples against the control here I am using following command for analysis.

cuffdiff \
  -p 10 \
  -o /home/pinky/Desktop/NGS/diff_out/son/ \
  --labels conditionA1, conditionB1, conditionC1 \
  -b /home/pinky/Desktop/NGS/conindex/rice_exc.fasta \
  -u /home/pinky/Desktop/NGS/cuffmerge/sonexp/transcripts.gtf \
  /home/pinky/Desktop/NGS/tophat_fasta/sample19/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample20/accepted_hits.bam /home/pinky/Desktop/NGS/tophat_fasta/sample21/accepted_hits.bam /home/pinky/Desktop/NGS/tophat_fasta/sample25/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample26/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample27/accepted_hits.bam /home/pinky/Desktop/NGS/tophat_fasta/sample31_1/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample32/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample33/accepted_hits.bam

Warning: Your version of Cufflinks is not up-to-date. It is recommended that you upgrade to Cufflinks v2.2.1 to benefit from the most recent features and bug fixes (http://cufflinks.cbcb.umd.edu).
Error: number of labels must match number of conditions

but instead of proper run I am getting this error. So I manipulated it in this way to resolve this.

cuffdiff \
  -p 10 \
  -o /home/pinky/Desktop/NGS/diff_out/son/ \
  --labels conditionA1, conditionA2, conditionA3, conditionB1, conditionB2, conditionB3, conditionC1, conditionC2, conditionC3 \
  -b /home/pinky/Desktop/NGS/conindex/rice_exc.fasta \
  -u /home/pinky/Desktop/NGS/cuffmerge/sonexp/transcripts.gtf /home/pinky/Desktop/NGS/tophat_fasta/sample19/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample20/accepted_hits.bam /home/pinky/Desktop/NGS/tophat_fasta/sample21/accepted_hits.bam /home/pinky/Desktop/NGS/tophat_fasta/sample25/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample26/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample27/accepted_hits.bam /home/pinky/Desktop/NGS/tophat_fasta/sample31_1/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample32/accepted_hits.bam, /home/pinky/Desktop/NGS/tophat_fasta/sample33/accepted_hits.bam

but still the problem remains same. so will you please guide or suggest something.

Thanks
Megha Varshney

ngs cufflinks RNA-Seq cuffdiff • 6.8k views
ADD COMMENT
0
Entering edit mode

Thanks @Goutham Atla and @ Daniel, I got solved my problem as previously I had made some mistakes.

ADD REPLY
6
Entering edit mode
9.0 years ago
Daniel ★ 4.0k

Remove the spaces between your commas, as cuffdiff is interpreting them as separate strings (it is expecting one comma separated list).

What you want to provide is the same number of strings as your conditions, so it looks like this but with full paths:

cuffdiff -p 10 -o diff_out/son/ --labels conditionA1,conditionB1,conditionC1 -b rice_exc.fasta -u transcripts.gtf sample1.bam,sample2.bam,sample3.bam sample4.bam,sample5.bam,sample6.bam sample7.bam,sample8.bam,sample9.bam

If it would make it easier to read, you can separate your string lines with a \ (press: space-backslash-newline) which only counts as a space, but separates your text, like so:

cuffdiff -p 10 -o diff_out/son/ --labels conditionA1,conditionB1,conditionC1 -b rice_exc.fasta -u transcripts.gtf \
sample1.bam,sample2.bam,sample3.bam \
sample4.bam,sample5.bam,sample6.bam \
sample7.bam,sample8.bam,sample9.bam
ADD COMMENT
2
Entering edit mode
9.0 years ago

From your command, it looks like that you have four groups of bam files (Each group is separated by space, within each groups replicates are separated by comma) and 8 conditions. Each condition separated by comma.

Basic format is

-L  A,B,C   SampleA_rep1.bam,SampleA_rep2.bam,SampleA_rep3.bam      SampleB_rep1.bam,SampleB_rep2.bam,SampleB_rep3.bam    SampleC_rep1.bam,SampleC_rep2.bam,SampleC_rep3.bam

The documentation may be misleading as it says "Specify a label for each sample, which will be included in various output files produced by Cuffdiff" but label should be given for each group.

ADD COMMENT

Login before adding your answer.

Traffic: 2182 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6