Genetic variance components analysis for binary phenotypes using generalized linear mixed models (GLMMs) and Gibbs sampling