plot_GrowthCurve: Fit and plot a dose-response curve for luminescence data...

In Luminescence: Comprehensive Luminescence Dating Data Analysis

Fit and plot a dose-response curve for luminescence data (Lx/Tx against dose)

Description

A dose-response curve is produced for luminescence measurements using a regenerative or additive protocol. This is a wrapper around fit_DoseResponseCurve and plot_DoseResponseCurve.

Usage

plot_GrowthCurve(
  object,
  mode = "interpolation",
  fit.method = "SSE",
  output.plot = TRUE,
  output.plotExtended = TRUE,
  plot_singlePanels = FALSE,
  verbose = TRUE,
  n.MC = 100,
  ...
)

Arguments

object	data.frame or a list of such objects (required): data frame with columns for Dose, LxTx, LxTx.Error and TnTx (optional). If these column names are used, then they can be passed in whatever order; otherwise columns are taken by position. If object is a list, the function is called on each of its elements. If fit.method = "OTORX" you have to provide the test dose in the same unit as the dose in a column called Test_Dose. The function searches explicitly for this column name. Only the first value will be used assuming a constant test dose over the measurement cycle.
mode	character (with default): selects calculation mode of the function. "interpolation" (default) calculates the De by interpolation, "extrapolation" calculates the equivalent dose by extrapolation (useful for MAAD measurements) and "alternate" calculates no equivalent dose and just fits the data points. Please note that for option "interpolation" the first point is considered as natural dose.
fit.method	character (with default): function used for fitting. Possible options are: LIN, QDR, SSE, ⁠SSE OR LIN⁠, SSE+LIN, DSE (not defined for extrapolation), GOK, OTOR and OTORX. See details.
output.plot	logical (with default): enable/disable the plot output.
output.plotExtended	logical (with default): If TRUE (default), 3 plots on one plot area are provided: the dose-response curve, a histogram from Monte Carlo error simulation and a test dose response plot. If FALSE, just the growth curve will be plotted.
plot_singlePanels	logical (with default): single plot output (TRUE/FALSE) to allow for plotting the results in single plot windows. Requires output.plotExtended = TRUE.
verbose	logical (with default): enable/disable output to the terminal.
n.MC	integer (with default): number of Monte Carlo simulations for error estimation.
...	Further arguments to fit_DoseResponseCurve (fit.force_through_origin, fit_weights, fit.includingRepeatedRegPoints, fit.NumberRegPoints, fit.NumberRegPointsReal, fit_bounds, txtProgressBar) and to plot_DoseResponseCurve (xlim, ylim, main, mtext, xlab, ylab, log, legend (TRUE/FALSE), legend.pos, reg_points_pch, density_polygon (TRUE/FALSE), density_polygon_col, density_rug (TRUE/FALSE), lwd_drc, col_drc, lty_drc, box (TRUE/FALSE)).

Value

Along with a plot (if wanted) the RLum.Results object produced by fit_DoseResponseCurve is returned invisibly.

Function version

1.2.3

How to cite

Kreutzer, S., Dietze, M., Colombo, M., 2026. plot_GrowthCurve(): Fit and plot a dose-response curve for luminescence data (Lx/Tx against dose). Function version 1.2.3. In: Kreutzer, S., Burow, C., Dietze, M., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J., Mercier, N., Philippe, A., Riedesel, S., Autzen, M., Mittelstrass, D., Gray, H.J., Galharret, J., Colombo, M., Steinbuch, L., de Boer, A., Bluszcz, A., 2026. Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 1.3.1. https://r-lum.github.io/Luminescence/

Author(s)

Sebastian Kreutzer, F2.1 Geophysical Parametrisation/Regionalisation, LIAG - Institute for Applied Geophysics (Germany)
Michael Dietze, GFZ Potsdam (Germany)
Marco Colombo, Institute of Geography, Heidelberg University (Germany) , RLum Developer Team

References

Berger, G.W., Huntley, D.J., 1989. Test data for exponential fits. Ancient TL 7, 43-46.

Guralnik, B., Li, B., Jain, M., Chen, R., Paris, R.B., Murray, A.S., Li, S.-H., Pagonis, P., Herman, F., 2015. Radiation-induced growth and isothermal decay of infrared-stimulated luminescence from feldspar. Radiation Measurements 81, 224-231.

Pagonis, V., Kitis, G., Chen, R., 2020. A new analytical equation for the dose response of dosimetric materials, based on the Lambert W function. Journal of Luminescence 225, 117333. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1016/j.jlumin.2020.117333")}

See Also

fit_DoseResponseCurve, plot_DoseResponseCurve

Examples

##(1) plot growth curve for a dummy dataset
data(ExampleData.LxTxData, envir = environment())
plot_GrowthCurve(LxTxData)

##(1b) horizontal plot arrangement
layout(mat = matrix(c(1,1,2,3), ncol = 2))
plot_GrowthCurve(LxTxData, plot_singlePanels = TRUE)

##(2) plot the growth curve with pdf output - uncomment to use
##pdf(file = "~/Desktop/Growth_Curve_Dummy.pdf", paper = "special")
plot_GrowthCurve(LxTxData)
##dev.off()

##(3) plot the growth curve with pdf output - uncomment to use, single output
##pdf(file = "~/Desktop/Growth_Curve_Dummy.pdf", paper = "special")
temp <- plot_GrowthCurve(LxTxData, plot_singlePanels = TRUE)
##dev.off()

##(4) plot resulting function for given interval x
x <- seq(1,10000, by = 100)
plot(
 x = x,
 y = eval(temp$Formula),
 type = "l"
)

##(5) plot using the 'extrapolation' mode
LxTxData[1,2:3] <- c(0.5, 0.001)
print(plot_GrowthCurve(LxTxData, mode = "extrapolation"))

##(6) plot using the 'alternate' mode
LxTxData[1,2:3] <- c(0.5, 0.001)
print(plot_GrowthCurve(LxTxData, mode = "alternate"))

Luminescence documentation built on Sept. 18, 2026, 9:07 a.m.