Inference for threshold models with variance components from the generalized linear mixed model perspective