Skip to contents

Reorder idents by the gene expression

Usage

srt_reorder(
  srt,
  features = NULL,
  reorder_by = NULL,
  layer = "data",
  assay = NULL,
  log = TRUE,
  distance_metric = "euclidean",
  verbose = TRUE
)

Arguments

srt

A Seurat object.

features

Features to plot: a character vector or a named list of assay gene names or numeric metadata columns.

reorder_by

Reorder groups instead of idents.

layer

Assay layer to use.

assay

Assay to use. NULL uses the default assay.

log

Whether log1p transformation needs to be applied.

distance_metric

Metric to compute distance.

verbose

Whether to print the message. Default is TRUE.

Examples

data(pancreas_sub)
pancreas_sub <- RunStandardWorkflow(pancreas_sub)
#>  [2026-08-30 05:58:54] Start standard processing workflow...
#>  [2026-08-30 05:58:54] Checking a list of <Seurat>...
#> ! [2026-08-30 05:58:54] Data 1/1 of the `srt_list` is "unknown"
#> Warning: Data 1/1 of the `srt_list` is "unknown"
#>  [2026-08-30 05:58:54] Perform `NormalizeData()` with `normalization.method = 'LogNormalize'` on 1/1 of `srt_list`...
#>  [2026-08-30 05:58:54] Perform `FindVariableFeatures()` on 1/1 of `srt_list`...
#>  [2026-08-30 05:58:54] Use the separate HVF from `srt_list`
#>  [2026-08-30 05:58:54] Number of available HVF: 2000
#>  [2026-08-30 05:58:55] Finished check
#>  [2026-08-30 05:58:55] Perform `ScaleData()`
#>  [2026-08-30 05:58:55] Perform pca linear dimension reduction
#>  [2026-08-30 05:58:55] Use stored estimated dimensions 1:23 for Standardpca
#>  [2026-08-30 05:58:55] Perform `Seurat::FindClusters()` with `cluster_algorithm = 'louvain'` and `cluster_resolution = 0.6`
#>  [2026-08-30 05:58:55] Reorder clusters...
#>  [2026-08-30 05:58:55] Skip `log1p()` because `layer = data` is not "counts"
#>  [2026-08-30 05:58:55] Perform umap nonlinear dimension reduction
#>  [2026-08-30 05:59:05] Standard processing workflow completed
pancreas_sub <- srt_reorder(
  pancreas_sub,
  reorder_by = "SubCellType",
  layer = "data"
)
#>  [2026-08-30 05:59:05] Skip `log1p()` because `layer = data` is not "counts"