-
Notifications
You must be signed in to change notification settings - Fork 270
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
Comments
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 |
I encounter the same error, and |
I encounter the same error, and --nolambda can't resolve it. How can I do? |
I also encounter the same error. How can I do? |
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.
The text was updated successfully, but these errors were encountered: