'Scores ngrams using the Poisson-Stirling measure.'

@classmethod
def poisson_stirling(cls, *marginals):

    exp = (_product(marginals[UNIGRAMS]) / (marginals[TOTAL] ** (cls._n - 1)))
    return (marginals[NGRAM] * (_log2((marginals[NGRAM] / exp)) - 1))
