Question regarding the support of the Geometric distribution in Julia

Hi all,

I am evaluating a likelihood involving the Geometric distribution, with a binary latent variable which indicates whether to set the parameter value p of the Geometric distribution to be equal to p=1.

As these terms when p=1 always evaluate to 1, it is not a strong requirement that I have to use something like logpdf to evaluate these observations, but I don’t understand why the support for the Geometric distribution within the Distributions package is 0<p<1 opposed to 0<p\leq1.

Thanks in advance.

Does it work if you set check_args=false?

I briefly looked at the source code and saw this:

But haven’t tested.