Description Usage Arguments Details Value Author(s) See Also Examples
This function takes a list object with potential end-member loadings and extracts those with modes in specified limits to describe them by mean and standard deviation and use these descriptions to propagate the uncertainties to end-member scores.
1 2 3 4 5 6 7 8 9 10 11 12 13 |
em |
|
limits |
|
classunits |
|
amount |
|
l |
|
mc_n |
|
type |
|
qt |
|
cores |
|
plot |
|
... |
Additional arguments passed to |
The function is used to extract potential end-member loadings based on their
mode positions and, optionally the height of the mode class, and use them to
infer mean and stanard deviation of all
end-members that match the group criteria defined by limits
. These
information are then used to model the uncertainty of the corresponding
end-member scores. The function uses input from two preceeding approaches.
In a compact protocol model.em
delivers these data in a predefined
way. In the extended protocol test.robustness
does this.
List
with statistic descriptions of end-member loadings
and scores.
Michael Dietze, Elisabeth Dietze
robust.loadings
, robust.scores
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 | ## Not run:
## load example data set
data(example_X)
## get weight transformation limit vector
l <- get.l(X = X)
## get minimum and maximum number of end-members
q <- get.q(X = X, l = l)
## get all potential model scenarios
EM_pot <- model.EM(X = X, q = q, plot = TRUE)
## define end-member mode class limits
limits <- cbind(c(61, 74, 95, 102),
c(64, 76, 100, 105))
## get robust end-members in the default way, with plot output
rem <- robust.EM(em = EM_pot,
limits = limits,
plot = TRUE)
## get robust end-members by only modelling uncertainty in loadings
robust_EM <- robust.EM(em = EM_pot,
limits = limits,
plot = TRUE)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.