Find shared and dataset-specific markers
Applies various filters to genes on the shared (W) and dataset-specific (V) components of the factorization, before selecting those which load most significantly on each factor (in a shared or dataset-specific way).
getFactorMarkers( object, dataset1 = NULL, dataset2 = NULL, factor.share.thresh = 10, dataset.specificity = NULL, log.fc.thresh = 1, pval.thresh = 0.05, num.genes = 30, print.genes = FALSE, verbose = TRUE )
object |
|
dataset1 |
Name of first dataset (default first dataset by order) |
dataset2 |
Name of second dataset (default second dataset by order) |
factor.share.thresh |
Use only factors with a dataset specificity less than or equalt to threshold (default 10). |
dataset.specificity |
Pre-calculated dataset specificity if available. Will calculate if not available. |
log.fc.thresh |
Lower log-fold change threshold for differential expression in markers (default 1). |
pval.thresh |
Upper p-value threshold for Wilcoxon rank test for gene expression (default 0.05). |
num.genes |
Max number of genes to report for each dataset (default 30). |
print.genes |
Print ordered markers passing logfc, umi and frac thresholds (default FALSE). |
verbose |
Print messages (TRUE by default) |
List of shared and specific factors. First three elements are dataframes of dataset1- specific, shared, and dataset2-specific markers. Last two elements are tables indicating the number of factors in which marker appears.
## Not run: # ligerex (liger object), factorization complete input markers <- getFactorMarkers(ligerex, num.genes = 10) # look at shared markers head(markers[[2]]) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.