Li Las
05/07/2025, 7:49 AMHello @channel, does anybody have an idea about this error please ? Caused by:
Process `pureclip_peak_call (PTBP1iCLIP_siMATR3_rep2_s)` terminated with an error exit status (134)
Command executed:
pureclip \
-i PTBP1iCLIP_siMATR3_rep2_s.dedup.bam \
-bai PTBP1iCLIP_siMATR3_rep2_s.dedup.bam.bai \
-g GRCh38.p13.genome.fa \
-nt 20 \
-bc 0 -dm 8 -iv 'all' \
-o "PTBP1iCLIP_siMATR3_rep2_s.sigxl.bed" \
-or "PTBP1iCLIP_siMATR3_rep2_s.8nt.peaks.bed"
pigz PTBP1iCLIP_siMATR3_rep2_s.sigxl.bed PTBP1iCLIP_siMATR3_rep2_s.8nt.peaks.bed
Command exit status:
134
Command output:
b0:0
mean:1
k:1
tp:0.00807885
initial transition frequencies:
0: 1 1 1 1
1: 1 1 1 1
2: 1 1 1 1
3: 1 1 1 1
Build HMM ...
*** Transition probabilitites ***
0: 0.25 0.25 0.25 0.25
1: 0.25 0.25 0.25 0.25
2: 0.25 0.25 0.25 0.25
3: 0.25 0.25 0.25 0.25
Baum-Welch ...
learn binomial parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
**** Convergence ! ****
learn gamma parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
**** Convergence ! ****
learn binomial parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
**** Convergence ! ****
learn gamma parameter
.. 0th iteration
computeEmissionProbs()
computeStatePosteriorsFB()
updateDensityParams()
**** Convergence ! ****
Apply learned parameters to whole dataset ...
chr1
Get preprocessed intervals:
Compute KDEs ...
Estiamte Ns ... 1
Apply learned HMM:
applyParameters
Command error:
INFO: Environment variable SINGULARITYENV_TMPDIR is set, but APPTAINERENV_TMPDIR is preferred
INFO: Environment variable SINGULARITYENV_NXF_DEBUG is set, but APPTAINERENV_NXF_DEBUG is preferred
ERROR: Can't find contigName in store.contigNameStore!all
terminate
called after throwing an instance of
'boost::exception_detail::clone_impl<boost::exception_detail::error_info_injector<std::domain_error>
>'
what(): Error in function
boost::math::binomial_distribution<long
double>::binomial_distribution: Success fraction argument is -nan,
but must be >= 0 and <= 1 !
.command.sh: line 9: 226464
Aborted (core dumped) pureclip -i
PTBP1iCLIP_siMATR3_rep2_s.dedup.bam -bai
PTBP1iCLIP_siMATR3_rep2_s.dedup.bam.bai -g GRCh38.p13.genome.fa -nt 20
-bc 0 -dm 8 -iv 'all' -o "PTBP1iCLIP_siMATR3_rep2_s.sigxl.bed" -or
"PTBP1iCLIP_siMATR3_rep2_s.8nt.peaks.bed"