Abstract Typically, analysis of data consisting of multiple observations on a cluster is complicated by within-cluster correlation. Estimating equations for generalized linear modeling of clustered data have recently received much attention. This article proposes an extension to the generalized estimating equation method proposed by Liang and Zeger, which treats within-cluster correlations as nuisance parameters. Using ideas from extended quasi-likelihood, estimating equations for regression and association parameters are provided simultaneously. The resulting estimators are proven to be asymptotically normal and consistent under certain conditions. The consistency of regression estimators allows incorrect modeling of the correlation among repeated responses. The method is illustrated with an analysis of data from a developmental toxicity study.