Advertisement

Statistical analysis of interfraction dose variations of HRCTV and OARs for cervical cancer HDR brachytherapy

Open AccessPublished:July 18, 2022DOI:https://doi.org/10.1016/j.adro.2022.101019

      Purpose

      High dose rate (HDR) brachytherapy for cervical cancer treatment includes significant uncertainties. This study was to quantify the interfraction dosimetric variation (IDV) of the high-risk clinical target volume (HRCTV) from the prescribed dose and the corresponding effect on organ at risk (OAR) dose based on a comprehensive statistical analysis.

      Methods and Materials

      Fifty cervical cancer patients treated with HDR intracavity brachytherapy (ICBT) from October 2019 to December 2020 were retrospectively analyzed. The OARs of interest were the rectum, bladder, sigmoid, and bowel. The dosimetric parameters evaluated for all patients was the dose absorbed by 90% of the HRCTV (D90) and the dose absorbed by 0.1 (D0.1cc) and 2 cm3 (D2cc) of each respective OAR. The HRCTV variations were from the prescribed dose and the OAR variations were from the corresponding tolerance dose. Distribution fitting of the HRCTV variations was determined to quantify the IDV. Comparative statistics of the HRCTV variations with the OAR variations were conducted to determine correlations.

      Results

      The mean HRCTV variation from the prescribed dose was -2.53% ± 8.74%. The HRCTV variations and OAR variations showed moderate to weak linear correlations despite the variations being relative to each other, with the bladder D2cc having the strongest correlation. There was a 30.0% (2.62% 95% confidence interval) probability of under-dosing the HRCTV (-5% variation from prescription) and a 23.3% (2.62% 95% CI) probability of overdosing the HRCTV (+5% variation from prescription). This tendency to under-dose the HRCTV was a consequence of HRCTV IDV not being normally distributed.

      Conclusions

      HRCTV dosimetric variations and OAR variations were complexly correlated with the bladder D2cc having the strongest correlation. HRCTV IDV were best described as a left skewed distribution that indicates a tendency of under-dosing the HRCTV. The clinical significance of such dose variations is expected and will be further investigated.

      Introduction

      Brachytherapy procedures are subject to varying levels of uncertainties, from source construction and calibration to delivery of clinical plans [1-3]. These uncertainties can result from technology or clinical procedures. Uncertainties associated with clinical procedures that clinicians have control over are called clinical uncertainties. Clinical uncertainties include the uncertainties of structure delineation and organ motion. [4-14]
      It is accepted that HRCTV delineation and OAR motion are the most significant components of the brachytherapy treatment uncertainty budget [15]. At our institution, HDR T&O ICBT treatments are online adaptive procedures: every fraction has a new CT scan and plan which results in interfraction dosimetric variations (IDV) of the HRCTV, especially, IDVs from the prescribed dose. Although dosimetric variations may not be considered an uncertainty in statistical terms, they are considered a form of uncertainty in brachytherapy [16]. IDVs from the prescribed dose may not be dominant among HDR T&O ICBT uncertainties, but they are important as they may have a significant impact on clinical outcomes. However, this type of uncertainty has not been studied with the same rigor as the aforementioned uncertainties.
      There have been studies of IDVs of OARs and the target volume or point. [12, 14, 17]. These studies acknowledge that IDVs are forms of uncertainties in brachytherapy, but mostly focus on deformable image registration (DIR) dosimetric parameters and their variation from dose-volume histogram (DVH) dosimetric parameters, not the variation of dose from the given prescription. Chakraborty et al. [18] and Jamema et al. [19] studied the effect of interfraction applicator position on OAR dose in cervical cancer brachytherapy in addition to the spatial change of the dosimetric parameters. However, neither focused on the IDV of the HRCTV nor its corresponding effect on OAR dose. Sharma et al studied Point A dosimetric variations from the given prescription in fractionated brachytherapy [20]. Despite the importance of continuing to use point doses in modern day cervical cancer brachytherapy, volumetric parameters such as the HRCTV have taken priority to Point A and other point dose parameters [21].
      To our knowledge, HRCTV dosimetric variations from the given prescription and the corresponding effect on OAR dose in HDR T&O ICBT have not been studied. Therefore, in this study we evaluate the IDV of the HRCTV from the prescribed dose and the corresponding effects on OAR dose in HDR T&O ICBT. Furthermore, we studied the distribution of IDVs from the prescribed dose to quantify the corresponding uncertainty.

      Methods and materials

      Data collection

      Fifty patients diagnosed with cancers of the uterine cervix and treated with HDR T&O ICBT from October 2019 to December 2020 were retrospectively analyzed. All patients were treated with prescriptions of 5 or 7 Gy per fraction for 2 to 5 fractions. There was a total of 188 fractions of HDR T&O ICBT evaluated. The delineation of structures followed the International Commission on Radiation Units report 89 (ICRU 89) and was conducted on CT images[22]. The HRCTV was delineated as the entire cervix, uterus, parametrium, and vagina. The OARs of interest were the rectum, bladder, sigmoid colon, and the bowel. Treatment planning for each fraction was conducted in the Varian Eclipse brachytherapy treatment planning system (TPS).

      Statistical analysis

      Dosimetric variation

      The calculation of dosimetric variations was done using the percent difference equation:
      %difference=DosedeliveredPrescriptionortolerancedosePrescriptionortolerancedose×100
      (1)


      where the dose delivered is D90 for the HRCTV, and dose in 0.1 cc (D0.1cc) and 2 cc (D2cc) for the OARs. The tolerance dose is 80% of the prescription for the OARs. The dosimetric variations were organized by structure for each patient. The mean dosimetric variation over a patient's course of HDR T&O ICBT constituted as a data point for each structure. Thus, for each structure there were 50 data points evaluated for the 50 patients in the study. This provided a description of each patient's IDV from the given prescription.
      Basic statistics
      Data analysis was performed in Python 3.7 via the use of the SciPy and DistFit packages. The mean, standard deviation, and median were calculated for each structure's dosimetric dataset. From the mean and standard deviation, the coefficient of variation (CV) was determined for each structure. The CV is a measurement of consistency in the data: the higher the CV the less consistent the data is, the lower the CV the more consistent the data is. The 95% confidence interval (CI) was calculated for all relevant parameters.
      Correlations
      Spearman's correlation coefficient (ρ) was calculated for each OAR dosimetric parameter against HRCTV D90. Spearman's correlation coefficient was used to determine if a non-linear, monotonic relationship existed between the relative variations in addition to linear relationship detection. As the correlation coefficient approaches -1 or 1, then the OAR dosimetric variation is correlated to the HRCTV dosimetric variation with a descending or ascending slope, respectively. Statistical significance for correlation coefficients was determined as p-values < 0.05.
      Comparison of median or mean variations using Wilcoxon signed-rank test or Student's t-test was not necessary for this study. This is due to the values of the OAR variations from tolerance having limited clinical significance, making the comparison of median or mean OAR variations to HRCTV variation obsolete.
      Data driven distribution fitting
      The HRCTV distribution was fitted to 89 different distributions using the DistFit function in Python 3.7 to find the best fit distribution. The histogram bin width can affect the fitted distribution. Therefore, distribution fitting was performed with limited dependence on histogram bin width to obtain an accurate fit [23]. To accurately fit a distribution to the HRCTV, the raw variation data was plotted as an empirical distribution, analogous to a line histogram. Each respective probability density function (PDF) was plotted along with the empirical distribution of the HRCTV. Distribution fits were ranked according to their residual sum of squares (RSS) score: the lower the score the better the fit. The RSS is the sum of squared distances from a given point on the empirical distribution curve to the corresponding point on the PDF curve. The RSS equation is shown in Equation (2):
      RSS=in(yif(xi))2
      (2)


      Where n is the maximum data point, yi is the ith point of the empirical distribution, and f(xi) is the ith point of the PDF. There were 50 evaluation points to fit and score the distributions.
      The RSS is only a relative measurement parameter and does not determine the statistical significance of a fitted distribution. Therefore, an Anderson-Darling (AD) test was performed to determine whether the fits determined by the RSS were statistically significant [24]. The AD test uses a distribution specific term to calculate the test statistic and serves mostly as a test of normality. AD test p-values > 0.05 or < 0.05 indicate that the data does or does not fit the distribution, respectively. Additionally, quantile-quantile (Q-Q) plots of the datasets were tabulated for visual interpretations of the AD test results.
      Probability of clinically significant HRCTV D90 variations from the given prescription calculations
      Each fitted distribution has a corresponding cumulative distribution function (CDF(x)) and survival function (SF(x) = 1-CDF(x)) as a function of dose variation of x. Evaluating CDFs at some desired value gives the probability of the variable obtaining a value less than or equal to “x”. The same can be said for calculating probabilities greater than or equal to some desired value using a distribution's survival function, which is one minus a distribution's CDF.
      The American Brachytherapy Society (ABS) estimated an HDR dose variation of ±0.25 Gy per fraction from a given prescription [25]. This equates to a ±5% variation for 5 Gy per fraction and a ±3.57% variation for 7 Gy per fraction. This is consistent with a clinically significant dose tolerance of ±5% for radiation therapy. Where -5% variations are under-dosed conditions and +5% variations are overdosed conditions. The probability of a patient's treatment course resulting in an under or overdose of the HRCTV was calculated using the best fitted distribution's CDF and survival function as shown in Equations (3) and (4):
      Underdoseprobability(%)=CDF(5%)×100
      (3)


      Overdoseprobability(%)=(1CDF(5%))×100
      (4)


      Robust distribution fitting
      The shaping parameters that define the fitted distributions will change with additional data. Thus, the fitted distributions determined in this study from the RSS score calculations are not robust to model different datasets. To compensate for this, we also used a non-parametric approach to model the data. The data were dichotomized, where 0 was indicative of an under or overdose incident and 1 a non-under or overdose incident. The proportion of incidents was used to determine the probability of under and overdosing the HRCTV as well as the probability of clinically significant variations (under-dose + overdose probability). This approach provided robust estimations of the probabilities of interest without imposing any distributional assumption on the data.
      Results
      Basic Statistics
      Fifty patients and a total of 188 fractions were analyzed for D90 variations from the given prescription and D2cc and D0.1cc variations from the corresponding tolerance dose. The mean variation for the HRCTV was -2.53% ± 8.74%, ranging from -25.4% to 8.90%, and the variations of up to > 20% (D2cc) and > 65% (D0.1cc) for the bladder were found, as indicated in Figure 1. The mean and CV for all structures are tabulated in Table 1. The D0.1cc CV is larger than the corresponding D2cc CV for all OARs except the bladder. The D0.1cc is accepted as a less robust parameter for dose reporting when compared to the D2cc [13,21]. The lower D0.1cc CV implies that the D0.1cc is more robust than the expected D2cc for the bladder when evaluated as mean course variations from tolerance. Distributions of all evaluated structures are displayed in Figure 1. By inspection, the OAR distributions appear to take different shapes than the HRCTV distributions despite the variations being relative to each other.
      Figure 1:
      Figure 1Histograms of all evaluated structures (bin width = 5%). Plot A is the HRCTV D90 distribution, plots B and C are the rectum D2cc and D0.1cc, plots D and E the bladder D2cc and D0.1cc, plots F and G the sigmoid D2cc and D0.1cc, and plots H and I the bowel D2cc and D0.1cc.
      Table 1The mean and coefficient of variation (CV) for HRCTV D90, and OAR D0.1cc and D2cc.
      Structure and parameterMean ± σ variation (%)Coefficient of variation (proportion)
      HRCTV D90-2.53 ± 8.743.45
      Rectum D2cc-25.9 ± 13.90.54
      Rectum D01.cc-1.55 ± 18.712.1
      Bladder D2cc-5.86 ± 12.02.05
      Bladder D0.1cc28.4 ± 19.30.68
      Sigmoid D2cc-28.3 ± 16.21.75
      Sigmoid D0.1cc-1.73 ± 22.312.9
      Bowel D2cc-49.8 ± 23.60.47
      Bowel D0.1cc-30.6 ± 34.61.13
      Nineteen of the 50 patients had their HDR T&O ICBT course result in an average under-dosing of the HRCTV with a mean variation of -11.6%. Ten of the 50 patients had their HDR T&O ICBT course result in an average overdosing of the HRCTV with a mean variation of 7.45%.
      Correlations
      The HRCTV D90 variations from the given prescription and OAR variations from the corresponding tolerance showed moderate to weak linear correlations. The rectum D2cc, bladder D2cc and D0.1cc showed statistically significant linear correlations (r = -0.305, -0.427, -0.373 and p-value = 0.033, 0.002, 0.008). Non-linear correlations also existed for both datasets. Table 2 has the tabulated Spearman's correlation coefficients (ρ) for both datasets. The non-linear correlations showed similar statistical strength as the linear correlations. The rectum D0.1cc and both bowel dosimetric parameters are the only variations consistently not correlated to HRCTV variations from prescription. The rectum D0.1cc was only moderately uncorrelated (ρ =-0.230, p-value = 0.111), whereas both bowel dosimetric parameters were strongly uncorrelated (p-value > 0.250). The bladder D2cc had the strongest correlation with HRCTV variations (ρ = -0.508, p < 0.001); the correlation is displayed in Figure 2. The sigmoid dosimetric parameters were not correlated HRCTV variations (p-value = 0.159 and p-value > 0.250 for D2cc and D0.1cc, respectively).
      Table 2Non-linear correlations.
      Structure and parameterSpearman's correlation coefficient (ρ)p-value
      Rectum D2cc-0.3020.035
      Rectum D01.cc-0.2300.111
      Bladder D2cc-0.508< 0.001
      Bladder D0.1cc-0.3880.006
      Sigmoid D2cc-0.2040.159
      Sigmoid D0.1cc-0.155> 0.250
      Bowel D2cc0.042> 0.250
      Bowel D0.1cc0.017> 0.250
      Figure 2:
      Figure 2Scatter plot of bladder D2cc vs. HRCTV D90. The HRCTV D90 variations from prescription are on the x-axes and the corresponding bladder D2cc variations from tolerance are on the y-axes.
      Fitted distributions
      HRCTV D90 variations from the prescription are not normally distributed (AD p-value < 0.001). Figure 3 displays the best fit distributions for the HRCTV D90 determined from the data driven analysis (RSS scores), and the corresponding normal distributions if the variations were normally distributed. The Generalized-extreme value (GEV) distribution is the best fit distribution for HRCTV D90 (AD p-value = 0.213). It is important to note that the distribution of HRCTV D90 variations from the prescription being left skewed and non-normal is more important than the actual fitted distribution, for the best fitted distribution may change with varying amounts of data.
      Figure 3:
      Figure 3Best fitted distributions and a fitted normal distribution for HRCTV D90 (A & C) determined from the data driven analysis (RSS score). The HRCTV D90 variations from prescriptions are on the x-axis and the corresponding probability density is on the y-axis. Q-Q plots and Anderson Darling (AD) test results are also tabulated (B & D) to show the statistical strength, or lack thereof, of the fitted distributions. The respective fitted distribution quantiles are on the x-axis and the HRCTV D90 variations from prescription quantiles are on the y-axis. AD p-values > 0.05 mean the distribution statistically fits the data. AD p-values < 0.05 mean the distributions do not statistically fit the data. From the Ad test and Q-Q plots, the variations are not normally distributed.
      The calculated probabilities from the GEV and normal distributions’ CDF and SF of under (-5% HRCTV D90 variation) and overdosing (+5% HRCTV D90 variation) the HRCTV, respectively, are tabulated in Table 3. The respective distributions’ corresponding mean, median, and standard deviations are also tabulated in Table 3. The GEV distribution had a higher probability of under-dosing the HRCTV (30.0%), when compared to overdosing the HRCTV (23.3%). Figure 4 provides a visual of the HRCTV under and overdosing probabilities calculated from the GEV's CDF.
      Table 3Probabilities and statistics from distributions. The probabilities were calculated using the Generalized extreme-value and normal CDF determined from their RSS scores.
      Under-dose probability (%)Overdose Probability (%)Significant variation probability (%)Mean ± Std (%)Median (%)
      Generalized-extreme Value30.023.353.3-2.32 ± 9.460.15
      Normal41.218.759.9-3.00 ± 9.00-3.00
      Figure 4:
      Figure 4Visualized probabilities of under and overdosing the HRCTV Generalized-extreme value CDF. HRCTV D90 variations from the given prescriptions are on the x-axis, and the probability of a corresponding variation from prescription is on the y-axis. Under-dosing the HRCTV was defined as HRCTV D90 variations less than -5%, overdoing the HRCTV was defined as HRCTV D90 variations greater than 5% [29].
      The non-parametric distributions are displayed in Figure 5. The under-dose, overdose, and clinically significant variation probability was 38.0%, 20.0%, and 58.0%, respectively. This resulted in an 8.0% difference for under-dose, 3.3% for overdose, and 4.7% for clinically significant variation when compared to the GEV distribution.
      Figure 5:
      Figure 5Discrete Bernoulli distributions for under-dose, overdosed, and clinically significant HRCTV D90 variations from the prescribed dose. On the x-axis, a value of 0 is indicative of an under or overdose incident, while a value of 1 is indicative of a non-under or overdose incident. On the y-axis, the probability mass function (PMF) gives the probability of the observed events. Plots A, B, and C are the under-dose, overdose, and clinically histograms.
      Tables 4 has this study's value-relevant parameter's 95% CI tabulated. The value-relevant parameters are the mean HRCTV D90 variation from the given prescription, the under-dose probability, and overdose probability. The clinically significant probability is implied from the under-dose and overdose probabilities. The value of these parameters holds clinical significance because they give the variation from the prescribed dose and the probability of significant variations from said dose. The value of the OAR variations holds limited clinical significance because dose is not prescribed to OARs, and it is a goal to limit OAR dose as much as possible. The 95% CI is tabulated for both distribution fitting techniques: data driven (RSS score) and non-parametric (robust).
      Table 495% Confidence intervals for HRCTV D90 variations from prescription.
      DatasetVariation from prescription
      HRCTV D90 mean variation-2.53% ± 2.42%
      Data driven under-dose probability30.0% ± 2.62%
      Data driven overdose probability23.3% ± 2.62%
      Non-parametric under-dose probability38.0% ± 2.42%
      Non-parametric overdose probability20.0% ± 2.42%
      Discussion
      The variation of HRCTV dose is an important issue as it may have potential effect on clinical outcomes. However, it has not been well addressed in the literature. Sharma et al studied IDVs of the target from the prescription given at Point A and the IDV of OAR dose for point dose parameters [20]. They found that the average IDV of Point A doses from the given prescription was 1.55% ± 1.07%. In this study, HRCTV dosimetric variations from the prescription and the corresponding effect on OAR dose have been successfully evaluated., Also, large HRCTV IDVs (-2.53% ± 8.74%, from -25.4% to 8.9%) and large OAR IDVs (e.g., up to > 20% (D2cc) and > 65% (D0.1cc) for the bladder) were obtained. Although the determination of clinical effect of IDVs should be based on clinical data and is beyond the scope of this paper, the clinical significance of the results in this study was anticipated and can be estimated based on certain models. Estimated using the dose response curves proposed by Tanderup, et al. [15], up to -9.1% change in local control and 12.4% change in morbidity could be caused by the aforementioned IDVs. This estimate may not be accurate, but at least indicates such IDVs may have significant effect on clinical outcomes. More thorough and systematic analysis will be performed based on clinical data in our future studies. Complex correlations between OAR doses and HRCTV D90 were found in this study. A linear correlation would indicate that a simple relationship between the respective variations is evident. That is, the cause of the OAR variations can simply be explained from the HRCTV variations. This is not the case. The variations show a stronger monotonic correlation than linear correlations. Monotonic functions are statistical functions with limited relevance clinically. Rather, it is stated that the relative variations are complexly correlated: there is more to the cause of the changing OAR dose than just the HRCTV dosimetric variation from prescription, despite the two being relative to each other.
      Uncertainties in brachytherapy are assumed to be random and, thus, normally distributed [16,29]. Nesvacil et al studied the simulated effect of systemic and random uncertainties on tumor control probability (TCP) and normal tissue complication probability (NTCP) models [16]. Systematic uncertainties were defined as consistent errors that are out of the control of clinicians, and random uncertainties were defined as dosimetric variations. They found that that TCP and NTCP models were generally robust to varying degrees of random uncertainties when combined with consistent systematic uncertainties. However, we have found that HDR T&O ICBT course dosimetric variations are not normally distributed, and, thus, cannot be assumed as a random uncertainty. The distribution of HRCTV variations is left-skewed, meaning there is a higher probability of under-dosing the HRCTV than overdosing the HRCTV. Assuming a normal distribution would result in either equal probabilities of under and overdosing the HRCTV, or overestimate the under-dose probability and underestimate the overdose probability as we have seen from the fitted normal distribution (Table 3). Both the data-driven distribution (GEV) determined from the RSS-score (30.0% ± 2.62%, 95% CI) and robust distribution (38.0% ± 2.42%, 95% CI) supports the claim that there is a tendency to under-dose the HRCTV throughout a patient's course of treatment, thus meaning the distribution of variations is left-skewed and non-normal. The observed non-random effect of dosimetric variations on TCP and NTCP models and clinical outcomes is recommended for future studies.
      The large IDVs found in this study indicated that delivering the prescribed dose to the target while sparing OARs is not always obtainable. Any techniques that can improve target coverage and OAR sparing should be encouraged to apply in clinical practice. Recently, clinical trials of hyaluronate gel injection spacers between the vagina and rectum have shown promising results in reducing rectum dose without sacrificing tumor coverage in GYN brachytherapy [26,27]. We recognize that different prescription doses and different fractionations may have different OAR dose constrains with different EQD2 values. However, the effect of differing dose per fraction is reduced in this study because of the use of tolerance doses relative to prescription doses. The 80% tolerance are generally conservative and provides a standard and uniform analysis [25,28,30]. The clinical implications of differing dose prescription will be studied extensively in future research.
      Interfraction contour variability and OAR motion may also affect dosimetric variations. However, interfraction contour variability and OAR motion were different uncertainties and not the focus of this study. In this study, we accepted the provided contours as the true anatomy and ignored possible OAR motion.
      In this study we only evaluated the correlations of HRCTV dosimetric variations on the D0.1cc and D2cc. However, for larger volume organs such as the sigmoid and the bowel, the D5cc and D10cc via dose surface histograms (DSH) are of clinical interest and are recommended for study purposes by ICRU 89 [22]. Volume coverage parameters such as the volume that receives 75% of the dose (V75) for OARs and the volume that receives 100% of the dose (V100) for HRCTV could also be used for evaluating HRCTV dosimetric variations from the given prescription [22]. Observations of the radiobiological effect of HRCTV dosimetric variations was not examined but will be examined in future studies. It is imperative to evaluate course variations using EQD2 to standardize the interfraction dosimetric variations in future studies. Doing so will strengthen the analysis as cumulative EBRT and brachytherapy doses are evaluated in this manner. A proper uncertainty analysis (adding uncertainties in quadrature) of the observed variations was not conducted in this study and will be included in future studies [2,31].
      Conclusion
      Dosimetric variations of the HRCTV from prescription and the corresponding effect on OAR dosimetric parameters were evaluated in this study. Complex correlations existed with HRCTV D90 variations from the given prescription and OAR dosimetric parameters. HRCTV D90 variations from the given prescription were well within the tolerance thresholds of ±5% in mean, but they formed a left skewed distribution best described by the Generalized-extreme value distribution that indicated an increased probability to exceed this tolerance with an increased tendency to under-dose the HRCTV. The clinical significance of such dose variations is expected and will be thoroughly and systematically investigated in future studies.
      Declaration of interests
      The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

      References

      [1] Rivard, M. J., Coursey, B. M., DeWerd, et al. (2004). Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. Medical physics31(3), 633-674.
      [2] DeWerd, L. A., Venselaar, J. L., Ibbott, G. S., et al. (2012). Overview on the dosimetric uncertainty analysis for photon-emitting brachytherapy sources, in the light of the AAPM Task Group No 138 and GEC-ESTRO report. Metrologia49(5), S253.
      [3] DeWerd, L. A., Ibbott, G. S., Meigooni, A. S., et al. (2011). A dosimetric uncertainty analysis for photon‐emitting brachytherapy sources: Report of AAPM Task Group No. 138 and GEC‐ESTRO. Medical physics38(2), 782-801.
      [4] Liu, H., Kinard, J., Maurer, J., et al. (2018). Evaluation of offline adaptive planning techniques in image‐guided brachytherapy of cervical cancer. Journal of applied clinical medical physics19(6), 316-322.
      [5] Duane, F. K., Langan, B., Gillham, C., et al. (2014). Impact of delineation uncertainties on dose to organs at risk in CT-guided intracavitary brachytherapy. Brachytherapy13(2), 210-218.
      [6] Hellebust, T. P., Tanderup, K., Lervåg, C., et al. (2013). Dosimetric impact of interobserver variability in MRI-based delineation for cervical cancer brachytherapy. Radiotherapy and Oncology107(1), 13-19.
      [7] Bell, L., Holloway, L., Bruheim, K., et al. (2020). Dose planning variations related to delineation variations in MRI-guided brachytherapy for locally advanced cervical cancer. Brachytherapy19(2), 146-153.
      [8] Saarnak, A. E., Boersma, M., van Bunningen, B. N., et al. (2000). Inter-observer variation in delineation of bladder and rectum contours for brachytherapy of cervical cancer. Radiotherapy and Oncology56(1), 37-42.
      [9] Arnesen, M. R., Bruheim, K., Malinen, E., et al. (2014). Spatial dosimetric sensitivity of contouring uncertainties in gynecological 3D-based brachytherapy. Radiotherapy and Oncology113(3), 414-419.
      [10] Petrič, P., Hudej, R., Rogelj, P., et al. (2013). Uncertainties of target volume delineation in MRI guided adaptive brachytherapy of cervix cancer: a multi-institutional study. Radiotherapy and Oncology107(1), 6-12.
      [11] Patel, S., Mehta, K. J., Kuo, H. C. G., et al. (2016). Do changes in interfraction organ at risk volume and cylinder insertion geometry impact delivered dose in high-dose-rate vaginal cuff brachytherapy?. Brachytherapy15(2), 185-190.
      [12] Kobayashi, K., Murakami, N., Wakita, A., et al. (2015). Dosimetric variations due to interfraction organ deformation in cervical cancer brachytherapy. Radiotherapy and Oncology117(3), 555-558.
      [13] Mazeron, R., Champoudry, J., Gilmore, J., et al. (2015). Intrafractional organs movement in three-dimensional image-guided adaptive pulsed-dose-rate cervical cancer brachytherapy: assessment and dosimetric impact. Brachytherapy14(2), 260-266.
      [14] Andersen, E. S., Noe, K. Ø., Sørensen, T. S., et al. (2013). Simple DVH parameter addition as compared to deformable registration for bladder dose accumulation in cervix cancer brachytherapy. Radiotherapy and Oncology107(1), 52-57.
      [15] Tanderup, K., Nesvacil, N., Pötter, R., et al. (2013). Uncertainties in image guided adaptive cervix cancer brachytherapy: impact on planning and prescription. Radiotherapy and Oncology107(1), 1-5.
      [16] Nesvacil, N., Tanderup, K., Lindegaard, J. C., et al. (2016). Can reduction of uncertainties in cervix cancer brachytherapy potentially improve clinical outcome?. Radiotherapy and Oncology120(3), 390-396.
      [17] van Heerden, L. E., van Wieringen, N., Koedooder, K., et al. (2018). Dose warping uncertainties for the accumulated rectal wall dose in cervical cancer brachytherapy. Brachytherapy17(2), 449-455.
      [18] Chakraborty, S., Patel, F. D., Patil, V. M., et al. (2014). Magnitude and implications of interfraction variations in organ doses during high dose rate brachytherapy of cervix cancer: a CT based planning study. International Scholarly Research Notices2014.
      [19] Jamema, S. V., Mahantshetty, U., Tanderup, K., et al. (2013). Inter-application variation of dose and spatial location of D2cm3 volumes of OARs during MR image based cervix brachytherapy. Radiotherapy and Oncology107(1), 58-62.
      [20] Brahmacharimayum Arunkumar Sharma, T., Singh, T., & Singh, L. J. (2011). Evaluation of variation in dose of organs at risk in intracavitary brachytherapy of cervical cancer–a prospective study. Journal of contemporary brachytherapy3(1), 23.
      [21] Mourya, A., Choudhary, S., Shahi, U. P., et al. (2021). A comparison between revised Manchester Point A and ICRU-89–recommended Point A definition absorbed-dose reporting using CT images in intracavitary brachytherapy for patients with cervical carcinoma. Brachytherapy20(1), 118-127.
      [22] The International Commission on Radiation Units and Measurements. (2016). Journal of the ICRU. Prescribing, Recording, and Reporting Brachytherapy for Cancer of the Cervix (No. 89). Oxford University Press. Volume 13 No 1-2 2013
      [23] Brown, L. D., & Hwang, J. G. (1993). How to approximate a histogram by a normal density. The American Statistician47(4), 251-255.
      [24] Engmann, S., & Cousineau, D. (2011). Comparing distributions: the two-sample Anderson-Darling test as an alternative to the Kolmogorov-Smirnoff test. Journal of applied quantitative methods6(3).
      [25] Nag, S., Erickson, B., Thomadsen, B., et al. (2000). The American Brachytherapy Society recommendations for high-dose-rate brachytherapy for carcinoma of the cervix. International Journal of Radiation Oncology* Biology* Physics48(1), 201-211.
      [28] Viswanathan, A. N., Thomadsen, B., & American Brachytherapy Society Cervical Cancer Recommendations Committee. (2012). American Brachytherapy Society consensus guidelines for locally advanced carcinoma of the cervix. Part I: general principles. Brachytherapy11(1), 33-46.
      [29] Tanderup, K., Hellebust, T. P., Lang, S., et al. (2008). Consequences of random and systematic reconstruction uncertainties in 3D image based brachytherapy in cervical cancer. Radiotherapy and Oncology89(2), 156-163.
      [30] Viswanathan, A. N., Beriwal, S., Jennifer, F., et al. (2012). American Brachytherapy Society consensus guidelines for locally advanced carcinoma of the cervix. Part II: high-dose-rate brachytherapy. Brachytherapy11(1), 47-52.
      [31] Kirisits, C., Rivard, M. J., Baltas, D., et al. (2014). Review of clinical brachytherapy uncertainties: analysis guidelines of GEC-ESTRO and the AAPM. Radiotherapy and oncology110(1), 199-212.