Skip to main content

Fitting NTCP models to bladder doses and acute urinary symptoms during post-prostatectomy radiotherapy



To estimate the radiobiological parameters of three popular normal tissue complication probability (NTCP) models, which describe the dose-response relations of bladder regarding different acute urinary symptoms during post-prostatectomy radiotherapy (RT). To evaluate the goodness-of-fit and the correlation of those models with those symptoms.


Ninety-three consecutive patients treated from 2010 to 2015 with post-prostatectomy image-guided intensity modulated radiotherapy (IMRT) were included in this study. Patient-reported urinary symptoms were collected pre-RT and weekly during treatment using the validated Prostate Cancer Symptom Indices (PCSI). The assessed symptoms were flow, dysuria, urgency, incontinence, frequency and nocturia using a Likert scale of 1 to 4 or 5. For this analysis, an increase by ≥2 levels in a symptom at any time during treatment compared to baseline was considered clinically significant. The dose volume histograms of the bladder were calculated. The Lyman-Kutcher-Burman (LKB), Relative Seriality (RS) and Logit NTCP models were used to fit the clinical data. The fitting of the different models was assessed through the area under the receiver operating characteristic curve (AUC), Akaike information criterion (AIC) and Odds Ratio methods.


For the symptoms of urinary urgency, leakage, frequency and nocturia, the derived LKB model parameters were: 1) D50 = 64.2Gy, m = 0.50, n = 1.0; 2) D50 = 95.0Gy, m = 0.45, n = 0.50; 3) D50 = 83.1Gy, m = 0.56, n = 1.00; and 4) D50 = 85.4Gy, m = 0.60, n = 1.00, respectively. The AUC values for those symptoms were 0.66, 0.58, 0.64 and 0.64, respectively. The differences in AIC between the different models were less than 2 and ranged within 0.1 and 1.3.


Different dose metrics were correlated with the symptoms of urgency, incontinence, frequency and nocturia. The symptoms of urinary flow and dysuria were poorly associated with dose. The values of the parameters of three NTCP models were determined for bladder regarding four acute urinary symptoms. All the models could fit the clinical data equally well. The NTCP predictions of urgency showed the best correlation with the patient reported outcomes.


The knowledge of dosimetric predictors for acute genito-urinary (GU) toxicity after intensity modulated radiotherapy (IMRT) for prostate cancer is largely lacking [1,2,3,4]. There is an increasing tendency to pay more attention to quality of life (QOL) issues, which is especially relevant in prostate cancer [5]. The incidence of acute moderate/severe GU toxicities, which play a major role on QOL, has increased in recent years, due to more aggressive treatments in terms of prescribed dose and fractionation schemes [1, 2].

Table 1 Dosimetric parameters in the patients with and without toxicity
Table 2 Summary of the best estimates and 95% confidence intervals of the parameters of the three examined normal tissue complication probability models for four acute bladder symptoms, respectively
Table 3 Summary of the results from the fit of the four normal tissue complication probability models for the different acute urinary symptoms during radiotherapy, respectively
Table 4 Summary of the cutoff values of gEUD / \( \overline{\overline{D}} \), which result in a statistically significant Odds Ratio (OR) larger than 1 for the different acute urinary symptoms
Fig. 1
figure 1

The dose-volume histograms (DVHs) of bladder for the acute symptom of urinary urgency. The DVHs of the patients with the symptom are represented by solid lines, whereas those without the symptoms by dashed lines. The thick lines represent the average DVHs, whereas the thin lines indicate the 68% confidence interval (one standard deviation)

Fig. 2
figure 2

The dose-response curves of bladder, for the acute symptoms of urinary urgency, urinary incontinence, frequency and nocturia for the three models. The unit on the dose axis is either the biologically effective uniform dose (\( \overline{\overline{D}} \) or BEUD), in the cases of the Relative Seriality (RS) model or the generalized equivalent uniform dose (gEUD) in the cases of the Lyman-Kutcher-Burman (LKB) and Logit models. The clinical response rates for different dose intervals (indicated by the horizontal error bars) are also plotted

