Marcus W. Beck, This email address is being protected from spambots. You need JavaScript enabled to view it.
Overview
The System Wide Monitoring Program (SWMP) was implemented by the National Estuarine Research Reserve System (NERRS) in 1995 to provide continuous monitoring data at over 140 continuous monitoring stations in 28 estuaries across the United States. SWMPr (pronounced "swamper") is an R package for retrieving, organizing, and analyzing estuary monitoring data from SWMP. Please cite the package as follows:
Beck MW. 2015. SWMPr: An R package for the National Estuarine Research Reserve System. Version 2.0.0. https://github.com/fawda123/SWMPr
All data obtained from the CDMO should be cited using the format:
National Estuarine Research Reserve System (NERRS). 2012. System-wide Monitoring Program. Data accessed from the NOAA NERRS Centralized Data Management Office website: http://cdmo.baruch.sc.edu/
Installing and loading the package
The package can be installed from CRAN (the archived package network for R) then loaded as follows:
install.packages('SWMPr')
library(SWMPr)
The development (unstable) version of the package can be installed from GitHub:
install.packages('devtools')
library('devtools')
install_github('fawda123/SWMPr')
library(SWMPr)
Using the package
A manuscript describing use of the package in detail is in preparation. For the time being, please see an old README file from the GitHub page. Below is a brief summary of the available functions. The help documentation for each function can be viewed in R after installing and loading the package (e.g., run '?all_params
' at the command line).
Retrieve
all_params |
Retrieve any number of records starting with the most recent at a given station, all parameters. Wrapper to exportAllParamsXMLNew function on web services. |
all_params_dtrng |
Retrieve records of all parameters within a given date range for a station. Optional argument for a single parameter. Wrapper to exportAllParamsDateRangeXMLNew . |
import_local |
Import files from a local path. The files must be in a specific format, specifically those returned from the CDMO using the zip downloads option for a reserve. |
import_remote |
Import SWMP site data from a remote independent server. These files have been downloaded from CDMO, processed using functions in this package, and uploaded to an Amazon server for quicker import into R. |
single_param |
Retrieve any number of records for a single parameter starting with the most recent at a given station. Wrapper to exportSingleParamXMLNew function on web services. |
Organize
comb.swmpr |
Combines swmpr objects to a common time series using setstep, such as combining the weather, nutrients, and water quality data for a single station. Only different data types can be combined. |
qaqc.swmpr |
Remove QAQC columns and remove data based on QAQC flag values for a swmpr object. Only applies if QAQC columns are present. |
qaqcchk.swmpr |
View a summary of the number of observations in a swmpr object that are assigned to different QAQC flags used by CDMO. The output is used to inform further processing but is not used explicitly. |
rem_reps.swmpr |
Remove replicate nutrient data that occur on the same day. The default is to average replicates. |
setstep.swmpr |
Format data from a swmpr object to a continuous time series at a given timestep. The function is used in comb.swmpr and can also be used with individual stations. |
subset.swmpr |
Subset by dates and/or columns for a swmpr object. This is a method passed to the generic `subset’ function provided in the base package. |
Analyze
aggreswmp.swmpr |
Aggregate swmpr objects for different time periods - years, quarters, months, weeks, days, or hours. Aggregation function is user-supplied but defaults to mean. |
aggremetab.swmpr |
Aggregate metabolism data from a swmpr object. This is primarily used within plot_metab but may be useful for simple summaries of raw daily data. |
ecometab.swmpr |
Estimate ecosystem metabolism for a combined water quality and weather dataset using the open-water method. |
decomp.swmpr |
Decompose a swmpr time series into trend, seasonal, and residual components. This is a simple wrapper to decompose . Decomposition of monthly or daily trends is possible. |
decomp_cj.swmpr |
Decompose a swmpr time series into grandmean, annual, seasonal, and events components. This is a simple wrapper to decompTs in the wq package. Only monthly decomposition is possible. |
hist.swmpr |
Plot a histogram for a swmpr object. |
lines.swmpr |
Add lines to an existing swmpr plot. |
na.approx.swmpr |
Linearly interpolate missing data (NA values) in a swmpr object. The maximum gap size that is interpolated is defined as a maximum number of records with missing data. |
plot.swmpr |
Plot a univariate time series for a swmpr object. The parameter name must be specified. |
plot_metab |
Plot ecosystem metabolism estimates after running ecometab on a swmpr object. |
plot_summary |
Create summary plots of seasonal/annual trends and anomalies for a water quality or weather parameter. |
smoother.swmpr |
Smooth swmpr objects with a moving window average. Window size and sides can be specified, passed to filter . |
Miscellaneous
calcKL |
Estimate the reaeration coefficient for air-sea gas exchange. This is only used within the ecometab function. |
map_reserve |
Create a map of all stations in a reserve using the ggmap package. |
metab_day |
Identify the metabolic day for each approximate 24 period in an hourly time series. This is only used within the ecometab function. |
param_names |
Returns column names as a list for the parameter type(s) (nutrients, weather, or water quality). Includes QAQC columns with ‘f_’ prefix. Used internally in other functions. |
parser |
Parses html returned from CDMO web services, used internally in retrieval functions. |
site_codes |
Metadata for all stations, wrapper to exportStationCodesXMLNew function on web services. |
site_codes_ind |
Metadata for all stations at a single site, wrapper to NERRFilterStationCodesXMLNew function on web services. |
swmpr |
Creates object of swmpr class, used internally in retrieval functions. |
time_vec |
Converts time vectors to POSIX objects with correct time zone for a site/station, used internally in retrieval functions. |