Skip to main content
Download PDF

Tracking of menstrual cycles and prediction of the fertile window via measurements of basal body temperature and heart rate as well as machine-learning algorithms

Reproductive Biology and Endocrinology volume 20, Article number: 118 (2022) Cite this article

Abstract

Background

Fertility awareness and menses prediction are important for improving fecundability and health management. Previous studies have used physiological parameters, such as basal body temperature (BBT) and heart rate (HR), to predict the fertile window and menses. However, their accuracy is far from satisfactory. Additionally, few researchers have examined irregular menstruators. Thus, we aimed to develop fertile window and menstruation prediction algorithms for both regular and irregular menstruators.

Methods

This was a prospective observational cohort study conducted at the International Peace Maternity and Child Health Hospital in Shanghai, China. Participants were recruited from August 2020 to November 2020 and followed up for at least four menstrual cycles. Participants used an ear thermometer to assess BBT and wore the Huawei Band 5 to record HR. Ovarian ultrasound and serum hormone levels were used to determine the ovulation day. Menstruation was self-reported by women. We used linear mixed models to assess changes in physiological parameters and developed probability function estimation models to predict the fertile window and menses with machine learning.

Results

We included data from 305 and 77 qualified cycles with confirmed ovulations from 89 regular menstruators and 25 irregular menstruators, respectively. For regular menstruators, BBT and HR were significantly higher during fertile phase than follicular phase and peaked in the luteal phase (all P < 0.001). The physiological parameters of irregular menstruators followed a similar trend. Based on BBT and HR, we developed algorithms that predicted the fertile window with an accuracy of 87.46%, sensitivity of 69.30%, specificity of 92.00%, and AUC of 0.8993 and menses with an accuracy of 89.60%, sensitivity of 70.70%, and specificity of 94.30%, and AUC of 0.7849 among regular menstruators. For irregular menstruators, the accuracy, sensitivity, specificity and AUC were 72.51%, 21.00%, 82.90%, and 0.5808 respectively, for fertile window prediction and 75.90%, 36.30%, 84.40%, and 0.6759 for menses prediction.

Conclusions

By combining BBT and HR recorded by the Huawei Band 5, our algorithms achieved relatively ideal performance for predicting the fertile window and menses among regular menstruators. For irregular menstruators, the algorithms showed potential feasibility but still need further investigation.

Trial registration

ChiCTR2000036556. Registered 24 August 2020.

Introduction

Accurate menstrual cycle tracking and identification of the fertile window are in high demand worldwide for determining the optimal time of intercourse among couples who desire to conceive. Wilcox et al. were among the first to identify the fertile window, the most fecund period of the menstrual cycle, which is comprised of the five days preceding ovulation and the day of ovulation [1]. Many biomarkers, including changes in basal body temperature (BBT), vaginal discharge, the position of cervix, and luteinizing hormone (LH) levels (which can be assessed in home urine tests), are used to detect ovulation and the fertile window to aid conception; among these, monitoring of BBT is one of the most convenient and least invasive [2]. The nadir (the lowest point) of BBT occurs approximately one day before ovulation [2]. In response to elevated levels of progesterone, BBT rises significantly approximately 2 days after the LH peak [3]. According to a cohort study, 21% of women have used BBT measurements to track their fecundity [4]. However, studies have shown that BBT measurements are not able to precisely predict the fertile window [5, 6]. In addition, traditional BBT records can only retrospectively confirm the fertile window by providing evidence of ovulation 2 or 3 days after its occurrence [3]. Even though there may be a BBT nadir before ovulation, this change does not feasibly predict the fertile window since not all nadirs coincide with ovulation [6].

Over the past few years, the use of websites and cellular phone applications (apps) to track menstrual cycles and predict the ovulation day or fertile window has become increasingly popular [7]. These apps are based on records of BBT, cervical fluid, cervix position and urine LH levels [8]. However, despite their convenience, the accuracy of their fertile window predictions is generally far from satisfactory. Current applications based on BBT have low prediction accuracy for ovulation since the algorithm is dependent on BBT data or previous user cycles [9]. Thus, more advanced app-based technology is urgently needed to aid couples who desire to conceive.

Advances in portable sensors and wearable technology provide continuous dynamic health information over the course of a day. In addition, machine learning based on health data has been increasingly applied to clinical medicine, including the field of fertility awareness. Pairing machine learning with health data from wearable devices has the potential to improve the practice of clinical medicine. Some studies have already used machine learning to predict menstruation, ovulation and the fertile window based on body temperature data recorded by wearables [10]. Since previous studies have shown that heart rate (HR) varies during different phases of the menstrual cycle, with a higher rate during ovulation [11], some studies have employed wearables that collect data on body temperature, HR, heart rate variability (HRV), respiratory rate, etc., to develop algorithms for the prediction of ovulation and the fertile window [12]. However, most of these studies failed to specifically focus on women with irregular menstrual cycles [13, 14]. Little is known about whether physiological data collected by wearables can be used to track menstrual cycles and predict the fertile window among women with irregular menstrual cycles. Furthermore, few studies have focused on the menstruation period of women with irregular cycles.

In this study, we aimed to evaluate alterations in BBT (measured with an ear thermometer) and HR (recorded by a wearable device, the Huawei Band 5) during different phases in the menstrual cycle. In addition, we further developed machine-learning algorithms that integrated BBT and HR data to predict the fertile window and menstruation days among women with regular menstrual cycles. Finally, we explored the feasibility of fertile window and menstruation day prediction among women with irregular menstrual cycles.

Methods

Study design and participants

