gamplot:

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

Usage

gamplot(x, y, sop = TRUE, pyhat = FALSE, eout = FALSE, xout = FALSE, outfun = out, plotit = TRUE, xlab = "X", ylab = "", zlab = "", theta = 50, phi = 25, expand = 0.5, scale = TRUE, ticktype = "simple")

Arguments

x
y
sop
pyhat
eout
xout
outfun
plotit
xlab
ylab
zlab
theta
phi
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, sop = TRUE, pyhat = FALSE, eout = FALSE, xout = FALSE, 
    outfun = out, plotit = TRUE, xlab = "X", ylab = "", zlab = "", 
    theta = 50, phi = 25, expand = 0.5, scale = TRUE, ticktype = "simple") 
{
    library(akima)
    library(mgcv)
    x <- as.matrix(x)
    np <- ncol(x)
    np1 <- np + 1
    if (ncol(x) > 4) 
        stop("x should have at most four columns of data")
    m <- elimna(cbind(x, y))
    x <- m[, 1:np]
    x <- as.matrix(x)
    y <- m[, np1]
    if (xout && eout) 
        stop("Can't have xout=eout=T")
    if (eout) {
        flag <- outfun(m)$keep
        m <- m[flag, ]
    }
    if (xout) {
        flag <- outfun(x, plotit = FALSE)$keep
        m <- m[flag, ]
    }
    x <- m[, 1:np]
    x <- as.matrix(x)
    y <- m[, np1]
    if (!sop) {
        if (ncol(x) == 1) 
            fitr <- fitted(gam(y ~ x[, 1]))
        if (ncol(x) == 2) 
            fitr <- fitted(gam(y ~ x[, 1] + x[, 2]))
        if (ncol(x) == 3) 
            fitr <- fitted(gam(y ~ x[, 1] + x[, 2] + x[, 3]))
        if (ncol(x) == 4) 
            fitr <- fitted(gam(y ~ x[, 1] + x[, 2] + x[, 3] + 
                x[, 4]))
    }
    if (sop) {
        if (ncol(x) == 1) 
            fitr <- fitted(gam(y ~ s(x[, 1])))
        if (ncol(x) == 2) 
            fitr <- fitted(gam(y ~ s(x[, 1]) + s(x[, 2])))
        if (ncol(x) == 3) 
            fitr <- fitted(gam(y ~ s(x[, 1]) + s(x[, 2]) + s(x[, 
                3])))
        if (ncol(x) == 4) 
            fitr <- fitted(gam(y ~ s(x[, 1]) + s(x[, 2]) + s(x[, 
                3]) + s(x[, 4])))
    }
    last <- fitr
    if (plotit) {
        if (ncol(x) == 1) {
            plot(x, fitr, xlab = xlab, ylab = ylab)
        }
        if (ncol(x) == 2) {
            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 = xlab, ylab = ylab, zlab = zlab, scale = scale, 
                ticktype = ticktype)
        }
    }
    if (!pyhat) 
        last <- "Done"
    last
  }

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