PSMC Illegal division by zero
1
1
Entering edit mode
19 months ago
yoce_pf ▴ 60

Dear all,

I'm using PSMC for many genomes, and only for two genomes, I received an error after performing the command psmc2history.pl

Error:
Illegal division by zero at /PATH/psmc_plot.pl line 119


I performed 100 PSMC bootstrap replicates and I concatenated these files to run the command: psmc2history.pl

I have seen the same error in other groups, but they have not figured out a solution.

I'm wondering if you could help me to fix and understand this error.

Best,
Yocelyn

genome software error sequence • 1.1k views
1
Entering edit mode

Hi,

I am having the same error. I was wondering if you found a solution to this problem? Thank you, Homa

0
Entering edit mode

@ _r_am have you managed to get the solution? I am getting the same error.

0
Entering edit mode

Same issue. 2/40 genomes failed. Any progress on a solution?

0
Entering edit mode

same error here, any solution?

0
Entering edit mode
4 months ago
Darija • 0

Do your .psmc bootstraped files have the same size. I had the same issue and the problem was that my .psmc files had different sizes, after I re ran the bootstrap and got the same sizes everything was ok. Also, when concatenating, be careful to include the correct files.. that was also one of the errors I had.

0
Entering edit mode

Can you update the code for performing bootstrap with PSMC? I'm having a problem in this line:

seq 100 | xargs -i echo psmc -N25 -t15 -r5 -b -p "4+25*2+4+6" \
-o round-{}.psmc split.fa | sh