This was a prospective observational cohort study conducted at the International Peace Maternity and Child Health Hospital (IPMCHH) in Shanghai, China. Participants were recruited from August 2020 to November 2020 and followed up between August 2020 and May 2021. Eligible participants were nonpregnant women aged 18–45 with natural menstrual cycles. For this study, we excluded women who had major diseases, had a history of pregnancy within the past 6 months, were breastfeeding, were taking medications that could interfere with the menstrual cycle (e.g., hormonal contraception or hormone replacement), traveled across time zones during the follow-up period or had a sleeping disorder. We divided women by menstrual cycle length, which was defined as the time in days from the first day of one menses to the first day of the next [15]. Since self-reported menstrual cycle length is proven reliable [15], women were queried about the usual length of the menstrual cycle without hormonal therapies in the past year [16,17,18]. Women with usual menstrual cycle lengths of 25–35 days were included in the regular group. Women with usual menstrual cycle lengths outside of this range were included in the irregular group [19, 20].

We recruited women via advertisements on WeChat APP of doctors at the IPMCHH. At enrollment, participants completed baseline questionnaires with items on age, weight, height, marital status, educational attainment, occupation, age at menarche, smoking status, alcohol consumption, and history of pregnancy and childbirth. The participants also received a wearable sensor (Huawei Band 5; Huawei Device Co, Ltd, Shenzhen, China), a smartphone (Huawei Mate 30; Huawei Device Co, Ltd, Shenzhen, China) and an ear thermometer (Braun IRT6520) to record essential physiological data. Participants were followed up at least four menstrual cycles and encouraged to retain in the cohort until they had completed 4 qualified menstrual cycles, which was defined as synced data for 80% of the cycle durations. Women were gifted the study materials (i.e., the sensor, smartphone, and thermometer) at the end of their follow-up.

Study protocol

Data collection

During the follow-up period, women were asked to wear the Huawei Band 5 every night while sleeping and to sync this data with their smartphones every morning. For data collection, the duration of continuous sleep had to exceed 4 h every night. The Huawei Band 5 measures HR and HRV continuously during sleep. It can also measure data related to sleep quantity and sleep quality. Since our research did not focus on sleep, these data were not included in the analysis. BBT was assessed daily with an ear thermometer in the morning upon waking, after participants had lain horizontally for 5 min. Participants were instructed to report menstruation on the smartphone by answering “yes” or “no” to the two questions (i.e., “Did your period start/end today?”) every day.

Determination of ovulation

Ovulation was determined by an ultrasound and detection of serum hormone levels [21]. In the first cycle, participants underwent transvaginal or abdominal ultrasound from cycle day 8–12 to the day that a follicle reached 17 mm [22]. Serum hormone levels, including LH, estradiol (E2), follicle-stimulating hormone (FSH), progesterone, testosterone, prolactin, and anti-Müllerian hormone (AMH), were measured on the day of large (i.e., at least 17 mm) follicle detection. According to hormone levels, another ultrasound was performed several days later to confirm ovulation (follicular enlargement followed by evidence of rupture). In cases where follicular enlargement was missed, the progesterone level was assayed. Based on the day of maximum follicular enlargement and levels of LH, E2, and progesterone as well as evidence of rupture, the day of ovulation was estimated.

During the next three menstrual cycles, the first day of follow-up was determined by a senior gynecologist according to the ovulation day in the first cycle. Basically, it would be one or two days before the previous ovulation day. All monitoring procedures were the same as those in the first cycle except that the blood samples were analyzed only for levels of LH, E2, FSH and progesterone. In the present study, data of cycles without determined ovulation were excluded.

According to the days of menstruation and ovulation, the menstrual cycle was divided into 4 phases. The duration of menstrual flow was defined as the menstrual phase reported by each woman [23]. The follicular phase (outside menses and the fertile window) was the time between the first day post-menses and 6 days before ovulation [24]. The fertile phase lasted from 5 days before ovulation to the ovulation day [1, 19, 25]. The luteal phase started post-ovulation and lasted until the day before menses [24].

Statistical analysis

Participant characteristics are presented with descriptive statistics. Continuous variables that were not normally distributed are displayed as medians and interquartile ranges. These variables were compared with the Wilcoxon test. Categorical variables are expressed as frequencies with proportions, with group differences assessed with the chi-square test and Fisher’s exact test.

Before analyzing the synced data from the wearable sensors, we excluded data collected during the first 30 min after sleep and the last 30 min before waking. HR was extracted from the photoplethysmographic signal. HRV indices, including the standard deviation of normal-to-normal intervals (SDNN) values and the low frequency (LF)/high frequency (HF) ratio, were obtained by analyzing R-R interval recordings. Missing values were replaced by the mean of values from two days before and after the day with missing data.

Data on qualified menstrual cycles with identified ovulation days were included in the analysis. Normality tests [26] indicated that physiological parameters BBT and HR were normally distributed, but SDNN and LF/HF ratio were skewed. Thus, SDNN and LF/HF ratio were natural log transformed (ln) to improve normality [27]. We did further analysis on the transformed data. Each participant contributed data from several menstrual cycles, which created a nested structure in our data. To estimate the effect of cycle phase on physiological parameters, we used a linear mixed model that estimated fixed and random effects [28] (R packages: lme4 and lmerTest). We specified cycle phase as a fixed effect and participants/menstrual cycle as a random intercept in the model. For multiple comparisons, the Bonferroni method was used to adjust the P value (R package: multcomp). We then set the group as a fixed effect and ran a series of models as mentioned above to investigate the group differences in physiological parameters. We tested all linear mixed models by examining their residuals.

All data analyses and visualizations were performed with R software (version 4.0.2, R Foundation for Statistical Computing, Vienna, Austria). A two-sided P value of less than 0.05 was considered statistically significant.

Construction of the algorithms

