SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
EdgeR - default offset drea11 Bioinformatics 3 08-18-2014 05:20 AM
Help with RSEM: how to change default phred33 to phred64 LearningBIOINFO RNA Sequencing 7 03-19-2014 10:21 PM
edgeR: How important is the FDR value? polsum Bioinformatics 23 02-07-2014 08:58 AM
edgeR: number of DGE genes and FDR?? SpreeFu Bioinformatics 4 03-12-2013 08:54 AM
Default Change in CASAVA / BCL->FASTQ skruglyak Bioinformatics 20 11-19-2012 12:51 PM

Reply
 
Thread Tools
Old 09-14-2015, 06:47 PM   #1
tinguzman
Member
 
Location: Philippines

Join Date: Aug 2014
Posts: 13
Default change default FDR on edgeR

Hi,

Could anyone please tell me how to change the FDR from default (0.05) to 0.0001 in edgeR?

et <- exactTest(y, pair=c("1","2"))
out.adjusted <- topTags(et, n=Inf, adjust.method="BH")
out.adjusted.keep <- out$table$FDR <= 0.0001

I'm not sure if i understand this correctly, topTags function with adjust.method will adjust the pvalues from et, using BH/FDR set at 5%, which is the default. then, the FDR column of the output file are the adjusted p.values, on this case retaining only those with <=0.0001 values.

appreciate your help!
tinguzman is offline   Reply With Quote
Old 09-15-2015, 01:10 AM   #2
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 310
Default

The BH/FDR method corrects your pvalues according to their number and distribution, there is no threshold to set there. Maybe the 5% you refer to is the usual threshold people use to filter results. So, if you want to keep genes with FDR < 0.0001 it seems to me that your code is fine (by the way <0.0001 is very stringent, if you have many genes maybe consider shareholding also on logFC and logCPM).
dariober is offline   Reply With Quote
Reply

Tags
degs, edger, fdr

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 04:02 AM.


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