pathview error
1
0
Entering edit mode
Fara • 0
@8e7ca809
Last seen 5 months ago
Japan

I am currently trying to conduct KEGG enrichment analysis using the Pathview package in R. However, R always responds with the following warning:

Warning: None of the genes or compounds mapped to the pathway! Argument gene.idtype or cpd.idtype may be wrong. 'select()' returned 1:1 mapping between keys and columns

Can you help me solve this problem?

This is the example of table. X1 is the gene name and X2 is the Log fold change. There are much more genes in the table.

#unlisting gene id and log fold change
Control_entrez_unlist_gene <-unlist(Control_entrez$NCBI.gene..formerly.Entrezgene..ID)
Control_entrez_unlist_fold = unlist(Control_entrez$Log.fold.change)

pathview_Control<- data.frame(matrix(ncol = 2, nrow = length(Control_entrez_unlist_gene)))

pathview_Control[1] = as.character(Control_entrez_unlist_gene)
pathview_Control[2] = Control_entrez_unlist_fold

pathview(gene.data = pathview_upregulated_PvsJControl, species = "dre", pathway.id = "dre04672", limit = list(gene=2, cpd=1))
KEGGgraph KEGG • 551 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 22 hours ago
United States

The obvious would be to check that you have genes that map to that pathway. It's unusual to have so many NA values for your NCBI gene IDs, so you might check that as well.

0
Entering edit mode

Thank you for the reply.

I already conducted the KEGG enrichment analysis and I got the pathway ID.

So there is no doubt it cannot be mapped.

Also, the reason why I have so many NA values is that I converted the Ensembl ID to Entrez ID.

Will removing the NA values help solve the problem?

ADD REPLY
0
Entering edit mode

Will removing the NA values help solve the problem?

Did you try?

And, FWIW, in the cases I had ensembl-based RNA-seq datasets I removed all genes that had no corresponding entrezid before doing enrichment analysis. For that I used the ensembl-based annoation libraries (EnsDb), and not the OrgDb.

ADD REPLY

Login before adding your answer.

Traffic: 667 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