sigma_a computation for binomial

[seekshreyas]

Hi @cbellei,

Thanks for an amazing load of work here. Really impressive.

While going over the code I saw this:

sigma_a_rvs = np.sqrt(a_rvs * (1 - a_rvs))
sigma_b_rvs = np.sqrt(b_rvs * (1 - b_rvs))

at: https://github.com/cbellei/abyes/blob/master/abyes/ab_exp.py#L148

Shouldn't the n_a in the formula as well?
As for binomial distribution
standard deviation (σx) is sqrt[ n * P * ( 1 - P ) ].

Thanks in advance. I'm happy to send a PR, but wanted to check first if I'm correct.

[cbellei]

Hello @seekshreyas , thanks for the nice words!

I haven't looked at this in a long time, but this should be the standard deviation for a Bernoulli distribution so I think it is ok. The binomial distribution is used in the examples in order to obtain one particular realization of a series of experiments drawn from a Bernoulli distribution.