pos.shapes(xy, pos.shp = c("range", "full", "circle", "xy", "range_axes",
"full_axes")[1], nb.shp = 12, nr.shp = 6, nc.shp = 5, circle.r.shp)range of xy,
full extent of the plane, circle as a rosewind,
on xy values provided, range_axes on the range of xy
but on the axes, full_axes same thing but on (0.85) range of the axes.
You can also directly pass a matrix (or a data.frame)
with columns named ("x", "y").