Patient-reported health-related QOL is gaining recognition as an important outcome measure after cancer treatment [6,7,8]. The traditional approach in recording normal tissue toxicity in radiation therapy has been the use of physician-assessed scoring systems. However, there are multiple recent reports indicating that this type of assessments do not agree with patient-based reports and tend to underestimate the frequency and severity of treatment-related symptoms [9, 10]. For example, a recent study analyzed patients, who received concurrent chemoradiation therapy for rectal cancer. Although these patients showed the same grade level of diarrhea as assessed using the Common Terminology Criteria for Adverse Events (CTCAE), they had wide variations in their patient-reported diarrhea severities using a validated QOL instrument [11]. Patient-reported symptoms after different times post prostate cancer RT have been recently studied. However, this is not the case for the symptoms that develop during the course of RT treatment, involving modern intensity-modulated radiotherapy (IMRT) with image-guidance.

Modelling the relationship between dosimetric data and development of acute symptoms during treatment can ultimately lead to future radiation planning dose constraints, which may minimize these acute issues, making radiation treatment even safer and well-tolerated for patients.

At the department of Radiation Oncology of the University of North Carolina, QOL were prospectively acquired using a validated questionnaire [12] during weekly treatment visits as part of the routine clinical work-flow. The purpose of this study is to analyze dose/volume/outcome data in order to estimate the values of commonly-used NTCP models and their ability to accurately represent the data.


Patient data

This retrospective analysis included 93 consecutive patients, who received post-prostatectomy external beam radiation therapy (RT) between 2010 and 2015. All the patients received post-prostatectomy RT with image guidance using CT on-rails or Tomotherapy cone beam imaging. All the patients were treated with intensity modulated radiation therapy (IMRT) or volumetric modulated arc therapy (VMAT) using 6 MV photons. The prescribed dose to the planning target volume (PTV) was 66.6 Gy at 1.8 Gy per fraction delivered in 37 fractions for 7 weeks. Plans were designed to deliver at least 95% of the prescription dose to the PTV.

Patients completed the Urinary Obstruction/Irritation scale of the validated Prostate Cancer Symptom Indices (PCSI) [6, 7] before starting RT (baseline) and weekly during treatment visits. For this study seven datasets were available: 1) before or during first week of RT (93 pt); 2) second week of RT (83 pt); 3) third week of RT (84 pt); 4) forth week of RT (80 pt); 5) fifth week of RT (76 pt); 6) sixth week of RT (84 pt); and 7) seventh or last week of RT (87 pt); The patients that did not miss more than two visits and were analyzed per symptom were 88 for urine flow, 85 for frequency, 86 for nocturia, 88 for dysuria, 90 for urgency and 87 for incontinence. The symptom questions were rated on a Likert scale of 1 to 4 or 5, where 1 signified “not at all” and 4 or 5 signified “frequently or very frequently.” In this study, a score difference (maximum from all the weekly recordings) from baseline (score before RT) of ≥2 represented clinical significance [13, 14]. Six urinary symptoms were examined: 1) Urine flow; 2) Frequency; 3) Nocturia; 4) Dysuria; 5) Urgency; and 6) Incontinence. Among 93 analyzed patients, 6, 14, 15, 14, 26 and 13 had a ≥ 2 point increase in acute urinary symptoms related to urine flow, frequency, nocturia, dysuria, urgency and incontinence, respectively.

We hypothesized that bladder is the organ at risk where dose received may be associated with development of measured acute GU symptoms. Bladder delineation is part of the standard clinical practice in prostate cancer radiotherapy. The bladder dose volume histogram (DVH) was calculated for each patient from their treatment plans. These DVHs were correlated with the patient reported outcome data of different acute urinary symptoms. Figure 1 illustrates the distribution of the DVHs of the patients with and without the symptom of urgency.

Radiobiological models

The doses in the DVHs were converted to equivalent doses of 2 Gy per fraction (EDQ2Gy) based on the linear quadratic model using an α/β value of 3 Gy [15, 16].

$$ {\mathrm{EQD}}_{2\mathrm{Gy}}=D\cdot \left(\frac{1+\frac{d}{\alpha /\beta }}{1+\frac{2}{\alpha /\beta }}\right) $$

where D is the physical dose, d is the dose per fraction. A dose distribution can be reduced to a single dose value using the generalized Equivalent Uniform Dose (gEUD) concept as follows [17, 18]:

