analogue (version 0.4-0)

RMSEP: Root mean square error of prediction

Description

Calculates or extracts the RMSEP from transfer function models.

Usage

RMSEP(object, ...)
## S3 method for class 'mat':
RMSEP(object, k, weighted = FALSE,
          \dots)
## S3 method for class 'bootstrap.mat':
RMSEP(object, type = c("birks1990", "standard"),
          ...)

Arguments

object

An R object.

k

numeric; the number of analogues to use in calculating the RMSEP. May be missing. If missing, k is extracted from the model using getK.

weighted

logical; Return the RMSEP for the weighted or unweighted model? Defaults is for the unweighted model.

type

The type of RMSEP to return/calculate.

...

Arguments passed to other methods.

Value

• A numeric vector of length 1 that is the RMSEP of object.

Details

There are two forms of RMSEP in common usage. Within palaeoecology, the RMSEP of Birks et al. (1990) is most familiar: $$\mathrm{RMSEP} = \sqrt{s_1^2 + s_2^2}$$ where where $s_1$ is the standard error of the out-of-bag (OOB) residuals and $s_2$ is the mean bias or the mean of the OOB residuals.

In the wider statistical literature, the following form of RMSEP is more commonly used: $$\mathrm{RMSEP} = \sqrt{\frac{\sum_{i=1}^n (y_i - \hat{y}_i)^2}{n}}$$ where $y_i$ are the observed values and $\hat{y}_i$ the transfer function predictions/fitted values.

The first form of RMSEP is returned by default or if type = "birks1990" is supplied. The latter form is returned if type = "standard" is supplied.

The RMSEP for objects of class "mat" is a leave-one-out cross-validated RMSEP, and is calculated as for type = "standard".

References

Birks, H.J.B., Line, J.M., Juggins, S., Stevenson, A.C. and ter Braak, C.J.F. (1990). Diatoms and pH reconstruction. Philosophical Transactions of the Royal Society of London; Series B, 327; 263--278.

See Also

mat, bootstrap

Examples

## continue the RLGH and SWAP example from ?join
example(join)

## fit the MAT model using the squared chord distance measure
swap.mat <- mat(swapdiat, swappH, method = "SQchord")

## Leave-one-out RMSEP for the MAT model
RMSEP(swap.mat)

## bootstrap training set
swap.boot <- bootstrap(swap.mat, n.boot = 100)
swap.boot

## extract the Birks et al (1990) RMSEP
RMSEP(swap.boot)

## Calculate the alternative formulation
RMSEP(swap.boot, type = "standard")