tsp1reg:

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

Usage

tsp1reg(x, y, plotit = FALSE, HD = FALSE, OPT = TRUE)

Arguments

x
y
plotit
HD
OPT

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, plotit = FALSE, HD = FALSE, OPT = TRUE) 
{
    temp <- matrix(c(x, y), ncol = 2)
    temp <- elimna(temp)
    x <- temp[, 1]
    y <- temp[, 2]
    ord <- order(x)
    xs <- x[ord]
    ys <- y[ord]
    vec1 <- outer(ys, ys, "-")
    vec2 <- outer(xs, xs, "-")
    v1 <- vec1[vec2 > 0]
    v2 <- vec2[vec2 > 0]
    if (!HD) 
        slope <- median(v1/v2, na.rm = TRUE)
    if (HD) 
        slope <- hd(v1/v2, na.rm = TRUE)
    if (OPT) {
        if (!HD) 
            coef <- median(y, na.rm = TRUE) - slope * median(x, 
                na.rm = TRUE)
        if (HD) 
            coef <- hd(y, na.rm = TRUE) - slope * hd(x, na.rm = TRUE)
    }
    if (!OPT) {
        if (!HD) 
            coef <- median(y - slope * x, na.rm = TRUE)
        if (HD) 
            coef <- hd(y - slope * x, na.rm = TRUE)
    }
    names(coef) <- "Intercept"
    coef <- c(coef, slope)
    if (plotit) {
        plot(x, y, xlab = "X", ylab = "Y")
        abline(coef)
    }
    res <- y - slope * x - coef[1]
    list(coef = coef, residuals = res)
  }

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