1 | rplot.res(x, y, pv = 1, est = tmean, scat = TRUE, fr = NA, plotit = TRUE, pyhat = FALSE, efr = 0.5, theta = 50, phi = 25, scale = TRUE, expand = 0.5, SEED = TRUE, varfun = pbvar, outfun = outpro, STAND = TRUE, nmin = 0, xout = FALSE, out = FALSE, eout = FALSE, xlab = "X", ylab = "Y", zscale = FALSE, zlab = " ", pr = TRUE, duplicate = "error", ticktype = "simple", LP = TRUE, ...)
|
x |
|
y |
|
pv |
|
est |
|
scat |
|
fr |
|
plotit |
|
pyhat |
|
efr |
|
theta |
|
phi |
|
scale |
|
expand |
|
SEED |
|
varfun |
|
outfun |
|
STAND |
|
nmin |
|
xout |
|
out |
|
eout |
|
xlab |
|
ylab |
|
zscale |
|
zlab |
|
pr |
|
duplicate |
|
ticktype |
|
LP |
|
... |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | ##---- 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, pv = 1, est = tmean, scat = TRUE, fr = NA, plotit = TRUE,
pyhat = FALSE, efr = 0.5, theta = 50, phi = 25, scale = TRUE,
expand = 0.5, SEED = TRUE, varfun = pbvar, outfun = outpro,
STAND = TRUE, nmin = 0, xout = FALSE, out = FALSE, eout = FALSE,
xlab = "X", ylab = "Y", zscale = FALSE, zlab = " ", pr = TRUE,
duplicate = "error", ticktype = "simple", LP = TRUE, ...)
{
xy = na.omit(cbind(x, y))
p = ncol(x)
p1 = p + 1
if (xout) {
flag = outfun(xy[, 1:p], plotit = FALSE, STAND = STAND,
...)$keep
xy = xy[flag, ]
}
x = xy[, 1:p]
y = xy[, p1]
res = y - rplot(x[, 1:2], y, est = est, scat = scat, varfun = varfun,
expand = expand, nmin = nmin, pyhat = TRUE, plotit = FALSE,
fr = fr, xout = FALSE)$yhat
outp = rplot(x[, pv], res, fr = fr, xout = FALSE, efr = efr,
theta = theta, phi = phi, scale = scale, SEED = SEED,
xlab = xlab, ylab = ylab, zlab = zlab, pr = FALSE, ticktype = ticktype,
LP = LP, ...)
outp
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.