This function calculates an optional residual end-member loading. It uses the modelled end-member loadings as input and evaluates the root of 1 minus the sum of all squared loadings to analyse the remaining variance, e.g. if not all (robust) EMs are included (cf. Dietze et al., 2012). Negative values are set to zero.

1 | ```
residual.EM(Vqn)
``` |

`Vqn` |
Numeric matrix with m robust end-member loadings. |

Numeric vector with residual end-member loading.

Michael Dietze, Elisabeth Dietze

Dietze E, Hartmann K, Diekmann B, IJmker J, Lehmkuhl F, Opitz S, Stauch G, Wuennemann B, Borchers A. 2012. An end-member algorithm for deciphering modern detrital processes from lake sediments of Lake Donggi Cona, NE Tibetan Plateau, China. Sedimentary Geology 243-244: 169-180.

`EMMA`

, `robust.EM`

1 2 3 4 5 6 7 8 9 10 11 12 13 | ```
## Some preparing steps to retrieve only robust end-members
## load example data, i.e. here TR
data(rEM, envir = environment())
## define mean robust end-member loadings
Vqn.rob <- rEM$Vqn.mean
## perform residual end-member loading calculation
Vqn.res <- residual.EM(Vqn.rob)
# Visualisation of the result
plot(NA, xlim = c(1, 80), ylim = c(0, 1))
for(i in 1:4) {lines(Vqn.rob[i,])}
lines(Vqn.res, col = 2)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.