## Not run:
# data(wheat)
# wheat$lat1 <- 69 - wheat$lat
# wheat$r <- factor(wheat$lat1)
# wheat$c <- factor(wheat$lon)
# wheat_sp <- wheat
# coordinates(wheat_sp) <- c("lon", "lat1")
# wheat_spg <- wheat_sp
# gridded(wheat_spg) <- TRUE
# wheat_spl <- as(wheat_spg, "SpatialPolygons")
# df <- as(wheat_spg, "data.frame")
# row.names(df) <- sapply(slot(wheat_spl, "polygons"),
# function(x) slot(x, "ID"))
# wheat <- SpatialPolygonsDataFrame(wheat_spl, data=df)
# ## End(Not run)
require(maptools)
wheat <- readShapeSpatial(system.file("etc/shapes/wheat.shp",
package="spdep")[1])
Run the code above in your browser using DataLab