
Last chance! 50% off unlimited learning
Sale ends in
qmplot
is the ggmap equivalent to the ggplot2
function qplot and allows for the quick plotting of maps
with data/models/etc. qmplot is still experimental.qmplot(x, y, ..., data, zoom, source = "stamen",
extent = "device", legend = "right", padding = 0.02,
darken = c(0, "black"), mapcolor = "color",
facets = NULL, margins = FALSE, geom = "auto",
stat = list(NULL), position = list(NULL),
xlim = c(NA, NA), ylim = c(NA, NA), main = NULL,
f = 0.05, xlab = deparse(substitute(x)),
ylab = deparse(substitute(y)))
get_map
get_map
facet_wrap
or facet_grid
depending on whether the formula is one sided or
two-sidedqmplot(lon, lat, data = crime)
# only violent crimes
violent_crimes <- subset(crime,
offense != 'auto theft' &
offense != 'theft' &
offense != 'burglary'
)
# rank violent crimes
violent_crimes$offense <-
factor(violent_crimes$offense,
levels = c('robbery', 'aggravated assault',
'rape', 'murder')
)
# restrict to downtown
violent_crimes <- subset(violent_crimes,
-95.39681 <= lon & lon <= -95.34188 &
29.73631 <= lat & lat <= 29.78400
)
theme_set(theme_bw())
qmplot(lon, lat, data = violent_crimes, colour = offense, darken = .5,
size = I(3.5), alpha = I(.6), legend = 'topleft')
qmplot(lon, lat, data = violent_crimes, geom = c('point','density2d'))
qmplot(lon, lat, data = violent_crimes) + facet_wrap(~ offense)
qmplot(lon, lat, data = violent_crimes, extent = 'panel') + facet_wrap(~ offense)
qmplot(lon, lat, data = violent_crimes, extent = 'panel', colour = offense) +
facet_wrap(~ month)
# doesn't quite work yet....
qmplot(long, lat, xend = long + delta_long,
yend = lat + delta_lat, data = seals, geom = 'segment')
library(scales)
library(grid)
options('device')$device(width = 4.98, height = 5.97)
qmplot(lon, lat, data = wind, size = I(.5), alpha = I(.5)) +
ggtitle('NOAA Wind Report Sites')
# thin down data set...
s <- seq(1, 227, 8)
thinwind <- subset(wind,
lon %in% unique(wind$lon)[s] &
lat %in% unique(wind$lat)[s]
)
# for some reason adding arrows to the following plot bugs
theme_set(theme_bw(18))
options('device')$device(width = 6.13, height = 7.04)
qmplot(lon, lat, data = thinwind, geom = 'tile', fill = spd, alpha = spd,
legend = 'bottomleft') +
geom_leg(aes(xend = lon + delta_lon, yend = lat + delta_lat)) +
scale_fill_gradient2('Wind Speed\nand\nDirection',
low = 'green', mid = muted('green'), high = 'red') +
scale_alpha('Wind Speed\nand\nDirection', range = c(.1, .75)) +
guides(fill = guide_legend(), alpha = guide_legend())
Run the code above in your browser using DataLab