def AKICc(N, rho, k):
from numpy import log, array<EOL>p = k<EOL>res = log(rho) + p/N/(N-p) + (<NUM_LIT>-(p+<NUM_LIT>)/N) * (p+<NUM_LIT:1.>) / (N-p-<NUM_LIT>)<EOL>return res<EOL>
r"""approximate corrected Kullback information

    .. math:: AKICc(k) = log(rho_k) + \frac{p}{N*(N-k)} + (3-\frac{k+2}{N})*\frac{k+1}{N-k-2}