# NOT RUN {
Npoints = 20
myPoints = SpatialPointsDataFrame(
cbind(
runif(Npoints, -15000000, 15000000),
runif(Npoints, -8000000, 8000000)),
data=data.frame(y1=c(NA, rnorm(Npoints-1)),
y2=c(sample(0:5, Npoints-1,replace=TRUE), NA)),
proj4string=moll(c(-100,0))
)
plot(myPoints)
gridlinesWrap(myPoints, lty=3, col='red')
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab