This is a test version of Biostars. For the public version, visit https://www.biostars.org.
reached iteration limit in SCTransform

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 warning not stop until I stop running the code. Does anyone face the same issue? Is there any solution?

What is alternative function to replace SCTransform?

Thanks

sctransform rna-seq seurat
• 4,413 views • link
updated by GenoMax GenoMax • written by re_raz • re_raz

1 answer

This error is related to a problem with parameter estimation (theta) on certain cells when fitting the negative binomial regression. I've run into this problem before, and I found a few issue responses from the authors of Seurat/SCTransform that these warnings can be ignored, because that parameter is not used downstream (due to parameter regularization). More info about this in the github issues here and here.

• 0 views • link
written by rpolicastro rpolicastro

Thank you for help, but the warning not stop and I have to stop the code

• 0 views • link
written by re_raz • re_raz

You can ignore the warnings. Let the code run and it should eventually finish.

• 0 views • link
written by rpolicastro rpolicastro

I will try that, Thanks

• 0 views • link
written by re_raz • re_raz

Not working for me, and I got a message Error: cannot allocate vector of size 18 GB although I am running my code on a remote server and memory.limit() give Inf.

• 0 views • link
written by re_raz • re_raz

Are you sure your remote server has no memory limit?

• 0 views • link
written by rpolicastro rpolicastro

I check by memory.limit() and give me Inf

• 0 views • link
written by re_raz • re_raz

Log in to answer this question.

More posts like this