IsoplotR (version 0.8)

radialplot: Visualise heteroscedastic data on a radial plot

Description

Implementation of a graphical device developed by Rex Galbraith to display several estimates of the same quantity that have different standard errors.

Usage

radialplot(x, ...)
"radialplot"(x, from = NA, to = NA, t0 = NA, transformation = "log", sigdig = 2, show.numbers = FALSE, pch = 21, bg = "white", k = 0, markers = NULL, ...)
"radialplot"(x, from = NA, to = NA, t0 = NA, transformation = "arcsin", sigdig = 2, show.numbers = FALSE, pch = 21, bg = "white", title = TRUE, markers = NULL, k = 0, exterr = TRUE, ...)
"radialplot"(x, from = NA, to = NA, t0 = NA, transformation = "log", type = 4, cutoff.76 = 1100, cutoff.disc = c(-15, 5), show.numbers = FALSE, pch = 21, bg = "white", markers = NULL, k = 0, exterr = TRUE, ...)
"radialplot"(x, from = NA, to = NA, t0 = NA, transformation = "log", show.numbers = FALSE, pch = 21, bg = "white", markers = NULL, k = 0, exterr = TRUE, ...)
"radialplot"(x, from = NA, to = NA, t0 = NA, transformation = "log", show.numbers = FALSE, pch = 21, bg = "white", markers = NULL, k = 0, ...)

Arguments

x

Either an nx2 matix of (transformed) values z and their standard errors s

OR

and object of class fissiontracks, UThHe, ArAr, or UPb

...

additional arguments to the generic points function

from

minimum age limit of the radial scale

to

maximum age limit of the radial scale

t0

central value

transformation

one of either log, linear or (if x has class fissiontracks)

sigdig

the number of significant digits of the numerical values reported in the title of the graphical output.

show.numbers

boolean flag (TRUE to show grain numbers)

pch

plot character (default is a filled circle)

bg

background colour of the plot character

k

number of mixture model components to fit.

markers

vector of ages of radial marker lines to add to the plot.

title

add a title to the plot?

exterr

propagate the external sources of uncertainty into the mixture model errors?

type

scalar indicating whether to plot the $^{207}$Pb/$^{235}$U age (type=1), the $^{206}$Pb/$^{238}$U age (type=2), the $^{207}$Pb/$^{206}$Pb age (type=3), the $^{207}$Pb/$^{206}$Pb-$^{206}$Pb/$^{238}$U age (type=4), or the (Wetherill) concordia age (type=5)

cutoff.76

the age (in Ma) below which the $^{206}$Pb/$^{238}$U and above which the $^{207}$Pb/$^{206}$Pb age is used. This parameter is only used if type=4.

cutoff.disc

two element vector with the maximum and minimum percentage discordance allowed between the $^{207}$Pb/$^{235}$U and $^{206}$Pb/$^{238}$U age (if $^{206}$Pb/$^{238}$U < cutoff.76) or between the $^{206}$Pb/$^{238}$U and $^{207}$Pb/$^{206}$Pb age (if $^{206}$Pb/$^{238}$U > cutoff.76). Set cutoff.disc=NA if you do not want to use this filter.

References

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.

Examples

data(examples)
radialplot(examples$FT1)
#fname <- system.file("FT1.csv",package="IsoplotR")
#FT <- read.data(fname,method="fissiontracks",format=1)