[R-sig-Geo] illustrate kernel density estimation as wireframe

Marco Helbich marco.helbich at gmx.at
Sat Oct 28 18:10:42 CEST 2006


hello,

how can i illustrate a kernel density estimation as an perspective plot or 
wireframe? enclosed you can find an kde example-code...

library(splancs)
data(uganda)
plot(uganda$poly, type = "n", xlab = "x", ylab = "y", main = "Bsp KDE: 
Craters in Uganda")
polygon(uganda$poly)
points(uganda)
events <- as.points(uganda$x, uganda$y)
kde <- kernel2d(events, uganda$poly, h0 = 200, 500, 500, kernel="quartic")
image(kde, add = TRUE)
polygon(uganda$poly)
points(uganda)

thanks for reply and best regards
marco




More information about the R-sig-Geo mailing list