Hi all,
I have a gene list with refseq names as input for DeSeq and EdgeR (using both), I want to add official gene names for the output as well. I tried the following in EdgeR but it did not work:
library(org.Hs.eg.db)
idfound <- y$genes$RefSeqID %in% mappedRkeys(org.Hs.egREFSEQ)
y <- y[idfound,]
dim(y)
Any good ways to do it?
I have a gene list with refseq names as input for DeSeq and EdgeR (using both), I want to add official gene names for the output as well. I tried the following in EdgeR but it did not work:
library(org.Hs.eg.db)
idfound <- y$genes$RefSeqID %in% mappedRkeys(org.Hs.egREFSEQ)
y <- y[idfound,]
dim(y)
Any good ways to do it?
Comment