I would like to know how to format the macs2 peaks to be used in diffbind package. I have the following files:
omodel_model.r
omodel_peaks.broadPeak
omodel_peaks.gappedPeak
omodel_peaks.xls
I can not use the xls as I am getting an error:
Error in Summary.factor(c(926L, 1085L, 1022L, 998L, 902L, 1360L, 677L, :
‘max’ not meaningful for factors
The contents of xls file are:
chr start end length pileup -log10(pvalue) fold_enrichment -log10(qvalue) name
Thanks in advance.