$$ {\mathrm{gEUD}}_{2\mathrm{Gy}}={\left(\sum \limits_i{\left({\mathrm{EQD}}_{2\mathrm{Gy}}\right)}_i^{1/n}\frac{V_i}{V_{\mathrm{tot}}}\right)}^n $$

where V i is the fractional subvolume of the organ being irradiated with a given dose and Vtot is the total volume of the organ. n is a parameter describing the volume dependence of the organ. gEUD2Gy is then used in the LKB and Logit models as follows [19]:

$$ \mathrm{NTCP}=\frac{1}{\sqrt{2\pi }}\underset{-\infty }{\overset{t}{\int }}{e}^{\frac{-{x}^2}{2}} dx $$


$$ t=\frac{{\mathrm{gEUD}}_{2 Gy}-{D}_{50}}{m\cdot {D}_{50}} $$

The Logit model is a logistic equation, which produces an analytical sigmoidal shaped curve commonly used in biology and it is defined in the following way [20]:

$$ \mathrm{NTCP}=\frac{1}{1+{\left(\frac{D_{50}}{gEUD_{2\mathrm{Gy}}}\right)}^k} $$

The Relative Seriality (RS) model used in this work applies a different approach to account for the volume dependence and it is mathematically expressed as follows [21, 22]:

$$ \mathrm{NTCP}={\left[1-\prod \limits_{i=1}^M{\left(1-P{\left({D}_i\right)}^s\right)}^{\Delta {v}_i}\right]}^{1/s} $$


$$ P\left({D}_i\right)=\exp \left[-{e}^{e\gamma -\left({EQD}_{2\mathrm{Gy}}^i/{D}_{50}\right)\cdot \left( e\gamma -\ln \ln 2\right)}\right] $$

where P(D i ) is the probability of response of an organ having the reference volume and being irradiated to dose D i . From the NTCP values the biologically effective uniform dose (\( \overline{\overline{D}} \)) can be derived by the following formula:

$$ P\left(\overrightarrow{D}\right)\equiv P\left(\overline{\overline{D}}\right)\Rightarrow \overline{\overline{D}}=\frac{e\gamma -\ln \left(-\ln \left(P\left(\overrightarrow{D}\right)\right)\right)}{e\gamma -\ln \left(\ln 2\right)} $$

The basic parameters of each model are: D50, which is the dose for a complication rate of 50%, the slope (gradient) of the dose response curve (m for LKB, k for Logit and γ for RS), and the parameter that accounts for the volume dependence of the organ (n for LKB and Logit, s for RS).

Statistical methods for fitting the NTCP models and evaluating the goodness-of-fit

The values of the parameters of the NTCP models and their 95% confidence intervals were determined using the maximum likelihood method. In this process, the predictions of the NTCP models were fitted to the clinical outcome results by changing the values of the model parameters until best estimates could be reached [23, 24]. The profile likelihood method was used to determine the confidence intervals of the model parameters. This method is often used when accurate interval estimates are difficult to obtain using standard methods (for example, when the log-likelihood function is highly non-normal in shape). For the 95% confidence region, the allowable difference from the maximum \( \ln {\overline{L}}_{\mathrm{max}}\left({\overline{x}}_i\right)-\ln {L}_{\mathrm{max}}\left({x}_i\right) \) is 1/2 × 3.84 = 1.92, for one degree of freedom.

The goodness-of-fit of the different NTCP models was assessed through the area under the receiver operating characteristic curve (AUC), maximum of the log-likelihood function, normal error distribution and Akaike information criterion (AIC) [25, 26]. More specifically, the AUC of a receiver operating characteristic (ROC) curve was computed and compared to the level of 0.5 (equivalent to a random predictor) [26]. The probability of achieving a worse fit was assessed by comparing the maximum log-likelihood value against the average log-likelihood value and its variance assuming a Gaussian distribution of the log-likelihood function [14, 23]. The Akaike information criterion was used to compare the fitting and complexity of the different models [25]. A lower Akaike number for a model indicates superiority of that model. Finally, the Odds Ratio (OR) method was applied to identify NTCP thresholds beyond which the risk of toxicity increases significantly.


