get_signatures-ConsensusPartition-method {cola} | R Documentation |
Get signature rows
## S4 method for signature 'ConsensusPartition' get_signatures(object, k, silhouette_cutoff = 0.5, fdr_cutoff = ifelse(identical(diff_method, "samr"), 0.05, 0.1), scale_rows = object@scale_rows, diff_method = c("Ftest", "ttest", "samr", "pamr", "one_vs_others"), anno = get_anno(object), anno_col = get_anno_col(object), internal = FALSE, show_row_dend = FALSE, show_column_names = FALSE, use_raster = TRUE, plot = TRUE, verbose = TRUE, seed = 888, ...)
object |
A |
k |
number of partitions. |
silhouette_cutoff |
cutoff for silhouette scores. Samples with values less than it are not used for finding signature rows. For selecting a proper silhouette cutoff, please refer to https://www.stat.berkeley.edu/~s133/Cluster2a.html#tth_tAb1. |
fdr_cutoff |
cutoff for FDR of the difference test between subgroups. |
scale_rows |
whether apply row scaling when making the heatmap. |
diff_method |
methods to get rows which are significantly different between subgroups, see 'Details' section. |
anno |
a data frame of annotations for the original matrix columns. By default it uses the annotations specified in |
anno_col |
a list of colors (color is defined as a named vector) for the annotations. If |
internal |
used internally. |
show_row_dend |
whether show row dendrogram. |
show_column_names |
whether show column names in the heatmap. |
use_raster |
internally used. |
plot |
whether to make the plot. |
verbose |
whether to print messages. |
seed |
random seed. |
... |
other arguments. |
Basically the function applies statistical test for the difference in subgroups for every row. There are following methods which test significance of the difference:
First it looks for the subgroup with highest mean value, compare to each of the other subgroups with t-test and take the maximum p-value. Second it looks for the subgroup with lowest mean value, compare to each of the other subgroups again with t-test and take the maximum p-values. Later for these two list of p-values take the minimal p-value as the final p-value.
use SAM (from samr package)/PAM (from pamr package) method to find significantly different rows between subgroups.
use F-test to find significantly different rows between subgroups.
diff_method
can also be a self-defined function. The function needs two arguments which are the matrix for the analysis
and the predicted classes. The function should returns a vector of FDR from the difference test.
A data frame with more than two columns:
which_row
:row index corresponding to the original matrix.
fdr
:the FDR.
the mean expression (depending rows are scaled or not) in each subgroup.
Zuguang Gu <z.gu@dkfz.de>
# There is no example NULL