Skip to main content

Estimating the causal effect of embryo transfer day on clinical in vitro fertilization outcomes using propensity score matching



For women undergoing in vitro fertilization (IVF), the clinical benefit of embryo transfer at the blastocyst stage (Day 5) versus cleavage stage (Day 3) remains controversial. The purpose of this study is to compare the implantation rate, clinical pregnancy rate and odds of live birth of Day 3 and Day 5 embryo transfer, and more importantly, to address the issue that patients were chosen to receive either transfer protocol due to their underlying clinical characteristics, i.e., confounding by indication.


We conducted a retrospective cohort study of 9,090 IVF cycles collected by Lee Women’s Hospital in Taichung, Taiwan from 1998 to 2014. We utilized the method of propensity score matching to mimic a randomized controlled trial (RCT) where each patient with Day 5 transfer was matched by another patient with Day 3 transfer with respect to other clinical characteristics. Implantation rate, clinical pregnancy rate, and odds of live birth were compared for women underwent Day 5 transfer and Day 3 transfer to estimate the causal effects. We further investigated the causal effects in subgroups by stratifying age and anti-Mullerian hormone (AMH).


Our analyses uncovered an evidence of a significant difference in implantation rate (p=0.04) favoring Day 5 transfer, and showed that Day 3 and Day 5 transfers made no difference in both odds of live birth (p=0.27) and clinical pregnancy rate (p=0.11). With the increase of gestational age, the trend toward non-significance of embryo transfer day in our result appeared to be consistent for subgroups stratified by age and AMH, while all analyses stratified by age and AMH were not statistically significant.


We conclude that for women without strong indications for Day 3 or Day 5 transfer, there is a small significant difference in implantation rate in favor of Day 5 transfer. However, the two protocols have indistinguishable outcomes on odds of live birth and clinical pregnancy rate.

Peer Review reports


There has been a long debate on whether the blastocyst stage (Day 5) embryo transfer should be routinely recommended to improve the implantation rate. [13] The introduction and improvements of advanced cell culture techniques have suggested blastocyst stage embryo transfer over cleavage stage (Day 3) embryo transfer for in vitro fertilization (IVF). Many published studies have suggested that Day 5 transfer provides higher implantation, pregnancy, and live birth rates [46]. However, others have argued that patients undergoing blastocyst culture are expected to have fewer transferable embryos available, and fewer embryos cryopreserved [7, 8]. Although several studies including randomized controlled trials (RCTs) have already been conducted, many of which were underpowered, leaving the margin of benefit between Day 3 and Day 5 transfers unclear.

Observational studies without randomization serve as an alternative study design to collect more participants to increase statistical power. However, the large sample size in observational studies pays the price of potentially introducing bias into the study. For example, Day 3 and Day 5 transfers were usually prescribed for patients with different clinical characteristics [9]. Without proper randomization, the difference in live birth rate or other associated outcomes is not necessarily attributable to the day of transfer. The phenomenon represents a typical confounding by indication [10]. The indication of receiving a certain treatment may affect the outcome of interest and thus exerts an undue influence on the association of the treatment and the outcome. Therefore, a powerful causal inference approach is imperative to determine the effectiveness of cleavage stage embryo transfer versus blastocyst stage embryo transfer.

While RCTs are considered as the gold standard approach for estimating the causal effects of treatments or exposures on outcomes, they are largely limited by ethical considerations and other practical limitations. More importantly, the practical constraints of RCTs make it difficult to distinguish between a null effect and a promising finding limited by the statistical power due to its small sample size. To integrate the advantages of RCT and observational study, here we adopt a propensity score matching approach. “Propensity score matching” section enables us to artificially construct an RCT from an observational study under certain assumptions [11]. The notion of propensity score was first introduced to be the probability of receiving a particular treatment conditional on observed baseline variables [12]. For the past few decades, different propensity score methods have been developed to remove the confounding effects under the potential outcome framework [1315]. “Propensity score matching” section entails forming matched pairs of treated and untreated subjects that share similar value of propensity scores. Once the matched sample has been formed, the causal effect of the treatment on an outcome can be easily accessed by standard statistical procedures. To answer our scientific question by implementing this particular causal method, we had dichotomous embryo transfer day (Day 3 vs. Day 5) as treatment, implantation rate, clinical pregnancy rate, and odds of live birth as outcomes.