According to best practices, BBT and HR waveforms were smoothed at the single cycle level. From the changes and analysis of physiological signs, biphasic rule and features were extracted from both waveforms and the fertile day was most likely located on the switch point of the waveforms. SDNN and LF/HF ratio were preprocessed in a similar way. More details were provided in supplementary methods. Based on the analysis and data preparation, we developed a probability function estimation model to estimate the probability of whether a day fell in the fertile window. Specifically, to estimate whether Day 1 was in the fertile window, data from the last consecutive period before Day 1 were used for model input. The construction of the menstruation prediction model is similar to that of the fertile window prediction model.

The dataset of enrolled patients was partitioned randomly into training and testing groups. In detail, data from 225 qualified cycles from 68 participants in the regular group were used as the training dataset; 80 qualified cycles from 21 participants in the regular group and 77 cycles from 25 participants in the irregular group were included in the two testing datasets respectively. Based on a training dataset, a fivefold cross validation method was used to determine the model. Firstly, the training data was divided into 5 parts. In every iteration, 4 parts were aggregated as training subsets and the last part as a testing subset. Then the hyperparameters were optimized and adjusted, and the model was preliminarily evaluated. Finally, the model performance was evaluated with the testing dataset. We calculated the accuracy, sensitivity, specificity, and area under the receiver operating characteristic (ROC) curve (AUC) to measure the algorithm performance. The algorithms were developed and visualized using Python software (version 3.10.4, The Python Software Foundation, USA).

Results

A total of 153 eligible women with regular menstrual cycles (n = 103) or irregular menstrual cycles (n = 50) were initially recruited for this study (Fig. 1). However, 3 women in the regular group withdrew early because they were no longer willing to participate, and one participant was lost to follow-up after the first cycle. Therefore, 149 women eventually completed all procedures. Of these participants, 10 women in the regular group and 14 women in the irregular group were excluded due to > 20% missing data. In addition, 11 women in the irregular group were excluded from the final analysis due to > 20% missing data and inability to determine ovulation day. Thus, by the end of data collection, 89 women in the regular group and 25 women in the irregular group completed qualified menstrual cycles, synced data on which was 80% complete and which had determined ovulation days. Data from these 114 women were included in the final analysis.

Fig. 1
figure 1

Flow chart of the study population

Full size image

The characteristics of all participants are summarized in Table S1. Table 1 presents the characteristics of the women included in the final analysis. Women in the irregular group were younger than those in the regular group (P = 0.005). However, body mass index (BMI), marital status, educational attainment, occupation, duration of menstruation, age at menarche, smoking status, alcohol consumption, and history of pregnancy and childbirth were comparable between the two groups. The AMH levels significantly differed between the two groups (P = 0.035); the irregular group was more likely to have a low AMH (< 0.7 ng/ml) (12.0% vs. 4.5%) or high AMH (> 7.5 ng/ml) (28.0% vs. 12.4%) than the regular group. During the entire follow-up period, women in the regular group had cycles with a median length of 29 (interquartile range [IQR] = 28–32) days, while women in the irregular group had a median cycle length of 33 (IQR = 27–40) days.

Table 1 Characteristics of participants included in the analysis
Full size table

The physiological parameters changed regularly throughout the menstrual cycle (Fig. 2 and Figure S1). For regular menstruators (Table 2), BBT was significantly lower in the follicular phase than in the menstrual phase (P < 0.001). In the fertile phase, BBT was 0.04 ℃ higher than in the follicular phase (P < 0.001). BBT in the luteal phase increased 0.18 ℃ compared with that in the fertile phase (P < 0.001). Unlike BBT, HR, ln(SDNN) and ln(LF/HF ratio) in the follicular phase were similar to those in the menstrual phase. Consistent with BBT, HR in the fertile phase was significantly higher (0.61 bpm; P < 0.001) than that in the previous stage. Furthermore, HR peaked in the luteal phase, increasing by 1.84 bmp (P < 0.001) compared with the fertile phase. In contrast to HR, ln(SDNN) in the fertile phase was lower (-0.04; P < 0.001) than that in the follicular phase. Ln(SDNN) reached its lowest value in the luteal phase. Compared to that in the follicular phase, ln(LF/HF ratio) in the fertile phase increased 0.03 (P < 0.001); this increase was maintained in the luteal phase. Throughout the whole cycle, the fertile phase always was the "inflection point" for the physiological parameters.

Fig. 2
figure 2

Physiological parameters of the regular group in the different phases of the menstrual cycle. Changes in BBT (A), HR (B), SDNN (C) and the LF/HF ratio (D) during the menstrual cycle of regular menstruators are depicted. The horizontal line represents the medians, boxes represent the values between 25–75%, and lines represent the values between 5–95%. BBT: basal body temperature; HR: heart rate; LF/HF: low frequency /high frequency ratio; SDNN: standard deviation of normal-to-normal intervals

Full size image
Table 2 The relationship between menstrual cycle and physiological parameters in the regular group
Full size table

The trend of changes in physiological parameters in the irregular group throughout the menstrual cycle was similar to that in the regular group (Table S2). For menstrual cycles with an ovulation day, the physiological parameters of all phases were similar between the regular and irregular groups (Table 3).

Table 3 Comparison of physiological parameters between the regular and irregular groups
Full size table

