EMMAgeo (version 0.9.7)

robust.loadings: Extract robust end-member loadings

Description

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.

Usage

robust.loadings(
  em,
  limits,
  classunits,
  amount,
  type = "mean",
  qt = c(0.25, 0.75),
  plot = FALSE,
  ...
)

Arguments

em

List of class "EMMAgeo_empot", i.e. the outout of model.EM() or test.robustness(), containing potential end-members, both in unscaled and rescaled version as well as further parameters.

limits

Numeric matrix with two columns, defining the class limits for the robust end-members to calculate. The first column defines the lower limits, the second column the upper limits. End-members are organised in rows.

classunits

Numeric vector, optional class units (e.g. micrometers or phi-units) of the same length as the number of (grain-size) classes per sample.

amount

Numeric matrix with two columns, defining the minimum and maximum amount of the modal class for each end-member.

type

Character scalar, type of statistics. One out of "mean" and "median". Default is "mean".

qt

Numeric vector of length two, quantiles to describe end-member loadings Default is c(0.25, 0.75) (i.e., the quartile range).

plot

Logical scalar, option to enable plot output. Default is FALSE.

…

Additional arguments passed to EMMA and plot.

Value

List with statistic descriptions of unscaled and scaled end-member loadings.

See Also

robust.EM, robust.scores

Examples

# NOT RUN {
## load example data set, potential end-members, output of model.EM()
data(example_EMpot)

## define limits for robust end-members
limits <- cbind(c(61, 74, 95, 102), 
                c(64, 76, 100, 105))

## get robust end-member loadings with plot output
robust_loadings <- robust.loadings(em = EMpot,
                                   limits = limits,
                                   plot = TRUE)
                    
# }