Skip to content

Running reactomeGSA on multiple resolution cluster separately #41

@Ankita-1211

Description

@Ankita-1211

Hi

Thank you for helping me running reactomeGSA for pathway analysis.

I need lil help in understanding this.

The seurat-result.rds file I am using is having multiple resolution cluster information.

my seurat@metadata look like this

<style> </style>
  orig.ident nCount_RNA nFeature_RNA percent.mt nCount_SCT nFeature_SCT SCT_snn_res.0.1 seurat_clusters SCT_snn_res.0.2 SCT_snn_res.0.3 SCT_snn_res.0.4 SCT_snn_res.0.5 SCT_snn_res.0.6 SCT_snn_res.0.7 SCT_snn_res.0.8 SCT_snn_res.0.9 SCT_snn_res.1.0 SCT_snn_res.1.1 SCT_snn_res.1.2 SCT_snn_res.1.3 SCT_snn_res.1.4 SCT_snn_res.1.5
AACACACTCGACGAGA SeuratProject 4462 2071 0.054459883 6629 2093 0 2 1 1 1 1 1 0 0 0 0 0 0 0 2 2
AACAGGGTCAGAATAG SeuratProject 19959 4112 0.091838268 8600 2649 0 3 1 1 1 1 1 0 0 0 0 0 0 0 4 4

how can I used this seurat cluster file for running reactomeGSA on each single resolution one by one.

So far when I try to set the indents as Idents(object = data) <- "SCT_snn_res.0.8"

and then running gsva_result <- analyse_sc_clusters(seuset, verbose = TRUE)
it gives me error Error: Only one identification found: 'integrated_snn_res.0.8'. Please ensure that cell / cluster ids are stored as the primary identification (Ident) of your Seurat object. Clustering has to be performed prior to this pathway analysis.

Could you pls guide me in running the reactomeGSA for each resolution cluster separately?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions