Previously I asked a question here about finding overlapping and unique peaks between two TFs from chip-seq experiments and I got really good answers for that. I have been trying to test one of the methods to achieve this task as was suggested. This method which can be used in R is DiffBind. The problem that I am facing with this is that I just cannot create the object to apply any of the functions from this package. I have peak files in bed format for the two TFs that look like this:
chr1 79092406 79094792 chr1-559 442.2 + chr1 198383685 198385392 chr1-2813 439.5 + chr1 200988993 200993982 chr1-141 400.5 + chr1 63836348 63842294 chr1-24549 296.8 +
and here's the DiffBind manual:
This is how im trying to create the object. I've converted my files to csv and then applied this:
TF1_dba = dba(sampleSheet="TF1.csv")
but it keeps giving me this error:
1 NA raw Error in file(file, "rt") : invalid 'description' argument
Can anyone tell me what I'm doing wrong?
Thanks a lot!!!