kgaps_stat function

Sufficient statistics for the KK-gaps model