R/plot.multiQQ_function.R

In hutchisonjtw/JNCCTools: A Package to Store Useful Functions for JNCC Staff

#' @export

plot.multiQQ <- function(x, y, ...) {
  plot(1, type = "n", xlab = "Theoretical quantiles", ylab = "Sample quantiles", xlim = c(0, x$xmax), ylim = c(0, x$xmax), main = x$main)
  lines(rbind(c(0, 0), c(x$xmax, x$xmax)))
  mapply(points, x$theoreticalProbs, col = 1:4, pch = 1:4)
  legend("topleft", x$distr, col = 1:4, pch = 1:4)
  invisible(NULL)
}