saezlab / liana

LIANA: a LIgand-receptor ANalysis frAmework
https://saezlab.github.io/liana/
GNU General Public License v3.0
181 stars 31 forks source link

Running GSEA for for mouse Liana x tensor-cell2cell #175

Open wgao688 opened 2 months ago

wgao688 commented 2 months ago

I am trying to run this code for a mouse dataset (https://ccc-protocols.readthedocs.io/en/latest/notebooks/ccc_R/QuickStart.html):

# Generate list with ligand-receptors pairs in DB
lr_pairs <- liana::select_resource('Consensus')[[1]] %>%
    select(ligand = source_genesymbol, receptor = target_genesymbol)
lr_list <- lr_pairs %>%
     unite('interaction', ligand, receptor, sep = '^') %>%
     pull(interaction)

# Specify the organism and pathway database to use for building the LR set
organism = "human"
pathwaydb = "KEGG"

# Generate ligand-receptor gene sets
lr_set <- c2c$external$generate_lr_geneset(lr_list = lr_list,
                                           complex_sep='_', # Separation symbol of the genes in the protein complex
                                           lr_sep='^', # Separation symbol between a ligand and a receptor complex
                                           organism=organism,
                                           pathwaydb=pathwaydb,
                                           readable_name=TRUE
                                           )

Change the organism = "mouse" does not work. I have been able to run the Liana x tensor-cell2cell up to this point. Thanks!

earmingol commented 1 month ago

Hi! using mouse doesn't support pathwaydb="KEGG", you have to change it to pathwaydb="GOBP" or pathwaydb="Reactome"