Table 1 presents a summary of the average mean doses to bladder for the patients with and without acute urinary symptoms. Additionally, the average values of the dose-volume metrics that correlated best with the outcome date are shown. Urine flow and dysuria were not included in this table because the correlation between their dosimetric and outcome data was poor. Consequently, the values of the parameters of the NTCP models were determined for the symptoms of urgency, incontinence, frequency and nocturia (Table 2). The corresponding dose-response curves for the examined models and structures are shown in Fig. 2 Broadly, speaking, the three examined models essentially show equivalent goodness-of-fit per symptom. The analytical results are shown in Table 3. For the models and symptoms shown in this table, statistically significant ORs were identified for a given gEUD or \( \overline{\overline{D}} \) cutoffs in every case (Table 4). More specifically, for the symptom of urgency, the biological doses range between 59.7–61.6 for corresponding statistically significant ORs of 7.4–11.4. For frequency, the biological doses range between 55.9–57.7 for ORs of 4.4–9.4. For nocturia, the biological doses range between 40.9–52.5 for ORs of 4.5–6.6. For the symptom of incontinence, no statistically significant dose cutoffs and OR were found. The NTCP predictions against the actual response rates were in agreement (differences ranged between 0.0–0.1%) for all the models and symptoms. Regarding the correlation of the NTCP prediction against the actual outcome data, the symptom of urgency shows slightly higher AUC values (0.66), but worse AIC values (108.6–109.9). On the other hand, the symptom of frequency shows the best combination of AUC and AIC values (0.63–0.64 and 78.6–79.0, respectively) (Table 3).


It is well-known that radiation therapy causes an impact on QOL and specifically patient-reported urinary issues. Many patients after radical prostatectomy require radiation treatment to the prostate bed; these patients are especially vulnerable to urinary symptoms during radiation treatment due to the prior surgery. However, the association between doses received to the bladder and development of patient-reported acute urinary issues in this specific patient population has not been well-studied. To our knowledge, this is the first study to address this important knowledge gap.

Many of these studies use physician-reported toxicity, but patient-reported outcomes (PRO) have become an increasingly important way to measure quality of life (QOL) after cancer treatment [6,7,8,9,10]. Reports indicate that PRO-QOL data are often not in-line with physician-reported acute toxicity for prostate cancer patients [13, 27,28,29,30,31,32]. Additionally, the acute toxicity of post-prostatectomy radiation and dosimetric correlations continues to be understudied. Although, there are many studies investigating frequency and severity issues of acute toxicity after prostate cancer radiotherapy, they are mostly based on physician-reported measures of quality-of-life [13, 27,28,29,30,31,32,33,34,35,36,37,38]. Additionally, not all patients treated in these studies were treated using modern treatment techniques. All the patients included in this study were treated with IMRT receiving a median dose of 66.6 Gy. They also underwent regular image guidance. This represents modern-era treatment techniques for this cohort, and should therefore provide an accurate picture of acute GU toxicity in current practice.

The model parameters based on the LENT/SOMA scoring system (≥ Grade 2) were D50 = 69.56 Gy, γ = 1.7 and s = 0.35 for the RS model and D50 = 78.68 Gy, m = 0.17 and n = 0.09 for the LKB model [39, 40]. Regarding the symptoms of bladder contracture and volume loss, the parameters of the RS model are D50 = 80.0 Gy, γ = 3.0 and s = 0.18, whereas those of the LKB model are D50 = 80.0 Gy, m = 0.11 and n = 0.5 [41,42,43]. All the aforementioned parameters were derived using physician-rated outcome scores. The last few years, there is a trend to evaluate radiation-induced toxicity using patient-reported outcome scores. For the symptom of urgency, the parameters for the LKB model are D50 = 150.0, m = 0.37 and n = 0.01 [44]. The values that were derived by the present study using a different outcome scoring system for the symptom of urgency during the period of RT indicate a higher radiosensitivity (lower D50). Regarding acute urinary symptoms, it has been recently reported that an impact of dosimetric parameters was found for most symptoms (frequency, intermittency, urgency and nocturia), confirming the existence of a dose–volume/surface effect for acute effects [45]. However, those data and similar ones from other studies have not been modelled yet. So, in the literature there is lack of model parameter values for acute urinary symptoms based on PRO-QOL scoring protocols.

