diffbind and macs2 output

Question: diffbind and macs2 output

4.0 years ago by

tony • 10

tony • 10 wrote:

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 • 2.8k views

ADD COMMENT • link •

modified 4.0 years ago by Gord Brown • 590 • written 4.0 years ago by tony • 10

Answer: diffbind and macs2 output

4.0 years ago by

Gord Brown • 590

United Kingdom

Gord Brown • 590 wrote:

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 written 4.0 years ago by Gord Brown • 590

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

ADD REPLY • link written 4.0 years ago by Gord Brown • 590

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 modified 2.5 years ago • written 2.5 years ago by ZheFrench • 10

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 written 12 months ago by 1692679247 • 0

Similar posts • Search »