We developed a series of predictive models (Table 4) for detecting the fertile window and menstruation using data from 114 women, including 305 qualified cycles with ovulation from 89 participants in the regular group and 77 qualified cycles with ovulation from 25 participants in the irregular group. The algorithm based on the BBT of the regular group achieved an accuracy of 86.65%, sensitivity of 68.30% and specificity of 91.30% in predicting the fertile window of women with regular menstrual cycles. Furthermore, the prediction performance improved to an accuracy of 87.46%, sensitivity of 69.30% and specificity of 92.00% when integrating HR collected by the wearable Huawei Band 5. The algorithm developed from data from regular menstruators was also applied to predict the fertile window of irregular menstruators and achieved an accuracy of 72.51%, sensitivity of 21.00%, and specificity of 82.90%. The algorithm trained with the BBT data from regular menstruators predicted menses with an accuracy of 87.80%, sensitivity of 66.10% and specificity of 93.10% in women with regular menstrual cycles. When we integrated HR data into the model, the prediction performance improved to an accuracy of 89.60%, sensitivity of 70.70% and specificity of 94.30%. The algorithm was also applied to predict menses of irregular menstruators and achieved an accuracy of 75.90%, sensitivity of 36.30%, and specificity of 84.40% for menstrual phase prediction. The ROC curves and the AUCs of different models are provided in Fig. 3. The fertile window prediction models developed among regular menstruators using BBT and BBT combined with HR had AUCs of 0.8986 and 0.8993, respectively. The latter one was applied to predict the fertile window among irregular menstruators and achieved an AUC of 0.5808. The menstruation prediction models developed among regular menstruators based on BBT and BBT combined with HR had AUCs of 0.7787 and 0.7849, respectively. The latter one was applied to predict menses among irregular menstruators and achieved an AUC of 0.6759.

Table 4 Accuracy, sensitivity and specificity of different models in fertile window and menstruation prediction
Full size table
Fig. 3
figure 3

Prediction performance of different models. ROC curves of fertile window (A) and menstruation (B) prediction models based on BBT and HR for regular and irregular groups. AUC: area under the receiver-operating characteristic curve; BBT: basal body temperature; HR: heart rate; ROC: receiver operating characteristic

Full size image

To further compare the performance of the current model in predicting fertile window with other algorithms, we tried to integrate SDNN and LF/HF ratio data into the model and use the data from the irregular group for training (Table S3). Overall, the models developed using four features had similar performance with the model based on BBT and HR. Among models based on BBT and HR, the model trained using BBT and HR data from the regular group still achieved the best accuracy, sensitivity and specificity among regular menstruators. It achieved the same sensitivity (35.8%), a slightly lower accuracy (78.78% vs 77.31%) and specificity (87.60% vs 85.80%) among irregular menstruators compared to the model based on data of the irregular group. In conclusion, the model which included only BBT and HR based on data from participants with regular menstrual cycles provided relatively robust and stable performance for fertile window prediction among both regular and irregular menstruators.

In order to improve the performance of fertile window prediction among irregular menstruators, we divided them into two subgroups based on the cycle length. The current model had better predictive accuracy and sensitivity in women with cycle length > 35 days and better predictive specificity in women with cycle length < 25 days (Table S4). Additionally, the current model was more suitable for women over 35 years old in predicting the fertile window among both regular menstruators and irregular menstruators (Table S5).

Discussion

Principal findings

To the best of our knowledge, this is the first study to apply wearable technology and machine learning methods to predict the fertile window and menstruation days for women with regular and irregular cycles. Here, we demonstrated that 1) the in-ear temperature and HR collected by the Huawei Band 5 followed a similar pattern in women with regular and irregular menstrual cycles; 2) the algorithm that used BBT data and HR data from the Huawei Band 5 outperformed the algorithm that used BBT data alone in predicting the fertile window and menstruation of women with regular menstrual cycles in terms of accuracy, sensitivity, specificity and AUC; 3) the fertile window and menstruation of women with irregular menstrual cycles can potentially be predicted using the algorithms but with lower accuracy, sensitivity, specificity and AUC.

According to the results, the tympanic membrane temperature during the luteal phase differed significantly from that of the menstrual phase, at approximately 0.2 ℃ higher; this difference was expected according to previous literature [3, 25]. Our finding of the pattern of changes in HR were consistent with some other studies that measured HR using electrocardiography. Matthew et al. examined HRV across the menstrual cycle in 13 women with normal menstruation, confirming an increase in HR and a decrease in ln(SDNN) after ovulation [11]. Similar to our findings, Bai et al. also found that the LF/HF ratio and HR increased from the follicular phase to the luteal phase [29]. These results demonstrate that data collected by the Huawei Band 5 were capable of reflecting phase-based changes in HR as accurately as an electrocardiogram. In other studies that used wearables to monitor HR, the Ava bracelet indicated an increase in nightly HR in the luteal phase, but the HRV ratio decreased in the luteal phase, which was different from our study [12]. Changes in HR and HRV may be explained by a decrease in parasympathetic activity after ovulation, in response to increased progesterone [11]. As both BBT and HR data reflected phase-based alterations during the menstrual cycle, it is reasonable that the incorporation of HR data in the algorithm improved the predictive power in our study.

In addition to the Huawei Band 5 used in our study, other wearables, such as the Ava bracelet, Oura ring, in-ear wearable thermometer and OvulaRing (a vaginal biosensor), have been used to record physiological parameters and thereby track menstrual cycles, with machine learning algorithms used to predict the ovulation day and fertile window [10, 12, 14, 30, 31]. Most of these studies exclusively recorded one parameter, body temperature; the Ava bracelet simultaneously monitored additional parameters, including HR, HRV, respiratory rate and skin perfusion. Therefore, it is reasonable that the Ava bracelet demonstrated better performance than other wearables in predicting the fertile window, with an accuracy of 90%, a specificity of 93%, and a sensitivity of 81% [12]. In our study, the accuracy, sensitivity, specificity and AUC of fertile window prediction and menstruation prediction also improved after including HR data collected with the Huawei Band 5. These findings indicate the importance of including various data collected by wearables to establish more reliable prediction models for the fertile window and menstruation. The higher sensitivity and specificity can better distinguish the actual days within the fertile window from other phases, thus helping couples optimize their chance of achieving pregnancy. Surprisingly, the performance of the algorithm did not improve a lot after inclusion data of SDNN and LF/HF ratio, with even lower accuracy, sensitivity and specificity for fertile window prediction among the regular group. This can be possibly explained by the greater dispersion of data of SDNN and LF/HF ratio than BBT and HR in our study (data not shown). Additionally, since the algorithm developed in this study will finally be applied into the software in the smartphone and wearable device, both functionality and simplicity should be taken into consideration. A robust performance based on relatively less features will increase the compatibility and stability of the software. What’s more, the definition of the fertile window in our study was in accordance with the gold standard [1, 32, 33]; therefore, our prediction results were more professional and reliable than those of some other studies [14, 31].

