data(chamois)
xy <- chamois$locs
## Try a number of neighbours between 6 and 18
(nn <- NNCH(xy, k = seq(6, 18, by=3)))
(i <- NNCH.k.area(nn, percent=100))
plot(i)
(i <- NNCH.der.k.area(nn, percent=100))
plot(i)
## Try a number of neighbours between 10 and 15
(nn <- NNCH(xy, k = 10:15))
opar <- par(mfrow=c(1,2))
(i <- NNCH.k.area(nn, percent=100))
plot(i)
(i <- NNCH.der.k.area(nn, percent=100))
plot(i)
par(opar)
## 12 neighbours seem to be a good choice
nn <- NNCH(xy, k = 12)
plot(nn)
Run the code above in your browser using DataLab