read_data | R Documentation |
Reads input seahore data from an excel Seahorse Wave File. It assumes your data is background normalized.
read_data(
rep_list,
norm = NULL,
sheet = 2,
delimiter = " ",
norm_column = "exp_group",
norm_method = "minimum"
)
rep_list |
A list of Seahorse Wave excel export files. One file per
replicate. If your data is in a directory called "seahorse_data", use
|
norm |
A csv file with the experimental groups and their normalization
values. Leave unset if normalization is not required. See |
sheet |
The number of the excel sheet containing the long-form Seahorse data. Default is 2 because the long-form output from Seahorse Wave is on sheet 2 |
delimiter |
The delimiter between the group name and the assay type in the Group column of the wave output. e.g. "Group1 MITO" would use a space character as delimiter. |
norm_column |
Whether to normalize by |
norm_method |
How to normalize each well or experimental group (specified by
See the |
Although ceas enables integration of multiple biological and/or technical replicates, previous work has reported high inter-plate variation (Yepez et. al 2018). If you don't want your replicate data combined, you can either:
make sure that the names of the common groups between the replicates are different.
in downstream analyses (get_energetics_summary
, bioscope_plot
,
rate_plot
, atp_plot
), use sep_reps = TRUE
to do all calculations and
plotting separately for each replicate.
NOTE: to maintain backwards compatibility sep_reps
is currently
FALSE
by default, but will be set to TRUE
in a future release.
a seahorse_rates table
YƩpez et al. 2018 OCR-Stats: Robust estimation and statistical testing of mitochondrial respiration activities using Seahorse XF Analyzer PLOS ONE 2018;13:e0199938. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1371/journal.pone.0199938")}
rep_list <- system.file("extdata", package = "ceas") |>
list.files(pattern = "*.xlsx", full.names = TRUE)
seahorse_rates <- read_data(rep_list, sheet = 2)
head(seahorse_rates, n = 10)
# normalization by well using raw cell count or protein quantity
norm_csv <- system.file("extdata", package = "ceas") |>
list.files(pattern = "well_norm.csv", full.names = TRUE)
seahorse_rates.norm <- read_data(
rep_list,
norm = norm_csv,
norm_column = "well",
norm_method = "self",
sheet = 2
)
head(seahorse_rates.norm, n = 10)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.