'Binomial variance function
Parameters
mu : array-like
mean parameters
Returns
variance : array
variance = mu/n * (1 - mu/n) * self.n'

def __call__(self, mu):

    p = self._clean((mu / self.n))
    return ((p * (1 - p)) * self.n)
