Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 error in data.frame (multiple treatments and multiple replicates) KYR Bioinformatics 3 03-19-2019 09:05 AM
DESeq2 without biol replicates sisterdot Bioinformatics 23 02-24-2019 08:01 AM
DESeq2 - paired formula design with replicates bryand RNA Sequencing 3 07-24-2015 02:22 AM
DESeq2 with no replicates - strange results frymor Bioinformatics 2 07-23-2014 02:17 AM
DESeq2 plotPCA combine replicates id0 Bioinformatics 2 06-17-2014 08:53 AM

Thread Tools
Old 07-13-2015, 12:57 PM   #1
Junior Member
Location: Eugene, OR

Join Date: Sep 2013
Posts: 3
Default DESeq2 with replicates no sig padj

So I am running data analysis on count data using DESeq2. I have three biological replicates for each condition. My results show significant p values but no padj values other than 1, even with very low pval. How can I troubleshoot? Truncated data set attached.

setwd("C:/Users/Melissa Randel/Desktop")
library( "DESeq" )
theta = 0.3
Fname_pos = "trunc_DESeq_out"

CountTable3L_pos = read.table("counts_table.txt")
Design3L_pos = data.frame( row.names = colnames( CountTable3L_pos ), condition = c( "untreated", "untreated", "untreated", "treated", "treated", "treated" ), libType = c( "single-end", "single-end", "single-end", "single-end", "single-end","single-end" ) )
singleSamples_pos = Design3L_pos$libType == "single-end"
countTable_pos = CountTable3L_pos[ , singleSamples_pos]
condition_pos = Design3L_pos$condition[ singleSamples_pos ]
full_cds_pos = newCountDataSet( countTable_pos, condition_pos )
rs_pos = rowSums ( counts ( full_cds_pos ))
use = (rs_pos > quantile(rs_pos, probs=theta))
cds_pos = full_cds_pos[use,]
cds_pos = estimateSizeFactors( cds_pos )
cds_pos = estimateDispersions( cds_pos )
res_pos = nbinomTest( cds_pos, "untreated", "treated" )
write.table(res_pos, file = Fname_pos, append = FALSE, quote = FALSE, sep = "\t", eol = "\n", na = "NA", dec = ".", row.names = TRUE, col.names = TRUE, qmethod = c("escape", "double"), fileEncoding = "")

Attached Files
File Type: zip (7.5 KB, 4 views)
mrandel is offline   Reply With Quote
Old 07-14-2015, 07:08 AM   #2
Senior Member
Location: California

Join Date: Jul 2014
Posts: 198

Your truncated dataset appears to (correctly) have no significant hits. I'd go through the steps in the vignette and take a look at the dispersion and MA plots. Relatedly, it looks like you're using DESeq - would recommend using DESeq2.
fanli is offline   Reply With Quote
Old 07-15-2015, 09:36 AM   #3
Michael Love
Senior Member
Location: Boston

Join Date: Jul 2013
Posts: 333

Small p-values are expected when you do many tests. One way to think about p-value correction is so you are not fooled by these small p-values which are expected under the null. Under the null, you would expect 1 p-value < 1/10,000 when you do 10,000 tests, 2 p-values less than 2/10,000, etc. If you run p-values like these (i.e. uniformly distributed) through an adjustment like the BH method, you will find there will generally be no set for which the FDR < x < 1. This is multiple test correction doing it's job. See the vignette for more description.

Also, as fanli points out, that is DESeq code, not DESeq2.
Michael Love is offline   Reply With Quote

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 08:38 AM.

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