Luminescence (version 0.2.1)

plot_RadialPlot: Produce a Galbraith's radial plot.

Description

A Galbraith radial plot is produced on a logarithmic or a linear De scale. The function is based on a rewritten S script of Rex Galbraith.

Usage

plot_RadialPlot(sample, 
                sample.groups, 
                sample.legend, 
                sample.lty = 1, 
                sample.pch = 1, 
                sample.col = "black", 
                sample.mtext = "default", 
                zscale.log = TRUE, 
                zaxis.scale, zaxis.group_circle = FALSE, yaxis.scale, 
                plot.2sigmaRange = TRUE, 
                plot.area_ratio = 4.5/6, 
                zlab = expression(paste(D[e], " [Gy]")), 
                main = expression(paste(D[e], " Distribution", sep = "")), 
                cex.global = 1, xscale_factor = 1.01)

Arguments

sample

data.frame (required): two column data frame with the input values, e.g. 'de' and 'se'.

sample.groups

list (optional): option to group the input data set like: sample.groups = list(c(1:14), c(15:26), c(27:40))

sample.legend

character (optional): character vector for a legend. This option is provided for the parameter sample.groups, but can also be used for one sample.

sample.lty

vector (with default): line type for the central value (see par. If the sample is grouped, a line type can be defined for each data set. )

sample.pch

vector (with default): point type for the presented data (see par. If the sample is grouped, a point type can be defined for each data.set

sample.col

vector or character (with default): colour of the points and lines (see colors). If the sample is grouped, a colour can be defined for each group.

sample.mtext

character (optional): mtext on the top of the plot. This option is only availabe if the grouping option is used. Otherwise information on the distribution is shown.

zscale.log

logical (with default): log De scale (TRUE/FALSE)

zaxis.scale

numeric (optional): option to set the z-scale manually. Example: zaxis.scale = seq(50,120, by = 10)

zaxis.group_circle

logical (with default): shows additional group circles for the 2-sigma uncertainties on the z-scale.

yaxis.scale

vector (optional): option for manual y-axis scaling. Example: yaxis.scale=c(-15,15)

plot.2sigmaRange

logical (with default): plot a grey polygon showing the 2-sigma range of the central value.

plot.area_ratio

vector (with default): option for manual plot ratio adjustments. Example plot.area_ratio=4.5/5.5

zlab

character (with default): z-axis (semi circle) label

main

character (with default): title of the plot

cex.global

numeric (with default): global scaling factor

xscale_factor

numeric (with default): scaling factor for the z-scale ticks.

Value

• Returns a plot.

Details

Details and the theoretical background on the radial plot are given in the cited literature.

References

Galbraith, R.F., 1988. Graphical Display of Estimates Having Differing Standard Errors. Technometrics, 30 (3), pp. 271-281.

Galbraith, R.F., 1990. The radial plot: Graphical assessment of spread in ages. International Journal of Radiation Applications and Instrumentation. Part D. Nuclear Tracks and Radiation Measurements, 17 (3), pp. 207-214.

Galbraith, R. & Green, P., 1990. Estimating the component ages in a finite mixture. International Journal of Radiation Applications and Instrumentation. Part D. Nuclear Tracks and Radiation Measurements, 17 (3), pp. 197-206.

Galbraith, R.F. & Laslett, G.M., 1993. Statistical models for mixed fission track ages. Nuclear Tracks And Radiation Measurements, 21 (4), pp. 459-470.

Galbraith, R.F., 1994. Some Applications of Radial Plots. Journal of the American Statistical Association, 89 (428), pp. 1232-1242.

Galbraith, R.F., 2010. On plotting OSL equivalent doses. Ancient TL, 28 (1), pp. 1-10.

Galbraith, R.F. & Roberts, R.G., 2012. Statistical aspects of equivalent dose and error calculation and display in OSL dating: An overview and some recommendations. Quaternary Geochronology, 11, pp.1-27.

See Also

plot, legend

Examples

##load data
data(ExampleData.DeValues)

##plot
plot_RadialPlot(ExampleData.DeValues,zscale.log = TRUE, zaxis.scale = seq(1995,4250, by = 500), 
                zlab = expression(paste(D[e], "[s]")))