Imports the raw channel data from Pyroscience Workbench output files. This allows "live" analyses while the trial is still running. This does not utilize the ".pyr" file, nor the text file that is created once the trial is finished. This utilizes the raw channel data found within the "ChannelData" folder that the software makes when the trial starts.
import_pyroscience_workbench(
folder,
o2_unit = "percent_a.s.",
sal = NULL,
keep_metadata = FALSE,
split_channels = FALSE,
merge_close_measurements = "min"
)A data frame (or list of data frames) is returned.
Date and time, POSIXct format. If split_channels = FALSE (default), then the timestamp is the average of all the measurements that were merged. For details, see merge_close_measurements.
Duration of measurement trial (minutes).
Oxygen measurement in desired unit as determined by o2_unit.
Temperature recorded or defined at beginning of measurement trial.
Salinity (psu). Only displayed if sal != NULL.
Warning or error messages from Pyroscience Workbench file.
Channel columns (CH_...) are repeated for each channel.
If keep_metadata = TRUE, then the following columns are appended to the returned data frame:
Phase recorded. Phase is inversely related to O2.
Intensity is an indicator of the quality of the signal.
Ambient light on the sensor. Expressed in mV.
Warning or error messages from Pyroscience Workbench file's temperature measurement.
Atmospheric pressure (mbar).
Warning or error messages from Pyroscience Workbench file's atmospheric pressure measurement.
If split_channels = TRUE, then "CH_X_" is removed from the column names and multiple data frames are returned in a named list.
a character string. The filepath to the parent folder (directory) which contains ChannelData.
a character string. The unit of O2 measurement to be output in the data frame. Options are described in conv_o2.
numeric. If o2_unit is a concentration rather than partial pressure, the salinity of the chamber in which measurements were made must be entered here.
logical. Should metadata from the file be returned as extra columns in the returned data frame? Default is FALSE.
logical. Should a list of data frames be returned with a separate data frame for each channel? Default is FALSE.
used only when split_channels = FALSE (the default). The frequency during which measurements are taken can be set uniquely for each channel in the Pyroscience Workbench software. When this happens, measurements may, at times, be nearly synchronized. When measurements are close together in time (even if not exactly at the same moment), it may be desirable to merge them together in the same row of the output dataframe and consider them to be the same timepoint. This parameter allows you to control whether that happens, and if so, how close is "close enough". Options are:
0: Do not merge close measurements no matter how close in time (even if 1 msec apart).
Merge measurements as close as the most frequently sampled channel (e.g. if channel 1 sampled every 5 seconds, channel 2 every 2 seconds, and channel 3 every 10 seconds, then any measurements within 2 seconds of each other will be merged on the same row in the output dataframe.)
Merge measurements as close as the least frequently sampled channel (e.g. if channel 1 sampled every 5 seconds, channel 2 every 2 seconds, and channel 3 every 10 seconds, then any measurements within 10 seconds of each other will be merged on the same row in the output dataframe. Warning: this will duplicate more frequent channels. Do not let your downstream statistics be altered by artificially raising the number of observations.)
A numeric value specifying how many seconds apart is "close enough" to merge measurements to the same timepoint. This may result in duplications of the same observations across multiple rows or not merging multiple observations as expected. Examine the output carefully.
Matthew A. Birk, matthewabirk@gmail.com
import_presens, import_witrox, conv_o2
if (FALSE) {
folder <- system.file('extdata/pyro_wb/', package = 'respirometry')
import_pyroscience_workbench(folder = folder, o2_unit = 'umol_per_l', sal = c(0, 35))
# I want each channel as a separate data frame.
data_list <- import_pyroscience_workbench(folder = folder, split_channels = TRUE)
data_list$CH_2 # here's the channel 2 data frame.
}
Run the code above in your browser using DataLab