Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Error message in cuffcompare Line RNA Sequencing 3 09-20-2012 12:05 PM
DEXSeq error xguo Bioinformatics 4 12-06-2011 05:31 AM
gatk python script m_elena_bioinfo Bioinformatics 7 11-20-2011 06:40 AM
HTSeq Script from DEXSeq Reports Assertion Fail in SAM file FuzzyCoder Bioinformatics 5 09-27-2011 08:52 AM
GEOquery error message akolman Bioinformatics 8 06-06-2011 07:26 AM

Thread Tools
Old 09-07-2011, 12:55 PM   #1
Location: Rockville, MD

Join Date: Jun 2010
Posts: 10
Default DEXSeq python script error message

When I try to run, I got a python error at line 29:
exons[f.iv] += ( f.attr['gene_id'], f.attr['transcript_id'] )
TypeError: unsupported operand type(s) for +=: 'StepVector' and 'tuple'
Anyone has any idea how to fix it?
shilez is offline   Reply With Quote
Old 09-07-2011, 11:18 PM   #2
Jane M
Senior Member
Location: Paris

Join Date: Aug 2011
Posts: 239

I have no idea, but I managed to run this script, so there is no error in it. Nevertheless, if I remeber well, a certain version of python is needed...
Jane M is offline   Reply With Quote
Old 09-08-2011, 05:27 AM   #3
Simon Anders
Senior Member
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994

Could it be that you have an outdated version of HTSeq installed?
Simon Anders is offline   Reply With Quote
Old 09-08-2011, 10:42 AM   #4
Location: Rockville, MD

Join Date: Jun 2010
Posts: 10

Thank you, guys.
Once I update the HTSeq, it works fine.
However, I do got the error message "Same name found on two chromosomes".
Does it mean we need to remove the duplicated genes from the annotation gtf file first before running

*and also an error from HTSeq
ValueError: start is larger than end

Last edited by shilez; 09-08-2011 at 10:57 AM.
shilez is offline   Reply With Quote
Old 09-08-2011, 10:55 PM   #5
Senior Member
Location: Heidelberg

Join Date: Aug 2010
Posts: 165

The first error can be solved by adding an additional ID tag to these genes. (e.g. if the gene is called A, you could write A_dup1 and A_dup2 foreach one so that the gene names become unique). They are usually not to many in the annotation files. (I have done it manually for the ensembl mouse gtf annotation).

The second I have also seen it, for some genes, the start value is higher to the end value. I guess this is a mistake in the gtf files. Can be corrected just by shifting the columns in this cases. Also, there are usually not to many.
areyes is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:31 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO