Dear R community,

I'm having hard time to understand the kde function in "ks" package. Let me
use a 3-dimensional kernel smooth example to explain my question using the
elevation data in geoR.

### here is what I did ###

library(ks)
require(geoR)
data(elevation)
elev.df <- data.frame(x = elevation$coords[,"x"], y =
elevation$coords[,"y"], z = elevation$data)
H<-Hpi(elev.df)
elev.smt<-kde(elev.df,H=H)
plot(elev.smt, drawpoints=TRUE)

If I understand it correctly, with the kde function, I'm using the two
coordinate variables x and y to estimate (or say smooth) elevation (z). Is
this correct?

With this kernel smooth, my goal is to identify peaks, which are defined as
areas that have estimated elevations >=950. Can someone show me how to do
this?

Thanks!


Best
Gary

	[[alternative HTML version deleted]]

