garrettgman / ggsubplot Goto Github PK

Error in get(x, envir = this, inherits = inh)(this, ...) : 
  could not find function "aesply"

Error in layout_base(data, vars, drop = drop) : 
  At least one layer must contain all variables used for facetting

var = ggplot(test)+ geom_subplot(ref=ref_box(fill="white",alpha=0.2), aes(Longitude, Latitude, group=City, subplot = geom_bar(aes(x = value, y = variable, fill=variable), stat="identity")))

library(maps)
library(ggplot2)
library(RColorBrewer)
library(ggsubplot)
## You will also need to have the mapproj package installed for coord_map()


#############################
###      Map layers       ###
#############################

# getbox by Heike Hoffman, 
# https://github.com/ggobi/paper-climate/blob/master/code/maps.r
getbox <- function (map, xlim, ylim) {
  # identify all regions involved
  small <- subset(map, (long > xlim[1]) & (long < xlim[2]) & (lat > ylim[1]) & (lat < ylim[2]))
  regions <- unique(small$region)
  small <- subset(map, region %in% regions)  

  # now shrink all nodes back to the bounding box
  small$long <- pmax(small$long, xlim[1])
  small$long <- pmin(small$long, xlim[2])
  small$lat <- pmax(small$lat, ylim[1])
  small$lat <- pmin(small$lat, ylim[2])

  # Remove slivvers
  small <- ddply(small, "group", function(df) {
    if (diff(range(df$long)) < 1e-6) return(NULL)
    if (diff(range(df$lat)) < 1e-6) return(NULL)
    df
  })

  small
}

## map layer
## adapted from map_nasa
# https://github.com/ggobi/paper-climate/blob/master/code/maps.r

# assembling data
world <- map_data("world")
states <- map_data("state")
states$group <- max(world$group) + states$group
both <- rbind(world, states)
americas <- getbox(both, xlim = c(-115, -55), ylim = c(-21.1, 36.6))
nasa$fahrenheit <- nasa$surftemp * 9 / 5 - 459.67

# building layers
map_americas <- list(
  geom_polygon(aes(long, lat, group = group), data = americas, fill = "grey70", 
    colour = "grey60", inherit.aes = FALSE, show_guide = FALSE),
  scale_x_continuous("", breaks = NULL, expand = c(0.02, 0)),
  scale_y_continuous("", breaks = NULL, expand = c(0.02, 0))) 

afghanistan <- getbox(world, c(60,75), c(28, 39))
map_afghan <- list(
  geom_polygon(aes(long, lat, group = group), data = afghanistan, 
    fill = "grey80", colour = "white", inherit.aes = FALSE, 
    show_guide = FALSE),
  scale_x_continuous("", breaks = NULL, expand = c(0.02, 0)),
  scale_y_continuous("", breaks = NULL, expand = c(0.02, 0)))  



###################################################
###      Figure 1: Afghanistan casualties       ###
###################################################

ggplot(casualties) + 
  map_afghan +
  annotate("text", label = "Kabul", x = 69.1667, y = 34.453,
    color = "white", size = 3.5) +
  annotate("text", label = "Kandahar", x = 65.7053, y = 31.4978,
    color = "white", size = 3.4) +
  geom_subplot2d(aes(lon, lat,
    subplot = geom_bar(aes(victim, ..count.., fill = victim), 
    color = rev(brewer.pal(5,"Blues"))[1], size = 1/4)),
    bins = c(15,12), ref = NULL, width = rel(0.8), height = rel(1)) + 
  ggtitle("Afghan War Diary - Casualty counts by type and location") +
  coord_map() +
  scale_fill_manual(values = rev(brewer.pal(5,"Blues"))[c(1,4,2,3)]) 


###################################################
###      Figure 2: Overplotted diamonds         ###
###################################################
ggplot(diamonds) +
  geom_subplot2d(aes(log10(carat), log10(price), 
    subplot = geom_bar(aes(color, fill = color), position = "dodge")), 
    bins = c(10, 14), y_scale = free, height.adjust = 0.8,
    width.adjust = 0.8, ref = ref_box(aes(color = length(color)))) +
  scale_color_gradient("Total\ncount", low = "grey70", high = "black") +
  ggtitle("Diamonds, carat vs. price")


###################################################
###      Figure 3: Glyphmap                     ###
###################################################
ggplot(nasa) + 
  map_americas + 
  geom_subplot(aes(long, lat, group = id, 
    subplot = geom_star(aes(r = surftemp, angle = date, 
    fill = mean(surftemp)), r.zero = FALSE))) + 
  scale_fill_gradient("Average\nTemperature (F)   ", 
    breaks = c(55,60,65,70,75,80), guide = guide_colorbar()) +
  ggtitle("Surface temperature fluctuations 1995 - 2001") +
  coord_map()


###################################################
###      Figure 3: Interactions                 ###
###################################################
ggplot(nasa) + 
  geom_subplot(aes(min(surftemp), max(surftemp), group = id, 
    subplot = geom_star(aes(r = surftemp, angle = date, 
    fill = mean(surftemp)), r.zero = FALSE)), width = 2, height = 2) + 
  scale_fill_gradient("Average\nTemperature (F)   ", 
    breaks = c(55,60,65,70,75,80), guide = guide_colorbar()) +
  scale_x_continuous("Minimum temperature") +
  scale_y_continuous("Maximum temperature") +
  ggtitle("Temperature extremes and seasonality\nby region 1995 - 2001") +
  coord_fixed()

ggplot(casualties, aes(lon, lat), size = 0) +
    geom_subplot2d(aes(lon, lat, 
      subplot = geom_bar(aes(victim, ..count.., fill = victim))), 
      breaks=list(c(60,65,70, 75), c(28, 33, 38)), ref = NULL, width = rel(0.8))

Error in layout_base(data, vars, drop = drop) : 
  At least one layer must contain all variables used for facetting

> R.version
               _                           
platform       x86_64-w64-mingw32          
arch           x86_64                      
os             mingw32                     
system         x86_64, mingw32             
status                                     
major          3                           
minor          1.0                         
year           2014                        
month          04                          
day            10                          
svn rev        65387                       
language       R                           
version.string R version 3.1.0 (2014-04-10)
nickname       Spring Dance 

> sessionInfo()
R version 3.1.0 (2014-04-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] ggsubplot_0.3.1 ggplot2_0.9.3.1 devtools_1.5   

loaded via a namespace (and not attached):
 [1] colorspace_1.2-4 digest_0.6.4     evaluate_0.5.5   grid_3.1.0       gtable_0.1.2     httr_0.3         labeling_0.2     MASS_7.3-32     
 [9] memoise_0.2.1    munsell_0.4.2    parallel_3.1.0   plyr_1.8.1       proto_0.3-10     Rcpp_0.11.1      RCurl_1.95-4.1   reshape2_1.4    
[17] scales_0.2.4     stringr_0.6.2    tools_3.1.0      whisker_0.3-2