'find_gene_modules' error in the Monocle 3 tutorial
2
0
Entering edit mode
14 months ago
vHelmholtz ▴ 20

Dear all,

I have just started the 'Monocle 3' and followed it using the same code and data in the tutorial here.

When I run the code as below,

pr_graph_test_res <- graph_test(cds_subset, neighbor_graph="knn", cores=8)
pr_deg_ids <- row.names(subset(pr_graph_test_res, morans_I > 0.01 & q_value < 0.05))

gene_module_df <- find_gene_modules(cds_subset[pr_deg_ids,], resolution=1e-3)

it returns the error as below.

*Error in uwot(X = X, n_neighbors = n_neighbors, n_components = n_components, :

n_neighbors must be smaller than the dataset size*

I will very much appreciate it if you give me any advice or suggestion about this error.

Thanks,

Joshua

Monocle 3 Monocle RNA-Seq Single cell • 926 views
ADD COMMENT
0
Entering edit mode
14 months ago
vHelmholtz ▴ 20

The code works. There is no error. By mistake, I skipped running the previous code. Sorry for the confusion.

ADD COMMENT
0
Entering edit mode
12 months ago
krovi138 • 0

Hello,

I'm encountering the same exact error after running all the previous lines of code. How did you go about rectifying this? I'd love any help you could provide! Thanks!

ADD COMMENT
0
Entering edit mode

Hi,

Most likely your rownames(cds@preprocess_aux$gene_loadings) are different gene ID type than rownames(cds). It could happen e.g. when you change rownames(cds) after running preprocess_cds(). Just modify your rownames(cds) and make sure that your gene subset you provide for find_gene_modules() is present in rownames(cds@preprocess_aux$gene_loadings).

Best, Michal

ADD REPLY

Login before adding your answer.

Traffic: 2398 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6