DiffBind Issue with coordinates
0
0
Entering edit mode
2.6 years ago
aboyd003 • 0

I am analyzing ChIP data. I called peaks with MACS2 run through a RECAP wrapper script. I am attempting differential binding analysis using DiffBind. Code shown below. Please note I have tried the suggestion in biocondutor support from Stark (the developer) to explicitly set the datatype as part of dba.report.

MADSmeta <- dba(sampleSheet="SampleSheet.csv")


The resulting output has coordinates outside of the contig sizes. Peak in the DiffBind results shows up is

7000000037415150 5009119 5009427

But the contig size (listed below) is substantially less than the start of the peak.

7000000037415150 46831

ChIP-Seq • 833 views
0
Entering edit mode

This problem also occurs when I use diffbind under windows, but when I use diffbind under Linux, it will calculate quickly.