Question

diffbind and macs2 output

1

Entering edit mode

tony ▴ 10

@tony-9402

Last seen 8.9 years ago

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.

diffbind macs2 • 5.3k views

ADD COMMENT • link updated 8.9 years ago by Gord Brown ▴ 670 • written 8.9 years ago by tony ▴ 10

score 1 · Answer 1 · 2015-12-21

1

Entering edit mode

Gord Brown ▴ 670

@gord-brown-5664

Last seen 3.9 years ago

United Kingdom

Hi,

When you use the "--broad" option, the .xls output format is slightly different from the usual MACS format (missing the "summit" column). DiffBind doesn't account for that, unfortunately (though I'll put it on the to-do list).

The "broadPeak" and "gappedPeak" formats are almost "bed" format; just use the first 6 columns only, and rename them to have a ".bed" suffix. This command will extract the first 6 columns:

> cut -f 1-6 something.broadPeak > something.bed

Hope this helps...

- Gord

ADD COMMENT • link 8.9 years ago Gord Brown ▴ 670

0

Entering edit mode

(And change your sample sheet to show that the files are 'bed', not 'macs' format)

ADD REPLY • link 8.9 years ago Gord Brown ▴ 670

0

Entering edit mode

Hi,

17 months later does Diffbind handles .xls format from MACS2 --broad option ? or do we still have to use your trick ?

Just to know , thanks !

Edit : Seems to work using macs in Peaks Column

ADD REPLY • link 7.5 years ago ZheFrench ▴ 60

0

Entering edit mode

hello!

if i want to analysis broadpeak ,can i use this file ---macs.xls? and select 'macs' format in the samplesheet now?

thanks

ADD REPLY • link 5.9 years ago 1692679247 • 0