Open lo2811 opened 3 years ago
Are you using the same version of Seurat as before? That could also be a reason why the same code is not working anymore.
In any case, without a reproducible example I cannot say what the problem is exactly. If you are sure changes in the sctransform package are the cause, you can also try to install an older version, i.e. a specific commit. For example remotes::install_github("ChristophH/sctransform@73e2e3e")
would install the develop version as of 30. March 2021
Other commits to the develop branch are listed here
Another option is to change the parameter estimation method that sctransform is using. E.g.
if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager")
BiocManager::install("glmGamPoi")
epi <- SCTransform(epi, batch_var="dataset", method = "glmGamPoi")
Thank you Christoph. I tried glmGamPoi before and it worked, but will try installing the specific sctransform and Seurat version as well.
Hi Christoph,
I am running SCTransform with batch_var and encounter this error: "Error in while ((it <- it + 1) < limit && abs(del) > eps) { : missing value where TRUE/FALSE needed". Here is the full error:
I was able to use this same code to run on the same object in April/May this year. Now I need to regenerate the Seurat object and run into this error. I'm using the development sctransform, just reinstalled it on my machine this week using: remotes::install_github("ChristophH/sctransform@develop").
The error goes away if I increase the n_cells to 10,000, but I would prefer to keep everything consistent with the previous times (n_cells=5000). It is strange that I have this error now, since I'm running the exact same code on the same samples.
Thank you! Linh