distributions.py 文件源码

def poisson_logpdf(k,l):
    '''
    Gives the log-pdf for a poisson distribution with rate l 
    evaluated at points k. k should be a vector of integers.

    Parameters
    ----------

    Returns
    -------
    '''
    # k,l = map(np.float128,(k,l))
    return k*slog(l)-l-np.array([scipy.special.gammaln(x+1) for x in k])
    #return k*slog(l)-l-np.array([log_factorial(x) for x in k])