maxmouchet

def penalty(β):
    """
    L: log-likelihood
    n: number of data points
    k: number of parameters
    """
    return lambda L, n, k: β(n)*k - 2*L

AIC = penalty(lambda n: 2)
BIC = penalty(lambda n: np.log(n))