alpha.beta.ion.range: Evaluation of ion alpha and beta LQ parameters

View source: R/biological-functions.R

alpha.beta.ion.rangeR Documentation

Evaluation of ion alpha and beta LQ parameters

Description

The evaluation is performed for a monoenergetic ion. It uses a C++ implementation for the evaluation. It accepts a single set of parameter depending on the choosen radiobiological model, or alternatively a full range of variability (min,max) for each parameter.

Usage

alpha.beta.ion.range(model = "MKM",
  model.parameters = data.frame(alpha0 = 0.1295, beta0 = 0.03085, rN = 4,
  rd = 0.31), cell.name = NULL, particle = "H", energies = NULL,
  lets = NULL, calculusType = "rapidMKM", precision = 0.5,
  ignore.stdout = TRUE, ignore.stderr = TRUE, get.raw.data = FALSE,
  remove.temp.files = TRUE)

Arguments

model

the name of the model (options: 'MKM', 'LEMI', 'LEMII', 'LEMIII')

model.parameters

a dataframe containing the "biological" parameter associated to a specific biological tissue/cell line. In the case of the MKM:

  • alpha0 LQ alpha parameter of the reference radiation [Gy^-1]

  • beta0 LQ beta parameter of the reference radiation [Gy^-2]

  • rN cell nucleus radius [um]

  • rd domain radius [um]

In the case of LEMI, LEMII, LEMIII:

  • alpha0 LQ alpha parameter of the reference radiation [Gy^-1]

  • beta0 LQ beta parameter of the reference radiation [Gy^-2]

  • rN cell nucleus radius [um]

  • Dt threshold dose [Gy]

Optionally, instead of the specific values, the range of the parameters can be specified in the same data.frame. In the case of MKM:

  • alpha0.min, alpha0.max, alpha0.N range of variability for parameter alpha0 and number of steps

  • beta0.min, beta0.max, beta0.N range of variability for parameter beta0 and number of steps

  • rN.min, rN.max, rN.N range of variability for parameter rN and number of steps

  • rd.min, rN.max domain radius [um]

particle

available particles: 'H', 'He', 'Li', 'Be', 'B,', 'C', 'N', 'O', 'F', 'Ne'.

energies

vector of energies for the particle [MeV]

lets

vector of LETs for the particle [keV/um]. It is used if energies is NULL.

calculusType

the type of the evaluations options are:

  • rapidKase implementation of MKM as described in (Kase2008)

  • rapidMKM more correct implementation for the MKM from the one described in (Kase2008)

  • newMKMKase "rapidKase" implementation + non poissionian correction for beta

  • newMKM "rapidMKM" implementation + non poissionian correction for beta

  • rapidScholz fast implementation of the LEMI, LEMII, LEMIII (Scholz2000)

  • rapidRusso fast implementation of the LEMI, LEMII, LEMIII, more accurate than the rapidScholz (Russo2010)

  • slow_alphaIon_betaIon slow Monte Carlo evaluation (compatible with MKM, LEMI, LEMII and LEMIII).

precision

used only for the Monte Carlo evaluation. If precision < 1.0 it represents the target relative standard deviation to be reached in the evaluations of the cell survivals. If precision >= 1 it represents the number of cells to be simulated.

get.raw.data

If calculusType=slow_alphaIon_betaIon, returns the raw data (survival vs dose for each simulated cell) from which alpha and beta parameters can be extracted, otherwhise it is ignored.

cellType

name of the tissue/cell line (optional)

Value

a data.frame containing all the information specified including the alpha and beta MKM evaluation (note, in the MKM implementation beta = betaX).

See Also

Other LEM/MKM Models: alpha.beta.mkm, alpha.fun.mkm


planit-group/Rplanit documentation built on Dec. 5, 2022, 11:10 p.m.