Compute a Neighborhood Matrix from Spatial Polygons
The function takes a SpatialPolygonsDataFrame
and computes
the neighbor penalty matrix that can be used to fit a Markov
random field, e.g., using the smooth constructor
smooth.construct.mrf.smooth.spec
.
## Compute the neighborhood matrix. neighbormatrix(x, type = c("boundary", "dist", "delaunay", "knear"), k = 1, id = NULL, nb = FALSE, names = NULL, ...) ## Plot neighborhood structure. plotneighbors(x, add = FALSE, ...)
x |
An object of class |
type |
Which type of neighborhood structure should be used,
|
k |
For |
id |
An identifier variable for which the penalty matrix should be computed. |
nb |
Should only the neighborhood object be returned. |
names |
Specifies the column where the regions names are provided in the data
slot in the |
add |
Should the neighborhood structure be added to an existing plot? |
... |
Arguments to be passed to function |
data("LondonFire") ## Compute polygon boundary based ## neighborhood matrix. nm <- neighbormatrix(LondonBoroughs) print(nm) ## Plot neighborhood structures. plotneighbors(LondonBoroughs) plotneighbors(LondonBoroughs, type = "delaunay") plotneighbors(LondonBoroughs, type = "dist", d1 = 0, d2 = 0.15)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.