Dear Konrad,

The last component is the omega part (>=1, also notice different weight for it). So that is fine. The other values > 1 suggest that there is something funny going on with discretization. Could you post the P and Q values for this beta?

HTH,
Sergei

Dear Konrad,

The discretization procedure fails to perform correctly because of the steep density gradient around 1. I'll fix the code (it needs to be in done in HyPhy source), and post a note here when it's done. In the meantime, try setting the minimum betaP and betaQ to 0.1

Cheers,
Sergei

Dear Konrad,

I have looked into how to resolve this discretization issue; unfortunately, there is not much that the current implementation lets you do - it discretizes the mixture, and because the gradient of the beta gets very steep around '1', the double precision is insufficient to resolve discretization intervals - one of them (incorrectly) spills over one, and you get silly values for the beta part. Note that if you remove the omega part, the beta distribution with same p and q will be discretized correctly, because of the [0,1] support coded into the definition of the category variable.

I'll actually add a new way to define mixtures (e.g. category mixed = (cat1,p1,cat2,p2,...,catN,(p_n)), which will allow to specify "hard" bounds on the components and avoid round-off error related issues.

In the meantime, I suggest inspecting the output for oddities like the one you found - this is a fairly rare case, and it suggests that beta+w is a poor choice to model the distribution of rates in your case.

Thanks for bringing this to my attention and stay tuned.

Cheers,
Sergei