Nothing
# ellipsoidellipsoidhull around partitions
".ellipsoidhull" <- function (x, ...) {
# x <- prt
cl <- match.call()
res <- vector("list", length = getK(x))
for (i in 1:getK(x)) {
xy <- as.matrix(coordinates(x)[partitioning(x) == i,])
if (nrow(xy) >= 3) {
res[[i]] <- cluster::ellipsoidhull(x = xy, ...)
} else {
if (ncol(xy) > 1) {
res[[i]] <- xy.coords(x = xy[, 1], y = xy[, 2])
} else {
res[[i]] <- xy.coords(x = t(xy)[, 1], y = t(xy)[, 2])
}
}
}
points(SpatialPointsVegsoup(x))
if (any(names(cl) == "col")) {
#stopifnot()
}
#
sapply(res, function (x) {
if (class(x) == "ellipsoid") {
lines(predict(x), ...) # col = col)
} else {
points(x, ...) # col = col, cex = 2)
}
})
return(invisible(cl))
}
# cluster defines:
# ellipsoidhull(x, tol=0.01, maxit=5000, ret.wt = FALSE,
# ret.sqdist = FALSE, ret.pr = FALSE)
#if (!isGeneric("ellipsoidhull")) {
setGeneric("ellipsoidhull",
function (x, tol = 0.01, maxit = 5000, ret.wt = FALSE,
ret.sqdist = FALSE, ret.pr = FALSE)
standardGeneric("ellipsoidhull")
)
#}
setMethod("ellipsoidhull",
signature(x = "VegsoupPartition"),
function (x, tol = 0.01, maxit = 5000, ret.wt = FALSE,
ret.sqdist = FALSE, ret.pr = FALSE)
.ellipsoidhull(x, tol, maxit, ret.wt, ret.sqdist, ret.pr)
)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.