| panel.hexloess {hexbin} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
panel.hexloess(bin, w = NULL, span = 2/3, degree = 1, family = c("symmetric", "gaussian"), evaluation = 50, lwd = add.line$lwd, lty = add.line$lty, col, col.line = add.line$col, ...)
bin |
~~Describe bin here~~ |
w |
~~Describe w here~~ |
span |
~~Describe span here~~ |
degree |
~~Describe degree here~~ |
family |
~~Describe family here~~ |
evaluation |
~~Describe evaluation here~~ |
lwd |
~~Describe lwd here~~ |
lty |
~~Describe lty here~~ |
col |
~~Describe col here~~ |
col.line |
~~Describe col.line here~~ |
... |
~~Describe ... here~~ |
~~ If necessary, more details than the description above ~~
~Describe the value returned If it is a LIST, use
comp1 |
Description of 'comp1' |
comp2 |
Description of 'comp2' |
...
....
~~further notes~~
~Make other sections like Warning with section{Warning }{....} ~
~~who you are~~
~put references to the literature/web site here ~
~~objects to See Also as help, ~~~
##---- 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 (bin, w = NULL, span = 2/3, degree = 1, family = c("symmetric",
"gaussian"), evaluation = 50, lwd = add.line$lwd, lty = add.line$lty,
col, col.line = add.line$col, ...)
{
x <- bin@xcm
y <- bin@ycm
if (is.null(w))
w <- bin@count
control <- loess.control(...)
notna <- !(is.na(x) | is.na(y))
new.x <- seq(min(x[notna]), max(x[notna]), length = evaluation)
family <- match.arg(family)
iterations <- if (family == "gaussian")
1
else control$iterations
fit <- stats:::simpleLoess(y, x, w, span, degree, FALSE,
FALSE, normalize = FALSE, "none", "interpolate", control$cell,
iterations, control$trace.hat)
kd <- fit$kd
z <- .C("loess_ifit", as.integer(kd$parameter), as.integer(kd$a),
as.double(kd$xi), as.double(kd$vert), as.double(kd$vval),
as.integer(evaluation), as.double(x), fit = double(evaluation),
PACKAGE = "stats")$fit
if (length(x) > 0) {
if (!missing(col)) {
if (missing(col.line)) {
col.line <- col
}
}
add.line <- trellis.par.get("add.line")
panel.lines(new.x, z, col = col.line, lty = lty, lwd = lwd)
}
}