SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsiste Newsbot! Literature Watch 0 08-09-2012 02:00 AM
0 value in RPKM (RNA-seq) cuekt RNA Sequencing 0 02-14-2012 10:54 AM
RNA-seq, RPKM and heatmap??? sdwy2008 RNA Sequencing 10 07-14-2011 10:41 PM

Reply
 
Thread Tools
Old 10-17-2013, 10:53 PM   #1
joyjane88
Junior Member
 
Location: CAU,Beijing,China

Join Date: Mar 2013
Posts: 4
Red face how to caculate rpkm for bacteria RNA-SEQ data?

hello, everyone! I am a new member here.
Recently, I am working with the bacterial RNA-SEQ data. My data was in strand-specific PE reads fq file. I have mapped them onto the reference genome sequence and have got the sorted and indexed bam file. Nextly,
I want to make a differential gene expression analysis. So, I want to know what tools is recommended by the hands-on. DEseq, EdgeR, Cuffdiff or any other one, which is better?
Any opinions or Suggestions will be appreciated. Thank you!
joyjane88 is offline   Reply With Quote
Old 10-18-2013, 12:33 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

I've generally been happy with DESeq2. As to which is better, there's really no single answer to that. They all use slightly different methods, so it's not unheard of for one package to model experiment A (but not B) better than the others.
dpryan is offline   Reply With Quote
Old 10-18-2013, 01:18 AM   #3
joyjane88
Junior Member
 
Location: CAU,Beijing,China

Join Date: Mar 2013
Posts: 4
Default

Quote:
Originally Posted by dpryan View Post
I've generally been happy with DESeq2. As to which is better, there's really no single answer to that. They all use slightly different methods, so it's not unheard of for one package to model experiment A (but not B) better than the others.
Do you mean that all these tools are suitable for my data and will produce very similar results? I have seen a paper shown that DEseq identify more DEGs than EdgeR when the authors dealing their data with these two pakages however they select the EdgeR result for subsequent study, can I come to a conclusion that DEseq is more senstive?
joyjane88 is offline   Reply With Quote
Old 10-18-2013, 01:25 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Quote:
Originally Posted by joyjane88 View Post
Do you mean that all these tools are suitable for my data and will produce very similar results? I have seen a paper shown that DEseq identify more DEGs than EdgeR when the authors dealing their data with these two pakages however they select the EdgeR result for subsequent study, can I come to a conclusion that DEseq is more senstive?
Yes, they'll all work for your data and probably produce pretty similar results. I would expect that DESeq2 will produce a little better results than the others (including DESeq1). I would say that DESeq2 is probably more sensitive, simply due to the shrinkage methods. Aside from that, they're all pretty similar.
dpryan is offline   Reply With Quote
Old 10-18-2013, 02:32 AM   #5
joyjane88
Junior Member
 
Location: CAU,Beijing,China

Join Date: Mar 2013
Posts: 4
Smile

Quote:
Originally Posted by dpryan View Post
Yes, they'll all work for your data and probably produce pretty similar results. I would expect that DESeq2 will produce a little better results than the others (including DESeq1). I would say that DESeq2 is probably more sensitive, simply due to the shrinkage methods. Aside from that, they're all pretty similar.
Thanks for your suggestions and explaination! I will undertake my work with DEseq2.
joyjane88 is offline   Reply With Quote
Reply

Tags
paired-ends, rna-seq bacterial, rna-seq counts, strand-specific

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:20 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO