Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

getting error when running gsea #178

Open
ChenyuZhang97 opened this issue Sep 26, 2024 · 0 comments
Open

getting error when running gsea #178

ChenyuZhang97 opened this issue Sep 26, 2024 · 0 comments

Comments

@ChenyuZhang97
Copy link

Hi, I am trying to run the code pasted below, and get this error:

LR_List <- liana::select_resource('CellChatDB')[[1]] %>%
            select(ligand = source_genesymbol,
                   receptor = target_genesymbol) %>%
            unite('interaction',
                  ligand,
                  receptor,
                  sep = '^') %>%
            pull(interaction)
          
          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 = 'mouse', #using mouse doesn't support pathwaydb="KEGG", you have to change it to pathwaydb="GOBP" or pathwaydb="Reactome"
                                                     pathwaydb = 'GOBP', #GOBP, KEGG,  Reactome
                                                     readable_name = TRUE)
         
          LR_Loadings <- factors[['interactions']] %>%
            column_to_rownames("lr")
          
          GSEA_Results <- c2c$external$run_gsea(loadings = LR_Loadings,
                                                lr_set = LR_Set,
                                                output_folder = 'D:\\LIANA\\Comparison\\Results',
                                                weight = as.integer(1),
                                                min_size = as.integer(5),
                                                permutations = as.integer(999),
                                                processes = as.integer(6),
                                                random_state = as.integer(6),
                                                significance_threshold = as.numeric(0.05))



  0%|          | 0/11 [00:00<?, ?it/s]2024-09-26 23:56:28,134 [WARNING] Duplicated values found in preranked stats: 0.17% of genes
The order of those genes will be arbitrary, which may produce unexpected results.
  0%|          | 0/11 [00:00<?, ?it/s]
Error in py_call_impl(callable, call_args$unnamed, call_args$named) : 
  Exception: Error parsing gmt parameter for gene sets

And by running this command reticulate::py_last_error(), I got:

Traceback (most recent call last):
  File "D:\Softwares\Anaconda\Lib\site-packages\cell2cell\external\gseapy.py", line 223, in run_gsea
    gseapy.prerank(rnk=test,
  File "D:\Softwares\Anaconda\Lib\site-packages\gseapy\__init__.py", line 357, in prerank
    pre.run()
  File "D:\Softwares\Anaconda\Lib\site-packages\gseapy\gsea.py", line 426, in run
    gmt = self.load_gmt(gene_list=dat2.index.values, gmt=self.gene_sets)
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\Softwares\Anaconda\Lib\site-packages\gseapy\base.py", line 206, in load_gmt
    genesets_dict = self.load_gmt_only(gmt)
                    ^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\Softwares\Anaconda\Lib\site-packages\gseapy\base.py", line 195, in load_gmt_only
    raise Exception("Error parsing gmt parameter for gene sets")
Exception: Error parsing gmt parameter for gene sets

── R Traceback ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
    ▆
 1. └─c2c$external$run_gsea(...)
 2.   └─reticulate:::py_call_impl(callable, call_args$unnamed, call_args$named)
 

I'm able to run the Liana x tensor-cell2cell R tutorials up to this point. Any suggestion for this? Thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant