# NOT RUN {
data(vipers)
data(grid)
plot(grid, cex=0.5, col='lightgrey', asp=1,
main="Vipers data", xlab="Longitude", ylab="Latitude")
points(vipers[,1:2], pch=vipers$lin)
legend(1, 38, legend=c("West", "South", "East"), pch=1:3, title="Lineages")
# }
Run the code above in your browser using DataLab