| library(neotoma) | |
| all_ds <- get_dataset(datasettype = "pollen", gpid = "United States") | |
| as_sites <- get_site(all_ds) | |
| site_elevations <- as_sites$elev | |
| # It's also within the datasets: | |
| site_elevations <- sapply(all_ds, function(x)x$site.data$elev) | |
| hist(site_elevations, breaks = seq(0, 4000, by = 100)) |