joinTreeData: joinTreeData: compiles tree data

View source: R/joinTreeData.R

joinTreeDataR Documentation

joinTreeData: compiles tree data

Description

This function combines tree location and visit data for measurements that have only 1 record per visit. Must run importData first. Abandoned plots are excluded from function.

Usage

joinTreeData(
  park = "all",
  from = 2006,
  to = as.numeric(format(Sys.Date(), "%Y")),
  QAQC = FALSE,
  locType = c("VS", "all"),
  panels = 1:4,
  status = c("all", "active", "live", "dead"),
  speciesType = c("all", "native", "exotic", "invasive"),
  canopyPosition = c("all", "canopy"),
  dist_m = NA,
  eventType = c("complete", "all"),
  output = "short"
)

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 2006 to current year

to

Year to stop analysis, ranging from 2006 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

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.

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.

status

Filter by live, dead, or all. Acceptable options are:

"all"

Default. Includes all trees with any status, including excluded or missing.

"active"

Includes all trees with an active monitoring status, including "DF".

"live"

live trees only

"dead"

dead trees only

speciesType

Allows you to filter on native, exotic or include all species.

"all"

Default. Returns all species.

"native"

Returns native species only

"exotic"

Returns exotic species only

"invasive"

Returns species on the Indicator Invasive List

canopyPosition

Allows you to filter on tree crown class

"all"

Returns all canopy positions

"canopy"

Returns only dominant, codominant, and intermediate crown classes. Since only live trees are assigned crown classes, this also only returns live trees.

dist_m

Filter trees by a distance that is less than or equal to the specified distance in meters of the tree to the center of the plot. If no distance is specified, then all trees will be selected. For example, to select an area of trees that is 100 square meters in area, use a distance of 5.64m.

eventType

Allows you to include only complete sampling events or all sampling events

"complete"

Default. Only include sampling events for a plot that are complete.

"all

Include all plot events with a record in tblCOMN.Event, including plots missing most of the data associated with that event (eg ACAD-029.2010). This feature is currently hard-coded in the function.

output

Allows you to return all columns or just the most important columns for analysis. Valid inputs are "short" and "verbose".

Value

returns a dataframe with plot-level and visit-level tree data. Returns records for all specified plots and events, even if no trees meet the specified arguments (eg dead or exotic trees), although all associated data (eg TagCode, DBH), will be NA for those plot/events. ScientificName will be "None present".

Examples

## Not run: 
importData()
# compile tree data in all parks for live trees only in cycle 3, excluding QAQC visits
live_trees <- joinTreeData(status = 'live', from = 2014, to = 2017)

# compile ACAD trees within 100m^2 circle in cycle 3
ACAD_100m <- joinTreeData(park = 'ACAD', from = 2014, to = 2017, dist_m = 5.64)

# compile dead trees in MABI for cycle 3
MABI_dead <- joinTreeData(park = 'MABI', from = 2014, to = 2017, status = 'dead')

# compile exotic trees in MIMA in all years
MIMA_exotic <- joinTreeData(park = 'MIMA', speciesType = 'exotic')

# compile all visits in ROVA for 2019, including QAQC visits
ROVA_trees <- joinTreeData(park = "ROVA", from = 2019, to = 2019, QAQC = TRUE)

## End(Not run)

KateMMiller/forestNETN documentation built on Nov. 11, 2024, 6:07 a.m.