I'm trying to load in peaksets from a BED file (the summits output from macs2 callpeaks), and I'm running into an error

Here's what the first line of my BED file looks like:

chr1    4785666 4785667 BMDM-RNAPIIpS2-CTR-1-ChIPseq-GTACACCT_S9_L003_R1_001_peak_1 18.38304

And my DiffBind call us

ctr <- dba.peakset(NULL, 
               peaks="BMDM-RNAPIIpS2-CTR-1-ChIPseq-GTACACCT_S9_L003_R1_001_summits.bed",
               peak.caller="bed",
               sampID="CTR-1",
               replicate=1)

My error reads "Error in if file.info(peaks)$size > 0) { : missing value where TRUE/FALSE needed"

What is the issue? I'm new to DiffBind so I'm not sure what's going on. Thank you!

You probably shouldn't be using the *_summits.bed files, which identify only the point of maximal pileup, but rather the *_peaks.xls files, which include the full enriched intervals.