Extract all protein-coding genes from TxDb.
Hi
I am trying to process the TxDb.Mmusculus.UCSC.mm10.knownGene annotation database in R.
I want to extract all the information of protein-coding genes from it and still keep it as a TxDb structure. I am not sure what the correct operation is?
I have tried:
txdb <- TxDb.Mmusculus.UCSC.mm10.knownGene
all_genes <- genes(txdb, columns = "CDSID")
# then filter out CDSID is NA and convert back to TxDb
It can work, but it will lose a lot of information, making Txdb no longer usable for my downstream.
• 650 views
•
link
written
by
octpus616 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
should I remove the transcripts from annotate peak region in ATAC-seq
written by QX •
-
How to use NCBI annotation instead of Ensembl in analysing Transcriptome from RNAseq data
written by Vojtěch •
-
locateVariants Fails to Give the Correct Gene ID
written by cisplatin1987 •
-
How to choose output when using a "By" function from a TxDb object; GenomicFeatures R
written by jon.klonowski
-
map subsetByOverlaps back to input coordinates
written by brtawe •
-
How to resolve NAs when annotating a diff gene list with org.Hs.eg.db terms?
written by TonyCN •
-
Merge columns adding conditions for every row
written by miriam.riquelmep •
-
WGCNA analysis - use only linkage coefficients for results analysis
written by hennet.lauriane •
-
locateVariants missing some mutation site
written by noeD
-
Problem with extracting information from txdb object!
written by Parham
