SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
No significant DE genes in cuffdiff RocheKermit RNA Sequencing 7 01-13-2014 03:02 AM
cummeRbund significant genes cnyh RNA Sequencing 0 03-06-2013 04:17 AM
Cufflinks2 analysis without significant genes rboettcher Bioinformatics 3 11-08-2012 05:58 AM
RNA-Seq: Length Bias Correction for RNA-seq Data in Gene Set Analyses. Newsbot! Literature Watch 0 01-22-2011 02:02 AM

Reply
 
Thread Tools
Old 03-26-2013, 07:30 AM   #1
fahim
Member
 
Location: louisville

Join Date: Oct 2010
Posts: 10
Default No significant genes in RNA-seq analyses

Hi all
I am doing RNA-seq analyses using tophat and cufflinks. I followed all the steps of nature protocol papers without any problem except the 'cummRbund' portion.

I am adding selected snippet results below.

Code:
> 
cuff_data <- readCufflinks("diff_out")
> cuff_data
CuffSet instance with:
	 2 samples
	 35210 genes
	 65828 isoforms
	 46226 TSS
	 24611 CDS
	 35210 promoters
	 46226 splicing
	 20203 relCDS
> csVolcano(genes(cuff_data), 'C1', 'C2')
Warning message:
Removed 6924 rows containing missing values (geom_point). 
> gene.diff[gene.diff$significant=='yes',]
 [1] gene_id          sample_1         sample_2         status           value_1          value_2          log2_fold_change test_stat       
 [9] p_value          q_value          significant     
<0 rows> (or 0-length row.names)
> getSig(cuff_data,alpha=0.05,level='genes')
character(0)
> mySigGeneIds<-getSig(cuff_data,alpha=0.5,level='genes')
> getSig(cuff_data,alpha=0.5,level='genes')
character(0)
> length(getSig(cuff_data,alpha=0.9,level='genes'))
[1] 6759
I am amazed not to see anything differentially expressed. I am wondering if I did any mistake during the entire run.
My 'volcanoplot' also does not show any significant point in the label, however, I see some points that are above 2.0 on y-axis, demonstrating a p-value of 0.001.

I read some other posts suggesting to use edgeR.


Thanks and help appreciated

Last edited by fahim; 03-26-2013 at 07:55 AM.
fahim is offline   Reply With Quote
Reply

Tags
cummerbund 2, rna-seq

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 02:30 PM.


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