kr__modcomp function

F-test and degrees of freedom based on Kenward-Roger approximation