This calculates the log fold change for a gimap dataset based on the annotation and metadata provided.
calc_crispr(.data = NULL, gimap_dataset, normalized = TRUE)
Data can be piped in with tidyverse pipes from function to function. But the data must still be a gimap_dataset
A special dataset structure that is setup using the `setup_data()` function.
Default is TRUE meaning that we should expect to look for normalized data in the gimap_dataset.
if (FALSE) {
gimap_dataset <- get_example_data("gimap")
# Highly recommended but not required
run_qc(gimap_dataset)
gimap_dataset <- gimap_dataset %>%
gimap_filter() %>%
gimap_annotate() %>%
gimap_normalize(
timepoints = "day",
replicates = "rep"
) %>%
calc_crispr()
# To see results
gimap_dataset$crispr_score
}