R/pvalVolcanoPlot.R

In fcros: A Method to Search for Differentially Expressed Genes and to Detect Recurrent Chromosomal Copy Number Aberrations

pvalVolcanoPlot <- function(af, thr = 0.05) {
    ratio <- af$FC2;
    pval <- af$p.value;
    de.idx <- (pval <= thr);
    x <- log2(ratio);
    y <- -log10(pval);
    plot(x, y, ylab = "-log10(p-Value)", xlab = "log2(FC)", col = "blue", main = "Volcano plot");
    lines(x[de.idx], y[de.idx], type = "p", col = "red");
}