@@ -4,7 +4,7 @@ author: "David LeBauer, Chris Black"
44date : " `r Sys.Date()`"
55output : md_document
66---
7- ``` {r setup, include=FALSE}
7+ ``` {r traits-03-mac-traits- setup, include=FALSE}
88knitr::opts_chunk$set(echo = FALSE, cache = TRUE)
99library(dplyr)
1010library(tidyr)
@@ -26,7 +26,7 @@ betykey <- "9999999999999999999999999999999999999999"
2626
2727query sites for season 2
2828
29- ``` {r}
29+ ``` {r traits-03-query-mac-sites }
3030sites <- betydb_query(
3131 table = "sites",
3232 city = "Maricopa", sitename = "~Season 2 range",
@@ -36,7 +36,7 @@ sites <- betydb_query(
3636
3737plot polygons on a map
3838
39- ``` {r}
39+ ``` {r traits-03-map-mac-polygons }
4040site_bounds <- (sites
4141 %>% rowwise()
4242 %>% do(boundaries = readWKT(text = .$geometry, id = .$id)))
@@ -48,7 +48,7 @@ leaflet() %>% addTiles() %>% addPolygons(data=site_bounds)
4848
4949## Cultivars
5050
51- ``` {r}
51+ ``` {r traits-03-mac-cultivars }
5252cultivars <- (
5353 betydb_query(
5454 table = "cultivars", limit = "none",
@@ -84,7 +84,7 @@ sites <- (sites
8484
8585First look up variables by name. Let's look for measurements related to canopy size:
8686
87- ``` {r}
87+ ``` {r traits-03-height-cover-ndvi }
8888variables <- betydb_query(
8989 table = "variables", name = "~^(NDVI|canopy_height|canopy_cover|)$",
9090 betyurl = betyurl, key = betykey, api_version = "beta")
@@ -96,7 +96,7 @@ Exercise: Why are there two variables named canopy_height, and what database fie
9696
9797Now retrieve all available measurements for each variable.
9898
99- ``` {r}
99+ ``` {r traits-03-get-variables }
100100vars_measures <- (variables
101101 %>% group_by(id, name)
102102 %>% do(traits = betydb_record( # Get full trait list by variable ID
@@ -113,7 +113,7 @@ traitdata <- left_join( # ensures we keep a blank row for any variables with no
113113
114114Add cultivar information so we can plot by ecotype
115115
116- ``` {r}
116+ ``` {r traits-03-cultivar-info }
117117traitdata <- (traitdata
118118 %>% rename(
119119 variable_name = name,
@@ -124,7 +124,7 @@ traitdata <- (traitdata
124124 %>% left_join(cultivars, by = "cultivar_id"))
125125```
126126
127- ``` {r}
127+ ``` {r traits-03-plots }
128128(ggplot(
129129 traitdata %>% filter(variable_name == "canopy_height"),
130130 aes(date, mean, group = site_id * cultivar_id))
0 commit comments