As indicated by the results of the goodness-of-fit, the three models studied here fitted the clinical data with similar accuracy. However, it should be pointed out that an agreement between those models is observed when fitting organs showing a parallel behavior, whereas discrepancies between them are observed when fitting organs of serial behavior. In the literature, we could not identify any other study deriving parameters for multiple NTCP models for acute urinary symptoms. In order to clinically validate the derived parameters, an independent cohort of patients with similar clinical characteristics should be available. In this study, the correlation of those parameters with the outcome data was mainly performed to identify the differences between the models regarding the goodness-of-fit process.

It should be noticed that for the symptoms of urgency and incontinence, the dose/volume metrics that correlated with the outcome data and the NTCP models had the volume expressed in percentage of the whole of volume of bladder, whereas frequency and nocturia had the best correlation when the volume was expressed in absolute units (cc). This means that for the symptoms of frequency and nocturia it is not the relative size of bladder receiving a given dose that is associated with them but the absolute volume of bladder. It is also clinically relevant that we did not find a clear association between bladder dose and urinary flow or dysuria. This is consistent with our clinical hypotheses. Radiation therapy can cause slower urinary flow due to swelling of the prostate which obstructs the urethra; but in patients who have previously undergone a radical prostatectomy, this mechanism does not apply. Further, clinically, dysuria is commonly deemed to be related to radiation irritation of the urethra, which was not assessed in this study because a Foley catheter is not routinely used in our treatment planning process. Therefore, the null findings in our analysis related to urine flow and dysuria provide further face validity regarding the results.

The current study does not account for other clinical factors that may have been associated with acute urinary toxicity, such as the use of medications (hormonal therapy), Prostate-specific antigen (PSA), Gleason score and volume of the gross tumor volume (GTV). However, few of these factors have consistently been found to impact acute toxicity [1, 46]. Most of those factors are not directly related to the function of bladder but mostly affect the treatment plan and consequently the dose distribution delivered to the patient. So, analyzing the dose to bladder against those symptoms should lead to better correlations. Moreover, given the small sample size of this study additional subset analyses could not performed.

A few points of caution should be mentioned regarding this study. First, it is difficult to compare the results of this study with previously reported findings due to the fact that a new patient reported outcome system (PCSI) was used for the acute urinary symptoms instead of the most common CTCAE scoring system. Second, although patient toxicity was based on the validated QOL instrument specifications, the definition of a significant toxicity was somewhat subjective (a two-point increase in the symptom score on the four-point symptom scale). However, other groups using patient reported outcomes have followed this approach because it acknowledges the importance of taking into account the baseline status [47,48,49]. Third, we had limited possibility of performing sub-sample test validation or k-fold validation of our models due to the size of our dataset. A future goal should be to collect new clinical data to separately validate the predictive ability of the models and their suitability to be used as constraints in treatment planning.


In this study it was shown that all the examined NTCP models (LKB, Logit and Relative Seriality) could fit the individual patient reported outcome data with very similar accuracy. The NTCP predictions of frequency correlated a little better with the outcome data than the rest of the acute symptoms. The values of the model parameters for the different acute symptoms could not be compared with previously published values, which is most likely attributed to the patient reported outcome system used instead of the CTCAE system. Further investigation with a larger patient cohort could verify the suitability of using the values of the NTCP model as an additional constraint in IMRT treatment plan optimization.


\( \overline{\overline{D}} \) or BEUD:

biologically effective uniform dose


Akaike information criterion


area under the curve


Common Terminology Criteria for Adverse Events


dose volume histogram

EDQ2Gy :

equivalent doses of 2 Gy per fraction


generalized Equivalent Uniform Dose


gross tumor volume




intensity modulated radiotherapy




normal tissue complication probability


Odds Ratio


Prostate Cancer Symptom Indices


patient-reported outcomes


Prostate-specific antigen


planning target volume


quality of life


receiver operating characteristic


Relative Seriality


radiation therapy


