avimedi {ade4}R Documentation

Fauna Table for Constrained Ordinations

Description

avimedi is a list containing the information about 302 sites :
frequencies of 51 bird species ; two factors (habitats and Mediterranean origin).

Usage

data(avimedi)

Format

This list contains the following objects:

fau
is a data frame 302 sites - 51 bird species.
plan
is a data frame 302 sites - 2 factors : reg with two levels Provence (Pr, South of France) and Corsica (Co) ; str with six levels describing the vegetation from a very low matorral (1) up to a mature forest of holm oaks (6).
nomesp
is a vector 51 latin names.

Source

Blondel, J., Chessel, D., & Frochot, B. (1988) Bird species impoverishment, niche expansion, and density inflation in mediterranean island habitats. Ecology, 69, 1899–1917.

Examples

## Not run: 
data(avimedi)
par(mfrow = c(2,2))
coa1 <- dudi.coa(avimedi$fau, scan = FALSE, nf = 3)
s.class(coa1$li,avimedi$plan$str:avimedi$plan$reg,
    sub = "Correspondences Analysis")
bet1 <- between(coa1, avimedi$plan$str, scan = FALSE)
s.class(bet1$ls, avimedi$plan$str,
    sub = "Between Analysis")
wit1 <- within(coa1, avimedi$plan$reg, scan=FALSE)
s.class(wit1$li, avimedi$plan$str,
    sub = "Within Analysis")
pcaiv1 <- pcaiv(coa1, avimedi$plan, scan = FALSE)
s.match(pcaiv1$li, pcaiv1$ls, clab = 0,
    sub = "Canonical Correspondences Analysis")
s.class(pcaiv1$li, avimedi$plan$str:avimedi$plan$reg, 
    add.plot = TRUE)
par(mfrow=c(1,1))

## End(Not run)

Worked out examples


> library(ade4)
> ### Name: avimedi
> ### Title: Fauna Table for Constrained Ordinations
> ### Aliases: avimedi
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(avimedi)
> par(mfrow = c(2,2))
> coa1 <- dudi.coa(avimedi$fau, scan = FALSE, nf = 3)
> s.class(coa1$li,avimedi$plan$str:avimedi$plan$reg,
+     sub = "Correspondences Analysis")
> bet1 <- between(coa1, avimedi$plan$str, scan = FALSE)
> s.class(bet1$ls, avimedi$plan$str,
+     sub = "Between Analysis")
> wit1 <- within(coa1, avimedi$plan$reg, scan=FALSE)
> s.class(wit1$li, avimedi$plan$str,
+     sub = "Within Analysis")
> pcaiv1 <- pcaiv(coa1, avimedi$plan, scan = FALSE)
> s.match(pcaiv1$li, pcaiv1$ls, clab = 0,
+     sub = "Canonical Correspondences Analysis")
> s.class(pcaiv1$li, avimedi$plan$str:avimedi$plan$reg, 
+     add.plot = TRUE)
> par(mfrow=c(1,1))
> 
> 
> 
> 

[Package ade4 Index]