betweencoinertia {ade4}R Documentation

Between-class coinertia analysis

Description

Performs a between-class analysis after a coinertia analysis

Usage

betweencoinertia(obj, fac, scannf = TRUE, nf = 2)
## S3 method for class 'betcoi':
plot(x, xax = 1, yax = 2, ...)
## S3 method for class 'betcoi':
print(x, ...)

Arguments

obj a coinertia analysis (object of class coinertia) obtained by the function coinertia
fac a factor partitioning the rows in classes
scannf a logical value indicating whether the eigenvalues barplot should be displayed
nf if scannf FALSE, an integer indicating the number of kept axes
xax the axis number to be plotted for the x-axis
yax the axis number to be plotted for the y-axis
x an object of class betcoi
... further arguments passed to or from other methods

Details

This analysis is equivalent to do a between-class analysis on each initial dudi, and a coinertia analysis on the two between analyses. This function returns additional outputs for the interpretation.

Value

An object of the class betcoi. Outputs are described by the print function

Author(s)

Stephane Dray dray@biomserv.univ-lyon1.fr and Jean Thioulouse jean.thioulouse@univ-lyon1.fr

References

Franquet E., Doledec S., and Chessel D. (1995) Using multivariate analyses for separating spatial and temporal effects within species-environment relationships. Hydrobiologia, 300, 425–431.

See Also

coinertia, between

Examples

data(meaudret)
pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 4)
pca2 <- dudi.pca(meaudret$fau, scal = FALSE, scan = FALSE, nf = 4)
   
bet1 <- between(pca1, meaudret$plan$sta, scan = FALSE, nf = 2)
bet2 <- between(pca2, meaudret$plan$sta, scan = FALSE, nf = 2)
coib <- coinertia(bet1, bet2, scannf = FALSE)

coi <- coinertia(pca1,pca2,scannf = FALSE,nf = 3)
coi.b <- betweencoinertia(coi,meaudret$plan$sta, scannf = FALSE)
## coib and coi.b are equivalent

plot(coi.b)

Worked out examples


> library(ade4)
> ### Name: betweencoinertia
> ### Title: Between-class coinertia analysis
> ### Aliases: betweencoinertia print.betcoi plot.betcoi
> ### Keywords: multivariate
> 
> ### ** Examples
> 
> data(meaudret)
> pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 4)
> pca2 <- dudi.pca(meaudret$fau, scal = FALSE, scan = FALSE, nf = 4)
> 
> bet1 <- between(pca1, meaudret$plan$sta, scan = FALSE, nf = 2)
> bet2 <- between(pca2, meaudret$plan$sta, scan = FALSE, nf = 2)
> coib <- coinertia(bet1, bet2, scannf = FALSE)
> 
> coi <- coinertia(pca1,pca2,scannf = FALSE,nf = 3)
> coi.b <- betweencoinertia(coi,meaudret$plan$sta, scannf = FALSE)
> ## coib and coi.b are equivalent
> 
> plot(coi.b)

> 
> 
> 
>
[Package ade4 Index]