By default, addKernelDensity() function pools all observations:
snowMap(add.title = FALSE)
addKernelDensity()
But you can also view the kernel densities of individual pump neighborhoods. Either by selecting a subset of pump neighborhood or by specifying which pumps should be considered (i.e., defining the “population” of pump neighborhoods).
snowMap(add.title = FALSE)
addKernelDensity(pump.subset = c(6, 8))
snowMap(add.title = FALSE)
addKernelDensity(pump.select = c(6, 8))