h_cullis function

Cullis heritability for lme4 models