Modeling heterogeneity for bivariate survival data by the compound Poisson distribution with random scale

We propose a bivariate Weibull regression model with heterogeneity (frailty or random effect) which is generated by compound Poisson distribution with random scale. We assume that the bivariate survival data follow bivariate Weibull of Hanagal (2004). There are some interesting situations like survival times in genetic epidemiology, dental implants of patients and twin births (both monozygotic and dizygotic) where genetic behavior (which is unknown and random) of patients follows a known frailty distribution. These are the situations which motivate us to study this particular model. We propose a two stage maximum likelihood estimation procedure for the parameters in the proposed model and develop large sample tests for testing significance of regression parameters.