In biometrical genetic analyses of binary traits, the use of family data overcomes some limitations of twin studies, particularly in terms of sample size and types of genetic or environmental factors that can be estimated. However, because of computational problems, recent methods in the application of generalized linear mixed models for family data structure have limited the ability to handle large data sets with general covariates. In this paper, we investigate the use of the hierarchical likelihood approach to the analysis of binary traits from family data. In a simulation study, the method is shown to be highly accurate for the estimation of both the variance components and fixed regression parameters, even for small family sizes. For illustration, we analyze a real data set of familial aggregation of preeclampsia, a pregnancy-induced hypertension. When possible, the analysis is compared with the exact maximum likelihood approach. Genet. Epidemiol. 30: 37-47, 2006. r 2005 Wiley-Liss, Inc.