Gene annotation outputs NA values for real genes

Hello everyone, I have noticed that I was removing with the na.omit() function many known genes that were not successfully annotated with annotationdbi.

res$symbol <- mapIds(org.Mm.eg.db, 
                                   keys=row.names(res),
                                   column="SYMBOL",
                                   keytype="ENSEMBL",
                                   multiVals="first")

do you have any ideas why is this happening, by any chance?

rna-seq annotation

• 1,226 views • link

written 6.1 years ago by Mozart

Can you provide some example gene IDs? - grazie mille.

• 0 views • link

written 6.1 years ago by Kevin Blighe

0 answers

No answers yet.

More posts like this

Get SYMBOL using mapIds with T2T genome ?
written 2.1 years ago by Picasa

Hi I use to retrieve symbol from ENSEMBL id using this function: resOrdered$symbol <- mapIds(org.Hs.eg.db, keys=row.names(resOrdered), column="SYMBOL", keytype="ENSEMBL", multiVals="first") but it only works with the …
NAs using org.Hs.eg.db for Ensembl ID to Gene Symbol annotation
written 3.8 years ago by Jakpa •

Hi Everyone, I did Ensembl ID to Gene Symbol annotation using `org.Hs.eg.db` with this code: annotn = mapIds(org.Hs.eg.db, keys = rownames(res), keytype = "ENSEMBL", column …
What should I do with the following error? None of the keys entered are valid keys for 'ENSEMBL'.
written 4.3 years ago by abradford •

> library("AnnotationDbi") > library("org.Hs.eg.db") > anno_df$gene_name <- mapIds(org.Hs.eg.db, + keys=anno_df$gene_id, + column="SYMBOL", + keytype="ENSEMBL", + multiVals="first") Error in .testForValidKeys(x, keys, keytype, fks) : None of …
Error while adding gene names
written 6.0 years ago by Kumar

Hi, I am getting following ERROR. Please correct me in my script. res$symbol <- mapIds(org.Hs.eg.db, keys=row.names(res), column="SYMBOL", keytype="ENSEMBL", multiVals="first") ERROR: Error in .testForValidKeys(x, keys, keytype, …
Can't fetch pathways by entrez id?
written 8.6 years ago by Tania

**Hello Everyone** I have some expression results from GeneChip Human Genome U133 Plus 2.0 Array. They are already generated, I need to read them and …
Problem with Ensembl version identifiers after running DESeq2
written 8.8 years ago by Lila M

Hi everybody, I have a problem with my Ensembl ID after running DESEq2 (I'm using hg38 genome): dds <- DESeq(ds_matrix) res <- results(dds) baseMean log2FoldChange …
AnnotationDbi using UCSC annotation
written 9.0 years ago by jadepinket •

Below is the error I get, any suggestions would be helpful res$symbol <- mapIds(TxDb.Hsapiens.UCSC.hg38.knownGene, keys=row.names(res), column="CDSNAME", keytype="TXID", multiVals="first") Error in .testForValidKeys(x, keys, keytype, fks) : …
AnnotationDbi using UCSC annotation
written 9.1 years ago by jadepinket •

I have used UCSC GTF file to annotate my reads, now I am having trouble doing the following library("AnnotationDbi") library("org.Hs.eg.db") columns(org.Hs.eg.db) res$symbol = mapIds(org.Hs.eg.db, keys=row.names(res), …
Annotate ouput file from Deseq2
written 10.9 years ago by BM

I am trying to annotate the results output file from Desq2 so it contains gene names and symbols. The RNA-seq count file I have used …
Deseq2 gene name annotation
written 11.0 years ago by BM

I am trying to annotate the Ensemble ID in the Deseq2 results file and add a column of Gene symbols and gene names. I have …