fixed fit stat calculations for NxN model

Added code to calculate k (# of parameters) for AIC, AIC.c, and BIC for
the NxN model. Also made AIC, AIC.c, and BIC elements of fit object
returned by fit.grt() function.