R/fdrBiCurve.R

Defines functions fdrBiCurve

Documented in fdrBiCurve

fdrBiCurve <-
function(fdr.table,maxLogP=5,...) {
	fdrt=fdr.table
	plot(dr~logp,fdrt,xlim=c(0.1,maxLogP),type="l",xaxt="n",mgp=c(2.3,1,0),xlab="pvalue",ylab="N(DEGs)",...)
	axis(1,labels=10^(-c(1:maxLogP)),at=c(1:maxLogP)) 
	lines(fdr~logp,fdrt,xlim=c(1,maxLogP),type="l",col="red")
	abline(v=-log(empiricalFDR(fdrt),10),lty=3)
	legend("topright",c("real","simul"),lty=1,col=c("black","red"),bty="n")
}

Try the empiricalFDR.DESeq2 package in your browser

Any scripts or data that you put into this service are public.

empiricalFDR.DESeq2 documentation built on May 30, 2017, 4:23 a.m.