Implementation of a graphical device developed by Rex Galbraith to display several estimates of the same quantity that have different standard errors.
radialplot(x, ...)# S3 method for default
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", sigdig = 2, show.numbers = FALSE, pch = 21,
levels = NA, bg = c("white", "red"), title = TRUE, k = 0,
markers = NULL, ...)
# S3 method for fissiontracks
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "arcsin", sigdig = 2, show.numbers = FALSE, pch = 21,
levels = NA, bg = c("white", "red"), title = TRUE, markers = NULL,
k = 0, exterr = TRUE, ...)
# S3 method for UPb
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, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
common.Pb = 0, ...)
# S3 method for PbPb
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
i2i = TRUE, ...)
# S3 method for ArAr
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
i2i = FALSE, ...)
# S3 method for UThHe
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, ...)
# S3 method for ReOs
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
i2i = TRUE, ...)
# S3 method for SmNd
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
i2i = TRUE, ...)
# S3 method for RbSr
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
i2i = TRUE, ...)
# S3 method for LuHf
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, exterr = TRUE,
i2i = TRUE, ...)
# S3 method for ThU
radialplot(x, from = NA, to = NA, t0 = NA,
transformation = "log", show.numbers = FALSE, pch = 21, levels = NA,
bg = c("white", "red"), markers = NULL, k = 0, i2i = TRUE, ...)
Either an nx2 matix of (transformed) values z and their standard errors s
OR
and object of class fissiontracks, UThHe,
ArAr, ReOs, SmNd, RbSr, LuHf,
ThU, PbPb or UPb
additional arguments to the generic points
function
minimum age limit of the radial scale
maximum age limit of the radial scale
central value
one of either log, linear or
(if x has class fissiontracks)
the number of significant digits of the numerical values reported in the title of the graphical output.
boolean flag (TRUE to show grain
numbers)
plot character (default is a filled circle)
a vector with additional values to be displayed as different background colours of the plot symbols.
a vector of two background colours for the plot symbols.
If levels=NA, then only the first colour will be
used. If levels is a vector of numbers, then bg
is used to construct a colour ramp.
add a title to the plot?
number of peaks to fit using the finite mixture models of
Galbraith and Green (1993). Setting k='auto'
automatically selects an optimal number of components based on
the Bayes Information Criterion (BIC). Setting k='min'
estimates the minimum value using a three parameter model
consisting of a Normal distribution truncated by a discrete
component.
vector of ages of radial marker lines to add to the plot.
propagate the external sources of uncertainty into the mixture model errors?
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)
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.
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.
apply a common lead correction using one of three methods:
1: use the isochron intercept as the initial Pb-composition
2: use the Stacey-Kramer two-stage model to infer the initial
Pb-composition
3: use the Pb-composition stored in
settings('iratio','Pb206Pb204') and
settings('iratio','Pb207Pb204')
`isochron to intercept': calculates the initial
(aka `inherited', `excess', or `common') \(^{40}\)Ar/\(^{36}\)Ar,
\(^{207}\)Pb/\(^{204}\)Pb, \(^{87}\)Sr/\(^{86}\)Sr,
\(^{143}\)Nd/\(^{144}\)Nd, \(^{187}\)Os/\(^{188}\)Os or
\(^{176}\)Hf/\(^{177}\)Hf ratio from an isochron
fit. Setting i2i to FALSE uses the default values
stored in settings('iratio',...) or zero (for the Pb-Pb
method). When applied to data of class ThU, setting
i2i to TRUE applies a detrital Th-correction.
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.
# NOT RUN {
data(examples)
radialplot(examples$FT1)
# }
Run the code above in your browser using DataLab