respirometry (version 1.2.1)

make_bins: Make time bins for MO2 calculations

Description

The width of time bins seems to be an under-appreciated consideration when calculating metabolic rates if PO2 or time are interesting covariates. The wider the bins, the higher the precision of your calculated MO2 value (more observations to average over), but at a loss of resolution of an interesting covariate. The narrower the bins, the higher the resolution of the PO2 or time covariate, but at a cost of lower precision. For Pcrit trials, I have found good success using bins of 1/10th the trial duration at the highest PO2s (where good precision is important) and 1/100th the trial duration at the lowest PO2s (where good resolution is important).

Usage

make_bins(
  o2,
  duration,
  good_data = TRUE,
  min_o2_width = 1/100,
  max_o2_width = 1/10,
  n_bins = 10
)

Arguments

o2

numeric vector of O2 observations.

duration

numeric vector of the timepoints for each observation (minutes).

good_data

logical vector of whether O2 observations are "good" measurements and should be included in analysis. Default is that all observations are TRUE.

min_o2_width

Default is 1/100th of the total "good" trial duration.

max_o2_width

Default is 1/10th of the total "good" trial duration.

n_bins

Default is 10.

Value

A data.frame with two columns is returned.

o2

The O2 value below which the corresponding bin width is applied.

width

The bin width at which all data below the corresponding O2 value will be binned.

See Also

calc_MO2

Examples

# NOT RUN {
# get O2 data
file <- system.file('extdata', 'witrox_file.txt', package = 'respirometry')
o2_data <- na.omit(import_witrox(file, split_channels = TRUE)$CH_4)

# Total trial duration is 21.783 minutes

make_bins(o2 = o2_data$O2, duration = o2_data$DURATION) # creates the default 10 bins. At the
# highest O2 levels, bin widths are 21.783/10 = 2.1783 mins and at the lowest O2 levels, bin
# widths are 0.21783 mins.

bins <- make_bins(o2 = o2_data$O2, duration = o2_data$DURATION, min_o2_width = 1/20,
max_o2_width = 1/3, n_bins = 5) # creates 5 bins. At the highest O2 levels, bin widths are
# 21.783/3 = 7.261 mins and at the lowest O2 levels, bin widths are 21.783/20 = 1.089 mins.

(mo2 <- calc_MO2(duration = o2_data$DURATION, o2 = o2_data$O2,
bin_width = bins, vol = 10, temp = o2_data$TEMP, sal = o2_data$SAL))

# }