def AIC_chisq(t_dof, t_chisq):
    """
        t_dof: int
            Number of degrees of freedom i.e. number of parameters of the model
        t_chisq: float
            \u03C7\u00B2 of the model.

        Retruns: float
            AIC

        Assumptions:
            * normally distributed errors

        This function computes
        $$
            AIC = \chi^2 + 2k
        $$

        Following Sz. Borsanyi et. al.
            doi.org/10.1126/science.1257050
    """

    return t_chisq + 2*t_dof