View Single Post
Old 01-04-2012, 10:55 PM   #1
idyll_ty
Junior Member
 
Location: LA

Join Date: Nov 2011
Posts: 5
Default Why edgeR ouput some p values larger than 1

I was using edgeR to do differential expression on a RNA-seq dataset, but when I check the p.value distribution of the output, I found ~700 genes got a p value larger than 1... Just wondering if anyone else got this problem before? And why I got such a weird a result. Here's my codes,

d <- DGEList(counts=countsTable, group=conds, lib.size = lib.sizes) #I used the raw read counts as instructed
d <- calcNormFactors(d, method="TMM")
d <- estimateCommonDisp(d)
et.commonDis <- exactTest(d)

hist(et.commonDis$table$p.value, breaks=100, col="skyblue", border="slateblue", main="pval distribution_edgeR")

I got a picture as attached...

May I have your ideas? Thanks.
Attached Images
File Type: png Screen shot 2012-01-04 at 10.53.27 PM.png (34.6 KB, 20 views)
idyll_ty is offline   Reply With Quote