rioja (version 0.9-6)

MR: Multiple regression

Description

Functions for reconstructing (predicting) environmental values from biological assemblages using multiple regression.

Usage

MR(y, x, check.data=TRUE, lean=FALSE, ...)

MR.fit(y, x, lean=FALSE)

## S3 method for class 'MR':
predict(object, newdata=NULL, sse=FALSE, nboot=100,
      match.data=TRUE, verbose=TRUE, \dots)

## S3 method for class 'MR':
crossval(object, cv.method="loo", verbose=TRUE, ngroups=10,
      nboot=100, h.cutoff=0, h.dist=NULL, \dots)

## S3 method for class 'MR':
performance(object, \dots)

## S3 method for class 'MR':
print(x, \dots)

## S3 method for class 'MR':
summary(object, full=FALSE, \dots)

## S3 method for class 'MR':
plot(x, resid=FALSE, xval=FALSE, xlab="", 
      ylab="", ylim=NULL, xlim=NULL, add.ref=TRUE,
      add.smooth=FALSE, \dots)

## S3 method for class 'MR':
residuals(object, cv=FALSE, \dots)

## S3 method for class 'MR':
coef(object, \dots)

## S3 method for class 'MR':
fitted(object, \dots)

Arguments

y

a data frame or matrix of biological abundance data.

x, object

a vector of environmental values to be modelled or an object of class wa.

newdata

new biological data to be predicted.

check.data

logical to perform simple checks on the input data.

match.data

logical indicate the function will match two species datasets by their column names. You should only set this to FALSE if you are sure the column names match exactly.

lean

logical to exclude some output from the resulting models (used when cross-validating to speed calculations).

full

logical to show head and tail of output in summaries.

resid

logical to plot residuals instead of fitted values.

xval

logical to plot cross-validation estimates.

xlab, ylab, xlim, ylim

additional graphical arguments to plot.wa.

add.ref

add 1:1 line on plot.

add.smooth

add loess smooth to plot.

cv.method

cross-validation method, either "loo", "lgo" or "bootstrap".

verbose

logical to show feedback during cross-validaton.

nboot

number of bootstrap samples.

ngroups

number of groups in leave-group-out cross-validation, or a vector contain leave-out group menbership.

h.cutoff

cutoff for h-block cross-validation. Only training samples greater than h.cutoff from each test sample will be used.

h.dist

distance matrix for use in h-block cross-validation. Usually a matrix of geographical distances between samples.

sse

logical indicating that sample specific errors should be calculated.

cv

logical to indicate model or cross-validation residuals.

...

additional arguments.

Value

• Function MR returns an object of class MR with the following named elements:
• coefficients
  species coefficients (the updated "optima").
• fitted.values
  fitted values for the training set.
• call
  original function call.
• x
  environmental variable used in the model.
• Function crossval also returns an object of class MR and adds the following named elements:
• predicted
  predicted values of each training set sample under cross-validation.
• residuals.cv
  prediction residuals.
• If function predict is called with newdata=NULL it returns the fitted values of the original model, otherwise it returns a list with the following named elements:
• fit
  predicted values for newdata.
• if sample specific errors were generated the list will also include:
• fit.boot
  mean of the bootstrap estimates of newdata.
• v1
  squared standard error of the bootstrap estimates for each new sample.
• v2
  mean squared error for the training set samples, across all bootstrap samples.
• SEP
  standard error of prediction, calculated as the square root of v1 + v2.
• Function performance returns a matrix of performance statistics for the MR model. See performance, for a description of the summary.

Details

Function MR performs multiple regrssion. It is a wrapper to lm. Function predict predicts values of the environemntal variable for newdata or returns the fitted (predicted) values from the original modern dataset if newdata is NULL. Variables are matched between training and newdata by column name (if match.data is TRUE). Use compare.datasets to assess conformity of two species datasets and identify possible no-analogue samples. MR has methods fitted and rediduals that return the fitted values (estimates) and residuals for the training set, performance, which returns summary performance statistics (see below), coef which returns the species coefficients, and print and summary to summarise the output. MR also has a plot method that produces scatter plots of predicted vs observed measurements for the training set.

See Also

WA, MAT, performance, and compare.datasets for diagnostics.

Examples

data(IK)
spec <- IK$spec
SumSST <- IK$env$SumSST
core <- IK$core

fit <- IKFA(spec, SumSST)
fit
# cross-validate model
fit.cv <- crossval(fit, cv.method="lgo")
# How many components to use?
screeplot(fit.cv)

#predict the core
pred <- predict(fit, core, npls=2)

#plot predictions - depths are in rownames
depth <- as.numeric(rownames(core))
plot(depth, pred$fit[, 2], type="b")

# fit using only factors 1, 2, 4, & 5
# and using polynomial terms
# as Imbrie & Kipp (1971)
fit2 <- IKFA(spec, SumSST, ccoef=c(1, 2, 4, 5), IsPoly=TRUE)
fit2.cv <- crossval(fit2, cv.method="lgo")
screeplot(fit2.cv)

# predictions with sample specific errors
# takes approximately 1 minute to run
pred <- predict(fit, core, sse=TRUE, nboot=1000)
pred