diff --git a/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R b/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R index 681559a..af4fd7d 100644 --- a/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R +++ b/Dim_reduction_clustering_DEanalysis/DE_analysis_functions.R @@ -9,7 +9,7 @@ library(data.table) # id2 = identity 2 # 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...), default is 1:13, for our lizard and turtle files +# 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){