Our main purpose of this study was therefore to identify the causal effect of embryo transfer day on multiple major clinical outcomes in IVF. Additionally, we examine whether the causal effect of embryo transfer day on these outcomes can be modified by age and anti-Mullerian hormone (AMH) levels.


Patient selection

A retrospective cohort analysis was conducted using a dataset collected from Lee Women’s Hospital in Taichung, Taiwan. This study was approved by the Institutional Review Board on Biomedical Science Research at Academia Sinica. The dataset we obtained consisted of 9,090 fresh cycles from 1998 to 2014. This is a protocol based upon the first fresh transfer after oocyte retrieval and cryopreservation is just considered between the covariates in order to describe better patients characteristics. The information recorded in this database for each cycle included patients’ age at cycle start, the date of oocyte retrieval, cycle status, cycle outcome, AMH, follicle stimulating hormone (FSH), estradiol (E2), along with other 64 variables. In order to explore the causal effects of embryo transfer day on odds of live birth and other associated outcomes, only patients with transfer day of Day 3 and Day 5 and without missing or incomplete outcome information were included. To minimize the potential bias caused by numerical outliers, cycles with outliers in certain variables, such as age, weight, height, total number of oocytes retrieved, total sperm count, sperm motility, normal sperm morphology, AMH, FSH, period day of retrieval, number of embryos transferred and peak E2 were excluded. Cycles with PGD/PGS (preimplantation genetic diagnosis/preimplantation genetic screening) cases were also excluded. Therefore, the final dataset contained 4,127 autologous IVF cycles, with known embryo transfer day, and confirmed information about live births, number of gestational sacs, and number of fetal heartbeats at 4 weeks of gestation. Details about the inclusion and exclusion criteria of the study were documented as a flowchart in Fig. 1. Controlled ovarian stimulation was achieved by a downregulation protocol in which patients were administered leuprolide acetate (Lupron, Takeda Chemical Industries, Ltd., Osaka, Japan) during the midluteal phase. Subsequently, patients received recombinant follicle stimulating hormone (Gonal-F; Serono, Bari, Italy) on day 3 for ovarian stimulation. The retrieved oocytes were cultured in Quinn’s Advantage Fertilization Medium (Sage Bio- Pharma, Inc., Trumbull, CT, USA) with 15% serum protein substitute (SPS, Sage BioPharma, Inc) at 37C with the gas mixture of 5% oxygen, 5% carbon dioxide and 90% of nitrogen. Following conventional insemination or intracytoplasmic sperm injection (ICSI), a fertilization medium (SAGE Biopharma, USA) with 15% serum protein substitute (SPS; SAGE Biopharma, USA) was used to further culture the embryos. 70 to 72 hours after insemination or ICSI, all cleaved embryos were group cultured in microdrops of a blastocyst medium (Sage BioPharma, Inc.) with 15% SPS. Further details about the clinical aspects of the IVF protocols or embryology for the study were documented in a previous paper [16].

Fig. 1
figure 1

CONSORT diagram for data processing

Statistical methods

Comparison of baseline characteristics and handling missing data

To observe the discrepancies of marginal characteristics of the clinical characteristics between patients undergoing Day 3 embryo transfer and Day 5 embryo transfer, a series of statistical tests were performed. We applied Welch’s t-tests for continuous variables, chi-square tests for binary variables, and Poisson regressions for count variables. The comparisons of baseline characteristics were conducted both before and after matching to evaluate the performance of our constructed propensity score. Because data were largely missing in peak E2, we utilized Generalized Association Plots (GAP) [17] to retrieve the important variables highly associated with peak E2. A single imputation for peak E2 was then performed using the number of embryos frozen, number of oocytes retrieved, number of follicles, and embryo transfer day. In order to avoid the confounding effect of prescribed dosage given by multiple physicians, we converted Gonal-F (follitropin alfa), Menotropin and recombinant luteinizing hormone (rLH) into dichotomous values, and combined Menotropin and rLH into a single variable (whether or not the patient was injected with either Menotropin or rLH throughout the cycle).

