lh3 / psmc

Implementation of the Pairwise Sequentially Markovian Coalescent (PSMC) model
Other
146 stars 60 forks source link

hmm error #51

Closed ndussex closed 1 year ago

ndussex commented 1 year ago

Hi,

I am running the final step of the psmc using a pseudo-diploid chromosome such as:

psmc -N25 -t30 -r5 -d -p "4+25*2+4+6" -o Pseudo_ChrX.psmc Pseudo_ChrX.psmcfa

and I get the following error, which I have never seen so far:

psmc: khmm.c:278: hmm_post_decode: Assertion `max_k >= 0' failed

Has anyone seen this before?

Cheers, Nic