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

Receptor-ligand analysis cannot be reproduced #922

Open
Amberzhang66 opened this issue Dec 2, 2024 · 0 comments
Open

Receptor-ligand analysis cannot be reproduced #922

Amberzhang66 opened this issue Dec 2, 2024 · 0 comments
Assignees

Comments

@Amberzhang66
Copy link

Dear authors,

Description

I'm following your tutorial:
https://squidpy.readthedocs.io/en/stable/notebooks/examples/graph/compute_ligrec.html

However, I cannot reproduce it.

When running:

sq.pl.ligrec(res, source_groups="Erythroid", alpha=0.005)

I encountered errors below:

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Cell In[7], line 1
----> 1 sq.pl.ligrec(res, source_groups="Erythroid", alpha=0.005)

File [~/.conda/envs/py39_2/lib/python3.9/site-packages/squidpy/pl/_ligrec.py:296](https://ondemand.carc.usc.edu/node/a03-10.hpc.usc.edu/49178/lab/tree/~/.conda/envs/py39_2/lib/python3.9/site-packages/squidpy/pl/_ligrec.py#line=295), in ligrec(adata, cluster_key, source_groups, target_groups, means_range, pvalue_threshold, remove_empty_interactions, remove_nonsig_interactions, dendrogram, alpha, swap_axes, title, figsize, dpi, save, **kwargs)
    293 pvals = pvals[pvals <= pvalue_threshold]
    295 if remove_empty_interactions:
--> 296     pvals, means = filter_values(pvals, means, mask=~(pd.isnull(means) | pd.isnull(pvals)), kind="NaN")
    297 if remove_nonsig_interactions and alpha is not None:
    298     pvals, means = filter_values(pvals, means, mask=pvals <= alpha, kind="non-significant")

File [~/.conda/envs/py39_2/lib/python3.9/site-packages/squidpy/pl/_ligrec.py:206](https://ondemand.carc.usc.edu/node/a03-10.hpc.usc.edu/49178/lab/tree/~/.conda/envs/py39_2/lib/python3.9/site-packages/squidpy/pl/_ligrec.py#line=205), in ligrec.<locals>.filter_values(pvals, means, mask, kind)
    203 means = means.loc[mask_rows]
    205 if pvals.empty:
--> 206     raise ValueError(f"After removing rows with only {kind} interactions, none remain.")
    208 mask_cols = mask.any(axis=0)
    209 pvals = pvals.loc[:, mask_cols]

ValueError: After removing rows with only NaN interactions, none remain.

Could you please help with the problem here? Thank you so much!

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

2 participants