R/ggd.qqplot.R

In holgerman/HT12ProcessoR: Preprocessing Illuminas HT12 v4 data

ggd.qqplot <- function(pvector, main = NULL, ...) {
    # 21/8 p = 0 abgefangen
    pvector[pvector == 0] = .Machine$double.xmin
    o = -log10(sort(pvector, decreasing = F))
    e = -log10(1:length(o)/length(o))
    plot(e, o, pch = 19, cex = 1, main = main, ..., xlab = expression(Expected ~ ~-log[10](italic(p))), ylab = expression(Observed ~ ~-log[10](italic(p))), xlim = c(0, 
        max(e)), ylim = c(0, max(o)))
    lines(e, e, col = "red")
}