added functions to compute goodness of fit for C-stat