##
data(NYdata)
head(NYdata)
# plots in NY map
NYsite<-unique(cbind(NYdata[,1:3]))
head(NYsite)
# map
library(maps)
map(database="state",regions="new york")
points(NYsite[,2:3],pch=19)
##Run the code above in your browser using DataLab