To the best of our knowledge, this is the first study to include both women with regular cycles and those with irregular cycles to predict the fertile window and menstruation with the aid of wearables and machine learning, and the application of algorithms to both groups supports the generalizability of our results. Menstrual trackers such as smartphone apps or wearable sensors which collect women’s health data have become increasingly popular, but most of these studies only included women with regular menstruation cycles [12, 13, 34]. In real-world scenarios, it is impossible for all women to have menstrual cycles that fall within a certain length. A recent study based on 32,595 women found that approximately 13% of participants had actual menstrual cycle lengths that did not fall within the range of 23–35 days. In addition, even in women with a typical 28-day cycle, these apps failed to precisely predict the fertile window [7]. Therefore, the development of more multifaceted and accurate apps to track menstrual cycles is needed for women with irregular menstrual cycles.

Our study showed that data derived from ovulatory women, regardless of their menstrual pattern, have the potential to predict the fertile window and menstruation of women with irregular menstruation cycles. Ovulatory women with a cycle length that either falls within or outside of 25–35 days exhibit similar trends in physiological parameters during the menstrual cycle. The physiological parameters in all phases were similar between regular menstruators and irregular menstruators. These findings establish a theoretical foundation for the development of a prediction model. Accordingly, our preliminary algorithm had a relatively high specificity, 82.90%, and a low sensitivity, 21.00%, for women with menstrual irregularity in fertile window prediction; the algorithm also had a relatively high specificity, 84.40%, and a low sensitivity, 36.30% in menstruation prediction. This indicates a high false-negative rate but a low false-positive rate, which is necessary to detect the actual fertile window and menstruation. We also tried to improve the performance through subgroup analysis based on the cycle length and age of participants, especially for the irregular group. Based on previous literatures, short or long cycles had different pathophysiological mechanisms. A short cycle length could be associated with luteal phase deficiency and a long cycle length could be associated with polycystic ovary syndrome [35, 36]. In addition, age has been reported to be related with physiological parameters during menstrual cycles [37, 38]. These could partly explain why some of the predictive score improved in certain subgroups, while the sample size was too small to draw a solid conclusion. Further studies with larger sample sizes and more variables are needed to optimize these results.

Another strength of our study is that we used a combination of transvaginal sonographic examinations and serum LH levels as our reference standard to determine ovulation. This provides the most reliable method of assessing and predicting ovulatory status. Previous studies estimated ovulation onset via a home urine LH test without resorting to confirmatory ultrasounds or serum hormone tests [12, 14]. Based on prior studies comparing urinary LH kits and endometrium biopsy, false-positive test results are frequent [39]. In addition, ovulation does not always occur one day after the surge in urinary hormone levels [40]. These factors could reduce the precision of ovulation day confirmation. The serial transvaginal follicular ultrasound used in this study serves as the gold standard for identifying the ovulation day [41], thus enhancing the reliability of our prediction algorithms.

Limitations

This study has some limitations regarding both the device and the algorithm. We only included two parameters, BBT and HR, in the current algorithm. If more variables are measured in the future, the performance of prediction models should be further improved. In addition, the tympanic membrane temperature had to be measured with an ear thermometer and recorded in the app daily, which may have decreased the compliance of participants and the accuracy of data. Our approach would be difficult to apply if a user failed to test or record her temperature in the morning. In the future, a wearable device capable of collecting real-time temperature measurements can improve the performance and generalizability of the machine learning model.

Conclusions

In conclusion, we demonstrated that BBT and indices related to HR changed in a specific pattern during the menstrual cycle for both regular and irregular menstruators. This pilot study also exhibited the efficacy of machine learning in predicting fertile windows and menstruation in regular menstruators using BBT and HR data collected by a wearable device. We further established the feasibility of fertile window and menstruation prediction among irregular menstruators. Since determining the days when sexual intercourse is more likely to result in successful conception is highly important in clinical practice, it is expected that continuous health data monitoring through a wearable device and a well-established algorithm will help women successfully conceive in the future.

Availability of data and materials

The data for this study is available from the corresponding author on reasonable request.

Abbreviations

AMH:

Anti-Müllerian hormone

Apps:

Applications

AUC:

Area under the receiver-operating characteristic curve

BBT:

Basal body temperature

BMI:

Body mass index

E2:

Estradiol

FSH:

Follicle-stimulating hormone

HF:

High frequency

HR:

Heart rate

HRV:

Heart rate variability

LF:

Low frequency

LH:

Luteinizing hormone

ROC:

Receiver operating characteristic

SDNN:

Standard deviation of normal-to-normal intervals

