R package for Tampa Bay Estuary Program functions. Please see the vignettes for a full description.
The package can be installed from r-universe. The source code is available on the tbep-tech GitHub group web page: https://github.com/tbep-tech/tbeptools. Note that tbeptools only needs to be installed once, but it needs to be loaded every new R session (i.e., library(tbeptools)
).
# enable repos
options(repos = c(
tbeptech = 'https://tbep-tech.r-universe.dev',
CRAN = 'https://cloud.r-project.org'))
# install tbeptools
install.packages('tbeptools')
# load tbeptools
library(tbeptools)
After the package is loaded, you can view the help files for each function by typing a question mark followed by the function name, e.g., ?read_importwq
, on the console. The help files provide a brief description of what each function does and the required arguments that are needed to run the function.
The vignettes are organized by topic and are an excellent place to start for understanding how to use the package. Currently, there are six vignettes available for tbeptools:
- Water Quality Data: Overview of functions for working with water quality data and the water quality report card
- Tampa Bay Nekton Index: Overview of functions to import, analyze, and plot results for the Tampa Bay Nekton Index
- Tampa Bay Benthic Index: Overview of functions to import data for Tampa Bay Benthic Index, under development
- Tidal Creeks Assessment: Overview of functions to import, analyze, and plot results for the assessment of tidal creeks in southwest Florida
- Seagrass Transect Data: Overview of functions to import, analyze, and plot results for the seagrass transect data collected in Tampa Bay
- Habitat Master Plan: Overview of functions to analyze and create a report card for the Tampa Bay Habitat Master Plan 2020 update
- Fecal Indicator Bacteria: Overview of functions to import, analyze, and plot results for Fecal Indicator Bacteria (FIB)
The core functions in tbeptools are in three categories based on mode of use. Each function is named using a prefix for the mode of use, followed by what the function does. The prefixes are:
-
read
: Import current data from the main site. -
anlz
: Analyze or summarize the imported data. -
show
: Create a plot of the analyzed data.
The functions can be easily found in RStudio after loading the package and typing the prefix at the command line. An autofill dialog box will pop up showing all functions that apply for the prefix. This eliminates the need for searching for individual functions if all you know is the category of function you need (e.g., read
, anlz
, or show
).
Each function also includes a semi-descriptive suffix that generally describes what category it applies to (e.g, water quality, seagrass) and what it does (e.g., imports, formats). These follow a loose convention that attempts to strike a balance between description and brevity. The optimal balance is often hard to achieve. To aid in understanding, we provide a brief description of suffixes that are used more than once.
Suffix descriptions:
attain
: Analyze functions that summarize data relative to attainment categories specific to bay segmentsave
,med
: Analyze functions that summarize data into averages or mediansbenthic
: Applies to benthic monitoring data used for the Tampa Bay Benthic Indexentero
: Applies to Enterococcus, a specific type of Fecal Indicator Bacteria, whenfib
functions are not applicablefib
: Applies generally to Fecal Indicator Bacteriafim
: Applies to data from the Fisheries Independent Monitoring program used for the Tampa Bay Nekton Indexform
: An intermediate function for formatting imported data for downstream analysishmp
: Functions that work with Habitat Master Plan dataimport
: A function used to import data from a source external to the packageindic
: A function that analyzes or plots individual tidal creek indicator values, as opposed to integrated creek scoresiwr
: Functions or data that apply to the Impaired Waters Rule (IWR) data maintained by the Florida Department of Environmental Protection used as source data for the tidal creek functionsmatrix
: A plotting function that creates a report card style matrixmet
: A function that analyses or plots individual metrics for integrated indices, e.g., TBBI, TBNIphyto
: Applies to phytoplankton data from the Hillsborough County Environmental Protection Commissionplotly
: A plotting function that returns an interactive plotly objectscr
: A function that analyses or plots summary scores for integrated indices, e.g., TBBI, TBNIseg
,site
: Functions that analyze or plot results relative to bay segments or individual monitoring sitestbbi
: Applies to the Tampa Bay Benthic Index (TBBI)tbni
: Applies to the Tampa Bay Nekton Index (TBNI)tdlcrk
: Applies to tidal creekstransect
: Applies to seagrass transect datawq
: Applies to water quality
The function reference page can also be viewed for a complete list of functions organized by category, a description of what they do, and links to the help files.
The following example demonstrates use of a subset of the functions for water quality data to read a file from the Hillsborough County Environmental Protection Commission long-term monitoring dataset (available from https://www.tampabay.wateratlas.usf.edu/), analyze monthly and annual averages by major bay segments of Tampa Bay, and plot an annual time series for one of the bay segments.
# load the package
library(tbeptools)
# read current data
wqdat <- read_importwq(xlsx = "wqdata.xlsx", download_latest = TRUE)
wqdat
## # A tibble: 26,611 x 22
## bay_segment epchc_station SampleTime yr mo
## <chr> <dbl> <dttm> <dbl> <dbl>
## 1 HB 6 2021-06-08 10:59:00 2021 6
## 2 HB 7 2021-06-08 11:13:00 2021 6
## 3 HB 8 2021-06-08 14:15:00 2021 6
## 4 MTB 9 2021-06-08 13:14:00 2021 6
## 5 MTB 11 2021-06-08 11:30:00 2021 6
## # ... with 26,606 more rows, and 17 more variables:
## # Latitude <dbl>, Longitude <dbl>, Total_Depth_m <dbl>,
## # Sample_Depth_m <dbl>, tn <dbl>, tn_q <chr>, sd_m <dbl>,
## # sd_raw_m <dbl>, sd_q <chr>, chla <dbl>, chla_q <chr>,
## # Sal_Top_ppth <dbl>, Sal_Mid_ppth <dbl>,
## # Sal_Bottom_ppth <dbl>, Temp_Water_Top_degC <dbl>,
## # Temp_Water_Mid_degC <dbl>, ...
# analyze monthly and annual means by bay segment
avedat <- anlz_avedat(wqdat)
avedat
## $ann
## # A tibble: 584 x 4
## yr bay_segment var val
## <dbl> <chr> <chr> <dbl>
## 1 1974 HB mean_chla 22.4
## 2 1974 LTB mean_chla 4.24
## 3 1974 MTB mean_chla 9.66
## 4 1974 OTB mean_chla 10.2
## 5 1975 HB mean_chla 27.9
## # ... with 579 more rows
##
## $mos
## # A tibble: 4,484 x 5
## bay_segment yr mo var val
## <chr> <dbl> <dbl> <chr> <dbl>
## 1 HB 1974 1 mean_chla 36.2
## 2 LTB 1974 1 mean_chla 1.75
## 3 MTB 1974 1 mean_chla 11.5
## 4 OTB 1974 1 mean_chla 4.4
## 5 HB 1974 2 mean_chla 42.4
## # ... with 4,479 more rows
# show annual time series of chlorophyll for Hillsborough bay segment
show_thrplot(wqdat, bay_segment = "HB", yrrng = c(1975, 2020))
Functions in tbeptools
also support the creation of content for interactive, online dashboards that can facilitate more informed decisions without requiring an intimate understanding of the R programming language or the methods for analysis. These dashboards include assessments for water quality, seagrasses, nekton communities, and tidal creeks.
Please report any issues and suggestions on the issues link for the repository. A guide to posting issues can be found here.
Please view our contributing guidelines for any changes or pull requests.