Following the Rsem section of the tximport vigenettes ended by having a dataset of "transcripts" and "counts". still one critical step.. converting these transcripts to genes .. I'm more into gene-level.. but it's not the case as the same gene has multiple transcripts.
Is it possible to convert transcript-level to gene-level dataset to use it in DESeq2 ?
Hello,
I have a dataset contains the following columns: transcript_id, gene_id, length, effective_length, expected_count, TPM, FPKM, and IsoPct.
And for the same gene_id more than two transript_id .. so it's a transcript-level dataset
I thought I need to use tximport to convert it to counts so I could use it for DESeq2, but I didn't really get how can I convert it exactly.
1 answer
This looks like RSEM output.
Follow the instructions here: https://bioconductor.org/packages/devel/bioc/vignettes/tximport/inst/doc/tximport.html#RSEM for the isoforms.results file.
Log in to answer this question.
More posts like this
-
The difference of RNA-seq data between TCGA and ENCODE
written by dao •
-
Extraction of information from stringtie output for selected gene ids
written by kuntalasb •
-
RSEM isoform.results outputs for edgeR
written by weichengz •
-
Stringtie output gtf file only contains STRG,no original annotations
written by yiren •
-
Deseq2 with using data from ENCODE project RNA-seq data
written by jrnick7 •
-
Warning: gene "NR_111945" (on chr1) has reference transcripts on both strands?
written by Mithil Gaikwad •
-
splitncigarreads + stringtie
written by sorrymouse
-
Sringtie gff output has "STRG." instead of genid
written by mkh •
-
RSEM output interpretation
written by komal.rathi
-
How to use RSEM to get the expression level for each gene?
written by moxu
