Frailty models are becoming increasingly popular in multivariate survival analysis. Shared frailty models in particular are often used despite their limitations. To overcome their disadvantages numerous correlated frailty models were established during the last decade. In the present study, we examine bivariate correlated frailty models, and especially the behavior of the parameter estimates when using different estimation strategies. We consider three different bivariate frailty models: the gamma model and two versions of the log-normal model. The traditional maximum likelihood procedure of parameter estimation in the gamma case with an explicit available likelihood function is compared with maximum likelihood methods based on numerical integration and a Bayesian approach using MCMC methods with the help of a comprehensive simulation study. We detected a strong dependence between the two parameter estimates (variance and correlation of frailties) in the bivariate correlated frailty model and analyzed this dependence in detail.