Propensity score matching

To estimate the causal effect of embryos transfer day on different clinical outcomes, we implemented the method of propensity score matching to mimic a randomized controlled trial. The propensity score of receiving Day 5 transfer was constructed by a logistic regression where the probability of receiving Day 5 embryo transfer was modeled through a logit link function. The formula and deviation of the propensity score are detailed in Supple-ment. Embryo transfer day was regressed on a series of observed baseline covariates to derive a propensity score for each IVF cycle. The baseline covariates selected in the propensity score model were based on the literature as well as their clinical implications [18]. The variables (Xji,j=1,...,J) contributed to the propensity score included: age (linear and quadratic terms), body mass index (BMI), AMH, tubal factor (divided into six categories: undecided, removed, patent, sticking, ligation and obstructed for both left and right side of Fallopian tubes), sperm motility, sperm count, normal sperm morphology, azoospermia, polycystic ovarian syndrome (PCOS), Gonal-F, menotropin, rLH, peak estradiol, day of trigger, number of follicles, number of oocytes retrieved, number of embryos transferred, number of embryos frozen, number of mature oocytes and ICSI vs. IVF. Furthermore, to adjust for the secular trend regarding changes or modifications of the protocols made by clinical professionals, the year of cycles performed was also included as an additional matching factor, i.e., we only paired cycles that the embryos were implanted in the same years. Matched sets of Day 3 and Day 5 transfer subjects were thus formed: for each cycle whose embryo transfer was Day 5, we selected another cycle with the nearest propensity score from those with an embryo transfer Day 3. The method of greedy matching was used to produce balanced matched samples without replacement [19]. Those without matched cycles from the other transfer protocol were excluded from the subsequent analyses.

Post-matching analyses

Using the matched dataset, the effect of embryo transfer day on odds of live birth was then estimated through logistic regression, adjusting for propensity score and the remaining covariates that appeared to be statistically different between Day 3 and Day 5 transfers after propensity score matching. On the other hand, the effect of embryo transfer day on the number of gestational sacs and the number of fetal heartbeats was then estimated by Poisson regression where the number of embryos transferred was treated as an offset. Similarly, the Poisson regression adjusted for propensity score and the remaining significant covariates to account for residual confounding.

Stratification by age and AMH

Additional analyses were conducted stratifying cycles by age and AMH. We stratified cycles based on the age categories recommended by SART (age <35, 35 ≤age<37, 37 ≤age≤38 years) and the quantiles of AMH (AMH =0, 0 ≤AMH≤1.39, 1.39 <AMH≤3.95, AMH >3.95 ng/ml). The method of propensity score matching, logistic regressions and Poisson regressions were applied repeatedly to each subgroup to identify the causal effects of embryo transfer day on associated IVF outcomes.


Crude analyses

The clinical characteristics, live birth rate and other IVF associated outcomes among patients with Day 3 and Day 5 transfers before “Propensity score matching” section, were shown in Table 1. The crude live birth rate was slightly higher among those with Day 5 transfer (34.2% vs. 32.1%, p = 0.23); Day 5 transfers also showed advantage over Day 3 transfers on the average number of sacs (0.84 vs. 0.70, p <0.001) and number of fetal heartbeats (0.75 vs. 0.65, p = 0.001). The average age of patients in Day 3 was higher than that in Day 5 (32.89 vs. 31.61, p<0.001). For patients who had their embryos transferred on Day 3, the average AMH was much lower than that in group Day 5 (1.87 vs. 4.35, p<0.001), and the number of oocytes retrieved was fewer in those with Day 3 transfer (11.17 vs. 20.65, p<0.001). Other significantly different characteristics include BMI, azoospermia, PCOS, ICSI, Gonal-F, menotropin/rLH, peak E2, right duct, number of follicles, number of matured oocytes, and number of embryos frozen. As shown in Table 1, 16 out of 23 clinical characteristics were statistically significantly different between women receiving Day 3 transfer and those with Day 5, suggesting the two groups of patients had markedly distinct underlying features.

Table 1 Demographic characteristics before and after matching

Post matching analyses

