hexagons(x, res = 100, ...)

require(sp)
  data(meuse)
    coordinates(meuse) <- ~x+y 

hex.polys <- hexagons(meuse, res=100)   
  plot(hex.polys)
    plot(meuse,pch=20,add=TRUE)

# Points intersecting hexagons  
hex.pts <- na.omit(over(meuse,hex.polys))
(hex.pts <- data.frame(PTID=rownames(hex.pts), hex.pts))