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

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