首页 | 本学科首页   官方微博 | 高级检索  
     检索      


Survival analysis for a new compounded bivariate failure time distribution in shock and competing risk models via an EM algorithm
Authors:Shirin Shoaee  Esmaile Khorram
Institution:1. Department of Actuarial Science, Faculty of Mathematical Sciences, Shahid Beheshti University, Tehran, Iran;2. Sh_Shoaee@sbu.ac.ir;4. Department of Statistics, Faculty of Mathematics and Computer Science, Amirkabir University of Technology, Tehran, Iran
Abstract:Abstract

In this paper, two bivariate models based on the proposed methods of Marshall and Olkin are introduced. In the first model, the new bivariate distribution is presented based on the proposed method of Marshall and Olkin (1967 Marshall, A. W., and I. Olkin. 1967. A multivariate exponential distribution. Journal of the American Statistical Association 62 (317):3044. doi: 10.2307/2282907.Taylor & Francis Online], Web of Science ®] Google Scholar]) which has natural interpretations, and it can be applied in fatal shock models or in competing risks models. In the second model, the proposed method of Marshall and Olkin (1997 Marshall, A. W., and I. Olkin. 1997. A new method of adding a parameter to a family of distributions with application to the exponential and weibull families. Biometrika 84 (3):64152. doi: 10.1093/biomet/84.3.641.Crossref], Web of Science ®] Google Scholar]) is generalized to bivariate case and a new bivariate distribution is introduced. We call these new distributions as the bivariate Gompertz (BGP) distribution and bivariate Gompertz-geometric (BGPG) distribution, respectively. Moreover, the BGP model can be obtained as a special case of the BGPG model. Then, we present various properties of the new bivariate models. In this regard, the joint and conditional density functions, the joint cumulative distribution function can be obtained in compact forms. Also, the aging properties and the bivariate hazard gradient are discussed. This model has five unknown parameters and the maximum likelihood estimators cannot be obtained in explicit form. We propose to use the EM algorithm to compute the maximum likelihood estimators of the unknown parameters, and it is computationally quite tractable. Also, Monte Carlo simulations are performed to investigate the effectiveness of the proposed algorithm. Finally, we analyze three real data sets for illustrative purposes.
Keywords:Bivariate model  competing risks model  expectation-maximization algorithm  gompertz distribution  shock model
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号