joinRegenData: joinRegenData: summarizes seedling and sapling data

In KateMMiller/forestNETN: Package for importing, joining and querying NETN forest database tables

joinRegenData: summarizes seedling and sapling data

Description

This function combines seedling and sapling data collected in microplots, and calculates the stocking index. Each row represents a species observed per visit. If no seedlings or saplings were observed, function returns "None present" for ScientificName and 0 for densities. Permanently Missing species tallies for species or whole plots have NAs. These are rare, have always been low canopy species initially recorded as shrubs, and mostly occur in data <2011. Note that the stocking index only includes saplings < 2.5cm DBH, but the sapling density returned is all saplings > 1cm and <10cm DBH. For the few plots with > 10 saplings of a given species in a microplot, their counts are included in the stocking index only if the average DBH of saplings measured is <=2.5 cm. This may underestimate the stocking index for those plots, but their index values are still way higher than most plots. Must run importData() or importCSV() first.

Usage

joinRegenData(
  park = "all",
  from = 2006,
  to = as.numeric(format(Sys.Date(), "%Y")),
  QAQC = FALSE,
  panels = 1:4,
  locType = c("VS", "all"),
  eventType = c("complete", "all"),
  speciesType = c("all", "native", "exotic", "invasive"),
  canopyForm = c("all", "canopy"),
  numMicros = 3,
  units = c("micro", "sq.m", "ha", "acres"),
  ...
)

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
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.
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.
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
canopyForm	Allows you to filter on species growth form "all" Default. Returns all species, including low canopy species. "canopy" Returns canopy-forming species only
numMicros	Allows you to select 1, 2, or 3 microplots of data to summarize
units	Calculates seedling and sapling densities based on different units. "micro" Default. Returns seedling and sapling densities per microplot. "sq.m" Returns seedling and sapling densities per square meter "ha" Returns seedling and sapling densities per hectare "acres" Returns densities per acre
...	Other arguments passed to function.

Value

returns a dataframe with seedling and sapling densities, and stocking index metrics for each species observed per visit.

Examples

## Not run: 
importCSV('./forest_csvs/')
# compile seedling and sapling data for all parks and all species in cycle 3
regen_data <- joinRegenData(canopyForm = 'all', from = 2014, to = 2017)

# compile regen data for only canopy-forming (default) and native species in SAGA for all years
SAGA_regen <- joinRegenData(park = 'SAGA', speciesType = 'native')

# compile only 1 microplot of data for ACAD native canopy-forming species for all but first year
ACAD_regen_m1 <- joinRegenData(park = 'ACAD', speciesType = 'native', numMicros = 1, from = 2007)

## End(Not run)

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