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

unlimited warning reached iteration limit in sctransform #3285

Closed
MWRR opened this issue Jul 19, 2020 · 2 comments
Closed

unlimited warning reached iteration limit in sctransform #3285

MWRR opened this issue Jul 19, 2020 · 2 comments

Comments

@MWRR
Copy link

MWRR commented Jul 19, 2020

Hi
I got hundreds of warning "reached iteration limit" when I run

Sproject <- SCTransform(Sproject, vars.to.regress = c("percent.mt", "percent.rps", "percent.rpl", "percent.rrna", "nCount_RNA", "nFeature_RNA"), verbose = FALSE, return.only.var.genes = FALSE)

and the run was killed, and the same happened when I changed return.only.var.genes to TRUE.
I checked all issues #1426, #1378, #25, but no solution work.

I stuck at this and I could not move to the next step in my code, any help, please.

@scideas
Copy link

scideas commented Jul 21, 2020

Although you won't be able to use the SCTransform algorithm, a bypass would be to run normal normlization and scaling prior to RunPCA()
The sequence would be:
NormalizeData()
ScaleData(Sproject, vars.to.regress= c("percent.mt", "percent.rps", "percent.rpl", "percent.rrna", "nCount_RNA", "nFeature_RNA") )
RunPCA()
...
This at least worked for me when I was having issues yesterday and would at least allow you to move forward to generating UMAPs in the meantime.

@timoast
Copy link
Collaborator

timoast commented Jul 24, 2020

This was addressed by @ChristophH in satijalab/sctransform#25. My guess would be that if you decreased the number of latent variables in the model you might see less of these cases where you can't reliably estimate theta, but Christoph or @yuhanH might be able to comment further

@timoast timoast closed this as completed Jul 24, 2020
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

3 participants