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

Issue caused by No significant gene #14

Open
Rong-ao opened this issue Aug 5, 2024 · 2 comments
Open

Issue caused by No significant gene #14

Rong-ao opened this issue Aug 5, 2024 · 2 comments

Comments

@Rong-ao
Copy link

Rong-ao commented Aug 5, 2024

I got an error Error in expr[i, ] : subscript out of bounds when I was trying to apply function rogue() to my cell clusters. I tried to traced what happened in this function and found an issue in function ent.toli(). This function was used in rogue() code:
tmp.res <- ent.toli(tmp.res, tmp.matr, span = span, r = r, n = remove.outlier.n)

In ent.toli(), the issue was caused by code below:
sig.gene <- ent %>% dplyr::filter(p.adj < 0.05) %>% dplyr::pull(Gene)
ng <- length(sig.gene)
expr <- expr[sig.gene,]
When I check tmp.res of my cell cluster, I found there is no gene with p.adj < 0.05, causing sig.gene is an empty list and make expr[sig.gene,] an error.

So I have some questions here:

  1. How should I understand no gene with p.adj < 0.05 in a cluster? (I think they are not informative genes...right? I am a little bit confused about 'informative gene' actually.)
  2. The function ent.toli() is to remove outlier cells before calculating ROGUE as you mentioned in code. Is that necessary?
  3. Should I skip ent.toli() step manually when I have this issue?

Thanks for your help! :)

@Felix-Xie
Copy link

getting the same error and confused as to what no gene with p.adj < 0.05 indicates

@Rong-ao
Copy link
Author

Rong-ao commented Sep 19, 2024

@Felix-Xie I tried to manually skip those 3 lines of codes in ent.toli() to avoid this bug, which means that using all genes rather than filtered genes to calculate ROGUE index. It also works.
Actually this issue only occured in one cluster with high heterogeneity in my scRNA data, so I guess p.adj < 0.05 indicate the gene expression pattern in one cluster is too noisy than expectation.

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