meteval: Evaluate metabolism results

View source: R/meteval.R

metevalR Documentation

Evaluate metabolism results

Description

Evaluate metabolism results before and after weighted regression

Usage

meteval(metab_in, ...)

## S3 method for class 'metab'
meteval(metab_in, all = TRUE, ...)

Arguments

metab_in

input metab object as returned from ecometab

...

additional arguments passed to other methods

all

logical indicating if all evaluation summaries are returned or just mean, sd, and percent anomalies

Details

This function provides summary statistics of metabolism results to evaluate the effectiveness of weighted regression. These estimates are mean production, standard deviation of production, percent of production estimates that were anomalous, mean respiration, standard deviation of respiration, percent of respiration estimates that were anomalous, correlation of dissolved oxygen with tidal height changes, correlation of production with tidal height changes, and the correlation of respiration with tidal height changes. The correlation estimates are based on an average of the correlations by each month in the time series from the raw data for dissolved oxygen and the daily results for the metabolic estimates. Dissolved oxygen is correlated directly with tidal height at each time step. The metabolic estimates are correlated with the tidal height ranges during the day for production and during the night for respiration. Tidal height ranges are estimated from the raw data during each diurnal period for each metabolic day.

In general, useful results for weighted regression are those that remove the correlation of dissolved oxygen, production, and respiration with tidal changes. Similarly, the mean estimates of metabolism should not change if a long time series is evaluated, whereas the standard deviation and percent anomalous estimates should decrease.

Tables 2 and 3 in Beck et al. 2015 were created using these methods.

Value

A two-element list of summary statistics for the complete period of record (cmp) and by month (mos). The complete record summary has columns named meanPg, sdPg, anomPg, meanRt, sdRt, anomRt. The monthly summary has DOcor, Pgcor, Rtcor for the correlations of each with the tidal cycle for the given month and anomPg and anomRt for the anomalous tallies of the metabolism estimates in each month. See the details above for a meaning of each.

References

Beck MW, Hagy III JD, Murrell MC. 2015. Improving estimates of ecosystem metabolism by reducing effects of tidal advection on dissolved oxygen time series. Limnology and Oceanography Methods. DOI: 10.1002/lom3.10062

See Also

ecometab

Examples

## Not run: 

# load library and sample data
# metab_obs and metab_dtd
 library(WtRegDO)
 data(metab_obs)
 data(metab_dtd)

 meteval(metab_obs)
 meteval(metab_dtd)

## End(Not run)

fawda123/WtRegDO documentation built on March 18, 2024, 9:04 p.m.