Skip to contents

This function retrieves counts from a number of methods that are available for normalisation, with the default of limma-voomLmFit.

Usage

compute_single_de(
  data = NULL,
  treatment_samples = NULL,
  control_samples = NULL,
  method = NULL,
  batch = 1,
  k = 2,
  spikes = NULL
)

Arguments

data

A tidyseurat object merged with metadata. Must contain columns "Well_ID", "Row", "Column"

treatment_samples

Value in the column "combined_id" representing replicates of treatment samples in the data

control_samples

Value in the column "combined_id" representing replicates of control samples in the data

method

One of "Seurat_wilcox", "DESeq2", "edgeR", "RUVg", "RUVs", "RUVr", "limma_voom"

batch

Either empty, a single value, or a vector corresponding to the number of samples

k

Parameter k for RUVSeq methods, check RUVSeq tutorial

spikes

List of genes to use as spike controls

Value

Data frame of DE counts

Examples

data(mini_mac)
treatment_samples="Staurosporine_0.1"
control_samples<-"DMSO_0"
top_table <- compute_single_de(mini_mac, treatment_samples, control_samples, method = "limma_voom")