[R] unimodal VS bimodal normal distribution - how to get a pvalue?

[Emmanuel Levy]

Dear All,

I have a distribution of values and I would like to assess the
uni/bimodality of the distribution.

I managed to decompose it into two normal distribs using Mclust, and
the BIC criteria is best for two parameters.
However, the problem is that the BIC criteria is not a P-value, which
I would need ideally.

I saw the diptest package but it is not for gaussian distributions.

Any hint at a package or way-around this would be greatly appreciated.

Best,

Emmanuel