qregsm:

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

Usage

qregsm(x, y, est = hd, qval = 0.5, sm = TRUE, plotit = TRUE, pyhat = FALSE, fr = 0.8, nboot = 40, xlab = "X", ylab = "Y")

Arguments

x
y
est
qval
sm
plotit
pyhat
fr
nboot
xlab
ylab

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, est = hd, qval = 0.5, sm = TRUE, plotit = TRUE, 
    pyhat = FALSE, fr = 0.8, nboot = 40, xlab = "X", ylab = "Y") 
{
    x <- as.matrix(x)
    X <- cbind(x, y)
    X <- elimna(X)
    np <- ncol(X)
    p <- np - 1
    x <- X[, 1:p]
    x <- as.matrix(x)
    y <- X[, np]
    vals <- matrix(NA, ncol = length(y), nrow = length(qval))
    for (i in 1:length(qval)) {
        if (sm) 
            vals[i, ] <- rplotsm(x, y, est = est, q = qval[i], 
                pyhat = TRUE, plotit = FALSE, fr = fr, nboot = nboot)$yhat
        if (!sm) 
            vals[i, ] <- rungen(x, y, est = est, q = qval[i], 
                pyhat = TRUE, plotit = FALSE, fr = fr)$output
    }
    if (p == 1) {
        if (plotit) {
            plot(x, y, xlab = xlab, ylab = ylab)
            for (i in 1:length(qval)) {
                sx <- sort(x)
                xorder <- order(x)
                sysm <- vals[i, ]
                lines(sx, sysm[xorder])
            }
        }
    }
    output <- "Done"
    if (pyhat) 
        output <- vals
    output
  }

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