FindMarkersCondition: Differential expression testing between conditions via...

In mgildea87/CVRCFunc: Functions for CVRC bioinformatics

Differential expression testing between conditions via pseudobulking and DESeq2

Description

Differential expression testing between conditions via pseudobulking and DESeq2

Usage

FindMarkersCondition(
  seurat,
  clus_ident,
  sample_ident,
  condition_ident,
  conditions,
  expfilt_counts = 1,
  expfilt_freq = 0.5,
  out_dir = "FindMarkersCondition_outs",
  alpha = 0.1
)

Arguments

seurat	A Seurat object
clus_ident	Identity for clusters. Normally 'seurat_clusters' but can be any identity
sample_ident	Sample identities. Identity class that indicates how to partition samples
condition_ident	Identity class for conditions to be tested
conditions	A vector of the 2 conditions within condition_ident to be included in the DESeq2 model.
expfilt_counts	genes with less than expfilt_counts in expfilt_freq * sample number will be removed from DESeq2 model. 1 by default.
expfilt_freq	genes that have greater than expfilt_counts in greater than expfilt_freq fraction of cells will be kept for the DESeq2 model. 0.5 by default
out_dir	Name of output directory
alpha	FDR adjusted p-value threshold for significance in plotting. 0.1 by default.

Value

.csv files with marker genes per clus_ident. .pdf files with diagnostic plots

mgildea87/CVRCFunc documentation built on Nov. 9, 2024, 7:39 p.m.