diffBind - bam not accessible
0
0
Entering edit mode
@eli7javasky-9060
Last seen 9.0 years ago
European Union

Hi,

I'm trying to analyze chip-seq histone modification data with diffBind. When I try running diffBind, it tells me that my bam files cannot be accessed. It was able to make the plot based on the peaks so it is able to access my peak files. This is the code I've been trying:

mod=dba(sampleSheet=paste0(args[1]))

plot(mod)  #this works

mod=dba.count(mod, minOverlap=2)

plot(mod)

 

My sampleSheet contains the following columns and all of the location paths are correct (it's a .csv file):

 

SampleID, Condition, Replicate, bamReads, ControlID, bamControl, Peaks, ScoreCol
 

Do you know why it cannot access my bam files? The bam.bai files are in the same location as the bams. Any help would be greatly appreciated. Thank you so much,

Elisheva J

diffbind chip-seq bam • 2.6k views
ADD COMMENT
0
Entering edit mode

Hi,

Do the paths in your sample sheet contain tildes, by any chance?  (e.g. ~/bam/my_file.bam)?  The release version of DiffBind doesn't process them correctly (it's fixed in devel).  If so, try replacing them with absolute paths.  Also, leading or trailing white space in the filenames will cause problems (also fixed in devel).  If it's neither of those, post or email the sample sheet and I'll take a look at it. 

Cheers,

 - Gord

ADD REPLY

Login before adding your answer.

Traffic: 673 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6