joinSoilLabData: joinSoilLabData: compile and QC soil chemistry data by...

View source: R/joinSoilLabData.R

joinSoilLabDataR Documentation

joinSoilLabData: compile and QC soil chemistry data by horizon.


This function verifies whether O and A horizons were named corrected based on For duplicate horizons on a plot, chemistry variables are corrected using weighted averages, with sample depth as the weight. Must run importData first. Note that Earthworms are summarized in joinStandData(). Only works for complete visits and plots that haven't been abandoned. Note that data starts at 2007 because 2006 methods were pretty different.


  park = "all",
  from = 2007,
  to = as.numeric(format(Sys.Date(), "%Y")),
  panels = 1:4,
  locType = c("VS", "all"),
  layer = c("all", "O", "A")



Combine data from all parks or one or more parks at a time. Valid inputs:


Includes all parks in the network


Acadia NP only


Marsh-Billings-Rockefeller NHP only


Minute Man NHP only


Morristown NHP only


Roosevelt-Vanderbilt NHS only


Saint-Gaudens NHS only


Saratoga NHP only


Weir Farm NHS only


Year to start analysis, ranging from 2007 to current year.


Year to stop analysis, ranging from 2007 to current year


Allows you to remove or include QAQC events.


Default. Only returns visits that are not QAQC visits


Returns all visits, including QAQC visits


Allows you to select individual panels from 1 to 4. Default is all 4 panels (1:4). If more than one panel is selected, specify by c(1, 3), for example.


Allows you to only include plots that are part of the GRTS sample design or include all plots, such as deer exclosures.


Only include plots that are part of the Vital Signs GRTS sample design


Include all plots, such as plots in deer exclosures or test plots.


Allows you to filter on soil horizons


Default. Includes O and A horizons.


Return only samples from the O horizon.


Return only samples from the A horizon.


returns a dataframe containing each plot and visit with soil chemistry data for each horizon on a plot Plots that weren't sampled during a given cycle are not returned. Horizon depths are averaged across samples. Note that horizons that were combined after lab QC may be > 10 cm deep.


## Not run: 
importData() #imports using default odbc
# join only O horizon data for most recent cycle in ACAD. Note soil chemistry data are typically a year behind plot data.
soil_ACAD_O <- joinSoilLabData(park = 'ACAD', from = 2015, to = 2018, layers = 'O')

# join all park data from all layers and all years
soil_df_all <- joinSoilLabData(from = 2007, to = 2018, layers = 'all')

## End(Not run)

KateMMiller/forestNETN documentation built on April 29, 2023, 10:36 p.m.