SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
cufflinks "waiting for 1 threads to complete" letusgo Bioinformatics 6 01-17-2013 11:07 AM
Variable length paired-end read processing using Cufflinks ZhengXia Bioinformatics 0 10-03-2011 08:23 AM
Is it normal for cufflinks produce different results in multiple run superligang Bioinformatics 3 03-03-2011 01:40 AM
How to get the multiple matches of a same read in the results biobee07 Bioinformatics 1 02-18-2011 04:42 PM
Can maq support multiple threads? genelab Bioinformatics 1 07-21-2010 12:00 AM

Reply
 
Thread Tools
Old 09-06-2012, 07:01 AM   #1
david.eby
Junior Member
 
Location: Massachusetts

Join Date: Sep 2012
Posts: 3
Default Variable results with cufflinks 2.0.2 using multiple threads

I've been doing some testing of the cufflinks 2.0.2 tools and have been getting variable results when using multiple threads. That is, running cufflinks on the same BAM file multiple times with (for example) '-p 8' gives different results in the output files on each run. Likewise for cuffmerge (I haven't checked cuffdiff yet - that's next). I'm testing on the data from the March 2012 Nature Protocols paper using the pre-built Linux binaries.

Running single-threaded seems to be fine and always give identical results (the cuffmerge abundance results vary slightly, though that seems to be in the realm of rounding differences).

Is this a known issue? From the website it looks like there was a similar issue in v2.0.1. Are there any workarounds other than running single threaded? The multithreaded code is so much faster.
david.eby is offline   Reply With Quote
Old 09-25-2012, 09:35 AM   #2
david.eby
Junior Member
 
Location: Massachusetts

Join Date: Sep 2012
Posts: 3
Default

To reply to my own question: after doing much more testing it looks like everything is fine. The differences that I was seeing were mostly related to a combination of ordering of data in the output files and the variable way in which the 'CUFF' identifiers are assigned.
david.eby is offline   Reply With Quote
Reply

Tags
cufflinks 2.0.2, multi-threading

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 07:44 AM.


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