select markers using mclust EM method

markers_mclust(
  top_markers,
  column = ".dot",
  prob = 0.99,
  s_thres = NULL,
  method = c("max.one", "remove.min"),
  plot = FALSE,
  ...
)

Arguments

top_markers: output of top_markers()
column: character, specify which column used as group label
prob: numeric, probability cutoff for 1st component classification
s_thres: NULL or numeric, only features with score > threshold will be returned, if NULL will use 2 * average probability as threshold
method: can be "max.one" or "remove.min", if to only keep features in 1st component or return features not in the last component
plot: logical, if to plot mixture density and hist
...: other params for mclust::densityMclust()

Value

a list of markers for each group

Examples

data <- matrix(rnorm(100), 10, dimnames = list(1:10))
top_n <- top_markers(data, label = rep(c("A", "B"), 5))
markers_mclust(top_n)
#> $A
#> [1] "X9"
#> 
#> $B
#> [1] "X1"
#>