dircor: Directional mantel correlation

Description Usage Arguments Details Value Note Author(s) References Examples

View source: R/dircor.R

Description

Given a two-dimensional vegetation data frame and the x- and y-coordinates of the releves (the rows in the data frame) in geographical space, mantel correlation (function mantel in the vegan package) is evaluated at regular intervals of direction. Direction versus correlation is plotted including 95 percent confidence interval.

Usage

1
2
3
4
5
6
7
dircor(veg, x.axis, y.axis, step,...)
dircor2(veg, x.axis, y.axis, step = 5)

## Default S3 method:
dircor(veg, x.axis, y.axis, step,...)
## S3 method for class 'dircor'
plot(x,...)

Arguments

veg

A data frame of vegetation releves (rows) by species (columns)

x.axis

This is the x-coordinate in geographical space

y.axis

This is the y-coordinate in geographical space

step

The step length in degrees, used to draw the above mentioned function

x

An object of class "dircor"

...

Further variables used for printing

Details

The method presently uses correlation as distance, as.dist((1-cor(t(sveg^2.0)))/2), as distance measure for vegetation releves.

Value

An output list of class "dircor" with at least the following intems:

steps

The explicit steps used in degrees, 0 - 180 degrees

mean.correlation

Mantel correlation at each step

lower.limit

The lower confidence limits

upper.limit

The upper confidence limits

Note

Computation time is fairly long because function mantel() evaluates the confidence limits for each directional step between 0 and 180 degrees. See also mantel in the vegan package.

Author(s)

Otto Wildi

References

Legendre, P. & Fortin, M.-J. 1989. Spatial analysis and ecological modeling. Vegetatio 80: 107–138.

Wildi, O. 2017. Data Analysis in Vegetation Ecology. 3rd ed. CABI, Oxfordshire, Boston.

Examples

1
2
3
4
# vegetation data is taken from sveg
# the x- and y-axes are stored in object ssit
o.dircor<- dircor(sveg,ssit$x.axis,ssit$y.axis,step=40)
plot(o.dircor)

Example output

Loading required package: cluster
Loading required package: labdsv
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-33. For overview type 'help("mgcv-package")'.
This is labdsv 2.0-1
convert existing ordinations with as.dsvord()

Attaching package:labdsvThe following object is masked frompackage:stats:

    density

Loading required package: vegan
Loading required package: permute
Loading required package: lattice
This is vegan 2.5-7
Loading required package: nnet

Attaching package:nnetThe following object is masked frompackage:mgcv:

    multinom

Loading required package: tree
Call:
dircor.default(veg = sveg, x.axis = ssit$x.axis, y.axis = ssit$y.axis, 
    step = 40)

dave documentation built on May 2, 2019, 1:30 p.m.