Can't group samples in StringTie for differential expression
0
0
Entering edit mode
8.6 years ago
CandiceChuDVM ★ 2.5k

Hi all,

I have performed StringTie for individual samples and --merge for all samples by implementing the following commands:

$ stringtie $PATH1/sample1.aligns.sorted.bam -G $ANNOT -p 20 -o $PATH2/sample1.gtf     
$ stringtie --merge -G $ANNOT stringtie_gtf_list -o $PATH2/merged.gtf

Now I am going to run StringTie using the -B/-b and -e options for each sample in order to estimate transcript abundances and generate read coverage tables for Ballgown.

However, I am not sure how to group my samples. For example, I would like to group all 9 samples into 3 different groups:

group1: sample1, sample2, sample3
group2: sample4, sample5, sample6
group3: sample7, sample8, sample9

Do I have to execute stringtie three times (one group at a time)? For example:

$ stringtie $PATH1/sample1.aligns.sorted.bam,$PATH1/sample2.aligns.sorted.bam,$PATH1/sample3.aligns.sorted.bam \
 -G $PATH2/merged.gtf -B -e -o group1.stringtie.out.gtf

Or can I run them all together by separating them with different labels as in Cuffdiff?

Thanks!

RNA-Seq StringTie differential expression • 2.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 1882 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