Skip to content

Instantly share code, notes, and snippets.

Embed
What would you like to do?
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))
HQC = penalty(lambda n: 2*np.log(np.log(n)))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment