Exception AssertionError: 'Lambd

Exception AssertionError: 'Lambda must > 0, however we got 0' in 'MACS2.Prob.poisson_cdf' ignored

1. Manually add pseudocount to your bedGraph file

$ macs2 bdgopt -i treat.bdg -m add -p 0.1 -o treat_0.1.bdg
$ macs2 bdgopt -i control.bdg -m add -p 0.1 -o control_0.1.bdg

2. Run bdgcmp

$ macs2 bdgcmp -t treat_0.1.bdg -c control_0.1.bdg -m ppois -o ppois.bdg

3. Call peaks on score track using a cutoff

$ macs2 bdgpeakcall -i ppois.bdg -c 1.301 -l 245 -g 100 -o peaks.bed

如果实验组和对照组测序深度不一样，在进行比较（macs2 bdgcmp）之前需要将二者测序深度进行归一化；如果bdg来自 macs2 callpeak 则不需要， macs2 callpeak 会自动进行这一步操作。

Reads去重之后，实验组和对照组reads总数分别是199583 和 199867； 二者之间可以根据比例（199583/199867=.99858）进行操作。

$ macs2 bdgopt -i local_bias_raw.bdg -m multiply -p .99858 -o local_lambda.bdg

参考：
Advanced: Call peaks using MACS2 subcommands
fold enrichment
ppois signal track problem