bety_sites.R

library(PEcAn.DB)
library(ggplot2)
dbparms <- list(driver = "PostgreSQL", dbname = "bety", user = "bety", host = "localhost", password = "bety")
con <- db.open(dbparms)
views <- db.query("select * from traits_and_yields_view_private", con = con)
db.close(con)

world <- map_data("world")
ggplot() +  
  geom_jitter(data = views, 
              aes(lon, lat), alpha = 0.1, position = position_jitter(width = 1, height = 1)) +
  geom_polygon(data = world, aes(long, lat, group = group), color = "grey80", alpha = 0) + 
  theme_bw() + scale_color_brewer(palette="Dark2")