fit_model function

Goodness of fit statistics for transition matrix data