Diffbind bioconductor, peaks overlapping
0
0
Entering edit mode
7.4 years ago

Hello

As I understood I should create a csv file with names of bed files with input, control and peaks bed files of peaksets.

Are all of these files required for analysis? Can I do only with files with peaks? I'm trying to do without initial files (just don't have them) writing a csv file as follows:

SampleID,Peaks
1,2A.peaks.bed
2,3A.peaks.bed
3,5A.peaks.bed


But using this in R I have a mistake:

> sets = dba(sampleSheet="sample.csv")
1     NA raw
Error in Summary.factor(c(1L, 1112L, 2223L, 3334L, 4445L, 5556L, 6592L,  :
max not meaningful for factors
Warning:
In read.table(samplesheet, sep = ",", stringsAsFactors = F, header = T) :

ChIP-Seq diffbind bioconductor • 4.8k views
1
Entering edit mode

Hi Marina,

I was facing similar problem. I think, unknowingly you got some spaces while writing the CSV file. Open the file in a texteditor (e.g Textedit for MAC or Notepad for Windows) and then check if there are any unwanted spaces.

The file should somewhat look like this

SampleID,XXX,YYY,AAA


Not like this

SampleID,XXX, YYY,AAA


Even check the end of a line to see any extra spaces, which is visually unrecognizable

Thanks

0
Entering edit mode

Hi, thank you for the reply, I managed to make it work )

0
Entering edit mode

Hi Marina-orlova, May you please tell how you made it work? I met the same problem now. Thanks.

0
Entering edit mode

Hi zwang,

You need to make sure that your peak file is formatted correctly. It must be in BED format and follow this form:

CHROMOSOME   START   END   SCORE