References

  1. Wilcox AJ, Weinberg CR, Baird DD. Timing of sexual intercourse in relation to ovulation. Effects on the probability of conception, survival of the pregnancy, and sex of the baby. N Engl J Med. 1995;333(23):1517–21. https://doi.org/10.1056/NEJM199512073332301 PMID: 7477165.

    Article  PubMed  CAS  Google Scholar 

  2. Su HW, Yi YC, Wei TY, Chang TC, Cheng CM. Detection of ovulation, a review of currently available methods. Bioeng Transl Med. 2017;2(3):238–46. https://doi.org/10.1002/btm2.10058 PMID: 29313033.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Moghissi KS. Prediction and detection of ovulation. Fertil Steril. 1980;34(2):89–98. https://doi.org/10.1016/s0015-0282(16)44888-0 PMID: 6773821.

    Article  PubMed  CAS  Google Scholar 

  4. Stanford JB, Willis SK, Hatch EE, Rothman KJ, Wise LA. Fecundability in relation to use of fertility awareness indicators in a North American preconception cohort study. Fertil Steril. 2019;112(5):892–9. https://doi.org/10.1016/j.fertnstert.2019.06.036 PMID: 31731946.

    Article  PubMed  Google Scholar 

  5. Bauman JE. Basal body temperature: unreliable method of ovulation detection. Fertil Steril. 1981;36(6):729–33. https://doi.org/10.1016/s0015-0282(16)45916-9 PMID: 7308516.

    Article  PubMed  CAS  Google Scholar 

  6. Barron ML, Fehring RJ. Basal body temperature assessment: is it useful to couples seeking pregnancy? MCN Am J Matern Child Nurs. 2005;30(5):290–6; quiz78. https://doi.org/10.1097/00005721-200509000-00004 PMID: 16132004.

    Article  PubMed  Google Scholar 

  7. Setton R, Tierney C, Tsai T. The Accuracy of Web Sites and Cellular Phone Applications in Predicting the Fertile Window. Obstet Gynecol. 2016;128(1):58–63. https://doi.org/10.1097/AOG.0000000000001341 PMID: 27275788.

    Article  PubMed  Google Scholar 

  8. Stanford JB, Willis SK, Hatch EE, Rothman KJ, Wise LA. Fecundability in relation to use of mobile computing apps to track the menstrual cycle. Hum Reprod. 2020;35(10):2245–52. https://doi.org/10.1093/humrep/deaa176 PMID: 32910202.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Stanford JB, White GL, Hatasaka H. Timing intercourse to achieve pregnancy: current evidence. Obstet Gynecol. 2002;100(6):1333–41. https://doi.org/10.1016/s0029-7844(02)02382-7 PMID: 12468181.

    Article  PubMed  Google Scholar 

  10. Luo L, She X, Cao J, Zhang Y, Li Y, Song PXK. Detection and Prediction of Ovulation From Body Temperature Measured by an In-Ear Wearable Thermometer. IEEE Trans Biomed Eng. 2020;67(2):512–22. https://doi.org/10.1109/TBME.2019.2916823 PMID: 31095472.

    Article  PubMed  Google Scholar 

  11. Tenan MS, Brothers RM, Tweedell AJ, Hackney AC, Griffin L. Changes in resting heart rate variability across the menstrual cycle. Psychophysiology. 2014;51(10):996–1004. https://doi.org/10.1111/psyp.12250 PMID: 24942292.

    Article  PubMed  Google Scholar 

  12. Goodale BM, Shilaih M, Falco L, Dammeier F, Hamvas G, Leeners B. Wearable Sensors Reveal Menses-Driven Changes in Physiology and Enable Prediction of the Fertile Window: Observational Study. J Med Internet Res. 2019;21(4):e13404. https://doi.org/10.2196/13404 PMID: 30998226.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Shilaih M, Goodale BM, Falco L, Kubler F, De Clerck V, Leeners B. Modern fertility awareness methods: wrist wearables capture the changes in temperature associated with the menstrual cycle. Biosci Rep. 2018;38(6). https://doi.org/10.1042/BSR20171279. PMID: 29175999.

  14. Maijala A, Kinnunen H, Koskimaki H, Jamsa T, Kangas M. Nocturnal finger skin temperature in menstrual cycle tracking: ambulatory pilot study using a wearable Oura ring. BMC Womens Health. 2019;19(1):150. https://doi.org/10.1186/s12905-019-0844-9 PMID: 31783840.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Jukic AM, Weinberg CR, Wilcox AJ, McConnaughey DR, Hornsby P, Baird DD. Accuracy of reporting of menstrual cycle length. Am J Epidemiol. 2008;167(1):25–33. https://doi.org/10.1093/aje/kwm265 PMID: 17928401.

    Article  PubMed  Google Scholar 

  16. Wei S, Schmidt MD, Dwyer T, Norman RJ, Venn AJ. Obesity and menstrual irregularity: associations with SHBG, testosterone, and insulin. Obesity (Silver Spring). 2009;17(5):1070–6. https://doi.org/10.1038/oby.2008.641 PMID: 19180069.

    Article  CAS  Google Scholar 

  17. Wang ET, Cirillo PM, Vittinghoff E, Bibbins-Domingo K, Cohn BA, Cedars MI. Menstrual irregularity and cardiovascular mortality. J Clin Endocrinol Metab. 2011;96(1):E114-8. https://doi.org/10.1210/jc.2010-1709 PMID: 20980429.

    Article  PubMed  CAS  Google Scholar 

  18. Lagana AS, Veronesi G, Ghezzi F, Ferrario MM, Cromi A, Bizzarri M, et al. Evaluation of menstrual irregularities after COVID-19 vaccination: Results of the MECOVAC survey. Open Med (Wars). 2022;17(1):475–84. https://doi.org/10.1515/med-2022-0452 PMID: 35350834.

    Article  CAS  Google Scholar 

  19. Harlow SD, Ephross SA. Epidemiology of menstruation and its relevance to women’s health. Epidemiol Rev. 1995;17(2):265–86. https://doi.org/10.1093/oxfordjournals.epirev.a036193 PMID: 8654511.

    Article  PubMed  CAS  Google Scholar 

  20. Yarde F, Voorhuis M, Dolleman M, Knauff EA, Eijkemans MJ, Broekmans FJ. Antimullerian hormone as predictor of reproductive outcome in subfertile women with elevated basal follicle-stimulating hormone levels: a follow-up study. Fertil Steril. 2013;100(3):831–8. https://doi.org/10.1016/j.fertnstert.2013.05.009 PMID: 23755952.

    Article  PubMed  CAS  Google Scholar 

  21. Vermesh M, Kletzky OA, Davajan V, Israel R. Monitoring techniques to predict and detect ovulation. Fertil Steril. 1987;47(2):259–64 PMID: 3817171.

    Article  CAS  Google Scholar 

  22. Kerin J. Ovulation detection in the human. Clin Reprod Fertil. 1982;1(1):27–54 PMID: 6821195.

    PubMed  CAS  Google Scholar 

  23. Hayaki J, Holzhauer CG, Epstein EE, Cook S, Gaba A, Lorenzo AC, et al. Menstrual cycle phase, alcohol consumption, alcohol cravings, and mood among women in outpatient treatment for alcohol use disorder. Psychol Addict Behav. 2020;34(6):680–9. https://doi.org/10.1037/adb0000576 PMID: 32250129.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Allen AM, McRae-Clark AL, Carlson S, Saladin ME, Gray KM, Wetherington CL, et al. Determining menstrual phase in human biobehavioral research: A review with recommendations. Exp Clin Psychopharmacol. 2016;24(1):1–11. https://doi.org/10.1037/pha0000057 PMID: 26570992.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Simmons RG, Jennings V. Fertility awareness-based methods of family planning. Best Pract Res Clin Obstet Gynaecol. 2020;66:68–82. https://doi.org/10.1016/j.bpobgyn.2019.12.003 PMID: 32169418.

    Article  PubMed  Google Scholar 

  26. Kim HY. Statistical notes for clinical researchers: assessing normal distribution (2) using skewness and kurtosis. Restor Dent Endod. 2013;38(1):52–4. https://doi.org/10.5395/rde.2013.38.1.52 PMID: 23495371.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Choi SW. Life is lognormal! What to do when your data does not follow a normal distribution. Anaesthesia. 2016;71(11):1363–6. https://doi.org/10.1111/anae.13666 PMID: 27734487.

    Article  PubMed  CAS  Google Scholar 

  28. Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Softw. 2015;67(1):1–48. https://doi.org/10.18637/jss.v067.i01.

  29. Bai X, Li J, Zhou L, Li X. Influence of the menstrual cycle on nonlinear properties of heart rate variability in young women. Am J Physiol Heart Circ Physiol. 2009;297(2):H765-74. https://doi.org/10.1152/ajpheart.01283.2008 PMID: 19465541.

    Article  PubMed  CAS  Google Scholar 

  30. Zhu TY, Rothenbuhler M, Hamvas G, Hofmann A, Welter J, Kahr M, et al. The Accuracy of Wrist Skin Temperature in Detecting Ovulation Compared to Basal Body Temperature: Prospective Comparative Diagnostic Accuracy Study. J Med Internet Res. 2021;23(6):e20710. https://doi.org/10.2196/20710 PMID: 34100763.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Regidor PA, Kaczmarczyk M, Schiweck E, Goeckenjan-Festag M, Alexander H. Identification and prediction of the fertile window with a new web-based medical device using a vaginal biosensor for measuring the circadian and circamensual core body temperature. Gynecol Endocrinol. 2018;34(3):256–60. https://doi.org/10.1080/09513590.2017.1390737 PMID: 29082805.

    Article  PubMed  CAS  Google Scholar 

  32. Faust L, Bradley D, Landau E, Noddin K, Farland LV, Baron A, et al. Findings from a mobile application-based cohort are consistent with established knowledge of the menstrual cycle, fertile window, and conception. Fertil Steril. 2019;112(3):450-7 e3. https://doi.org/10.1016/j.fertnstert.2019.05.008 PMID: 31272722.

    Article  PubMed  Google Scholar 

  33. Ecochard R, Duterque O, Leiva R, Bouchard T, Vigil P. Self-identification of the clinical fertile window and the ovulation period. Fertil Steril. 2015;103(5):1319-25 e3. https://doi.org/10.1016/j.fertnstert.2015.01.031 PMID: 25724738.

    Article  PubMed  Google Scholar 

  34. Johnson S, Stanford JB, Warren G, Bond S, Bench-Capon S, Zinaman MJ. Increased Likelihood of Pregnancy Using an App-Connected Ovulation Test System: A Randomized Controlled Trial. J Womens Health (Larchmt). 2020;29(1):84–90. https://doi.org/10.1089/jwh.2019.7850 PMID: 31483187.

    Article  PubMed Central  Google Scholar 

  35. Boutzios G, Karalaki M, Zapanti E. Common pathophysiological mechanisms involved in luteal phase deficiency and polycystic ovary syndrome. Impact on fertility. Endocrine. 2013;43(2):314–7. https://doi.org/10.1007/s12020-012-9778-9 PMID: 22930247.

    Article  PubMed  CAS  Google Scholar 

  36. Dishi M, Enquobahrie DA, Abetew DF, Qiu C, Rudra CB, Williams MA. Age at menarche, menstrual cycle characteristics and risk of gestational diabetes. Diabetes Res Clin Pract. 2011;93(3):437–42. https://doi.org/10.1016/j.diabres.2011.07.001 PMID: 21816498.

    Article  PubMed  Google Scholar 

  37. Tatsumi T, Sampei M, Saito K, Honda Y, Okazaki Y, Arata N, et al. Age-Dependent and Seasonal Changes in Menstrual Cycle Length and Body Temperature Based on Big Data. Obstet Gynecol. 2020;136(4):666–74. https://doi.org/10.1097/AOG.0000000000003910 PMID: 32925608.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Vallejo M, Marquez MF, Borja-Aburto VH, Cardenas M, Hermosillo AG. Age, body mass index, and menstrual cycle influence young women’s heart rate variability –a multivariable analysis. Clin Auton Res. 2005;15(4):292–8. https://doi.org/10.1007/s10286-005-0272-9 PMID: 16032384.

    Article  PubMed  Google Scholar 

  39. McGovern PG, Myers ER, Silva S, Coutifaris C, Carson SA, Legro RS, et al. Absence of secretory endometrium after false-positive home urine luteinizing hormone testing. Fertil Steril. 2004;82(5):1273–7. https://doi.org/10.1016/j.fertnstert.2004.03.070 PMID: 15533341.

    Article  PubMed  Google Scholar 

  40. Behre HM, Kuhlage J, Gassner C, Sonntag B, Schem C, Schneider HP, et al. Prediction of ovulation by urinary hormone measurements with the home use ClearPlan Fertility Monitor: comparison with transvaginal ultrasound scans and serum hormone measurements. Hum Reprod. 2000;15(12):2478–82. https://doi.org/10.1093/humrep/15.12.2478 PMID: 11098014.

    Article  PubMed  CAS  Google Scholar 

  41. Stanford JB. Revisiting the fertile window. Fertil Steril. 2015;103(5):1152–3. https://doi.org/10.1016/j.fertnstert.2015.02.015 PMID: 25772771.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the doctor who performed the ovarian ultrasound for the participants and the nurse who collected the blood samples. In addition, we thank the engineers of the Information Department of the IPMCHH for setting up the menstrual cycle monitoring project in the hospital outpatient system, which greatly facilitated monitoring of ovulation.

