Displaying 1 result from an estimated 1 matches for "k1000".
Did you mean:
1000
2009 Aug 18
0
kernel density estimation for univariate data using splancs
...and cells.dim). Here is my code:
coords <- S at polygons[[1]]@Polygons[[1]]@coords
coord <- as.points(coords)
cases1.xy <- as.points(mergedis$jit.x.x, mergedis$pre.fsa.shp.Y2)
grd1 <- GridTopology(cellcentre.offset=c(50281627.2, 580082), cellsize=c(0.2, 0.2), cells.dim=c(75,100))
k1000 <- spkernel2d(cases1.xy, coord, h0=1000, grd1)
k5000 <- spkernel2d(cases1.xy, coord, h0=5000, grd1)
if (.sp_lt_0.9()) {
df <- AttributeList(list(k1000=k1000, k5000=k5000))
} else {
df <- data.frame(k1000=k1000, k5000=k5000)
}
kernels <- SpatialGridDataFrame(grd1, data=df)
spplot(...