On Sun, 24 Jul 2005, jack wrote: > Sorry Allin, I must admit I wasn't clear at all. If we really want to go > with the others, the following patch should be applied to CVS: > > int k = var->ncoeff; > > var->ll = -(g * n / 2.0) * (LN_2_PI + 1) - (n / 2.0) * var->ldet; > - var->AIC = (-2.0 / g) * var->ll + 2.0 * k; > - var->BIC = (-2.0 / g) * var->ll + k * log(n); > + var->AIC = (-2.0 * var->ll + 2.0 * k * g) / n; > + var->BIC = (-2.0 * var->ll + log(n) * k * g) / n;
I'm willing to do that, but it seems a bit odd to divide the whole thing by the number of observations. Is there an obvious rationale? Allin.
