-
Notifications
You must be signed in to change notification settings - Fork 8
/
Analyse_SAR.OSLdata.Rd
161 lines (135 loc) · 6.25 KB
/
Analyse_SAR.OSLdata.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/Analyse_SAR.OSLdata.R
\name{Analyse_SAR.OSLdata}
\alias{Analyse_SAR.OSLdata}
\title{Analyse SAR CW-OSL measurements (deprecated)}
\usage{
Analyse_SAR.OSLdata(
input.data,
signal.integral,
background.integral,
position,
run,
set,
dtype,
keep.SEL = FALSE,
info.measurement = "unknown measurement",
output.plot = FALSE,
plot_singlePanels = FALSE,
cex.global = 1,
...
)
}
\arguments{
\item{input.data}{\linkS4class{Risoe.BINfileData} (\strong{required}):
input data from a Risø BIN file, produced by the function \link{read_BIN2R}.}
\item{signal.integral}{\link{vector} (\strong{required}):
channels used for the signal integral, e.g. \code{signal.integral=c(1:2)}}
\item{background.integral}{\link{vector} (\strong{required}):
channels used for the background integral, e.g. \code{background.integral=c(85:100)}}
\item{position}{\link{vector} (\emph{optional}):
reader positions that want to be analysed (e.g. \code{position=c(1:48)}.
Empty positions are automatically omitted. If no value is given all
positions are analysed by default.}
\item{run}{\link{vector} (\emph{optional}):
range of runs used for the analysis. If no value is given the range of the
runs in the sequence is deduced from the \code{Risoe.BINfileData} object.}
\item{set}{\link{vector} (\emph{optional}):
range of sets used for the analysis. If no value is given the range of the
sets in the sequence is deduced from the \code{Risoe.BINfileData} object.}
\item{dtype}{\link{character} (\emph{optional}):
allows to further limit the curves by their data type (\code{DTYPE}),
e.g., \code{dtype = c("Natural", "Dose")} limits the curves to this two data types.
By default all values are allowed.
See \linkS4class{Risoe.BINfileData} for allowed data types.}
\item{keep.SEL}{\link{logical} (default):
option allowing to use the \code{SEL} element of the \linkS4class{Risoe.BINfileData} manually.
\strong{NOTE:} In this case any limitation provided by \code{run}, \code{set} and \code{dtype}
are ignored!}
\item{info.measurement}{\link{character} (\emph{with default}):
option to provide information about the measurement on the plot
output (e.g. name of the BIN or BINX file).}
\item{output.plot}{\link{logical} (\emph{with default}):
plot output (\code{TRUE/FALSE})}
\item{plot_singlePanels}{\link{logical} (\emph{with default}):
single plot output (\code{TRUE/FALSE}) to allow for plotting the results in
single plot windows. Requires \code{output.plot = TRUE}.}
\item{cex.global}{\link{numeric} (\emph{with default}):
global scaling factor.}
\item{...}{further arguments that will be passed to the function
\link{calc_OSLLxTxRatio} (supported: \code{background.count.distribution}, \code{sigmab},
\code{sig0}; e.g., for instrumental error) and can be used to adjust the plot.
Supported" \code{mtext}, \code{log}}
}
\value{
A plot (\emph{optional}) and \link{list} is returned containing the
following elements:
\item{LnLxTnTx}{\link{data.frame} of all calculated Lx/Tx values including signal, background counts and the dose points.}
\item{RejectionCriteria}{\link{data.frame} with values that might by used as rejection criteria. NA is produced if no R0 dose point exists.}
\item{SARParameters}{\link{data.frame} of additional measurement parameters obtained from the BIN file, e.g. preheat or read temperature
(not valid for all types of measurements).}
}
\description{
The function analyses SAR CW-OSL curve data and provides a summary of the
measured data for every position. The output of the function is optimised
for SAR OSL measurements on quartz.
The function works only for standard SAR protocol measurements introduced by
Murray and Wintle (2000) with CW-OSL curves. For the calculation of the
Lx/Tx value the function \link{calc_OSLLxTxRatio} is used.
\strong{Provided rejection criteria}
\verb{[recyling ratio]}: calculated for every repeated regeneration dose point.
\verb{[recuperation]}: recuperation rate calculated by comparing the \code{Lx/Tx} values of the zero
regeneration point with the \code{Ln/Tn} value (the \code{Lx/Tx} ratio of the natural
signal). For methodological background see Aitken and Smith (1988)
\verb{[IRSL/BOSL]}: the integrated counts (\code{signal.integral}) of an
IRSL curve are compared to the integrated counts of the first regenerated
dose point. It is assumed that IRSL curves got the same dose as the first
regenerated dose point. \strong{Note:} This is not the IR depletion ratio
described by Duller (2003).
}
\note{
Rejection criteria are calculated but not considered during the
analysis to discard values.
\strong{The analysis of IRSL data is not directly supported}. You may want to
consider using the functions \link{analyse_SAR.CWOSL} or
\link{analyse_pIRIRSequence} instead.
\strong{The development of this function will not be continued. We recommend to use the function \link{analyse_SAR.CWOSL} or instead.}
}
\section{Function version}{
0.2.18
}
\examples{
##load data
data(ExampleData.BINfileData, envir = environment())
##analyse data
suppressWarnings( # silence the deprecation warning
output <- Analyse_SAR.OSLdata(input.data = CWOSL.SAR.Data,
signal.integral = c(1:5),
background.integral = c(900:1000),
position = c(1:1),
output.plot = TRUE)
)
##combine results relevant for further analysis
output.SAR <- data.frame(Dose = output$LnLxTnTx[[1]]$Dose,
LxTx = output$LnLxTnTx[[1]]$LxTx,
LxTx.Error = output$LnLxTnTx[[1]]$LxTx.Error)
output.SAR
}
\references{
Aitken, M.J. and Smith, B.W., 1988. Optical dating: recuperation
after bleaching. Quaternary Science Reviews 7, 387-393.
Duller, G., 2003. Distinguishing quartz and feldspar in single grain
luminescence measurements. Radiation Measurements, 37 (2), 161-165.
Murray, A.S. and Wintle, A.G., 2000. Luminescence dating of quartz using an
improved single-aliquot regenerative-dose protocol. Radiation Measurements
32, 57-73.
}
\seealso{
\link{calc_OSLLxTxRatio}, \linkS4class{Risoe.BINfileData}, \link{read_BIN2R}, \link{fit_DoseResponseCurve}
}
\author{
Sebastian Kreutzer, Institute of Geography, Heidelberg University (Germany)\cr
Margret C. Fuchs, HZDR, Freiberg (Germany)
}
\keyword{datagen}
\keyword{dplot}