Retrieve Data from MacLeish Field Station

Download data from the Ada and Archibald MacLeish Field Station in Whately, MA. The Ada and Archibald MacLeish Field Station is a 260-acre patchwork of forest and farmland located in West Whately, MA that provides opportunities for faculty and students to pursue environmental research, outdoor education, and low-impact recreation (see < http://www.smith.edu/ceeds/macleish.php> for more information). This package contains weather data over several years, and spatial data on various man-made and natural structures.


The Ada and Archibald MacLeish Field Station is a 260-acre patchwork of forest and farmland located in West Whately, MA that provides opportunities for faculty and students to pursue environmental research, outdoor education, and low-impact recreation. More information can be found at (http://www.smith.edu/ceeds/macleish.php)

This R package allows you to download and process weather data using the ETL framework from the MacLeish Field Station. It also contains shapefiles for contextualizing spatial information.

devtools::install_github("beanumber/macleish")

Weather data from 2015 is available immediately from both the Whately and Orchard weather stations.

library(macleish)
glimpse(whately_2015)
## Observations: 52,560
## Variables: 8
## $ when            (time) 2015-01-01 00:00:00, 2015-01-01 00:10:00, 201...
## $ temperature     (dbl) -9.32, -9.46, -9.44, -9.30, -9.32, -9.34, -9.3...
## $ wind_speed      (dbl) 1.399, 1.506, 1.620, 1.141, 1.223, 1.090, 1.16...
## $ wind_dir        (dbl) 225.4, 248.2, 258.3, 243.8, 238.4, 241.7, 242....
## $ rel_humidity    (dbl) 54.55, 55.38, 56.18, 56.41, 56.87, 57.25, 57.7...
## $ pressure        (int) 985, 985, 985, 985, 984, 984, 984, 984, 984, 9...
## $ solar_radiation (dbl) 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
## $ rainfall        (int) 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
glimpse(orchard_2015)
## Observations: 52,547
## Variables: 9
## $ when         (time) 2015-01-01 00:00:00, 2015-01-01 00:10:00, 2015-0...
## $ temperature  (dbl) -9.62, -9.45, -9.52, -9.68, -9.56, -9.59, -9.51, ...
## $ wind_speed   (dbl) 0.490, 0.517, 0.576, 0.522, 0.296, 0.382, 0.497, ...
## $ wind_dir     (dbl) 216.3, 209.6, 215.5, 199.3, 186.4, 183.9, 195.6, ...
## $ rel_humidity (dbl) 61.97, 59.08, 58.69, 58.59, 57.32, 57.58, 56.97, ...
## $ pressure     (dbl) 1020, 1019, 1019, 1019, 1019, 1019, 1019, 1018, 1...
## $ par_density  (dbl) 0.548, 0.548, 0.548, 0.548, 0.548, 0.548, 0.548, ...
## $ par_total    (dbl) 0.088, 0.088, 0.088, 0.088, 0.088, 0.088, 0.088, ...
## $ rainfall     (int) 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...

Weather readings are logged every 10 minutes. Current and historical (dating back to 1/3/2012 for whately and 6/27/2014 for orchard) meteorological readings are available through the ETL framework. Please see the documentation for that package for more information about how this works.

macleish <- etl("macleish")
macleish %>%
  etl_update()
whately <- macleish %>%
  tbl("whately")
whately %>%
  mutate(the_year = strftime('%Y', when)) %>%
  group_by(the_year) %>%
  summarize(N = n(), begin = min(when), end = max(when), avg_temp = mean(temperature))
orchard <- macleish %>%
  tbl("orchard")
orchard %>%
  mutate(the_year = strftime('%Y', when)) %>%
  group_by(the_year) %>%
  summarize(N = n(), begin = min(when), end = max(when), avg_temp = mean(temperature))
daily <- whately %>%
  mutate(the_date = date(when)) %>%
  group_by(the_date) %>%
  summarize(N = n(), avgTemp = mean(temperature)) %>%
  collect()
library(ggplot2)
ggplot(data = daily, aes(x = as.Date(the_date), y = avgTemp)) +
  geom_point() + geom_smooth()

Spatial data is available through the macleish_layers data object.

class(macleish_layers)
## [1] "list"
names(macleish_layers)
##  [1] "landmarks"         "forests"           "streams"          
##  [4] "challenge_courses" "buildings"         "wetlands"         
##  [7] "boundary"          "research"          "soil"             
## [10] "trails"
library(leaflet)
leaflet() %>%
  addTiles() %>%
  addPolygons(data = macleish_layers[["boundary"]], 
              weight = 1, fillOpacity = 0.1) %>%
  addPolygons(data = macleish_layers[["buildings"]], 
              weight = 1, popup = ~ name) %>%
  addPolylines(data = macleish_layers[["trails"]], 
               weight = 1, color = "brown",
               popup = ~ name) %>%
  addPolylines(data = macleish_layers[["streams"]], 
               weight = 2) %>%
  addMarkers(data = subset(macleish_layers[["landmarks"]], grepl("Met", Label)), 
             popup = ~Label)

News

macleish 0.3.0

  • Added a NEWS.md file to track changes to the package.

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.

install.packages("macleish")

0.3.0 by Ben Baumer, a year ago


Browse source code at https://github.com/cran/macleish


Authors: Ben Baumer [aut, cre], Rose Goueth [aut], Wencong Li [aut], Weijia Zhang [aut], Nicholas Li [aut]


Documentation:   PDF Manual  


CC0 license


Imports dplyr, DBI, lubridate, readr

Depends on etl

Suggests ggplot2, RSQLite, broom, knitr, ggmap, leaflet, rmarkdown, clifro, testthat

Enhances sp, rgdal, rgeos


Suggested by mdsr.


See at CRAN