Funding

This research is supported by National Key Research and Development Program of China (2021YFC2700701), CAMS Innovation Fund for Medical Sciences (2019-I2M-5-064), National Natural Science Foundation of China (81661128010, 82171686), the International Science and Technology Collaborative Fund of Shanghai (18410711800), Program of Shanghai Academic Research Leader (20XD1424100), Clinical Research Plan of Shanghai Shenkang Hospital Development Center (SHDC12018X17, SHDC12019107, SHDC2020CR1008A), Science and Technology Innovation Fund of Shanghai Jiao Tong University (YG2019GD04, YG2020YQ29), National Natural Science Foundation of China (82088102), Collaborative Innovation Program of Shanghai Municipal Health Commission (2020CXJQ01), Technology Innovation Action Plan Hong Kong, Macao and Taiwan Science and Technology Cooperation Project (19410760100), Shanghai Frontiers Science Research Base of Reproduction and Development. Huawei (Huawei Device Co, Ltd, Shenzhen, China) provided the smartphones (Huawei Mate 30), smart bands (Huawei Band 5), and thermometers (Braun IRT6520) for the study.

Author information

Author notes

  1. Jia-Le Yu and Yun-Fei Su contributed equally to this work.

Authors and Affiliations

  1. The International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China

    Jia-Le Yu, Yun-Fei Su & Lu-Ting Chen

  2. Shanghai Key Laboratory of Embryo Original Diseases, Shanghai, China

    Jia-Le Yu, Yun-Fei Su & Lu-Ting Chen

  3. Shanghai Municipal Key Clinical Speciality, Shanghai, China

    Jia-Le Yu, Yun-Fei Su & Lu-Ting Chen

  4. Obstetrics and Gynecology Hospital, Institute of Reproduction and Development, Fudan University, No. 419, Fangxie Rd, Shanghai, 200011, China

    Chen Zhang, Li Jin, Xian-Hua Lin, He-Feng Huang & Yan-Ting Wu

  5. Research Units of Embryo Original Diseases, Chinese Academy of Medical Sciences (No. 2019RU056), Shanghai, China

    He-Feng Huang & Yan-Ting Wu

  6. Key Laboratory of Reproductive Genetics (Ministry of Education), Department of Reproductive Endocrinology, Women’s Hospital, Zhejiang University School of Medicine, Hangzhou, China

    He-Feng Huang

