Changes

observed <- sort(pvals) lobs <- -(log10(observed)) expected <- c(1:length(observed)) lexp <- -(log10(expected / (length(expected)+1))) plot(c(0,7), c(0,7), col&eq;"red", lwd&eq;3, type&eq;"l", xlab&eq;"Expected (-logP)", ylab&eq;"Observed (-logP)", xlim&eq;c(0,max(lobs)), ylim&eq;c(0,max(lobs)), las&eq;1, xaxs&eq;"i", yaxs&eq;"i", bty&eq;"l", main &eq; title) points(lexp, lobs, pch&eq;23, cex&eq;.4, bg&eq;"black") }