calculate geometric mean of mass of points
The maptools package has a function for computing the centroids of
polygons - install that and try help.search("centroid").
x <- read.shape(system.file("shapes/sids.shp", package="maptools")
mycentroid <- get.Pcent(x)