[R-sig-Geo] spatialpoints: each dot represents 100 individuals?

rubenfcasal rubenfcasal at gmail.com
Wed Apr 27 23:13:42 CEST 2016

Alternatively, you might also consider data binning (implemented in 
several packages: KernSmooth, ks, sm, npsp ,...). This technique is 
commonly used in nonparametric statistics to reduce the computational 
time (see e.g. Wand, M. P. (1994), Fast Computation of Multivariate 
Kernel Estimators, Journal of Computational and Graphical Statistics, 3, 

For instance, using the npsp package (maintained by me...), you could do 
something like this:


bin <- binning(earthquakes[, c("lon", "lat")], nbin = c(50,50))

# ‘bin$binw’ will contain the binning weights (aggregations) at 
locations ‘coords(bin)’


Additionally, you could estimate (nonparametrically) the spatial density:

h <- h.cv(bin, ncv = 2)$h

den <- np.den(bin, h = h)

plot(den, log = FALSE, main = 'Estimated density')

Best regards,


El 25/04/2016 a las 13:35, Juta Kawalerowicz escribió:
> Hi,
> I have a dataset with couple of million of points (individuals) and
> would like to do some mapping (I have the coordinates of each point)
> but given the number of observation I think it may be usuful to plot
> dots which represent 100 individuals (of a given group). Does anyone
> know a good way to aggregate up spatialpoints? Any suggestions would
> be much appreciated!
> Best wishes,
> Juta
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

	[[alternative HTML version deleted]]

More information about the R-sig-Geo mailing list