DR.Warkentin: Function that calculates TCP according Warkentin (Poisson)...
In robertogattabs/RadAgent: Package for handling DICOM RT files

Description Usage Arguments Details Value References

View source: R/Dose.Response.R View source: R/modeling.R

This function calculates the Tumor Control Probability according the Warkentin model.

1	DR.Warkentin(doses, TD50 = 45, gamma50 = 1.5, a = 1)

`doses`	Either a `dvhmatrix` class object or a vector with nominal doses
`TD50`	The value of dose that gives the 50% of probability of outcome
`gamma50`	The slope of dose/response curve at 50% of probability
`a`	Value for parallel-serial correlation in radiobiological response

This model is the equivalent of the Poisson model where the covariates and their coefficients have been reported as function of TD_{50} and γ_{50}. In Warkentin paper model computation starts with formula:

TCP=e^{-Np_{s}(D)}

where N is the initial number of clonogens, p_{s}(D) is the survival fraction after the dose D. The previous equation can be rewritten as function of TD_{50} and γ_{50}:

TCP=≤ft (\frac{1}{2} \right )^{e^{≤ft [2γ_{50}(1-D/D_{50})/ln2 \right ]}}

Using the assumption of independent subvolumes, for the case of heterogeneous irradiation, the overall probability of tumor control is the product of the probabilities of killing all clonogens in each tumor subvolume described by the differential DVH:

TCP=∏_{i}TCP(D_{i},v_{i})

Thus, for a given differential DVH ≤ft \{ D_{i},v_{i} \right \}, the TCP can be calculated using the following two-parameter TCP formula:

TCP=≤ft (\frac{1}{2} \right )^{∑_{i}v_{i}e^{≤ft [2γ_{50}(1-D/D_{50})/ln2 \right ]}}

In the model equation D can be either the nominal dose or the EUD as calculated by DVH.eud function.

A vector with TCP calculated according Warkentin (Poisson) model.

Warkentin B, Stavrev P, Stavreva N, Field C, Fallone BG. A TCP-NTCP estimation module using DVHs and known radiobiological models and parameter sets. J Appl Clin Med Phys. 2004 Winter;5(1):50-63. Epub 2004 Jan 1. PubMed PMID: 15753933.

robertogattabs/RadAgent documentation built on June 30, 2018, 12:02 a.m.