I am using Genesum package for estimating gene abundnace from transcript abundance. As input, I give sailfish generated expression file "Quant.sf" and gene annotation "GTF" file. But i get the following error
**Parsing input expression file
terminate called after throwing an instance of 'std::invalid_argument'
what(): stoi
Aborted (core dumped)**
Can someone help me solve this issue ??
TIA
1 answer
Hi Pam,
The reason for this is that Sailfish (& Salmon) have since changed their default output format (actually, making them simpler to read with standard tsv parsing functions etc.), and I have not yet updated GeneSum to keep pace. However, I should say that I'd actually recommend tximport, which solves the same problem while offering some more options than GeneSum. Also, it's worth noting that Sailfish and Salmon now actually have built-in support for aggregating their expression estimates to the gene level (though I'd still probably recommend using tximport).
Hi Rob, I have used tximport and now I have counts from each sample for DEG analysis. I need TPM values(summarised for gene) for downstream analysis like heatmaps etc. I enabled "scaledTPM" in tximport but i get in countsfromabundance column just "scaledTPM" string and not any values ??!!
Hi Pam,
Tx-import only generates counts as its output. The "countsfromabundance" field just described the method used to generate the counts from the input abundances. That said, given counts, computing TPM should be very simple. The "length" of a gene should be represented as the abundance-weighted combination of the lengths of its isoforms, and the count provided by Tx-import gives the read count for the gene. Perhaps Mike Love (the main tx-import author) might even have a function to perform this transformation and get back gene-level TPM. I'll point him here on twitter!
Check ?tximport and scan to the section Value. This section describes the object that is returned by tximport().
The return object is a list which contains three matrices, one of which is "abundance", this is the TPM summarized to gene-level.
Log in to answer this question.
-
i got an error message in bowtie2 run
written by carlosgonzalezcruz327 •
-
File write failed (wrong size) when filtering a bam file with interset
written by ManuelDB
-
Error while using bowtie2
written by Jimpix •
-
freebayes running error on public data
written by anithanagaraj93 •
-
Freebayes error message
written by ocortes •
-
Error in Cuffdiff
written by harshraje19 •
-
terminate called after throwing an instance of 'std::bad_alloc?
written by star
-
Problem with bedtools subtract
written by zion22
-
Kallisto run error
written by shoujun.gu
-
Control-FREEC error message
written by 2nelly
