bkreg:

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

Usage

bkreg(x, y, kerfun = akerd, pyhat = FALSE, plotit = TRUE, xlab = "X", ylab = "Y", zlab = "Z", xout = FALSE, outfun = outpro, pr = TRUE, theta = 50, phi = 25, duplicate = "error", expand = 0.5, SCALE = FALSE, ticktype = "simple", ...)

Arguments

x
y
kerfun
pyhat
plotit
xlab
ylab
zlab
xout
outfun
pr
theta
phi
duplicate
expand
SCALE
ticktype
...

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, kerfun = akerd, pyhat = FALSE, plotit = TRUE, 
    xlab = "X", ylab = "Y", zlab = "Z", xout = FALSE, outfun = outpro, 
    pr = TRUE, theta = 50, phi = 25, duplicate = "error", expand = 0.5, 
    SCALE = FALSE, ticktype = "simple", ...) 
{
    x = as.matrix(x)
    p = ncol(x)
    p1 = p + 1
    xx <- elimna(cbind(x, y))
    x <- xx[, 1:p]
    y <- xx[, p1]
    if (xout) {
        m <- cbind(x, y)
        flag <- outfun(x, plotit = FALSE, ...)$keep
        m <- m[flag, ]
        x <- m[, 1:p]
        y <- m[, p1]
    }
    x = as.matrix(x)
    flag <- (y == 1)
    mv = sum(flag)
    nv = sum(!flag)
    phat <- NA
    fhat <- kerfun(x[flag, ], pyhat = TRUE, plotit = FALSE, pts = x)
    ghat <- kerfun(x[!flag, ], pyhat = TRUE, plotit = FALSE, 
        pts = x)
    phat <- mv * fhat/(mv * fhat + nv * ghat)
    if (p == 1) {
        if (plotit) {
            plot(x, y, xlab = xlab, ylab = ylab)
            flag2 <- order(x)
            lines(x[flag2], phat)
        }
    }
    if (p == 2) {
        if (plotit) {
            library(akima)
            if (pr) {
                if (!scale) 
                  print("With dependence, suggest using scale=T")
            }
            fitr <- phat
            iout <- c(1:length(fitr))
            nm1 <- length(fitr) - 1
            for (i in 1:nm1) {
                ip1 <- i + 1
                for (k in ip1:length(fitr)) if (sum(x[i, ] == 
                  x[k, ]) == 2) 
                  iout[k] <- 0
            }
            fitr <- fitr[iout >= 1]
            mkeep <- x[iout >= 1, ]
            fit <- interp(mkeep[, 1], mkeep[, 2], fitr, duplicate = duplicate)
            persp(fit, theta = theta, phi = phi, expand = expand, 
                scale = scale, xlab = xlab, ylab = ylab, zlab = zlab, 
                ticktype = ticktype)
        }
    }
    if (!pyhat) 
        phat <- "Done"
    phat
  }

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