qhdsm:

In musto101/wilcox_R: What the package does (short line)

Usage

qhdsm(x, y, qval = 0.5, q = NULL, pr = FALSE, xout = FALSE, outfun = outpro, plotit = TRUE, xlab = "X", ylab = "Y", zlab = "Z", pyhat = FALSE, fr = NULL, LP = TRUE, theta = 50, phi = 25, ticktype = "simple", nmin = 0, scale = TRUE, pr.qhd = TRUE, pch = ".", ...)

Arguments

x
y
qval
q
pr
xout
outfun
plotit
xlab
ylab
zlab
pyhat
fr
LP
theta
phi
ticktype
nmin
scale
pr.qhd
pch
...

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, y, qval = 0.5, q = NULL, pr = FALSE, xout = FALSE, 
    outfun = outpro, plotit = TRUE, xlab = "X", ylab = "Y", zlab = "Z", 
    pyhat = FALSE, fr = NULL, LP = TRUE, theta = 50, phi = 25, 
    ticktype = "simple", nmin = 0, scale = TRUE, pr.qhd = TRUE, 
    pch = ".", ...) 
{
    if (!is.null(q)) 
        qval = q
    x <- as.matrix(x)
    X <- cbind(x, y)
    X <- elimna(X)
    np <- ncol(X)
    p <- np - 1
    if (p > 1 & length(q) > 1) 
        print("Only first quantile specified can be plotted")
    x <- X[, 1:p]
    x <- as.matrix(x)
    y <- X[, np]
    if (xout) {
        x <- as.matrix(x)
        flag <- outfun(x, plotit = FALSE, ...)$keep
        x <- x[flag, ]
        y <- y[flag]
        x <- as.matrix(x)
    }
    if (p == 1) {
        if (is.null(fr)) 
            fr = 0.8
        ord = order(x)
        x = sort(x)
        y = y[ord]
        est = matrix(NA, ncol = 3, nrow = length(qval))
        dimnames(est) = list(NULL, c("q", "Inter", "Slope"))
        x <- as.matrix(x)
        qest = matrix(NA, ncol = length(qval), nrow = length(y))
        for (j in 1:length(qval)) {
            rmd = NA
            for (i in 1:length(x)) rmd[i] <- hd(y[near(x, x[i], 
                fr)], q = qval[j])
            if (LP) 
                rmd = lplot(x, rmd, plotit = FALSE, pyhat = TRUE, 
                  pr = FALSE)$yhat
            qest[, j] = rmd
        }
        if (plotit) {
            plot(x, y, xlab = xlab, ylab = ylab, pch = pch)
            for (j in 1:ncol(qest)) lines(x, qest[, j])
        }
        if (!pyhat) 
            qest = "DONE"
    }
    if (p > 1) {
        if (is.null(fr)) 
            fr = 1
        if (p == 2) {
            if (pr.qhd) {
                if (!scale) 
                  print("scale=F is specified. If there is dependence, might want to use scale=TRUE")
            }
        }
        qest = rplot(x, y, est = hd, q = qval[1], fr = fr, plotit = plotit, 
            pyhat = pyhat, theta = theta, phi = phi, scale = scale, 
            SEED = FALSE, varfun = pbvar, xlab = xlab, ylab = ylab, 
            zlab = zlab, ticktype = ticktype, nmin = nmin, pr = pr)
        if (!pyhat) 
            qest = "DONE"
        if (pyhat) 
            qest = qest$yhat
    }
    qest
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.