Description Usage Arguments Value
Test for differential gene expression along tree using binomial test
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | aucprTestByFactor(
object,
cells.1,
cells.2,
label,
groups,
log.effect.size = 0.25,
auc.factor = 1,
min.auc.thresh = 0.1,
max.auc.thresh = Inf,
frac.must.express = 0.1,
frac.min.diff = 0,
genes.use = NULL,
min.groups.to.mark = 1,
report.debug = F
)
|
object |
An URD object |
cells.1 |
(Character vector) Cells in population of interest |
cells.2 |
(List) Cells to compare against. If a list, each entry is a population to compare against, and markers must beat all populations. |
label |
(Character) Label to use to split cells. Must be a column of |
groups |
(List) List of label values to include in each group. |
log.effect.size |
(Numeric) Minimum fold-change difference in mean expression to be considered differential (default, |
auc.factor |
(Numeric) The precision-recall AUC is determined for a random classifier is determined based on the size of populations. To be considered differential, genes must have an AUC this factor multiplied by the expected AUC of a random classifier. |
frac.must.express |
(Numeric) Fraction of cells of interest that must express the gene to consider it differential. |
frac.min.diff |
(Numeric) Minimum difference in fraction of cells expressing a gene to consider it differential. |
genes.use |
(Character vector) Which genes to test for differential expression (default |
min.groups.to.mark |
(Numeric) For how many of the groups must a gene be a marker in order to be considered differential? |
report.debug |
(Logical) If |
min.auc.threshold |
(Numeric) This acts as a lower bound for the AUC threshold, no matter how unbalanced the populations are. |
max.auc.threshold |
(Numeric) This acts as an upper bound for how high the AUC must be for a gene to be considered differential. |
a data.frame of gene expression results if report.debug=F
or a list with entries diff.exp
and stats
if report.debug=T
.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.