SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Error with GTF file when using htseq-count MDonlin Bioinformatics 13 01-13-2015 08:29 AM
Error when using HTSeq count shocker8786 Bioinformatics 2 02-19-2014 10:10 PM
htseq-count error message roll Bioinformatics 3 01-30-2014 02:30 AM
Strange error when using htseq-count shhuang Bioinformatics 13 11-19-2012 12:40 AM
htseq-count error sissi Bioinformatics 0 03-20-2012 11:40 PM

Reply
 
Thread Tools
Old 09-03-2014, 10:28 AM   #1
blakeoft
Member
 
Location: Connecticut

Join Date: Oct 2013
Posts: 79
Default DiffBind: error with dba.count

Hello,

I've been using DiffBind to analyze some ATAC-seq peak data. Briefly, peaks in these datasets represent sites of open chromatin. DiffBind seems to do a very good job, but I've encountered an error in my most recent comparison. Here's the R code:

Code:
> library(DiffBind)
> diffpeaks <- dba(sampleSheet = "metadata.csv")
# suppressed output
> diffpeaks <- dba.count(diffpeaks)
# suppressed output
Error in if (sum(tokeep) < length(tokeep)) { : 
  missing value where TRUE/FALSE needed
In addition: Warning message:
In dba.multicore.init(DBA$config) :
  Parallel execution unavailable: executing serially.
I've bolded the actual error. I've looked through the dba.count() code, and I can't tell where the error is coming from. My best guess is that it's from the pv.count() function, but I can't seem to find anything about it anywhere. This is really a mystery to me (and I thought I knew what was going on in R too )...

Iíve selected subsets of the metadata with varying sizes--sometimes it works, and sometimes I get the same error. This leads me to believe that it isnít a memory issue.

Does anyone have any ideas about this?

Thanks for any help,

Blake
blakeoft is offline   Reply With Quote
Old 09-09-2014, 11:03 AM   #2
rory
Member
 
Location: Cambridge, UK

Join Date: Aug 2008
Posts: 28
Default

It appears that DiffBind can't handle the "-" character (minus sign/hyphen) in the file names. Remove them and it should work!

Cheers-
Rory
rory is offline   Reply With Quote
Old 09-22-2014, 11:27 AM   #3
lime gin
Junior Member
 
Location: missoula, mt

Join Date: Jul 2012
Posts: 3
Default

Hi Blake, just curious if this has been resolved. I'm getting the same exact error with my data; while I've pinpointed the exact subset where the error is coming from, my file names all have alphanumeric characters.
lime gin is offline   Reply With Quote
Old 09-22-2014, 11:36 AM   #4
blakeoft
Member
 
Location: Connecticut

Join Date: Oct 2013
Posts: 79
Default

Unfortunately, I have not resolved the issue yet. I'm just curious, what kind of machine are you running DiffBind on? I've used it several times without error on a Unix based HPC. But ever since I've tried using it on my windows computer, I get this error. The HPC has been down for maintenance, so I haven't been able to see if running it there solves the problem. I'll post my results whenever the cluster is back up.
blakeoft is offline   Reply With Quote
Old 09-22-2014, 11:42 AM   #5
lime gin
Junior Member
 
Location: missoula, mt

Join Date: Jul 2012
Posts: 3
Default

Thanks for your response. I'm running it on my Windows computer through RStudio. We do have a Unix based HPC, but I've never tried running R there. I'll give it a shot and see how it goes.
lime gin is offline   Reply With Quote
Old 09-23-2014, 11:02 AM   #6
lime gin
Junior Member
 
Location: missoula, mt

Join Date: Jul 2012
Posts: 3
Default

Update: I ran the same set of data on our Unix based HPC and it worked!
lime gin is offline   Reply With Quote
Old 09-23-2014, 11:05 AM   #7
blakeoft
Member
 
Location: Connecticut

Join Date: Oct 2013
Posts: 79
Default

That's good to hear. I really don't have any sort of an idea of why DiffBind (specifically dba.count()) does not seem to work on Windows. I look forward to when I can use our HPC again! Thanks for the update.
blakeoft is offline   Reply With Quote
Reply

Tags
diffbind

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 11:37 AM.


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