After we matched our data using propensity score, only 577 patients were matched in each group, indicating that many patients had strong clinical indications for either Day 3 or Day 5 transfer and thus may not be matched with those receiving the other transfer protocol. Only one variable, the number of embryos frozen, out of 23 variables was still significantly different across the two groups after matching (Table 1). The distributions of the propensity scores across two groups before and after matching are shown in Fig. 2.

Fig. 2
figure 2

Propensity score distribution: a. Day 3 before matching; b. Day 5 before matching; c. Day 3 after matching; d. Day 5 after matching

A logistic regression was used to access the causal effect of embryo transfer day on odds of live birth, adjusting for the number of embryos frozen (Table 2). The analyses showed that embryo transfer day was not associated with odds of live birth (odds ratio 0.87 [95% confidence interval: 0.68, 1.11], p = 0.27). Our Poisson regression analyses also showed that embryo transfer day was not associated with clinical pregnancy rate (number of fetal heartbeats per embryo transferred; p = 0.11). Nevertheless, our analysis revealed a risk ratio (RR) of 1.14 in favor of Day 5 transfer in implantation rate (number of gestational sacs per embryos transferred; RR 1.14 [95% confidence interval: 1.0004, 1.30], p=0.04).

Table 2 Regression model estimating associated outcomes after matching

Stratified analyses

We further investigated the association between day of embryo transfer and IVF outcomes among subgroups of age and AMH, again using propensity score matching. Unlike our main analyses in Table 3, where blastocyst stage transfer is significantly associated with implantation rate, the “Stratified analyses” section by age (Table 3) and AMH (Table 4) showed that there was no evidence for association of embryo transfer day and IVF outcomes across all subgroups.

Table 3 Regression model estimating associated outcomes after matching stratified by age
Table 4 Regression model estimating associated outcomes after matching stratified by AMH


The principal finding of this study was that, comparing the women who shared similar baseline characteristics, there was evidence of a small but significant difference in implantation rate favoring Day 5 embryo transfer and no difference in the probability of live birth and clinical pregnancy.

The outcome of implantation, which is defined as the number of gestational sacs per embryo transferred, is often considered as the earliest primary endpoint in the process of IVF, and may affect the success of delivery. Our result showed that the blastocyst stage group had a higher rate of implantation, a relative risk of 1.14, than the group of cleavage stage transfer. This finding is in support of the previous studies that reported Day 5 transfers have purported advantages based on its morphological features [20, 21]. As the consequence of self selection, only the most viable embryos are expected to develop into the 64-cell blastocysts in the in vitro environment, eliminating those who had chromosomal abnormalities at an early age [22, 23]. Meanwhile, studies have also shown that the predictive morphological criteria of the embryos expected to be transferred on Day 3 are limited [2426]. In addition, it is recognized that premature exposure of early stage embryos to the uterine environment may induce homeostatic stress, reducing the potential of successful implantation [27]. This finding is also consistent with numerous trials reporting higher implantation rates for blastocyst transfer [2830].

Although blastocyst transfer is an effective procedure of increasing implantation rate, our result revealed that the transfer of embryos on Day 5 rather than on Day 3 did not change the overall probabilities of pregnancy and live birth. One of the potential cause for this phenomenon may be the existence of the post-treatment confounding similar to that of the lack of blinding effect in an RCT [31]. Blinding refers to maintaining unawareness of the assigned interventions among trial participants, healthcare providers, and assessors. However, even the goal of our method is to mimic an RCT in the context of an observational study, none of the people mentioned above in our study were blinded, which may leave room for potential bias due to residual confounding. For example, the clinicians may reflect their attitudes, paying extra attention or making differential decisions on the management, based on the allocation. Likewise, the patients may be more likely to comply with the clinical care after knowing they are about to or already had been assigned to Day 3 embryo transfer. Our propensity score was constructed based on the variables before implantation. However, numerous post-implantation factors may still affect the live birth rate. Even though blinding may reduced bias, the length of culture and the different procedures required based on the day of embryo transfer make it impossible to blind which group the patients were in for the technicians and clinicians.

