diff --git a/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R b/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R index af4fd7d..546d8ed 100644 --- a/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R +++ b/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R @@ -10,7 +10,8 @@ library(data.table) # min.pct = minimum percentage of cells expressing the gene # thresh.use = minimum AVERAGE DIFFERENCE for filtering, in log scale (base e), e.g. for 1.5 fold change, type log(1.5) # data.info.col = column number for columns in data.info that contain metadata (not clusters...) - + + FindMarkers.MAST <- function(object, id1, id2, min.pct, thresh.use,data.info.col = 1:13){ cells.1 <- names(object@ident[object@ident == id1])