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

add_embeddings error #36

Open
Jay-Leung opened this issue Jan 18, 2022 · 0 comments
Open

add_embeddings error #36

Jay-Leung opened this issue Jan 18, 2022 · 0 comments

Comments

@Jay-Leung
Copy link

Jay-Leung commented Jan 18, 2022

Hi,

I am using Seurat for my single cell analysis and so exported the loom file with Seurat UMAP embeddings using loom <- build_loom("Data/loom/RHB251019P_RG_RL_ScopeRembed.loom", dgem=exprMat, default.embedding=umap_embeddings, default.embedding.name = "umap_embeddings") to run pyscenic CLI.

However, after I added the SCENIC tSNE and UMAP using add_visualization.py and imported the loom into R for downstream visualization (http://htmlpreview.github.io/?https://github.com/aertslab/SCENIC/blob/master/Tutorials_JupyterNotebooks/SCENIC_tutorial_2-ExploringOutput.html), I realized the Seurat UMAP embeddings (umap embeddings) got replaced by SCENIC UMAP embedding values:

embeddings

I then tried to add the Seurat UMAP embeddings back manually:
add_embedding(loom = loom, embedding = seurat_embed, name = "Seurat_UMAP"), where the embedding was acquired through seurat_embed <- as.data.frame(Seurat::Embeddings(seuratobj, reduction = 'umap')) , but I got the following error:

[1] "Adding embedding Seurat_UMAP..."
Error in if (ca_extra_embeddings_dims[1] != ca_extra_embeddings_dims[2]) { :
missing value where TRUE/FALSE needed
In addition: Warning messages:
1: In FUN(X[[i]], ...) : NAs introduced by coercion
2: In FUN(X[[i]], ...) : NAs introduced by coercion

I have checked that the dims of X and Y embeddings are equal. Could you help with this problem?

Thanks!

Regards,
Jay

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

1 participant