volumetric modulated arc therapy


  1. Fiorino C, Valdagni R, Rancati T, et al. Dose–volume effects for normal tissues in external radiotherapy: pelvis. Radiother Oncol. 2009;93:153–67.

    Article  PubMed  Google Scholar 

  2. Budaus L, Bolla M, Bossi A, et al. Functional outcome and complications following radiation therapy for prostate cancer: a critical analysis of the literature. Eur Urol. 2012;61:112–27.

    Article  PubMed  Google Scholar 

  3. Viswanathan AN, Yorke ED, Marks LB, et al. Radiation dose volume effects of the urinary bladder. Int J Radiat Oncol Biol. 2010;76:S123–9.

    Article  Google Scholar 

  4. Roach M III, Nam J, Gagliardi G, et al. Radiation dose–volume effects of the penile bulb. Int J Radiat Oncol Biol Phys. 2010;76:S130–4.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Ferlay J, Parkin DM, Steliarova-Foucher E. Estimates of cancer incidence and mortality in Europe in 2008. Eur J Cancer. 2010;46:765–81.

    CAS  Article  PubMed  Google Scholar 

  6. Clark JA, Talcott JA. Symptom indexes to assess outcomes of treatment for early prostate cancer. Med Care. 2001;39:1118–30.

    CAS  Article  PubMed  Google Scholar 

  7. Chen RC, Basak R, Meyer AM, et al. Association between choice of radical prostatectomy, external beam radiotherapy, brachytherapy, or active surveillance and patient-reported quality of life among men with localized prostate cancer. JAMA. 2017;317:1141–50.

    Article  PubMed  Google Scholar 

  8. Wagner LI, Wenzel L, Shaw E, et al. Patient-reported outcomes in phase II cancer clinical trials: lessons learned and future directions. J Clin Oncol. 2007;25:5058–62.

    Article  PubMed  Google Scholar 

  9. Chen RC, Mamon HJ, Chen YH, et al. Patient-reported acute gastrointestinal symptoms during concurrent chemoradiation treatment for rectal cancer. Cancer. 2010;116:1879–86.

    Article  PubMed  Google Scholar 

  10. Flores LT, Bennett AV, Law EB, et al. Patient-reported outcomes vs. clinician symptom reporting during Chemoradiation for rectal cancer. Gastrointest Cancer Res. 2012;5:119–24.

    PubMed  PubMed Central  Google Scholar 

  11. Chen RC, Zhang Y, Chen MH, et al. Patient-reported quality of life during radiation treatment for localized prostate cancer: results from a prospective phase II trial. BJU Int. 2012;110:1690–5.

    Article  PubMed  Google Scholar 

  12. Talcott J, Clark J, Manola J, et al. Bringing prostate cancer quality of life research back to the bedside: translating numbers into a format that patients can understand. J Urol. 2006;176:1558–63.

    Article  PubMed  Google Scholar 

  13. Matzinger O, Duclos F, van den Bergh A, et al. Acute toxicity of curative radiotherapy for intermediate- and high-risk localised prostate cancer in the EORTC trial 22991. Eur J Cancer. 2009;45:2825–34.

    CAS  Article  PubMed  Google Scholar 

  14. Mavroidis P, Price A, Fried D, et al. Dose-volume toxicity modeling for de-intensified chemo-radiation therapy for HPV-positive oropharynx cancer. Radiother Oncol. 2017;124:240–7.

    Article  PubMed  Google Scholar 

  15. Fowler JF. Brief summary of radiobiological principles in fractionated radiotherapy. Semin Radiat Oncol. 1992;2:16–21.

    Article  Google Scholar 

  16. Fowler JF. Sensitivity analysis of parameters in linear-quadratic radiobiologic modeling. Int J Radiat Oncol Biol Phys. 2009;73:1532–7.

    Article  PubMed  Google Scholar 

  17. Niemierko A. A generalized concept of equivalent uniform dose. Med Phys. 1999;26:1100.

    Google Scholar 

  18. Kwa S, Lebesque J, Theuws JC, et al. Radiation pneumonitis as a function of mean lung dose: an analysis of pooled data of 540 patients. Int J Radiat Oncol Biol Phys. 1998;42:1–9.

    CAS  Article  PubMed  Google Scholar 

  19. Seppenwoolde Y, Lebesque J, Jaeger K, et al. Comparing different NTCP models that predict the incidence of radiation pneumonitis. Int J Radiat Oncol Biol Phys. 2003;55:724–35.

    Article  PubMed  Google Scholar 

  20. Jackson A, Ten Haken RK, Robertson JM, et al. Analysis of clinical complication data for radiation hepatitis using a parallel architecture model. Int J Radiat Oncol Biol Phys. 1995;31:883–91.

    CAS  Article  PubMed  Google Scholar 

  21. Lind BK, Mavroidis P, Hyödynmaa S, et al. Optimization of the dose level for a given treatment plan to maximize the complication free tumor cure. Acta Oncol. 1999;38:787–98.

    CAS  Article  PubMed  Google Scholar 

  22. Källman P, Ågren AK, Brahme A. Tumor and normal tissue responses to fractionated non uniform dose delivery. Int J Radiat Biol. 1992;62:249–62.

    Article  PubMed  Google Scholar 

  23. Mavroidis P, Laurell G, Kraepelien T, et al. Determination and clinical verification of dose-response parameters for esophageal stricture from head and neck radiotherapy. Acta Oncol. 2003;42:865–81.

    Article  PubMed  Google Scholar 

  24. Herring DF. Methods for extracting dose-response curves from radiation therapy data, I: a unified approach. Int J Radiat Oncol Biol Phys. 1980;6:225–32.

    CAS  Article  PubMed  Google Scholar 

  25. Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.

    Article  Google Scholar 

  26. Wang K, Eblan MJ, Deal AM, et al. Cardiac toxicity after radiotherapy for stage III non-small-cell lung cancer: pooled analysis of dose-escalation trials delivering 70 to 90 Gy. J Clin Oncol. 2017;35:1387–94.

    Article  PubMed  Google Scholar 

  27. De Meerleer G, Vakaet L, Meersschout S, et al. Intensity-modulated radiotherapy as primary treatment for prostate cancer: acute toxicity in 114 patients. Int J Radiat Oncol Biol Phys. 2004;60:777–87.

    Article  PubMed  Google Scholar 

  28. Guerrero Urbano T, Khoo V, Staffurth J, et al. Intensity-modulated radiotherapy allows escalation of the radiation dose to the pelvic lymph nodes in patients with locally advanced prostate cancer: preliminary results of a phase I dose escalation study. Clin Oncol. 2010;22:236–44.

    CAS  Article  Google Scholar 

  29. Muller AC, Lutjens J, Alber M, et al. Toxicity and outcome of pelvic IMRT for node-positive prostate cancer. Strahlenther Onkol. 2012;188:982–9.

    Article  PubMed  Google Scholar 

  30. Goldner G, Wachter-Gerstner N, Wachter S, et al. Acute side effects during 3-D-planned conformal radiotherapy of prostate cancer. Differences between patient's self-reported questionnaire and the corresponding doctor's report. Strahlenther Onkol. 2003;179:320–7.

    Article  PubMed  Google Scholar 

  31. Stephens RJ, Hopwood P, Girling DJ, et al. Randomized trials with quality of life endpoints: are doctors' ratings of patients' physical symptoms interchangeable with patients' self-ratings? Qual Life Res. 1997;6:225–36.

    CAS  Article  PubMed  Google Scholar 

  32. Clauser SB, Ganz PA, Lipscomb J, et al. Patient-reported outcomes assessment in cancer trials: evaluating and enhancing the payoff to decision making. J Clin Oncol. 2007;25:5049–50.

    Article  PubMed  Google Scholar 

  33. Carillo V, Cozzarini C, Rancati T, et al. Relationships between bladder dose–volume/surface histograms and acute urinary toxicity after radiotherapy for prostate cancer. Radiother Oncol. 2014;111:100–5.

    Article  PubMed  Google Scholar 

  34. Chua B, Min M, Wood M, et al. Implementation of an image guided intensity-modulated protocol for post-prostatectomy radiotherapy: planning data and acute toxicity outcomes. J Med Imaging Radiat Oncol. 2013;57:482–9.

    Article  PubMed  Google Scholar 

  35. Myrehaug S, Chan G, Craig T, et al. A treatment planning and acute toxicity comparison of two pelvic nodal volume delineation techniques and delivery comparison of intensity-modulated radiotherapy versus volumetric modulated arc therapy for hypofractionated high-risk prostate cancer radiotherapy. Int J Radiat Oncol Biol Phys. 2012;82:e657–62.

    Article  PubMed  Google Scholar 

  36. Vesprini D, Catton C, Jacks L, et al. Inverse relationship between biochemical outcome and acute toxicity after image-guided radiotherapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2012;83:608–16.

    Article  PubMed  Google Scholar 

  37. Gill S, Thomas J, Fox C, et al. Acute toxicity in prostate cancer patients treated with and without image-guided radiotherapy. Radiat Oncol. 2011;6:145.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Tsai HK, Manola J, Abner A, et al. Patient-reported acute gastrointestinal toxicity in men receiving 3-dimensional conformal radiation therapy for prostate cancer with or without neoadjuvant androgen suppression therapy. Urol Oncol. 2005;23:230–7.

    CAS  Article  PubMed  Google Scholar 

  39. Zhu J, Simon A, Haigron P, et al. The benefit of using bladder sub-volume equivalent uniform dose constraints in prostate intensity-modulated radiotherapy planning. Onco Targets Ther. 2016;9:7537–44.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Mavroidis P, Komisopoulos G, Buckey C, et al. Radiobiological evaluation of prostate cancer IMRT and conformal-RT plans using different treatment protocols. Phys Med. 2017;40:33–41.

    Article  PubMed  Google Scholar 

  41. Ågren AK: Quantification of the response of heterogeneous tumors and organized normal tissues to fractionated radiotherapy. Ph.D. Thesis. Stockholm: Stockholm University, 1995.

  42. Burman C, Kutcher GJ, Emami B, et al. Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys. 1991;21:123–35.

    CAS  Article  PubMed  Google Scholar 

  43. Emami B, Lyman J, Brown A, et al. Tolerance of normal tissue to therapeutic irradiation. Int J Radiat Oncol Biol Phys. 1991;21:109–22.

    CAS  Article  PubMed  Google Scholar 

  44. Thor M, Olsson C, Oh JH, et al. Urinary bladder dose-response relationships for patient-reported genitourinary morbidity domains following prostate cancer radiotherapy. Radiother Oncol. 2016;119:117–22.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Cozzarini C, Rancati T, Carillo V, et al. Multi-variable models predicting specific patient-reported acute urinary symptoms after radiotherapy for prostate cancer: results of a cohort study. Radiother Oncol. 2015;116:185–91.

    Article  PubMed  Google Scholar 

  46. Palorini F, Rancati T, Cozzarini C, et al. Multi-variable models of large international prostate symptom score worsening at the end of therapy in prostate cancer radiotherapy. Radiother Oncol. 2016;118:92–8.

    Article  PubMed  Google Scholar 

  47. Scher HI, Morris MJ, Stadler WM, et al. Trial design and objectives for castration-resistant prostate cancer: updated recommendations from the prostate cancer clinical trials working group 3. J Clin Oncol. 2016;34:1402–18.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Beetz I, Schilstra C, Van der Schaaf A, et al. NTCP models for patient-rated xerostomia and sticky saliva after treatment with intensity modulated radiotherapy for head and neck cancer: the role of dosimetric and clinical factors. Radiother Oncol. 2012;105:101–6.

    Article  PubMed  Google Scholar 

  49. Beetz I, Schilstra C, Van Luijk P, et al. External validation of three dimensional conformal radiotherapy based NTCP models for patient-rated xerostomia and sticky saliva among patients treated with intensity modulated radiotherapy. Radiother Oncol. 2012;105:94–100.

    Article  PubMed  Google Scholar 

Download references


Not applicable


Not applicable

Availability of data and materials

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

Author information

Authors and Affiliations



PM and RCC conceived the idea and designed the study. KAP, AZW, JS and SS collected and reviewed patient data. PM, JD and SKD performed the statistical modeling and analyzed the data. All authors participated in drafting and revising the manuscript. All authors have given their final approval of the manuscript.

Corresponding author

Correspondence to Panayiotis Mavroidis.

Ethics declarations

Ethics approval and consent to participate

There is an IRB approval for this study by the Ethics committee of the University of North Carolina. Additionally, this is retrospective analysis so no patient consent required.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mavroidis, P., Pearlstein, K.A., Dooley, J. et al. Fitting NTCP models to bladder doses and acute urinary symptoms during post-prostatectomy radiotherapy. Radiat Oncol 13, 17 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Radiobiological parameters
  • Acute urinary symptoms
  • Post-prostatectomy radiotherapy
  • NTCP
  • LKB
  • Logit
  • Relative seriality