hadley wickham wrote:
Here's a quick start using ggplot2 - one trick to make the map look
nicer is to use a bigger range when calling map, and then adjust the
limits of the plot to zoom in on the region of interest.
library(maps)
library(ggplot2)
langren <-
read.csv("http://euclid.psych.yorku.ca/SCS/Gallery/Private/langrens.csv";,
header=TRUE)
bbox <- c( 38.186, -9.184,
43.692, 28.674 )
bbox <- matrix(bbox, 2, 2, byrow=TRUE)
borders <- as.data.frame(map("world", plot = F,
xlim = expand_range(bbox[,2], 0.2),
ylim = expand_range(bbox[,1], 0.2))[c("x", "y")])
data(world.cities)
cities <- subset(world.cities,
name %in% c("Rome", "Toledo") & country.etc %in% c("Spain", "Italy"))
ggplot(langren, aes(long, lat)) +
geom_path(aes(x, y), borders, colour = "grey60") +
geom_point(y = 40) +
geom_text(aes(label = name), y = 40.1, angle = 90, hjust = 0, size = 3) +
geom_point(data = cities, colour = "red", size = 2) +
coord_cartesian(xlim=bbox[,2], ylim=bbox[,1])
Thanks very much, Hadley. However, it doesn't work for me, using
ggplot2_0.7, even when adding
langren$lat <- 39.68
> ggplot(langren, aes(long, lat)) +
+ geom_path(aes(x, y), borders, colour = "grey60") +
+ geom_point(y = 40) +
+ geom_text(aes(label = name), y = 40.1, angle = 90, hjust = 0, size =
3) +
+ geom_point(data = cities, colour = "red", size = 2) +
+ coord_cartesian(xlim=bbox[,2], ylim=bbox[,1])
Error in get("new", env = CoordCartesian, inherits =
TRUE)(CoordCartesian, :
unused argument(s) (xlim = c(-9.184, 28.674), ylim = c(38.186, 43.692))
Eliminating the last line (coord_caertesian) gives a different error,
and just the first gives
> ggplot(langren, aes(long, lat))
Error: No layers in plot
> sessionInfo()
R version 2.8.0 (2008-10-20)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
[1] splines grid stats graphics grDevices utils
datasets methods
[9] base
other attached packages:
[1] ggplot2_0.7 MASS_7.2-44 RColorBrewer_1.0-2
proto_0.3-8
[5] reshape_0.8.2 plyr_0.1.3 rimage_0.5-7
foreign_0.8-29
[9] maps_2.0-40
...
--
Michael Friendly Email: friendly AT yorku DOT ca
Professor, Psychology Dept.
York University Voice: 416 736-5115 x66249 Fax: 416 736-5814
4700 Keele Street http://www.math.yorku.ca/SCS/friendly.html
Toronto, ONT M3J 1P3 CANADA
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
