glmmkin function

Fit generalized linear mixed model with known relationship matrices