Joint modeling of success and treatment discontinuation in in vitrofertilization programs: a retrospective cohort study
© Troude et al.; licensee BioMed Central Ltd. 2012
Received: 18 November 2011
Accepted: 27 July 2012
Published: 3 August 2012
As discontinuation in in vitro fertilization (IVF) programs has been associated with a poor prognosis, one hypothesis is that some couple-specific predictive factors in IVF may be shared with opposite effect by both success (i.e. live birth) and treatment discontinuation processes. Our objective was to perform a joint analysis of these two processes to examine the hypothesis of a link between the two processes.
Analyses were conducted on a retrospective cohort of 3,002 women who began IVF between 1998 and 2002 in two French IVF centers: a Parisian center and a center in a medium-sized city in central France. A shared random effects model based on a joint modelization of IVF treatment success and discontinuation was used to study the link between the two processes.
Success and discontinuation processes were significantly linked in the medium-sized city center, whereas they were not linked in the Parisian center. The center influenced risk of treatment discontinuation but not chance of success. The well-known inverse-J relation between the woman’s age and chance of success was observed, as expected. Risk of discontinuation globally increased as the woman’s age increased.
The link between success and discontinuation processes could depend on the fertility center. In particular, the woman’s decision to pursue or to discontinue IVF in a particular center could depend on the presence of other IVF centers in the surrounding area.
Discontinuation of in vitro fertilization (IVF) treatment is common, whatever the technical or social management of patients [1–3]. During recent years, studies conducted in various countries such as the Netherlands, Sweden and the United Kingdom have reported very high rates of discontinuation in IVF programs, with 25% to 50% of couples discontinuing treatment after the first or second attempt [1, 2, 4–6]. In countries where the financial costs of IVF must be mainly borne by the couple, such as the United Kingdom and the United States, the decision to continue or to discontinue IVF treatment is probably strongly influenced by financial issues . In France, where couples receive financial support for IVF treatment (limited to 4 attempts for one pregnancy), discontinuation rates also appear to be high, with estimated cumulative discontinuation rates of 46% to 58% before the fourth attempt [7, 8]. Treatment discontinuation is thus not only a financial matter and may be motivated by the heavy psychological or physical burden of IVF treatment and/or by a poor prognosis [2, 3].
In an English cohort study of 2,056 couples, the characteristics of patients who discontinued IVF treatment after the first attempt were compared with those who had a second IVF attempt . Among couples who discontinued IVF treatment, more women were aged >35 years, had five or less oocytes retrieved at the first attempt and two or less embryos available at the first attempt . All these factors were associated with a lower chance of successful IVF. Other studies also demonstrated common factors associated with success and treatment interruption [4, 6, 9, 10]. Some predictive factors in IVF treatment thus appeared to be shared by both success and treatment discontinuation processes, having an opposite impact on the two processes. For example, as the woman’s age increases, the chance of success decreases and, at the same time, the risk of treatment interruption increases. Such a link between success and treatment interruption is confirmed by studies investigating reasons for IVF treatment discontinuation. Even if the reason for treatment interruption is difficult to assess (it is probably a multifactorial decision, and one that is often studied several years after interruption), these studies showed that one in four couples considered that their treatment interruption was due to a poor prognosis [2, 3]. Such results demonstrate the importance of studying discontinuation and success processes together.
Our aim was to investigate conjointly treatment discontinuation and success processes in two IVF centers and to examine the hypothesis that there is a link between the two processes in each center. We also aimed to explore the effects of the woman’s age and of the IVF center on the success process and on the discontinuation process. For this purpose, a shared random effects model was used.
Design and subjects
The study was conducted in two French IVF centers: a center located in Paris (Cochin) and a center located in a medium-sized city in central France (Clermont-Ferrand).
All women having their first aspiration in one of the two centers between 1998 and 2002 were included in the study (n = 3,037). Thirty-five women were excluded because the result of the first attempt was unknown, leaving 3,002 women. This study received approval from the French Data Protection Authority in September 2005 (authorization number 05-1334).
Information was collected from medical records for all aspirations undergone by the couples in the IVF center, as well as data on frozen embryo transfer (FET), up to 2005. The couples’ characteristics collected included the woman’s date of birth, date of aspiration, number of oocytes retrieved, IVF technique used, number of fresh embryos transferred, number of frozen embryos and the result of transfer (pregnancy, delivery). As the French social security system reimburses IVF treatment up to four aspirations, data collection was discontinued after the woman’s fourth aspiration.
The success of IVF was measured by a live birth after one attempt. The live birth rate was defined as a delivery resulting from fresh or frozen embryo transfer among women who had undergone one IVF attempt. Treatment discontinuation was defined as no treatment for at least two years in the IVF center (whatever the reason for discontinuation, e.g. maternal age, financial resources, move to another area, seeking IVF treatment elsewhere…) [8, 11, 12]. Discontinuation rate was defined as discontinuation among women who had not obtained a live birth after the IVF attempt.
Women’s characteristics at the first attempt were compared according to inclusion center, using the chi 2 test. Live birth rate and discontinuation rate at each attempt were also compared according to inclusion center.
Shared random effects model
These models included baseline factors ( and ). They included two observed determinants called “fixed effects”:
the IVF center (Parisian center/medium-sized city center) denoted and
and the woman’s age at the first attempt divided into five classes (age <25/25-29/30-34/35-39/≥40 years) denoted and .
It also included a “random effect” (f i ) that is common to both mixed logistic regression models. This shared “random effect” represents all the couple-specific and non-explicitly identified factors explaining both processes. The coefficient associated with this shared random effect (λ), could be interpreted as the “link” between the two processes: if λ equals zero, it means that there are no couple-specific factors (other than the observed covariates included in the model) that simultaneously explain IVF success and treatment discontinuation. On the contrary, if λ differs from zero, the interpretation is that there are couple-specific non-observed factors that determine both IVF success and treatment discontinuation. When λ is strictly negative, it means that these factors have an opposite relative impact on the two processes. In our model, a specific λ was considered for each center (λParis and λMedium-sized city) in order to allow different degrees of “links” between the two processes in the two IVF centers. Detailed information on this model is given in the Additional file 1.
The shared random effects model was fitted with Bayesian computational methods using Markov Chain Monte Carlo (MCMC)  and implemented in WinBUGS . We ran two independent MCMC chains (using different initial values for the parameters) of 300,000 simulations with a burn-in period of 50,000 and kept every 100th to reduce autocorrelation in the MCMC samples. Our results are therefore based on thinned samples of size 5,000. Convergence of the MCMC run was assessed by graphical inspection of the chains and by computing the Gelman-Rubin statistics as modified by Brooks and Gelman [17, 18] and intra-chains autocorrelations. Credible intervals (95%) were estimated. The λParis and λMedium-sized city parameters were tested according to Bayesian statistical theory methods using partial Bayes factors (BF) : no evidence that λ ≠ 0 if L = 2logBF is greater than -2, λ significantly different from 0 if L ≤ -2 and strongly different from 0 if L ≤ -6.
We performed a Bayes factor sensitivity analysis to prior choice of parameter distribution  and checked the stability of our Bayesian method . We computed partial Bayes factors from more or less informative priors (and only one MCMC run) by performing a split test sample analysis as described elsewhere . We considered three ways to split the sample of 3,002 women so as to define a learning sample (i.e., 2,702, 1,502 and 502 women) that would be more or less informative on the remaining test sample (i.e., 300, 1,500, 2,500 women).
Characteristics of study population at the first aspiration ( N = 3,002)
Parisian center (Cochin)
Medium-sized city center (Clermont-Ferrand)
Patient age (years)
Fresh embryos transferred
Observed live birth and discontinuation rates in the study population ( N = 3,002)
Parisian center (Cochin)
Medium-sized city center (Clermont-Ferrand)
Live birth rate
Estimations of the links between the success and treatment discontinuation processes for the Parisian and the Medium-sized city center
Posterior mean (posterior standard deviation)a
95% credible interval
Test of λb
L = 2.70 c - NS
λ Medium-sized city
L = -5.26 d - S
To study jointly success and treatment discontinuation in IVF programs, a shared random effects model was built and used to analyze data from two French IVF centers: a Parisian center and a medium-sized city center. We found no evidence of a link between success and discontinuation processes in the Parisian center, whereas we did find one in the medium-sized city center. In the medium-sized city center, the negative link observed between the two processes meant that women who discontinued treatment in this center had a lower probability of success. The direction of the link was expected and is in agreement with the literature, as previous studies have reported poorer prognostic factors of IVF success among women who discontinued treatment [1, 4, 10]. However, it is interesting that no evidence of such a link was observed in the Parisian center.
Woman’s age: effect on success and treatment discontinuation
By defining success as a live birth during the entire IVF program (first to fourth IVF attempts), we demonstrated an inverse-J relationship between the woman’s age and success. Such a relationship has already been shown . We also found that the probability of discontinuation varied inversely to the probability of success according to the woman’s age, and that the probability of discontinuation globally increased with increasing female age. In our study, using women aged 30-34 years as reference, OR of treatment discontinuation was 2.9 ([2.2-3.7]) for women older than 40 years. A few studies have already shown that women who discontinued treatment were generally older that women who persevered . In our study, we assessed how the cumulative risk of discontinuation alters with the woman’s age and we demonstrated a J-relationship between the woman’s age and cumulative risk of discontinuation during an IVF program.
IVF center: variability in success and treatment discontinuation
In our study, the crude cumulative success rate differed according to center (34% vs 41%), but after controlling for treatment discontinuation and the woman’s age, the probability of success no longer differed between centers (Figure 2). This result is similar to that of a study carried out in two centers in the Netherlands . Investigators found that the crude cumulative live birth rates differed between the two centers. However, this difference was not due to differences in success rates at each attempt, but rather to different discontinuation rates in each center. These observations showed that comparison of IVF centers should not be done on crude success rate and that treatment discontinuation is an important factor that should be taken into account.
Conversely, we found that treatment discontinuation rates differed between centers, being lower in the medium-sized city center. Moreover, a negative link was found between success and discontinuation in the medium-sized city center, whereas there was no evidence of such a link in the Parisian center. One explanation could be that the characteristics of women who discontinued treatment differed between the two centers. In the medium-sized city center, the negative link meant that women who discontinued had poorer prognostic factors. In the Parisian center, there was no reason why women with poorer prognostic factors did not discontinue treatment, but the higher level of discontinuation could indicate that women with good prognostic factors also discontinued IVF treatment in the Parisian center. This hypothesis of a more mixed population could explain the lack of significant link between success and discontinuation in the Parisian center. Differences between fertility centers may be linked to various factors such as patient selection, medical staff, or management practices (i.e. choice of IVF vs ICSI, number of embryos transferred). However, one major difference between the two centers in our study is their geographical environment: the Parisian center is surrounded by 23 other IVF centers (9 in Paris itself and 14 in the suburbs), whereas the medium-sized city center is the only one in this administrative area and the nearest other center is in the city of Lyon, a 2-hour drive away. Consequently, the medium-sized city center could be defined as a monopoly center, whereas the Parisian one competes with several other fertility centers. When there are several fertility centers close to the woman’s place of residence (a competition situation between centers), the population of women who discontinue is probably mixed, consisting of both patients with a poorer prognosis and patients who merely change IVF center, whatever their prognosis. On the contrary, in a monopoly center, as women cannot easily discontinue in order to begin another IVF program elsewhere, most treatment discontinuations are linked to poorer prognostic factors.
To the best of our knowledge, this is the first time that the hypothesis of an association between treatment discontinuation and a monopoly/competition situation of the IVF center has emerged in the literature. However, an association between success and a monopoly situation or competition between IVF centers has already been considered. Indeed, some studies have tested the association between an increasing number of multiple pregnancies and competition between IVF centers, the underlying hypothesis being that a greater number of embryos are transferred in centers that compete against various others, in order to maximize the chance of success [24, 25]. Recently, a large American study, conducted in clinics performing ART between 1995 and 2001 (n = 2374 clinic-years), has examined the relationship between competition and clinic-level ART outcomes and practice patterns . Defining competition as the number of clinics within a 20-mile radius (32.19 km) of a given clinic, they found no evidence of a significant relationship between competition and birth rates in multivariate models. Moreover, they found a lower, rather than a higher, rate of multiple births per ART cycle for clinics in highly competitive areas, as has been suggested in one previous study using another definition for competition . Our results are in agreement with the American study, showing no difference between the two centers, one being in a monopoly situation and the other in a competition situation, with regard to chance of success.
In our model, we included only female age and center as fixed effects. In the context of growing interest in understanding differences between IVF success rates according to center [28, 29], some studies have explored to what extent such differences may be linked to differences in patients’ characteristics. An English study has explored the influence of patients’ characteristics on live birth rate per cycle started . The authors demonstrated the impact of non-IVF related patient characteristics on the success rate and concluded that using a “standard patient group and outcome” did not improve validity of comparisons between centers. More recently, using IVF and ICSI treatment data from 11 IVF centers in the Netherlands, Lintsen et al. studied how differences in IVF success rates between centers could be explained by patient characteristics and concluded that only 17% of the variation between centers could be explained by patient mix . Thus, there is currently no clear evidence that other patient characteristics should be taken into account in our multivariate model. However, our shared random effects model could be extended by including temporal effects that could describe, for instance, the patient’s level of discouragement due to psychological and physical burden. Such a temporal effect could also be included in the model to test if the link between the success and treatment discontinuation processes may also depend on the IVF attempt. Obviously, it would be of great interest to conduct such analysis on a greater number of centers to better understand how the center’s situation impacts on treatment discontinuation.
Despite the increasing interest in understanding differences in IVF success rates between IVF centers, the reasons explaining such differences remain rather unclear. It is likely that differences in IVF centers success rates are a combination of patient and center characteristics . Treatment discontinuation rate could be one of the factors impacting on the center success rate but it has scarcely been investigated. In our study, we observed two French centers with different crude success rates. After controlling for the woman’s age and for the impact of discontinuation on success rate in a shared random effects model, success rates between the two centers no longer differed. Our results showed that discontinuation may be a very important factor in explaining success rate differences between centers, and it needs to be better understood.
Our study also enabled us to explore treatment discontinuation. Our main result was that discontinuation appeared very dissimilar in the two centers. The center strongly influenced the risk of treatment discontinuation (unlike the chance of success). An important perspective of this work will be to explore further the discontinuation process and differences between centers in a larger number of centers. Based on our results, a very promising hypothesis would be to explore the possible influence of the IVF center situation (monopoly/non-monopoly) on the probability of treatment discontinuation. Our hypothesis is that the probability of treatment discontinuation decreased in centers that were in a monopoly situation.
Based on IVF treatment, we hypothesize that a better understanding of the treatment discontinuation process in relation to competition may be very helpful in understanding differences in success rates between IVF centers. It would be very interesting to explore if a similar hypothesis would be pertinent for treatments other than IVF, when patients have the possibility of deciding to pursue their treatment in another medical center.
In vitro fertilization
Acknowledgments and funding
We thank Béatrice Ducot for helpful comments on an earlier version of the manuscript. This project was supported by the Agence Nationale de la Recherche [n°ANR-06-BLAN-0221-01] and by a post-doctoral grant from the Région Ile-de-France and Paris-Sud XI University [n°F-08-1220/R].
- Sharma V, Allgar V, Rajkhowa M: Factors influencing the cumulative conception rate and discontinuation of in vitro fertilization treatment for infertility. Fertil Steril. 2002, 78 (1): 40-46. 10.1016/S0015-0282(02)03160-6.View ArticlePubMedGoogle Scholar
- Olivius C, Friden B, Borg G, Bergh C: Why do couples discontinue in vitro fertilization treatment? a cohort study. Fertil Steril. 2004, 81 (2): 258-261. 10.1016/j.fertnstert.2003.06.029.View ArticlePubMedGoogle Scholar
- Rajkhowa M, McConnell A, Thomas GE: Reasons for discontinuation of IVF treatment: a questionnaire study. Hum Reprod. 2006, 21 (2): 358-363.View ArticlePubMedGoogle Scholar
- Land JA, Courtar DA, Evers JL: Patient dropout in an assisted reproductive technology program: implications for pregnancy rates. Fertil Steril. 1997, 68 (2): 278-281. 10.1016/S0015-0282(97)81515-4.View ArticlePubMedGoogle Scholar
- Olivius K, Friden B, Lundin K, Bergh C: Cumulative probability of live birth after three in vitro fertilization/intracytoplasmic sperm injection cycles. Fertil Steril. 2002, 77 (3): 505-510. 10.1016/S0015-0282(01)03217-4.View ArticlePubMedGoogle Scholar
- Malizia BA, Hacker MR, Penzias AS: Cumulative live-birth rates after in vitro fertilization. N Engl J Med. 2009, 360 (3): 236-243. 10.1056/NEJMoa0803072.View ArticlePubMedGoogle Scholar
- de La Rochebrochard E, Soullier N, Peikrishvili R, Guibert J, Bouyer J: High in vitro fertilization discontinuation rate in France. Int J Gynaecol Obstet. 2008, 103 (1): 74-75. 10.1016/j.ijgo.2008.05.009.View ArticlePubMedGoogle Scholar
- Soullier N, Bouyer J, Pouly JL, Guibert J, de La Rochebrochard E: Estimating the success of an in vitro fertilization programme using multiple imputation. Hum Reprod. 2008, 23 (1): 187-192.View ArticlePubMedGoogle Scholar
- Stolwijk AM, Wetzels AM, Braat DD: Cumulative probability of achieving an ongoing pregnancy after in-vitro fertilization and intracytoplasmic sperm injection according to a woman's age, subfertility diagnosis and primary or secondary subfertility. Hum Reprod. 2000, 15 (1): 203-209. 10.1093/humrep/15.1.203.View ArticlePubMedGoogle Scholar
- Verberg MF, Eijkemans MJ, Heijnen EM, Broekmans FJ, de Klerk C, Fauser BC, Macklon NS: Why do couples drop-out from IVF treatment? A prospective cohort study. Hum Reprod. 2008, 23 (9): 2050-2055. 10.1093/humrep/den219.View ArticlePubMedGoogle Scholar
- de La Rochebrochard E, Quelen C, Peikrishvili R, Guibert J, Bouyer J: Long-term outcome of parenthood project during in vitro fertilization and after discontinuation of unsuccessful in vitro fertilization. Fertil Steril. 2009, 92 (1): 149-156. 10.1016/j.fertnstert.2008.05.067.View ArticlePubMedGoogle Scholar
- Soullier N, Bouyer J, Pouly JL, Guibert J, de La Rochebrochard E: Effect of the woman's age on discontinuation of IVF treatment. Reprod Biomed Online. 2011, 22 (5): 496-500. 10.1016/j.rbmo.2011.01.013.View ArticlePubMedGoogle Scholar
- Gao S: A shared random effect parameter approach for longitudinal dementia data with non-ignorable missing data. Stat Med. 2004, 23 (2): 211-219. 10.1002/sim.1710.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu MC, Carroll RJ: Estimation and comparison of changes in the presence of informative right censoring by modeling the censoring process. Biometrics. 1988, 44 (1): 175-188. 10.2307/2531905.View ArticleGoogle Scholar
- Robert CP, Casella G: Monte Carlo Statistical Methods. 2004, New York: Springer, 2View ArticleGoogle Scholar
- Lunn DJ, Thomas A, Best N, Spiegelhalter D: WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility. Stat Comput. 2000, 10 (4): 325-337. 10.1023/A:1008929526011.View ArticleGoogle Scholar
- Gelman A: Prior distributions for variance parameters in hierarchical models. Bayesian Analysis. 2006, 1 (3): 515-533.Google Scholar
- Brooks SP, Gelman A: Alternative methods for monitoring convergence of iterative simulations. J Comput Graph Stat. 1998, 7: 434-455.Google Scholar
- Kass RE, Raftery AE: Bayes Factors. J Am Stat Assoc. 1995, 90 (430): 773-795. 10.1080/01621459.1995.10476572.View ArticleGoogle Scholar
- Ancelet S, Etienne M-P, Benoît H, Parent E: Modelling spatial zero-inflated continuous data with an exponentially compound Poisson process. Environ Ecol Stat. 2010, 17: 347-376. 10.1007/s10651-009-0111-6.View ArticleGoogle Scholar
- Templeton A, Morris JK, Parslow W: Factors that affect outcome of in-vitro fertilisation treatment. Lancet. 1996, 348 (9039): 1402-1406. 10.1016/S0140-6736(96)05291-9.View ArticlePubMedGoogle Scholar
- Pearson KR, Hauser R, Cramer DW, Missmer SA: Point of failure as a predictor of in vitro fertilization treatment discontinuation. Fertil Steril. 2009, 91 (4 Suppl): 1483-1485.View ArticlePubMedGoogle Scholar
- Witsenburg C, Dieben S, Van der Westerlaken L, Verburg H, Naaktgeboren N: Cumulative live birth rates in cohorts of patients treated with in vitro fertilization or intracytoplasmic sperm injection. Fertil Steril. 2005, 84 (1): 99-107. 10.1016/j.fertnstert.2005.02.013.View ArticlePubMedGoogle Scholar
- Fisk NM, Trew G: Two's company, three's a crowd for embryo transfer. Lancet. 1999, 354 (9190): 1572-1573. 10.1016/S0140-6736(99)00290-1.View ArticlePubMedGoogle Scholar
- Jones HW, Schnorr JA: Multiple pregnancies: a call for action. Fertil Steril. 2001, 75 (1): 11-13. 10.1016/S0015-0282(00)01612-5.View ArticlePubMedGoogle Scholar
- Henne MB, Bundorf MK: The effects of competition on assisted reproductive technology outcomes. Fertil Steril. 2010, 93 (6): 1820-1830. 10.1016/j.fertnstert.2008.02.159.View ArticlePubMedGoogle Scholar
- Steiner AZ, Paulson RJ, Hartmann KE: Effects of competition among fertility centers on pregnancy and high-order multiple gestation rates. Fertil Steril. 2005, 83 (5): 1429-1434. 10.1016/j.fertnstert.2004.10.048.View ArticlePubMedGoogle Scholar
- Lass A, Brinsden P: How do patients choose private in vitro fertilization treatment? A customer survey in a tertiary fertility center in the United Kingdom. Fertil Steril. 2001, 75 (5): 893-897. 10.1016/S0015-0282(01)01751-4.View ArticlePubMedGoogle Scholar
- Marcus HJ, Marcus DM, Marcus SF: How do infertile couples choose their IVF centers? An Internet-based survey. Fertil Steril. 2005, 83 (3): 779-781. 10.1016/j.fertnstert.2004.11.003.View ArticlePubMedGoogle Scholar
- Johnson A, El-Toukhy T, Sunkara SK, Khairy M, Coomarasamy A, Ross C, Bora S, Khalaf Y, Braude P: Validity of the in vitro fertilisation league tables: influence of patients' characteristics. BJOG. 2007, 114 (12): 1569-1574. 10.1111/j.1471-0528.2007.01539.x.View ArticlePubMedGoogle Scholar
- Lintsen AM, Braat DD, Habbema JD, Kremer JA, Eijkemans MJ: Can differences in IVF success rates between centres be explained by patient characteristics and sample size?. Hum Reprod. 2010, 25 (1): 110-117. 10.1093/humrep/dep358.View ArticlePubMedGoogle Scholar
- Sharif K, Afnan M: The IVF league tables: time for a reality check. Hum Reprod. 2003, 18 (3): 483-485. 10.1093/humrep/deg134.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2393/12/77/prepub