Skip to content

Instantly share code, notes, and snippets.

@mbk0asis

mbk0asis/volcano_plot_R

Last active September 25, 2015 07:35
Show Gist options
  • Save mbk0asis/fbaf7839cd9b5e74ed9e to your computer and use it in GitHub Desktop.
Save mbk0asis/fbaf7839cd9b5e74ed9e to your computer and use it in GitHub Desktop.
Download ZIP
Raw
volcano_plot_R
with(SCNT_TSA, plot(log2FoldChange,-log10(pval),pch=20,main="Volcano plot",
col=rgb(0.6,0.6,0.6, alpha = 0.1)))
with(subset(SCNT_TSA, pval<.05 & abs(log2FoldChange)>1),
points(log2FoldChange,-log10(pval),pch=20,col=rgb(1,0.5,0, alpha = 0.4)))
abline(v=c(-1,1),col="blue",lwd=1,lty=5)
abline(v=c(0,0),col="grey40",lwd=1,lty=2)
abline(c(-log10(0.05),0),col="blue",lwd=1,lty=5)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment