Bivariate logistic regression: modelling the association of small for gestational age births in twin gestations

Clustered binary responses, such as disease status in twins, frequently arise in perinatal health and other epidemiologic applications. The scientific objective involves modelling both the marginal mean responses, such as the probability of disease, and the within-cluster association of the multivariate responses. In this regard, bivariate logistic regression is a useful procedure with advantages that include (i) a single maximization of the joint probability distribution of the bivariate binary responses, and (ii) modelling the odds ratio describing the pairwise association between the two binary responses in relation to several covariates. In addition, since the form of the joint distribution of the bivariate binary responses is assumed known, parameters for the regression model can be estimated by the method of maximum likelihood. Hence, statistical inferences may be based on likelihood ratio tests and profile likelihood confidence intervals. We apply bivariate logistic regression to a perinatal database comprising 924 twin foetuses resulting from 462 pregnancies to model obstetric and clinical risk factors for the association of small for gestational age births in twin gestations.