robust.loadings {EMMAgeo}

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.

Author(s)

Michael Dietze, Elisabeth Dietze

See Also

robust.EM, robust.scores

Examples

## 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)