Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

paired-end MeRIP-seq #496

Open
BenxiaHu opened this issue Jan 5, 2022 · 4 comments
Open

paired-end MeRIP-seq #496

BenxiaHu opened this issue Jan 5, 2022 · 4 comments

Comments

@BenxiaHu
Copy link

BenxiaHu commented Jan 5, 2022

Hello,
I have paired-end MeRIP-seq and input RNA-seq. the input files are bam files generated by Hisat2. When I ran MACS3 with -f BAM, it works. However, when setting -f to BAMPE, I got the following errors:
INFO @ Wed, 05 Jan 2022 08:14:44: #3 Call peaks...
Traceback (most recent call last):
File "/macs3", line 727, in
main()
File "/macs3", line 48, in main
run( args )
File "/MACS3/Commands/callpeak_cmd.py", line 249, in
run
peakdetect.call_peaks()
File "MACS3/Signal/PeakDetect.pyx", line 112, in MACS3.Signal.PeakDetect.PeakDetect.call_peaks
File "MACS3/Signal/PeakDetect.pyx", line 184, in MACS3.Signal.PeakDetect.PeakDetect.__call_peaks_w_control
AssertionError: 1000 can't be smaller than 2024.6768798828125!

Does it mean that MACS3 analyzes paired-end MeRIP-seq using -f BAM?
Best

Additional context
Add any other context or screenshots about your use case here.

@taoliu
Copy link
Contributor

taoliu commented Jan 7, 2022

The error is caused by the fact that the average fragment length in the sample is 2024bps, and the default setting for estimating background noise with 'small local' (--slocal) of 1000bps is less than 2024. As for MeRIP-seq (or RNA data), I don't feel the estimation of the local noise is necessary, therefore, you can disable it with --nolambda. Hopefully, you won't see such error.

@AZREALH
Copy link

AZREALH commented Nov 7, 2022

The error is caused by the fact that the average fragment length in the sample is 2024bps, and the default setting for estimating background noise with 'small local' (--slocal) of 1000bps is less than 2024. As for MeRIP-seq (or RNA data), I don't feel the estimation of the local noise is necessary, therefore, you can disable it with --nolambda. Hopefully, you won't see such error.

I encounter the same error, and --nolambda can't resolve it. How can I do?

@hlldhy
Copy link

hlldhy commented Mar 15, 2023

The error is caused by the fact that the average fragment length in the sample is 2024bps, and the default setting for estimating background noise with 'small local' (--slocal) of 1000bps is less than 2024. As for MeRIP-seq (or RNA data), I don't feel the estimation of the local noise is necessary, therefore, you can disable it with --nolambda. Hopefully, you won't see such error.

I encounter the same error, and --nolambda can't resolve it. How can I do?

I encounter the same error, and --nolambda can't resolve it. How can I do?

@q1134269149
Copy link

I also encounter the same error. How can I do?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

5 participants