Recent paper has revealed that serum AMH level is positively correlated with ovarian responsiveness, embryo developmental competence, and cumulative live birth rate [32]. Therefore, we further explored the effect of embryo transfer day on live birth and other associated outcomes by categorizing patients’ age and AMH. Interestingly, we were able to identify a trend that is consistent with our main finding: an increase in p values of embryo transfer day on IVF outcomes as the gestational age increases from 14 days (implantation rate) to week 4 (clinical pregnancy). For instance, the group where patients’ age were under 35, the p values of embryo transfer day increased from 0.23 to 0.48 as the outcome progressed from implantation rate to clinical pregnancy rate; as well as the group of patients aged between 37 and 38, the p values increased from 0.19 to 0.33. Although this trend was preserved in both our main result and stratified results, there was no evidence of significant difference in the later.

A recently published systematic review has performed meta-analysis on 27 RCTs to determine whether blastocyst transfers improve live birth and other associated IVF outcomes compared with cleavage stage transfers [33]. The study has revealed very similar results on clinical pregnancy rates with ours, both did not find difference between Day 3 transfer and Day 5 transfer. Yet, the review did report an evidence of a significant difference in live birth rate between the two groups, favoring Day 5 transfer. It is important to note that, although there is discrepancy between the two results on live birth, the intervention between the two studies are not the same. The previous study compared Day 2 and Day 3 transfers with Day 5 and Day 6 transfers, whereas our study only compared Day 3 transfers with Day 5 transfers.

Our results showing no evidence of increasing odds of live birth in Day 5 embryo transfer is consistent with a recent cohort conducted in the UK [34]. Although our studies share similar study design, the patient population between the two studies are different. A growing body of existing literature have indicated the racial and ethnic disparities that may appear in assisted reproductive technology outcomes. Current evidence have demonstrated the predisposition of Asian women having worse IVF outcomes, which may be resulting from fundamental biological or genetic differences [35, 36]. For example, studies have shown the distribution of FSH receptor allelic variants varies among Asians and Caucasians [37]. Alternatively, behavioral and environmental differences may have also caused the decreased pregnancy rate in Asian women [38, 39]. To our knowledge, our study is the biggest cohort in Asian population in assessing the differences between blastocyst stage transfer and cleavage stage transfer.

The goal of propensity score matching is to mimic an RCT from an observational study. The validity of this approach depends on a few assumptions. First, the propensity of receiving Day 5 transfer was properly modeled. The tendency of physician’s preference in prescribing a certain transfer protocol is unknown. On the other hand, to reduce the potential multi-dimensional confounding effect into a one-dimensional propensity score (PS) itself is a strong assumption. The assumption states that the unknown tendency is fully captured by the logistic model (model (1) in Supplement). It requires that all relevant variables are included and the feature of the variables (linear, quadratic, categorical, cross-product between variables, etc.) is correctly specified. However, limitations may arise in practice when we excluded certain highly correlated variables, such as endometriosis and ovarian cancer, in order to adjust the confounding effects and produce a homogeneous dataset. This may threaten the statistical generalizability of this study. Second, the tolerance of imperfect matching did not result in residual confounding. Numerically, we may not find a perfect match with an identical PS and systematic imperfect matches may lead to confounding. We address this issue by additionally adjust for propensity score in the “Post-matching analyses” section.

As an RCT require a carefully selected participants to ensure internal validity, its generalizability to population that does not satisfy its inclusion criteria is of concern. A similar issue applies in the PS matching approach. As shown in Fig. 2., the propensity scores between Day 3 and 5 transfers were quite different before matching, indicating that there are patients that were very likely to receive Day 3 transfer and unlikely to receive Day 5 transfer and vice versa. After matching, we only preserve those whose tendency of receiving either transfer is not that high. This phenomenon reflects the clinical equipoise, a common basis for clinical trials. Therefore, focusing on the patients who are equivocal for their tendency of receiving either transfer is a consequence of mimicking an RCT. However, what inherits from the equipoise principle is its difficulty in generalize the result to the patients not included in our matched data. For example, the null effect of Day 3 and Day 5 transfers may not be applicable to patients who have strong clinical indications of receiving Day 3 or Day 5 transfer. We stress that our matched analyses results need to be cautiously interpreted.


