logadr:

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

Usage

logadr(x, y, est = mean, iter = 10, pyhat = FALSE, plotit = TRUE, fr = 0.8, xout = FALSE, eout = xout, outfun = out, theta = 50, phi = 25, expand = 0.5, STAND = TRUE, ticktype = "simple", scale = FALSE, ...)

Arguments

x
y
est
iter
pyhat
plotit
fr
xout
eout
outfun
theta
phi
expand
STAND
ticktype
scale
...

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 = mean, iter = 10, pyhat = FALSE, plotit = TRUE, 
    fr = 0.8, xout = FALSE, eout = xout, outfun = out, theta = 50, 
    phi = 25, expand = 0.5, STAND = TRUE, ticktype = "simple", 
    scale = FALSE, ...) 
{
    m <- elimna(cbind(x, y))
    x <- as.matrix(x)
    p <- ncol(x)
    p1 <- p + 1
    y <- m[, p1]
    x <- m[, 1:p]
    x <- as.matrix(x)
    if (STAND) {
        for (ip in 1:p) x[, ip] <- (x[, ip] - mean(x[, ip]))/sqrt(var(x[, 
            ip]))
    }
    if (xout) {
        keepit <- rep(T, nrow(x))
        flag <- outfun(x, plotit = FALSE)$out.id
        keepit[flag] <- F
        x <- x[keepit, ]
        y <- y[keepit]
    }
    x <- as.matrix(x)
    if (p == 1) 
        val <- logrsm(x[, 1], y, pyhat = T, plotit = plotit, 
            fr = fr, ...)$output
    if (p > 1) {
        np <- p + 1
        x <- m[, 1:p]
        y <- m[, np]
        fhat <- matrix(NA, ncol = p, nrow = length(y))
        fhat.old <- matrix(NA, ncol = p, nrow = length(y))
        res <- matrix(NA, ncol = np, nrow = length(y))
        dif <- 1
        for (i in 1:p) fhat.old[, i] <- logrsm(x[, i], y, pyhat = TRUE, 
            plotit = FALSE, fr = fr)$output
        eval <- NA
        for (it in 1:iter) {
            for (ip in 1:p) {
                res[, ip] <- y
                for (ip2 in 1:p) {
                  if (ip2 != ip) 
                    res[, ip] <- res[, ip] - fhat.old[, ip2]
                }
                fhat[, ip] = logrsm(x[, ip], y, pyhat = TRUE, 
                  plotit = FALSE, fr = fr)$output
            }
            eval[it] <- sum(abs(fhat/sqrt(sum(fhat^2)) - fhat.old/sqrt(sum(fhat.old^2))))
            if (it > 1) {
                itm <- it - 1
                dif <- abs(eval[it] - eval[itm])
            }
            fhat.old <- fhat
            if (dif < 0.01) 
                break
        }
        val <- apply(fhat, 1, sum)
        aval <- est(y - val, ...)
        val <- val + aval
        flag = (val < 0)
        val[flag] = 0
        flag = (val > 1)
        val[flag] = 1
        if (plotit && p == 2) {
            fitr <- val
            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, ]
            fitr <- interp(mkeep[, 1], mkeep[, 2], fitr)
            persp(fitr, theta = theta, phi = phi, expand = expand, 
                xlab = "x1", ylab = "x2", zlab = "", scale = scale, 
                ticktype = ticktype)
        }
    }
    if (!pyhat) 
        val <- "Done"
    val
  }

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