Converting raw counts or DEseq2 output to TxDb container
Interested in taking raw counts or DEseq2 output and importing into a TxDb container for downstream comparisons. Is there a straightforward way to go about this? Thank you.
• 3,032 views
•
link
updated
by
tarek.mohamed
•
written
by
rbronste
•
written
by
rbronste
1 answer
Hi, you can extract normalized raw counts via
dds<-DESeq(dds)
NormCounts<-counts(dds,normalized=TRUE)
• 0 views
•
link
written
by
tarek.mohamed
I do not get your point, can you clarify what do you mean. Thanks
• 0 views
•
link
written
by
tarek.mohamed
Log in to answer this question.
More posts like this
-
Quantification of RNAseq data
written by ahmedsalahtvi •
-
Is this featureCounts output normal, and how can I process it for DESeq2 analysis?
written by DdogBoss •
-
Can I normalize several samples with DESeq2 and then do pairwise comparisons?
written by bioinfo •
-
Inter-species RNAseq normalisation
written by pprice3 •
-
DESeq multiple comparisons post normalization
written by zapsep •
-
Consensus motif of closely related TFs
written by rbronste
-
Filter HT-Seq output for gene subclasses
written by rbronste
-
Official gene symbol list to BED file
written by rbronste
-
Calculating FPKM values from HT-counts
written by rbronste
-
RNA Seq analysis of virally infected cells
written by divya.nandakumar •
Could you give more details about what you want to do? a TxDB container contains transcript information for a given organism. You can't place count data in the container.
So I am trying to limit the transcripts I run my ChIP-seq peaks against in the ChIPseeker package to those that are actually expressed in our system, instead of running the peaks and annotating against all transcripts.
Okay I hope I've understood correctly, try the following:
You can then supply this new TxDB object to ChIPseeker.
Thats exactly what I wanted to do, thanks!
No problem, glad we solved the issue!
What would be your suggestion for how to best filter the GRanges object I get when I pull all transcripts from TxDb.Mmusculus.UCSC.mm10.knownGene against a DESeq output? Thanks again.
For each gene in your DESeq2 output, get its corresponding Entrez ID (using the mapIds function). Then map this ID to the corresponding TX NAME and just filter with this list.