In “Conclusion” section, this study provides evidence that there is difference in implantation rate favoring Day 5 transfer compared to Day 3 transfer. However, considering later IVF outcomes, such as clinical pregnancy and live birth delivery, we suggest that advising women who share equivocal baseline characteristic at consultation to initiate an IVF cycle leading to a Day 5 transfer does not appear to increase the ongoing clinical pregnancy rate and odds of live birth compared with initiating an IVF cycle leading to Day 3 transfer.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.


Day 3 embryo transfer:

Cleavage stage embryo transferred on day 3

Day 5 embryo transfer:

Blastocyst stage embryo transferred on day 5


Randomized controlled trial


Anti-Mullerian hormone


Body mass index


In vitro fertilization


Follicle stimulating hormone




Preimplantation genetic diagnosis/preimplantation genetic screening


Generalized Association Plots


Follitropin alfa


Recombinant luteinizing hormone


Polycystic ovarian syndrome


Propensity score


Confidence interval


Odds ratio


  1. Jones G, Trounson A. Blastocyst stage transfer: pitfalls and benefits: The benefits of extended. Hum Reprod (Oxford. Print). 1999; 14(6):1405–8.

    Article  CAS  Google Scholar 

  2. Guerif F, Lemseffer M, Bidault R, Gasnier O, Saussereau M, Cadoret V, Jamet C, Royère D. Single day 2 embryo versus blastocyst-stage transfer: a prospective study integrating fresh and frozen embryo transfers. Hum Reprod. 2009; 24(5):1051–8.

    Article  CAS  PubMed  Google Scholar 

  3. Schoolcraft W, Gardner D. Blastocyst versus day 2 or 3 transfer. In: Seminars in Reproductive Medicine: 2001. p. 259–68. Copyright 2001 by Thieme Medical Publishers, Inc., 333 Seventh Avenue, New ….

  4. Gardner D, Vella P, Lane M, Wagley L, Schlenker T, Schoolcraft W. Culture and transfer of human blastocysts increases implantation rates and reduces the need for multiple embryo transfers. Fertil Steril. 1998; 69(1):84–8.

    Article  CAS  PubMed  Google Scholar 

  5. Papanikolaou E, Kolibianakis E, Tournaye H, Venetis C, Fatemi H, Tarlatzis B, Devroey P. Live birth rates after transfer of equal number of blastocysts or cleavage-stage embryos in ivf. a systematic review and meta-analysis. Hum Reprod. 2008; 23(1):91–9.

    Article  PubMed  Google Scholar 

  6. Gardner D, Schoolcraft W, Wagley L, Schlenker T, Stevens J, Hesla J. A prospective randomized trial of blastocyst culture and transfer in in-vitro fertilization. Hum Reprod (Oxford, England). 1998; 13(12):3434–40.

    Article  CAS  Google Scholar 

  7. Blake D, Proctor M, Johnson N, Olive D. Cleavage stage versus blastocyst stage embryo transfer in assisted conception. Cochrane Database Syst Rev. 2002; 2.

  8. Tsirigotis M. Blastocyst stage transfer: pitfalls and benefits. too soon to abandon current practice?. Hum Reprod (Oxford, England). 1998; 13(12):3285–9.

    Article  CAS  Google Scholar 

  9. Braga D, Setti A, Rita de Cássia SF, Machado R, Iaconelli A, Borges E. Patient selection criteria for blastocyst transfers in extended embryo culture programs. J Assist Reprod Genet. 2012; 29(12):1357–62.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Milki A, Hinckley M, Fisch J, Dasig D, Behr B. Comparison of blastocyst transfer with day 3 embryo transfer in similar patient populations. Fertil Steril. 2000; 73(1):126–9.

    Article  CAS  PubMed  Google Scholar 

  11. Rubin D. Estimating causal effects from large data sets using propensity scores. Annals of internal medicine. 1997; 8_Part_2:757–63.

    Article  Google Scholar 

  12. Rosenbaum P, Rubin D. The central role of the propensity score in observational studies for causal effects. Biometrika. 1983; 70(1):41–55.

    Article  Google Scholar 

  13. Rosenbaum P. Model-based direct adjustment. J Am Stat Assoc. 1987; 82(398):387–94.

    Article  Google Scholar 

  14. Austin P, Mamdani M. A comparison of propensity score methods: a case-study estimating the effectiveness of post-ami statin use. Stat Med. 2006; 25(12):2084–106.

    Article  PubMed  Google Scholar 

  15. Rubin D. Estimating causal effects of treatments in randomized and nonrandomized studies. J Educ Psychol. 1974; 66(5):688.

    Article  Google Scholar 

  16. Chen H-H, Huang C-C, Cheng E-H, Lee T-H, Chien L-F, Lee M-S. Optimal timing of blastocyst vitrification after trophectoderm biopsy for preimplantation genetic screening. PloS one. 2017; 12(10):0185747.

    Google Scholar 

  17. Chen C. Generalized association plots: information visualization via iteratively generated correlation matrices. Stat Sin. 2002; 12(1):7–29. Retrieved August 8, 2021.

    Google Scholar 

  18. Goldman R, Kaser D, Missmer S, Srouji S, Farland L, Racowsky C. Building a model to increase live birth rate through patient-specific optimization of embryo transfer day. J Assist Reprod Genet. 2016; 33(11):1525–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Austin P. A comparison of 12 algorithms for matching on the propensity score. Stat Med. 2014; 33(6):1057–69.

    Article  PubMed  Google Scholar 

  20. Huisman G, Alberda A, Leerentveld R, Verhoeff A, Zeilmaker G. A comparison of in vitro fertilization results after embryo transfer after 2, 3, and 4 days of embryo culture. Fertil Steril. 1994; 61(5):970.

    Article  CAS  PubMed  Google Scholar 

  21. Huisman G, Fauser B, Eijkemans M, Pieters M. Implantation rates after in vitro fertilization and transfer of a maximum of two embryos that have undergone three to five days of culture. Fertil Steril. 2000; 73(1):117–22.

    Article  CAS  PubMed  Google Scholar 

  22. Scott L, Alvero R, Leondires M, Miller B. The morphology of human pronuclear embryos is positively related to blastocyst development and implantation. Hum Reprod. 2000; 15(11):2394–403.

    Article  CAS  PubMed  Google Scholar 

  23. Magli M, Jones G, Gras L, Gianaroli L, Korman I, Trounson A. Chromosome mosaicism in day 3 aneuploid embryos that develop to morphologically normal blastocysts in vitro. Hum Reprod. 2000; 15(8):1781–6.

    Article  CAS  PubMed  Google Scholar 

  24. Rijnders P, Jansen C. The predictive value of day 3 embryo morphology regarding blastocyst formation, pregnancy and implantation rate after day 5 transfer following in-vitro fertilization or intracytoplasmic sperm injection. Hum Reprod (Oxford, England). 1998; 13(10):2869–73.

    Article  CAS  Google Scholar 

  25. Graham J, Han T, Porter R, Levy M, Stillman R, Tucker M. Day 3 morphology is a poor predictor of blastocyst quality in extended culture. Fertil Steril. 2000; 74(3):495–7.

    Article  CAS  PubMed  Google Scholar 

  26. Neuber E, Mahutte N, Arici A, Sakkas D. Sequential embryo assessment outperforms investigator-driven morphological assessment at selecting a good quality blastocyst. Fertil Steril. 2006; 85(3):794–6.

    Article  PubMed  Google Scholar 

  27. Gardner D, Lane M, Calderon I, Leeton J. Environment of the preimplantation human embryo in vivo: metabolite analysis of oviduct and uterine fluids and metabolism of cumulus cells. Fertil Steril. 1996; 65(2):349–53.

    Article  CAS  PubMed  Google Scholar 

  28. Papanikolaou E, Camus M, Kolibianakis E, Van Landuyt L, Van Steirteghem A, Devroey P. In vitro fertilization with single blastocyst-stage versus single cleavage-stage embryos. N Engl J Med. 2006; 354(11):1139–46.

    Article  CAS  PubMed  Google Scholar 

  29. Frattarelli J, Leondires M, McKeeby J, Miller B, Segars J. Blastocyst transfer decreases multiple pregnancy rates in in vitro fertilization cycles: a randomized controlled trial. Fertil Steril. 2003; 79(1):228–30.

    Article  PubMed  Google Scholar 

  30. Papanikolaou E, D’haeseleer E, Verheyen G, Van de Velde H, Camus M, Van Steirteghem A, Devroey P, Tournaye H. Live birth rate is significantly higher after blastocyst transfer than after cleavage-stage embryo transfer when at least four embryos are available on day 3 of embryo culture. a randomized prospective study. Hum Reprod. 2005; 20(11):3198–203.

    Article  PubMed  Google Scholar 

  31. Schulz K, Grimes D. Blinding in randomised trials: hiding who got what. Lancet. 2002; 359(9307):696–700.

    Article  PubMed  Google Scholar 

  32. Ezoe K, Ni X, Kobayashi T, Kato K. Anti-müllerian hormone is correlated with cumulative live birth in minimal ovarian stimulation with clomiphene citrate: a retrospective cohort study. BMC Pregnancy and Childbirth. 2020; 20(1):1–9.

    Article  Google Scholar 

  33. Glujovsky D, Farquhar C, Retamar A, Sedo C, Blake D. Cleavage stage versus blastocyst stage embryo transfer in assisted reproductive technology. Cochrane Database Syst Rev. 2016; 6.

  34. Cameron N, Bhattacharya S, McLernon D. Cumulative live birth rates following blastocyst-versus cleavage-stage embryo transfer in the first complete cycle of ivf: a population-based retrospective cohort study. Hum Reprod. 2020; 35(10):2365–74.

    Article  CAS  PubMed  Google Scholar 

  35. Purcell K, Schembri M, Frazier L, Rall M, Shen S, Croughan M, Grainger D, Fujimoto V. Asian ethnicity is associated with reduced pregnancy outcomes after assisted reproductive technology. Fertil Steril. 2007; 87(2):297–302.

    Article  PubMed  Google Scholar 

  36. Langen E, Shahine L, Lamb J, Lathi R, Milki A, Fujimoto V, Westphal L. Asian ethnicity and poor outcomes after in vitro fertilization blastocyst transfer. Obstet Gynecol. 2010; 115(3):591–6.

    Article  PubMed  Google Scholar 

  37. Simoni M, Nieschlag E, Gromoll J. Isoforms and single nucleotide polymorphisms of the fsh receptor gene: implications for human reproduction. Hum Reprod Update. 2002; 8(5):413–21.

    Article  CAS  PubMed  Google Scholar 

  38. Sharara F, Seifer D, Flaws J. Environmental toxicants and female reproduction. Fertil Steril. 1998; 70(4):613–22.

    Article  CAS  PubMed  Google Scholar 

  39. Choy C, Lam C, Cheung L, Briton-Jones C, Cheung L, Haines C. Infertility, blood mercury concentrations and dietary seafood consumption: a case–control study. BJOG: Int J Obstet Gynaecol. 2002; 109(10):1121–5.

    CAS  Google Scholar 

Download references


Not applicable.


The study was supported by Academia Sinica (06T-1070830-1Q). The funding body played no role in the study design, data collection, data analysis or interpretation of findings.

Author information

Authors and Affiliations



HCH conducted the statistical analysis and prepared the manuscript under the supervision of YTH1. EYL, YTH2 & JYS assisted with the statistical analysis. CIL,WLZ, CHC, CCH, PYL, MSL and ML were involved in hypothesis generation, subjects recruitment, data management, result interpretation. All authors have read and approved the final version of the manuscript. (YTH1: Yen-Tsung Huang; YTH2: Yi-Ting Huang)

Corresponding author

Correspondence to Yen-Tsung Huang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board on Biomedical Science Research at Academia Sinica (Case ID: AS-IRB-BM-18078). Since this is a retrospective study, it has been granted an exemption from informed consent by IRB committee. All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

Deviation and formula for propensity score.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hsieh, HC., Lee, CI., Lai, EY. et al. Estimating the causal effect of embryo transfer day on clinical in vitro fertilization outcomes using propensity score matching. BMC Pregnancy Childbirth 21, 554 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: