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
GSVA + SingleCellExperiment: Still experimental? #88
Comments
hi Lluís!! Thanks for reporting this. In principle, there is no g <- gsva(sc3_clusters_h, geneIds(gsc)) but I have just pushed a fix that allows the software to deal with the library(Seurat)
library(SeuratData)
InstallData("pbmc3k")
data(pbmc3k)
sce <- as.SingleCellExperiment(pbmc3k)
gs1 <- GeneSet(SymbolIdentifier(), geneIds=rownames(sce)[1:10], setName="GS1")
gs2 <- GeneSet(SymbolIdentifier(), geneIds=rownames(sce)[11:20], setName="GS2")
gsc <- GeneSetCollection(gs1, gs2)
g <- gsva(sce, gsc)
Estimating GSVA scores for 2 gene sets.
Estimating ECDFs with Gaussian kernels
|======================================================================| 100%
Warning message:
In .local(expr, gset.idx.list, ...) :
Using 'SingleCellExperiment' objects as input is still in an experimental stage.
assay(g)[, 1:2]
AAACATACAACCAC AAACATTGAGCTAC
GS1 0.4827375 0.36579394
GS2 -0.1440312 -0.09209655 Depending on the features of your data, you may encounter errors running the default GSVA method. In such a case, please try calling GSVA with The CZI project you mention is alive and kicking, we're just slower than I initially planned due to some delays in the recruitment of the taskforce for the project. On April 1st @axelklenk started working on the robustness and single-cell part (btw, Axel, please sync your devel version of GSVA with this bugfix), and on July 1st the second recruit will start working on the spatial transcriptomics part. We expect to fully support single-cell data during the coming one or two devel-cycles, and spatial transcriptomics in the following two or three. One of the first things that Axel has done has been to clean up some CRAN-specific warnings and notes for the latest April release of GSVA, so I recommend you to switch to the current release 3.17 (GSVA 1.48.1). In any case, since you're based in BCN as well, we can discuss in person any problems or suggestions you may have for GSVA, just DM me on Slack or drop me an email. |
I've seen that there is now a method for SingleCellExperiment but one get's a warning that it is experimental (but released at least since 3.16, which I am using.).
I am finding an error using it with GSVA 1.46.0:
I hope this is enough to try to troubleshoot it (sorry for not producing a reproducible example, I can make one if needed).
I think there was a project to make GSVA work with single-cell and spatial transcriptomic data. Not sure if it has yet to start or this might be relevant for those that will work in it.
Thanks for maintaining this package!
The text was updated successfully, but these errors were encountered: