Nothing
# Copyright 2001-2022 by Nicholas Lewin-Koh and Roger Bivand
#
relativeneigh <- function(coords, nnmult=3) {
if (inherits(coords, "SpatialPoints")) {
if (!is.na(is.projected(coords)) && !is.projected(coords)) {
warning("relativeneigh: coordinates should be planar")
}
coords <- coordinates(coords)
} else if (inherits(coords, "sfc")) {
if (!inherits(coords, "sfc_POINT"))
stop("Point geometries required")
if (attr(coords, "n_empty") > 0L)
stop("Empty geometries found")
if (!is.na(sf::st_is_longlat(coords)) && sf::st_is_longlat(coords))
warning("relativeneigh: coordinates should be planar")
coords <- sf::st_coordinates(coords)
} else if (inherits(coords, "data.frame")) {
coords <- as.matrix(coords)
}
x <- coords
if (!is.matrix(x)) stop("Data not in matrix form")
if (any(is.na(x))) stop("Data cannot include NAs")
np <- nrow(x)
if(ncol(x)!=2) stop("Only works in 2d")
ngaballoc <- np*nnmult
g1<-g2<-rep(0,ngaballoc)
nogab <- 0
storage.mode(x) <- "double"
z <- .C("compute_relative", np=as.integer(np), from=as.integer(g1),
to=as.integer(g2), nedges=as.integer(nogab),
ngaballoc=as.integer(ngaballoc), x=x[,1],
y=x[,2], PACKAGE="spdep")
z$from<-z$from[1:z$nedges]
z$to<-z$to[1:z$nedges]
attr(z, "call") <- match.call()
class(z)<-c("Graph","relative")
z
}
plot.relative<-function(x, show.points=FALSE, add=FALSE, linecol=par(col),...){
if(!add) plot(x$x,x$y,type='n')
segments(x$x[x$from],x$y[x$from],
x$x[x$to],x$y[x$to],col=linecol)
if(show.points) points(x$x,x$y,...)
}
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.