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

Description Usage Arguments Value Examples

View source: R/joinSoilLabData.R

Description

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.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
joinSoilLabData(
  park = "all",
  from = 2007,
  to = 2021,
  QAQC = FALSE,
  panels = 1:4,
  locType = c("VS", "all"),
  layer = c("all", "O", "A"),
  ...
)

Arguments

park

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

"all"

Includes all parks in the network

"ACAD"

Acadia NP only

"MABI"

Marsh-Billings-Rockefeller NHP only

"MIMA"

Minute Man NHP only

"MORR"

Morristown NHP only

"ROVA"

Roosevelt-Vanderbilt NHS only

"SAGA"

Saint-Gaudens NHS only

"SARA"

Saratoga NHP only

"WEFA"

Weir Farm NHS only

from

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

to

Year to stop analysis, ranging from 2007 to current year

QAQC

Allows you to remove or include QAQC events.

FALSE

Default. Only returns visits that are not QAQC visits

TRUE

Returns all visits, including QAQC visits

panels

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.

locType

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

"VS"

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

"all"

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

layer

Allows you to filter on soil horizons

"all"

Default. Includes O and A horizons.

"O"

Return only samples from the O horizon.

"A"

Return only samples from the A horizon.

Value

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.

Examples

1
2
3
4
5
6
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')

KateMMiller/forestNETN documentation built on May 30, 2021, 6:39 a.m.