The seawaveQ R package provides functionality and help to fit a parametric regression model, SEAWAVE-Q, to pesticide concentration data from stream-water samples to assess trends. The model incorporates the strong seasonality and high degree of censoring common in pesticide data, and users can incorporate numerous ancillary variables such as streamflow anomalies. The model is fitted to pesticide data using maximum likelihood methods for censored data and is robust in terms of pesticide, stream location, and degree of censoring of the concentration data. This R package standardizes this methodology for trend analysis, documents the code, and provides help and tutorial information.
In previous investigations, the SEAWAVE-Q model assumed a linear trend across the period analyzed. For short trend periods, this assumption of a linear trend is adequate. However, as the period of record analyzed becomes longer, the assumption of linearity is problematic because of changes in pesticide regulation and use, some of which can be abrupt. In this update to the model, a restricted cubic spline option was added for long trend periods. This option allows for more flexibility in the time component of the model. Bootstrap functionality is included to determine statistical significance. Model results with the new restricted cubic spline option are compared to the linear trend option for two pesticide-site combinations.
|Title||seawaveQ—An R package providing a model and utilities for analyzing trends in chemical concentrations in streams with a seasonal wave (seawave) and adjustment for streamflow (Q) and other ancillary variables, version 2.0.0|
|Authors||Karen R. Ryberg, Benjamin C. York|
|Publication Subtype||USGS Numbered Series|
|Series Title||Open-File Report|
|Record Source||USGS Publications Warehouse|
|USGS Organization||WMA - Earth System Processes Division|