Binomial Likelihood Function

example.py

#
# Example from: http://warnercnr.colostate.edu/~gwhite/fw663/BinomialLikelihood.PDF
#

from math import factorial

# Number of flips
n = 11

# Number of heads
y = 7

# Binomial likelihood function
L = lambda p: C(n, y) *  p**y * (1-p)**(n-y)

# Combinations
C = lambda n, y: factorial(n)/(factorial(y) * factorial(n-y))