Doubly truncated data are found in astronomy, econometrics and survival analysis literature. They arise when an individual is potentially observed only if its failure time lies within a certain interval, unique to that individual. In this paper, we consider the pseudo-value approach to estimate regression coefficients for doubly truncated data. The pseudo-values generated from the nonparametric maximum likelihood estimates (NPMLE) of the survival function are used as response variables in a generalized estimating equation (GEE) to obtain estimates of the model parameters. The performance of the proposed estimators is investigated through a simulation study.