R/qcomhdMC.R

In mkanai/WRShd: An Experimental Refactor of the Harrell-Davis Estimate Functions from R. R. Wilcox' Robust Statistics Package "WRS"

#' Compare quantiles via the Harrell-Davis estimator
#' 
#' @param x a numeric \code{vector}
#' @param y a numeric \code{vector}
#' @param q a desired quantile
#' @param nboot a number of bootstraps
#' @param plotit a logical indicating whether plotting the result
#' @param SEED a logical indicating whether the seed should be set
#' @param xlab a string of x-axis label
#' @param ylab a string of y-axis label
#' @param alpha a significance level
#' @param cores a number of cores used for computation
#' @export
qcomhdMC <- function(x, y, q = c(.1, .25, .5, .75, .9), nboot = 2000, plotit = TRUE, SEED = TRUE, xlab = "Group 1", ylab = "Est.1-Est.2", alpha = .05, cores = -1) {
    require(parallel)
    if(cores <= 0){cores <- parallel::detectCores()}
    qcomhd(x, y, q, nboot, plotit, SEED, xlab, ylab, alpha, cores)
}