Negative binomial layer

@johnbb thanks for direction on SpecialFunctions.jl, than it can be written with loggamma fce as:

nll = @. loggamma(n) + loggamma(y + 1) - loggamma(n + y) - n * log( p) - y * log(1 - p)

what is also faster solution.
Lubo