Abstract A maximum likelihood estimation procedure is developed for a bivariate frailty regression model, in which dependence is generated by a gamma distribution. It is assumed that the random lifetimes follow a bivariate Weibull distribution proposed in Hanagal [Economic Quality Control 11:...