Authors
  1. Jia-Le Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yun-Fei Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chen Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Li Jin
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xian-Hua Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Lu-Ting Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. He-Feng Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Yan-Ting Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

He-Feng Huang and Yan-Ting Wu designed the study and were responsible for quality control. Jia-Le Yu and Yun-Fei Su conducted the statistical analysis and drafted the manuscript. Chen Zhang provided oversight for the statistical analysis and in drafting the manuscript. Jia-Le Yu and Lu-Ting Chen were responsible for data collection. He-Feng Huang, Yan-Ting Wu, Li Jin and Xian-Hua Lin critically revised the manuscript. All authors provided administrative, technical, and material support and agreed with the final version of the article.

Corresponding authors

Correspondence to He-Feng Huang or Yan-Ting Wu.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Medical Research Ethics Committee of the IPMCHH in Shanghai, China (GKLW2020-06). All research activities were carried out in line with the principles of the Declaration of Helsinki. Written informed consent was obtained from each participant before the initiation of any study procedures.

Consent for publication

All co-authors have read and approved the final version of this manuscript and have agreed to its submission for publication. All participants consented having their data published.

Competing interests

None.

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: FigureS1.

Physiological parameters of the irregular group in the different phases of the menstrual cycle. Table S1. Characteristics of all participants. Table S2. The relationship between menstrual cycle and physiological parameters in the irregular group. Table S3. Accuracy, sensitivity and specificity of different models in fertile window prediction. Table S4. Performance of prediction model developed using BBT and HR data of regular group in fertile window prediction in different irregular subgroups. Table S5. Performance of prediction model developed using BBT and HR data of regular group in fertile window prediction in different age subgroups. Supplementary Methods.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Yu, JL., Su, YF., Zhang, C. et al. Tracking of menstrual cycles and prediction of the fertile window via measurements of basal body temperature and heart rate as well as machine-learning algorithms. Reprod Biol Endocrinol 20, 118 (2022). https://doi.org/10.1186/s12958-022-00993-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12958-022-00993-4

Keywords

Reproductive Biology and Endocrinology

ISSN: 1477-7827

Contact us