I'm having some trouble using the following function to load the gene id's from UCSC to annotate my RNAseq results - getting the following error when i try to implement it:
Error in .testForValidKeys(x, keys, keytype, fks) :
None of the keys entered are valid keys for 'TXNAME'. Please use the keys method to see a listing of valid arguments.
The function:
transcriptToSymbolAndDescription = function(ucsc_transcriptID)
{
df = select(TxDb.Mmusculus.UCSC.mm10.knownGene,
keys=ucsc_transcriptID,
columns='GENEID',
keytype='TXNAME')
symbolAndDescription = select(org.Mm.eg.db, df$GENEID, c("SYMBOL", "GENENAME"))
row.names(symbolAndDescription) = ucsc_transcriptID
return(symbolAndDescription)
}
any help or suggestions would be greatly appreciated! :)