Use absorbance data to correct inner-filter effect in FEEM objects and collections of them.
feemife(x, …)
# S3 method for feem
feemife(x, absorbance, abs.path = 1, …)
# S3 method for feemcube
feemife(x, absorbance, abs.path, …)
# S3 method for list
feemife(x, absorbance, abs.path, …)A FEEM object, a FEEM data cube, or a list of them.
If x is a FEEM object: a two-column matrix-like object suitable
for xy.coords containing the absorbance spectrum of the
sample: the wavelengths in the first column and the unitless
absorbance values in the second column.
Otherwise, this could be a list of such objects or a multi-column
matrix-like objects. If x contains names of the samples
(is a named list or had names specified when calling
feemcube), absorbance is a named list or has
named columns, and all samples from x can be looked up in
absorbance, results of this lookup are used. If name lookup
fails but (given \(N\)-sample x) absorbance has
exactly \(N+1\) columns or is an \(N\)-element list, absorbances
are supposed to be present in the same order as the samples in
x.
If x is a FEEM object, a number specifying the length of the
optical path used when measuring the absorbance, cm.
Otherwise, a vector of them of exactly same length as the number of
FEEMs in x.
If not set, assumed to be \(1\).
No parameters besides those described above are allowed.
An object of the same kind as x, with inner filter effect
corrected.
The formula used is:
$$ I_{corr}(\lambda_{em}, \lambda_{ex}) = % I_{orig}(\lambda_{em}, \lambda_{ex}) 10^{ % \frac{A(\lambda{em}) + A(\lambda{ex})}{2 L_{abs}} % } $$
J.R. Lakowicz, Principles of Fluorescence Spectroscopy, 3rd ed., Springer US, 2006. https://www.springer.com/la/book/9780387312781
D.N. Kothawala, K.R. Murphy, C.A. Stedmon, G.A. Weyhenmeyer, L.J. Tranvik, Inner filter correction of dissolved organic matter fluorescence: Correction of inner filter effects, Limnology and Oceanography: Methods. 11 (2013) 616-630. 10.4319/lom.2013.11.616
# NOT RUN {
data(feems)
dataset <- feemcube(feems, FALSE)
str(dataset)
str(absorp)
plot(feemife(dataset,absorp) / dataset)
# }
Run the code above in your browser using DataLab