vegan (version 1.0-1)

wascores: Weighted Averages Scores for Species

Description

Computes Weighted Averages scores of species for ordination configuration or for environmental variables.

Usage

wascores(x, w, expand=FALSE)

Arguments

x

Environmental variables or ordination scores.

w

Weights: species abundances.

expand

Not yet implemented.

Value

• Returns a matrix where species define rows and environmental variables define columns.

Details

Computes weighted averages.

See Also

isoMDS

Examples

library(MASS)
library(mva)
data(varespec)
data(varechem)
vare.dist <- vegdist(wisconsin(varespec))
vare.mds <- isoMDS(vare.dist)
vare.points <- postMDS(vare.mds$points, vare.dist)
vare.wa <- wascores(vare.points, varespec)
eqscplot(vare.points, pch="+")
text(vare.wa, rownames(vare.wa), cex=0.8, col="blue")
## Omit rare species (frequency <= 4)
freq <- apply(varespec>0, 2, sum)
eqscplot(vare.points, pch="+")
text(vare.wa[freq > 4,], rownames(vare.wa)[freq > 4], cex=0.8, col="blue")